Inverse problems for the spectral fractional Laplacian with inhomogeneous Dirichlet boundary data

Aquí presento un resumen técnico detallado del artículo "Inverse Problems for the Spectral Fractional Laplacian with Inhomogeneous Dirichlet Boundary Data" (Problemas inversos para el Laplaciano fraccional espectral con datos de frontera de Dirichlet inhomogéneos), escrito por Ravi Shankar Jaiswal, Pu-Zhao Kow y Suman Kumar Sahoo.

---

1. Planteamiento del Problema

El artículo aborda un problema inverso de tipo Calderón para el Laplaciano fraccional espectral con datos de frontera de Dirichlet inhomogéneos. El objetivo principal es determinar un potencial desconocido $q$ en un dominio acotado $\Omega \subset \mathbb{R}^n$ ($n \ge 2$) a partir de mediciones en la frontera.

El problema directo se formula como la ecuación de Schrödinger fraccional:
$$\begin{cases} ((-\Delta_D)^s - q) u = 0 & \text{en } \Omega, \\ u|_{\partial\Omega} = g, \end{cases}$$
donde:

• $(-\Delta_D)^s$ es el Laplaciano fraccional espectral con condiciones de frontera de Dirichlet inhomogéneas (definido a través de la descomposición espectral del Laplaciano clásico con condiciones de Dirichlet homogéneas).
• $s \in (1/2, 1)$ es el orden de la fracción.
• $g$ es el dato de frontera.
• $q \in L^\infty(\Omega)$ es el potencial a recuperar.

Se define el Mapa Dirichlet-a-Neumann (DN) asociado, denotado como $\Lambda_q^s$, que mapea el dato de frontera $g$ a una derivada normal generalizada (o traza de Neumann) de la solución. El problema inverso consiste en determinar si el mapa $\Lambda_q^s$ determina unívocamente el potencial $q$.

2. Metodología

Los autores emplean una combinación de análisis funcional, teoría espectral de operadores diferenciales y técnicas de ecuaciones en derivadas parciales (EDP) avanzadas:

• Aproximación de Born y Linealización: Dado que el problema es no lineal en $q$, los autores estudian la linealización del mapa DN alrededor de $q=0$. Definen la derivada de Fréchet $d\Lambda^s[q]$ (conocida como la aproximación de Born), que es un operador lineal que mapea perturbaciones del potencial a perturbaciones en el mapa de frontera.
• Identidad de Alessandrini: Se utiliza una identidad integral que relaciona la diferencia entre dos mapas DN con la diferencia de los potenciales integrados contra productos de soluciones. Esto permite transformar el problema de recuperación de $q$ en un problema de determinar una función a partir de sus transformadas de Fourier ponderadas.
• Soluciones de Óptica Geométrica Compleja (CGO): La herramienta central para la prueba de unicidad es la construcción de soluciones especiales de alta frecuencia (CGO) para la ecuación de Laplace ($\Delta u = 0$) y la ecuación biarmónica ($\Delta^2 u = 0$). Estas soluciones tienen la forma $u(x) = e^{\alpha \cdot x}(1 + r(x))$, donde $\alpha$ es un vector complejo y $r(x)$ es un término de resto pequeño.
• Análisis Asintótico: Se utiliza un parámetro pequeño $h$ (semiclásico) para expandir las soluciones y analizar los términos dominantes en el límite $h \to 0$. Esto permite extraer información sobre los coeficientes del operador a partir de la identidad integral.
• Análisis de Dimensiones: Se trata por separado el caso de dimensión $n \ge 3$ y el caso crítico $n=2$, requiriendo técnicas ligeramente diferentes debido a la estructura de las soluciones CGO y las simetrías en dos dimensiones.

3. Contribuciones Clave

El artículo presenta dos contribuciones principales:

A. Unicidad y Estabilidad para el Potencial Escalar (Teorema 1.1)

Los autores demuestran que el mapa DN linealizado $d\Lambda^s$ es inyectivo. Específicamente:

