Infinite-dimensional Lipschitz stability in the Calderón problem and general Zernike bases

1. Il Problema: Stabilità nel Problema di Calderón

Il problema di Calderón (o problema inverso di conduttività) consiste nel determinare il coefficiente di conduttività elettrica $\gamma(x)$ all'interno di un dominio $\Omega$ (in questo caso la palla unitaria $B \subset \mathbb{R}^d$ con $d \ge 2$) misurando solo dati al bordo.
I dati sono forniti dalla mappa Dirichlet-to-Neumann (DtN), $\Lambda(\gamma)$, che associa un potenziale elettrico al bordo ($u|_{\partial B}$) alla corrente normale uscente ($\nu \cdot \gamma \nabla u|_{\partial B}$).

Stato dell'arte e sfida:

• Per classi generali di conduttività (infinite-dimensionali), la stabilità è nota essere logaritmica (condizionata), il che implica che piccole errori nei dati di bordo possono causare grandi errori nella ricostruzione interna.
• La stabilità Lipschitz (dove l'errore nella soluzione è limitato linearmente dall'errore nei dati) è stata finora dimostrata solo per classi di conduttività finite-dimensionali.
• L'obiettivo del paper è dimostrare l'esistenza di classi infinite-dimensionali di conduttività per le quali vale la stabilità Lipschitz, superando il limite logaritmico.

2. Metodologia e Strumenti Matematici

Gli autori adottano un approccio basato sull'analisi funzionale e sulla teoria degli operatori, strutturato in tre fasi principali:

A. Riduzione alla Determinazione al Bordo

Il punto di partenza è un risultato classico di Sylvester e Uhlmann (Teorema 1.1), che garantisce la stabilità Lipschitz per la determinazione del coefficiente sulla frontiera ($\partial B$).
Il paper dimostra che se la differenza tra due conduttività, $\gamma_1 - \gamma_2$, può essere controllata dalla sua restrizione al bordo (cioè se la norma $L^\infty$ interna è limitata dalla norma $L^\infty$ al bordo), allora si ottiene la stabilità Lipschitz globale.
Ciò è possibile se $\gamma_1 - \gamma_2$ soddisfa un principio del massimo debole (es. è soluzione di un'equazione ellittica o è armonica).

B. Costruzione di Basi Zernike Generalizzate

Per costruire classi infinite-dimensionali con questa proprietà, gli autori introducono una nuova famiglia di funzioni di base per $L^2(B)$, generalizzando le polinomi di Zernike (tradizionalmente usati in 2D per l'ottica e la tomografia) a qualsiasi dimensione spaziale $d \ge 2$.
La costruzione avviene tramite:

1. Armoniche Sferiche: Una base ortonormale $\{f_{\ell,m}\}$ per $L^2(\partial B)$.
2. Polinomi Radiali: Risoluzione di un problema agli autovalori di Sturm-Liouville singolare per ottenere una base ortonormale $\{R_{\ell,k}\}$ nello spazio pesato $L^2_{r^{d-1}}([0,1])$.
3. Combinazione: Le funzioni di base sono $\psi_{\ell,m,k}(r\theta) = R_{\ell,k}(r)f_{\ell,m}(\theta)$.

C. Condizioni di Decadimento e Vincoli

Il cuore della dimostrazione risiede nell'imporre condizioni specifiche sui coefficienti dell'espansione di $\gamma_1 - \gamma_2$ nella base Zernike.
Affinché la traccia al bordo sia ben definita e la stabilità Lipschitz valga, i coefficienti $(c_{n,k})$ devono appartenere a uno spazio pesato $\ell^1$ (decadimento sufficientemente rapido) e soddisfare una condizione di "non cancellazione" nelle direzioni radiali (Equazione 1.2).
Essenzialmente, la libertà nella direzione radiale è limitata: non è permesso che i termini si cancellino eccessivamente vicino al bordo, il che garantirebbe che la norma interna sia controllata dal valore al bordo.

3. Risultati Chiave

1. Teorema di Stabilità Lipschitz Infinito-Dimensionale (Teorema 1.6):
   Gli autori dimostrano che esiste un'infinità di classi infinite-dimensionali di conduttività per le quali vale la stabilità Lipschitz. In particolare, se la differenza $\gamma_1 - \gamma_2$ appartiene a uno spazio $A$ definito da una base $\{\psi_{n,k}\}$ e i coefficienti soddisfano una specifica disuguaglianza di decadimento, allora:
   $$\|\gamma_1 - \gamma_2\|_{L^2(B)} \le C \|\Lambda(\gamma_1) - \Lambda(\gamma_2)\|$$
   dove $C$ è una costante indipendente dalla dimensione della classe.

2. Costruzione delle Basi Zernike in $d$ dimensioni (Sezione 2):
   Viene fornita una costruzione esplicita e rigorosa delle basi Zernike per $d \ge 2$. Si dimostra che queste basi soddisfano tutte le ipotesi tecniche necessarie (Assunzione 1.4), inclusa la limitatezza delle funzioni vicino al bordo e il decadimento dei coefficienti.

   • Si forniscono formule esplicite per $d=2$ e $d=3$, rilevanti per applicazioni pratiche.
   • Si dimostra che i polinomi radiali $R_{\ell,k}$ sono polinomi di grado $\ell+2k$.

3. Corollari Pratici (Corollario 1.7):
   Viene mostrato che la stabilità vale anche per sottospazi infiniti dove si fissa l'indice radiale $k$ (variazioni con un profilo radiale fisso ma con tutte le frequenze angolari libere), purché i coefficienti decadano abbastanza velocemente.

4. Contributi e Significatività

• Superamento del Limite Logaritmico: Il lavoro rompe il paradigma secondo cui la stabilità Lipschitz è possibile solo per modelli parametrici a parametri finiti. Dimostra che la stabilità può essere mantenuta anche in spazi di dimensione infinita, a patto di imporre vincoli strutturali (decadimento dei coefficienti) che riflettano la fisica del problema (comportamento al bordo).
• Generalizzazione delle Basi Zernike: L'estensione delle basi di Zernike a dimensioni arbitrarie $d \ge 2$ è un contributo matematico significativo, fornendo strumenti analitici per problemi inversi in 3D e oltre, non solo in 2D.
• Implicazioni per la Ricostruzione Numerica: La stabilità Lipschitz è cruciale per la convergenza dei metodi iterativi di ricostruzione (come il metodo di Newton o algoritmi di discesa del gradiente). Questo risultato suggerisce che algoritmi di ricostruzione basati su espansioni in basi Zernike (con coefficienti opportunamente regolarizzati) potrebbero convergere in modo stabile e veloce, anche per problemi non lineari, senza bisogno di discretizzazioni troppo grossolane.
• Applicabilità: I risultati sono rilevanti per la Tomografia a Impedenza Elettrica (EIT) e altri problemi di imaging medico o industriale dove la conduttività interna deve essere stimata da misure di superficie.

In sintesi, il paper fornisce una prova elementare ma potente che la "maledizione" della stabilità logaritmica nel problema di Calderón può essere evitata per classi di funzioni realistiche e infinite-dimensionali, aprendo la strada a metodi di inversione più robusti e teoricamente fondati.