Resolvent estimates for the Schrödinger operator with $L^\infty$ electric and magnetic potentials and applications to the local energy decay

1. Problématique et Contexte

L'objectif principal de cet article est d'étendre les résultats récents concernant la décroissance de l'énergie locale pour les solutions de l'équation des ondes associée à un opérateur de Schrödinger magnétique.

• Opérateur étudié : L'opérateur $P = (i\nabla + b(x))^2 + V(x)$ défini sur $L^2(\Omega)$, où $\Omega = \mathbb{R}^d \setminus \mathcal{O}$ ($\mathcal{O}$ étant un obstacle borné ou vide).
• Potentiels : Les auteurs considèrent des potentiels électriques $V \in L^\infty(\mathbb{R}^d; \mathbb{R})$ et magnétiques $b \in L^\infty(\mathbb{R}^d; \mathbb{R}^d)$.
• Hypothèses géométriques :
  • Cas a) : $\Omega = \mathbb{R}^d$ ($d \ge 3$) avec un potentiel magnétique non trivial.
  • Cas b) : $\Omega$ est le complémentaire d'un obstacle non piégeant (non-trapping) avec $b \equiv 0$.
• Limitation des travaux précédents : L'article [13] (référence clé) traitait de potentiels à décroissance exponentielle ($|V| + |b| \le Ce^{-c|x|}$), permettant une continuation analytique du résolvant pondéré.
• Nouvelle contribution : Cet article généralise ces résultats à une classe beaucoup plus large de potentiels à décroissance sous-exponentielle (ou plus rapide que toute puissance, mais plus lente que l'exponentielle pure dans certains cas), caractérisée par une fonction de poids $\Theta$.

2. Méthodologie

La stratégie repose sur l'obtention d'estimations précises pour le résolvant pondéré et ses dérivées par rapport au paramètre spectral $\lambda$, uniformes en $\lambda$ et en l'ordre de dérivation $k$.

A. Classe de potentiels admissibles

Les potentiels sont supposés satisfaire une condition de décroissance de la forme :
$$|V(x)| + |b(x)| \le C \Theta(c\langle x \rangle)$$
où $\langle x \rangle = (1+|x|^2)^{1/2}$ et $\Theta$ est une fonction décroissante satisfaisant des conditions de régularité de type Gevrey.

• Exemple typique : $\Theta(r) = e^{-(r+1)^s}$ avec $0 < s \le 1$.
• La fonction $\Theta$ doit satisfaire des inégalités impliquant une suite croissante $m_k$ (liée aux dérivées de $\Theta$) :
  • Condition de croissance des dérivées : $|\partial_r^j \Theta| \le \tilde{C} \Theta$.
  • Condition de stabilité sous produit : $\Theta(r_1)\Theta(r_2) \le C_1 \Theta(C_2(r_1+r_2))$.
  • Conditions sur la suite $m_k$ (liées à la factorielle et aux inégalités de type Gevrey) : $k! + \tilde{m}_k + m_{k+2} + \sqrt{m_{2k}} \le C_3^{k+1} m_k$ et une inégalité de type Leibniz pour les produits.

B. Estimations du Résolvant (Théorème 1.1)

Le cœur de la preuve consiste à établir des bornes pour les dérivées du résolvant pondéré $\mu(P-\lambda^2)^{-1}\mu$, où $\mu(x) = \sqrt{\Theta(c\langle x \rangle)}$.

• Résultat clé : Pour tout $k \ge 0$ et $\lambda \in \mathbb{R}$ ($|\lambda| \ge \delta$), on a :
  $$\left\| \frac{d^k}{d\lambda^k} \mu \nabla^\ell (P-\lambda^2)^{-1} \mu \right\| \le C^{k+1} m_k (|\lambda|+1)^{\ell-1}$$
  où $\ell \in \{0, 1\}$.
• Techniques utilisées :
  1. Résolvant libre : Estimation des dérivées du résolvant de l'opérateur de Laplace pondéré en utilisant la formule de représentation intégrale et le principe de Huygens (pour les dimensions impaires) ou des propriétés de noyaux analytiques (pour les dimensions paires).
  2. Identités de résolvant : Utilisation d'identités de type "resolvent identity" pour exprimer le résolvant perturbé $P$ en fonction du résolvant libre $P_0$ et des termes de perturbation.
  3. Lemmes d'algèbre d'opérateurs : Démonstration que la classe de fonctions à valeurs opérateurs satisfaisant ces bornes de dérivées (liées à $m_k$) est stable par produit et inversion (Lemmes 3.1 et 3.2). Cela permet de passer du résolvant libre au résolvant perturbé.
  4. Principe d'absorption limite : Utilisation des estimations de haute fréquence (hypothèse non-piégeante) et du principe d'absorption limite pour traiter les basses fréquences et assurer la continuité jusqu'à l'axe réel.

C. Décroissance de l'énergie locale (Théorème 1.2)

Une fois les estimations du résolvant établies, la décroissance temporelle est obtenue par une transformation de Fourier inverse et une optimisation sur l'ordre de dérivation $k$ en fonction du temps $t$.

• On introduit une fonction de coupure $\psi$ dans l'espace spectral, construite à partir de fonctions de classe Gevrey.
• En choisissant $k(t)$ comme le plus grand entier tel que $m_k^{1/k} \le t$, on obtient une décroissance de la forme $e^{-k(c_0 t)}$.

3. Résultats Principaux

1. Extension des estimations de résolvant : Les auteurs prouvent que les estimations uniformes sur les dérivées du résolvant pondéré, valables pour les potentiels exponentiels, s'étendent aux potentiels à décroissance $\Theta$ satisfaisant les conditions (1.7)-(1.10). Cela inclut des décroissances de type $e^{-|x|^s}$ ($0 < s < 1$).
2. Décroissance de l'énergie locale : Pour les potentiels satisfaisant $|V(x)| + |b(x)| \le C e^{-c|x|^s}$ avec $0 < s < 1$, la décroissance de l'énergie locale est de la forme :
   $$E_{loc}(t) \le C_0 e^{-c_0 t^s}$$
   Cela généralise le résultat exponentiel ($s=1$) obtenu dans [13] à des décroissances plus lentes (sous-exponentielles).
3. Cas des dimensions paires et impaires :
   • En dimension impaire, sous l'hypothèse de non-piégeage et d'absence de résonance à zéro, les estimations valent pour tout $\lambda \in \mathbb{R}$.
   • En dimension paire, les résultats sont obtenus sur le revêtement de Riemann du logarithme ou avec des coupures spectrales appropriées.
4. Obstacles non-piégeants : Le résultat est étendu au cas d'un obstacle non-piégeant en dimension $d \ge 2$, à condition que le potentiel magnétique soit nul ($b \equiv 0$).

4. Signification et Impact

• Généralisation théorique : Ce travail brise la barrière de la décroissance exponentielle stricte. Il montre que la méthode de continuation analytique (ou ses équivalents via les estimations Gevrey) fonctionne pour une classe beaucoup plus large de potentiels, reliant la régularité de la décroissance spatiale des potentiels à la régularité Gevrey du résolvant.
• Lien Géométrie-Dynamique : Il renforce le lien entre la géométrie de l'obstacle (condition non-piégeante) et la vitesse de dissipation de l'énergie, même en présence de perturbations magnétiques et électriques complexes.
• Outils techniques : La démonstration des lemmes de stabilité (produit et inversion) pour les classes de fonctions à dérivées contrôlées par une suite $m_k$ fournit un outil puissant applicable à d'autres problèmes d'EDP où les coefficients ne sont pas analytiques mais de classe Gevrey.
• Applications physiques : Les résultats sont pertinents pour la modélisation de systèmes quantiques ou d'ondes dans des milieux où les potentiels décroissent de manière sous-exponentielle (par exemple, certains modèles de plasmas ou de milieux désordonnés), offrant des prédictions précises sur le taux de dissipation de l'énergie.

En résumé, cet article fournit une théorie robuste pour l'analyse spectrale et la dynamique des ondes dans des milieux perturbés par des potentiels à décroissance lente, en établissant un pont précis entre la vitesse de décroissance spatiale des coefficients et la vitesse de décroissance temporelle de l'énergie.