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

1. Il Problema Studiato

Il lavoro si concentra sullo studio di problemi inversi per l'operatore Laplaciano frazionario spettrale ($(-\Delta_D)^s$) in un dominio limitato $\Omega \subset \mathbb{R}^n$ ($n \ge 2$) con bordo liscio, considerando dati al contorno di Dirichlet non omogenei.

Il problema diretto è formulato come un'equazione di Schrödinger frazionaria:
$$\begin{cases} ((-\Delta_D)^s - q) u = 0 & \text{in } \Omega, \\ u|_{\partial\Omega} = g, \end{cases}$$
dove $q \in L^\infty(\Omega)$ è un potenziale incognito e $g$ è il dato al contorno. L'obiettivo è determinare il potenziale $q$ (o i coefficienti di un'equazione differenziale associata) misurando la mappa Dirichlet-to-Neumann (DN) associata, che lega il dato al contorno $g$ alla sua "derivata normale" generalizzata.

A differenza dei lavori precedenti che trattano il Laplaciano frazionario in senso integrale o con dati al contorno esterni, questo studio utilizza la definizione spettrale basata sugli autovalori del Laplaciano classico con condizioni di Dirichlet omogenee, estendendola a dati non omogenei.

2. Metodologia

Gli autori adottano un approccio analitico sofisticato che combina teoria spettrale, analisi funzionale e soluzioni di ottica geometrica complessa (CGO).

• Definizione della Mappa DN: Viene definita rigorosamente la mappa $\Lambda_q$ che mappa i dati di Dirichlet $g \in H^{s-1/2}(\partial\Omega)$ in una quantità che rappresenta la derivata normale generalizzata. La definizione si basa sulla scomposizione della soluzione $u$ in una parte armonica (rispetto all'operatore frazionario) e una parte corretta dal potenziale.
• Linearizzazione (Approssimazione di Born): Per risolvere il problema inverso, gli autori studiano la linearizzazione della mappa DN attorno al potenziale nullo ($q=0$). La derivata di Fréchet $d\Lambda$ è un operatore lineare che mappa il potenziale $q$ in una perturbazione della mappa DN. L'obiettivo principale è dimostrare l'iniettività di questo operatore lineare.
• Identità di Alessandrini: Viene derivata un'identità fondamentale (analoga a quella classica per l'equazione di Schrödinger) che collega la differenza tra le mappe DN di due potenziali all'integrale del prodotto dei potenziali per soluzioni dell'equazione omogenea.
• Soluzioni CGO (Complex Geometric Optics): La parte centrale della dimostrazione consiste nella costruzione di soluzioni speciali (CGO) per l'equazione $\Delta^2 u = 0$ (nel caso di recupero dei coefficienti) e per l'equazione frazionaria. Queste soluzioni sono della forma $e^{\phi/h} a(x, h)$, dove $h$ è un parametro semiclassico piccolo.
  • Per $n \ge 3$, vengono costruite soluzioni asintotiche basate su operatori di trasporto.
  • Per $n = 2$, viene utilizzata una tecnica di fase stazionaria (Lemma 4.5) con funzioni di fase specifiche nel piano complesso.
• Analisi Asintotica: Sostituendo le soluzioni CGO nell'identità di Alessandrini e facendo tendere $h \to 0$, si ottengono equazioni integrali che permettono di recuperare i coefficienti incogniti punto per punto o fino a un gauge naturale.

3. Contributi Chiave e Risultati Principali

Il documento presenta due risultati principali:

A. Unicità e Stabilità per il Potenziale (Teorema 1.1)

Gli autori dimostrano che il potenziale $q$ è univocamente determinato dalla mappa DN linearizzata $d\Lambda$ (derivata di Fréchet a $q=0$).

• Ipotesi: Il dominio è $C^{1,\alpha}$, $1/2 < s < 1$, e i potenziali sono limitati in norma $L^2$.
• Risultato: Se $d\Lambda[q^{(1)}] = d\Lambda[q^{(2)}]$, allora $q^{(1)} \equiv q^{(2)}$.
• Stabilità: Viene stabilita una stima di stabilità logaritmica (tipica dei problemi inversi mal posti):
  $$\|q^{(1)} - q^{(2)}\|_{H^{-1}} \le C w(\|d\Lambda[q^{(1)}] - d\Lambda[q^{(2)}]\|)$$
  dove il modulo di continuità $w(t)$ è del tipo $|\log t|^{-1}$. Questo conferma che il problema è ben posto in senso debole ma instabile.

B. Recupero dei Coefficienti fino al Gauge Naturale (Teorema 1.4)

Il secondo risultato affronta un problema inverso più generale, cercando di recuperare coefficienti $\theta_2$ (matrice), $\theta_1$ (vettore) e $\theta_0$ (scalare) che appaiono in un'identità di integrazione per parti derivata dall'operatore Laplaciano.

• Contesto: Si studia la recovering di coefficienti nell'identità $\int_\Omega (\theta_2 : \nabla^{\otimes 2} u + \theta_1 \cdot \nabla u + \theta_0 u) v \, dx = 0$ per funzioni armoniche e soluzioni frazionarie.
• Risultato per $n \ge 3$: I coefficienti sono determinati univocamente fino a un gauge naturale definito da una funzione scalare $w$ tale che $w = \partial_\nu w = 0$ sul bordo. In particolare:
  $$\theta_2 = w I_n, \quad \theta_1 = 2\nabla w, \quad \theta_0 = \Delta w.$$
  Questo significa che non si possono distinguere i coefficienti da una trasformazione di gauge generata da $w$.
• Risultato per $n = 2$: Sotto l'ipotesi aggiuntiva che $\theta_0 = 0$ e $\text{tr}(\theta_2) = 0$, si dimostra che $\theta_1 = 0$ e $\theta_2 = 0$. Il caso bidimensionale è più complesso e richiede condizioni aggiuntive per ottenere l'unicità completa.

4. Significato e Impatto

• Estensione della Teoria Classica: Il lavoro estende i risultati classici di Calderón e della teoria dei problemi inversi per l'equazione di Schrödinger al contesto frazionario spettrale, un'area di ricerca relativamente recente ma in rapida crescita.
• Gestione dei Dati Non Omogenei: A differenza di molti studi sul Laplaciano frazionario che si basano su dati al contorno esterni (problema di Cauchy esterno), questo articolo affronta direttamente il caso di dati di Dirichlet non omogenei sul bordo fisico del dominio, che è più rilevante per le applicazioni fisiche reali.
• Nuovi Strumenti Analitici: La costruzione di soluzioni CGO per l'operatore spettrale frazionario con dati al contorno non omogenei e la successiva analisi asintotica rappresentano un avanzamento tecnico significativo.
• Limiti Dimensionali: Il lavoro chiarisce le differenze sostanziali tra il caso $n \ge 3$ e $n=2$, mostrando come la struttura del gauge naturale cambi e richieda ipotesi aggiuntive in due dimensioni.

In sintesi, il paper fornisce una soluzione rigorosa al problema inverso linearizzato per il Laplaciano frazionario spettrale, stabilendo unicità e stabilità logaritmica, e risolvendo un problema di recupero di coefficienti differenziali fino a un gauge naturale, contribuendo significativamente alla teoria matematica dei problemi inversi non locali.