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:2608.00341 · Advances in atomic, molecular, and optical physics

Phase-Space Methods for Many-Body Quantum Optics

Edgar Guardiola-Navarrete, Silvia Cardenas-Lopez, Ana Asenjo-Garcia

quant-phSubmitted Jul 31, 2026• Announced Aug 4, 2026• 2 citations
arXivHTMLPDF

Abstract

Many-body quantum-optical systems, where a collection of emitters interacts through a common electromagnetic reservoir, exhibit rich out-of-equilibrium behavior and hold promise for applications in quantum technologies. However, exact numerical simulations of their dynamics quickly become unfeasible due to the exponential growth of the Hilbert space with system size. Semiclassical, phase-space approaches -- such as the Truncated Wigner approximation (TWA) -- provide computationally efficient alternatives by capturing leading-order quantum fluctuations. In this paper, we present a comprehensive overview of how to tackle problems in many-body quantum optics using phase-space methods. We derive the exact partial differential equation governing many-body dissipative evolution in any phase-space representation and discuss the approximations that yield the dissipative TWA proposed by Mink and Fleischhauer [SciPost Phys. 15, 233 (2023)]. We find that $P$ and $Q$ distributions are generally suboptimal for many-body quantum optics. Additionally, we extend the formalism to calculate multi-time correlation functions, thereby broadening the scope of phase-space simulations of open spin systems to include coherence and spectral properties, as well as directional correlations of collectively radiating emitters. These developments provide valuable tools for investigating exotic light sources driven by collective dissipation, driven-dissipative phase transitions, and a wealth of many-body phenomena arising in state-of-the-art experimental platforms.

Journal reference
Adv. At. Mol. Opt. Phys. 74, 87 (2025)
DOI
10.1016/bs.aamop.2025.04.003
Comments
51 pages, 12 figures

Related papers

Truncated Wigner approximation for spins in continuous phase space

Jens Hartmann, Tom Schlegel, Viktoria Noel +2 more

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.

quant-pharXiv:2607.17295•Jul 21, 2026•0 citations

User-Friendly Truncated Wigner Approximation for Dissipative Spin Dynamics

Hossein Hosseinabadi, Oksana Chelpanova, Jamir Marino

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.

arXiv:2503.17443•PRX Quantum•Aug 1, 2025•13 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

Collective radiative interactions in the discrete truncated Wigner approximation

Christopher D. Mink, Michael Fleischhauer

Interfaces of light and matter serve as a platform for exciting many-body physics and photonic quantum technologies. Due to the recent experimental realization of atomic arrays at sub-wavelength spacings, collective interaction effects such as superradiance have regained substantial interest. Their analytical and numerical treatment is however quite challenging. Here we develop a semiclassical approach to this problem that allows to describe the coherent and dissipative many-body dynamics of interacting spins while taking into account lowest-order quantum fluctuations. For this purpose we extend the discrete truncated Wigner approximation, originally developed for unitarily coupled spins, to include collective, dissipative spin processes by means of truncated correspondence rules. This maps the dynamics of the atomic ensemble onto a set of semiclassical, numerically inexpensive stochastic differential equations. We benchmark our method with exact results for the case of Dicke decay, which shows excellent agreement. For small arrays we compare to exact simulations, again showing good agreement at early times and at moderate to strong driving, and to a second order cumulant expansion. We conclude by studying the radiative properties of a spatially extended three-dimensional, coherently driven gas and compare the coherence of the emitted light to experimental results.

arXiv:2305.19829•SciPost Physics•Dec 8, 2023•32 citations

Phase-space methods for simulating the dissipative many-body dynamics of collective spin systems

Julian Huber, Peter Kirton, Peter Rabl

We describe an efficient numerical method for simulating the dynamics and steady states of collective spin systems in the presence of dephasing and decay. The method is based on the Schwinger boson representation of spin operators and uses an extension of the truncated Wigner approximation to map the exact open system dynamics onto stochastic differential equations for the corresponding phase space distribution. This approach is most effective in the limit of very large spin quantum numbers, where exact numerical simulations and other approximation methods are no longer applicable. We benchmark this numerical technique for known superradiant decay and spin-squeezing processes and illustrate its application for the simulation of non-equilibrium phase transitions in dissipative spin lattice models.

