An inverse problem for semilinear wave equations on metric tree graphs

1. Il Problema

Il lavoro affronta un problema inverso per un'equazione delle onde semilineare definita su grafi metrici ad albero (metric tree graphs). L'obiettivo è determinare la struttura geometrica del grafo (connettività e lunghezze degli spigoli) e i coefficienti incogniti dell'equazione partendo da dati di misura al contorno.

L'equazione modello considerata è:
$$(\partial_t^2 - \partial_x^2)u(x, t) + q(x)u(x, t) + a(x, t)u^3(x, t) = 0$$
definita su un grafo metrico $\Omega$ fino a un tempo $T > 0$.

• $q(x)$ è un potenziale indipendente dal tempo.
• $a(x, t)$ è un coefficiente dipendente dal tempo associato al termine non lineare cubico.
• Le condizioni al contorno sono di Dirichlet su un sottoinsieme di vertici $\Gamma_0$ (foglie del grafo) e condizioni di Kirchhoff-Neumann (continuità della soluzione e somma nulla delle derivate normali) sui vertici interni.
• Su una singola foglia $\gamma_0$, viene imposta una condizione di Dirichlet omogenea ($u=0$), mentre su tutte le altre foglie $\Gamma_0 = \Gamma \setminus \{\gamma_0\}$ si applicano sorgenti di controllo e si misurano le derivate normali.

L'obiettivo inverso: Data la mappa Dirichlet-to-Neumann (DtN) $\Lambda_T$ definita su $\Gamma_0$ fino al tempo $T$, è possibile recuperare:

1. La struttura del grafo metrico (fino all'isometria, ignorando i vertici di grado 2).
2. Il potenziale $q(x)$.
3. Il coefficiente non lineare $a(x, t)$ in un dominio temporale spaziale specifico.

2. Metodologia

Gli autori combinano tecniche di analisi funzionale, teoria delle onde geometriche e metodi di linearizzazione per risolvere il problema. La strategia si articola in diversi passaggi fondamentali:

A. Linearizzazione e Riduzione al Caso Lineare

Il primo passo cruciale è sfruttare la struttura della non linearità. Poiché il termine non lineare è cubico ($u^3$), gli autori utilizzano una linearizzazione multipla (derivata terza rispetto a parametri di scala delle sorgenti).

• Dimostrano che la mappa DtN per l'equazione semilineare determina univocamente la mappa DtN per l'equazione lineare associata ($a=0$).
• Questo permette di separare il problema: la struttura del grafo e il potenziale $q$ possono essere recuperati utilizzando risultati esistenti per le equazioni lineari (metodo del "leaf-peeling" o "leaf-cleaning").
• Una volta noti la topologia e $q$, il focus si sposta esclusivamente sul recupero del coefficiente $a(x, t)$.

B. Soluzioni di Ottica Geometrica

Per isolare il coefficiente $a(x, t)$, gli autori costruiscono soluzioni ad alta frequenza (onde geometriche) della forma $e^{i h^{-1}\theta(t)} \chi(t)$, dove $h$ è un parametro piccolo.

• Analizzano il comportamento di queste onde quando si propagano attraverso i vertici del grafo ad albero.
• Derivano formule asintotiche esplicite per le onde riflesse e trasmesse ai vertici, tenendo conto delle condizioni di Kirchhoff-Neumann. In particolare, calcolano i coefficienti di riflessione e trasmissione in funzione del grado del vertice.
• Questo permette di tracciare le traiettorie delle onde e prevedere dove le diverse componenti (onde incidenti e riflesse) si intersecheranno.

C. Identità Integrali e Metodo di Intersezione

Il cuore della dimostrazione per il recupero di $a(x, t)$ si basa su un'identità integrale ottenuta tramite integrazione per parti (simile all'identità di Green o al metodo di Calderón):
$$\int_0^T \int_\Omega 6 a v_0 v_1 v_2 v_3 \, dx dt = \int_0^T \int_{\Gamma_0} f_0 \partial_n w \, dt$$
dove $v_1, v_2, v_3$ sono onde lineari inviate dalle sorgenti e $v_0$ è un'onda ausiliaria inviata "all'indietro" nel tempo.

