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$.

12 pages, 5 figures. Accepted in Phys. Rev. Research

References in corpus (7)

• SIESTA: recent developments and applications
• 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
• 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