Entire solutions to a strongly competitive nonlinear Schrödinger system

1. Problématique et Contexte

L'article étudie l'existence de solutions entières (définies sur tout l'espace $\mathbb{R}^N$) pour un système d'équations de Schrödinger non linéaires en régime de compétition forte. Le système considéré est :

$$\begin{cases} -\Delta u_1 + u_1 = u_1^p - \Lambda u_1^{a_1} u_2^{a_2} & \text{dans } \mathbb{R}^N \\ -\Delta u_2 + u_2 = u_2^p - \Lambda u_1^{b_1} u_2^{b_2} & \text{dans } \mathbb{R}^N \end{cases}$$

où $N \ge 2$, $p > 1$, les exposants $a_i, b_i > 0$, et $\Lambda > 0$ est un paramètre de compétition.

Contexte scientifique :

• Régime de compétition ($\Lambda \to +\infty$) : Ce régime modélise la ségrégation de phases, où les supports des composantes $u_1$ et $u_2$ deviennent asymptotiquement disjoints. Ce phénomène a été largement étudié pour des solutions radiales ou sur des domaines bornés.
• Structure variationnelle : Sauf pour des choix très spécifiques d'exposants (comme dans le système de Bose-Einstein où $p=3$ et $a_1=b_2=a_2-1=b_1-1=1$), ce système n'admet pas de structure variationnelle, rendant les méthodes variationnelles classiques inapplicables.
• Question ouverte : La question centrale (Q2) est de savoir si l'on peut construire des solutions positives non radiales $(u_1, u_2)$ qui ressemblent, lorsque $\Lambda \to +\infty$, aux parties positive et négative d'une solution signée $W^*$ de l'équation scalaire limite $-\Delta W + W = |W|^{p-1}W$.

2. Méthodologie

Les auteurs utilisent une approche perturbative basée sur la réduction de Lyapunov-Schmidt, combinée à une analyse fine des interactions entre les pics de concentration.

A. Construction de l'Ansatz (Approximation)

Les auteurs construisent une solution approchée $(U_1^*, U_2^*)$ composée de la somme de plusieurs copies de la solution radiale positive unique $U$ de l'équation scalaire $-\Delta U + U = U^p$.

• Configuration géométrique : Les pics de concentration sont disposés selon une structure géométrique spécifique :
  • Deux polygones réguliers concentriques à $k$ côtés.
  • Les pics de $U_1$ sont situés sur les arêtes du polygone extérieur et alternés avec ceux de $U_2$.
  • Les pics de $U_2$ sont situés sur les sommets du polygone intérieur et sur les rayons reliant les sommets des deux polygones.
• Paramètres : La configuration dépend de trois entiers $(m, n, k)$ et de paramètres continus liés aux distances entre les pics ($\ell, \bar{\ell}, \bar{\ell}'$), qui tendent vers l'infini avec $\Lambda$.

B. Conditions d'équilibre (Balancing Conditions)

Pour que l'approximation soit valide, les forces d'interaction entre les pics doivent s'annuler à l'ordre principal. Les auteurs dérivent deux types de conditions :

1. Condition d'équilibre interne : Équilibre des forces entre les pics d'une même composante (gouverné par le principe de superposition).
2. Condition d'équilibre externe (Nouvelle contribution) : Équilibre des forces entre les pics de composantes différentes. Contrairement au cas variationnel, ici les interactions non linéaires ($a_2, b_1 > 1$) ne se décomposent pas simplement par superposition. Les auteurs doivent effectuer des développements asymptotiques délicats pour montrer que les termes dominants de ces interactions sont cohérents avec une superposition formelle, malgré l'absence de linéarité.

Ces conditions permettent de déterminer les relations entre les distances $\ell, \bar{\ell}, \bar{\ell}'$ et le paramètre $\Lambda$.

C. Réduction de Lyapunov-Schmidt

Le problème est réduit à la résolution d'une équation finie-dimensionnelle pour les paramètres de translation des pics ($\alpha_j, \beta_h, \gamma_h$).

1. Problème linéaire : On inverse l'opérateur linéarisé autour de l'ansatz dans un espace de fonctions pondérées, en imposant des conditions d'orthogonalité par rapport au noyau de l'opérateur linéarisé (lié aux invariances par translation).
2. Estimation de l'erreur : On contrôle la norme de l'erreur résiduelle $E$, en tenant compte des termes d'interaction non linéaires complexes.
3. Point fixe : On résout le système non linéaire résultant pour les paramètres de translation en utilisant un argument de point fixe (contraction), montrant qu'il existe une solution exacte proche de l'ansatz pour $\Lambda$ suffisamment grand.

3. Résultats Principaux

Théorème 1.1 :
Sous des hypothèses techniques sur les exposants (sub-critique pour $p$, et conditions spécifiques sur $a_i, b_i$ garantissant que la configuration n'est pas critique, voir hypothèse (1.7)), et pour tout triplet admissible $(m, n, k)$ (où $k$ est suffisamment grand), il existe $\Lambda_0 > 0$ tel que pour tout $\Lambda > \Lambda_0$, le système (1.1) admet une solution positive non radiale $u = (u_1, u_2)$.

• Multiplicité : Puisqu'il existe une infinité de triplets admissibles $(m, n, k)$, le théorème garantit l'existence d'une infinité de solutions non radiales.
• Comportement asymptotique : Lorsque $\Lambda \to +\infty$, les composantes $u_1$ et $u_2$ convergent vers les parties positive et négative d'une solution signée non radiale de l'équation scalaire limite, dont les pics sont distribués selon la configuration géométrique décrite.

4. Contributions Clés et Innovations Techniques

1. Première construction non radiale : C'est la première preuve d'existence de solutions non radiales positives pour des systèmes de Schrödinger fortement compétitifs sur l'espace entier $\mathbb{R}^N$.
2. Gestion des interactions non linéaires : Le défi majeur était l'absence de structure variationnelle et la non-linéarité des termes de couplage ($a_2, b_1 > 1$). Les auteurs ont réussi à développer des expansions asymptotiques précises pour les forces d'interaction inter-composantes, démontrant que le terme dominant se comporte comme si une superposition linéaire était possible, ce qui permet de fermer le schéma de preuve.
3. Configuration géométrique complexe : L'utilisation de deux polygones concentriques avec une alternance de pics offre une flexibilité géométrique inédite pour générer des solutions non radiales, dépassant les configurations radiales ou à symétrie simple.
4. Conditions sur les exposants : Les auteurs identifient des conditions précises (hypothèse 1.7) sous lesquelles la méthode fonctionne. Ils notent que le cas critique du système de Bose-Einstein ($a_1+1=b_2, a_2-1=b_1$) échappe à leur méthode directe car l'approximation n'est pas assez précise sans terme de correction supplémentaire.

5. Signification et Impact

Ce travail ouvre de nouvelles perspectives dans l'étude des systèmes elliptiques non variationnels.

• Il démontre que la ségrégation de phases peut se produire dans des configurations géométriques complexes et non radiales, même sans structure variationnelle.
• Il fournit un cadre robuste pour traiter les interactions fortes dans des systèmes couplés non linéaires, avec des applications potentielles en physique mathématique (condensats de Bose-Einstein, modèles de Lotka-Volterra généralisés).
• La méthode développée pour contrôler les interactions non linéaires sans superposition directe pourrait être adaptée à d'autres problèmes de systèmes compétitifs ou coopératifs dans des régimes limites.

En résumé, cet article résout un problème ouvert majeur en construisant une famille infinie de solutions non radiales pour un système de Schrödinger compétitif, en surmontant les difficultés techniques liées à l'absence de structure variationnelle et à la complexité des interactions non linéaires.