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:2403.11327

Quantum dynamics in the self-consistent quadratic approximation

Frank Ernesto Quintela Rodriguez

Published Mar 17, 2024• 0 citations
arXivHTMLPDF

Abstract

A self-consistent quadratic theory is presented to account for nonlinear contributions in quantum dynamics. Evolution equations are shown to depend on higher-order gradients of the Hamiltonian, which are incorporated via their equations of motion or through perturbative calculations. The dynamics is proven trace-preserving, with the Hamiltonian acting as a constant of motion for initial Gaussian states. Nonlinear response functions are calculated perturbatively, and sufficient conditions are provided for the existence of their classical limit.

DOI
10.48550/arxiv.2403.11327

Related papers

Nonlinear differential dynamics of Gaussian States

Julio A. López-Saldívar, Margarita A. Man’ko, Vladimir I. Man’ko

The dynamics of Gaussian states for systems with Hamiltonians, quadratic in the position and momentum operators, gives rise to the definition of a system of nonlinear coupled differential equations for the density matrix parameters of the system states. In this work, we show that, using the Gaussian-state covariance-matrix evolution, it is possible to solve the system of these nonlinear equations. Some examples of applications of the method developed are given for one- and two-dimensional quadratic-system Gaussian states. It is noted that this formalism can also be used to find new constants of motion related to the covariance matrices of quadratic systems.

AIP conference proceedings•Jan 1, 2021•2 citations

Quantum Dynamics

Uri Peskin

Explicitly time-dependent Hamiltonians are ubiquitous in applications of quantum theory. It is therefore necessary to solve the time-dependent Schrödinger equation directly. The system’s dynamics is associated with a unitary time-evolution operator (a propagator), formally given as an infinite Dyson series. Time-dependent observables are invariant under unitary time-dependent transformations, where it is sometimes useful to transform the time-evolution from the states into the corresponding operators. This is carried out in part (in full) by transforming to the interaction (Heisenberg) picture. The corresponding equations of motion for the time-dependent operators are introduced. For quadratic potential energy functions, the time evolution of quantum expectation values coincides with the corresponding classical dynamics. This is demonstrated and analyzed in detail for Gaussian wave packets and a coherent state. Finally, we derive exact and approximate expressions for time-dependent transition probabilities and transition rates between quantum states. The validity of time-dependent perturbation theory is analyzed by comparison to exact dynamics.

Cambridge University Press eBooks•May 11, 2023•0 citations

Chaos in time-dependent variational approximations to quantum dynamics

Fred Cooper, John F. Dawson, Salman Habib +1 more

Dynamical chaos has recently been shown to exist in the Gaussian approximation in quantum mechanics and in the self-consistent mean field approach to studying the dynamics of quantum fields. In this study, we first note that any variational approximation to the dynamics of a quantum system based on the Dirac action principle leads to a classical Hamiltonian dynamics for the variational parameters. Since this Hamiltonian is generically nonlinear and nonintegrable, the dynamics thus generated can be chaotic, in distinction to the exact quantum evolution. We then restrict our attention to a system of two biquadratically coupled quantum oscillators and study two variational schemes, the leading order large-$N$ (four canonical variables) and Hartree (six canonical variables) approximations. The chaos seen in the approximate dynamics is an artifact of the approximations: this is demonstrated by the fact that its onset occurs on the same characteristic time scale as the breakdown of the approximations when compared to numerical solutions of the time-dependent Schr\"odinger equation.

arXiv:quant-ph/9610013•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Feb 1, 1998•46 citations

A semiclassical approximation to quantum dynamics

Yitzhak Weissman, Joshua Jortner

In this paper we advance a new approximation scheme for the dynamics of an N particle system, which is based on the expansion of the Heisenberg equations of motion in powers of h/ in the coherent state representation. The time evolution of the physical observables is generated classically by a Hamiltonian Hs, which is derived from the original Hamiltonian by adding 2N2 virtual particles to the system. These virtual particles are related to the dispersion of the quantum mechanical wave function, and reflect the (time dependent) sensitivity of the instantaneous positions and momenta xi, pi to the variation of the initial conditions. The approximate equations of motion preserve the commutation relations between the coordinates and the momenta and the total (quantum) energy. The present approximation may be adequate for the study of the dynamics of coupled anharmonic oscillators in the quasiperiodic regime.

