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:1705.06203 · Physical Review A

Nonequilibrium dynamics of spin-boson models from phase-space methods

Asier Piñeiro Orioli, Arghavan Safavi-Naini, Michael L. Wall, Ana María Rey

Published Sep 5, 2017• 46 citations
arXivHTMLPDF

Abstract

An accurate description of the nonequilibrium dynamics of systems with coupled spin and bosonic degrees of freedom remains theoretically challenging, especially for large system sizes and in higher than one dimension. Phase-space methods such as the truncated Wigner approximation (TWA) have the advantage of being easily scalable and applicable to arbitrary dimensions. In this work we adapt the TWA to generic spin-boson models by making use of recently developed algorithms for discrete phase spaces [J. Schachenmayer, A. Pikovski, and A. M. Rey, Phys. Rev. X 5, 011022 (2015)]. Furthermore we go beyond the standard TWA approximation by applying a scheme based on the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy of equations to our coupled spin-boson model. This allows us, in principle, to study how systematically adding higher-order corrections improves the convergence of the method. To test various levels of approximation we study an exactly solvable spin-boson model, which is particularly relevant for trapped-ion arrays. Using TWA and its BBGKY extension we accurately reproduce the time evolution of a number of one- and two-point correlation functions in several dimensions and for an arbitrary number of bosonic modes.

DOI
10.1103/physreva.96.033607

Related papers

Simulation of quantum spin dynamics by phase space sampling of Bogoliubov-Born-Green-Kirkwood-Yvon trajectories

Lorenzo Pucci, Analabha Roy, Michael Kästner

A numerical method, suitable for the simulation of the time evolution of quantum spin models of arbitrary lattice dimension, is presented. The method combines sampling of the Wigner function with evolution equations obtained from the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy. Quantum correlations are treated through both the Wigner function sampling and the BBGKY evolution, bringing about highly accurate estimates of correlation functions. The method is particularly suitable for long-range interacting systems, and we demonstrate its power by comparing with exact results as well as other numerical methods. As an application we compute spin squeezing in a frustrated two-dimensional lattice with antiferromagnetic power-law interactions and a transverse field, making predictions pertinent to future ion trap experiments.

arXiv:1510.03768•Physical review. B./Physical review. B•May 4, 2016•51 citations

The truncated Wigner approximation on discrete phase spaces

Johannes Schachenmayer, Alexander Pikovski, Ana María Rey

A phase-space method, suitable to study the many-body spin dynamics in generic spin lattice models is introduced, discussed, and tested. This technique, which we refer to as the dTWA, is based on the discrete sampling of the Wigner function instead of the continuous sampling in traditional Truncated Wigner Approximation (TWA) methods. The power of the method is demonstrated by comparisons with analytical and numerically exact calculations, which show a significant improvement compared to the TWA, e.g. by capturing revivals of spin-observables. The dTWA is demonstrated to not only correctly reproduce dynamics of one-particle observables but also two-point correlations and spin squeezing, a property relying on entanglement. This method thus opens the possibility to study the quantum dynamics in long-range spin models, accessible to recent state-of-the-art experiments, in regimes where other numerical techniques are inapplicable.

•Aug 19, 2014•0 citations

Dynamics of correlations in two-dimensional quantum spin models with long-range interactions: a phase-space Monte-Carlo study

Johannes Schachenmayer, Alexander Pikovski, A M Rey

Interacting quantum spin models are remarkably useful for describing different types of physical, chemical, and biological systems.Significant understanding of their equilibrium properties has been achieved to date, especially for the case of spin models with short-range couplings.However, progress toward the development of a comparable understanding in long-range interacting models, in particular out-of-equilibrium, remains limited.In a recent work, we proposed a semiclassical numerical method to study spin models, the discrete truncated Wigner approximation (DTWA), and demonstrated its capability to correctly capture the dynamics of one-and two-point correlations in one-dimensional (1D) systems.Here we go one step forward and use the DTWA method to study the dynamics of correlations in two-dimensional (2D) systems with many spins and different types of long-range couplings, in regimes where other numerical methods are generally unreliable.We compute spatial and time-dependent correlations for spin-couplings that decay with distance as a power-law and determine the velocity at which correlations propagate through the system.Sharp changes in the behavior of those velocities are found as a function of the power-law decay exponent.Our predictions are relevant for a broad range of systems including solid state materials, atom-photon systems and ultracold gases of polar molecules, trapped ions, Rydberg, and magnetic atoms.We validate the DTWA predictions for small 2D systems and 1D systems, but ultimately, in the spirt of quantum simulation, experiments will be needed to confirm our predictions for large 2D systems.

