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

Sintesi Tecnica: Determinazione Stabile della Perturbazione del Primo Ordine dell'Operatore Biharmonico da Dati Parziali

Enunciato del Problema
Questo articolo affronta un problema inverso di valore al contorno per l'operatore biarmonico perturbato in un dominio limitato $\Omega \subseteq \mathbb{R}^n$ ($n \ge 3$) con un bordo liscio $\partial\Omega$. L'operatore è definito come $L_{A,q} = \Delta^2 + A(x) \cdot D + q(x)$, dove $D = i^{-1}\nabla$, $A$ è un campo vettoriale in $W^{1,\infty}(\Omega, \mathbb{C}^n)$, e $q$ è un potenziale in $L^\infty(\Omega, \mathbb{C})$. Lo studio si concentra sul contesto dei dati parziali: il recupero dei coefficienti $A$ e $q$ da misurazioni effettuate solo su sottoinsiemi aperti arbitrariamente piccoli $\gamma_1, \gamma_2 \subset \partial\Omega$.

Il problema assume che i coefficienti siano a priori noti in un intorno $\omega_0$ del bordo. I dati sono forniti da una mappa Dirichlet-to-Neumann (DtN) parziale $\Lambda_{\gamma_1, \gamma_2}^{A,q}$, che mappa i dati di Dirichlet $(f, g)$ supportati su $\gamma_1$ nelle derivate normali $(\partial_\nu u, \partial_\nu(\Delta u))$ misurate su $\gamma_2$. L'obiettivo primario è stabilire stime di stabilità (quantificando come gli errori nelle misurazioni influenzano la ricostruzione dei coefficienti) sia per la perturbazione del primo ordine $A$ che per la perturbazione del termine di ordine zero $q$.

Metodologia
Gli autori utilizzano il metodo delle soluzioni di Geometria Complessa Ottica (CGO), una tecnica standard per i problemi inversi di operatori ellittici. La strategia di prova coinvolge diversi componenti tecnici chiave:

1. Costruzione di Soluzioni CGO: Gli autori costruiscono soluzioni $u_1$ e $u_2$ per le equazioni adiugata e diretta, rispettivamente, della forma $u(x) = e^{ix\cdot\zeta/h}(a(x) + r(x))$, dove $\zeta$ è un vettore complesso soddisfacente $\zeta \cdot \zeta = 0$, $h$ è un parametro semiclassico, e $r(x)$ è un termine di resto che svanisce per $h \to 0$.
2. Identità Integrale e Formula di Green: Applicando la formula di Green alla differenza di due operatori e utilizzando la mappa DtN parziale, gli autori derivano un'identità integrale che mette in relazione la differenza dei coefficienti ($A_1 - A_2$ e $q_1 - q_2$) con la differenza delle mappe DtN.
3. Unicità di Continuazione Quantitativa Locale: Un'innovazione cruciale in questo lavoro è la derivazione di una stima di unicità di continuazione quantitativa locale per la soluzione dell'equazione biarmonica non omogenea con condizioni al contorno omogenee. Ciò viene ottenuto utilizzando stime di Carleman per l'operatore biarmonico semiclassico $h^4\Delta^2$. Questo passaggio permette agli autori di controllare le norme delle soluzioni specificamente su un sottoinsieme di misurazione $\gamma_2$ piuttosto che sull'intero bordo, il che è essenziale per il contesto dei dati parziali.
4. Decomposizione dei Campi Vettoriali: Per recuperare il campo vettoriale $A$ dalla sua derivata esterna $dA$ (che si ottiene naturalmente tramite la trasformata di Fourier dell'identità integrale), gli autori utilizzano una decomposizione di tipo Helmholtz $A = A_{sol} + \nabla \phi$. Stimano la parte solenoide $A_{sol}$ usando il limite su $dA$ e la parte potenziale $\nabla \phi$ risolvendo un problema al contorno per $\phi$ e applicando le disuguaglianze di embedding di Sobolev e di interpolazione.
5. Ottimizzazione dei Parametri: Le stime di stabilità vengono derivate ottimizzando il parametro semiclassico $h$ e il cutoff di frequenza $\rho$ nel dominio di Fourier, bilanciando i termini di errore derivanti dai residui CGO e i pesi esponenziali nelle stime di Carleman.

Contributi Chiave e Risultati
L'articolo stabilisce stime di stabilità di tipo logaritmico e di tipo doppio-logaritmico, che sono caratteristiche dei problemi inversi mal posti.

• Stabilità per la Perturbazione del Primo Ordine ($A$):
  Il Teorema 1.1 stabilisce una stima di stabilità di tipo log per il campo vettoriale $A$. Nello specifico, la norma $L^\infty$ della differenza $A_1 - A_2$ è limitata da una somma di una potenza della differenza DtN e di una potenza negativa del logaritmo della differenza DtN:
  $$\|A_1 - A_2\|_{L^\infty} \lesssim \|\Lambda_1 - \Lambda_2\|^{\mu_1} + |\log \|\Lambda_1 - \Lambda_2\||^{-\mu_2}$$
  Questo risultato integra i precedenti risultati di unicità fornendo un limite quantitativo di stabilità.

• Stabilità per la Perturbazione del Termine di Ordine Zero ($q$) con Termine del Primo Ordine:
  Il Teorema 1.2 stabilisce una stima di stabilità di tipo log-log per il potenziale $q$ in presenza del termine del primo ordine $A$:
  $$\|q_1 - q_2\|_{H^{-1}} \lesssim |\log |\log \|\Lambda_1 - \Lambda_2\|||^{-\mu'}$$
  Il Corollario 1.3 estende questo risultato a una stima $L^\infty$ sotto ulteriori assunzioni di regolarità su $q$.

• Stabilità per la Perturbazione del Termine di Ordine Zero ($q$) senza Termine del Primo Ordine:
  Il Teorema 1.4 mostra che in assenza della perturbazione del primo ordine ($A=0$), la stabilità per $q$ migliora da 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)}$$
  Ciò evidenzia la specifica difficoltà introdotta dal termine del primo ordine nell'operatore biarmonico.

Significato e Rivendicazioni
Gli autori affermano che questo lavoro fornisce i primi risultati di stabilità per la perturbazione del primo ordine dell'operatore biarmonico utilizzando dati parziali. Sebbene l'unicità per questo problema sia stata precedentemente stabilita (citando [42]), l'analisi di stabilità è significativamente più complessa.

L'articolo sottolinea due principali novità tecniche:

1. Unicità di Continuazione Locale: La derivazione di una stima di unicità di continuazione quantitativa locale per operatori ellittici di ordine superiore (specificamente l'operatore biarmonico) utilizzando stime di Carleman. Gli autori notano che questa è la prima tale risultante per operatori di ordine superiore in questo specifico contesto di dati parziali.
2. Gestione dei Dati Parziali: La capacità di derivare stime di stabilità utilizzando misurazioni su sottoinsiemi aperti arbitrariamente piccoli del bordo, anziché richiedere dati sull'intero bordo o configurazioni geometriche specifiche (come parti inaccessibili piatte).

Gli autori posizionano i loro risultati come un complemento ai teoremi di unicità esistenti e un passo avanti nell'analisi di stabilità di problemi inversi di ordine superiore, tracciando paralleli con i metodi utilizzati per gli operatori di Schrödinger e di Schrödinger magnetico, pur affrontando le sfide uniche poste dalla struttura dell'operatore biarmonico e dalla mancanza di invarianza di gauge per la sua perturbazione del primo ordine.