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

1. Problématique et Contexte

L'article se concentre sur l'analyse des solutions faibles du système d'hydrodynamique quantique (QHD) unidimensionnel, incluant un terme de collision, posé sur un tore $T = \mathbb{R}/\mathbb{Z}$. Le système est décrit par les équations (1.1) :
$$\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_{xx}^2 \sqrt{\rho}}{\sqrt{\rho}}\right) - \frac{1}{\tau} J \\ -\partial_{xx}^2 V = \rho - C(x) \end{cases}$$
où $\rho$ est la densité de masse, $J$ la densité de quantité de mouvement, $V$ le potentiel électrique auto-cohérent (satisfaisant l'équation de Poisson), $p(\rho)$ la pression isentropique, et $\tau > 0$ le temps de relaxation de la quantité de mouvement.

Objectif principal : Établir rigoureusement la limite de temps de relaxation lorsque $\tau \to 0$. Formellement, ce système converge vers l'équation de dérive-diffusion quantique (QDD) (équation 1.5) :
$$\partial_t \bar{\rho} + \partial_x \left[ \frac{1}{2}\bar{\rho}\partial_x \left(\frac{\partial_{xx}^2 \sqrt{\bar{\rho}}}{\sqrt{\bar{\rho}}}\right) - \partial_x p(\bar{\rho}) - \bar{\rho}\partial_x \bar{V} \right] = 0$$

Défis majeurs :

• Les bornes d'énergie standard ne suffisent pas à assurer la compacité nécessaire pour passer à la limite dans le terme de tension quantique non linéaire (le terme de type $(\partial_x \sqrt{\rho})^2$), qui peut subir des oscillations fortes ou des concentrations.
• La présence de régions de vide ($\rho = 0$) complique la définition des solutions faibles et la régularité nécessaire.
• Les travaux antérieurs nécessitaient souvent des données initiales « bien préparées » (well-prepared) ou la régularité des solutions de l'équation limite, ce qui limite la généralité des résultats.

2. Méthodologie

Les auteurs développent une approche combinant plusieurs outils avancés d'analyse mathématique :

1. Solutions GCP (Generalized Chemical Potential) :
   Pour surmonter le manque de compacité des solutions faibles classiques, les auteurs introduisent une classe de solutions plus régulières, appelées solutions GCP. Ces solutions sont caractérisées par la finitude d'un fonctionnel d'ordre supérieur $I(t)$, lié au potentiel chimique $\mu$ :
   $$I(t) = \int_T \frac{\rho}{2} (\mu^2 + \sigma^2) dx$$
   où $\mu$ est le potentiel chimique et $\sigma = \partial_t \log \sqrt{\rho}$. Cette fonctionnel fournit des estimations de régularité supplémentaires (type $H^2$ pour $\sqrt{\rho}$) essentielles pour la compacité.

2. Lifting de la fonction d'onde et Factorisation Polaire :
   L'analyse repose sur l'équivalence formelle et rigoureuse entre le système hydrodynamique (1.1) et une équation de Schrödinger-Langevin (1.21) :
   $$i\partial_t \psi + \frac{1}{2}\partial_{xx}^2 \psi = f'(|\psi|^2)\psi + \frac{1}{\tau} S \psi + V \psi$$
   Grâce aux techniques de lifting de la fonction d'onde et de factorisation polaire (développées dans des travaux précédents des auteurs), ils établissent une correspondance biunivoque entre les solutions hydrodynamiques à densité strictement positive et les solutions de l'équation de Schrödinger-Langevin. Cela permet d'utiliser la structure de l'équation de Schrödinger pour prouver l'existence globale.

3. Estimations d'Entropie et Dissipation :
   Les auteurs utilisent une entropie physique $H(\rho) = \int \rho \log(\rho/M_0) dx$ et une fonctionnelle combinée $F_\tau(t) = H(\rho_\tau) + E_\tau(t) + c_1 I_\tau(t)$. Ils démontrent que pour $\tau$ suffisamment petit et sous des hypothèses spécifiques sur l'énergie interne et le profil de dopage, cette fonctionnelle décroît exponentiellement, garantissant la stabilité et la dissipation du système.

4. Méthode d'Entropie Relative :
   Pour obtenir un taux de convergence explicite, ils utilisent la méthode de l'entropie relative $H(\rho_\tau | \bar{\rho})$ entre la solution du système QHD et la solution limite QDD.

3. Résultats Clés

L'article établit les résultats théoriques suivants :

• Existence Globale et Bien-posé (Théorème 2 & 4) :
  Les auteurs prouvent l'existence et l'unicité de solutions globales en temps (GCP solutions) pour le système QHD avec des données initiales à densité strictement positive. Ils démontrent également la bien-posé globale de l'équation de Schrödinger-Langevin associée dans l'espace $H^2$. Ces solutions satisfont des lois de conservation de l'énergie et des bornes uniformes sur le fonctionnel $I(t)$.

• Limite de Temps de Relaxation (Théorème 6) :
  C'est le résultat central. Sans hypothèse de données initiales « bien préparées » (c'est-à-dire même si la condition initiale de la quantité de mouvement ne correspond pas exactement à la loi de comportement limite), les auteurs prouvent que la suite de solutions $(\rho_\tau, J_\tau)$ converge vers une solution faible $\bar{\rho}$ de l'équation QDD.

  • Taux de convergence : Ils établissent un taux de convergence explicite en $O(\tau)$ :
    $$\|\rho_\tau - \bar{\rho}\|_{L^\infty_t L^2_x} + \|\sqrt{\rho_\tau}(\partial_{xx}^2 \log \sqrt{\rho_\tau} - \partial_{xx}^2 \log \sqrt{\bar{\rho}})\|_{L^2_{t,x}} \leq C \tau$$
  • La convergence de la densité de quantité de mouvement est valable pour $t > 0$ (phénomène de couche initiale admis).

