Numerical integration of quantum time evolution in a curved manifold

Abstract

The numerical integration of the Schrödinger equation by discretization of time is explored for the curved manifolds arising from finite representations based on evolving basis states. In particular, the unitarity of the evolution is assessed, in the sense of the conservation of mutual scalar products in a set of evolving states, and with them the conservation of orthonormality and particle number. Although the adequately represented equation is known to give rise to unitary evolution in spite of curvature, discretized integrators easily break that conservation, thereby deteriorating their stability. The Crank Nicolson algorithm, which offers unitary evolution in Euclidian spaces independent of time-step size $\mathrm{d}t$, can be generalised to curved manifolds in different ways. Here we compare a previously proposed algorithm that is unitary by construction, albeit integrating the wrong equation, with a faithful generalisation of the algorithm, which is, however, not strictly unitary for finite $\mathrm{d}t$.

arXiv comment12 pages, 5 figures. Accepted in Phys. Rev. Research

References in corpus (14)

• The SIESTA method for ab initio order-N materials simulation
• Efficient index handling of multidimensional periodic boundary conditions
• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• QuantumATK: An integrated platform of electronic and atomic-scale modelling tools
• SIESTA: recent developments and applications
• Calculation of the optical response of C60 and Na8 using time-dependent density functional theory and local orbitals
• Electronic stopping power in gold: The role of d electrons and the H/He anomaly
• Electronic stopping power in a narrow band gap semiconductor from first principles
• The CECAM Electronic Structure Library and the modular software development paradigm
• Anisotropy of electronic stopping power in graphite
• Quantum mechanics in an evolving Hilbert space
• Exponential Integrators in Time-Dependent Density Functional Calculations
• Efficient ab initio calculation of electronic stopping in disordered systems via geometry pre-sampling: application to liquid water
• Real-Time Time-Dependent Density Functional Theory within FHI-aims

Cited by in corpus (1)

• Manifold curvature and Ehrenfest forces with a moving basis