The initial-to-final-state inverse problem with critically-singular potentials

1. Problématique et Contexte

L'article s'intéresse à un problème inverse pour l'équation de Schrödinger en dimensions élevées ($n \ge 2$). Le cadre est le suivant :

• Équation : On considère l'équation de Schrödinger avec un potentiel indépendant du temps $V(x)$ :
  $$i\partial_t u = -\Delta u + V u, \quad u(0, x) = f(x)$$
  où $u \in C([0, T]; L^2(\mathbb{R}^n))$.
• Données : On dispose de l'application "état initial vers état final" (initial-to-final-state map), notée $U_T$, qui envoie l'état initial $f \in L^2(\mathbb{R}^n)$ à l'instant $t=0$ vers la solution $u(T, \cdot)$ à un temps fixe $T$.
• Question centrale : L'application $U_T$ détermine-t-elle de manière unique le potentiel $V$ (et donc l'hamiltonien $-\Delta + V$) ?

Contexte antérieur :
Des résultats de unicité existaient déjà pour des potentiels dépendant du temps, mais sous des hypothèses de décroissance très fortes (décroissance super-exponentielle à l'infini). Pour les potentiels indépendants du temps, des travaux précédents (notamment de Caro et Ruiz) avaient établi l'unicité pour des potentiels bornés avec une décroissance polynomiale.

Objectif de ce papier :
Affaiblir considérablement les hypothèses sur le potentiel $V$. Les auteurs visent à prouver l'unicité pour des potentiels indépendants du temps qui peuvent présenter :

1. Des singularités de type $L^q$ sur des ensembles de mesure finie.
2. Une intégrabilité $L^1$ (décroissance $L^1$) à l'infini, sans nécessiter de décroissance super-exponentielle ni même polynomiale forte.

2. Résultats Principaux

Le résultat principal est énoncé dans le Théorème 1.1 :

Soient $V_1, V_2$ deux potentiels indépendants du temps appartenant à l'espace $L^1(\mathbb{R}^n) \cap L^q(\mathbb{R}^n)$, avec :

• $q > 1$ si $n = 2$,
• $q \ge n/2$ si $n \ge 3$.

Si les applications état initial-état final coïncident ($U_T^1 = U_T^2$), alors les potentiels sont égaux presque partout :
$$V_1 = V_2 \quad \text{dans } \mathbb{R}^n.$$

Points clés de l'amélioration :

• Le cas critique $q = n/2$ (pour $n \ge 3$) est inclus. C'est un cas difficile car les estimées de résolvant classiques ne fournissent pas de décroissance en la fréquence d'énergie $\lambda$.
• La condition de décroissance à l'infini est réduite à une simple intégrabilité $L^1$, ce qui est une amélioration majeure par rapport aux hypothèses de décroissance super-linéaire ou super-exponentielle requises dans les cas dépendants du temps.

3. Méthodologie et Outils Techniques

La preuve repose sur une stratégie en plusieurs étapes, combinant l'analyse harmonique, la théorie des équations aux dérivées partielles et l'analyse fonctionnelle.

A. Relation d'orthogonalité de type Alessandrini

L'approche standard pour les problèmes inverses consiste à utiliser une relation d'orthogonalité. Si $U_T^1 = U_T^2$, alors pour toute paire de solutions $u_1, v_2$ associées respectivement aux potentiels $V_1$ et $V_2$, on a :
$$\int_{\Sigma} (V_1 - V_2) u_1 v_2 \, dx dt = 0$$
où $\Sigma = (0, T) \times \mathbb{R}^n$.
Le défi majeur est que cette relation est initialement valable pour des solutions "physiques" (dans $C([0, T]; L^2)$), mais pour reconstruire $V$, il faut tester cette relation avec des solutions oscillantes spécifiques (états stationnaires) qui ne sont pas dans $L^2$ sur l'espace entier. Les auteurs doivent donc étendre cette relation d'orthogonalité à ces solutions non-physiques.

B. Construction d'États Stationnaires (Stationary States)

Les auteurs construisent des solutions spéciales de la forme :
$$\psi_{\lambda, \omega}(t, x) = e^{-i\lambda^2 t} \left( e^{-i\lambda \omega \cdot x} + w_{cor}(x) \right)$$
où $e^{-i\lambda \omega \cdot x}$ est une solution de l'équation libre et $w_{cor}$ est un terme de correction dû au potentiel $V$.
La construction de $w_{cor}$ repose sur l'inversion d'un opérateur de la forme $(I - P_\lambda \circ V)$, où $P_\lambda$ est l'opérateur de résolvant de l'équation de Helmholtz $(\Delta + \lambda^2)u = f$.

C. Gestion du Cas Critique et Estimations de Résolvant

C'est ici que réside l'innovation technique principale du papier :

1. Cas non-critique ($q > n/2$) : L'estimation de résolvant classique de Kenig-Ruiz-Sogge suffit. Elle garantit que l'opérateur $P_\lambda \circ V$ devient petit lorsque la fréquence $\lambda \to \infty$, permettant l'inversion par série de Neumann.
2. Cas critique ($q = n/2$) : L'estimation classique de Kenig-Ruiz-Sogge ne fournit pas de décroissance en $\lambda$ pour l'opérateur $P_\lambda \circ V$ dans les espaces standards.
   • Solution : Les auteurs introduisent une nouvelle échelle d'espaces de Banach ($X_\lambda$ et son dual $X^*_\lambda$) qui combine des espaces $L^p$ avec des poids dépendant de $\lambda$.
   • Ils établissent une estimation de résolvant raffinée (une amélioration de l'estimation de Kenig-Ruiz-Sogge) dans ces espaces. Cela permet de récupérer un facteur de décroissance en $\lambda$ même au point critique, rendant l'inversion de l'opérateur possible.

D. Extension de la Relation d'Orthogonalité

Pour utiliser les états stationnaires (qui ne sont pas dans $L^2$) dans la relation d'orthogonalité, les auteurs développent des identités d'intégration par parties (Lemmes 4.1 à 4.3). Ils utilisent des approximations par des fonctions à support compact et des régularisations pour justifier le passage à la limite, en exploitant la structure oscillatoire des solutions (phases imaginaires pures) qui permet d'éviter les hypothèses de décroissance forte requises dans le cas dépendant du temps (où l'on utilise des solutions d'optique géométrique complexe à croissance exponentielle).

4. Preuve de l'Unicité

Une fois la relation d'orthogonalité établie pour les états stationnaires $\psi_1, \psi_2$ :
$$\int_{\Sigma} (V_1 - V_2) \psi_1 \psi_2 = 0$$
Les auteurs choisissent des vecteurs de fréquence $\omega_1, \omega_2$ dépendant d'un vecteur $\xi$ arbitraire et d'un paramètre $\lambda$ grand. En développant les termes, l'intégrale se sépare en :

1. Un terme principal contenant la transformée de Fourier de la différence des potentiels $F = V_1 - V_2$ : $\hat{F}(\xi)$.
2. Des termes d'erreur impliquant les corrections $w_{cor}$.

Grâce aux estimées de décroissance des termes de correction (obtenues via les espaces $X^*_\lambda$ dans le cas critique et les espaces $L^p$ standards dans le cas non-critique), les termes d'erreur tendent vers zéro lorsque $\lambda \to \infty$.
Il en résulte que $\hat{F}(\xi) = 0$ pour tout $\xi$, ce qui implique $V_1 = V_2$ presque partout.

5. Signification et Impact

• Avancée théorique : Ce travail repousse les limites de la théorie des problèmes inverses pour l'équation de Schrödinger en éliminant le besoin de décroissance super-exponentielle pour les potentiels statiques. Il montre que la structure temporelle (indépendance du temps) est suffisante pour compenser le manque de décroissance spatiale forte.
• Outils nouveaux : L'introduction d'espaces de Banach adaptés pour traiter le cas critique $q=n/2$ et la refinement des estimées de résolvant de Kenig-Ruiz-Sogge constituent des outils puissants qui pourraient être appliqués à d'autres problèmes de diffusion ou d'ondes avec des singularités critiques.
• Contraste avec le cas dépendant du temps : L'article met en lumière une différence fondamentale : dans le cas dépendant du temps, les solutions d'optique géométrique complexe croissent dans certaines directions, imposant une décroissance très forte du potentiel pour contrôler les erreurs. Dans le cas statique, les solutions sont purement oscillatoires (unimodulaires), permettant une analyse plus fine avec des potentiels beaucoup moins réguliers à l'infini.

En résumé, ce papier résout un problème d'unicité ouvert en affaiblissant drastiquement les hypothèses de régularité et de décroissance des potentiels, en s'appuyant sur une analyse fine des opérateurs de résolvant et une ingénierie subtile des espaces fonctionnels.