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

Representation of quantum dynamics of interacting systems through classical trajectories

Anatoli Polkovnikov

Published May 20, 2009• 1 citation
OpenAlex

Abstract

We discuss a phase space representation of quantum dynamics of systems with many degrees of freedom. This representation is based on a perturbative expansion in quantum fluctuations around one of the classical limits. We explicitly analyze expansions around three such limits: (i) corpuscular or Newtonian limit in the coordinate-momentum representation, (ii) wave or Gross-Pitaevskii limit for interacting bosons in the coherent state representation, and (iii) Bloch limit for the spin systems. We discuss both the semiclassical (truncated Wigner) approximation and further quantum corrections appearing in the form of either stochastic quantum jumps along the classical trajectories or the nonlinear response to such jumps. We also discuss how quantum jumps naturally emerge in the analysis of non-equal time correlation functions. This representation of quantum dynamics is closely related to the phase space methods based on the Wigner-Weyl quantization and to the Keldysh technique. We show how such concepts as the Wigner function, Weyl symbol, Moyal product, Bopp operators, and others automatically emerge from the Feynmann's path integral representation of the evolution in the Heisenberg representation. We illustrate the applicability of this expansion with various examples mostly in the context of cold atom systems including sine-Gordon model, one- and two-dimensional Bose Hubbard model, Dicke model and others.

Related papers

Quantum corrections to the dynamics of interacting bosons: Beyond the truncated Wigner approximation

Anatoli Polkovnikov

We develop a consistent perturbation theory in quantum fluctuations around the classical evolution of a system of interacting bosons. The zero-order approximation gives the classical Gross-Pitaevskii equations. In the next order we recover the truncated Wigner approximation, where the evolution is still classical but the initial conditions are distributed according to the Wigner transform of the initial density matrix. Further corrections can be characterized as quantum scattering events, which appear in the form of a nonlinear response of the observable to an infinitesimal displacement of the field along its classical evolution. At the end of the paper we give a few numerical examples to test the formalism.

arXiv:cond-mat/0303628•Physical Review A•Nov 4, 2003•132 citations

Moyal quantum mechanics: the semiclassical Heisenberg dynamics

T. A. Osborn, Molzahn, F H

The Moyal description of quantum mechanics, based on the Wigner-Weyl isomorphism between operators and symbols, provides a comprehensive phase space representation of dynamics. The Weyl symbol image of the Heisenberg picture evolution operator is regular in ħ and so presents a preferred foundation for semiclassical analysis. Its semiclassical expansion "coefficients," acting on symbols that represent observables, are simple, globally defined (phase space) differential operators constructed in terms of the classical flow. The first of two presented methods introduces a cluster-graph expansion for the symbol of an exponentiated operator, which extends Groenewold′s formula for the Weyl product of two symbols and has ħ as its natural small parameter. This Poisson bracket based cluster expansion determines the Jacobi equations for the semiclassical expansion of "quantum trajectories." Their Green function solutions construct the regular ħ ⇓ 0 asymptotic series for the Heisenberg-Weyl evolution map. The second method directly substitutes such a series into the Moyal equation of motion and determines the ħ coefficients recursively. In contrast to the WKB approximation for propagators, the Heisenberg-Weyl description of evolution involves no essential singularity in ħ, no Hamilton-Jacobi equation to solve for the action, and no multiple trajectories, caustics, or Maslov indices.

•Jan 1, 1994•46 citations

Phase space path integrals in Monte Carlo quantum dynamics

Stavros Caratzoulas, Philip Pechukas

In this article we present two phase space path integrals in terms of coherent states. The first one is derived in a standard fashion but using a nonstandard resolution of the identity in terms of coherent states with different width parameters. The second path integral emerges from a novel phase space representation in terms of coherent states distributed on n-dimensional manifolds embedded in the 2n-dimensional phase space of an n-degree-of-freedom system. These states are shown to form locally complete basis sets since we show that fairly smooth and localized functions can be expanded in terms of them in a unique way. In this representation the time evolution operator can be cast in the form of a phase space path integral. Both path integrals can be evaluated by straightforward implementation of Monte Carlo methods. In both cases the probability amplitude between two phase points turns out to be proportional to the average of the phase, eih∫(p dq−H dt), over a Markov process of phase space paths consisting of classical trajectories interrupted at intervals by Gaussian ‘‘quantum jumps.’’ The numerical evaluation of these expressions through importance sampling is demonstrated.

The Journal of Chemical Physics•Apr 22, 1996•10 citations

