[论文解读] Family of Gaussian wavepacket dynamics methods from the perspective of a nonlinear Schrödinger equation

Jiří Vaníček|arXiv (Cornell University)|Feb 20, 2023

Advanced Fiber Laser Technologies被引用 4

一句话总结

本文将多种高斯波包动力学方法——如Heller的非冻结高斯、变分高斯及单Hessian近似——统一于一个基于状态依赖二次势的非线性薛定谔方程框架下。推导出一般运动方程，证明了规范性、能量和辛结构的守恒，并提出一种新的‘单四次’变分方法，在保持低计算成本的同时，相较于局部立方近似显著提升了精度，且保留了关键几何性质。

ABSTRACT

Many approximate solutions of the time-dependent Schrödinger equation can be formulated as exact solutions of a nonlinear Schrödinger equation with an effective Hamiltonian operator depending on the state of the system. We show that several well-known Gaussian wavepacket dynamics methods, such as Heller's original thawed Gaussian approximation or Coalson and Karplus's variational Gaussian approximation, fit into this framework if the effective potential is a quadratic polynomial with state-dependent coefficients. We study such a nonlinear Schrödinger equation in general: in particular, we derive general equations of motion for the Gaussian's parameters, demonstrate the time reversibility and norm conservation, and analyze conservation of the energy, effective energy, and symplectic structure. We also describe efficient geometric integrators of arbitrary even orders of accuracy in the time step for the numerical solution of this nonlinear Schrödinger equation. The general presentation is illustrated by examples of this family of Gaussian wavepacket dynamics, including the variational and nonvariational thawed and frozen Gaussian approximations, and special limits of these methods based on the global harmonic, local harmonic, single-Hessian, local cubic, and single quartic approximations for the potential energy. Without substantially increasing the cost, the proposed single quartic variational thawed Gaussian wavepacket dynamics improves the accuracy over the local cubic approximation and, at the same time, conserves both the effective energy and symplectic structure, in contrast to the much more expensive local quartic approximation. Most results are presented in both Heller's and Hagedorn's parametrizations of the Gaussian wavepacket.

研究动机与目标

• 将多种高斯波包动力学方法统一于一个基于状态依赖二次势的非线性薛定谔方程的理论框架下。
• 推导高斯波包参数的一般运动方程，并证明基本守恒律：规范性、有效能量与辛结构。
• 开发高阶几何积分器，以高效且精确地求解非线性薛定谔方程。
• 提出并分析一种新型的‘单四次’变分高斯近似，其精度优于局部立方方法，同时保持能量与辛结构守恒。
• 通过将其应用于变分与非变分的非冻结与冻结高斯近似（包括全局谐振、局部谐振、单Hessian、局部立方与局部四次极限），展示该框架的普适性。

提出的方法

• 将高斯波包动力学表述为一个非线性时变薛定谔方程的精确解，其有效哈密顿量的参数（二次势系数）依赖于波包的瞬时状态。
• 通过将时变变分原理应用于高斯试探波函数，推导出高斯波包中心、宽度矩阵、相位与动量的一般运动方程。
• 通过分析非线性系统的基本几何与代数性质，证明时间可逆性、规范性守恒，以及有效能量与辛结构的守恒。
• 引入高阶几何积分器，以保持辛结构与有效能量，确保数值模拟中长期的稳定性和精度。
• 将该框架应用于多种近似方法：变分与非变分的非冻结与冻结高斯，以及通过谐振、立方与四次势近似得到的极限形式。
• 通过在局部立方近似中增加一个四阶导数项，提出一种新型的‘单四次’变分方法，保持计算效率的同时提升精度，并维持几何不变量。

实验结果

研究问题

• RQ1 多种高斯波包动力学方法（如Heller的非冻结高斯与Coalson-Karplus的变分高斯）能否统一于一个具有状态依赖势的非线性薛定谔方程框架下？
• RQ2 在此非线性框架中，高斯波包参数的一般运动方程为何？其守恒的几何性质（如规范性、能量、辛结构）为何？
• RQ3 如何构造高阶几何积分器，以在非线性薛定谔方程的数值求解中保持辛结构与有效能量？
• RQ4 能否设计一种新的变分高斯近似，使其精度优于局部立方近似，同时不牺牲能量或辛结构守恒？
• RQ5 在高斯波包动力学中，包含更高阶势能导数（如四阶导数）对精度与守恒性质有何影响？

主要发现

• 在非线性薛定谔方程框架下，有效能量被精确守恒，该框架推广并统一了多种高斯波包动力学方法。
• 运动方程保持了辛结构，确保了数值模拟中长期的稳定性和精度。
• 所提出的‘单四次’变分高斯近似在保持有效能量与辛结构守恒的同时，显著提升了相对于局部立方近似的精度。
• 该框架揭示，尽管单Hessian近似形式简单，却能精确守恒有效能量，并在扩展相空间中具有哈密顿结构。
• 该通用理论在Heller与Hagedorn的高斯波包参数化形式下均成立，使其在不同形式间具有广泛适用性。
• 推导出高阶几何积分器，可保持辛结构与有效能量，从而实现高斯波包的高效且精确的长时间演化。