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

Discrete-phase-space method for driven-dissipative dynamics of strongly interacting bosons in optical lattices

Kazuma Nagao, Ippei Danshita, Seiji Yunoki

Published Dec 11, 2024• 5 citations
arXivHTMLPDF

Abstract

We develop a discrete truncated Wigner method to analyze the real-time evolution of dissipative $\mathrm{SU}(\mathcal{N})$ spin systems coupled with a Markovian environment. This semiclassical approach is not only numerically efficient but also particularly capable of accurately capturing local loss processes due to its local linearity in the dynamical equations. We apply the method to a state-of-the-art experiment involving an analog quantum simulator of a three-dimensional dissipative Bose-Hubbard model in a strongly interacting regime. Our numerical results show good agreement with experimental data, specifically capturing the continuous quantum Zeno effect in the dynamics subjected to a gradual change of the ratio between the hopping amplitude and the onsite interaction across the superfluid Mott insulator crossover. Furthermore, we present comparative analyses with the continuous truncated Wigner method, derived as an effective Fokker-Planck equation for $\mathrm{SU}(\mathcal{N})$ classical spin variables, showing that the discrete method outperforms the continuous one in simulating the long-time dynamics of SU(2) and SU(3) spin models. The discrete phase-space framework offers a versatile and powerful tool for exploring a wide range of open quantum many-body systems in dimensions higher than one dimension, where numerically exact methods are impractical due to the exponential growth of the Hilbert-space dimension.

DOI
10.1103/physreva.110.063310

Related papers

Hybrid discrete-continuous truncated Wigner approximation for driven, dissipative spin systems

Christopher D. Mink, David Petrosyan, Michael Fleischhauer

We present a systematic approach for the semiclassical treatment of many-body dynamics of interacting, open spin systems. Our approach overcomes some of the shortcomings of the recently developed discrete truncated Wigner approximation (DTWA) based on Monte Carlo sampling in a discrete phase space that improves the classical treatment by accounting for lowest-order quantum fluctuations. We provide a rigorous derivation of the DTWA by embedding it in a continuous phase space, thereby introducing a hybrid discrete-continuous truncated Wigner approximation . We derive a set of operator-differential mappings that yield an exact equation of motion (EOM) for the continuous SU(2) Wigner function of spins. The standard DTWA is then recovered by a systematic neglection of specific terms in this exact EOM. The hybrid approach permits us to determine the validity conditions and to gain a detailed understanding of the quality of the approximation, paving the way for systematic improvements. Furthermore, we show that the continuous embedding allows for a straightforward extension of the method to open spin systems subject to dephasing, losses, and incoherent drive, while preserving the key advantages of the discrete approach, such as a positive definite Wigner distribution of typical initial states. We derive exact stochastic differential equations for processes which cannot be described by the standard DTWA due to the presence of nonclassical noise. We illustrate our approach by applying it to the dissipative dynamics of Rydberg excitation of one-dimensional arrays of laser-driven atoms and compare it to exact results for small systems.

arXiv:2203.17120•Physical Review Research•Nov 28, 2022•28 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

SU(3) truncated Wigner approximation for strongly interacting Bose gases

Kazuma Nagao, Yosuke Takasu, Yoshiro Takahashi +1 more

We develop and utilize the SU(3) truncated Wigner approximation (TWA) in order to analyze far-from-equilibrium quantum dynamics of strongly interacting Bose gases in an optical lattice. Specifically, we explicitly represent the corresponding Bose-Hubbard model at an arbitrary filling factor with restricted local Hilbert spaces in terms of SU(3) matrices. Moreover, we introduce a discrete Wigner sampling technique for the SU(3) TWA and examine its performance as well as that of the SU(3) TWA with the Gaussian approximation for the continuous Wigner function. We directly compare outputs of these two approaches with exact computations regarding dynamics of the Bose-Hubbard model at unit filling with a small size and that of a fully connected spin-1 model with a large size. We show that both approaches can quantitatively capture quantum dynamics on a timescale of $\ensuremath{\hbar}/(Jz)$, where $J$ and $z$ denote the hopping energy and the coordination number. We apply the two kinds of SU(3) TWA to dynamical spreading of a two-point correlation function of the Bose-Hubbard model on a square lattice with a large system size, which has been measured in recent experiments. Noticeable deviations between the theories and experiments indicate that proper inclusion of effects of the spatial inhomogeneity, which is not straightforward in our formulation of the SU(3) TWA, may be necessary.

