Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1301.7596

Adaptive propagation of quantum few-body systems with time-dependent Hamiltonians

Jonas C. Cremon

Published Jan 31, 2013• 0 citations
arXivHTMLPDF

Abstract

In this study, a variety of methods are tested and compared for the numerical solution of the Schrödinger equation for few-body systems with explicitely time-dependent Hamiltonians, with the aim to find the optimal one. The configuration interaction method, generally applied to find stationary eigenstates accurately and without approximations to the wavefunction's structure, may also be used for the time-evolution, which results in a large linear system of ordinary differential equations. The large basis sizes typically present when the configuration interaction method is used calls for efficient methods for the time evolution. Apart from efficiency, adaptivity (in the time domain) is the other main focus in this study, such that the time step is adjusted automatically given some requested accuracy. A method is suggested here, based on an exponential integrator approach, combined with different ways to implement the adaptivity, which was found to be faster than a broad variety of other methods that were also considered.

DOI
10.48550/arxiv.1301.7596

Related papers

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

Dynamical Simulations of Schrödinger's Equation via Rank-Adaptive Tensor Decompositions

N. Anders Petersson, Chase Hodges-Heilmann, Stefanie Günther

We study low-rank tensor methods for the numerical solution of Schrödinger's equation with time-independent and explicitly time-dependent Hamiltonians, motivated by large-scale simulations of many-body quantum systems and quantum computing devices subject to time-dependent control pulses. We outline the recent application of the "basis update and Galerkin" (BUG) method for tensor trains, and describe the established TDVP and TDVP-2 algorithms based on the time-dependent variational principle. For comparison, we also consider the BUG method in the Tucker format. All these approaches enable memory efficient representations of partially entangled quantum states and thereby mitigate the exponential cost of conventional state-vector formulations. The rank-adaptivity relies on the truncated singular value decomposition, in which the rank of a matrix is reduced by setting its smallest singular values to zero, based on a threshold parameter that controls the truncation error. Numerical experiments on representative time-independent and time-dependent Hamiltonian models quantify the tradeoff between accuracy and compression across methods, with particular attention to the interplay between the time-step and the truncation threshold, and how the computational effort scales with the number of sub-systems in the quantum system.

quant-pharXiv:2603.13990•Mar 17, 2026•0 citations

Approximate Exponential Integrators for Time-Dependent Equation-of-Motion Coupled Cluster Theory

David B. Williams‐Young, Stephen H. Yuwono, A. Eugene DePrince +1 more

High Resolution Image Download MS PowerPoint Slide With a growing demand for time-domain simulations of correlated many-body systems, the development of efficient and stable integration schemes for the time-dependent Schrödinger equation is of keen interest in modern electronic structure theory. In this work, we present two approaches for the formation of the quantum propagator for time-dependent equation-of-motion coupled cluster theory based on the Chebyshev and Arnoldi expansions of the complex, nonhermitian matrix exponential, respectively. The proposed algorithms are compared with the short-iterative Lanczos method of Cooper et al. [ J. Phys. Chem. A 2021 125, 5438–5447], the fourth-order Runge–Kutta method, and exact dynamics for a set of small but challenging test problems. For each of the cases studied, both of the proposed integration schemes demonstrate superior accuracy and efficiency relative to the reference simulations.

arXiv:2305.07592•Journal of Chemical Theory and Computation•Dec 12, 2023•5 citations

Efficient time-evolution of matrix product states using average Hamiltonians

Belal Abouraya, Jirawat Saiphet, Fedor Jelezko +1 more

Simulating quantum many-body systems (QMBS) is one of the long-standing, highly non-trivial challenges in condensed matter physics and quantum information due to the exponentially growing size of the system's Hilbert space. To date, tensor networks have been an essential tool for studying such quantum systems, owing to their ability to efficiently capture the entanglement properties of the systems they represent. One of the well-known tensor network architectures, namely matrix product states (MPS), is the standard method for simulating one-dimensional QMBS. Here, we propose a simple, yet efficient, method to augment the already available MPS algorithms to simulate the dynamics of time-dependent Hamiltonians with better accuracy and a faster convergence rate, giving a second-order convergence compared to the first-order convergence of the standard method. We apply our proposed method to simulate the dynamics of a chain of single spins associated with nitrogen-vacancy color centers in diamonds, which has potential applications for practical and scalable quantum technologies, and find that our method improves the average error for a system of few NV centers by a factor of about 1000 for moderate step sizes. Our work paves the way for efficient simulation of QMBS under the influence of time-dependent Hamiltonians.

quant-pharXiv:2602.04955•Feb 6, 2026•0 citations

Numerical quantum propagation with time-dependent Hamiltonian

Wusheng Zhu, Xinsheng Zhao

A numerical method was proposed to propagate the quantum system with a time-dependent Hamiltonian. The propagator is composed of the time-independent-Hamiltonian-type multisplit forms, which becomes possible by realizing that the evolution operator has a close link with the time-shift operator. The present algorithm can achieve a high order of accuracy through the recurrence of multisplit forms with predetermined coefficients, and it has the advantage that only the Hamiltonian itself at different time instants is involved, while other currently available schemes need to perform either a time integration or a time differential of the Hamiltonian. This essential difference has made our proposal highly competitive. A one-dimensional model system is considered to verify the accuracy and efficiency of the numerical scheme.

The Journal of Chemical Physics•Dec 1, 1996•8 citations

The Operator Method of Solution of the Non-Relativistic Schrödinger Equation for Few-Body Systems with a Pair-Wise Interaction

V. A. Sidorovich