The Journal of Chemical Physics•Nov 1, 1979•9 citations

Gaussian Stochastic Linearization for Open Quantum Systems Using Quadratic Approximation of Hamiltonians

Igor G. Vladimirov, Ian R. Petersen

This paper extends the energy-based version of the stochastic linearization method, known for classical nonlinear systems, to open quantum systems with canonically commuting dynamic variables governed by quantum stochastic differential equations with non-quadratic Hamiltonians. The linearization proceeds by approximating the actual Hamiltonian of the quantum system by a quadratic function of its observables which corresponds to the Hamiltonian of a quantum harmonic oscillator. This approximation is carried out in a mean square optimal sense with respect to a Gaussian reference quantum state and leads to a self-consistent linearization procedure where the mean vector and quantum covariance matrix of the system observables evolve in time according to the effective linear dynamics. We demonstrate the proposed Hamiltonian-based Gaussian linearization for the quantum Duffing oscillator whose Hamiltonian is a quadro-quartic polynomial of the momentum and position operators. The results of the paper are applicable to the design of suboptimal controllers and filters for nonlinear quantum systems.

arXiv:1202.0946•Feb 5, 2012•7 citations

Tangent equations of motion for nonlinear response functions

Atsushi Ono

Nonlinear response functions, formulated as multipoint correlation functions or Volterra kernels, encode the dynamical and spectroscopic properties of physical systems and underpin a wide range of nonlinear transport and optical phenomena. However, their evaluation rapidly becomes prohibitive at high orders because of combinatorial (often factorial) scaling or severe numerical errors. Here, we establish a systematic and efficient framework to compute nonlinear response functions directly from real-time dynamics, without explicitly constructing multipoint correlators or relying on numerically unstable finite-difference methods for order-resolved extraction. Our approach is based on the Gateaux derivative with respect to the external field in function space, which yields a closed hierarchy of tangent equations of motion (TEOM). Propagating the TEOM alongside the original dynamics isolates each perturbative order with high accuracy, providing a term-by-term decomposition of physical contributions. The computational cost scales exponentially with response order in the fully general setting and reduces to polynomial complexity when all perturbation directions are identical; both regimes avoid the factorial scaling of explicit multipoint-correlator evaluations. We demonstrate the power of TEOM by computing frequency-resolved fifth-order response functions for a solid-state electron model and by obtaining nonlinear response functions up to the 49th order with controlled accuracy in a classical Duffing oscillator. We further show that our time-evolution formulation allows optical conductivities to be evaluated directly while remaining numerically stable even near zero frequency. TEOM can be incorporated seamlessly into existing real-time evolution methods, yielding a general framework for computing nonlinear response functions in quantum and classical dynamical systems.

