The nonlocal attraction-repulsion transport equation with power kernels

Résumé Technique : L'équation de transport à attraction-répulsion non locale avec noyaux de puissance

Énoncé du problème
Cet article étudie l'équation de continuité non locale sur $\mathbb{R}^d$ décrivant l'évolution d'une densité de probabilité $\phi_t$ régie par la compétition entre une attraction vers une mesure de fond prescrite $\omega$ et une auto-répulsion entre les particules. La dynamique est gouvernée par l'équation :
$$\partial_t \phi_t = \text{div}\left(\phi_t (\nabla V - K_r * \phi_t)\right), \quad t > 0,$$
où le potentiel de fond est $V = \psi_a * \omega$, et les noyaux d'interaction sont des lois de puissance définies par $\psi_s(x) = |x|^{1+s}$ et $K_s(x) = \nabla \psi_s(x) = (1+s)|x|^{s-1}x$. Les exposants $a$ (attraction) et $r$ (répulsion) appartiennent à l'intervalle $[0, 1)$.

L'étude se concentre sur le régime « discordant » où $a \neq r$, ainsi que sur le cas à exposants égaux $a=r$ avec une masse $\omega(\mathbb{R}^d)$ arbitraire. L'équation est formellement le flot de gradient de Wasserstein de l'énergie d'interaction :
$$E(\phi) = \int_{\mathbb{R}^d} (\psi_a * \omega) d\phi - \frac{1}{2} \iint_{\mathbb{R}^d \times \mathbb{R}^d} \psi_r(x-y) d\phi(x) d\phi(y).$$
Lorsque $a=r$ et que $\omega$ est une mesure de probabilité, cette énergie correspond à la Discrépance Moyenne Maximale (MMD) avec un noyau de distance négative. Les principaux défis mathématiques proviennent de la singularité des noyaux à l'origine (pour $s=0$) et de l'absence de bornes lipschitziennes globales, ce qui empêche l'application directe des théories de flot de gradient standard nécessitant la $\lambda$-convexité.

Méthodologie
Les auteurs emploient une approche analytique multi-facettes :

1. Régularisation par le rayon au carré : Pour traiter les singularités et le manque de continuité lipschitzienne, les auteurs introduisent une régularisation $\psi_{s,\epsilon}(x) = (|x|^2 + \epsilon)^{(1+s)/2}$. Contrairement à la mollification par convolution, celle-ci préserve exactement la convexité et les estimations de signe unidirectionnelles. Ils établissent la bien-posée globale pour le système régularisé et passent à la limite $\epsilon \to 0$ en utilisant des arguments de compacité.
2. Cadre Lagrangien : Les solutions sont construites et analysées comme des images directes (push-forwards) des données initiales par une application de flux $X_t$, satisfaisant l'EDO caractéristique $\dot{X}_t = -\nabla V(X_t) + K_r * \phi_t(X_t)$.
3. Caractérisation par frontière libre : Pour les états stationnaires, les auteurs utilisent l'analyse de Fourier et la théorie des distributions homogènes. Ils inversent les opérateurs de convolution en utilisant les laplaciens fractionnaires pour réduire la condition de stationnarité $\phi \nabla F_\phi = 0$ à un problème de frontière libre impliquant un problème de Dirichlet extérieur fractionnaire.
4. Estimations d'énergie et de moments : Des bornes $L^\infty$ uniformes et des estimations de moments sont dérivées via des arguments de type Gronwall sur le système régularisé. Celles-ci sont cruciales pour prouver la compacité des solutions et le confinement du support dans les régimes attractifs.

Contributions clés et résultats