arXiv:1501.06593•New Journal of Physics•Jun 10, 2015•69 citations

Two-dimensional correlation propagation dynamics with a cluster discrete phase-space method

Kazuma Nagao, Seiji Yunoki

Nonequilibrium dynamics of highly controlled quantum systems is a challenging issue in statistical physics and quantum many-body physics, relevant to recent experimental developments of analog and digital quantum simulations. In this work, we develop a discrete phase-space approach for general $\mathrm{SU}(N$) spin systems that utilizes cluster mean-field equations, which capture nontrivial quantum correlations inside each cluster, beyond the capability of the standard discrete truncated Wigner approximation for individual classical spins. Our formalism, based on a cluster phase-point operator, enables efficient numerical samplings of cluster phase-space variables, where the total number of noise variables for a direct product state is independent of the specific way in which the entire system is divided into multiple equally sized finite clusters. We numerically demonstrate that the cluster discrete truncated Wigner approximation (C-dTWA) method can reproduce key results in a recent experiment on correlation propagation dynamics in a two-dimensional Bose-Hubbard system. We further compare the results of C-dTWA for clusters of $2\ifmmode\times\else\texttimes\fi{}2$ sites with those from a two-dimensional tensor network method and discuss that both approaches agree very well in the short-time region, where the energy conservation is well maintained in the tensor network simulations. Since we formulate the C-dTWA method in a general form, it has the potential for application to various dynamical problems in isolated and open quantum systems, even in higher dimensions.

arXiv:2404.18594•Physical Review A•Feb 12, 2025•4 citations

Path integral approach to the truncated Wigner approximation of driven-dissipative spins

Viktoria Noel, Michael Fleischhauer, Igor Lesanovsky

Phase-space approaches such as the truncated Wigner approximation (TWA) provide an efficient semiclassical framework for performing approximate simulations of the dynamics of open quantum many-body systems outside the reach of exact numerical methods but beyond the mean-field level. For bosonic systems, TWA is known to be equivalent to a Keldysh path-integral formulation truncated at second order in the so-called quantum fluctuations. This semiclassical approach provides an alternative transparent route towards approximate stochastic equations of motion, which can be efficiently solved. Here we establish the corresponding path-integral formulation for interacting open spin-$1/2$ systems using the continuous $\mathrm{SU}(2)$ phase space. We show, in particular, that a consistent treatment of dissipation requires correctly mapping operator products onto the curved spin phase space, leading to stochastic equations that coincide with those obtained from the continuous TWA formulation and thus reproduce the exact dynamics of a single dissipative spin. Our results provide a unified field-theoretic foundation for the TWA to dissipative spin dynamics and offer a systematic starting point for extensions beyond the semiclassical approximation.

quant-pharXiv:2607.26029•Jul 29, 2026•0 citations

A generalized phase space approach for solving quantum spin dynamics

Bihui Zhu, Ana Maria Rey, Johannes Schachenmayer

