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

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

基本信息

  • 论文ID: 2511.03347
  • 标题: Reversibility, covariance and coarse-graining for Langevin dynamics: On the choice of multiplicative noise
  • 作者: Mario Ayala (TU München), Nicolas Dirr (Cardiff University), Grigorios A. Pavliotis (Imperial College London), Johannes Zimmer (TU München)
  • 分类: math.PR (概率论)
  • 发表时间: 2025年11月6日 (arXiv预印本)
  • 论文链接: https://arxiv.org/abs/2511.03347

摘要

本文研究随机微分方程(SDEs)中可逆性、几何结构与乘性噪声选择(特别是Itô、Stratonovich、Klimontovich解释)之间的相互作用。基于统一的几何框架,作者推导了扩散过程关于黎曼流形上Gibbs测度可逆的代数条件。该条件连续依赖于参数λ∈0,1,它在Itô(λ=0)、Stratonovich(λ=1/2)和Klimontovich(λ=1)约定之间插值。对于具有块对角扩散结构的可逆慢-快系统,作者利用Dirichlet形式理论证明,可逆性和Klimontovich噪声解释在粗粒化下都被保持。特别地,通过投影到低维流形获得的慢变量的有效动力学保留了Klimontovich解释,并关于边缘Gibbs测度/自由能保持可逆。

研究背景与动机

核心问题

本文要解决的核心问题是:在乘性噪声(状态依赖噪声)情形下,如何选择随机积分的解释(Itô、Stratonovich或Klimontovich)以保证扩散过程关于给定Gibbs测度的可逆性?

问题重要性

  1. 统计力学基础: 可逆性是平衡态统计力学的核心概念,它不仅要求不变测度存在,还要求满足细致平衡条件,对应零熵产生率
  2. 采样算法: 在MCMC方法和随机梯度下降等算法中,可逆Langevin动力学是基础工具
  3. 物理建模: 乘性噪声自然出现在涨落流体力学、活性物质模型等物理系统中
  4. 收敛速度优化: 通过优化扩散矩阵σ可以加速平衡态收敛,这在采样问题中至关重要

现有方法的局限性

  1. 加性噪声的限制: 经典的过阻尼Langevin方程 dXt = -∇V(Xt)dt + √2dWt 仅适用于常数波动矩阵
  2. Itô解释的困境: 在Itô解释下,乘性噪声几乎总是破坏可逆性,除非满足极强的调和坐标条件
  3. 缺乏统一框架: 不同噪声解释之间的关系及其对可逆性的影响缺乏系统研究
  4. 粗粒化保持性未知: 在多尺度系统中,噪声解释和可逆性在粗粒化后是否保持尚不清楚

研究动机

作者的出发点是建立一个统一的几何框架,系统地刻画不同噪声解释下可逆性的代数条件,并证明Klimontovich解释在粗粒化过程中具有特殊的保持性质,这为建模具有非平凡几何和噪声结构的可逆动力学提供了理论基础。

核心贡献

  1. 统一的可逆性条件: 推导了参数化的代数条件 (2λ-1)∇^c·(σσ^T) = 2λσ(∇^c·σ^T),它完全刻画了SDE关于Gibbs测度的可逆性,其中λ∈0,1插值不同噪声解释
  2. 三种噪声的精确刻画:
    • Itô噪声(λ=0): 要求协变散度消失 ∇^c·(σσ^T)=0,等价于调和坐标条件,极其严格
    • Stratonovich噪声(λ=1/2): 要求 σ∇^c·σ^T=0,是一种噪声不可压缩性条件
    • Klimontovich噪声(λ=1): 在一维情形总是满足,高维情形对对角或块对角σ自动满足
  3. 粗粒化保持定理: 证明了对于具有独立噪声源的慢-快系统,Klimontovich噪声解释和可逆性在平均化极限下同时保持,有效动力学仍关于边缘Gibbs测度可逆
  4. Dirichlet形式方法: 利用Kuwae-Shioya的Mosco收敛理论,为多尺度系统的粗粒化提供了严格的变分框架,这是对经典生成元方法的优雅替代
  5. 几何视角: 将扩散矩阵M=σσ^T解释为黎曼度量g=M^{-1},在此几何框架下统一处理噪声修正和几何修正

