Bifurcations in Isoperimetric Problems with Nonlocal Interactions

1. Problématique et Contexte

L'article étudie un problème variationnel isopérimétrique modélisé sur le modèle de la goutte liquide, mais étendu à des dimensions $N \ge 2$ et incluant des interactions non locales. L'objectif est d'analyser les points critiques de l'énergie fonctionnelle $E(\Omega)$ sous contrainte de volume fixe $|\Omega| = m$.

L'énergie est définie comme la somme de deux termes :
$$E(\Omega) = \text{Per}(\Omega) + D(\Omega)$$
où :

• $\text{Per}(\Omega)$ est le périmètre (surface de la frontière), favorisant la compacité (tension de surface).
• $D(\Omega)$ est un terme d'interaction non locale de type Coulombien :
  • Pour $N \ge 3$ et $0 < \lambda < N-1$ : $D(\Omega) = \frac{1}{2} \iint_{\Omega \times \Omega} \frac{1}{|x-y|^\lambda} dx dy$.
  • Pour $N=2$ (cas logarithmique) : $D(\Omega) = \frac{1}{2\pi} \iint_{\Omega \times \Omega} \log \left(\frac{1}{|x-y|}\right) dx dy$.

Contexte physique :

• En dimension $N=3$ avec $\lambda=1$, ce modèle décrit les noyaux atomiques (protons et neutrons). Le terme de périmètre représente la tension de surface, tandis que le terme non local représente la répulsion coulombienne entre protons.
• Il est bien établi que les boules sont des points critiques naturels. Pour de petits volumes, la boule est le minimiseur global. Pour de grands volumes, la répulsion domine et les boules cessent d'être des minimiseurs globaux (et même locaux au-delà d'un certain seuil).

Question centrale : Au-delà de la simple perte de stabilité des boules, existe-t-il des branches de solutions non sphériques qui bifurquent de la famille des boules ? L'article vise à prouver l'existence d'une suite infinie de rayons critiques où ces bifurcations se produisent.

2. Méthodologie

Les auteurs utilisent une approche analytique rigoureuse basée sur la théorie de la bifurcation et l'analyse fonctionnelle.

1. Paramétrisation des domaines :
   Les domaines $\Omega$ sont supposés étoilés par rapport à l'origine et paramétrés par une fonction $\phi : S^{N-1} \to \mathbb{R}$ :
   $$\Omega_\phi = \{ x \in \mathbb{R}^N : |x| < \phi(x/|x|) \}$$
   On pose $\phi(z) = R + u(z)$, où $R$ est le rayon de la boule de référence et $u$ est une petite perturbation dans l'espace $C^{2,\alpha}(S^{N-1})$.

2. Réduction à une équation sur la sphère :
   En utilisant les conditions d'Euler-Lagrange pour les points critiques sous contrainte de volume, le problème est reformulé comme la recherche de zéros d'un opérateur non linéaire $\Phi(R, u) = 0$ défini sur la sphère unité $S^{N-1}$.

   • La solution triviale est $u \equiv 0$ (la boule).
   • L'objectif est de trouver des solutions non triviales ($u \not\equiv 0$) bifurquant de $(R, 0)$.

3. Théorème de Crandall-Rabinowitz :
   L'outil principal est le théorème de bifurcation de Crandall-Rabinowitz. Pour l'appliquer, les auteurs doivent vérifier :

   • La régularité $C^\infty$ de l'opérateur $\Phi$.
   • Que le noyau (kernel) de l'opérateur linéarisé $D_u \Phi(R, 0)$ est de dimension 1 (ou de dimension finie contrôlée par symétrie).
   • La condition de transversalité : la dérivée par rapport au paramètre de bifurcation (le rayon $R$) de l'opérateur linéarisé appliqué au vecteur propre du noyau ne doit pas appartenir à l'image (range) de l'opérateur.

4. Analyse des Harmoniques Sphériques :
   L'analyse du spectre de l'opérateur linéarisé repose sur la décomposition en harmoniques sphériques. Les auteurs exploitent les propriétés de symétrie des harmoniques sphériques pour réduire l'espace fonctionnel et isoler les modes de bifurcation spécifiques (pairs et impairs).

5. Réduction de Lyapunov-Schmidt :
   Pour montrer l'absence de bifurcation en dehors des rayons critiques, ils utilisent une réduction de Lyapunov-Schmidt pour prouver que toute solution proche d'une boule non critique est nécessairement triviale (uniquement des translations ou des homothéties).

