Very accurate time propagation of coupled Schrödinger equations for femto- and attosecond physics and chemistry, with C++ source code
arXiv:2307.02479 · doi:10.1016/j.cpc.2023.108839
Abstract
In this article, I present a very fast and high-precision (up to 33 decimal places) C++ implementation of the semi-global time propagation algorithm for a system of coupled Schrödinger equations with a time-dependent Hamiltonian. It can be used to describe time-dependent processes in molecular systems after excitation by femto- and attosecond laser pulses. It also works with an arbitrary user supplied Hamiltonian and can be used for nonlinear problems. The semi-global algorithm is briefly presented, the C++ implementation is described and five sample simulations are shown. The accompanying C++ source code package is included. The high precision benchmark (long double and float128) shows the estimated calculation costs. The presented method turns out to be faster and more accurate than the global Chebyshev propagator.
arXiv comment22 pages, 9 figures, 8 tables, 4 code listings
References in corpus (11)
- Numerical Solution of the Gross-Pitaevskii Equation for Bose-Einstein Condensation
- Massive 3-loop Feynman diagrams reducible to SC* primitives of algebras of the sixth root of unity
- WavePacket: A Matlab package for numerical quantum dynamics. I: Closed quantum systems and discrete variable representations
- WavePacket: A Matlab package for numerical quantum dynamics. II: Open quantum systems, optimal control, and model reduction
- Semi-global approach for propagation of the time-dependent Schrödinger equation for time-dependent and nonlinear problems
- A Chebychev propagator with iterative time ordering for explicitly time-dependent Hamiltonians
- Extended Hylleraas three-electron integral
- New, Highly Accurate Propagator for the Linear and Nonlinear Schrödinger Equation
- Long-range asymptotics of exchange energy in the hydrogen molecule
- Implementation of high-precision computation capabilities into the open-source dynamic simulation framework YADE
- Optimization of the Femtosecond Laser Impulse for Excitation and the Spin-Orbit Mediated Dissociation in the NaRb Dimer