方法详解

任务定义

输入:

  • 光滑约束势函数 V: ℝ^d → ℝ (满足 lim_{|x|→∞}V(x)=+∞)
  • 光滑可逆矩阵场 σ: ℝ^d → GL(d) (波动率矩阵)
  • 噪声参数 λ∈0,1

输出: 判断SDE 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 是否关于Gibbs测度 G(dx)=1ZVe−V(x)volM(dx)G(dx) = \frac{1}{Z_V}e^{-V(x)}\text{vol}_M(dx) 可逆,其中 vol_M(dx)=√(ω_M(x))dx 是黎曼体积元,ω_M(x)=det(M^{-1}(x))。

约束: σ需满足一致椭圆性条件,即存在C>0使得 ξ^T M(x)ξ ≥ C‖ξ‖² 对所有x和ξ成立。

理论框架

1. 几何设置

将扩散矩阵 M(x)=σ(x)σ^T(x) 诱导的逆 g(x)=M^{-1}(x) 视为黎曼度量,使(ℝ^d,g)成为黎曼流形。在此流形上:

  • Levi-Civita联络: 唯一的无挠且与度量相容的联络,满足 ∇^c_ℓ g_=0
  • Christoffel符号: Γijk=12gkℓ(∂igjℓ+∂jgiℓ−∂ℓgij)Γ^k_{ij} = \frac{1}{2}g^{kℓ}(∂_i g_{jℓ} + ∂_j g_{iℓ} - ∂_ℓ g_{ij})
  • 协变散度: 对矩阵场A,定义行协变散度为 (∇c⋅A)j=∂iAji+ΓikiAjk(∇^c·A)^j = ∂_i A^{ji} + Γ^i_{ik}A^{jk}

2. 生成元分析

SDE (8)的无穷小生成元为: LMf=ΔMf+B⋅∇f\mathcal{L}_M f = Δ_M f + B·∇f 其中:

  • Laplace-Beltrami算子: ΔMf=1ωM∂i(ωMMij∂jf)Δ_M f = \frac{1}{\sqrt{ω_M}}∂_i(\sqrt{ω_M}M^{ij}∂_j f)
  • 修正漂移: 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}

第二项是噪声修正(依赖于λ),第三项是几何修正。

3. 可逆生成元

Gibbs测度G(dx)的可逆生成元(参考BGL13)为: LGf=ΔMf+Γ(−V,f)\mathcal{L}_G f = Δ_M f + Γ(-V,f) 其中 Γ(-V,f)=-M^{ij}(∂_j V)∂_i f 是carré-du-champ算子。利用Christoffel符号的收缩性质: ΓikjMik=(∇c⋅(σσT))jΓ^j_{ik}M^{ik} = (∇^c·(σσ^T))^j 可将可逆生成元写为: 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

主定理(Theorem 3.1)

陈述: G(dx)关于X_t可逆(即生成元在L²(G)中自伴)当且仅当: (2λ−1)∇c⋅(σσT)=2λσ(∇c⋅σT)(2λ-1)∇^c·(σσ^T) = 2λσ(∇^c·σ^T)

证明思路:

  1. 第一步: 写出可逆生成元L_G的坐标表示(如上)
  2. 第二步: 利用噪声转换公式(56),将L写为包含λ修正的形式: 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. 第三步: 利用度量相容性(17)证明几何修正项等于协变散度: −ΓikjMik=(∇c⋅(σσT))j-Γ^j_{ik}M^{ik} = (∇^c·(σσ^T))^j
  4. 第四步: 比较L=L_G的系数,得到可逆性等价于: 2λ[(∇c⋅(σσT))−σ(∇c⋅σT)]=(∇c⋅(σσT))2λ[(∇^c·(σσ^T)) - σ(∇^c·σ^T)] = (∇^c·(σσ^T)) 整理即得定理条件。

