On the Calderón problem with piecewise polynomial anisotropic conductivities and many-flat-face interfaces

Resumen Técnico: Unicidad para el Problema de Calderón con Conductividades Anisotrópicas Polinomiales por Partes

Enunciado del Problema
El artículo aborda el problema de la conductividad inversa (problema de Calderón) en dimensiones $n \ge 3$. El objetivo es determinar un campo de matrices simétricas y uniformemente elípticas $\sigma(x)$ que representa una conductividad anisotrópica dentro de un dominio acotado $\Omega \subset \mathbb{R}^n$. La determinación se basa en un mapa de Dirichlet-a-Neumann (DN) local, $\Lambda_{\Omega, \Sigma}$, definido en un subconjunto relativamente abierto $\Sigma \subset \partial \Omega$.

En el caso anisotrópico general, se conjetura que la unicidad solo se cumple módulo de las difeomorfismos que fijan la frontera (libertad de gauge), ya que diferentes conductividades relacionadas por tales transformaciones producen mediciones de frontera idénticas. Para superar este gauge, el artículo impone supuestos estructurales específicos:

1. Partición Conocida: El dominio se subdivide en una colección finita conocida de celdas Lipschitz disjuntas $\{D_\alpha\}$.
2. Estructura Polinomial por Partes: En cada celda $D_\alpha$, la conductividad $\sigma$ es la restricción de una matriz de valores polinomiales simétricos $P_\alpha$ de grado conocido.
3. Discontinuidad: La conductividad puede ser discontinua a través de las interfaces entre las celdas.

El problema inverso consiste en recuperar unívocamente los coeficientes de las matrices polinomiales $\{P_\alpha\}$ a partir del mapa DN local en $\Sigma$.

Metodología
La prueba de unicidad se basa en un argumento de "despojo de capas" (layer-stripping) combinado con un teorema de recuperación local para una sola celda. La metodología procede en dos etapas principales:

1. Recuperación de Celda Local (Determinación de Caras Planas):
   El artículo establece que si la conductividad en una sola celda es un polinomio de grado como máximo $N$, y si el mapa DN es conocido en suficientes caras planas de esa celda, el polinomio queda determinado de forma unívoca.

   • Formas de Schur: Los datos de frontera determinan la "forma de Schur conormal" de la conductividad en cada interfaz plana. Esta forma es una cantidad invariante de gauge derivada del símbolo principal del mapa DN.
   • Propagación Polinomial: Dado que la conductividad es polinomial, su forma de Schur en una cara plana también es polinomial. El conocimiento de la forma de Schur en un pequeño parche abierto de una cara se extiende a todo el hiperplano de soporte mediante el teorema de identidad para polinomios.
   • Recuperación Puntual: En los puntos donde tres hiperplanos de soporte se intersectan con conormales linealmente independientes, las tres formas de Schur correspondientes determinan unívocamente la matriz de conductividad en ese punto.
   • Argumento Algebraico: Al requerir que las formas de Schur coincidan en un número suficiente de caras planas genéricas (específicamente $N+3$ caras que satisfacen una condición de genericidad respecto a sus intersecciones), el artículo utiliza un argumento de divisibilidad para demostrar que la diferencia entre dos candidatas de conductividad polinomial debe anularse identicamente.

2. Recuperación Global (Despojo de Capas y Extensión Interior):
   La unicidad global se logra recuperando las celdas secuencialmente de acuerdo con un "orden de recuperación de caras planas admisible".

   • Admisibilidad: El orden requiere que, en cada paso, la siguiente celda a recuperar tenga acceso a suficientes caras planas genéricas (ya sean interfaces exteriores o internas con celdas previamente recuperadas) y que la región recuperada permanezca conectada por caras a la frontera inicial medida.
   • Extensión Interior: Una vez recuperada una región de celdas, esta actúa como una "región conocida". El artículo utiliza un argumento de aproximación de Runge (basado en la continuación única y el teorema de Holmgren) para mostrar que el mapa DN local en la frontera original $\Sigma$ determina el mapa DN local en las interfaces internas de la región no recuperada.
   • Propagación: Esto permite aplicar el teorema de recuperación local de forma iterativa a la siguiente celda en la secuencia, propagando la unicidad desde la frontera hacia el interior.

Resultados Clave

• Teorema 1.4 (Unicidad): Bajo los supuestos de una subdivisión Lipschitz plana finita conocida, conductividades polinomiales por partes y la existencia de un orden de recuperación de caras planas admisible, el mapa DN local en $\Sigma$ determina unívocamente la conductividad $\sigma$ en $\Omega$. Específicamente, si $\Lambda_{\Omega, \Sigma}^{\sigma(1)} = \Lambda_{\Omega, \Sigma}^{\sigma(2)}$, entonces $\sigma(1) = \sigma(2)$ en todo el dominio.
• Corolario 1.5 (Estabilidad): Debido a la naturaleza de dimensión finita del espacio de parámetros (los coeficientes de los polinomios), el artículo deriva un estimador de estabilidad de Hölder. En conjuntos compactos de parámetros admisibles, la diferencia en los coeficientes polinomiales está acotada por una potencia de Hölder de la diferencia en los mapas DN.

Significancia y Reivindicaciones
El artículo afirma proporcionar un resultado de unicidad riguroso para un problema de Calderón anisotrópico de dimensión finita que incluye conductividades discontinuas. Su significancia radica en:

• Superar el Gauge: Al combinar supuestos estructurales (polinomios por partes) con condiciones geométricas (muchas caras planas), el trabajo elimina el gauge de difeomorfismo que típicamente impide la unicidad en problemas anisotrópicos.
• Generalización de Trabajos Previos: Extiende los resultados existentes para conductividades constantes por partes a la clase más general de conductividades polinomiales por partes.
• Hipótesis Geométrica: El artículo introduce una condición geométrica específica ("orden de recuperación de caras planas admisible") que involucra intersecciones genéricas de caras planas. Esta condición es análoga a las hipótesis de fronteras curvas en otros trabajos, pero está adaptada a modelos poliédricos o de tipo malla donde las interfaces planas son intrínsecas.
• Estabilidad: La derivación de un estimador de estabilidad de Hölder es una consecuencia directa del entorno de dimensión finita y la analiticidad del mapa directo, proporcionando una medida cuantitativa de la robustez de la recuperación.

La estrategia de la prueba separa el problema algebraico de recuperar un polinomio a partir de sus formas de Schur en hiperplanos del problema analítico de propagar los datos de frontera a través de regiones conocidas, utilizando la aproximación de Runge y las propiedades de continuación única de las ecuaciones elípticas con coeficientes analíticos.