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:2008.13481 · Physical Review Research

Performance evaluation of the discrete truncated Wigner approximation for quench dynamics of quantum spin systems with long-range interactions

Masaya Kunimi, Kazuma Nagao, Shimpei Goto, Ippei Danshita

Published Jan 19, 2021• 16 citations
arXivHTMLPDF

Abstract

The authors study a validity timescale of the discrete truncated Wigner approximation for several quantum-spin systems with long-range interactions of step-function form by using the Bogoliubov-Born-Green-Kirkwood-Yvon equation.
DOI
10.1103/physrevresearch.3.013060

Related papers

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

Asier Piñeiro Orioli, Arghavan Safavi-Naini, Michael L. Wall +1 more

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.

arXiv:1705.06203•Physical Review A•Sep 5, 2017•46 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

Quenches near criticality of the quantum Ising chain—power and limitations of the discrete truncated Wigner approximation

Stefanie Czischek, Martin Gärttner, Markus K. Oberthaler +2 more

Abstract The semi-classical discrete truncated Wigner approximation has recently been proposed as a simulation method for spin-1/2 systems. While it appears to provide a powerful approach which shows promising results in higher dimensions and for systems with long-range interactions, its performance is still not well understood in general. Here we perform a systematic benchmarking on the one-dimensional transverse-field Ising model and point to limitations of the approximation arising after sudden quenches into the quantum critical regime. Our procedure allows to identify the limitations of the semi-classical simulations and with that to determine the regimes and questions where quantum simulators can provide information which is inaccessible to semi-classics.

arXiv:1805.05221•Quantum Science and Technology•Sep 25, 2018•18 citations

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

Discrete truncated Wigner approach to dynamical phase transitions in Ising models after a quantum quench

Reyhaneh Khasseh, Angelo Russomanno, Markus Schmitt +2 more

By means of the discrete truncated Wigner approximation, we study dynamical phase transitions arising in the steady state of transverse-field Ising models after a quantum quench. Starting from a fully polarized ferromagnetic initial condition, these transitions separate a phase with nonvanishing magnetization along the ordering direction from a disordered symmetric phase upon increasing the transverse field. We consider two paradigmatic cases, a one-dimensional long-range model with power-law interactions $\ensuremath{\propto}1/{r}^{\ensuremath{\alpha}}$ decaying algebraically as a function of distance $r$ and a two-dimensional system with short-range nearest-neighbor interactions. In the former case, we identify dynamical phase transitions for $\ensuremath{\alpha}\ensuremath{\lesssim}2$ and we extract the critical exponents from a data collapse of the steady-state magnetization for up to 1200 lattice sites. We find identical exponents for $\ensuremath{\alpha}\ensuremath{\lesssim}0.5$, suggesting that the dynamical transitions in this regime fall into the same universality class as the nonergodic mean-field limit. The two-dimensional Ising model is believed to be thermalizing, which we also confirm using exact diagonalization for small system sizes. Thus, the dynamical transition is expected to correspond to the thermal phase transition, which is consistent with our data upon comparing to equilibrium quantum Monte Carlo simulations. We further test the accuracy of the discrete truncated Wigner approximation by comparing against numerically exact methods such as exact diagonalization, tensor network, as well as artificial neural network states and we find good quantitative agreement on the accessible time scales. Finally, our work provides an additional contribution to the understanding of the range and the limitations of qualitative and quantitative applicability of the discrete truncated Wigner approximation.

arXiv:2004.09812•Physical review. B./Physical review. B•Jul 6, 2020•36 citations

Truncated Wigner dynamics of biclique quantum spin glasses

Dries Sels

Quantum spin glasses are often considered testbeds for studying quantum optimization algorithms and as such have been the subject of various quantum advantage claims. Here we investigate the near adiabatic dynamics of biclique quantum spin glasses within the (discrete) truncated Wigner approximation (TWA). Benchmarks on small systems show that TWA recovers sample-to-sample fluctuations of the Edwards-Anderson order parameter, over a wide range of annealing times, with increasing fidelity when the system size increases. We extract critical exponents from the Binder cumulant in line with theoretical expectations, reproducing recent quantum experiments. The computational cost of the method is minimal and it can easily be applied to tens of thousands of qubits.

cond-mat.dis-nnCross-listedarXiv:2606.20187•Jun 19, 2026•0 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.
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.
Abstract The semi-classical discrete truncated Wigner approximation has recently been proposed as a simulation method for spin-1/2 systems. While it appears to provide a powerful approach which shows promising results in higher dimensions and for systems with long-range interactions, its performance is still not well understood in general. Here we perform a systematic benchmarking on the one-dimensional transverse-field Ising model and point to limitations of the approximation arising after sudden quenches into the quantum critical regime. Our procedure allows to identify the limitations of the semi-classical simulations and with that to determine the regimes and questions where quantum simulators can provide information which is inaccessible to semi-classics.
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.
By means of the discrete truncated Wigner approximation, we study dynamical phase transitions arising in the steady state of transverse-field Ising models after a quantum quench. Starting from a fully polarized ferromagnetic initial condition, these transitions separate a phase with nonvanishing magnetization along the ordering direction from a disordered symmetric phase upon increasing the transverse field. We consider two paradigmatic cases, a one-dimensional long-range model with power-law interactions $\ensuremath{\propto}1/{r}^{\ensuremath{\alpha}}$ decaying algebraically as a function of distance rrr and a two-dimensional system with short-range nearest-neighbor interactions. In the former case, we identify dynamical phase transitions for $\ensuremath{\alpha}\ensuremath{\lesssim}2$ and we extract the critical exponents from a data collapse of the steady-state magnetization for up to 1200 lattice sites. We find identical exponents for $\ensuremath{\alpha}\ensuremath{\lesssim}0.5$, suggesting that the dynamical transitions in this regime fall into the same universality class as the nonergodic mean-field limit. The two-dimensional Ising model is believed to be thermalizing, which we also confirm using exact diagonalization for small system sizes. Thus, the dynamical transition is expected to correspond to the thermal phase transition, which is consistent with our data upon comparing to equilibrium quantum Monte Carlo simulations. We further test the accuracy of the discrete truncated Wigner approximation by comparing against numerically exact methods such as exact diagonalization, tensor network, as well as artificial neural network states and we find good quantitative agreement on the accessible time scales. Finally, our work provides an additional contribution to the understanding of the range and the limitations of qualitative and quantitative applicability of the discrete truncated Wigner approximation.
Quantum spin glasses are often considered testbeds for studying quantum optimization algorithms and as such have been the subject of various quantum advantage claims. Here we investigate the near adiabatic dynamics of biclique quantum spin glasses within the (discrete) truncated Wigner approximation (TWA). Benchmarks on small systems show that TWA recovers sample-to-sample fluctuations of the Edwards-Anderson order parameter, over a wide range of annealing times, with increasing fidelity when the system size increases. We extract critical exponents from the Binder cumulant in line with theoretical expectations, reproducing recent quantum experiments. The computational cost of the method is minimal and it can easily be applied to tens of thousands of qubits.