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.

Domain decomposition algorithms for the two dimensional nonlinear Schr{ö}dinger equation and simulation of Bose-Einstein condensates

Super-localised wave function approximation of Bose-Einstein condensates

New optimized Schwarz algorithms for one dimensional Schrödinger equation with general potential

An efficient and conservative compact finite difference scheme for the coupled Gross-Pitaevskii equations describing spin-1 Bose-Einstein condensate.

Domain Decomposition Method for Poisson--Boltzmann Equations based on Solvent Excluded Surface

Computation of excited states for the nonlinear Schr{ö}dinger equation: numerical and theoretical analysis

A two level approach for simulating Bose-Einstein condensates by localized orthogonal decomposition

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

An overlapping domain decomposition splitting algorithm for stochastic nonlinear Schroedinger equation

A Schur Complement Preconditioner for Scalable Parallel Fluid Simulation

Quasi-optimal domain decomposition method for neural network-based computation of the time-dependent Schrödinger equation

On Deterministic Numerical Methods for the Quantum Boltzmann-Nordheim Equation. I. Spectrally Accurate Approximations, Bose-Einstein Condensation, Fermi-Dirac Saturation

Domain and Subband decomposition approach for 2 D simulation of quantum transport phenomena ∗

On an interior-exterior nonoverlapping domain decomposition method for the Poisson--Boltzmann equation

Quantum kinetic theory: modelling and numerics for Bose-Einstein condensation

A Regularized Newton Method for Computing Ground States of Bose–Einstein Condensates

Numerical Solution to Coupled Nonlinear Schrödinger Equations on Unbounded Domains

Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations with absorption

Newton-based alternating methods for the ground state of a class of multi-component Bose-Einstein condensates