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.

Frozen Gaussian Sampling for Scalar Wave Equations

Frozen Gaussian approximation for the fractional Schrödinger equation

Frozen Gaussian approximation for high frequency wave propagation in periodic media

Convergence of frozen Gaussian approximation for high frequency wave propagation

Seismic modeling using the frozen Gaussian approximation

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

Global Geometrical Optics Method

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

Gaussian Wave Packet Transform Based Numerical Scheme for the Semi-Classical Schrödinger Equation with Random Inputs

Newton-type Gauss–Seidel Lax–Friedrichs High-Order Fast Sweeping Methods for Solving Generalized Eikonal Equations at Large-Scale Discretization

Numerical Simulation of Non-Gaussian Wave Elevation and Kinematics Based on Two-Dimensional Fourier Transform

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

Solving Caldeira-Leggett Model by Inchworm Method with Frozen Gaussian Approximation

A Grid-Free Fluid Solver based on Gaussian Spatial Representation

Gaussian Beam Methods for the Schrodinger Equation in the Semi-classical Regime

Extended WKB analysis for the linear vectorial wave equation in the high-frequency regime

Mitigating the Cycle-Skipping of Full-Waveform Inversion by Random Gradient Sampling

A Fast and Accurate Solver for the Fractional Fokker-Planck Equation with Dirac-Delta Initial Conditions

Wavefield Separation Algorithm of Helmholtz Theory-Based Discontinuous Galerkin Method Using Unstructured Meshes on GPU

Fourier pseudospectral methods for the spatial variable-order fractional wave equations

Non-Gaussian random wave simulation by two-dimensional Fourier transform and linear oscillator responselo Morison force