arXiv:2011.10049•SciPost Physics•Feb 22, 2021•37 citations

Phase-space methods for simulating the dissipative many-body dynamics of collective spin systems

Julian Huber, Peter Kirton, Peter Rabl

We describe an efficient numerical method for simulating the dynamics and steady states of collective spin systems in the presence of dephasing and decay. The method is based on the Schwinger boson representation of spin operators and uses an extension of the truncated Wigner approximation to map the exact open system dynamics onto stochastic differential equations for the corresponding phase space distribution. This approach is most effective in the limit of very large spin quantum numbers, where exact numerical simulations and other approximation methods are no longer applicable. We benchmark this numerical technique for known superradiant decay and spin-squeezing processes and illustrate its application for the simulation of non-equilibrium phase transitions in dissipative spin lattice models.

•Feb 22, 2021•34 citations
Many-body quantum-optical systems, where a collection of emitters interacts through a common electromagnetic reservoir, exhibit rich out-of-equilibrium behavior and hold promise for applications in quantum technologies. However, exact numerical simulations of their dynamics quickly become unfeasible due to the exponential growth of the Hilbert space with system size. Semiclassical, phase-space approaches -- such as the Truncated Wigner approximation (TWA) -- provide computationally efficient alternatives by capturing leading-order quantum fluctuations. In this paper, we present a comprehensive overview of how to tackle problems in many-body quantum optics using phase-space methods. We derive the exact partial differential equation governing many-body dissipative evolution in any phase-space representation and discuss the approximations that yield the dissipative TWA proposed by Mink and Fleischhauer [SciPost Phys. 15, 233 (2023)]. We find that PPP and QQQ distributions are generally suboptimal for many-body quantum optics. Additionally, we extend the formalism to calculate multi-time correlation functions, thereby broadening the scope of phase-space simulations of open spin systems to include coherence and spectral properties, as well as directional correlations of collectively radiating emitters. These developments provide valuable tools for investigating exotic light sources driven by collective dissipation, driven-dissipative phase transitions, and a wealth of many-body phenomena arising in state-of-the-art experimental platforms.
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 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.
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.
Interfaces of light and matter serve as a platform for exciting many-body physics and photonic quantum technologies. Due to the recent experimental realization of atomic arrays at sub-wavelength spacings, collective interaction effects such as superradiance have regained substantial interest. Their analytical and numerical treatment is however quite challenging. Here we develop a semiclassical approach to this problem that allows to describe the coherent and dissipative many-body dynamics of interacting spins while taking into account lowest-order quantum fluctuations. For this purpose we extend the discrete truncated Wigner approximation, originally developed for unitarily coupled spins, to include collective, dissipative spin processes by means of truncated correspondence rules. This maps the dynamics of the atomic ensemble onto a set of semiclassical, numerically inexpensive stochastic differential equations. We benchmark our method with exact results for the case of Dicke decay, which shows excellent agreement. For small arrays we compare to exact simulations, again showing good agreement at early times and at moderate to strong driving, and to a second order cumulant expansion. We conclude by studying the radiative properties of a spatially extended three-dimensional, coherently driven gas and compare the coherence of the emitted light to experimental results.
We describe an efficient numerical method for simulating the dynamics and steady states of collective spin systems in the presence of dephasing and decay. The method is based on the Schwinger boson representation of spin operators and uses an extension of the truncated Wigner approximation to map the exact open system dynamics onto stochastic differential equations for the corresponding phase space distribution. This approach is most effective in the limit of very large spin quantum numbers, where exact numerical simulations and other approximation methods are no longer applicable. We benchmark this numerical technique for known superradiant decay and spin-squeezing processes and illustrate its application for the simulation of non-equilibrium phase transitions in dissipative spin lattice models.
We describe an efficient numerical method for simulating the dynamics and steady states of collective spin systems in the presence of dephasing and decay. The method is based on the Schwinger boson representation of spin operators and uses an extension of the truncated Wigner approximation to map the exact open system dynamics onto stochastic differential equations for the corresponding phase space distribution. This approach is most effective in the limit of very large spin quantum numbers, where exact numerical simulations and other approximation methods are no longer applicable. We benchmark this numerical technique for known superradiant decay and spin-squeezing processes and illustrate its application for the simulation of non-equilibrium phase transitions in dissipative spin lattice models.