Inverse Problems for the Monge--Ampère Equation: Linearization and Nonlinear Recovery

Technische Zusammenfassung: Inverse Probleme für die Monge–Ampère-Gleichung: Linearisierung und nichtlineare Rekonstruktion

Problemstellung
Diese Arbeit untersucht inverse Randwertprobleme für die voll nichtlineare Monge–Ampère-Gleichung:
$$\det D^2 u = a(x, u, \nabla u) \quad \text{in } \Omega,$$
wobei $\Omega \subset \mathbb{R}^n$ eine beschränkte Domäne ist und $a$ eine unbekannte Nichtlinearität darstellt. Motiviert durch das klassische Calderón-Problem (elektrische Impedanztomographie), untersuchen die Autoren, ob der „Jet“ der Nichtlinearität $a$ entlang einer gegebenen streng konvexen Hintergrundlösung $\Phi$ aus Randmessungen rekonstruiert werden kann. Konkret ist das Ziel, die Ableitungen von $a$ bezüglich seiner Argumente $(x, z, p)$ entlang der Trajektorie $(x, \Phi(x), \nabla \Phi(x))$ zu bestimmen, unter Verwendung des nichtlinearen Cauchy-Datensatzes.

Im Gegensatz zu semilinearen oder quasilinearen Gleichungen besitzen voll nichtlineare Gleichungen im Allgemeinen keine kanonische Dirichlet-zu-Neumann-Abbildung (DN-Abbildung) und lassen sich nicht ohne Weiteres mit linearen inversen Problemen verknüpfen. Die Monge–Ampère-Gleichung stellt aufgrund ihrer voll nichtlinearen Natur eine besondere Herausforderung dar, doch sie besitzt eine Divergenzstruktur über die Cofaktor-Matrix der Hesse-Matrix, welche die Autoren ausnutzen, um einen systematischen Rahmen zu etablieren.

Methodik
Die Autoren entwickeln einen systematischen Rahmen für die höherwertige Linearisierung um eine streng konvexe Hintergrundlösung $\Phi$. Der Ansatz erfolgt in drei Hauptphasen:

1. Lokale Wohldefiniertheit und glatte Abhängigkeit:
   Die Autoren etablieren zunächst die lokale Wohldefiniertheit des Dirichlet-Problems für Randdaten nahe $\Phi|_{\partial \Omega}$. Sie beweisen, dass der Lösungsoperator $S: \eta \mapsto u_\eta$ in Bezug auf die Randdaten glatt ($C^k$) ist. Dies ermöglicht die Ableitung expliziter Formeln für die höherwertigen Fréchet-Ableitungen des Lösungsoperators.

2. Linearisierung der Cauchy-Daten:
   Ein nichtlinearer Cauchy-Datensatz wird als die Menge der Paare $(\eta, \partial_\nu u_\eta|_{\partial \Omega})$ für Randdaten nahe $\Phi|_{\partial \Omega \text{}}$ definiert. Die Autoren zeigen, dass die erste Variation dieses nichtlinearen Datensatzes mit dem Cauchy-Datensatz eines linearen elliptischen Operators $L_\Phi$ übereinstimmt, der die linearisierte Gleichung regelt:
   $$L_\Phi v = \text{div}(C_\Phi \nabla v) - b_\Phi \cdot \nabla v - q_\Phi v = 0,$$
   wobei die Koeffizienten durch die Hintergrundlösung und die Nichtlinearität bestimmt werden:

• $C_\Phi = \text{cof } D^2 \Phi$ (die Cofaktor-Matrix der Hesse-Matrix),
• $b_\Phi = \nabla_p a(x, \Phi, \nabla \Phi)$,
• $q_\Phi = \partial_z a(x, \Phi, \nabla \Phi)$.

Entscheidend ist, dass der Hauptteil dieses linearisierten Operators anisotrop ist, bestimmt durch die Cofaktor-Matrix $C_\Phi$.

3. Reduktion und höherwertige Identitäten:
   Unter Verwendung einer Alessandrini-Typ-Identität etablieren die Autoren ein Reduktionsprinzip: Die Gleichheit der nichtlinearen Cauchy-Datensätze impliziert die Gleichheit der entsprechenden linearisierten DN-Abbildungen. Dies reduziert das nichtlineare inverse Problem auf ein anisotropes Calderón-artiges inverses Problem für den Operator $L_\Phi$.

