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.

Convergence of frozen Gaussian approximation for high frequency wave propagation

Frozen Gaussian approximation for high frequency wave propagation in periodic media

Frozen Gaussian approximation for the fractional Schrödinger equation

Frozen Gaussian Sampling for Scalar Wave Equations

Frozen Gaussian Series Representation of the Imaginary Time Propagator Theory and Numerical Tests

Seismic modeling using the frozen Gaussian approximation

Gaussian beam approach for the boundary value problem of high frequency Helmholtz equation

Global Geometrical Optics Method

Finite-Temperature Correlation Functions Obtained from Combined Real- and Imaginary-Time Propagation of Variational Thawed Gaussian Wavepackets

Global Geometrical Optics Approximation to the High Frequency Helmholtz Equation with Discontinuous Media

A Convergent Multiscale Gaussian-Beam Parametrix for the Wave Equation

A Gaussian Beam Method for High Frequency Solution of Symmetric Hyperbolic Systems with Polarized Waves.

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

The Eulerian Gaussian Beam Method for High Frequency Wave Propagation in the Reduced Momentum Space

An Error Bound for the Time-Sliced Thawed Gaussian Propagation Method

A super-convergent quasi-discontinuous Galerkin method for approximating inertia-gravity and Rossby waves in geophysical flows

Strong convergence rate of finite difference approximations for stochastic cubic Schrödinger equations

Asymptotically Compatible Approximations of Linear Nonlocal Conservation Laws with Variable Horizon

Computing High Frequency Solutions of Symmetric Hyperbolic Systems with Polarized Waves

Superconvergence Analysis of Linear FEM Based on Polynomial Preserving Recovery for Helmholtz Equation with High Wave Number.

A Surface Hopping Gaussian Beam Method for High-Dimensional Transport Systems