• Dissipation Exponentielle (Théorème 7) :
  Pour des énergies internes spécifiques (de type $f(\rho) = (\rho - M_0)^{2n}$) et un profil de dopage constant, ils montrent que les solutions du système résolu (avec $\tau$ petit) possèdent une dissipation d'entropie et d'énergie exponentielle en temps, généralisant les résultats connus pour l'équation DLSS (Derrida-Lebowitz-Speer-Spohn).

4. Contributions et Significativité

• Rigueur sans hypothèses restrictives : Contrairement à la littérature précédente, ce travail ne nécessite ni la régularité des solutions de l'équation limite, ni des données initiales bien préparées. Cela rend le résultat applicable à des situations physiques plus réalistes où l'équilibre initial n'est pas parfait.
• Nouvelle classe de solutions : L'introduction et l'exploitation des solutions GCP permettent de contourner les difficultés de compacité liées aux termes quantiques non linéaires, offrant un cadre robuste pour l'analyse des limites hydrodynamiques quantiques.
• Taux de convergence explicite : La fourniture d'un taux de convergence linéaire en $\tau$ est une avancée significative par rapport aux résultats qualitatifs antérieurs.
• Lien Schrödinger-Hydrodynamique : La preuve rigoureuse de l'équivalence entre le système QHD avec collision et l'équation de Schrödinger-Langevin pour des densités non nulles ouvre la voie à l'utilisation d'outils d'analyse d'équations d'ondes pour étudier des systèmes hydrodynamiques quantiques dissipatifs.

En résumé, ce papier fournit une justification mathématique complète et rigoureuse de la transition du modèle hydrodynamique quantique avec collision vers le modèle de dérive-diffusion quantique, en surmontant les obstacles analytiques liés à la régularité et à la compacité grâce à une approche innovante combinant la factorisation polaire, les fonctionnels d'ordre supérieur et l'entropie relative.