Eigenstructure of the linearized electrical impedance tomography problem under radial perturbations

1. Il Problema

Il paper si concentra sul problema inverso della Tomografia ad Impedenza Elettrica (EIT), specificamente sull'analisi della derivata di Fréchet dell'applicazione di Dirichlet-Neumann (ND) linearizzata.

• Contesto Fisico: Si considera un mezzo conduttivo descritto dall'equazione di conducibilità $-\nabla \cdot (\gamma \nabla u) = 0$ all'interno di una palla unitaria $B \subset \mathbb{R}^d$ ($d \ge 2$).
• Mappatura: L'obiettivo è studiare la mappa non lineare $\gamma \mapsto \Lambda(\gamma)$, che associa una distribuzione di conducibilità $\gamma$ alla mappa ND (che lega la densità di corrente di Neumann al potenziale di Dirichlet sul bordo).
• Linearizzazione: L'analisi si focalizza sulla derivata di Fréchet $F = D\Lambda(1)$ valutata in $\gamma=1$ (conduttività omogenea) e applicata a una perturbazione $\eta$. La relazione fondamentale è data da:
  $$\langle (F\eta)f, g \rangle_{L^2(\partial B)} = -\int_B \eta \nabla u_f \cdot \nabla u_g \, dx$$
  dove $u_f$ e $u_g$ sono potenziali armonici con dati al bordo di Neumann $f$ e $g$.
• Sfida: Sebbene la continuità di $F: L^\infty(B) \to \mathcal{L}(L^2_\diamond(\partial B))$ sia nota, le proprietà di approssimabilità tramite operatori di rango finito e la compattezza in topologie più forti (come $L^2(B)$) non erano state pienamente caratterizzate, specialmente in dimensioni $d > 2$.

2. Metodologia

L'autore adotta un approccio basato sulla simmetria radiale e sull'analisi spettrale:

• Geometria e Coordinate: Il dominio è la palla unitaria. Si utilizzano coordinate polari $(r, \theta)$ per sfruttare la simmetria.
• Sottospazio delle Perturbazioni: L'analisi è limitata al sottospazio $\mathcal{R} \subset L^2(B)$ delle perturbazioni radialmente simmetriche (indipendenti dall'angolo $\theta$). Questo permette di separare le variabili.
• Basi Ortogonali:
  • Angolare: Si utilizzano gli armonici sferici $\{f_{\ell,m}\}$ come base ortonormale per lo spazio $L^2(\partial B)$.
  • Radiale: Si utilizzano i polinomi di Jacobi (shiftati e normalizzati) $\{P_k(r)\}$ come base per lo spazio radiale pesato $L^2_{r^{d-1}}((0,1))$.
• Strumenti Analitici:
  • Sviluppo in serie delle perturbazioni $\eta$ nella base di Jacobi.
  • Calcolo esplicito degli elementi di matrice dell'operatore $F\eta$ nella base degli armonici sferici.
  • Uso di identità combinatorie e funzioni Beta per valutare somme di fattoriali e coefficienti.
  • Applicazione della disuguaglianza di Cauchy-Schwarz per stimare la decadenza degli autovalori.

3. Contributi Chiave e Risultati Principali

A. Struttura degli Autovalori (Teorema 2.1)

Il risultato centrale è la caratterizzazione completa della struttura spettrale dell'operatore linearizzato $F\eta$ quando $\eta \in \mathcal{R}$.

• Diagonalizzazione: L'operatore $F\eta$ è diagonale nella base degli armonici sferici. Gli armonici sferici $f_{\ell,m}$ sono autovettori di $F\eta$.
• Formula Esplicita: Gli autovalori $\lambda_\ell(\eta)$ sono indipendenti dall'ordine $m$ e dipendono solo dal grado $\ell$ e dai coefficienti di Fourier-Jacobi della perturbazione $\eta$. La formula è:
  $$\lambda_\ell(\eta) = \sum_{k=0}^{2\ell-2} (-1)^{k+1} a_k(\eta) \sqrt{2k+d} \frac{1}{\ell} \frac{(2\ell-2+d)! (2\ell-2)!}{(2\ell-2+d+k)! (2\ell-2-k)!}$$
  dove $a_k(\eta)$ sono i coefficienti di espansione di $\eta$ in polinomi di Jacobi.
• Decadimento Uniforme: Viene dimostrato che gli autovalori decadono uniformemente rispetto al grado $\ell$. In particolare:
  $$|\lambda_\ell(\eta)| \le C_d \|\eta\|_{L^2(B)} \ell^{-1/2}$$
  Questo decadimento è cruciale per le proprietà di compattezza.

B. Compattezza e Approssimabilità (Corollario 2.2)

• Estensione a $L^2(B)$: Il lavoro estende la definizione della derivata di Fréchet al sottospazio infinito-dimensionale $\mathcal{R} \subset L^2(B)$, andando oltre lo spazio standard $L^\infty(B)$.
• Compattezza: Grazie al decadimento degli autovalori ($\ell^{-1/2}$), l'operatore $F$, ristretto a $\mathcal{R}$, è un operatore compatto.
• Approssimazione di Rango Finito: L'operatore può essere approssimato nella norma operatoriale da una sequenza di operatori di rango finito. Questo è dimostrato costruendo l'approssimazione tronca $F_N$ (somma dei primi $N$ termini della decomposizione spettrale) e mostrando che $\|F - F_N\|_{op} \to 0$ quando $N \to \infty$.

4. Significato e Implicazioni

1. Avanzamento Teorico: Questo è il primo risultato che dimostra la compattezza e l'approssimabilità di rango finito della derivata di Fréchet dell'EIT linearizzato in dimensioni $d \ge 2$ per perturbazioni in $L^2(B)$ (invece che solo in $L^\infty$).
2. Implicazioni Numeriche:
   • La possibilità di approssimare $F$ con operatori di rango finito è fondamentale per gli algoritmi numerici (es. metodi di regolarizzazione, algoritmi di Newton).
   • Conferma che l'informazione contenuta nelle perturbazioni radiali è "compatta", il che giustifica l'uso di tecniche di compressione o troncamento spettrale nelle ricostruzioni EIT.
3. Generalizzazione: Mentre studi precedenti (es. Garde e Hyvönen) avevano trattato il caso bidimensionale, questo lavoro generalizza i risultati a qualsiasi dimensione spaziale $d \ge 2$.
4. Struttura Matematica: La scoperta che le armoniche sferiche sono autovettori per perturbazioni radiali semplifica enormemente l'analisi del problema inverso in geometrie sferiche, permettendo di trattare il problema come una serie di problemi scalari indipendenti per ogni grado $\ell$.

In sintesi, il paper fornisce un quadro analitico rigoroso per la linearizzazione dell'EIT in domini sferici, stabilendo proprietà di compattezza essenziali per la stabilità e l'efficienza dei metodi numerici di ricostruzione delle immagini.