The time-relaxation limit for weak solutions to the quantum hydrodynamics system

1. Problem Statement

The paper investigates the time-relaxation limit ($\tau \to 0$) of the Quantum Hydrodynamics (QHD) system with a collisional damping term, posed on a one-dimensional torus ($\mathbb{T} = \mathbb{R}/\mathbb{Z}$).

The QHD system (1.1) describes the evolution of mass density $\rho$ and momentum density $J$:
$$\begin{cases} \partial_t \rho + \partial_x J = 0 \\ \partial_t J + \partial_x \left(\frac{J^2}{\rho}\right) + \partial_x p(\rho) + \rho \partial_x V = \frac{1}{2}\rho \partial_x \left( \frac{\partial_x^2 \sqrt{\rho}}{\sqrt{\rho}} \right) - \frac{1}{\tau} J \\ -\partial_x^2 V = \rho - C(x) \end{cases}$$
where $V$ is the self-consistent electric potential (Poisson equation), $C(x)$ is the background doping profile, $p(\rho)$ is the isentropic pressure, and $\tau > 0$ is the momentum relaxation time.

The goal is to rigorously prove that as $\tau \to 0$, the solutions of the QHD system converge to the solutions of the Quantum Drift-Diffusion (QDD) equation (1.5):
$$\partial_t \bar{\rho} + \partial_x \left[ \frac{1}{2}\bar{\rho}\partial_x \left( \frac{\partial_x^2 \sqrt{\bar{\rho}}}{\sqrt{\bar{\rho}}} \right) - \partial_x p(\bar{\rho}) - \bar{\rho}\partial_x \bar{V} \right] = 0$$
Key Challenges:

• Compactness: Standard energy bounds are insufficient to control the weak limit of the nonlinear quantum stress tensor term $(\partial_x \sqrt{\rho})^2$, which may suffer from oscillations and concentration.
• Vacuum: The presence of vacuum regions ($\rho=0$) creates singularities in the chemical potential and phase functions.
• Initial Data: Previous works often required "well-prepared" initial data (where the initial momentum matches the limiting drift-diffusion profile). This paper aims to handle general initial data, allowing for an initial layer.

2. Methodology

The authors employ a multi-faceted approach combining hydrodynamic analysis, wave function lifting, and entropy methods.

A. GCP Solutions (Generalized Chemical Potential)

To overcome the lack of compactness in standard weak solutions, the authors introduce a class of GCP solutions. These are weak solutions with finite energy and a bounded generalized chemical potential functional $I(t)$:
$$I(t) = \int_{\mathbb{T}} \frac{\rho}{2} (\mu^2 + \sigma^2) dx$$
where $\mu$ is the chemical potential and $\sigma = \partial_t \log \sqrt{\rho}$.

• Assumption: The analysis restricts to solutions where the density remains strictly positive ($\rho \ge \delta > 0$) to avoid vacuum singularities.
• Role: The finiteness of $I(t)$ provides higher-order regularity (specifically $H^2$ bounds on $\sqrt{\rho}$) necessary to pass to the limit in the quantum stress tensor.

B. Wave Function Lifting and Polar Factorization

The authors exploit the formal equivalence between the QHD system and the Schrödinger-Langevin equation (1.21):
$$i\partial_t \psi + \frac{1}{2}\partial_x^2 \psi = f'(|\psi|^2)\psi + \frac{1}{\tau} S \psi + V \psi$$
where $\psi = \sqrt{\rho} e^{iS}$.

• They use wave function lifting and polar factorization techniques (developed in their previous works) to rigorously establish the equivalence between the hydrodynamic variables $(\rho, J)$ and the wave function $\psi$ for strictly positive densities.
• This allows them to prove global well-posedness for the QHD system by first proving it for the Schrödinger-Langevin equation in $H^2$.

C. Rescaling and Entropy Estimates

The system is rescaled using $t' = \tau t$ and $(\rho_\tau, J_\tau) = (\rho, \frac{1}{\tau}J)$.

• Physical Entropy: The authors utilize the physical entropy $H(\rho) = \int \rho \log(\rho/M_0) dx$.
• Dissipation: For small $\tau$ and specific internal energy forms (e.g., $f(\rho) = (\rho - M_0)^{2n}$), they prove that the total energy and the functional $I(t)$ decay exponentially. This establishes the dissipative nature of the collisional QHD system.

D. Relative Entropy Method

To prove the convergence rate, the authors use the relative entropy method:
$$H(\rho_\tau | \bar{\rho}) = \int_{\mathbb{T}} \left( g(\rho_\tau) - g(\bar{\rho}) - g'(\bar{\rho})(\rho_\tau - \bar{\rho}) \right) dx$$
where $g(s) = s \log(s/M_0)$. By deriving a differential inequality for the time derivative of the relative entropy and controlling the remainder terms using the uniform bounds of GCP solutions, they establish the convergence rate.

3. Key Contributions and Results

1. Global Well-Posedness of GCP Solutions (Theorem 2 & 4)

• Proved the existence and uniqueness of global-in-time GCP solutions for the QHD system with strictly positive initial density.
• Established the equivalence between the QHD system and the Schrödinger-Langevin equation for these solutions.
• Derived uniform bounds on the higher-order functional $I(t)$, ensuring the density remains bounded away from vacuum.

2. Rigorous Relaxation-Time Limit (Theorem 6)

• Proved that as $\tau \to 0$, the sequence of GCP solutions $(\rho_\tau, J_\tau)$ converges to a weak solution $\bar{\rho}$ of the QDD equation.
• Novelty: The result holds for general initial data (not necessarily well-prepared), meaning an initial layer is allowed. The convergence of momentum occurs only for $t > 0$.
• Explicit Convergence Rate: The paper provides an explicit convergence rate of order $O(\tau)$:
  $$\|\rho_\tau - \bar{\rho}\|_{L^\infty_t L^2_x} + \|\sqrt{\rho_\tau}(\partial_x^2 \log \sqrt{\rho_\tau} - \partial_x^2 \log \sqrt{\bar{\rho}})\|_{L^2_{t,x}} \le C \tau$$

3. Dissipation and Exponential Decay (Theorem 7)

• Demonstrated that for small $\tau$ and specific internal energy/doping profiles, the GCP solutions exhibit exponential decay of the physical entropy, total energy, and the functional $I(t)$.
• This highlights the dissipative mechanism introduced by the collision term, breaking the time-reversibility of the standard quantum mechanics.

4. Existence of QDD Solutions

• As a by-product, the analysis constructs global weak solutions (and dissipative solutions under convexity of $f$) to the limiting QDD equation as strong limits of the QHD solutions.

4. Significance

• Mathematical Rigor: This work provides one of the first rigorous justifications of the relaxation limit for weak solutions of quantum hydrodynamic systems without assuming the smoothness of the limiting solution or well-prepared initial data.
• Handling Nonlinearity: By introducing the GCP framework and utilizing the wave function lifting, the authors successfully control the highly nonlinear quantum stress tensor term, which is the main obstacle in such limits.
• Physical Relevance: The results are crucial for modeling carrier transport in semiconductor devices at the nanoscale, where quantum effects and collisional relaxation are both significant. The explicit convergence rate allows for better error estimation in numerical simulations and asymptotic modeling.
• Dissipative Structure: The identification of the exponential decay of entropy for the collisional QHD system bridges the gap between conservative quantum dynamics and dissipative classical/semi-classical transport models.

In summary, the paper establishes a robust mathematical framework for the transition from Quantum Hydrodynamics to Quantum Drift-Diffusion, overcoming significant analytical hurdles related to compactness, vacuum, and initial layers.