Infinite-order multisoliton solutions to the Benjamin--Ono equation and soliton resolution

1. Problématique et Contexte

L'article s'intéresse au comportement asymptotique en temps long des solutions de l'équation de Benjamin–Ono (BO) sur la droite réelle :
$$\partial_t u - \partial_x |D_x|u + \partial_x(u^2) = 0, \quad u(0, x) = u_0(x).$$
Cette équation, intégrable, modélise les ondes de gravité internes unidirectionnelles dans un fluide à deux couches de profondeur infinie.

Contexte antérieur :
Il a été démontré récemment (notamment dans [GGM26]) que pour des données initiales $u_0$ suffisamment décroissantes à l'infini (en particulier $x u_0 \in L^2(\mathbb{R})$), la solution se décompose asymptotiquement en une superposition finie de solitons (ondes solitaires) plus un terme de rayonnement qui tend vers zéro. Le nombre de solitons est fini grâce à la décroissance spatiale de la donnée initiale.

Objectif de l'article :
Les auteurs visent à étendre ce résultat de « résolution en solitons » à une classe de données initiales où le nombre de solitons est infini. Ces données, appelées « multisolitons d'ordre infini », appartiennent à $L^2(\mathbb{R})$ mais ne satisfont pas la condition de décroissance rapide ($x u_0 \notin L^2(\mathbb{R})$). Le défi principal est de prouver que, malgré l'absence de décroissance rapide et la présence d'une infinité de modes, la solution se décompose toujours en une superposition infinie de solitons indépendants sans terme de rayonnement résiduel significatif.

2. Méthodologie

La preuve repose sur une combinaison d'analyse spectrale, de théorie des opérateurs et d'une formule explicite récente pour l'équation de Benjamin–Ono.

A. Données initiales et Définition

Les auteurs définissent une classe de données initiales $u_0$ comme une somme infinie de profils de solitons élémentaires $R_p$ :
$$u_0(x) = \sum_{j=1}^\infty \frac{2\text{Im}(p_j)}{|x + p_j|^2},$$
où $(p_j)_{j \ge 1}$ est une suite de points dans le demi-plan supérieur $\mathbb{C}^+$. La condition de convergence imposée est :
$$\sum_{j,k} \frac{\text{Im}(p_j)}{|p_j - p_k|^2} < \infty.$$
Cette condition assure que $u_0 \in L^2(\mathbb{R})$ mais implique que $u_0 \notin L^1(\mathbb{R})$ et $x u_0 \notin L^2(\mathbb{R})$.

B. Outils Spectraux et Opérateurs de Lax

L'approche utilise la structure de paire de Lax de l'équation. L'opérateur de Lax associé à $u_0$ est défini sur l'espace de Hardy $L^2_+(\mathbb{R})$ par :
$$L_{u_0} = -i \frac{d}{dx} - T_{u_0},$$
où $T_{u_0}$ est l'opérateur de Toeplitz.
Un résultat clé de la section 3 est la démonstration que les fonctions propres $\phi_j$ de $L_{u_0}$ (associées aux valeurs propres négatives $\lambda_j$) appartiennent au domaine de l'opérateur $X^*$ (adjoint de la multiplication par $x$), ce qui est crucial pour manipuler les termes de dérive dans la formule explicite.

C. Formule Explicite

Les auteurs utilisent une formule explicite dérivée par P. Gérard ([Gé23]) reliant la solution $u(t,x)$ à l'opérateur de Lax :
$$\Pi u(t, z) = \frac{1}{2i\pi} I_+ \left( (X^* - 2t L_{u_0} - z)^{-1} \Pi u_0 \right),$$
où $\Pi$ est le projecteur de Riesz sur $L^2_+$, $z \in \mathbb{C}^+$, et $I_+$ désigne la limite de la transformée de Fourier en $0^+$. Cette formule permet d'analyser la dynamique temporelle via le spectre de $L_{u_0}$.

D. Stratégie de Preuve

