Asymptotic behavior of solutions to elliptic problems with Robin boundary conditions

Titel: Asymptotisches Verhalten von Lösungen elliptischer Probleme mit Robin-Randbedingungen

Autoren: Mengyao Chen, Massimo Grossi, Qi Li

---

1. Problemstellung

Die Autoren untersuchen das asymptotische Verhalten positiver Lösungen $u_\beta$ des semilinearen elliptischen Problems mit Robin-Randbedingungen, wenn der Parameter $\beta \to 0$ geht. Das betrachtete Problem lautet:

$$\begin{cases} -\Delta u = u^p & \text{in } \Omega, \\ u > 0 & \text{in } \Omega, \\ \frac{\partial u}{\partial \nu} + \beta u = 0 & \text{auf } \partial\Omega, \end{cases}$$

wobei:

• $\Omega \subset \mathbb{R}^N$ ($N \ge 2$) ein beschränktes, glattes Gebiet ist.
• $\nu$ die äußere Einheitsnormale auf $\partial\Omega$ bezeichnet.
• $p \ge 0$ ein Exponent ist.
• $\beta > 0$ ein Parameter ist, der die Robin-Bedingung steuert.

Für den Spezialfall $p=1$ wird das Problem durch das zugehörige Eigenwertproblem ersetzt:
$$-\Delta u = \lambda_{1,\beta} u \quad \text{in } \Omega, \quad \frac{\partial u}{\partial \nu} + \beta u = 0 \quad \text{auf } \partial\Omega.$$

Die Robin-Bedingung interpoliert zwischen der Dirichlet-Bedingung ($\beta \to \infty$) und der Neumann-Bedingung ($\beta = 0$). Der Fokus liegt hier auf dem Grenzübergang $\beta \to 0$, der dem Neumann-Fall entspricht.

---

2. Methodik

Die Analyse stützt sich auf eine Kombination aus Variationsmethoden, elliptischer Regularitätstheorie und einer spezifischen Transformation der Lösung.

• Variationsformulierung: Für $0 \le p \le \frac{N+2}{N-2}$ werden Lösungen als kritische Punkte des Funktionals
  $$J_\beta(u) = \frac{1}{2}\int_\Omega |\nabla u|^2 dx + \frac{\beta}{2}\int_{\partial\Omega} u^2 dS - \frac{1}{p+1}\int_\Omega u_+^{p+1} dx$$
  charakterisiert.
• Hilfsfunktion und Transformation: Ein zentrales technisches Werkzeug ist die Einführung einer Hilfsfunktion $v_\beta$, definiert durch:
  $$v_\beta = u_\beta - d_\beta, \quad \text{wobei } d_\beta = \frac{c_\beta}{\beta |\partial\Omega|} \text{ und } c_\beta = \int_\Omega u_\beta^p dx.$$
  Diese Transformation führt zu einem neuen Problem für $v_\beta$ mit der Eigenschaft $\int_{\partial\Omega} v_\beta dS = 0$. Dies ermöglicht es, das Verhalten von $u_\beta$ durch das Verhalten des konstanten Terms $d_\beta$ und des Restterms $v_\beta$ zu analysieren.
• Skalierung und Konvergenz: Für den Fall $p > 1$ wird eine Skalierung $\tilde{u}_\beta = u_\beta / \|u_\beta\|_\infty$ verwendet, um Widersprüche zu erzeugen und das Verschwinden der Lösungen zu beweisen.
• Existenz im superkritischen Fall: Für den Fall, dass $\Omega$ eine Kugel ist und $p > \frac{N+2}{N-2}$ (superkritisch), nutzen die Autoren bekannte Ergebnisse über radiale Lösungen im ganzen Raum $\mathbb{R}^N$ (Serrin-Zou-Lösungen) und passen diese durch Skalierung an die Randbedingung an.

---

3. Hauptergebnisse

Die Arbeit liefert eine vollständige Klassifizierung des asymptotischen Verhaltens von $u_\beta$ für $\beta \to 0$, abhängig vom Exponenten $p$.

A. Eindeutigkeit (Theorem 1.1)

• Für $0 \le p \le \frac{N+2}{N-2}$ ($p \neq 1$) existiert für hinreichend kleines $\beta$ genau eine positive Lösung.
• Für $p=1$ existiert für jedes $\beta$ genau eine positive Eigenfunktion (normiert auf $\|u\|_\infty = 1$).
• Für $p > \frac{N+2}{N-2}$ existiert für hinreichend kleines $\beta$ höchstens eine positive Lösung.