3. Résultats Principaux

Les auteurs établissent l'existence d'une suite non bornée de rayons $(R_n)_{n \ge 2}$ tels que des branches de solutions non sphériques bifurquent de la boule de rayon $R_n$.

A. Cas de la dimension $N=3$ (Modèle de la goutte liquide, $\lambda=1$)

• Rayons de bifurcation : $R_n = \left( \frac{3(2n+1)(n+2)}{8\pi} \right)^{1/3}$ pour $n \ge 2$.
  • Pour $n=2$, on retrouve le résultat connu de Frank ($m=10$).
• Structure des solutions :
  • Si $n$ est pair : Il existe une unique branche de solutions bifurquant de $R_n$. Les solutions sont invariantes sous le groupe $O(2) \times O(1)$ (symétrie axiale). La forme du domaine passe de prolate à oblate selon que le rayon est inférieur ou supérieur à $R_n$. La bifurcation est transcritique (le volume change à l'ordre linéaire).
  • Si $n$ est impair : Il existe plusieurs branches de solutions (pour $m = \frac{n+1}{2}, \dots, n$). Les solutions possèdent des symétries plus complexes (invariance sous des rotations d'angle $2\pi/m$). La bifurcation est sous-critique ou supercritique (le volume est conservé à l'ordre linéaire).
• Unicité locale : Pour tout $R \neq R_n$, il n'existe pas de solutions non triviales proches de la boule.

B. Cas de la dimension $N=2$ (Interaction logarithmique)

• Rayons de bifurcation : $R_n = (2n(n+1))^{1/3}$ pour $n \ge 2$.
• Structure des solutions :
  • Pour chaque $n$, il existe une branche de solutions bifurquant de $R_n$.
  • Les solutions sont de la forme $u_t(\phi) \approx t \sin(n\phi)$.
  • Contrairement au cas $N=3$ pair, la bifurcation préserve le volume à l'ordre linéaire (et probablement à l'ordre quadratique), quelle que soit la parité de $n$.
  • Le nombre de branches de solutions ne croît pas avec $n$ (dimension du noyau constante).

C. Cas général $N \ge 4$ (Interaction avec exposant $\lambda$)

• Rayons de bifurcation : Une formule explicite dépendant de $N$ et $\lambda$ est donnée (équation 1.11).
• Complexité accrue :
  • Pour $n$ pair, il existe des branches de solutions invariantes sous $O(k) \times O(N-k)$ pour $k=1, \dots, N-1$. Cela révèle une structure de bifurcation beaucoup plus riche que dans les dimensions inférieures, avec des types de symétrie géométrique supplémentaires.
  • Pour $n$ impair, des branches avec des symétries spécifiques (liées aux harmoniques sphériques de type "produit") sont également construites.

4. Contributions Clés et Signification

1. Généralisation du résultat de Frank : L'article étend le résultat de Frank (2014) qui ne traitait que du premier point de bifurcation ($n=2$) en dimension 3. Les auteurs prouvent l'existence d'une suite infinie de points de bifurcation pour toutes les dimensions et types d'interactions considérés.
2. Classification des symétries : L'article fournit une description précise des symétries des solutions bifurquées, distinguant nettement les cas pairs et impairs, et montrant comment la dimension de l'espace ambiant enrichit la structure des solutions (notamment pour $N \ge 4$).
3. Comportement du volume : L'analyse fine du développement asymptotique des solutions révèle des différences fondamentales entre les cas pairs et impairs concernant la variation de volume au voisinage du point de bifurcation (bifurcation transcritique vs sous/supercritique).
4. Rigueur mathématique : La preuve combine des techniques avancées d'analyse non linéaire (théorème de bifurcation, régularité des opérateurs intégraux, théorie spectrale des harmoniques sphériques) pour établir des résultats d'existence et d'unicité locale.

Conclusion :
Ce travail démontre que la structure des équilibres dans les problèmes isopérimétriques avec interactions non locales est bien plus riche que la simple compétition entre sphères et formes toroïdales. Il établit l'existence d'une hiérarchie infinie de formes d'équilibre non sphériques, apparaissant à des échelles de volume spécifiques, ce qui a des implications potentielles pour la compréhension de la fission nucléaire et de la formation de motifs dans les systèmes biologiques ou physiques.