Stable determination of the first order perturbation of the biharmonic operator from partial data

Resumen Técnico: Determinación Estable de la Perturbación de Primer Orden del Operador Biharmónico a partir de Datos Parciales

Planteamiento del Problema
Este artículo aborda un problema de valor en la frontera inverso para el operador biarmónico perturbado en un dominio acotado $\Omega \subseteq \mathbb{R}^n$ ($n \ge 3$) con una frontera $\partial\Omega$ suave. El operador se define como $L_{A,q} = \Delta^2 + A(x) \cdot D + q(x)$, donde $D = i^{-1}\nabla$, $A$ es un campo vectorial en $W^{1,\infty}(\Omega, \mathbb{C}^n)$, y $q$ es un potencial en $L^\infty(\Omega, \mathbb{C})$. El estudio se centra en el escenario de datos parciales: la recuperación de los coeficientes $A$ y $q$ a partir de mediciones tomadas únicamente en subconjuntos abiertos arbitrariamente pequeños $\gamma_1, \gamma_2 \subset \partial\Omega$.

El problema asume que los coeficientes son conocidos a priori en una vecindad $\omega_0$ de la frontera. Los datos son proporcionados por un mapa de Dirichlet-a-Neumann (DtN) parcial $\Lambda_{\gamma_1, \gamma_2}^{A,q}$, que mapea los datos de Dirichlet $(f, g)$ soportados en $\gamma_1$ hacia las derivadas normales $(\partial_\nu u, \partial_\nu(\Delta u))$ medidas en $\gamma_2$. El objetivo principal es establecer estimaciones de estabilidad (cuantificar cómo los errores en las mediciones afectan la reconstrucción de los coeficientes) tanto para la perturbación de primer orden $A$ como para la perturbación de orden cero $q$.

Metodología
Los autores emplean el método de soluciones de Óptica Geométrica Compleja (CGO, por sus siglas en inglés), una técnica estándar para problemas inversos de operadores elípticos. La estrategia de prueba involucra varios componentes técnicos clave:

1. Construcción de Soluciones CGO: Los autores construyen soluciones $u_1$ y $u_2$ de las ecuaciones adjunta y directa, respectivamente, de la forma $u(x) = e^{ix\cdot\zeta/h}(a(x) + r(x))$, donde $\zeta$ es un vector complejo que satisface $\zeta \cdot \zeta = 0$, $h$ es un parámetro semiclasico, y $r(x)$ es un término de residuo que se anula cuando $h \to 0$.
2. Identidad Integral y Fórmula de Green: Mediante la aplicación de la fórmula de Green a la diferencia de dos operadores y el uso del mapa DtN parcial, los autores derivan una identidad integral que relaciona la diferencia de los coeficientes ($A_1 - A_2$ y $q_1 - q_2$) con la diferencia de los mapas DtN.
3. Continuidad Única Cuantitativa Local: Una innovación crucial en este trabajo es la derivación de una estimación de continuidad única cuantitativa local para soluciones de la ecuación biarmónica no homogénea con condiciones de frontera homogéneas. Esto se logra utilizando estimaciones de Carleman para el operador biarmónico semiclasico $h^4\Delta^2$. Este paso permite a los autores controlar las normas de la solución específicamente sobre el subconjunto de medición $\gamma_2$ en lugar de sobre toda la frontera, lo cual es esencial para el escenario de datos parciales.
4. Descomposición de Campos Vectoriales: Para recuperar el campo vectorial $A$ a partir de su derivada exterior $dA$ (que se obtiene naturalmente mediante la transformada de Fourier de la identidad integral), los autores utilizan una descomposición de tipo Helmholtz $A = A_{sol} + \nabla \phi$. Estiman la parte solenoidal $A_{sol}$ utilizando la cota sobre $dA$ y la parte de potencial $\nabla \phi$ resolviendo un problema de valor en la frontera para $\phi$ y aplicando desigualdades de incrustación de Sobolev e interpolación.
5. Optimización de Parámetros: Las estimaciones de estabilidad se derivan optimizando el parámetro semiclasico $h$ y el corte de frecuencia $\rho$ en el dominio de Fourier, equilibrando los términos de error provenientes de los residuos CGO y los pesos exponenciales en las estimaciones de Carleman.

Contribuciones Clave y Resultados
El artículo establece estimaciones de estabilidad de tipo logarítmico y de tipo doble logarítmico, las cuales son características de los problemas inversos mal condicionados.

• Estabilidad para la Perturbación de Primer Orden ($A$):
  El Teorema 1.1 establece una estimación de estabilidad de tipo log para el campo vectorial $A$. Específicamente, la norma $L^\infty$ de la diferencia $A_1 - A_2$ está acotada por una suma de una potencia de la diferencia del DtN y una potencia negativa del logaritmo de la diferencia del DtN:
  $$\|A_1 - A_2\|_{L^\infty} \lesssim \|\Lambda_1 - \Lambda_2\|^{\mu_1} + |\log \|\Lambda_1 - \Lambda_2\||^{-\mu_2}$$
  Este resultado complementa los resultados de unicidad previos al proporcionar una cota de estabilidad cuantitativa.

• Estabilidad para la Perturbación de Orden Cero ($q$) con Término de Primer Orden:
  El Teorema 1.2 establece una estimación de estabilidad de tipo log-log para el potencial $q$ cuando el término de primer orden $A$ está presente:
  $$\|q_1 - q_2\|_{H^{-1}} \lesssim |\log |\log \|\Lambda_1 - \Lambda_2\|||^{-\mu'}$$
  El Corolario 1.3 extiende esto a una estimación $L^\infty$ bajo supuestos adicionales de regularidad para $q$.

• Estabilidad para la Perturbación de Orden Cero ($q$) sin Término de Primer Orden:
  El Teorema 1.4 muestra que en ausencia de la perturbación de primer orden ($A=0$), la estabilidad para $q$ mejora de tipo log-log a tipo log:
  $$\|q_1 - q_2\|_{H^{-1}} \lesssim \|\Lambda_1 - \Lambda_2\|^{2/3} + |\log \|\Lambda_1 - \Lambda_2\||^{-2/(n+2)}$$
  Esto resalta la dificultad específica introducida por el término de primer orden en el operador biarmónico.

Significado y Pretensiones
Los autores afirman que este trabajo proporciona los primeros resultados de estabilidad para la perturbación de primer orden del operador biarmónico utilizando datos parciales. Si bien la unicidad para este problema ya había sido establecida (citando [42]), el análisis de estabilidad es significativamente más desafiante.

El artículo enfatiza dos innovaciones técnicas principales:

1. Continuidad Única Local: La derivación de una estimación de continuidad única cuantitativa local para operadores elípticos de orden superior (específicamente el operador biarmónico) utilizando estimaciones de Carleman. Los autores señalan que este es el primer resultado de este tipo para operadores de orden superior en este contexto específico de datos parciales.
2. Manejo de Datos Parciales: La capacidad de derivar estimaciones de estabilidad utilizando mediciones en subconjuntos abiertos arbitrariamente pequeños de la frontera, en lugar de requerir datos en toda la frontera o configuraciones geométricas específicas (como partes inaccesibles planas).

Los autores posicionan sus resultados como un complemento a los teoremas de unicidad existentes y un paso adelante en el análisis de estabilidad de problemas inversos de orden superior, trazando paralelismos con los métodos utilizados para los operadores de Schrödinger y de Schrödinger magnético, mientras abordan los desafíos únicos planteados por la estructura del operador biarmónico y la falta de invarianza de gauge para su perturbación de primer orden.