The Rabinowitz continuum of subcritical Gelfand problems and free boundary-type equations arising in plasma physics

1. Problématique et Contexte

L'article s'intéresse au comportement qualitatif du continuum non borné de Rabinowitz associé au problème de Gelfand sous-critique :
$$\begin{cases} -\Delta v = \mu(1 + v)^p & \text{dans } \Omega \\ v > 0 & \text{dans } \Omega \\ v = 0 & \text{sur } \partial\Omega \end{cases} \quad (Q)_\mu$$
où $\Omega \subset \mathbb{R}^N$ ($N \ge 2$) est un domaine borné de classe $C^3$, et $1 < p < p_N$ (avec $p_N = \infty$ si $N=2$ et $p_N = \frac{N}{N-2}$ si $N \ge 3$).

État de l'art et lacunes :

• Il est bien connu qu'il existe un paramètre critique $\mu^*$ tel que pour $\mu < \mu^*$, il existe au moins deux solutions (une solution minimale et une solution de type "mountain-pass"), une solution unique pour $\mu = \mu^*$, et aucune pour $\mu > \mu^*$.
• La branche de solutions minimales est lisse et monotone en fonction de $\mu$.
• Cependant, une description précise et globale de la connectivité et du comportement qualitatif des solutions non minimales (en particulier la monotonie de la norme $L^\infty$ le long de la branche globale) est absente pour des domaines généraux, sauf dans le cas spécifique des boules (où les solutions sont radiales) ou de certains domaines convexes en dimension 2.
• L'objectif est de développer une méthode générale pour décrire ce continuum, y compris les solutions non minimales, au-delà des résultats classiques limités aux solutions minimales.

2. Méthodologie : Une Approche par Problèmes à Frontière Libre

Les auteurs proposent une nouvelle approche en reliant le problème de Gelfand $(Q)_\mu$ à un problème contraint de type frontière libre $(P)_\lambda$, qui apparaît en physique des plasmas (équations de type Grad-Shafranov) :

$$\begin{cases} -\Delta \psi = [\alpha + \lambda \psi]_+^p & \text{dans } \Omega \\ \int_\Omega [\alpha + \lambda \psi]_+^p = 1 \\ \psi = 0 & \text{sur } \partial\Omega \end{cases} \quad (P)_\lambda$$

L'idée centrale :

1. Paramétrisation : Si $(\alpha_\lambda, \psi_\lambda)$ résout $(P)_\lambda$ avec $\alpha_\lambda > 0$, alors $v(\lambda) = \frac{\lambda}{\alpha_\lambda}\psi_\lambda$ résout $(Q)_\mu$ avec $\mu = \mu_\lambda = \lambda \alpha_\lambda^{p-1}$.
2. Remplacement de la norme : Au lieu d'étudier la monotonie de la norme $L^\infty$ de la solution (qui est difficile à contrôler pour les solutions non minimales), les auteurs utilisent l'énergie du plasma $E_\lambda = \frac{1}{2} \int_\Omega |\nabla \psi_\lambda|^2$ comme paramètre de suivi.
3. Unicité et Stabilité : Le cœur de la méthode repose sur la résolution d'un problème ouvert concernant l'unicité des solutions de $(P)_\lambda$ et l'analyse spectrale de l'opérateur linéarisé associé.

3. Contributions Clés et Résultats Principaux

A. Résolution d'un problème ouvert sur l'unicité (Théorème 1.1)

Les auteurs résolvent un problème ouvert de longue date concernant l'unicité des solutions des équations de type Grad-Shafranov.

• Ils introduisent une première valeur propre pondérée $\sigma_1(\alpha_\lambda, \psi_\lambda)$ associée à l'opérateur linéarisé non-local de $(P)_\lambda$.
• Ils définissent $\lambda^*(\Omega, p)$ comme le supremum des $\lambda$ pour lesquels $\sigma_1 > 0$ (solutions "stables").
• Résultat : Pour tout $\lambda \in [0, \lambda^*(\Omega, p))$, il existe une solution unique de $(P)_\lambda$. Cette solution forme une courbe $C^1$ (voire analytique) $G^*(\Omega)$ où $\alpha_\lambda$ est strictement décroissant et l'énergie $E_\lambda$ est strictement croissante.
• Cela établit que les solutions positives connues précédemment sont en fait les seules solutions dans cet intervalle.