SU(3) Semiclassical Representation of Quantum Dynamics of Interacting Spins

Shainen M. Davidson, Anatoli Polkovnikov

We present a formalism for simulating quantum dynamics of lattice spin-1 systems by first introducing local hidden variables and then doing semiclassical (truncated Wigner) approximation in the extended phase space. In this way, we exactly take into account the local on-site Hamiltonian and approximately treat spin-spin interactions. In particular, we represent each spin with eight classical SU(3) variables. Three of them represent the usual spin components and five others are hidden variables representing local spin-spin correlations. We compare our formalism with exact quantum dynamics of fully connected spin systems and find very good agreement. As an application, we discuss quench dynamics of a Bose-Hubbard model near the superfluid-insulator transition for a 3D lattice system consisting of 1000 sites.

arXiv:1406.1192•Physical Review Letters•Jan 29, 2015•52 citations

Coherent state path integrals in the Weyl representation

Luana Carolina Santos, Marcus A. M. de Aguiar

We construct a representation of the coherent state path integral using the Weyl symbol of the Hamiltonian operator. This representation is very different from the usual path integral forms suggested by Klauder and Skagerstam (1985 Coherent States: Applications in Physics and Mathematical Physics (Singapore: World Scientific)), which involve the normal or the antinormal ordering of the Hamiltonian. These different representations, although equivalent quantum mechanically, lead to different semiclassical limits. We show that the semiclassical limit of the coherent state propagator in the Weyl representation involves classical trajectories that are independent of the width of coherent states. This propagator is also free from the phase corrections found in Baranger et al ( 2001 J. Phys. A: Math. Gen. 34 7227 ) for the two Klauder forms and provides an explicit connection between the Wigner and the Husimi representations of the evolution operator.

arXiv:quant-ph/0607136•Journal of Physics A Mathematical and General•Oct 10, 2006•16 citations

Phase-space algorithm for simulating quantum nonlinear response functions of bosons using stochastic classical trajectories

Benoit Palmieri, Yuki Nagata, Shaul Mukamel

Using the positive P-representation of the density matrix, we develop an algorithm for calculating the quantum many-body nonlinear response functions of a system of bosons driven impulsively by external fields. The formalism maps the quantum time evolution of N boson degrees of freedom into a stochastic dynamics of 4N classical degrees of freedom. The first- and the third-order response functions are calculated by propagating the parameters of the P-representation using a set of coupled Langevin equations with multiplicative noise. These parameters serve as classical variables. Two classical ways for computing the response functions are presented. In the nonequilibrium method, an observable is calculated for weak impulsive pulses, and the response functions are obtained by taking its derivatives with respect to the pulse amplitudes. In the alternative, equilibrium simulation, the response functions are expressed in terms of time-correlation functions involving the P-representation parameters and stability matrices representing the perturbation of the trajectories. The stability matrices can be propagated simultaneously with the Langevin equations for the parameters. The formalism is generalized for a many-body boson system coupled to a harmonic bath.

