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.

Approximate Solution of the Classical Liouville Equation Using Gaussian Phase Packet Dynamics: Application to Enhanced Equilibrium Averaging and Global Optimization

Lyapunov Function Method for Analyzing Stability of Quasi-Hamiltonian Systems under Combined Gaussian and Poisson White Noise Excitations

Simulated Annealing Using Coarse-Grained Classical Dynamics - Smoluchowski Dynamics in the Gaussian Density Approximation

Simulated Annealing Using the Classical Density Distribution

Application of Grey System Theory in Approximate Calculation of Wave Packet Evolution

Approximate Exponential Integrators for Time-Dependent Equation-of-Motion Coupled Cluster Theory

Linearized Semiclassical Initial Value Time Correlation Functions Using the Thermal Gaussian Approximation: Applications to Condensed Phase Systems.

Stationary State Distribution and Efficiency Analysis of the Langevin Equation Via Real or Virtual Dynamics.

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

Augmented Variational Superposed Gaussian Approximation for Langevin Equations with Rational Polynomial Functions

Stochastic Integration of the Cahn-Hilliard Phase Field Equations

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

Liouville models of particle-laden flow

Lagrangian Liouville models of multiphase flows with randomly forced inertial particles

Path Integral Liouville Dynamics for Thermal Equilibrium Systems.

An approach for generating trajectory-based dynamics which conserves the canonical distribution in the phase space formulation of quantum mechanics. I. Theories

An Energetic Variational Approach for the Cahn-Hilliard Equation with Dynamic Boundary Conditions: Derivation and Analysis

Time-Dependent Quantum Transport: An Efficient Method Based On Liouville-Von-Neumann Equation For Single-Electron Density Matrix

Treatment of long-range interactions arising in the Enskog-Vlasov description of dense fluids

Diffuse Interface Methods for Multiple Phase Materials: An Energetic Variational Approach

GNN-Assisted Phase Space Integration with Application to Atomistics