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

1. Il Problema

Il lavoro affronta un problema inverso per l'equazione di Schrödinger in dimensioni elevate ($n \ge 2$). L'obiettivo è determinare se la mappa che associa uno stato iniziale $f \in L^2(\mathbb{R}^n)$ (al tempo $t=0$) allo stato finale $u(T, \cdot)$ (al tempo fisso $t=T$) determina univocamente il potenziale $V(x)$ che genera l'evoluzione temporale.

L'equazione considerata è:
$$\begin{cases} i\partial_t u = -\Delta u + V(x)u & \text{in } (0, T) \times \mathbb{R}^n, \\ u(0, x) = f(x) & \text{in } \mathbb{R}^n. \end{cases}$$
La mappa iniziale-finale è definita come $U_T: f \mapsto u(T, \cdot)$. La domanda centrale è: se $U_T^1 = U_T^2$ per due potenziali $V_1$ e $V_2$, ne consegue che $V_1 = V_2$ quasi ovunque?

2. Contesto e Stato dell'Arte

Prima di questo lavoro, la letteratura sui potenziali dipendenti dal tempo richiedeva condizioni di decadimento molto forti (decadimento super-esponenziale all'infinito) per garantire l'unicità, anche per potenziali limitati o con singolarità locali.
Per i potenziali indipendenti dal tempo, lavori precedenti (es. [10]) avevano ottenuto risultati di unicità sotto assunzioni di decadimento super-lineare.
Il contributo di questo paper è rilassare ulteriormente queste ipotesi, permettendo potenziali che possono presentare:

1. Singolarità critiche: Appartenenti a $L^q(\mathbb{R}^n)$ con $q = n/2$ (per $n \ge 3$) o $q > 1$ (per $n=2$).
2. Decadimento debole: Appartenenti a $L^1(\mathbb{R}^n)$ all'infinito, senza richiedere decadimento super-esponenziale.

3. Metodologia e Strumenti Tecnici

La prova si basa su una strategia di tipo "Calderón", adattata all'equazione di Schrödinger, che coinvolge i seguenti passaggi chiave:

A. Relazione di Ortogonalità di Alessandrini

Si parte dall'assunzione $U_T^1 = U_T^2$. Questo implica una relazione di ortogonalità (di tipo Alessandrini) per le soluzioni fisiche $u_1, v_2$ associate ai potenziali $V_1, V_2$:
$$\int_{\Sigma} (V_1 - V_2) u_1 v_2 \, dx dt = 0.$$
Tuttavia, le soluzioni fisiche appartengono a $C([0, T]; L^2)$ e non sono sufficienti per recuperare informazioni puntuali sul potenziale. È necessario estendere questa relazione a soluzioni speciali chiamate stati stazionari.

B. Costruzione di Stati Stazionari

Gli stati stazionari sono soluzioni dell'equazione di Schrödinger della forma:
$$\psi(t, x) = e^{-i\lambda^2 t} (e^{-i\lambda \omega \cdot x} + w_{cor}(x)),$$
dove $e^{-i\lambda \omega \cdot x}$ è un'onda piana libera e $w_{cor}$ è un termine di correzione dovuto al potenziale.
La costruzione di $w_{cor}$ richiede l'inversione di un operatore del tipo $(I - P_\lambda \circ V)$, dove $P_\lambda$ è l'operatore risolutore dell'equazione di Helmholtz $(\Delta + \lambda^2)u = f$.

C. Stime del Risolutore e Casi Critici

La difficoltà principale risiede nel controllo dell'operatore di correzione per grandi valori dell'energia $\lambda$.

• Caso non-endpoint ($q > n/2$): Si utilizzano le stime classiche di Kenig-Ruiz-Sogge per il risolutore di Helmholtz, che garantiscono un decadimento in $\lambda$, permettendo di invertire l'operatore tramite serie di Neumann.
• Caso endpoint critico ($q = n/2$): Le stime classiche non forniscono decadimento in $\lambda$. Per superare questo ostacolo, gli autori introducono una nuova scala di spazi di Banach ($X_\lambda$ e $X^*_\lambda$) che combina stime $L^p$ e $L^q$. All'interno di questo framework, dimostrano una versione raffinata delle stime di Kenig-Ruiz-Sogge, che recupera un fattore di decadimento anche nel caso critico. Questo è il contributo analitico principale del lavoro.

D. Estensione della Relazione di Ortogonalità

Un altro passo cruciale è estendere la relazione di ortogonalità (2) dalle soluzioni fisiche ($L^2$) agli stati stazionari, che non appartengono a $L^2(\mathbb{R}^n)$ a causa della parte oscillante. Gli autori utilizzano tecniche di integrazione per parti e approssimazione con funzioni di taglio per giustificare rigorosamente l'uso di queste soluzioni come funzioni di prova.

4. Risultati Principali

Il Teorema 1.1 stabilisce l'unicità sotto le seguenti condizioni:
Siano $V_1, V_2 \in L^1(\mathbb{R}^n) \cap L^q(\mathbb{R}^n)$ potenziali indipendenti dal tempo, dove:

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

Se le mappe iniziale-finale coincidono ($U_T^1 = U_T^2$), allora $V_1 = V_2$ quasi ovunque.

La prova procede in due fasi:

1. Caso non-endpoint ($q > n/2$): L'uso delle stime standard permette di mostrare che la trasformata di Fourier della differenza $F = V_1 - V_2$ tende a zero per $\lambda \to \infty$.
2. Caso endpoint ($q = n/2$): L'uso degli spazi $X_\lambda$ e delle stime raffinate permette di ottenere lo stesso risultato di decadimento, dimostrando che $\hat{F}(\xi) = 0$ per ogni $\xi$, e quindi $F=0$.

5. Significato e Contributi

Questo lavoro rappresenta un significativo avanzamento nella teoria dei problemi inversi per l'equazione di Schrödinger:

• Rilassamento delle condizioni di decadimento: A differenza dei risultati precedenti per potenziali dipendenti dal tempo che richiedevano decadimento super-esponenziale, questo studio dimostra che per potenziali indipendenti dal tempo è sufficiente un'integrabilità $L^1$ all'infinito.
• Gestione delle singolarità critiche: Il paper risolve il problema per potenziali nella classe critica $L^{n/2}$, un caso che le stime standard non riescono a gestire senza tecniche avanzate.
• Innovazione Analitica: Lo sviluppo di una nuova scala di spazi di Banach e il raffinamento delle stime del risolutore di Helmholtz offrono strumenti potenti che potrebbero essere applicati ad altri problemi di scattering e inversi con potenziali singolari.
• Semplicità del setting temporale: Il lavoro evidenzia come l'indipendenza temporale del potenziale permetta di evitare l'uso di soluzioni ottiche geometriche complesse (che crescono esponenzialmente), sostituendole con stati stazionari puramente oscillanti, semplificando così le ipotesi di regolarità richieste.

In sintesi, il paper chiude un gap importante nella teoria dell'unicità per potenziali singolari, dimostrando che la conoscenza dello stato finale è sufficiente a identificare un potenziale anche in presenza di singolarità forti e decadimento lento, purché il potenziale sia indipendente dal tempo.