arXiv:2008.09900•Physical Review Research•Nov 1, 2021•9 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

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

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
We develop a discrete truncated Wigner method to analyze the real-time evolution of dissipative SU(N)\mathrm{SU}(\mathcal{N})SU(N) spin systems coupled with a Markovian environment. This semiclassical approach is not only numerically efficient but also particularly capable of accurately capturing local loss processes due to its local linearity in the dynamical equations. We apply the method to a state-of-the-art experiment involving an analog quantum simulator of a three-dimensional dissipative Bose-Hubbard model in a strongly interacting regime. Our numerical results show good agreement with experimental data, specifically capturing the continuous quantum Zeno effect in the dynamics subjected to a gradual change of the ratio between the hopping amplitude and the onsite interaction across the superfluid Mott insulator crossover. Furthermore, we present comparative analyses with the continuous truncated Wigner method, derived as an effective Fokker-Planck equation for SU(N)\mathrm{SU}(\mathcal{N})SU(N) classical spin variables, showing that the discrete method outperforms the continuous one in simulating the long-time dynamics of SU(2) and SU(3) spin models. The discrete phase-space framework offers a versatile and powerful tool for exploring a wide range of open quantum many-body systems in dimensions higher than one dimension, where numerically exact methods are impractical due to the exponential growth of the Hilbert-space dimension.
We present a systematic approach for the semiclassical treatment of many-body dynamics of interacting, open spin systems. Our approach overcomes some of the shortcomings of the recently developed discrete truncated Wigner approximation (DTWA) based on Monte Carlo sampling in a discrete phase space that improves the classical treatment by accounting for lowest-order quantum fluctuations. We provide a rigorous derivation of the DTWA by embedding it in a continuous phase space, thereby introducing a hybrid discrete-continuous truncated Wigner approximation . We derive a set of operator-differential mappings that yield an exact equation of motion (EOM) for the continuous SU(2) Wigner function of spins. The standard DTWA is then recovered by a systematic neglection of specific terms in this exact EOM. The hybrid approach permits us to determine the validity conditions and to gain a detailed understanding of the quality of the approximation, paving the way for systematic improvements. Furthermore, we show that the continuous embedding allows for a straightforward extension of the method to open spin systems subject to dephasing, losses, and incoherent drive, while preserving the key advantages of the discrete approach, such as a positive definite Wigner distribution of typical initial states. We derive exact stochastic differential equations for processes which cannot be described by the standard DTWA due to the presence of nonclassical noise. We illustrate our approach by applying it to the dissipative dynamics of Rydberg excitation of one-dimensional arrays of laser-driven atoms and compare it to exact results for small systems.
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.
We develop and utilize the SU(3) truncated Wigner approximation (TWA) in order to analyze far-from-equilibrium quantum dynamics of strongly interacting Bose gases in an optical lattice. Specifically, we explicitly represent the corresponding Bose-Hubbard model at an arbitrary filling factor with restricted local Hilbert spaces in terms of SU(3) matrices. Moreover, we introduce a discrete Wigner sampling technique for the SU(3) TWA and examine its performance as well as that of the SU(3) TWA with the Gaussian approximation for the continuous Wigner function. We directly compare outputs of these two approaches with exact computations regarding dynamics of the Bose-Hubbard model at unit filling with a small size and that of a fully connected spin-1 model with a large size. We show that both approaches can quantitatively capture quantum dynamics on a timescale of $\ensuremath{\hbar}/(Jz)$, where JJJ and zzz denote the hopping energy and the coordination number. We apply the two kinds of SU(3) TWA to dynamical spreading of a two-point correlation function of the Bose-Hubbard model on a square lattice with a large system size, which has been measured in recent experiments. Noticeable deviations between the theories and experiments indicate that proper inclusion of effects of the spatial inhomogeneity, which is not straightforward in our formulation of the SU(3) TWA, may be necessary.
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.
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.
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.