[论文解读] Stochastic Simulation of Nonadiabatic Dynamics at Long Time

Daniel A. Uken, Alessandro Sergi+1 more|arXiv (Cornell University)|Mar 16, 2010

Quantum chaos and dynamical systems参考文献 1被引用 3

一句话总结

本文提出了一种广义随机采样方案，用于在Wigner-Heisenberg表示下进行非绝热动力学模拟，显著降低了长时间演化下的统计误差增长。通过使用定制的转移概率权重，过滤掉导致能量剧烈波动的轨迹，该方法将可访问的积分时间相比原始采样方案延长了整整一个数量级，如自旋-玻色子模型模拟所示，长时间下的误差条显著减小。

ABSTRACT

Using a generalized energy-conserving transition probability, it is shown how nonadiabatic calculations, within the Wigner-Heisenberg representation of quantum mechanics, can be reliably extended to far longer times than those allowed by a primitive sampling scheme. Tackling the spin-boson model as a paradigmatic example, substantial numerical evidence is provided that effective integration of the dynamics can be achieved for a wide range of temperatures and friction.

研究动机与目标

• 解决在长时间时间区间内非绝热量子动力学随机模拟中长期存在的统计误差增长问题。
• 通过改进量子跃迁的随机采样，提高非绝热动力学模拟的数值稳定性。
• 在原始采样方案的极限之外，实现对强系统- bath 耦合体系中非绝热过程的可靠模拟。
• 提供一种过滤不物理解轨迹的方法，这些轨迹在随机传播过程中会引发虚假的能量波动。
• 在自旋-玻色子模型中，展示广义采样方法在一系列阻尼和耦合参数下的有效性。

提出的方法

• 该方法在非绝热跃迁的随机采样中采用广义转移概率权重，其来源于线性化的Wigner-Heisenberg传播器。
• 该权重函数被设计用于抑制导致大能量波动的轨迹，从而有效过滤掉不物理解的随机跃迁。
• 该方法基于分段绝热轨迹表示，其中子系统能级之间的跃迁通过修改后的概率分布进行随机采样。
• 该采样方案在动量跳跃（MJ）近似框架内实现，确保了混合量子-经典动力学中能量守恒和反作用效应的一致性。
• 广义采样被应用于自旋-玻色子模型，时间演化通过轨迹的系综平均计算得出。
• 该方法使用阈值参数 $ c_{ m E} $ 控制高能波动事件的过滤，确保数值稳定性。

实验结果

研究问题

• RQ1 广义采样方案是否能降低长时间非绝热动力学模拟中的统计误差增长？
• RQ2 对高能波动轨迹的过滤如何影响随机非绝热动力学的长时间稳定性？
• RQ3 与原始采样相比，所提出的采样方法在多大程度上能够延长积分时间？
• RQ4 广义采样是否在非绝热动力学中保持能量守恒和物理一致性？
• RQ5 该方法是否能在自旋-玻色子模型的广泛阻尼和耦合参数范围内有效应用？

主要发现

• 广义采样方案显著降低了统计误差增长，在 $ t = 20 $ 之前误差条保持较小，而原始采样在 $ t = 10 $ 时已出现不可接受的大误差。
• 对于 $ eta = 3 $，$ eta = 0.25 $，和 $ eta = 1.0 $，广义采样在 $ t > 16 $ 时仍保持低误差水平，而原始采样在 $ t = 10 $ 后即失效。
• 在广义采样下，零跃迁轨迹（$ n=0 $）的比例随时间增加，表明不物理解跃迁受到抑制。
• 具有一个或两个跃迁的轨迹数量保持稳定且符合物理实际，而原始采样在跃迁统计上表现出不受控的波动。
• 该方法使积分时间相比以往方案延长了一个数量级，具体取决于数值参数。
• 结果为广义采样在包括高阻尼和强耦合区域在内的广泛参数范围内显著降低统计误差提供了强有力的数值证据。