Large-Norm Solutions and the Relaxation-Time Limit for Quantum Hydrodynamics on the Two-Dimensional Torus

Technical Summary: Large-Norm Solutions and the Relaxation-Time Limit for Quantum Hydrodynamics on the Two-Dimensional Torus

Problem Statement
This paper addresses the mathematical analysis of the collisional quantum hydrodynamic (QHD) system arising in semiconductor device modeling on the two-dimensional torus ($T^2$). The system couples the continuity and momentum equations with a Poisson equation for the self-consistent electric potential and includes a quantum Bohm potential term. The primary challenges tackled are:

1. Global Well-Posedness for Large Norms: Establishing the existence and uniqueness of global weak solutions for initial data with arbitrarily large norms, provided the mass is sufficiently large and the density remains strictly positive. Previous results often relied on smallness assumptions for initial data.
2. Relaxation-Time Limit: Rigorously justifying the limit as the relaxation time parameter $\tau \to 0$, where the QHD system converges to the quantum drift-diffusion (QDD) system. A specific difficulty is handling general initial data that are not "well-prepared" (i.e., where the initial momentum does not satisfy the constitutive relation of the limiting system), which leads to the formation of an initial layer.

Methodology
The analysis relies on a functional framework centered on Generalized Chemical Potential (GCP) solutions, a concept introduced in the authors' previous work [7]. The core methodology involves:

• Wave Function Lifting and Polar Factorization: The authors establish a rigorous equivalence between the hydrodynamic variables $(\rho, J)$ and a complex wave function $\psi$ via the Madelung transformation ($\rho = |\psi|^2, J = \text{Im}(\bar{\psi}\nabla\psi)$). This allows the QHD system to be analyzed through a corresponding Schrödinger–Langevin equation (1.17).
• Combined Functionals: The stability analysis utilizes a combination of the total energy functional $E(t)$, the higher-order GCP functional $I(t)$, and the physical entropy $H(\rho)$. The functional $I(t)$, defined via the $L^2$ norms of the chemical potential $\mu$ and the time derivative of the log-density $\sigma$, controls higher-order regularity.
• Logarithmic Sobolev-Type Inequality: A critical technical ingredient is a logarithmic Sobolev-type inequality (Lemma 17) adapted to the torus. This inequality controls the $L^\infty$ oscillations of the density in terms of the energy and the GCP functional, which is essential for preventing vacuum formation in higher dimensions where standard energy estimates are insufficient.
• Relative Entropy Method: For the relaxation-time limit, the authors employ the relative entropy method to derive explicit convergence rates. This approach compares the solution of the rescaled QHD system with the solution of the limiting QDD system without requiring the limiting solution to be smooth or the initial data to be well-prepared.

Key Contributions and Results

1. Global Well-Posedness (Theorem 1):
   The paper proves the global existence and uniqueness of GCP weak solutions with strictly positive density ($\inf \rho \ge \delta > 0$).

   • Large Data: The result does not require small initial norms. Instead, it imposes a condition (1.14) relating the initial energy $E_0$ and GCP functional $I_0$ to the total mass $M_0$. As $M_0 \to \infty$, the allowable size of $E_0$ and $I_0$ grows, accommodating large-norm solutions.
   • Stability and Decay: The combined functional $F(t) = H(\rho) + E(t) + c_1 I(t)$ is shown to decay exponentially for small $\tau$ (Theorem 4), establishing stability and convergence to equilibrium.

2. Schrödinger–Langevin Well-Posedness (Theorem 3):
   As an auxiliary result, the paper establishes the global well-posedness of the Schrödinger–Langevin equation (1.17) in $H^2(T^2)$ for initial data with strictly positive modulus. This is achieved via a Picard iteration scheme and the wave function lifting technique.

3. Relaxation-Time Limit (Theorem 6):
   The authors rigorously justify the convergence of the rescaled QHD system to the QDD system as $\tau \to 0$.

   • General Initial Data: Unlike previous works, this result does not assume well-prepared initial data. Consequently, an initial layer is allowed to form, and convergence of the momentum density is established for $t' > 0$.
   • Explicit Rate: The paper provides an explicit convergence rate of order $O(\tau)$ in the $L^\infty(0, T; L^2)$ norm for the density and a corresponding rate for the second-order derivatives of the log-density.
   • Compactness: The proof relies on compactness techniques derived from uniform bounds on the physical entropy and the GCP functional, avoiding the need for smoothness assumptions on the limiting solution.

Significance
The paper claims significance in extending the analysis of quantum hydrodynamic systems to the two-dimensional torus with large initial data and strictly positive densities. By introducing the GCP framework and utilizing the wave function lifting method, the authors overcome the difficulties associated with vacuum formation and strong oscillations that typically plague energy-based compactness arguments in higher dimensions. Furthermore, the rigorous justification of the relaxation-time limit with an explicit convergence rate, without requiring well-prepared initial data, provides a more robust theoretical foundation for the transition from quantum hydrodynamic to quantum drift-diffusion models in semiconductor physics. The results also highlight the utility of the Schrödinger–Langevin equation as a tool for analyzing the underlying hydrodynamic system.