cond-mat.str-elCross-listedarXiv:2603.22137•Physical Review Research•Mar 24, 2026•0 citations
A self-consistent quadratic theory is presented to account for nonlinear contributions in quantum dynamics. Evolution equations are shown to depend on higher-order gradients of the Hamiltonian, which are incorporated via their equations of motion or through perturbative calculations. The dynamics is proven trace-preserving, with the Hamiltonian acting as a constant of motion for initial Gaussian states. Nonlinear response functions are calculated perturbatively, and sufficient conditions are provided for the existence of their classical limit.
The dynamics of Gaussian states for systems with Hamiltonians, quadratic in the position and momentum operators, gives rise to the definition of a system of nonlinear coupled differential equations for the density matrix parameters of the system states. In this work, we show that, using the Gaussian-state covariance-matrix evolution, it is possible to solve the system of these nonlinear equations. Some examples of applications of the method developed are given for one- and two-dimensional quadratic-system Gaussian states. It is noted that this formalism can also be used to find new constants of motion related to the covariance matrices of quadratic systems.
Explicitly time-dependent Hamiltonians are ubiquitous in applications of quantum theory. It is therefore necessary to solve the time-dependent Schrödinger equation directly. The system’s dynamics is associated with a unitary time-evolution operator (a propagator), formally given as an infinite Dyson series. Time-dependent observables are invariant under unitary time-dependent transformations, where it is sometimes useful to transform the time-evolution from the states into the corresponding operators. This is carried out in part (in full) by transforming to the interaction (Heisenberg) picture. The corresponding equations of motion for the time-dependent operators are introduced. For quadratic potential energy functions, the time evolution of quantum expectation values coincides with the corresponding classical dynamics. This is demonstrated and analyzed in detail for Gaussian wave packets and a coherent state. Finally, we derive exact and approximate expressions for time-dependent transition probabilities and transition rates between quantum states. The validity of time-dependent perturbation theory is analyzed by comparison to exact dynamics.
Dynamical chaos has recently been shown to exist in the Gaussian approximation in quantum mechanics and in the self-consistent mean field approach to studying the dynamics of quantum fields. In this study, we first note that any variational approximation to the dynamics of a quantum system based on the Dirac action principle leads to a classical Hamiltonian dynamics for the variational parameters. Since this Hamiltonian is generically nonlinear and nonintegrable, the dynamics thus generated can be chaotic, in distinction to the exact quantum evolution. We then restrict our attention to a system of two biquadratically coupled quantum oscillators and study two variational schemes, the leading order large-NNN (four canonical variables) and Hartree (six canonical variables) approximations. The chaos seen in the approximate dynamics is an artifact of the approximations: this is demonstrated by the fact that its onset occurs on the same characteristic time scale as the breakdown of the approximations when compared to numerical solutions of the time-dependent Schr\"odinger equation.
In this paper we advance a new approximation scheme for the dynamics of an N particle system, which is based on the expansion of the Heisenberg equations of motion in powers of h/ in the coherent state representation. The time evolution of the physical observables is generated classically by a Hamiltonian Hs, which is derived from the original Hamiltonian by adding 2N2 virtual particles to the system. These virtual particles are related to the dispersion of the quantum mechanical wave function, and reflect the (time dependent) sensitivity of the instantaneous positions and momenta xi, pi to the variation of the initial conditions. The approximate equations of motion preserve the commutation relations between the coordinates and the momenta and the total (quantum) energy. The present approximation may be adequate for the study of the dynamics of coupled anharmonic oscillators in the quasiperiodic regime.
This paper extends the energy-based version of the stochastic linearization method, known for classical nonlinear systems, to open quantum systems with canonically commuting dynamic variables governed by quantum stochastic differential equations with non-quadratic Hamiltonians. The linearization proceeds by approximating the actual Hamiltonian of the quantum system by a quadratic function of its observables which corresponds to the Hamiltonian of a quantum harmonic oscillator. This approximation is carried out in a mean square optimal sense with respect to a Gaussian reference quantum state and leads to a self-consistent linearization procedure where the mean vector and quantum covariance matrix of the system observables evolve in time according to the effective linear dynamics. We demonstrate the proposed Hamiltonian-based Gaussian linearization for the quantum Duffing oscillator whose Hamiltonian is a quadro-quartic polynomial of the momentum and position operators. The results of the paper are applicable to the design of suboptimal controllers and filters for nonlinear quantum systems.
Nonlinear response functions, formulated as multipoint correlation functions or Volterra kernels, encode the dynamical and spectroscopic properties of physical systems and underpin a wide range of nonlinear transport and optical phenomena. However, their evaluation rapidly becomes prohibitive at high orders because of combinatorial (often factorial) scaling or severe numerical errors. Here, we establish a systematic and efficient framework to compute nonlinear response functions directly from real-time dynamics, without explicitly constructing multipoint correlators or relying on numerically unstable finite-difference methods for order-resolved extraction. Our approach is based on the Gateaux derivative with respect to the external field in function space, which yields a closed hierarchy of tangent equations of motion (TEOM). Propagating the TEOM alongside the original dynamics isolates each perturbative order with high accuracy, providing a term-by-term decomposition of physical contributions. The computational cost scales exponentially with response order in the fully general setting and reduces to polynomial complexity when all perturbation directions are identical; both regimes avoid the factorial scaling of explicit multipoint-correlator evaluations. We demonstrate the power of TEOM by computing frequency-resolved fifth-order response functions for a solid-state electron model and by obtaining nonlinear response functions up to the 49th order with controlled accuracy in a classical Duffing oscillator. We further show that our time-evolution formulation allows optical conductivities to be evaluated directly while remaining numerically stable even near zero frequency. TEOM can be incorporated seamlessly into existing real-time evolution methods, yielding a general framework for computing nonlinear response functions in quantum and classical dynamical systems.