• Bien-posée globale : L'article étabète la bien-posée globale lagrangienne pour des données initiales arbitraires dans $L^1 \cap L^\infty$ avec des moments finis. La solution est unique dans la classe des solutions lagrangiennes distributionnelles satisfaisant des bornes $L^\infty$ et de moments uniformes. De plus, les auteurs prouvent la propagation de la régularité $W^{n,\infty}$ pour les données initiales lisses.
• Compacité uniforme du support : Dans les régimes dominés par l'attraction ($a > r$, ou $a=r$ avec $\omega(\mathbb{R}^d) > 1$), les auteurs prouvent que si les données initiales sont à support compact, le support de la solution reste uniformément borné pour tout temps. Cela est obtenu en montrant que le diamètre du support ne peut croître indéfiniment en raison de la dominance de la force attractive sur la répulsion aux grandes distances. Un contre-exemple est fourni montrant que cette propriété échoue si $a=r > 1$.
• Caractérisation des états stationnaires : Pour la plage non-quadratique dominée par l'attraction ($0 \le r \le a < 1$), l'article caractérise les états stationnaires à flux nul comme des solutions d'un problème de frontière libre. Plus précisément, la densité stationnaire $\phi$ est exprimée par $\phi = Q_{a,r} - L_r F$, où $L_r$ est un opérateur de convolution inverse lié au laplacien fractionnaire, et $F$ satisfait un problème de Dirichlet extérieur fractionnaire.
• Exemples stationnaires explicites : En utilisant la caractérisation dérivée, les auteurs construisent des mesures stationnaires explicites pour des géométries spécifiques :
  • 1D : Solutions d'intervalles impliquant des potentiels de Riesz.
  • 2D : Solutions de disques utilisant le noyau de Green Boggio-Kelvin pour le laplacien fractionnaire.
  • 3D : Solutions pour des fonds de type boule uniforme et gaussien, où la caractérisation se réduit à des opérateurs différentiels locaux (Laplacien) en raison des exposants spécifiques ($a=r=0$).
• Convergence à long terme : L'article prouve que toute solution globale possédant une énergie bornée et des bornes de moments uniformes possède une sous-suite convergeant vers un état stationnaire à flux nul. Si l'ensemble $\omega$-limite contient un état stationnaire unique, la convergence totale est établie.

Signification et affirmations
Les auteurs positionnent ce travail comme une extension des résultats unidimensionnels et à exposants égaux vers des dimensions arbitraires $d \ge 1$ et le régime discordant $a \neq r$.

• Avancement théorique : L'article fournit un cadre rigoureux pour analyser les flux d'attraction-répulsion singuliers où l'énergie n'est pas $\lambda$-convexe, une condition qui empêche généralement les théories d'existence de flot de gradient standard. L'utilisation de la régularisation par le rayon au carré est mise en évidence comme une innovation technique clé préservant les propriétés structurelles nécessaires.
• Aperçu de modélisation : Les résultats offrent un prototype de modèle pour l'« accès inégal » dans la distribution des ressources. Dans le cas dominé par l'attraction, la population d'équilibre $\phi_\infty$ est attirée vers la distribution de ressources $\omega$ mais confinée par la répulsion interne, ce qui résulte en un état stationnaire qui ne coïncide pas avec $\omega$. La caractérisation explicite de $\phi_\infty$ permet de quantifier l'inégalité dans la distribution de la richesse ou des ressources.
• Lien avec la MMD : L'article clarifie la relation entre ces équations de transport et la Discrépance Moyenne Maximale (MMD) utilisée en apprentissage automatique. Alors que des travaux précédents traitaient le flot de gradient de la MMD pour des noyaux lisses ou des cas spécifiques, cet article étend l'analyse au noyau de distance négative singulier dans des dimensions supérieures, établissant la bien-posée et le comportement asymptotique pour le cas $a=r$ avec $\omega(\mathbb{R}^d) \ge 1$.
• Validation numérique : Les profils stationnaires théoriques dérivés de la caractérisation de la frontière libre sont vérifiés par rapport à des simulations de particules numériques, mont'rant un accord en une et deux dimensions.

Les auteurs notent explicitement que leurs résultats sur le confinement du support uniforme ne s'appliquent pas au cas $a=r$ avec $\omega(\mathbb{R}^d)=1$ si $\omega$ n'est pas à support compact, car le support est censé croître pour correspondre à $\omega$. Ils reconnaissent également que le cas $a=r$ avec $\omega(\mathbb{R}^d)=1$ est traité dans un travail parallèle de Rosenzweig, Slepčev et Wang, se concentrant sur les taux de convergence de champ moyen.