[论文解读] The strong convergence of operator-splitting methods for the Langevin dynamics model

Adam Telatovich, Xiantao Li|arXiv (Cornell University)|Jun 13, 2017

Theoretical and Computational Physics参考文献 34被引用 6

一句话总结

本文针对带有加性噪声的朗之万动力学模型中的算子分裂方法，提供了严格的强收敛性分析。通过利用Kunita的解表示形式与伊藤-泰勒展开，作者推导出新型高阶算子分裂格式，其强收敛阶数最高可达3阶，表明标准分裂格式（如SVV和朴素分裂）仅能达到1阶。

ABSTRACT

We study the strong convergence of some operator-splitting methods for the Langevin dynamics model with additive noise. It will be shown that a direct splitting of deterministic and random terms, including the symmetric splitting methods, only offers strong convergence of order 1. To improve the order of strong convergence, a new class of operator-splitting methods based on Kunita's solution representation are proposed. We present stochastic algorithms with strong orders up to 3. Both mathematical analysis and numerical evidence are provided to verify the desired order of accuracy.

研究动机与目标

• 分析广泛使用的算子分裂方法在带有加性噪声的朗之万动力学模型中的强收敛行为。
• 揭示现有分裂格式（如随机速度Verlet（SVV）和朴素分裂）的根本局限性，表明其仅能达到强收敛阶1。
• 基于伊藤SDE的Kunita解表示，提出一类新型高阶算子分裂格式。
• 通过数学分析与数值实验验证理论收敛阶。
• 研究所提出的高阶方法在过阻尼与欠阻尼极限下的行为。

提出的方法

• 作者采用Kunita对伊藤SDE的解表示形式，将解表达为包含漂移与扩散微分算子对易子的算子指数形式。
• 通过系统性地截断Kunita展开至更高阶次，构造高阶数值格式，其系数通过反复对漂移与扩散项应用微分算子得到。
• 该方法涉及计算至三阶的多重伊藤积分，包括类似 $ \bigint_{0}^{\tau} \bigint_{0}^{s_3} \bigint_{0}^{s_2} dW_{s_1} ds_2 ds_3 $ 的项，这些项被精确解析计算以保证精度。
• 通过多指标 $ \boldsymbol{\nu} $ 对项进行分组，格式阶数由总阶数 $ l(\boldsymbol{\nu}) + n(\boldsymbol{\nu}) $ 决定，其中 $ l $ 表示维纳积分的个数，$ n $ 表示微分算子应用的次数。
• 对于加性噪声，强收敛阶3与阶2.5的格式完全一致，因为所有额外的系数函数 $ g_{\boldsymbol{\nu}} $（当 $ \boldsymbol{\nu} \notin \tilde{\boldsymbol{\nu}}_{2.5} $ 时）由于常数或零导数而消失。
• 所得到的格式为显式、可积的子步骤，既保持了经典分裂方法的实现简便性，又实现了更高的强收敛阶。

实验结果

研究问题

• RQ1 为何像SVV和朴素分裂这样的标准算子分裂方法尽管在分子动力学模拟中广泛应用，却仅能达到强收敛阶1？
• RQ2 在带有加性噪声的朗之万动力学模型中，能否实现算子分裂格式的更高阶强收敛？
• RQ3 基于Kunita解表示与伊藤-泰勒展开构造高阶格式的理论基础是什么？
• RQ4 通过强收敛误差范数评估时，新格式在精度上与经典方法相比如何？
• RQ5 新提出的高阶格式是否在过阻尼与欠阻尼极限下保持正确的物理行为？

主要发现

• 尽管SVV和朴素分裂方法在子步骤中精确且广受欢迎，但其强收敛阶仅为1阶。
• 通过应用Kunita的解表示并截断伊藤-泰勒展开至更高阶次，作者构造出强收敛阶最高可达3阶的新算子分裂格式。
• 对于加性噪声，强收敛阶3与阶2.5的格式完全相同，因为所有额外的系数函数 $ g_{\boldsymbol{\nu}} $（当 $ \boldsymbol{\nu} \notin \tilde{\boldsymbol{\nu}}_{2.5} $ 时）由于 $ g_{(1,0,0)} $ 与 $ g_{(1,0,0,0)} $ 在 $ L^1 $ 与 $ L^0 $ 运算下为常数或零而消失。
• 新格式在数学上被证明可达到预期的强收敛阶，数值结果也证实了理论预测。
• 该方法在过阻尼与欠阻尼两种极限下均适用，且在两种情形下均表现出一致的收敛行为。
• 该方法避免了对力的高阶导数的需求（与伊藤-泰勒方法不同），也避免了多次力评估（与龙格-库塔方法不同），因此在分子动力学模拟中具有实际可行性。