© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2312.01115 · The Journal of Physical Chemistry A

Simple and General Unitarity Conserving Numerical Real-Time Propagators of the Time-Dependent Schrödinger Equation Based on Magnus Expansion

Taner M. Ture, Seogjoo Jang

Published Apr 2, 2024• 3 citations

Abstract

Magnus expansion (ME) provides a general way to expand the real-time propagator of a time-dependent Hamiltonian within the exponential such that the unitarity is satisfied at any order. We use this property and explicit integration of Lagrange interpolation formulas for the time-dependent Hamiltonian within each time interval and derive approximations that preserve unitarity for the differential time evolution operators of general time-dependent Hamiltonians. The resulting second-order approximation is the same as using the average of Hamiltonians for two end points of time. We identify three fourth-order approximations involving commutators of Hamiltonians at different times and also derive a sixth-order expression. A test of these approximations along with other available expressions for a two-state time-dependent Hamiltonian with sinusoidal time dependences provides information on the relative performance of these approximations and suggests that the derived expressions can serve as useful numerical tools for time evolution in time-resolved spectroscopy, quantum control, quantum sensing, real-time ab initio quantum dynamics, and open system quantum dynamics.

DOI: 10.1021/acs.jpca.3c07866

Related papers

The Magnus Expansion for Non-Hermitian Hamiltonians

Yair Mulian

The Magnus expansion offers a powerful geometric method to express a time-ordered exponential operator solution, typical of time-dependent linear differential equations, as an ordinary operatorial exponential. This representation possesses advantageous theoretical and qualitative features, such as preserving the underlying Lie algebraic structure of the system while remaining an exact solution to the original differential equation. For any finite-dimensional Hermitian Hamiltonian governing a quantum system, the standard Magnus expansion guarantees a manifestly unitary representation at every finite order of truncation. However, this foundational property is no longer preserved if the underlying Hamiltonian is non-Hermitian. In this work, we derive a generalized, unitarized expansion framework that systematically restores and maintains the property of geometric unitarity for all bounded finite-dimensional non-Hermitian Hamiltonians.

arXiv:2505.09559•Modern mathematical physics•Jul 24, 2026•0 citations

A pedagogical approach to the Magnus expansion

Sergio Blanes, Fernando Casas, José Ángel Oteo +1 more

Time-dependent perturbation theory as a tool to compute approximate solutions of the Schrödinger equation does not preserve unitarity. Here we present, in a simple way, how the Magnus expansion (also known as exponential perturbation theory ) provides such unitary approximate solutions. The purpose is to illustrate the importance and consequences of such a property. We suggest that the Magnus expansion may be introduced to students in advanced courses of quantum mechanics.

European Journal of Physics•Jun 16, 2010•80 citations

EXPONENTIAL PROPAGATORS (INTEGRATORS) FOR THE TIME-DEPENDENT SCHRÖDINGER EQUATION

André D. Bandrauk, Huizhong Lu

The time-dependent Schrödinger Equation (TDSE) is a parabolic partial differential equation (PDE) comparable to a diffusion equation but with imaginary time. Due to its first order time derivative, exponential integrators or propagators are natural methods to describe evolution in time of the TDSE, both for time-independent and time-dependent potentials. Two splitting methods based on Fer and/or Magnus expansions allow for developing unitary factorizations of exponentials with different accuracies in the time step △t. The unitary factorization of exponentials to high order accuracy depends on commutators of kinetic energy operators with potentials. Fourth-order accuracy propagators can involve negative or complex time steps, or real time steps only but with gradients of potentials, i.e. forces. Extending the propagators of TDSE's to imaginary time allows to also apply these methods to classical many-body dynamics, and quantum statistical mechanics of molecular systems.

Journal of Theoretical and Computational Chemistry•Apr 23, 2013•70 citations

Propagators for Quantum-Classical Models: Commutator-Free Magnus Methods

Adrián Gómez Pueyo, Sergio Blanes, Alberto Castro