• Scegliendo opportunamente i tempi di invio delle onde ($t_0, s, \dots$) e i supporti delle funzioni di taglio, gli autori fanno sì che i supporti delle quattro onde ($v_0, v_1, v_2, v_3$) si intersechino in una regione spaziale-temporale molto piccola (approssimabile a un punto).
• In questo limite, l'integrale a sinistra converge al valore puntuale di $a(x, t)$ in quel punto.
• Variando i parametri di ritardo, è possibile "scansionare" il dominio e ricostruire $a(x, t)$.

D. Strategia Ricorsiva su Alberi (Sheaf Decomposition)

Per gestire grafi complessi, gli autori utilizzano una decomposizione ricorsiva basata sui fasci di spigoli (sheaves):

1. Identificano un sottografo a stella (un vertice interno collegato a diverse foglie).
2. Recuperano $a$ sugli spigoli incidenti alle foglie controllate e sull'ultimo spigolo che porta alla foglia non controllata, utilizzando le onde riflesse dalle altre foglie.
3. Dimostrano che le onde riflesse da parti più profonde dell'albero non interferiscono con la regione di interesse se i tempi sono scelti correttamente (evitando sovrapposizioni indesiderate).
4. Procedono ricorsivamente verso l'interno dell'albero, trattando ogni sotto-albero come un nuovo problema su un grafo più piccolo, fino a coprire l'intero albero.

3. Risultati Principali

Il teorema principale (Teorema 1) stabilisce quanto segue:

Sia $\Omega$ un albero metrico finito e $\Gamma$ l'insieme delle sue foglie. Se si dispone della mappa DtN su $\Gamma_0 = \Gamma \setminus \{\gamma_0\}$ per un tempo $T > 2D(\gamma_0)$ (dove $D(\gamma_0)$ è la massima distanza da $\gamma_0$ alle altre foglie), allora:

1. La struttura metrica dell'albero (connettività e lunghezze degli spigoli) è univocamente determinata, a meno di vertici di grado 2.
2. Il potenziale $q(x)$ è univocamente determinato su tutto l'albero.
3. Il coefficiente $a(x, t)$ è univocamente determinato nel dominio:
   $$\{(x, t) \in \Omega \times [0, T] : D(\gamma_0) + d(x, \gamma_0) \le t \le T - D(\gamma_0) + d(x, \gamma_0)\}$$
   Questo dominio rappresenta la regione raggiungibile dalle onde che partono dalle foglie controllate, si riflettono e tornano indietro, entro il tempo $T$.

4. Contributi Chiave

• Primo risultato su equazioni non lineari: Questo è uno dei primi lavori che risolve un problema inverso per equazioni differenziali non lineari su grafi metrici. La maggior parte della letteratura precedente si era concentrata su equazioni lineari.
• Estensione del metodo di linearizzazione: Gli autori adattano con successo il metodo di linearizzazione sviluppato per varietà Riemanniane (in particolare il lavoro di Lassas et al.) al contesto discreto e complesso dei grafi, affrontando le sfide poste dalle riflessioni multiple ai vertici.
• Gestione delle riflessioni complesse: Dimostrano come isolare le interazioni delle onde su alberi complessi, mostrando che è possibile "pulire" l'albero foglia per foglia senza che le riflessioni da parti distanti dell'albero oscurino l'informazione locale.
• Robustezza rispetto alla topologia: Il metodo funziona per alberi di qualsiasi configurazione, purché si abbia controllo su tutte le foglie tranne una.

5. Significato e Implicazioni

• Applicazioni Scientifiche: I grafi metrici sono modelli fondamentali in fisica (es. guide d'onda, fibre ottiche, reti di nanotubi di carbonio) e ingegneria. La capacità di ricostruire parametri non lineari da dati di superficie è cruciale per il monitoraggio strutturale e l'identificazione di difetti in materiali complessi.
• Avanzamento Teorico: Il lavoro colma un vuoto significativo nella teoria dei problemi inversi, estendendo la comprensione dai sistemi lineari a quelli non lineari in geometrie non lisce (grafi).
• Metodologia Innovativa: La combinazione di ottica geometrica, linearizzazione di ordine superiore e decomposizione ricorsiva offre un nuovo paradigma per affrontare problemi inversi su strutture reticolari, con potenziali estensioni future a grafi con cicli (non solo alberi).

In sintesi, il paper fornisce una soluzione costruttiva e rigorosa per un problema inverso altamente non banale, dimostrando che la conoscenza della risposta dinamica al contorno è sufficiente per ricostruire sia la geometria che i parametri fisici non lineari di un sistema a rete.