Abstract Numerical techniques to efficiently model out-of-equilibrium dynamics in interacting quantum many-body systems are key for advancing our capability to harness and understand complex quantum matter. Here we propose a new numerical approach which we refer to as generalized discrete truncated Wigner approximation (GDTWA). It is based on a discrete semi-classical phase space sampling and allows to investigate quantum dynamics in lattice spin systems with arbitrary S ≥ 1/2. We show that the GDTWA can accurately simulate dynamics of large ensembles in arbitrary dimensions. We apply it for S > 1/2 spin-models with dipolar long-range interactions, a scenario arising in recent experiments with magnetic atoms. We show that the method can capture beyond mean-field effects, not only at short times, but it also can correctly reproduce long time quantum-thermalization dynamics. We benchmark the method with exact diagonalization in small systems, with perturbation theory for short times, and with analytical predictions made for models which feature quantum-thermalization at long times. We apply our method to study dynamics in large S > 1/2 spin-models and compute experimentally accessible observables such as Zeeman level populations, contrast of spin coherence, spin squeezing, and entanglement quantified by single-spin Renyi entropies. We reveal that large S systems can feature larger entanglement than corresponding S = 1/2 systems. Our analyses demonstrate that the GDTWA can be a powerful tool for modeling complex spin dynamics in regimes where other state-of-the art numerical methods fail.

