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

Aquí presento un resumen técnico detallado del artículo "Estabilidad Lipschitz de dimensión infinita en el problema de Calderón y bases Zernike generales" de Henrik Garde, Markus Hirvensalo y Nuutti Hyvönen.

1. El Problema: Inestabilidad Logarítmica vs. Estabilidad Lipschitz

El problema inverso de conductividad de Calderón busca determinar una conductividad interna $\gamma(x)$ en un dominio $\Omega$ (en este caso, la bola unitaria $B \subset \mathbb{R}^d$) a partir de mediciones en la frontera, representadas por el operador de Dirichlet-a-Neumann $\Lambda(\gamma)$.

• Estado del arte: Para clases generales e infinitas de conductividades, la estabilidad de este problema es logarítmica condicional. Esto significa que pequeños errores en las mediciones de frontera pueden provocar grandes errores en la reconstrucción interna, lo que hace que el problema sea altamente inestable.
• Limitación actual: La estabilidad Lipschitz (donde el error en la reconstrucción es proporcional linealmente al error en la medición, lo cual es deseable para algoritmos de reconstrucción) solo se ha demostrado previamente para clases de conductividades de dimensión finita.
• Objetivo del trabajo: Demostrar la existencia de clases infinitas de conductividades para las cuales se cumple la estabilidad Lipschitz en el dominio completo, no solo en la frontera.

2. Metodología y Enfoque

Los autores desarrollan un marco teórico que combina análisis funcional, teoría de ecuaciones diferenciales parciales y polinomios ortogonales.

A. Reducción a la Frontera

El punto de partida es el teorema clásico de Sylvester-Uhlmann, que establece que la conductividad está determinada Lipschitz-mente en la frontera ($\partial B$) a partir de $\Lambda(\gamma)$. El desafío es extender esta estabilidad al interior del dominio.

• Los autores proponen que si la diferencia entre dos conductividades, $\gamma_1 - \gamma_2$, pertenece a ciertas clases funcionales específicas, su norma en el interior puede acotarse por su restricción en la frontera.
• Utilizan el principio del máximo débil para funciones que satisfacen ciertas ecuaciones elípticas, lo que permite acotar el máximo global por el máximo en la frontera.

B. Construcción de Bases Zernike Generales

Para definir estas clases infinitas, los autores construyen una base ortonormal para $L^2(B)$ en cualquier dimensión $d \geq 2$, generalizando las bases Zernike (usualmente conocidas en 2D para óptica y tomografía).

• Separación de variables: Descomponen la función en coordenadas polares $x = r\theta$.
• Parte angular: Utilizan armónicos esféricos $\{f_{\ell,m}\}$ en la esfera unitaria $\partial B$.
• Parte radial: Resuelven un problema de autovalores de Sturm-Liouville singular para obtener polinomios radiales $R_{\ell,k}(r)$. Estos son eigenfunciones de un operador Zernike generalizado en $d$ dimensiones.
• Definición de la base: Las funciones base son $\psi_{\ell,m,k}(r\theta) = R_{\ell,k}(r)f_{\ell,m}(\theta)$.

C. Condiciones de Decaimiento y Restricción Radial

El núcleo de la contribución teórica reside en imponer condiciones sobre los coeficientes de la expansión de la conductividad en esta base:

1. Espacio ponderado $\ell^1$: Se requiere que los coeficientes de la expansión de la diferencia de conductividades decaigan lo suficientemente rápido (en un espacio $\ell^1$ ponderado por la norma $L^\infty$ de las funciones base cerca de la frontera). Esto asegura que la traza en la frontera esté bien definida dentro de la serie.
2. Restricción de cancelación radial: Se introduce una condición técnica (Ecuación 1.2) que limita la "libertad" en la dirección radial. Específicamente, prohíbe que los coeficientes de diferentes órdenes radiales ($k$) se cancelen excesivamente entre sí para un mismo modo angular. Esto implica que la conductividad no puede variar arbitrariamente en el radio si se desea estabilidad.

3. Contribuciones Clave y Resultados

A. Teorema de Estabilidad Lipschitz Infinita (Teorema 1.6)

Los autores demuestran que si la diferencia de conductividades $\gamma_1 - \gamma_2$ pertenece a un subespacio $A$ definido por la base Zernike y satisface la condición de no-cancelación radial, entonces:
$$\|\gamma_1 - \gamma_2\|_{L^2(B)} \leq C \|\Lambda(\gamma_1) - \Lambda(\gamma_2)\|$$
donde $C$ es una constante independiente de las conductividades. Esto establece estabilidad Lipschitz en dimensión infinita.

B. Construcción y Propiedades de las Bases Zernike $d$-dimensionales

• Derivan explícitamente los polinomios radiales $R_{\ell,k}(r)$ en términos de funciones hipergeométricas de Gauss y coeficientes binomiales.
• Demuestran que estas funciones forman una base ortonormal completa en $L^2(B)$ con peso $r^{d-1}$.
• Establecen cotas uniformes y el comportamiento en la frontera ($r=1$), demostrando que $R_{\ell,k}(1) > 0$ y acotado inferiormente, lo cual es crucial para la constante de estabilidad.

C. Casos Específicos (2D y 3D)

• 2D: Recuperan las bases Zernike clásicas utilizadas en óptica y tomografía computarizada, mostrando cómo se integran en este marco de estabilidad.
• 3D: Proporcionan la formulación explícita para la dimensión 3, relevante para aplicaciones médicas y geofísicas, utilizando armónicos esféricos de Laplace.

D. Corolario de Subespacios Fijos

Si se restringe la conductividad a un subespacio donde solo varía un índice radial fijo $k$ (es decir, $\gamma_1 - \gamma_2 \in A_k$), la condición de estabilidad se simplifica significativamente, garantizando la estabilidad Lipschitz con una constante explícita que depende de $k$ y $d$.

4. Significado e Impacto

1. Ruptura de la barrera de dimensión: Este trabajo es pionero al demostrar que la estabilidad Lipschitz no es exclusiva de problemas de dimensión finita. Abre la puerta a la reconstrucción estable de conductividades con una complejidad mucho mayor (infinita) que las aproximaciones por elementos finitos o parámetros constantes.
2. Fundamento para algoritmos iterativos: La estabilidad Lipschitz es un requisito fundamental para garantizar la convergencia de métodos de reconstrucción iterativos (como Newton o métodos de gradiente). Sin ella, estos métodos pueden divergir o estancarse.
3. Compromiso entre libertad y estabilidad: El artículo ilustra un compromiso matemático preciso: se puede tener un espacio de búsqueda infinito, pero la estabilidad Lipschitz requiere restricciones en la "suavidad" o el comportamiento de los coeficientes radiales (decaimiento rápido y falta de cancelación destructiva).
4. Generalización dimensional: Al proporcionar una construcción válida para cualquier $d \geq 2$, el trabajo unifica el tratamiento del problema de Calderón en diferentes dimensiones espaciales, algo que no estaba claro en la literatura previa (especialmente para $d > 2$).

En resumen, el artículo proporciona una prueba elemental pero rigurosa de que, bajo ciertas condiciones estructurales en la expansión de la conductividad (específicamente usando bases Zernike generalizadas), el problema inverso de Calderón deja de ser inestable y se vuelve Lipschitz-estable, incluso en espacios de dimensión infinita.