关键引理(Proposition 2.3): 证明了 [∇c⋅(σσT)]j−[σ(∇c⋅σT)]j=[∇⋅(σσT)]j−[σ(∇⋅σT)]j[∇^c·(σσ^T)]^j - [σ(∇^c·σ^T)]^j = [∇·(σσ^T)]^j - [σ(∇·σ^T)]^j 即Christoffel符号项在该组合中完全抵消,这使得可以用欧氏散度表达相同的条件。

技术创新点

  1. 行协变散度的引入: 标准协变散度对(2,0)-张量M=σσ^T因度量相容性恒为零,作者定义的行协变散度捕捉了正确的几何类比
  2. 统一插值框架: 通过单一参数λ统一三种噪声解释,使可逆性条件连续依赖于λ,揭示了它们之间的内在联系
  3. 几何与随机分析的桥梁: 将随机积分的选择(分析问题)与黎曼几何中的协变运算(几何问题)联系起来
  4. Graham协变形式的澄清: 指出在较弱条件σ∇^c·σ^T=0下,Stratonovich SDE就等价于Itô SDE,改进了Diósi2024的结果

粗粒化理论(Section 4)

慢-快系统设置

考虑二维系统(推广到高维直接):

dX^n_t = -σ²_1(X^n_t,Y^n_t)∂_x V(X^n_t,Y^n_t)dt + \sqrt{2}σ_1(X^n_t,Y^n_t)◦_K dW¹_t \\ dY^n_t = -nσ²_2(X^n_t,Y^n_t)∂_y V(X^n_t,Y^n_t)dt + \sqrt{2n}σ_2(X^n_t,Y^n_t)◦_K dW²_t \end{cases}$$ 其中: - X是慢变量,Y是快变量 - 噪声源W¹和W²独立(对应块对角扩散矩阵) - 均采用Klimontovich解释(λ=1) - 联合Gibbs测度为 μ(dxdy) = Z^{-1}_V e^{-V(x,y)}dxdy ### 平均化极限 当n→∞时,经典结果[Kha68]表明X^n_t收敛到有效过程: $$d\bar{X}_t = \bar{b}(\bar{X}_t)dt + \sqrt{2}\bar{σ}_1(\bar{X}_t)dW_t$$ 其中有效系数为: $$\bar{b}(x) = \frac{1}{Z_V(x)}\int_ℝ [∂_x σ²_1(x,y) - σ²_1(x,y)∂_x V(x,y)]e^{-V(x,y)}dy$$ $$\bar{σ}²_1(x) = \frac{1}{Z_V(x)}\int_ℝ σ²_1(x,y)e^{-V(x,y)}dy$$ 边缘测度为 μ_∞(dx) ∝ ∫_ℝ e^{-V(x,y)}dy。 ### 保持性定理 **关键计算**: 通过微分关系 $$∂_x(\bar{σ}²_1(x)) = \bar{b}(x) + \bar{σ}²_1(x)\frac{∂_x Z_V(x)}{Z_V(x)}$$ 可将有效SDE改写为Klimontovich形式: $$d\bar{X}_t = -\bar{σ}²_1(\bar{X}_t)∂_x\ln\left(\int_ℝ e^{-V(\bar{X}_t,y)}dy\right)dt + \sqrt{2}\bar{σ}_1(\bar{X}_t)◦_K dW_t$$ 这直接证明了:**Klimontovich噪声解释在平均化下被保持**,且有效动力学关于边缘Gibbs测度μ_∞可逆。 ### Dirichlet形式方法 #### Mosco收敛框架 作者采用Kuwae-Shioya[2003]的理论,在变化Hilbert空间上研究Dirichlet形式的收敛: 1. **预极限空间**: $$H¹_n(μ) = \{f∈L²(ℝ²,μ) : \int_{ℝ²}∑_{ij}M^{ij}_n(∂_i f)(∂_j f)μ(dxdy)<∞\}$$ 其中 M_n = diag(σ²_1, nσ²_2) 2. **极限空间**: $$H¹(μ_∞) = \{f∈L²(ℝ,μ_∞) : \int_ℝ \bar{σ}²_1(x)(∂_x f)²μ_∞(dx)<∞\}$$ 3. **嵌入映射**: Φ_n: H¹(μ_∞)→H¹_n(μ) 定义为 (Φ_n f)(x,y)=f(x) **关键性质(Remark 4.1)**: Φ_n是等距嵌入,且保持Dirichlet形式: $$\mathcal{E}_n(Φ_n f) = \mathcal{E}(f), \quad ∀f∈H¹(μ_∞)$$ **主要结果(Theorem 4.1)**: Dirichlet形式序列{(E_n,D(E_n))}在Mosco意义下收敛到(E,D(E)),这保证了半群和预解式的收敛,从而严格证明了平均化原理。 ### 高维推广 对于(d+m)维系统,只要噪声源独立(块对角结构),所有结果直接推广。有效扩散参数为矩阵值Bochner积分: $$\bar{σ}_1(x)\bar{σ}^T_1(x) = \frac{1}{Z_V(x)}\int_{ℝ^m}σ_1(x,y)σ^T_1(x,y)e^{-V(x,y)}dy$$ 对于σ_1(x,y)=UΛ(x,y)U^T(U正交常数,Λ对角)的情形,作者通过显式计算(Section 4.7)验证了Klimontovich条件(55)在平均化后仍然成立。 ## 实验设置 **注**: 本文是纯理论数学论文,不包含数值实验。所有结果都是通过严格的数学证明得到的定理和命题。 ### 理论验证方法 1. **一维案例分析**(Section 1.1): 通过显式计算生成元,验证只有Klimontovich解释(λ=1)能保证任意σ(x)下的可逆性 2. **几何例子**: - 对角矩阵σ=diag(σ_1,...,σ_d)自动满足Klimontovich条件 - 块对角结构保持可逆性 - 常正交变换σ=UDU^T(U常数)保持条件 3. **反例构造**: 若U=U(x)依赖位置,则Klimontovich条件通常失效,说明噪声方向的旋转引入几何扭曲 ## 相关工作 ### 随机微分方程理论 1. **Itô积分** [Itô44]: 保证积分是鞅,是弱解和收敛性理论的基础 2. **Stratonovich积分** [Str66]: 保持经典微分链式法则,物理建模中自然 3. **Klimontovich积分** [Kli90, Hän78]: 与Fick扩散定律一致的Fokker-Planck方程[ER25] ### 可逆性与非平衡统计力学 1. **细致平衡条件** [BGL13, Pav14]: 自伴生成元等价于零熵产生 2. **熵产生理论** [MN03, CP23]: 不可逆成分的投影决定熵产生率 3. **涨落关系** [JQQ04]: Gallavotti-Cohen型涨落定理 ### 几何随机分析 1. **Graham协变形式** [Gra77, Gra85]: 首次提出Fokker-Planck方程的协变形式 2. **流形上的扩散** [Hsu02, Wei72]: 黎曼流形上的Laplace-Beltrami算子 3. **Diósi的协变条件** [Dió24]: 本文推广并改进了其Stratonovich-Itô等价性结果 ### 采样与优化 1. **预条件Langevin动力学** [Che+24]: 与Fisher-Rao梯度流的联系 2. **谱隙优化** [Lel+25]: 优化σ以最大化收敛速度 3. **随机梯度下降** [LTE19, DZ20]: Fokker-Planck视角 ### 多尺度与粗粒化 1. **平均化原理** [Kha68, PS08]: 经典的慢-快系统理论 2. **Dirichlet形式方法** [BR14]: 通过Dirichlet形式研究平均化 3. **Mosco收敛** [KS03, Kol06]: 变化Hilbert空间上的谱收敛理论 4. **两尺度收敛** [DDP23]: 周期势中的布朗运动 ### 数据驱动方法 1. **学习扩散矩阵** [Zhu+25]: 从轨迹数据学习可逆/不可逆成分 2. **发现演化方程** [Li+19, Hua+25]: 从粒子数据学习宏观算子 ## 结论与讨论 ### 主要结论 1. **噪声选择的决定性作用**: 在乘性噪声情形下,随机积分的解释(Itô、Stratonovich、Klimontovich)根本性地影响可逆性: - Itô解释过于严格(需调和坐标) - Klimontovich解释最为宽松(一维总成立,高维对常见结构成立) - Stratonovich介于两者之间 2. **Klimontovich噪声的优越性**: - 在粗粒化过程中自然保持 - 有效动力学自动保持可逆结构 - 边缘Gibbs测度保持Gibbs形式 3. **几何视角的统一性**: 将扩散矩阵视为黎曼度量,协变散度条件提供了内在的几何刻画,独立于坐标选择 4. **Dirichlet形式的优雅性**: Mosco收敛框架为多尺度分析提供了变分视角,避免了传统生成元方法的技术复杂性 ### 局限性 1. **块对角结构的要求**: 粗粒化保持性定理要求慢-快变量的噪声源独立(块对角M_n),这排除了耦合噪声的情形 2. **光滑性假设**: 要求σ和V是C²光滑的,对于非光滑或退化情形(如边界条件)理论不适用 3. **一致椭圆性**: 需要扩散矩阵一致正定,排除了部分自由度"冻结"的情形 4. **Gibbs测度的特殊性**: 理论专门针对Gibbs形式的测度,对于更一般的不变测度(如非平衡稳态)需要修改 5. **高维Klimontovich条件**: 在高维情形,Klimontovich条件(35)仍是非平凡的代数约束,不是对所有σ都成立 ### 未来方向 1. **耦合噪声的粗粒化**: 研究非块对角扩散矩阵情形下的平均化,可能需要修正的Klimontovich解释 2. **非平衡推广**: 将理论扩展到包含外力或循环流的非平衡稳态 3. **数值算法**: 基于Klimontovich离散化开发新的数值格式,保持离散层面的可逆性 4. **学习框架**: 从数据中同时学习σ、V和最优λ,实现数据驱动的可逆建模 5. **无穷维推广**: 将理论扩展到SPDE和无穷维流形(如概率测度空间) 6. **量子类比**: 探索与量子力学中路径积分表示的联系 ## 深度评价 ### 优点 1. **理论深度与原创性**: - 首次给出统一的参数化可逆性条件,揭示三种噪声解释的连续谱 - Klimontovich噪声在粗粒化下的保持性是全新发现,具有深刻的物理意义 - 行协变散度的引入巧妙解决了标准协变散度恒为零的困境 2. **数学严格性**: - 所有定理都有完整的证明,引用文献准确 - Mosco收敛的应用展示了高超的泛函分析技巧 - 度量相容性等几何性质的利用非常优雅 3. **框架的统一性**: - 将随机分析、黎曼几何、非平衡统计力学和多尺度分析有机结合 - 从一维到高维,从单一尺度到多尺度,理论框架一以贯之 - Dirichlet形式方法提供了与经典方法互补的视角 4. **物理直觉**: - 将σ解释为"噪声景观",小σ区域抑制涨落,大σ区域噪声主导 - 协变散度条件的几何意义清晰(调和坐标、不可压缩性等) - 粗粒化保持性符合热力学第二定律的精神 5. **写作质量**: - 结构清晰,从动机到理论到应用层层递进 - 一维例子(Section 1.1)提供了直观的入门 - 附录完整,包含噪声转换公式和Mosco收敛背景 ### 不足 1. **应用范围的限制**: - 缺乏数值实验验证理论预测(虽然是数学论文,但简单算例会增强说服力) - 块对角假设限制了实际系统的适用性(许多物理系统有耦合噪声) - 未讨论边界条件和约束系统 2. **技术细节的复杂性**: - 对非专家而言,协变导数和Christoffel符号的计算可能构成障碍 - Mosco收敛的抽象性较高,与物理直觉的联系不够直接 - 某些证明步骤(如Proposition 2.3)的动机不够清晰 3. **与物理文献的对话不足**: - Graham的工作[Gra77,Gra85]虽被引用,但对比分析不够深入 - 与涨落流体力学[MRZ25]的具体联系仅一笔带过 - 缺少与实验可观测量(如熵产生率测量)的联系 4. **开放问题的讨论**: - 对于不满足Klimontovich条件的σ,是否存在其他保持粗粒化的噪声解释? - 非块对角情形的粗粒化是否有部分保持性? - 离散时间近似的可逆性如何? ### 影响力评估 1. **对数学领域的贡献**: - 为随机分析与黎曼几何的交叉提供了新工具(行协变散度) - 推广了Dirichlet形式在多尺度问题中的应用 - 可能激发关于一般张量场散度理论的研究 2. **对统计力学的贡献**: - 澄清了不同噪声解释的热力学意义 - 为非平衡态理论提供了新的几何视角 - Klimontovich噪声的自然性得到数学支持 3. **对采样算法的潜在影响**: - 可指导设计保持可逆性的预条件Langevin算法 - 多尺度采样可利用粗粒化保持性 - 为自适应MCMC提供理论基础 4. **实用价值**: - 对于需要设计可逆动力学的应用(如分子动力学、材料科学),提供了明确的设计准则 - 粗粒化理论可应用于复杂系统的降维建模 - 数据驱动方法[Zhu+25]可利用本文框架 5. **可复现性**: - 纯理论结果,可复现性极高 - 所有定理都有证明,读者可独立验证 - 引用文献充分,便于追溯 ### 适用场景 1. **分子动力学模拟**: 设计保持Gibbs分布的Langevin热浴,特别是非均匀温度场情形 2. **贝叶斯采样**: 在复杂后验分布上设计高效MCMC算法,通过优化σ加速收敛 3. **多尺度建模**: 从微观随机模型推导宏观确定性或随机方程,保持热力学一致性 4. **活性物质**: 建模自驱动粒子系统的涨落流体力学,保持能量-动量守恒 5. **机器学习**: 随机优化算法的理论分析,理解不同离散化方案对收敛的影响 6. **金融数学**: 设计保持无套利条件的随机波动率模型 ## 参考文献(关键文献) 1. **[BGL13]** Bakry, Gentil, Ledoux. *Analysis and geometry of Markov diffusion operators*. Springer, 2013. (可逆扩散的经典参考) 2. **[Pav14]** Pavliotis. *Stochastic processes and applications*. Springer, 2014. (随机过程教材,本文多处引用) 3. **[KS03]** Kuwae & Shioya. "Convergence of spectral structures". *Comm. Anal. Geom.* 2003. (Mosco收敛理论基础) 4. **[Gra77]** Graham. "Covariant formulation of non-equilibrium statistical thermodynamics". *Z. Phys. B* 1977. (首次提出协变Fokker-Planck方程) 5. **[Kha68]** Khasminskij. "On the principle of averaging". *Kybernetika* 1968. (平均化原理经典文献) 6. **[CP23]** Costa & Pavliotis. "The entropy production of stationary diffusions". *J. Phys. A* 2023. (熵产生的现代理论) 7. **[Lel+25]** Lelièvre et al. "Optimizing the diffusion coefficient". *Math. Comp.* 2025. (扩散矩阵优化) 8. **[Dió24]** Diósi. "The covariant Langevin equation". *Rep. Math. Phys.* 2024. (本文推广的前期工作) --- **总体评价**: 这是一篇高质量的理论数学论文,在随机分析与黎曼几何的交叉领域做出了实质性贡献。主要创新在于统一刻画了不同噪声解释下的可逆性条件,并首次证明了Klimontovich噪声在粗粒化下的保持性。理论严格,框架优雅,对统计力学和采样算法有潜在的深远影响。主要局限在于应用范围(块对角假设)和缺乏数值验证。建议未来工作包括数值实现、耦合噪声推广和与实验的联系。论文适合发表在概率论或数学物理的顶级期刊。