Abstract:As a basic principle, benefits of adaptive discretisations are an improved balance between required accuracy and efficiency as well as an enhancement of the reliability of numerical computations. In this work, the capacity of locally adaptive space and time discretisations for the numerical solution of low-dimensional nonlinear Schrödinger equations is investigated. The considered model equation is related to the time-dependent Gross-Pitaevskii equation arising in the description of Bose-Einstein condensates in dilute gases. The performance of the Fourier-pseudo spectral method constrained to uniform meshes versus the locally adaptive finite element method and of higher-order exponential operator splitting methods with variable time stepsizes is studied. Numerical experiments confirm that a local time stepsize control based on a posteriori local error estimators or embedded splitting pairs, respectively, is effective in different situations with an enhancement either in efficiency or reliability. As expected, adaptive time-splitting schemes combined with fast Fourier transform techniques are favourable regarding accuracy and efficiency when applied to Gross-Pitaevskii equations with a defocusing nonlinearity and a mildly varying regular solution. However, the numerical solution of nonlinear Schrödinger equations in the semi-classical regime becomes a demanding task. Due to the highly oscillatory and nonlinear nature of the problem, the spatial mesh size and the time increments need to be of the size of the decisive parameter [Formula: see text], especially when it is desired to capture correctly the quantitative behaviour of the wave function itself. The required high resolution in space constricts the feasibility of numerical computations for both, the Fourier pseudo-spectral and the finite element method. Nevertheless, for smaller parameter values locally adaptive time discretisations facilitate to determine the time stepsizes sufficiently small in order that the numerical approximation captures correctly the behaviour of the analytical solution. Further illustrations for Gross-Pitaevskii equations with a focusing nonlinearity or a sharp Gaussian as initial condition, respectively, complement the numerical study.

Adaptive Time Stepping for a Two-Time Integro-Differential Equation in Non-Equilibrium Quantum Dynamics

Using Dynamic Mode Decomposition to Predict the Dynamics of a Two-Time Non-Equilibrium Green's Function

Using dynamic mode decomposition to predict the dynamics of a two-time non-equilibrium Green’s function

Computing the many-body Green's function with adaptive variational quantum dynamics

Non-equilibrium Green's function approach to inhomogeneous quantum many-body systems using the Generalized Kadanoff Baym Ansatz

Dynamic Mode Decomposition for Extrapolating Nonequilibrium Green?s-function Dynamics

Adaptive time-stepping for aggregation-shattering kinetics

Making Trotterization adaptive and energy-self-correcting for NISQ devices and beyond

Excitations and spectra from equilibrium real-time Green's functions

Analyzing and predicting non-equilibrium many-body dynamics via dynamic mode decomposition

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

Reversible time-step adaptation for the integration of few-body systems

Telescopic projective integration for kinetic equations with multiple relaxation times

Quantum trajectories for time-local non-Lindblad master equations

Asymptotic Behaviour of Time Stepping Methods for Phase Field Models

Time-dependent resonant tunneling transport: Keldysh and Kadanoff-Baym nonequilibrium Green's functions in an analytically soluble problem

Single-seed generation of Brownian paths and integrals for adaptive and high order SDE solvers

A numerical study of adaptive space and time discretisations for Gross-Pitaevskii equations

Few-particle quantum dynamics---comparing Nonequilibrium Green's functions with the generalized Kadanoff-Baym ansatz to density operator theory

Heisenberg formulation of adiabatic elimination for open quantum systems with two time-scales

Solving the Kadanoff-Baym equations for inhomogenous systems: Application to atoms and molecules