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:2305.19829 · SciPost Physics

Collective radiative interactions in the discrete truncated Wigner approximation

Christopher D. Mink, Michael Fleischhauer

Published Dec 8, 2023• 32 citations
arXivHTMLPDF

Abstract

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.

DOI
10.21468/scipostphys.15.6.233

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

Phase-Space Methods for Many-Body Quantum Optics

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

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.

quant-pharXiv:2608.00341•Advances in atomic, molecular, and optical physics•Aug 4, 2026•2 citations

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

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

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

Light-Induced Collective Interactions in Arrays of Quantum Emitters

Oriol Rubies-Bigordà

The interaction between light and matter has captivated physicists for centuries, from early studies of vision and refraction in ancient Greece to the development of quantum mechanics and quantum electrodynamics in the past century. While the response of a single quantum emitter to light is well understood, the radiative properties of an ensemble of closely spaced emitters are far more intricate. Coupling to a shared electromagnetic environment induces coherent and dissipative interactions between emitters, giving rise to a collective response that cannot be captured by treating them independently. In the regime of few excitations, the system hosts delocalized subradiant states, that is, coherent superpositions that are largely decoupled from the electromagnetic field and thus decay at suppressed rates. While this weak coupling makes subradiant states attractive for quantum technologies, it also renders them difficult to manipulate. At higher excitation densities, the intrinsic nonlinearity of emitters and the exponential growth of the Hilbert space make theoretical and numerical descriptions of the system and its dynamics increasingly challenging. This thesis explores two fundamental questions: How can subradiant and dark states be selectively accessed and harnessed for practical applications in quantum technologies? And how can interacting ensembles of quantum emitters be efficiently simulated to uncover their many-body physics? The first part of the thesis develops protocols for controlling and addressing dark states in free-space and waveguide-coupled atomic arrays, demonstrating their utility in quantum storage and the deterministic generation of entangled photonic states. Incorporating atomic motion, we further show that collective subradiant states can enhance cooling in dense atomic arrays, offering new avenues for controlling motional dynamics. In the second part, we introduce cumulant expansions of the equations of motion as a powerful tool to analytically and numerically investigate collective decay in the many-body regime. We first examine the collective decay of fully excited atomic arrays in free space, characterizing the onset and scaling of the superradiant burst across different geometries. In collaboration with experiments on ultracold erbium atoms in optical lattices, we provide the first direct observations of many-body collective effects in free-space ordered arrays, including early-time superradiant bursts, late-time subradiant tails, and the emergence of atomic correlations throughout the dynamics. Finally, we theoretically and numerically explore the transient formation of multi-excitation subradiant states, and demonstrate how the existence of multiple dissipation channels suppresses steady-state superradiance in extended arrays.

•May 1, 2025•0 citations
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 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.
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 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.
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.
The interaction between light and matter has captivated physicists for centuries, from early studies of vision and refraction in ancient Greece to the development of quantum mechanics and quantum electrodynamics in the past century. While the response of a single quantum emitter to light is well understood, the radiative properties of an ensemble of closely spaced emitters are far more intricate. Coupling to a shared electromagnetic environment induces coherent and dissipative interactions between emitters, giving rise to a collective response that cannot be captured by treating them independently. In the regime of few excitations, the system hosts delocalized subradiant states, that is, coherent superpositions that are largely decoupled from the electromagnetic field and thus decay at suppressed rates. While this weak coupling makes subradiant states attractive for quantum technologies, it also renders them difficult to manipulate. At higher excitation densities, the intrinsic nonlinearity of emitters and the exponential growth of the Hilbert space make theoretical and numerical descriptions of the system and its dynamics increasingly challenging. This thesis explores two fundamental questions: How can subradiant and dark states be selectively accessed and harnessed for practical applications in quantum technologies? And how can interacting ensembles of quantum emitters be efficiently simulated to uncover their many-body physics? The first part of the thesis develops protocols for controlling and addressing dark states in free-space and waveguide-coupled atomic arrays, demonstrating their utility in quantum storage and the deterministic generation of entangled photonic states. Incorporating atomic motion, we further show that collective subradiant states can enhance cooling in dense atomic arrays, offering new avenues for controlling motional dynamics. In the second part, we introduce cumulant expansions of the equations of motion as a powerful tool to analytically and numerically investigate collective decay in the many-body regime. We first examine the collective decay of fully excited atomic arrays in free space, characterizing the onset and scaling of the superradiant burst across different geometries. In collaboration with experiments on ultracold erbium atoms in optical lattices, we provide the first direct observations of many-body collective effects in free-space ordered arrays, including early-time superradiant bursts, late-time subradiant tails, and the emergence of atomic correlations throughout the dynamics. Finally, we theoretically and numerically explore the transient formation of multi-excitation subradiant states, and demonstrate how the existence of multiple dissipation channels suppresses steady-state superradiance in extended arrays.