2025-11-23T09:16:16.957000

Reversibility, covariance and coarse-graining for Langevin dynamics: On the choice of multiplicative noise

Ayala, Dirr, Pavliotis et al.
We study the interplay between reversibility, geometry, and the choice of multiplicative noise (in particular Itô, Stratonovich, Klimontovich) in stochastic differential equations (SDEs). Building on a unified geometric framework, we derive algebraic conditions under which a diffusion process is reversible with respect to a Gibbs measure on a Riemannian manifold. The condition depends continuously on a parameter $λ\in [0,1]$ which interpolates between the conventions of Itô ($λ= 0$), Stratonovich ($λ= \frac 1 2$) and Klimontovich ($λ= 1$). For reversible slow-fast systems of SDEs with a block-diagonal diffusion structure, we show, using the theory of Dirichlet forms, that both reversibility and the Klimontovich noise interpretation are preserved under coarse-graining. In particular, we prove that the effective dynamics for the slow variables, obtained via projection onto a lower-dimensional manifold, retain the Klimontovich interpretation and remain reversible with respect to the marginal Gibbs measure/free energy. Our results provide a flexible variational framework for modeling coarse-grained reversible dynamics with nontrivial geometric and noise structures.
academic

Réversibilité, covariance et renormalisation pour la dynamique de Langevin : Sur le choix du bruit multiplicatif

Informations fondamentales

  • ID de l'article: 2511.03347
  • Titre: Reversibility, covariance and coarse-graining for Langevin dynamics: On the choice of multiplicative noise
  • Auteurs: Mario Ayala (TU München), Nicolas Dirr (Cardiff University), Grigorios A. Pavliotis (Imperial College London), Johannes Zimmer (TU München)
  • Classification: math.PR (Théorie des probabilités)
  • Date de publication: 6 novembre 2025 (prépublication arXiv)
  • Lien de l'article: https://arxiv.org/abs/2511.03347

Résumé

Cet article étudie l'interaction entre la réversibilité, les structures géométriques et le choix du bruit multiplicatif (en particulier les interprétations d'Itô, Stratonovich et Klimontovich) dans les équations différentielles stochastiques (EDS). Sur la base d'un cadre géométrique unifié, les auteurs déduisent les conditions algébriques pour que les processus de diffusion soient réversibles par rapport à la mesure de Gibbs sur une variété riemannienne. Cette condition dépend continûment d'un paramètre λ∈0,1 qui interpole entre les conventions d'Itô (λ=0), Stratonovich (λ=1/2) et Klimontovich (λ=1). Pour les systèmes lent-rapide réversibles avec structure de diffusion bloc-diagonale, les auteurs démontrent à l'aide de la théorie des formes de Dirichlet que la réversibilité et l'interprétation du bruit de Klimontovich sont préservées lors de la renormalisation. En particulier, la dynamique effective des variables lentes obtenue par projection sur une variété de dimension inférieure préserve l'interprétation de Klimontovich et reste réversible par rapport à la mesure de Gibbs marginale/énergie libre.

Contexte et motivation de la recherche

Problème fondamental

