Solving the Schrödinger eigenvalue problem by the imaginary time propagation technique using splitting methods with complex coefficients

Abstract

The Schrödinger eigenvalue problem is solved with the imaginary time propagation technique. The separability of the Hamiltonian makes the problem suitable for the application of splitting methods. High order fractional time steps of order greater than two necessarily have negative steps and can not be used for this class of diffusive problems. However, there exist methods which use fractional complex time steps with positive real parts which can be used with only a moderate increase in the computational cost. We analyze the performance of this class of schemes and propose new methods which outperform the existing ones in most cases. On the other hand, if the gradient of the potential is available, methods up to fourth order with real and positive coefficients exist. We also explore this case and propose new methods as well as sixth-order methods with complex coefficients. In particular, highly optimized sixth-order schemes for near integrable systems using positive real part complex coefficients with and without modified potentials are presented. A time-stepping variable order algorithm is proposed and numerical results show the enhanced efficiency of the new methods.

12 pages of RevTex4-1, as submitted to journal, revised version

References in corpus (3)

• New families of symplectic splitting methods for numerical integration in dynamical astronomy
• Solving the Schrödinger eigenvalue problem by the imaginary time propagation technique using splitting methods with complex coefficients
• Any order imaginary time propagation method for solving the Schrodinger equation

Cited by in corpus (21)

• Solving the Schrödinger eigenvalue problem by the imaginary time propagation technique using splitting methods with complex coefficients
• Nonlinear dynamics in a synthetic momentum state lattice
• Momentum-Space Josephson Effects
• Solving fractional Schroedinger-type spectral problems: Cauchy oscillator and Cauchy well
• Iterative approaches to the self-consistent nuclear energy density functional problem. Heavy ball dynamics and potential preconditioning
• Non-Hermitian coupled-mode theory for incoherently pumped exciton-polariton condensates
• Iterative Power Algorithm for Global Optimization with Quantics Tensor Trains
• Imaginary-time time-dependent density functional theory and its application for robust convergence of electronic states
• Polarization control of attosecond pulses using bi-chromatic elliptically polarized laser
• A Jacobi spectral method for computing eigenvalue gaps and their distribution statistics of the fractional Schrödinger operator
• Linked and knotted synthetic magnetic fields
• Splitting and composition methods with embedded error estimators
• Diffusion of Quantum Vortices
• Competing role of interactions in synchronization of exciton-polariton condensates
• The Wigner Function of Ground State and One-Dimensional Numerics
• Numerical model of the Gross-Pitaevskii equation for rotating Bose-Einstein condensates using smoothed-particle hydrodynamics
• Nonlocally-induced (quasirelativistic) bound states: Harmonic confinement and the finite well
• Dispersion Managed Elliptical Atomtronics for Interferometry
• Characterizing Fano Resonances During Recollision
• Reduced temporal convergence rates in high-order splitting schemes
• Controlling resonant enhancement in higher-order harmonic generation