We consider the numerical propagation of models that combine both quantum and classical degrees of freedom, usually, electrons and nuclei, respectively. We focus, in our computational examples, on the case in which the quantum electrons are modeled with time-dependent density-functional theory, although the methods discussed below can be used with any other level of theory. Often, for these so-called quantum-classical molecular dynamics models, one uses some propagation technique to deal with the quantum part and a different one for the classical equations. While the resulting procedure may, in principle, be consistent, it can however spoil some of the properties of the methods, such as the accuracy order with respect to the time step or the preservation of the geometrical structure of the equations. Few methods have been developed specifically for hybrid quantum-classical models. We propose using the same method for both the quantum and classical particles, in particular, one family of techniques that proves to be very efficient for the propagation of Schrödinger-like equations: the (quasi)-commutator free Magnus expansions. These have been developed, however, for linear systems, yet our problem is nonlinear: formally, the full quantum-classical system can be rewritten as a nonlinear Schrödinger equation, i.e., one in which the Hamiltonian depends on the system itself. The Magnus expansion algorithms for linear systems require the application of the Hamiltonian at intermediate points in a given propagating interval. For nonlinear systems, this poses a problem as this Hamiltonian is unknown due to its dependence on the state. We approximate it by employing a higher order extrapolation using previous steps as input. The resulting technique can then be regarded as a multistep technique or, alternatively, as a predictor corrector formula.

Journal of Chemical Theory and Computation•Jan 30, 2020•6 citations

Exponential Perturbative Expansions and Coordinate Transformations

Ana Arnal, Fernando Casas, Cristina Chiralt

We propose a unified approach for different exponential perturbation techniques used in the treatment of time-dependent quantum mechanical problems, namely the Magnus expansion, the Floquet–Magnus expansion for periodic systems, the quantum averaging technique, and the Lie–Deprit perturbative algorithms. Even the standard perturbation theory fits in this framework. The approach is based on carrying out an appropriate change of coordinates (or picture) in each case, and it can be formulated for any time-dependent linear system of ordinary differential equations. All of the procedures (except the standard perturbation theory) lead to approximate solutions preserving by construction unitarity when applied to the time-dependent Schrödinger equation.

Mathematical and Computational Applications•Aug 13, 2020•6 citations

Exponential perturbative expansions and coordinate transformations

Ana Arnal, Fernando Casas, Cristina Chiralt

We propose a unified approach for different exponential perturbation techniques used in the treatment of time-dependent quantum mechanical problems, namely the Magnus expansion, the Floquet--Magnus expansion for periodic systems, the quantum averaging technique and the Lie--Deprit perturbative algorithms. Even the standard perturbation theory fits in this framework. The approach is based on carrying out an appropriate change of coordinates (or picture) in each case, and can be formulated for any time-dependent linear system of ordinary differential equations. All the procedures (except the standard perturbation theory) lead to approximate solutions preserving by construction unitarity when applied to the time-dependent Schrödinger equation.

arXiv:2401.12955•Jan 23, 2024•0 citations

Magnus expansion (ME) provides a general way to expand the real-time propagator of a time-dependent Hamiltonian within the exponential such that the unitarity is satisfied at any order. We use this property and explicit integration of Lagrange interpolation formulas for the time-dependent Hamiltonian within each time interval and derive approximations that preserve unitarity for the differential time evolution operators of general time-dependent Hamiltonians. The resulting second-order approximation is the same as using the average of Hamiltonians for two end points of time. We identify three fourth-order approximations involving commutators of Hamiltonians at different times and also derive a sixth-order expression. A test of these approximations along with other available expressions for a two-state time-dependent Hamiltonian with sinusoidal time dependences provides information on the relative performance of these approximations and suggests that the derived expressions can serve as useful numerical tools for time evolution in time-resolved spectroscopy, quantum control, quantum sensing, real-time ab initio quantum dynamics, and open system quantum dynamics.