1. Construction spectrale : Ils montrent que pour ces données initiales, l'opérateur $L_{u_0}$ possède une infinité de valeurs propres négatives $(\lambda_j)_{j \ge 1}$.
2. Formule de Plancherel déformée : Ils établissent une identité liant la norme $L^2$ de la donnée initiale à la somme des valeurs propres : $\sum 2\pi|\lambda_k| = \|\Pi u_0\|_{L^2}^2$.
3. Analyse asymptotique : En utilisant la formule explicite et des estimations de résolvantes, ils prouvent que le terme associé à chaque valeur propre $\lambda_j$ converge vers un soliton se déplaçant à la vitesse $c_j = 1/\text{Im}(p_j^\infty) = 2|\lambda_j|$.
4. Contrôle du reste : Ils démontrent que l'erreur entre la solution exacte et la somme tronquée des $N$ premiers solitons tend vers zéro lorsque $N \to \infty$ et $t \to \infty$.

3. Résultats Principaux

Le théorème principal (Théorème 1.2) établit la résolution en solitons pour ces données d'ordre infini :

1. Décomposition asymptotique : Pour toute donnée initiale $u_0$ de la forme (1.3) satisfaisant la condition (1.4), il existe une suite de paramètres $(p_j^\infty)_{j \ge 1}$ (avec des parties imaginaires strictement croissantes) telle que la solution $u(t, \cdot)$ se décompose asymptotiquement en une somme infinie de solitons :
   $$\lim_{N \to \infty} \limsup_{t \to \infty} \left\| u(t, \cdot) - \sum_{j=1}^N R_{p_j^\infty}(\cdot - c_{p_j^\infty} t) \right\|_{L^2} = 0.$$
2. Absence de rayonnement : Contrairement aux cas génériques où un terme de rayonnement persiste, ici le terme de reste tend vers zéro dans des normes plus fortes ($L^\infty$ et $\dot{H}^s$ pour $s \ge 1/2$) lorsque la somme est prise sur l'ensemble infini des solitons.
3. Régularité : La solution $u(t, \cdot)$ appartient à tous les espaces de Sobolev $H^s(\mathbb{R})$ pour tout $s \in \mathbb{R}$, et la somme infinie de solitons $u_{sol}(t, x)$ définit une fonction continue tendant vers 0 à l'infini.
4. Carte spectrale : L'article définit une application entre la suite des paramètres initiaux $(p_j)$ et la suite des paramètres asymptotiques $(p_j^\infty)$. La caractérisation exacte de l'image de cette application reste un problème ouvert (lié à un problème inverse spectral).

4. Contributions Clés

• Extension à l'infini : C'est la première construction rigoureuse de solutions multisolitons d'ordre infini pour l'équation de Benjamin–Ono sur la ligne, dépassant le cadre des données à décroissance rapide ($x u_0 \in L^2$).
• Nouvelle classe de données : Introduction d'une classe de données initiales dans $L^2$ mais non dans $L^1$ (ni à moment fini), élargissant considérablement le domaine de validité de la théorie de la résolution en solitons.
• Propriétés des fonctions propres : Preuve que les fonctions propres de l'opérateur de Lax associé à ces potentiels singuliers appartiennent au domaine de l'opérateur $X^*$, une étape technique nécessaire pour utiliser la formule explicite.
• Formule de Plancherel pour l'ordre infini : Établissement d'une relation de type Plancherel reliant la norme $L^2$ de la somme infinie de solitons à la somme des valeurs propres de l'opérateur de Lax.

5. Signification et Impact

Ce travail est significatif pour plusieurs raisons :

• Théorie des systèmes intégrables : Il démontre la robustesse du phénomène de résolution en solitons, même dans des régimes où la densité d'énergie est faible à l'infini et où le nombre de modes est infini.
• Problèmes inverses : Il soulève des questions profondes sur les problèmes inverses spectraux pour les opérateurs de Lax associés à des potentiels non décroissants, reliant la dynamique non linéaire à la théorie spectrale.
• Méthodologie : L'utilisation de la formule explicite de Gérard combinée à l'analyse spectrale fine offre un modèle pour l'étude d'autres équations intégrables avec des données initiales peu régulières ou à décroissance lente.

En résumé, Gassot et Gérard réussissent à prouver que la structure de soliton est intrinsèque à la dynamique de l'équation de Benjamin–Ono, même pour des configurations initiales complexes et infiniment riches, confirmant que la radiation est un phénomène transitoire qui disparaît complètement au profit d'une superposition infinie d'ondes solitaires.