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

Novel Phase-Space Methods to Simulate Strongly-Interacting Many-Body Quantum Dynamics

Shainen M. Davidson

Published Jan 1, 2017• 2 citations
OpenAlex

Abstract

Understanding the collective behaviour of many-body quantum systems is an important subject in many areas of physics. With advances in ultra-cold gas experiments, the dynamics of strongly-interacting systems can now be studied in the lab. However, there is a paucity of theoretical techniques available to simulate such systems. One technique is phase-space methods, often known as the Truncated Wigner Approximation; however, its applicability in its naive form is limited. In this work, we expound on techniques to expand the regimes in which it can be effective. This involves creating a novel phase-space that is tailored to the problem at hand, and associated classical equations of motion. We show techniques for lattice systems with local finite Hilbert spaces, for fermionic systems, and for many-body localized systems. In all cases, we benchmark the accuracy of the approximation against exact results.

Related papers

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

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

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

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

The truncated Wigner approximation on discrete phase spaces

Johannes Schachenmayer, Alexander Pikovski, Ana María Rey

A phase-space method, suitable to study the many-body spin dynamics in generic spin lattice models is introduced, discussed, and tested. This technique, which we refer to as the dTWA, is based on the discrete sampling of the Wigner function instead of the continuous sampling in traditional Truncated Wigner Approximation (TWA) methods. The power of the method is demonstrated by comparisons with analytical and numerically exact calculations, which show a significant improvement compared to the TWA, e.g. by capturing revivals of spin-observables. The dTWA is demonstrated to not only correctly reproduce dynamics of one-particle observables but also two-point correlations and spin squeezing, a property relying on entanglement. This method thus opens the possibility to study the quantum dynamics in long-range spin models, accessible to recent state-of-the-art experiments, in regimes where other numerical techniques are inapplicable.

•Aug 19, 2014•0 citations

Exact and rigorous methods in quantum many body physics

Zhiyuan Wang

Quantum many body physics is an exciting research area, involving novel phases of matter with fundamentally new properties, but is also notoriously hard due to complexity of interacting quantum systems. Some popular approaches involve approximation techniques and numerical simulations, which are known to fail in several important cases. In contrast, rigorous mathematical tools, such as exact solutions and operator inequalities, have a narrower range of applicability, but provide guaranteed results and insights into the underlying physical mechanisms. This research aims to develop new techniques in this direction and use them to explore novel phases of matter. My first direction is to construct toy models where exact solutions are possible. Such models are important as they prove that certain physical phenomena are theoretically possible in nature, and often lead to the discovery of new phases. Such insights are provided by three exactly solvable models I have discovered: (1) a family of 1D quantum spin models hosting free parastatistical quasiparticles (an exotic type of identical particles beyond fermions, bosons, and anyons), proving for the first time that parastatistics is theoretically possible as an emergent phenomenon; (2) a 3D classical Ising model whose phases are characterized by topological features of certain loop observables, suggesting existence of previously unknown classical phases and phase transitions with topological order parameters; and (3) a family of models with exact p-wave superconducting ground states demonstrating the existence of Majorana quasiparticles and non-Abelian statistics in particle number-conserving systems. My second direction is to derive rigorous bounds and exact constraints on physical observables, which are applicable to large families of quantum many-body systems. I present three directions of progress: (1) a method that dramatically improves the upper bounds on the speed of information propagation in locally-interacting systems, which significantly extends the scope of these bounds and enables new applications; (2) bounds on finite-size errors in numerical simulations of lattice systems, including quench dynamics and gapped ground states; and (3) a locality bound on gapped ground states of power-law interacting systems, which leads to a generalization of the aforementioned error bounds to such systems. These error bounds have important theoretical implications such as proving the existence of the thermodynamic limit and stability of phases, and are practically useful in determining the validity of finite-size numerical simulations.

•Aug 1, 2022•0 citations
Understanding the collective behaviour of many-body quantum systems is an important subject in many areas of physics. With advances in ultra-cold gas experiments, the dynamics of strongly-interacting systems can now be studied in the lab. However, there is a paucity of theoretical techniques available to simulate such systems. One technique is phase-space methods, often known as the Truncated Wigner Approximation; however, its applicability in its naive form is limited. In this work, we expound on techniques to expand the regimes in which it can be effective. This involves creating a novel phase-space that is tailored to the problem at hand, and associated classical equations of motion. We show techniques for lattice systems with local finite Hilbert spaces, for fermionic systems, and for many-body localized systems. In all cases, we benchmark the accuracy of the approximation against exact results.
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.
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 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.
A phase-space method, suitable to study the many-body spin dynamics in generic spin lattice models is introduced, discussed, and tested. This technique, which we refer to as the dTWA, is based on the discrete sampling of the Wigner function instead of the continuous sampling in traditional Truncated Wigner Approximation (TWA) methods. The power of the method is demonstrated by comparisons with analytical and numerically exact calculations, which show a significant improvement compared to the TWA, e.g. by capturing revivals of spin-observables. The dTWA is demonstrated to not only correctly reproduce dynamics of one-particle observables but also two-point correlations and spin squeezing, a property relying on entanglement. This method thus opens the possibility to study the quantum dynamics in long-range spin models, accessible to recent state-of-the-art experiments, in regimes where other numerical techniques are inapplicable.
Quantum many body physics is an exciting research area, involving novel phases of matter with fundamentally new properties, but is also notoriously hard due to complexity of interacting quantum systems. Some popular approaches involve approximation techniques and numerical simulations, which are known to fail in several important cases. In contrast, rigorous mathematical tools, such as exact solutions and operator inequalities, have a narrower range of applicability, but provide guaranteed results and insights into the underlying physical mechanisms. This research aims to develop new techniques in this direction and use them to explore novel phases of matter. My first direction is to construct toy models where exact solutions are possible. Such models are important as they prove that certain physical phenomena are theoretically possible in nature, and often lead to the discovery of new phases. Such insights are provided by three exactly solvable models I have discovered: (1) a family of 1D quantum spin models hosting free parastatistical quasiparticles (an exotic type of identical particles beyond fermions, bosons, and anyons), proving for the first time that parastatistics is theoretically possible as an emergent phenomenon; (2) a 3D classical Ising model whose phases are characterized by topological features of certain loop observables, suggesting existence of previously unknown classical phases and phase transitions with topological order parameters; and (3) a family of models with exact p-wave superconducting ground states demonstrating the existence of Majorana quasiparticles and non-Abelian statistics in particle number-conserving systems. My second direction is to derive rigorous bounds and exact constraints on physical observables, which are applicable to large families of quantum many-body systems. I present three directions of progress: (1) a method that dramatically improves the upper bounds on the speed of information propagation in locally-interacting systems, which significantly extends the scope of these bounds and enables new applications; (2) bounds on finite-size errors in numerical simulations of lattice systems, including quench dynamics and gapped ground states; and (3) a locality bound on gapped ground states of power-law interacting systems, which leads to a generalization of the aforementioned error bounds to such systems. These error bounds have important theoretical implications such as proving the existence of the thermodynamic limit and stability of phases, and are practically useful in determining the validity of finite-size numerical simulations.