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:1406.1192 · Physical Review Letters

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

Shainen M. Davidson, Anatoli Polkovnikov

Published Jan 29, 2015• 52 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevlett.114.045701

Related papers

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

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

Kazuma Nagao, Ippei Danshita, Seiji Yunoki

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.

arXiv:2307.16170•Physical Review A•Dec 11, 2024•5 citations

Representation of quantum dynamics of interacting systems through classical trajectories

Anatoli Polkovnikov

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.

•May 20, 2009•1 citation

Semiclassical representation of the Hubbard model

Yuki Yamasaki, Hidemaro Suwa, Cristian D. Batista +1 more

By revisiting the path-integral formulation of the Hubbard model, we propose a theoretical approach based on a semiclassical approximation employing an unconventional coherent-state representation. Within this framework, a subset of the dynamical variables is treated as static, yielding a nonperturbative scheme that is applicable at finite temperature, incorporates intersite correlations, and can be naturally extended to multiorbital systems. We assess the validity of the approximation by comparing its results with exact solutions for one- and two-site systems, focusing in particular on the particle number, double occupancy, hopping amplitude, and spin correlations, and find that the present approach qualitatively reproduces the exact behavior. Quantitatively, deviations arise, which is associated with the continuum (non-discretized) character of the underlying density of states. Furthermore, we derive the exact transformation associated with the coherent-state construction, thereby providing additional insight into the representation of the Hubbard model.

arXiv:2604.02769•Apr 3, 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

Quantum dynamics of Bose-Hubbard Hamiltonians beyond the Hartree-Fock-Bogoliubov approximation: The Bogoliubov back-reaction approximation

Igor Tikhonenkov, J. R. Anglin, Amichay Vardi

We formulate a method for studying the quantum field dynamics of ultracold Bose gases confined within optical lattice potentials, within the lowest Bloch-band Bose-Hubbard model. Our formalism extends the two-sites results of Vardi and Anglin [Phys. Rev. Lett. 86, 568 (2001)] to the general case of $M$ lattice sites. The methodology is based on mapping the Bose-Hubbard Hamiltonian to an $\mathrm{SU}(M)$ pseudospin problem and truncating the resulting hierarchy of dynamical equations for correlation functions, up to pair correlations between $\mathrm{SU}(M)$ generators. Agreement with few-site exact many-particle calculations is consistently better than the corresponding Hartree-Fock-Bogoliubov approximation. Moreover, our approximation compares favorably with a more elaborate two-particle irreducible effective action formalism, at a fraction of the analytic and numerical effort.

Physical Review A•Jan 10, 2007•36 citations
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 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 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 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.
By revisiting the path-integral formulation of the Hubbard model, we propose a theoretical approach based on a semiclassical approximation employing an unconventional coherent-state representation. Within this framework, a subset of the dynamical variables is treated as static, yielding a nonperturbative scheme that is applicable at finite temperature, incorporates intersite correlations, and can be naturally extended to multiorbital systems. We assess the validity of the approximation by comparing its results with exact solutions for one- and two-site systems, focusing in particular on the particle number, double occupancy, hopping amplitude, and spin correlations, and find that the present approach qualitatively reproduces the exact behavior. Quantitatively, deviations arise, which is associated with the continuum (non-discretized) character of the underlying density of states. Furthermore, we derive the exact transformation associated with the coherent-state construction, thereby providing additional insight into the representation of the Hubbard model.
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.
We formulate a method for studying the quantum field dynamics of ultracold Bose gases confined within optical lattice potentials, within the lowest Bloch-band Bose-Hubbard model. Our formalism extends the two-sites results of Vardi and Anglin [Phys. Rev. Lett. 86, 568 (2001)] to the general case of MMM lattice sites. The methodology is based on mapping the Bose-Hubbard Hamiltonian to an SU(M)\mathrm{SU}(M)SU(M) pseudospin problem and truncating the resulting hierarchy of dynamical equations for correlation functions, up to pair correlations between SU(M)\mathrm{SU}(M)SU(M) generators. Agreement with few-site exact many-particle calculations is consistently better than the corresponding Hartree-Fock-Bogoliubov approximation. Moreover, our approximation compares favorably with a more elaborate two-particle irreducible effective action formalism, at a fraction of the analytic and numerical effort.