[论文解读] The projection of a nonlocal mechanical system onto the irreversible generalized Langevin equation, II: Numerical simulations

Marc Vogt, Rigoberto Hernandez|arXiv (Cornell University)|Jan 14, 2005

Spectroscopy and Quantum Chemical Studies被引用 10

一句话总结

该论文表明，不可逆广义朗之万方程（iGLE）——用于在随机动力学中建模非平稳摩擦的方程——可作为确定性开放哈密顿系统的投影而物理实现。数值模拟证实，哈密顿系统中的时间依赖能量和溶剂力相关性与iGLE精确匹配，从而验证了iGLE作为物理基础模型而非现象学构造的有效性，其中非局部耗散项在此一致性的形成中起着关键作用。

ABSTRACT

The irreversible generalized Langevin equation (iGLE) contains a nonstationary friction kernel that in certain limits reduces to the GLE with space-dependent friction. For more general forms of the friction kernel, the iGLE was previously shown to be the projection of a mechanical system with a time-dependent Hamiltonian. [R. Hernandez, J. Chem. Phys. 110, 7701 (1999)] In the present work, the corresponding open Hamiltonian system is further explored. Numerical simulations of this mechanical system illustrate that the time dependence of the observed total energy and the correlations of the solvent force are in precise agreement with the projected iGLE.

研究动机与目标

• 通过将iGLE从具有时间依赖哈密顿量的确定性开放哈密顿系统中投影，建立其作为物理可实现模型。
• 验证iGLE的非平稳摩擦核和涨落-耗散关系（FDR）可自然地从具有时间依赖哈密顿量的机械系统中产生。
• 证明iGLE中非局部耗散项在重现底层哈密顿系统动力学方面具有必要性。
• 通过全机械系统的数值积分验证iGLE的随机行为，显示其速度自相关函数和力相关性与模拟结果一致。
• 将非平稳随机动力学的物理基础从现象学模型扩展至与确定性哈密顿框架等价的体系。

提出的方法

• iGLE被推导为具有时间依赖哈密顿量的机械系统的投影，其中摩擦核由时间依赖函数g(t)调制。
• 系统由一个标记粒子与200个谐振子组成的浴组成，总哈密顿量被显式构建以重现iGLE的动力学。
• iGLE中的非局部耗散项源自浴的时间演化耦合，其显式计算对准确模拟至关重要。
• 使用辛积分器对全哈密顿方程进行数值积分，以在长时间内保持能量和动力学的守恒。
• 在每个时间步长通过关系式ξ(t) = v̇ + γ₀g(t)z₂(t) − F(t)计算随机力ξ(t)，以确保符合FDR。
• 直接比较哈密顿系统与投影后的iGLE之间的速度自相关函数和力相关函数。

实验结果

研究问题

• RQ1 不可逆广义朗之万方程（iGLE）能否作为具有时间依赖哈密顿量的确定性机械系统的物理投影而实现？
• RQ2 iGLE中的非局部耗散项是否在重现底层哈密顿系统正确动力学方面具有必要性？
• RQ3 哈密顿系统中的速度自相关函数和溶剂力相关性是否在定量上与iGLE的预测完全一致？
• RQ4 即使在非平稳条件下，哈密顿系统中随机力的涨落-耗散关系（FDR）是否仍然成立？
• RQ5 摩擦的初始强度和g(t)的时间演化如何影响哈密顿系统与iGLE之间的一致性？

主要发现

• 数值模拟显示，在所有情况和时间尺度下，哈密顿系统中标记粒子的速度自相关函数与iGLE的预测完全一致。
• 哈密顿系统中溶剂力的自相关函数在所有时间点和所有参数选择下均满足涨落-耗散关系（FDR）。
• 非局部耗散项至关重要：若省略该部分，模拟结果将导致哈密顿系统与iGLE之间的一致性被破坏。
• 在自由粒子情况下，早期时间（t = -4.0）的动力学为弹道运动，表现为线性速度自相关函数，与iGLE预测一致。
• 在谐振子势阱情况下，由于束缚效应，系统对摩擦变化的敏感性降低，但当包含非局部项时，iGLE仍能准确重现动力学。
• 所有情况在长时间（t = 4.0）均达到相同平衡态，证实iGLE正确捕捉了系统的长时统计行为。