Segregated solutions for a critical Choquard system with a small interspecies repulsive force

1. Problématique et Contexte

L'article étudie un système d'équations elliptiques couplées avec des non-linéarités non locales de type Choquard, défini dans l'espace $\mathbb{R}^4$. Le système considéré est :

$$\begin{cases} -\Delta u = (|x|^{-\mu} * |u|^{\frac{4-\mu}{2}})|u|^{\frac{2-\mu}{2}}u + \beta uv^2 & \text{dans } \mathbb{R}^4 \\ -\Delta v = (|x|^{-\mu} * |v|^{\frac{4-\mu}{2}})|v|^{\frac{2-\mu}{2}}v + \beta u^2v & \text{dans } \mathbb{R}^4 \end{cases}$$

où :

• $*$ désigne la convolution standard.
• $0 < \mu \le 4$.
• $\beta < 0$ représente un paramètre de couplage répulsif (régime compétitif) entre les deux composantes $u$ et $v$.
• Le système est critique au sens de l'inégalité de Hardy-Littlewood-Sobolev (HLS). L'exposant critique supérieur est $2^*_\mu = \frac{2N-\mu}{N-2}$, qui pour $N=4$ correspond à la puissance non linéaire présente dans le terme de Choquard.

Contexte physique et mathématique :
Les équations de type Choquard (ou Schrödinger-Newton) modélisent des phénomènes physiques tels que les polarons, les plasmas à un composant, ou les étoiles à bosons. La difficulté majeure réside dans la nature non locale du terme d'interaction et le caractère critique du système, qui empêche la compacité des opérateurs d'immersion de Sobolev. L'objectif est de construire des solutions où les composantes se séparent (ségrégation) en présence d'une force répulsive faible ($|\beta|$ petit).

2. Méthodologie

L'auteur utilise la méthode de réduction de Lyapunov-Schmidt (ou méthode de réduction finie dimensionnelle) pour construire des solutions par perturbation. La démarche se décompose en plusieurs étapes clés :

A. Ansatz (Approximation initiale)

La solution recherchée $(u, v)$ est approchée par la somme de solutions fondamentales :

• La première composante $u$ est proche de la solution radiale positive $U$ de l'équation de Choquard critique isolée.
• La deuxième composante $v$ est approchée par une superposition de $k$ "bubbles" (solutions concentrées) $U_{\lambda, \xi_j}$, centrées aux sommets $\xi_j$ d'un polygone régulier dans le plan $(x_1, x_2)$.
  $$u \approx U + \phi_1, \quad v \approx \sum_{j=1}^k U_{\lambda, \xi_j} + \phi_2$$
  où $\lambda$ est un paramètre d'échelle très grand (concentration) et $\phi_1, \phi_2$ sont des termes de reste petits.

B. Espace fonctionnel et Symétries

Pour gérer la non-compacité et les interactions entre les pics, l'auteur travaille dans des espaces pondérés $X$ définis par des normes spécifiques ($\|\cdot\|_*$ et $\|\cdot\|_{**}$). Ces espaces imposent des symétries :

• Parité par rapport aux variables $x_2, x_3, x_4$.
• Invariance par rotation de $2\pi/k$ dans le plan $(x_1, x_2)$.
• Invariance par transformation de Kelvin (cruciale pour le cas critique en dimension 4).

C. Réduction Linéaire et Inversibilité

Le problème est reformulé en une équation linéaire perturbée :
$$L(\phi_1, \phi_2) = E + N(\phi_1, \phi_2)$$
où $L$ est l'opérateur linéarisé autour de l'ansatz, $E$ est le terme d'erreur (résidu de l'approximation) et $N$ les termes non linéaires quadratiques et supérieurs.

• L'auteur démontre l'inversibilité de l'opérateur linéaire $L$ restreint au sous-espace orthogonal au noyau de l'opérateur linéarisé de l'équation limite.
• Une estimation précise du terme d'erreur $E$ est obtenue, montrant qu'il est de l'ordre de $O(|\beta|/\lambda)$.

D. Problème Réduit (Finite Dimensional Reduction)

La résolution du système complet se ramène à la résolution d'une équation scalaire pour le paramètre $\lambda$ (la taille des pics). En projetant l'équation sur le noyau de l'opérateur linéarisé, on obtient une fonctionnelle réduite $c(\lambda)$.
L'objectif est de trouver $\lambda$ tel que $c(\lambda) = 0$. L'analyse asymptotique de $c(\lambda)$ révèle une compétition entre deux termes :

1. Un terme provenant de l'interaction entre les $k$ pics de la composante $v$ (terme de répulsion géométrique).
2. Un terme provenant du couplage $\beta$ entre $u$ et $v$.

3. Résultats Principaux

Le théorème principal (Théorème 1.2) établit l'existence de solutions pour tout $k \ge 2$ et pour $\beta$ négatif avec $|\beta|$ suffisamment petit.

Énoncé du résultat :
Il existe $\beta_0 < 0$ tel que pour tout $\beta \in [\beta_0, 0)$, le système admet une solution $(u, v)$ de la forme :
$$u = U + \phi_1, \quad v = \sum_{j=1}^k U_{\lambda, \xi_j} + \phi_2$$
avec les propriétés suivantes :

• Les erreurs $\phi_1, \phi_2$ tendent vers zéro dans la norme appropriée.
• Le paramètre d'échelle $\lambda$ satisfait une relation asymptotique précise :
  $$\frac{1}{\lambda} = O\left( e^{-\frac{c}{|\beta|}} \right)$$
  Cela signifie que lorsque la force répulsive $|\beta|$ tend vers 0, les pics de la composante $v$ deviennent extrêmement concentrés ($\lambda \to \infty$).

Structure de la solution :

• La composante $u$ reste globale et ressemble à la solution fondamentale de l'équation de Choquard seule.
• La composante $v$ présente un comportement de "blow-up" (explosion) en $k$ points distincts formant un polygone régulier.
• Les deux composantes sont ségrégées : là où $v$ est très concentrée (les sommets du polygone), $u$ est très faible, et vice-versa, en raison de la répulsion $\beta < 0$.

4. Contributions Clés et Signification

1. Première étude de ce type : À la connaissance de l'auteur, c'est le premier travail à combiner la théorie des systèmes compétitifs (ségrégation) avec l'équation de Choquard critique. La littérature précédente traitait soit des systèmes de Schrödinger couplés, soit des équations de Choquard isolées.
2. Gestion de la criticité et de la non-localité : L'article surmonte les difficultés techniques majeures liées à la combinaison de l'exposant critique de Sobolev-HLS et de la non-localité. L'utilisation de la transformation de Kelvin et des espaces pondérés adaptés est essentielle pour obtenir les estimées nécessaires.
3. Rôle du paramètre $\beta$ : L'article démontre que même une très faible force répulsive suffit à induire une ségrégation complète et une concentration extrême (blow-up) de l'une des composantes, suivant une loi exponentielle en fonction de $1/|\beta|$.
4. Limitation dimensionnelle : L'auteur note que la méthode ne fonctionne que pour $N=4$. Pour $N \ge 5$, la croissance du terme de couplage devient sous-linéaire par rapport à la croissance critique, ce qui empêche l'application de la méthode de réduction standard utilisée ici.

Conclusion

Ce travail apporte une avancée significative dans l'analyse des systèmes elliptiques critiques non locaux. Il prouve que la structure géométrique des solutions (polygones réguliers) et le phénomène de ségrégation persistent même en présence de termes non locaux complexes, offrant ainsi un nouveau cadre pour l'étude des interactions compétitives dans les modèles de physique mathématique (comme les étoiles à bosons ou les plasmas).