Durch die Verwendung von Integralidentitäten, die Produkte von Lösungen der linearisierten Gleichung und deren Adjungierten involvieren, etablieren die Autoren unter geeigneten Dichtebedingungen für diese Produkte die Möglichkeit, höhere Ableitungen von $a$ entlang des Hintergrund-Jets zu rekonstruieren.

Wesentliche Ergebnisse

• Rekonstruktion erster Ordnung: Wenn zwei Nichtlinearitäten $a_1$ und $a_2$ dieselben nichtlinearen Cauchy-Daten nahe einer Hintergrundlösung $\Phi$ erzeugen, dann fallen die DN-Abbildungen ihrer assoziierten linearisierten Operatoren zusammen. Folglich sind die Größen $\partial_z a(x, \Phi, \nabla \Phi)$ und $\nabla_p a(x, \Phi, \nabla \Phi)$ entlang des Hintergrund-Jets bestimmt, bis auf die natürliche Gauge-Invarianz des anisotropen Calderón-Problems (randfixierende Diffeomorphismen).
• Höherwertige Rekonstruktion: Wenn die Hintergrundlösung fixiert ist ($\Phi_1 = \Phi_2 = \Phi$) und die Ableitungen niedrigerer Ordnung bereits als identisch bekannt sind, ermöglichen höherwertige Linearisierungsidentitäten die Rekonstruktion höherer Ableitungen $d^k_{(z,p)} a(x, \Phi, \nabla \Phi)$. Dies beruht auf einer Dichte-Eigenschaft von Produkten von Lösungen der linearisierten Gleichung und deren Adjungierten.
• Spezifische Anwendungen:
  • Semilinearer Fall: Wenn $a$ unabhängig von $\nabla u$ ist (semilinear), reduziert sich der linearisierte Operator auf eine Schrödinger-Gleichung. In Dimensionen $n \ge 3$ liefert dies Unizität für die Ableitung erster Ordnung $\partial_z a$ entlang quadratischer Hintergründe ($\Phi(x) = |x|^2/2$).
  • Optimal Transport: Für Nichtlinearitäten, die aus dem optimalen Transport stammen, speziell in der Gleichgewichtskonfiguration $\Phi(x) = |x|^2/2$, reduziert sich die linearisierte Gleichung auf eine standardmäßige isotrope Leitfähigkeitsgleichung. Dies stellt eine direkte Verbindung zum klassischen Calderón-Problem her und ermöglicht die Rekonstruktion von Transportdichten.
  • Inverse Quellprobleme: Im Falle von $a(x, u, \nabla u) = f(x)$ reduziert sich das Problem auf die Rekonstruktion des Quellterms $f(x)$, was mit den Starrigkeitseigenschaften der Klasse der Cofaktor-Matrizen streng konvexer Funktionen verknüpft ist.

Bedeutung und Behauptungen
Das Paper beansprucht, die erste systematische Untersuchung inverser Probleme für voll nichtlineare Monge–Ampère-Gleichungen aus der Perspektive von Calderón-artigen Randmessungen zu liefern. Sein primärer Beitrag besteht darin, zu zeigen, dass das inverse Problem trotz der voll nichtlinearen Natur der Gleichung systematisch auf eine Familie linearer anisotroper inverser Probleme durch höherwertige Linearisierung reduziert werden kann.

Die Autoren betonen, dass ihr Ansatz die Lücke zwischen voll nichtlinearen PDEs und linearen inversen Problemen schließt, indem er die spezifische Divergenzstruktur des Monge–Ampère-Operators nutzt. Während die Rekonstruktion der vollen Nichtlinearität von der Lösbarkeit des zugehörigen anisotropen Calderón-Problems abhängt (welches in allgemeinen Dimensionen ohne zusätzliche geometrische Annahmen weiterhin offen ist), stellt das Paper fest, dass die nichtlinearen Daten die linearisierten Koeffizienten bis auf die natürliche Gauge bestimmen. Darüber hinaus zeigen sie auf, dass unter spezifischen geometrischen Einstellungen (wie etwa quadratischen Hintergründen) und strukturellen Annahmen (Semilinearitäten oder Gleichgewichtskonfigurationen des optimalen Transports) unbedingte Unizitätsresultate erzielt werden können, die die vollständige Taylor-Entwicklung der Nichtlinearität entlang der Hintergrundtrajektorie rekonstruieren.