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.
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 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.
We review the truncated Wigner approximation (TWA) for spins as a computationally inexpensive numerical approximation method to describe interacting and / or dissipative many-body spin systems. Using the Wigner-Moyal mapping from Hilbert space to a suitable phase space, the many-body density matrix is represented by a c-number distribution, the Wigner function. The gauge freedom in continuous phase space can be exploited to find positive Wigner functions for a large class of spin states, including entangled ones. Employing different sets of correspondence rules, we derive equations of motion for the Wigner function, which, applying controlled approximations, can be mapped to stochastic differential equations. This allows a computationally inexpensive simulation of expectation values. Using a phase-space analog of the quantum regression theorem also multi-time correlations and spectra can be obtained. To illustrate the potential of the method, we benchmark the TWA for spins with some exactly solvable problems of interacting, dissipative spin systems, and then discuss its application to collective processes, such as the superradiant emission of light. Extending the TWA to imaginary time furthermore provides a tool to approximately calculate thermal and ground states of spin Hamiltonians. Finally, we show that the TWA stochastic equations can equivalently be derived within a path-integral approach, provided that the operator products in the dissipator are rigorously mapped onto the curved phase space.
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.
The truncated Wigner approximation (TWA) enables us to investigate bosonic quantum many-body dynamics, including open quantum systems described by the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) equation. In the TWA, the Weyl-Wigner transformation, a way of mapping from quantum-mechanical operators to c -numbers, of the GKSL equation leads to the Fokker–Planck equation, which we calculate by reducing it to the corresponding stochastic differential equations. However, the Fokker–Planck equation is not always reduced to the stochastic differential equations depending on details of jump operators. In this work, we clarify the condition for obtaining the stochastic differential equations from the Fokker–Planck equation and derive analytical expressions of these equations for a system with an arbitrary Hamiltonian with jump operators that do not couple different states. This result leads to remarkable simplification of the complicated calculations required in the conventional procedures of the TWA. In the course of the derivation, we formulate the GKSL equation by using the path-integral representation based on the Weyl-Wigner transformation, which gives us a clear interpretation of the relation between the TWA and quantum fluctuations and allows us to calculate the non-equal time correlation functions in the TWA. In the benchmark calculations, we numerically confirm that the relaxation dynamics of physical quantities including the non-equal time correlation functions obtained in our formulation agrees well with the exact ones in the numerically solvable models. Since our formulation is widely applicable to bosonic systems such as the Bose-Hubbard model, it facilitates studies of the open quantum many-body phenomenon.
Long-range interactions can protect dissipative order against fluctuations, but determining when they do so in two dimensions remains out of reach for existing many-body methods. We develop a continuous SU(2)-covariant truncated-Wigner framework that maps the problem onto stochastic dynamics with one Bloch sphere per site. Derived systematically from the exact Stratonovich-Weyl evolution of the Lindblad equation, it generates both drift and reservoir noise for arbitrary spin S. Benchmarks against exact and tensor-network calculations validate its description of non-equilibrium dynamics and steady states. Applied to 400-site two-dimensional dissipative anisotropic XY lattices, our method reveals physics missed by mean-field and perturbative approaches: vortex-antivortex proliferation destroys the predicted staggered-XY phase for nearest-neighbour coupling, whereas algebraically decaying interactions stabilize an extended Z2-ordered phase. We obtain the first beyond-mean-field steady-state phase diagram of the long-range interacting model, revealing reentrant loss of order at strong coupling and closure of the ordered region beyond α≈3.5. Our new method provides a scalable route to strongly correlated driven-dissipative systems with long-range interactions, higher-dimensional geometries, spin-boson couplings, and spins beyond one half, relevant to modern experiments with magnetic atoms and polar molecules in optical lattices and tweezer arrays.
We put forward a user-friendly framework of the truncated Wigner approximation (TWA) for dissipative quantum many-body systems. Our approach is computationally affordable and it features a straightforward implementation. The leverage of the method can be ultimately traced to an intimate connection between the TWA and the semiclassical limit of the quantum Langevin equation, which we unveil by resorting to a path-integral representation of the Lindbladian. Our approach allows us to explore dynamics from early to late times in a variety of models at the core of modern atomic, molecular, and optical research, including lasing, central-spin models, driven arrays of Rydbergs, and correlated emission in free space. Notably, our TWA approach outperforms the cumulant-expansion method in certain models and performs comparably well in others, all while offering significantly lower computational costs and a much simpler formulation of the dynamical equations. We therefore argue that TWA could, in the near future, become a primary tool for a fast and efficient first exploration of driven-dissipative many-body dynamics on consumer-grade computers.
We present a semiclassical phase-space method to calculate thermal and ground states of large interacting spin systems. To this end, we extend the recently developed truncated Wigner approximation for spins (TWA) to the imaginary time, termed iTWA. The evolution of the canonical density matrix in imaginary time is mapped to a partial differential equation of its Wigner function. Truncation at the Fokker-Planck level leads to a set of stochastic differential equations, which can be efficiently simulated even for large systems. We show that for general Ising Hamiltonians the approximation becomes exact for large imaginary times subject only to sampling errors. Thus the iTWA is ideal to determine the ground state of spin glasses or to find solutions to quadratic unconstrained binary optimization problems (QUBO) on a controlled approximation level. We illustrate this for MaxCut on random, unweighted 3-regular graphs, encoded in an anti-ferromagnetic Ising Hamiltonian, for which finding the exact ground state and even approximations to it beyond a certain accuracy is know to be NP hard. Furthermore, in order to assess the quality of the method also for general spin models, we analyze the ground-state quantum phase transition of the transverse-field Ising model in one and two spatial dimensions, finding reasonably good agreement with the exact behavior.