arXiv:1905.08782•New Journal of Physics•Jul 25, 2019•59 citations
An accurate description of the nonequilibrium dynamics of systems with coupled spin and bosonic degrees of freedom remains theoretically challenging, especially for large system sizes and in higher than one dimension. Phase-space methods such as the truncated Wigner approximation (TWA) have the advantage of being easily scalable and applicable to arbitrary dimensions. In this work we adapt the TWA to generic spin-boson models by making use of recently developed algorithms for discrete phase spaces [J. Schachenmayer, A. Pikovski, and A. M. Rey, Phys. Rev. X 5, 011022 (2015)]. Furthermore we go beyond the standard TWA approximation by applying a scheme based on the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy of equations to our coupled spin-boson model. This allows us, in principle, to study how systematically adding higher-order corrections improves the convergence of the method. To test various levels of approximation we study an exactly solvable spin-boson model, which is particularly relevant for trapped-ion arrays. Using TWA and its BBGKY extension we accurately reproduce the time evolution of a number of one- and two-point correlation functions in several dimensions and for an arbitrary number of bosonic modes.
A numerical method, suitable for the simulation of the time evolution of quantum spin models of arbitrary lattice dimension, is presented. The method combines sampling of the Wigner function with evolution equations obtained from the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy. Quantum correlations are treated through both the Wigner function sampling and the BBGKY evolution, bringing about highly accurate estimates of correlation functions. The method is particularly suitable for long-range interacting systems, and we demonstrate its power by comparing with exact results as well as other numerical methods. As an application we compute spin squeezing in a frustrated two-dimensional lattice with antiferromagnetic power-law interactions and a transverse field, making predictions pertinent to future ion trap experiments.
A phase-space method, suitable to study the many-body spin dynamics in generic spin lattice models is introduced, discussed, and tested. This technique, which we refer to as the dTWA, is based on the discrete sampling of the Wigner function instead of the continuous sampling in traditional Truncated Wigner Approximation (TWA) methods. The power of the method is demonstrated by comparisons with analytical and numerically exact calculations, which show a significant improvement compared to the TWA, e.g. by capturing revivals of spin-observables. The dTWA is demonstrated to not only correctly reproduce dynamics of one-particle observables but also two-point correlations and spin squeezing, a property relying on entanglement. This method thus opens the possibility to study the quantum dynamics in long-range spin models, accessible to recent state-of-the-art experiments, in regimes where other numerical techniques are inapplicable.
Interacting quantum spin models are remarkably useful for describing different types of physical, chemical, and biological systems.Significant understanding of their equilibrium properties has been achieved to date, especially for the case of spin models with short-range couplings.However, progress toward the development of a comparable understanding in long-range interacting models, in particular out-of-equilibrium, remains limited.In a recent work, we proposed a semiclassical numerical method to study spin models, the discrete truncated Wigner approximation (DTWA), and demonstrated its capability to correctly capture the dynamics of one-and two-point correlations in one-dimensional (1D) systems.Here we go one step forward and use the DTWA method to study the dynamics of correlations in two-dimensional (2D) systems with many spins and different types of long-range couplings, in regimes where other numerical methods are generally unreliable.We compute spatial and time-dependent correlations for spin-couplings that decay with distance as a power-law and determine the velocity at which correlations propagate through the system.Sharp changes in the behavior of those velocities are found as a function of the power-law decay exponent.Our predictions are relevant for a broad range of systems including solid state materials, atom-photon systems and ultracold gases of polar molecules, trapped ions, Rydberg, and magnetic atoms.We validate the DTWA predictions for small 2D systems and 1D systems, but ultimately, in the spirt of quantum simulation, experiments will be needed to confirm our predictions for large 2D systems.
Nonequilibrium dynamics of highly controlled quantum systems is a challenging issue in statistical physics and quantum many-body physics, relevant to recent experimental developments of analog and digital quantum simulations. In this work, we develop a discrete phase-space approach for general SU(N\mathrm{SU}(NSU(N) spin systems that utilizes cluster mean-field equations, which capture nontrivial quantum correlations inside each cluster, beyond the capability of the standard discrete truncated Wigner approximation for individual classical spins. Our formalism, based on a cluster phase-point operator, enables efficient numerical samplings of cluster phase-space variables, where the total number of noise variables for a direct product state is independent of the specific way in which the entire system is divided into multiple equally sized finite clusters. We numerically demonstrate that the cluster discrete truncated Wigner approximation (C-dTWA) method can reproduce key results in a recent experiment on correlation propagation dynamics in a two-dimensional Bose-Hubbard system. We further compare the results of C-dTWA for clusters of $2\ifmmode\times\else\texttimes\fi{}2$ sites with those from a two-dimensional tensor network method and discuss that both approaches agree very well in the short-time region, where the energy conservation is well maintained in the tensor network simulations. Since we formulate the C-dTWA method in a general form, it has the potential for application to various dynamical problems in isolated and open quantum systems, even in higher dimensions.
Phase-space approaches such as the truncated Wigner approximation (TWA) provide an efficient semiclassical framework for performing approximate simulations of the dynamics of open quantum many-body systems outside the reach of exact numerical methods but beyond the mean-field level. For bosonic systems, TWA is known to be equivalent to a Keldysh path-integral formulation truncated at second order in the so-called quantum fluctuations. This semiclassical approach provides an alternative transparent route towards approximate stochastic equations of motion, which can be efficiently solved. Here we establish the corresponding path-integral formulation for interacting open spin-1/21/21/2 systems using the continuous SU(2)\mathrm{SU}(2)SU(2) phase space. We show, in particular, that a consistent treatment of dissipation requires correctly mapping operator products onto the curved spin phase space, leading to stochastic equations that coincide with those obtained from the continuous TWA formulation and thus reproduce the exact dynamics of a single dissipative spin. Our results provide a unified field-theoretic foundation for the TWA to dissipative spin dynamics and offer a systematic starting point for extensions beyond the semiclassical approximation.
Abstract Numerical techniques to efficiently model out-of-equilibrium dynamics in interacting quantum many-body systems are key for advancing our capability to harness and understand complex quantum matter. Here we propose a new numerical approach which we refer to as generalized discrete truncated Wigner approximation (GDTWA). It is based on a discrete semi-classical phase space sampling and allows to investigate quantum dynamics in lattice spin systems with arbitrary S ≥ 1/2. We show that the GDTWA can accurately simulate dynamics of large ensembles in arbitrary dimensions. We apply it for S > 1/2 spin-models with dipolar long-range interactions, a scenario arising in recent experiments with magnetic atoms. We show that the method can capture beyond mean-field effects, not only at short times, but it also can correctly reproduce long time quantum-thermalization dynamics. We benchmark the method with exact diagonalization in small systems, with perturbation theory for short times, and with analytical predictions made for models which feature quantum-thermalization at long times. We apply our method to study dynamics in large S > 1/2 spin-models and compute experimentally accessible observables such as Zeeman level populations, contrast of spin coherence, spin squeezing, and entanglement quantified by single-spin Renyi entropies. We reveal that large S systems can feature larger entanglement than corresponding S = 1/2 systems. Our analyses demonstrate that the GDTWA can be a powerful tool for modeling complex spin dynamics in regimes where other state-of-the art numerical methods fail.