[论文解读] On the Convergence of Time Splitting Methods for Quantum Dynamics in the Semiclassical Regime

François Golse, Shi Jin+1 more|arXiv (Cornell University)|Jun 8, 2019

Numerical methods for differential equations参考文献 24被引用 21

一句话总结

本文为半经典 regime 下 von Neumann 方程的时间分裂方法——Lie-Trotter 与 Strang 方法——建立了关于 Planck 常数的统一误差估计。通过采用一种类似于 2 阶 Wasserstein 距离的伪度量，比较量子密度算符与经典相空间概率密度，证明了数值收敛性与 $\hbar$ 无关，从而在 $\hbar \to 0$ 时仍可使用 $\hat{\Delta t} = O(1)$ 的时间步长。关键结果是保持在 $\hbar$ 上一致的显式误差界，解决了半经典量子动力学模拟中长期存在的挑战。

ABSTRACT

By using the pseudo-metric introduced in [F. Golse, T. Paul: Archive for Rational Mech. Anal. 223 (2017) 57-94], which is an analogue of the Wasserstein distance of exponent $2$ between a quantum density operator and a classical (phase-space) density, we prove that the convergence of time splitting algorithms for the von Neumann equation of quantum dynamics is uniform in the Planck constant $\\hbar$. We obtain explicit uniform in $\\hbar$ error estimates for the first order Lie-Trotter, and the second order Strang splitting methods.

研究动机与目标

• 为解决量子动力学模拟中解析小 de Broglie 波长 ($\sim \hbar$) 的计算挑战。
• 严格证明时间分裂谱方法中使用 $O(1)$ 时间步长的合理性，此前该方法仅在形式层面得到支持。
• 建立在 Planck 常数 $\hbar$ 上一致的误差估计，使半经典 regime 下的可靠数值模拟成为可能。
• 通过使用能捕捉混合量子态半经典行为的相空间度量，弥合量子动力学与经典极限之间的鸿沟。

提出的方法

• 作者采用一种受 2 阶 Wasserstein 距离启发的伪度量，该度量定义在量子密度算符与经典相空间概率密度之间。
• 该度量用于比较通过时间分裂方法获得的量子态的 Wigner 变换与经典输运解之间的差异。
• 分析依赖于 Wigner 变换，将量子演化与相空间中的经典 Liouville 动力学联系起来。
• 通过结合迹类与 Hilbert-Schmidt 算子范数，并利用量子与经典测度之间的耦合技术，推导出误差估计。
• 该方法利用已知的指数分裂格式（Lie-Trotter 与 Strang）的算子范数界，并将其应用于混合态的 von Neumann 方程。
• 一个关键技术工具是使用相干态 $|q,p\rangle$ 作为构建量子与经典分布之间耦合的基。

实验结果

研究问题

• RQ1 时间分裂方法在 von Neumann 方程中能否在半经典 regime 下实现与 $\hbar$ 无关的时间步长的收敛？
• RQ2 是否存在一种合适的度量，使得在量子动力学模拟中可实现 $\hbar$-一致的误差估计？
• RQ3 当 $\hbar \to 0$ 时，Lie-Trotter 与 Strang 分裂格式的误差表现如何？
• RQ4 Wigner 变换框架能否严格证明时间分裂谱方法中使用 $O(1)$ 时间步长的合理性？
• RQ5 在半经典极限下，误差对 $\hbar$、$\Delta t$ 与初始数据的依赖关系如何？

主要发现

• 一阶 Lie-Trotter 分裂方法的误差界为 $\Delta t + \min(\frac{\Delta t}{\hbar}, \sqrt{\hbar})$，当 $\hbar = \Delta t^{4/3}$ 时简化为 $\Delta t + \Delta t^{1/3}$，且该界在 $\hbar$ 上一致。
• 二阶 Strang 分裂方法的误差界为 $\Delta t + \min(\frac{\Delta t^2}{\hbar}, \sqrt{\hbar})$，当 $\hbar = \Delta t^{4/3}$ 时简化为 $\Delta t + \Delta t^{2/3}$，且该界在 $\hbar$ 上一致。
• 误差估计在 $\hbar$ 上一致，意味着收敛速率不会随 $\hbar \to 0$ 而恶化，这对半经典模拟至关重要。
• 所使用的伪度量是 2 阶 Wasserstein 距离的量子-经典类比，使量子密度算符与经典相空间密度之间得以比较。
• 分析表明，若使用标准的 $L^2$ 或 Sobolev 范数，误差将呈现 $O((\Delta t / \hbar)^m)$ 的尺度，这会排除使用 $O(1)$ 时间步长，而所选度量避免了此问题。
• 结果验证了时间分裂谱方法在半经典 regime 下的经验成功，并为其使用 $\Delta t = O(1)$ 提供了严格的理论基础。