Physical Review E•Oct 12, 2010•2 citations
We discuss a phase space representation of quantum dynamics of systems with many degrees of freedom. This representation is based on a perturbative expansion in quantum fluctuations around one of the classical limits. We explicitly analyze expansions around three such limits: (i) corpuscular or Newtonian limit in the coordinate-momentum representation, (ii) wave or Gross-Pitaevskii limit for interacting bosons in the coherent state representation, and (iii) Bloch limit for the spin systems. We discuss both the semiclassical (truncated Wigner) approximation and further quantum corrections appearing in the form of either stochastic quantum jumps along the classical trajectories or the nonlinear response to such jumps. We also discuss how quantum jumps naturally emerge in the analysis of non-equal time correlation functions. This representation of quantum dynamics is closely related to the phase space methods based on the Wigner-Weyl quantization and to the Keldysh technique. We show how such concepts as the Wigner function, Weyl symbol, Moyal product, Bopp operators, and others automatically emerge from the Feynmann's path integral representation of the evolution in the Heisenberg representation. We illustrate the applicability of this expansion with various examples mostly in the context of cold atom systems including sine-Gordon model, one- and two-dimensional Bose Hubbard model, Dicke model and others.
We develop a consistent perturbation theory in quantum fluctuations around the classical evolution of a system of interacting bosons. The zero-order approximation gives the classical Gross-Pitaevskii equations. In the next order we recover the truncated Wigner approximation, where the evolution is still classical but the initial conditions are distributed according to the Wigner transform of the initial density matrix. Further corrections can be characterized as quantum scattering events, which appear in the form of a nonlinear response of the observable to an infinitesimal displacement of the field along its classical evolution. At the end of the paper we give a few numerical examples to test the formalism.
The Moyal description of quantum mechanics, based on the Wigner-Weyl isomorphism between operators and symbols, provides a comprehensive phase space representation of dynamics. The Weyl symbol image of the Heisenberg picture evolution operator is regular in ħ and so presents a preferred foundation for semiclassical analysis. Its semiclassical expansion "coefficients," acting on symbols that represent observables, are simple, globally defined (phase space) differential operators constructed in terms of the classical flow. The first of two presented methods introduces a cluster-graph expansion for the symbol of an exponentiated operator, which extends Groenewold′s formula for the Weyl product of two symbols and has ħ as its natural small parameter. This Poisson bracket based cluster expansion determines the Jacobi equations for the semiclassical expansion of "quantum trajectories." Their Green function solutions construct the regular ħ ⇓ 0 asymptotic series for the Heisenberg-Weyl evolution map. The second method directly substitutes such a series into the Moyal equation of motion and determines the ħ coefficients recursively. In contrast to the WKB approximation for propagators, the Heisenberg-Weyl description of evolution involves no essential singularity in ħ, no Hamilton-Jacobi equation to solve for the action, and no multiple trajectories, caustics, or Maslov indices.
In this article we present two phase space path integrals in terms of coherent states. The first one is derived in a standard fashion but using a nonstandard resolution of the identity in terms of coherent states with different width parameters. The second path integral emerges from a novel phase space representation in terms of coherent states distributed on n-dimensional manifolds embedded in the 2n-dimensional phase space of an n-degree-of-freedom system. These states are shown to form locally complete basis sets since we show that fairly smooth and localized functions can be expanded in terms of them in a unique way. In this representation the time evolution operator can be cast in the form of a phase space path integral. Both path integrals can be evaluated by straightforward implementation of Monte Carlo methods. In both cases the probability amplitude between two phase points turns out to be proportional to the average of the phase, eih∫(p dq−H dt), over a Markov process of phase space paths consisting of classical trajectories interrupted at intervals by Gaussian ‘‘quantum jumps.’’ The numerical evaluation of these expressions through importance sampling is demonstrated.
We present a formalism for simulating quantum dynamics of lattice spin-1 systems by first introducing local hidden variables and then doing semiclassical (truncated Wigner) approximation in the extended phase space. In this way, we exactly take into account the local on-site Hamiltonian and approximately treat spin-spin interactions. In particular, we represent each spin with eight classical SU(3) variables. Three of them represent the usual spin components and five others are hidden variables representing local spin-spin correlations. We compare our formalism with exact quantum dynamics of fully connected spin systems and find very good agreement. As an application, we discuss quench dynamics of a Bose-Hubbard model near the superfluid-insulator transition for a 3D lattice system consisting of 1000 sites.
We construct a representation of the coherent state path integral using the Weyl symbol of the Hamiltonian operator. This representation is very different from the usual path integral forms suggested by Klauder and Skagerstam (1985 Coherent States: Applications in Physics and Mathematical Physics (Singapore: World Scientific)), which involve the normal or the antinormal ordering of the Hamiltonian. These different representations, although equivalent quantum mechanically, lead to different semiclassical limits. We show that the semiclassical limit of the coherent state propagator in the Weyl representation involves classical trajectories that are independent of the width of coherent states. This propagator is also free from the phase corrections found in Baranger et al ( 2001 J. Phys. A: Math. Gen. 34 7227 ) for the two Klauder forms and provides an explicit connection between the Wigner and the Husimi representations of the evolution operator.
Using the positive P-representation of the density matrix, we develop an algorithm for calculating the quantum many-body nonlinear response functions of a system of bosons driven impulsively by external fields. The formalism maps the quantum time evolution of N boson degrees of freedom into a stochastic dynamics of 4N classical degrees of freedom. The first- and the third-order response functions are calculated by propagating the parameters of the P-representation using a set of coupled Langevin equations with multiplicative noise. These parameters serve as classical variables. Two classical ways for computing the response functions are presented. In the nonequilibrium method, an observable is calculated for weak impulsive pulses, and the response functions are obtained by taking its derivatives with respect to the pulse amplitudes. In the alternative, equilibrium simulation, the response functions are expressed in terms of time-correlation functions involving the P-representation parameters and stability matrices representing the perturbation of the trajectories. The stability matrices can be propagated simultaneously with the Langevin equations for the parameters. The formalism is generalized for a many-body boson system coupled to a harmonic bath.