A new operator approach to the solution of the non-relativistic Schrödinger equation for few-body quantum systems with pair-wise interaction is presented. The problem is solved on the basis of the sequential account of configurations (SAC) method proposed. The wave-function of a few-body system determined with allowance for the configuration interaction (CI-wavefunction) is expressed in terms of the one-configuration wavefunctions with use of the configuration mixing operator (±) introduced in the paper. The examination is carried out with an example of a few-electron atomic system. An analytical expression for the (±) -operator has been represented and an integral equation for the CI-wavefunction has been obtained. It is shown that: (I) the correlation interaction inside an isolated configuration does not destroy the independent particle approximation; (ii) in the framework of the present approach the self-action is excluded. The present approach has enabled to receive an analytical expression for the (±) -operator without applying the non-physical procedure of 'switching-on' the interconfiguration interaction.

Physica Scripta•Sep 1, 1999•0 citations
In this study, a variety of methods are tested and compared for the numerical solution of the Schrödinger equation for few-body systems with explicitely time-dependent Hamiltonians, with the aim to find the optimal one. The configuration interaction method, generally applied to find stationary eigenstates accurately and without approximations to the wavefunction's structure, may also be used for the time-evolution, which results in a large linear system of ordinary differential equations. The large basis sizes typically present when the configuration interaction method is used calls for efficient methods for the time evolution. Apart from efficiency, adaptivity (in the time domain) is the other main focus in this study, such that the time step is adjusted automatically given some requested accuracy. A method is suggested here, based on an exponential integrator approach, combined with different ways to implement the adaptivity, which was found to be faster than a broad variety of other methods that were also considered.
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 study low-rank tensor methods for the numerical solution of Schrödinger's equation with time-independent and explicitly time-dependent Hamiltonians, motivated by large-scale simulations of many-body quantum systems and quantum computing devices subject to time-dependent control pulses. We outline the recent application of the "basis update and Galerkin" (BUG) method for tensor trains, and describe the established TDVP and TDVP-2 algorithms based on the time-dependent variational principle. For comparison, we also consider the BUG method in the Tucker format. All these approaches enable memory efficient representations of partially entangled quantum states and thereby mitigate the exponential cost of conventional state-vector formulations. The rank-adaptivity relies on the truncated singular value decomposition, in which the rank of a matrix is reduced by setting its smallest singular values to zero, based on a threshold parameter that controls the truncation error. Numerical experiments on representative time-independent and time-dependent Hamiltonian models quantify the tradeoff between accuracy and compression across methods, with particular attention to the interplay between the time-step and the truncation threshold, and how the computational effort scales with the number of sub-systems in the quantum system.
High Resolution Image Download MS PowerPoint Slide With a growing demand for time-domain simulations of correlated many-body systems, the development of efficient and stable integration schemes for the time-dependent Schrödinger equation is of keen interest in modern electronic structure theory. In this work, we present two approaches for the formation of the quantum propagator for time-dependent equation-of-motion coupled cluster theory based on the Chebyshev and Arnoldi expansions of the complex, nonhermitian matrix exponential, respectively. The proposed algorithms are compared with the short-iterative Lanczos method of Cooper et al. [ J. Phys. Chem. A 2021 125, 5438–5447], the fourth-order Runge–Kutta method, and exact dynamics for a set of small but challenging test problems. For each of the cases studied, both of the proposed integration schemes demonstrate superior accuracy and efficiency relative to the reference simulations.
Simulating quantum many-body systems (QMBS) is one of the long-standing, highly non-trivial challenges in condensed matter physics and quantum information due to the exponentially growing size of the system's Hilbert space. To date, tensor networks have been an essential tool for studying such quantum systems, owing to their ability to efficiently capture the entanglement properties of the systems they represent. One of the well-known tensor network architectures, namely matrix product states (MPS), is the standard method for simulating one-dimensional QMBS. Here, we propose a simple, yet efficient, method to augment the already available MPS algorithms to simulate the dynamics of time-dependent Hamiltonians with better accuracy and a faster convergence rate, giving a second-order convergence compared to the first-order convergence of the standard method. We apply our proposed method to simulate the dynamics of a chain of single spins associated with nitrogen-vacancy color centers in diamonds, which has potential applications for practical and scalable quantum technologies, and find that our method improves the average error for a system of few NV centers by a factor of about 1000 for moderate step sizes. Our work paves the way for efficient simulation of QMBS under the influence of time-dependent Hamiltonians.
A numerical method was proposed to propagate the quantum system with a time-dependent Hamiltonian. The propagator is composed of the time-independent-Hamiltonian-type multisplit forms, which becomes possible by realizing that the evolution operator has a close link with the time-shift operator. The present algorithm can achieve a high order of accuracy through the recurrence of multisplit forms with predetermined coefficients, and it has the advantage that only the Hamiltonian itself at different time instants is involved, while other currently available schemes need to perform either a time integration or a time differential of the Hamiltonian. This essential difference has made our proposal highly competitive. A one-dimensional model system is considered to verify the accuracy and efficiency of the numerical scheme.
A new operator approach to the solution of the non-relativistic Schrödinger equation for few-body quantum systems with pair-wise interaction is presented. The problem is solved on the basis of the sequential account of configurations (SAC) method proposed. The wave-function of a few-body system determined with allowance for the configuration interaction (CI-wavefunction) is expressed in terms of the one-configuration wavefunctions with use of the configuration mixing operator (±) introduced in the paper. The examination is carried out with an example of a few-electron atomic system. An analytical expression for the (±) -operator has been represented and an integral equation for the CI-wavefunction has been obtained. It is shown that: (I) the correlation interaction inside an isolated configuration does not destroy the independent particle approximation; (ii) in the framework of the present approach the self-action is excluded. The present approach has enabled to receive an analytical expression for the (±) -operator without applying the non-physical procedure of 'switching-on' the interconfiguration interaction.