The Magnus expansion offers a powerful geometric method to express a time-ordered exponential operator solution, typical of time-dependent linear differential equations, as an ordinary operatorial exponential. This representation possesses advantageous theoretical and qualitative features, such as preserving the underlying Lie algebraic structure of the system while remaining an exact solution to the original differential equation. For any finite-dimensional Hermitian Hamiltonian governing a quantum system, the standard Magnus expansion guarantees a manifestly unitary representation at every finite order of truncation. However, this foundational property is no longer preserved if the underlying Hamiltonian is non-Hermitian. In this work, we derive a generalized, unitarized expansion framework that systematically restores and maintains the property of geometric unitarity for all bounded finite-dimensional non-Hermitian Hamiltonians.

The time-dependent Schrödinger Equation (TDSE) is a parabolic partial differential equation (PDE) comparable to a diffusion equation but with imaginary time. Due to its first order time derivative, exponential integrators or propagators are natural methods to describe evolution in time of the TDSE, both for time-independent and time-dependent potentials. Two splitting methods based on Fer and/or Magnus expansions allow for developing unitary factorizations of exponentials with different accuracies in the time step △t. The unitary factorization of exponentials to high order accuracy depends on commutators of kinetic energy operators with potentials. Fourth-order accuracy propagators can involve negative or complex time steps, or real time steps only but with gradients of potentials, i.e. forces. Extending the propagators of TDSE's to imaginary time allows to also apply these methods to classical many-body dynamics, and quantum statistical mechanics of molecular systems.

We consider the numerical propagation of models that combine both quantum and classical degrees of freedom, usually, electrons and nuclei, respectively. We focus, in our computational examples, on the case in which the quantum electrons are modeled with time-dependent density-functional theory, although the methods discussed below can be used with any other level of theory. Often, for these so-called quantum-classical molecular dynamics models, one uses some propagation technique to deal with the quantum part and a different one for the classical equations. While the resulting procedure may, in principle, be consistent, it can however spoil some of the properties of the methods, such as the accuracy order with respect to the time step or the preservation of the geometrical structure of the equations. Few methods have been developed specifically for hybrid quantum-classical models. We propose using the same method for both the quantum and classical particles, in particular, one family of techniques that proves to be very efficient for the propagation of Schrödinger-like equations: the (quasi)-commutator free Magnus expansions. These have been developed, however, for linear systems, yet our problem is nonlinear: formally, the full quantum-classical system can be rewritten as a nonlinear Schrödinger equation, i.e., one in which the Hamiltonian depends on the system itself. The Magnus expansion algorithms for linear systems require the application of the Hamiltonian at intermediate points in a given propagating interval. For nonlinear systems, this poses a problem as this Hamiltonian is unknown due to its dependence on the state. We approximate it by employing a higher order extrapolation using previous steps as input. The resulting technique can then be regarded as a multistep technique or, alternatively, as a predictor corrector formula.

We propose a unified approach for different exponential perturbation techniques used in the treatment of time-dependent quantum mechanical problems, namely the Magnus expansion, the Floquet–Magnus expansion for periodic systems, the quantum averaging technique, and the Lie–Deprit perturbative algorithms. Even the standard perturbation theory fits in this framework. The approach is based on carrying out an appropriate change of coordinates (or picture) in each case, and it can be formulated for any time-dependent linear system of ordinary differential equations. All of the procedures (except the standard perturbation theory) lead to approximate solutions preserving by construction unitarity when applied to the time-dependent Schrödinger equation.

We propose a unified approach for different exponential perturbation techniques used in the treatment of time-dependent quantum mechanical problems, namely the Magnus expansion, the Floquet--Magnus expansion for periodic systems, the quantum averaging technique and the Lie--Deprit perturbative algorithms. Even the standard perturbation theory fits in this framework. The approach is based on carrying out an appropriate change of coordinates (or picture) in each case, and can be formulated for any time-dependent linear system of ordinary differential equations. All the procedures (except the standard perturbation theory) lead to approximate solutions preserving by construction unitarity when applied to the time-dependent Schrödinger equation.