Le problème fondamental que cet article résout est : Comment choisir l'interprétation de l'intégrale stochastique (Itô, Stratonovich ou Klimontovich) dans le cas du bruit multiplicatif (bruit dépendant de l'état) pour assurer la réversibilité du processus de diffusion par rapport à une mesure de Gibbs donnée ?

Importance du problème

  1. Fondements de la mécanique statistique: La réversibilité est un concept fondamental de la mécanique statistique d'équilibre. Elle ne requiert pas seulement l'existence d'une mesure invariante, mais aussi la satisfaction de la condition de bilan détaillé, correspondant à un taux de production d'entropie nul.
  2. Algorithmes d'échantillonnage: Dans les méthodes MCMC et la descente de gradient stochastique, la dynamique de Langevin réversible est un outil fondamental.
  3. Modélisation physique: Le bruit multiplicatif apparaît naturellement dans l'hydrodynamique des fluctuations, les modèles de matière active et autres systèmes physiques.
  4. Optimisation de la vitesse de convergence: L'optimisation de la matrice de diffusion σ peut accélérer la convergence vers l'équilibre, ce qui est crucial pour les problèmes d'échantillonnage.

Limitations des approches existantes

  1. Restriction au bruit additif: L'équation de Langevin sursamortie classique dX_t = -∇V(X_t)dt + √2dW_t s'applique uniquement aux matrices de volatilité constantes.
  2. Dilemme de l'interprétation d'Itô: Sous l'interprétation d'Itô, le bruit multiplicatif détruit presque toujours la réversibilité, sauf sous des conditions extrêmement restrictives de coordonnées harmoniques.
  3. Absence de cadre unifié: La relation entre les différentes interprétations du bruit et leur impact sur la réversibilité manque d'étude systématique.
  4. Préservation lors de la renormalisation inconnue: Dans les systèmes multi-échelles, il est encore incertain si l'interprétation du bruit et la réversibilité sont préservées après renormalisation.

Motivation de la recherche

Le point de départ des auteurs est d'établir un cadre géométrique unifié qui caractérise systématiquement les conditions algébriques de réversibilité sous différentes interprétations du bruit, et de démontrer que l'interprétation de Klimontovich possède des propriétés de préservation particulières lors de la renormalisation. Cela fournit une base théorique pour la modélisation de dynamiques réversibles avec des géométries et des structures de bruit non triviales.

Contributions principales

  1. Conditions de réversibilité unifiées: Déduction de la condition algébrique paramétrée (2λ-1)∇^c·(σσ^T) = 2λσ(∇^c·σ^T), qui caractérise complètement la réversibilité de l'EDS par rapport à la mesure de Gibbs, où λ∈0,1 interpole les différentes interprétations du bruit.
  2. Caractérisation précise des trois bruits:
    • Bruit d'Itô (λ=0): Requiert que la divergence covariante s'annule ∇^c·(σσ^T)=0, équivalent à la condition de coordonnées harmoniques, extrêmement restrictif.
    • Bruit de Stratonovich (λ=1/2): Requiert σ∇^c·σ^T=0, une condition d'incompressibilité du bruit.
    • Bruit de Klimontovich (λ=1): Toujours satisfait en dimension un, automatiquement satisfait en dimension supérieure pour σ diagonal ou bloc-diagonal.
  3. Théorème de préservation lors de la renormalisation: Démonstration que pour les systèmes lent-rapide avec sources de bruit indépendantes, l'interprétation du bruit de Klimontovich et la réversibilité sont simultanément préservées dans la limite de moyenne, et la dynamique effective reste réversible par rapport à la mesure de Gibbs marginale.
  4. Méthode des formes de Dirichlet: Utilisation de la théorie de convergence de Mosco de Kuwae-Shioya, fournissant un cadre variationnel rigoureux pour la renormalisation des systèmes multi-échelles, alternative élégante aux méthodes classiques de générateur.
  5. Perspective géométrique: Interprétation de la matrice de diffusion M=σσ^T comme métrique riemannienne g=M^{-1}, unifiant le traitement des corrections de bruit et des corrections géométriques dans ce cadre géométrique.

Détails de la méthode

Définition de la tâche

Entrées:

  • Fonction potentiel lisse et inversible V: ℝ^d → ℝ (satisfaisant lim_{|x|→∞}V(x)=+∞)
  • Champ de matrice lisse et inversible σ: ℝ^d → GL(d) (matrice de volatilité)
  • Paramètre de bruit λ∈0,1

Sorties: Déterminer si l'EDS dXt=−σ(Xt)σT(Xt)∇V(Xt)dt+2σ(Xt)◦λdWtdX_t = -σ(X_t)σ^T(X_t)∇V(X_t)dt + \sqrt{2}σ(X_t) ◦_λ dW_t est réversible par rapport à la mesure de Gibbs G(dx)=1ZVe−V(x)volM(dx)G(dx) = \frac{1}{Z_V}e^{-V(x)}\text{vol}_M(dx) où vol_M(dx)=√(ω_M(x))dx est l'élément de volume riemannien, ω_M(x)=det(M^{-1}(x)).

Contraintes: σ doit satisfaire la condition d'ellipticité uniforme, c'est-à-dire qu'il existe C>0 tel que ξ^T M(x)ξ ≥ C‖ξ‖² pour tous x et ξ.

Cadre théorique

1. Configuration géométrique

Interprétation de l'inverse de la matrice de diffusion M(x)=σ(x)σ^T(x) induite par g(x)=M^{-1}(x) comme métrique riemannienne, rendant (ℝ^d,g) une variété riemannienne. Sur cette variété :

  • Connexion de Levi-Civita: Connexion unique sans torsion et compatible avec la métrique, satisfaisant ∇^c_ℓ g_=0
  • Symboles de Christoffel: Γijk=12gkℓ(∂igjℓ+∂jgiℓ−∂ℓgij)Γ^k_{ij} = \frac{1}{2}g^{kℓ}(∂_i g_{jℓ} + ∂_j g_{iℓ} - ∂_ℓ g_{ij})
  • Divergence covariante: Pour un champ de matrice A, la divergence covariante des lignes est définie comme (∇c⋅A)j=∂iAji+ΓikiAjk(∇^c·A)^j = ∂_i A^{ji} + Γ^i_{ik}A^{jk}

2. Analyse du générateur

Le générateur infinitésimal de l'EDS (8) est : LMf=ΔMf+B⋅∇f\mathcal{L}_M f = Δ_M f + B·∇f où :

  • Opérateur de Laplace-Beltrami: ΔMf=1ωM∂i(ωMMij∂jf)Δ_M f = \frac{1}{\sqrt{ω_M}}∂_i(\sqrt{ω_M}M^{ij}∂_j f)
  • Dérive corrigée: Bi=−Mij(∇V)j+2λ(∂jMji−σiℓ∂kσkℓ)+MkjΓkjiB^i = -M^{ij}(∇V)_j + 2λ(∂_j M^{ji} - σ^{iℓ}∂_k σ^{kℓ}) + M^{kj}Γ^i_{kj}

Le deuxième terme est la correction du bruit (dépendant de λ), le troisième terme est la correction géométrique.

3. Générateur réversible

Le générateur réversible de la mesure de Gibbs (référence BGL13) est : LGf=ΔMf+Γ(−V,f)\mathcal{L}_G f = Δ_M f + Γ(-V,f) où Γ(-V,f)=-M^{ij}(∂_j V)∂_i f est l'opérateur carré-du-champ. Utilisant la propriété de contraction des symboles de Christoffel : ΓikjMik=(∇c⋅(σσT))jΓ^j_{ik}M^{ik} = (∇^c·(σσ^T))^j on peut écrire le générateur réversible comme : LGf=Mij∂ij2f+(∇c⋅(σσT))j∂jf−Mij(∂iV)∂jf\mathcal{L}_G f = M^{ij}∂²_{ij}f + (∇^c·(σσ^T))^j∂_j f - M^{ij}(∂_i V)∂_j f

Théorème principal (Théorème 3.1)

Énoncé: G(dx) est réversible par rapport à X_t (c'est-à-dire que le générateur est auto-adjoint dans L²(G)) si et seulement si : (2λ−1)∇c⋅(σσT)=2λσ(∇c⋅σT)(2λ-1)∇^c·(σσ^T) = 2λσ(∇^c·σ^T)

Esquisse de la preuve:

  1. Première étape: Écriture de la représentation en coordonnées du générateur réversible L_G (comme ci-dessus)
  2. Deuxième étape: Utilisation de la formule de transformation du bruit (56), écriture de L sous la forme incluant la correction λ : Lf=Mij∂ij2f+2λ[(∇c⋅(σσT))−σ(∇c⋅σT)]j∂jf−Mij(∂iV)∂jf\mathcal{L}f = M^{ij}∂²_{ij}f + 2λ[(∇^c·(σσ^T)) - σ(∇^c·σ^T)]^j∂_j f - M^{ij}(∂_i V)∂_j f
  3. Troisième étape: Utilisation de la compatibilité métrique (17) pour prouver que le terme de correction géométrique égale la divergence covariante : −ΓikjMik=(∇c⋅(σσT))j-Γ^j_{ik}M^{ik} = (∇^c·(σσ^T))^j
  4. Quatrième étape: Comparaison des coefficients de L=L_G, obtenant que la réversibilité est équivalente à : 2λ[(∇c⋅(σσT))−σ(∇c⋅σT)]=(∇c⋅(σσT))2λ[(∇^c·(σσ^T)) - σ(∇^c·σ^T)] = (∇^c·(σσ^T)) La réorganisation donne la condition du théorème.

Lemme clé (Proposition 2.3): Preuve que [∇c⋅(σσT)]j−[σ(∇c⋅σT)]j=[∇⋅(σσT)]j−[σ(∇⋅σT)]j[∇^c·(σσ^T)]^j - [σ(∇^c·σ^T)]^j = [∇·(σσ^T)]^j - [σ(∇·σ^T)]^j c'est-à-dire que les termes des symboles de Christoffel s'annulent complètement dans cette combinaison, permettant d'exprimer la même condition en utilisant la divergence euclidienne.

Points d'innovation technique

  1. Introduction de la divergence covariante des lignes: La divergence covariante standard sur le tenseur (2,0) M=σσ^T s'annule identiquement en raison de la compatibilité métrique. La divergence covariante des lignes définie par l'auteur capture l'analogue géométrique correct.
  2. Cadre d'interpolation unifié: Unification des trois interprétations du bruit par un seul paramètre λ, rendant la condition de réversibilité continûment dépendante de λ, révélant les connexions intrinsèques entre elles.
  3. Pont entre analyse géométrique et stochastique: Connexion du choix de l'intégrale stochastique (problème d'analyse) avec les opérations covariantes en géométrie riemannienne (problème géométrique).
  4. Clarification de la forme covariante de Graham: Démonstration que sous la condition plus faible σ∇^c·σ^T=0, l'EDS de Stratonovich est équivalente à l'EDS d'Itô, améliorant le résultat de Diósi 2024.

Théorie de la renormalisation (Section 4)

Configuration du système lent-rapide

Considération d'un système bidimensionnel (généralisation à dimension supérieure directe) :

undefined