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.

High-order geometric integrators for the local cubic variational Gaussian wavepacket dynamics

High-order geometric integrators for the local cubic variational Gaussian wavepacket dynamics

High-order geometric integrators for the variational Gaussian approximation

Family of Gaussian wavepacket dynamics methods from the perspective of a nonlinear Schrödinger equation

Multi-symplectic Variational Integrators for Nonlinear Schrödinger Equations with Variable Coefficients

Arbitrary-order Energy-Preserving Exponential Integrators for the Cubic Schrödinger Equation

Collocation in Conjunction with the Local Variational Iteration Method

Surrogate Lagrangians for Variational Integrators: High Order Convergence with Low Order Schemes

Efficient integration of the variational equations of multi-dimensional Hamiltonian systems: Application to the Fermi-Pasta-Ulam lattice

New Way to Construct High Order Hamiltonian Variational Integrators

High-order conservative and accurately dissipative numerical integrators via auxiliary variables

Nonstandard Finite Difference Variational Integrators for Nonlinear Schrödinger Equation with Variable Coefficients

Degenerate Variational Integrators for Magnetic Field Line Flow and Guiding Center Trajectories

High-order linearly implicit structure-preserving exponential integrators for the nonlinear Schr\"odinger equation

Efficient evaluation of two-center Gaussian integrals in periodic systems

Efficient Calculation of Derivatives of Integrals in a Basis of Non-Separable Gaussians Through Exploitation of Sparsity

High order variational integrators in the optimal control of mechanical systems

General Local Energy-Preserving Integrators for Solving Multi-Symplectic Hamiltonian PDEs.

Geometric integrators and the Hamiltonian Monte Carlo method

Nonlinear Geometric Optics Method Based Multi-Scale Numerical Schemes for Highly-Oscillatory Transport Equations

An accelerated linear method for optimizing non-linear wavefunctions in variational Monte Carlo