B. Paramétrisation Globale du Continuum de Rabinowitz (Théorème 1.2)

En exploitant l'unicité et la stabilité, les auteurs construisent une paramétrisation globale du continuum de solutions de $(Q)_\mu$.

• La branche $R_+ = \{(\mu_\lambda, v(\lambda))\}$ est une courbe analytique.
• Monotonie de l'énergie : L'énergie $E_\lambda$ est strictement croissante le long de toute la branche, y compris pour les solutions non minimales. Cela généralise la propriété de monotonie classique des solutions minimales.
• Comportement asymptotique : Lorsque $\lambda$ approche une valeur critique $\lambda_+$ (où $\alpha_\lambda \to 0^+$), on observe que $\mu_\lambda \to 0^+$ et $\|v(\lambda)\|_\infty \to +\infty$.

C. Estimations Spectrales et Inégalités (Théorème 1.4)

Une question ouverte majeure était de savoir si $\lambda^*(\Omega, p) > \lambda_+(\Omega, p)$ (où $\lambda_+$ est la valeur où $\alpha_\lambda$ s'annule).

• Les auteurs prouvent que pour une boule unité $D_N$, on a $\lambda^*(D_N, p) \ge \lambda_+(D_N, p)$.
• Cette inégalité est obtenue par une preuve par contradiction utilisant des estimées de régularité fines de $\alpha_\lambda$ et un lemme de pliage (bending lemma) de type Crandall-Rabinowitz.
• Cela garantit que la branche unique de solutions stables couvre tout le continuum de solutions positives jusqu'à la singularité.

D. Profil en Cloche (Bell-Shaped) et Comportement de $\mu$ (Théorème 1.5)

Pour le cas de la boule unité $D_N$, les auteurs caractérisent complètement le comportement de $\mu_\lambda$ en fonction de l'énergie.

• La courbe $(\mu_\lambda, E_\lambda)$ présente un profil en cloche (bell-shaped).
• Il existe un unique point de retournement $\lambda_t$ tel que :
  • $\frac{d\mu_\lambda}{d\lambda} > 0$ pour $\lambda < \lambda_t$.
  • $\frac{d\mu_\lambda}{d\lambda} < 0$ pour $\lambda > \lambda_t$.
• Cela confirme que la branche globale part de $(0,0)$, atteint un $\mu_{max}$, puis redescend vers $\mu=0$ tout en voyant l'énergie et la norme $L^\infty$ exploser. Ce résultat généralise les connaissances classiques (limitées aux solutions minimales ou aux cas radiaux) à l'ensemble du continuum incluant les solutions non minimales.

4. Signification et Impact

1. Nouvelle Paramétrisation Globale : L'article fournit la première description globale et monotone (en termes d'énergie) du continuum de Rabinowitz pour des domaines généraux, dépassant la simple comptabilité des solutions.
2. Unicité des Équations de Grad-Shafranov : La preuve de l'unicité des solutions pour $(P)_\lambda$ dans la plage sous-critique résout un problème mathématique indépendant important en physique des plasmas.
3. Extension aux Solutions Non Minimales : La méthode permet de suivre les solutions non minimales (issues du théorème du point selle) avec la même rigueur que les solutions minimales, en utilisant la stabilité spectrale pondérée comme critère de sélection.
4. Généralisation du Profil en Cloche : Bien que le profil en cloche soit connu pour les boules, cette approche offre une nouvelle perspective basée sur l'énergie et la contrainte, suggérant des voies pour étendre ces résultats à des domaines convexes non radiaux (notamment en dimension 2).

En résumé, ce travail établit un pont rigoureux entre les problèmes de Gelfand et les équations de frontière libre de la physique des plasmas, offrant un outil puissant pour analyser la structure globale des solutions non linéaires elliptiques.