An algebraic method for determining Magnus propagators of arbitrary order is described. These propagators preserve the wronskian and therefore the linear independence of the solutions and can be used for numerical solutions of multichannel time-dependent and time-independent Schrödinger equations. As an example the propagator for the single-channel Schrödinger equation is derived analytically for quadratic reference potentials up to seventh order in the stepsize. Recent results obtained by purely numerical procedures are verified, corrected, improved and considerably extended.
An algebraic method for determining Magnus propagators of arbitrary order is described. These propagators preserve the wronskian and therefore the linear independence of the solutions and can be used for numerical solutions of multichannel time-dependent and time-independent Schrödinger equations. As an example the propagator for the single-channel Schrödinger equation is derived analytically for quadratic reference potentials up to seventh order in the stepsize. Recent results obtained by purely numerical procedures are verified, corrected, improved and considerably extended.
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 consider the numerical integration of the Schrödinger equation with a time-dependent Hamiltonian given as the sum of the kinetic energy and a time-dependent potential. Commutator-free (CF) propagators are exponential propagators that have shown to be highly efficient for general time-dependent Hamiltonians. We propose new CF propagators that are tailored for Hamiltonians of the said structure, showing a considerably improved performance. We obtain new fourth- and sixth-order CF propagators as well as a novel sixth-order propagator that incorporates a double commutator that only depends on coordinates, so this term can be considered as cost-free. The algorithms require the computation of the action of exponentials on a vector similar to the well-known exponential midpoint propagator, and this is carried out using the Lanczos method. We illustrate the performance of the new methods on several numerical examples.
We derive and discuss a formula, due to Magnus, for the exponential representation of the operator solution to Schrödinger's equation when the Hamiltonian is time dependent. The formula gives a unitary time-displacement operator in every order of approximation. We study the usefulness of the first- and second-order approximations for the kind of problem posed by the semiclassical theory of inelastic collisions, basing our discussion on two exactly soluble two-state problems. The algebraic structure of the Magnus formula is in itself useful; to illustrate this, we solve exactly the problems of the linearly forced harmonic oscillator and the harmonic oscillator with time-dependent force constant.
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.