• Si dos potenciales $q^{(1)}$ y $q^{(2)}$ (acotados en norma $L^2$) producen el mismo mapa DN linealizado para todos los datos de frontera suaves, entonces $q^{(1)} \equiv q^{(2)}$.
• Se establece un resultado de estabilidad logarítmica (tipo Hölder logarítmico):
  $$\|q^{(1)} - q^{(2)}\|_{H^{-1}(\mathbb{R}^n)} \le C w(\|d\Lambda^s[q^{(1)}] - d\Lambda^s[q^{(2)}]\|_*),$$
  donde la función de continuidad es $w(t) \sim |\log t|^{-1}$. Esto confirma que el problema es mal planteado (ill-posed), pero estable bajo perturbaciones controladas.

B. Recuperación de Coeficientes de Orden Superior (Teorema 1.4)

El segundo resultado aborda un problema más general donde el operador no es solo un potencial escalar, sino que involucra términos de primer y segundo orden derivados de una función de gauge $w$.

• Se considera la identidad integral:
  $$\int_\Omega (\theta_2 : \nabla^{\otimes 2} u + \theta_1 \cdot \nabla u + \theta_0 u) v \, dx = 0,$$
  para funciones $u$ armónicas ($\Delta^2 u = 0$) y $v$ tales que $(-\Delta_D)^s v = 0$.
• Resultado: Bajo ciertas condiciones de soporte y regularidad en los coeficientes $\theta_2, \theta_1, \theta_0$, se demuestra que estos coeficientes pertenecen a una clase de gauge específica definida por una función $w$ (es decir, $\theta_2 = w I_n$, $\theta_1 = 2\nabla w$, $\theta_0 = \Delta w$).
• Caso $n=2$: Bajo la suposición adicional de que $\theta_0 = 0$ y la traza de $\theta_2$ es cero, se demuestra que los coeficientes son idénticamente nulos ($\theta_1 = 0, \theta_2 = 0$).

4. Resultados Técnicos Destacados

• Definición Rigurosa del Mapa DN: Se proporciona una definición precisa del mapa DN para el Laplaciano fraccional espectral con datos inhomogéneos, utilizando el marco de [APR18] y operadores de Poisson generalizados.
• Propiedad de Aproximación de Runge: Se utiliza implícitamente la propiedad de aproximación de Runge (una forma cuantitativa de la Propiedad de Continuación Única - UCP) para asegurar que las soluciones de la ecuación fraccional pueden aproximar funciones armónicas, lo cual es crucial para conectar el problema fraccional con soluciones CGO clásicas.
• Construcción de Soluciones CGO para el Laplaciano Fraccional: En el Apéndice A, los autores refinan la construcción de soluciones CGO para la ecuación $\Delta^k u = 0$ en dominios acotados, adaptando los resultados de [SS23] para manejar los términos de transporte y los restos en normas semiclásicas $H^m_{[h]}$.
• Diferencia entre Dimensiones: Se destaca la dificultad técnica en $n=2$, donde la identificación exacta del gauge requiere condiciones adicionales sobre la traza del tensor $\theta_2$, a diferencia de $n \ge 3$ donde la estructura algebraica de las soluciones permite una recuperación más directa.

5. Significado e Impacto

Este trabajo es significativo por varias razones:

1. Avance en Problemas Inversos Fraccionales: Extiende la teoría de problemas inversos (clásicamente desarrollada para el Laplaciano entero) al régimen fraccional espectral con condiciones de frontera inhomogéneas, un escenario que es físicamente relevante pero matemáticamente más complejo que el caso de datos exteriores (donde se usa el Laplaciano fraccional integral).
2. Conexión entre Operadores Locales y No Locales: El resultado demuestra que, a pesar de la naturaleza no local del operador fraccional, la información en la frontera es suficiente para recuperar el potencial interior, aprovechando la fuerte propiedad de continuación única.
3. Generalización de Resultados Previos: El Teorema 1.4 generaliza resultados recientes sobre la recuperación de coeficientes en ecuaciones de cuarto orden (como en [SS23]) al contexto de operadores fraccionales, mostrando cómo la estructura de gauge se mantiene o se modifica en este marco.
4. Aplicaciones Potenciales: Los resultados tienen implicaciones para la imagen médica y la geofísica, donde los modelos fraccionales pueden describir mejor la difusión anómala o la propagación de ondas en medios heterogéneos que los modelos clásicos.

En resumen, el artículo establece las bases teóricas sólidas para la resolución de problemas inversos para el Laplaciano fraccional espectral, proporcionando pruebas de unicidad, estabilidad y métodos constructivos basados en soluciones de óptica geométrica compleja.