B. Asymptotisches Verhalten nach Exponenten

1. Sublinearer Fall ($0 \le p < 1$): Aufblähen (Blow-up)
   Die Lösung $u_\beta$ konvergiert nicht gegen eine endliche Konstante, sondern bläht sich gleichmäßig im gesamten Gebiet $\Omega$ auf.
   $$u_\beta \sim |\Omega|^{\frac{1}{1-p}} |\partial\Omega|^{-\frac{1}{1-p}} \beta^{-\frac{1}{1-p}} \quad \text{für } \beta \to 0.$$
   Das Wachstum ist also polynomiell in $\beta^{-1}$.

2. Linearer Fall ($p = 1$): Konvergenz gegen eine Konstante
   Im Eigenwertproblem konvergiert die normierte Eigenfunktion $u_\beta$ gleichmäßig gegen die Konstante 1.
   $$u_\beta \to 1 \quad \text{gleichmäßig in } \Omega.$$
   Dies steht im Einklang damit, dass der erste Eigenwert $\lambda_{1,\beta} \to 0$ geht.

3. Superlinearer Fall ($p > 1$): Konvergenz gegen Null
   Die Lösung $u_\beta$ konvergiert gleichmäßig gegen 0.
   $$u_\beta \sim |\Omega|^{-\frac{1}{p-1}} |\partial\Omega|^{\frac{1}{p-1}} \beta^{\frac{1}{p-1}} \quad \text{für } \beta \to 0.$$
   Das Verhalten ist also durch eine Potenz von $\beta$ bestimmt, die gegen 0 geht.

C. Existenz im superkritischen Fall ($p \ge \frac{N+2}{N-2}$)

Für allgemeine Gebiete ist die Existenz von Lösungen im superkritischen Fall ein offenes Problem. Die Autoren beweisen jedoch für den Spezialfall, dass $\Omega$ eine Kugel ist:

• Theorem 1.4: Wenn $\Omega$ eine Kugel ist, $p > \frac{N+2}{N-2}$ und $0 < \beta < \frac{2}{p-1}$, dann existiert eine radiale positive Lösung. Der Beweis nutzt die Existenz von Lösungen im $\mathbb{R}^N$ mit langsamer Abklingrate und passt den Skalierungsfaktor so an, dass die Robin-Randbedingung erfüllt ist.

---

4. Bedeutung und Beitrag

• Vollständige Charakterisierung: Das Paper schließt eine Lücke in der Literatur, indem es das asymptotische Verhalten für alle $p \ge 0$ im Grenzübergang $\beta \to 0$ einheitlich beschreibt. Bisherige Arbeiten konzentrierten sich oft nur auf spezifische Exponentenbereiche oder das Dirichlet/Neumann-Limit.
• Unterscheidung der Regime: Es wird deutlich gezeigt, wie kritisch der Exponent $p$ ist:
  • $p < 1$: Die Nichtlinearität ist schwach genug, dass die Lösung "entartet" und unendlich wird, um die Randbedingung zu erfüllen.
  • $p = 1$: Das System ist im Gleichgewicht (Eigenwertproblem).
  • $p > 1$: Die Nichtlinearität dominiert und zwingt die Lösung, gegen Null zu kollabieren, da die Robin-Bedingung bei kleinem $\beta$ kaum noch eine "Wand" darstellt.
• Neue Existenzresultate: Die Arbeit liefert den ersten Existenzbeweis für radiale Lösungen im superkritischen Fall ($p > 2^*-1$) unter Robin-Randbedingungen auf einer Kugel, ein Bereich, in dem Standard-Variationsmethoden aufgrund des Fehlens der Kompaktheit der Sobolev-Einbettung versagen.
• Technische Innovation: Die Einführung der Hilfsfunktion $v_\beta$ und die genaue Abschätzung des Terms $d_\beta$ erlauben eine präzise Bestimmung der führenden Ordnung (leading-order asymptotics) für alle Fälle, was über qualitative Aussagen (wie "Blow-up" oder "Konvergenz") hinausgeht.

Zusammenfassend liefert das Paper eine tiefgehende Analyse, wie die Wechselwirkung zwischen der Nichtlinearität $u^p$ und der schwachen Robin-Randbedingung ($\beta \to 0$) die Struktur der Lösungen fundamental verändert.