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

Fully Quantum Scalable Description of Driven-Dissipative Lattice Models

Piotr Deuar, Ferrier, A, Matuszewski, M, Giuliano Orso, Marzena Hanna Szymańska

Published Apr 1, 2021• 25 citations
OpenAlex

Abstract

Methods for modeling large driven-dissipative quantum systems are becoming increasingly urgent due to recent experimental progress in a number of photonic platforms. We demonstrate the positive-P method to be ideal for this purpose across a wide range of parameters, focusing on the archetypal driven-dissipative Bose-Hubbard model. Notably, these parameters include intermediate regimes where interactions and dissipation are comparable, and especially cases with low occupations for which common semiclassical approximations can break down. The presence of dissipation can alleviate instabilities in the method that are known to occur for closed systems, allowing the simulation of dynamics up to and including the steady state. Throughout the parameter space of the model, we determine the magnitude of dissipation that is sufficient to make the method useful and stable, finding its region of applicability to be complementary to that of the truncated Wigner method. We then demonstrate its use in a number of examples with nontrivial quantum correlations, including a demonstration of solving the urgent open problem of large and highly nonuniform systems with tens of thousands of sites.

Related papers

Fully Quantum Scalable Description of Driven-Dissipative Lattice Models

Piotr Deuar, Alex Ferrier, Michał Matuszewski +2 more

The problem of modeling large driven dissipative quantum systems is tackled with a phase space method, opening the way to progress on a myriad of experimental platforms described by a Bose Hubbard model, including highly non-uniform systems with tens of thousands of sites.

arXiv:2012.02014•PRX Quantum•Feb 5, 2021•37 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

Truncated correlation hierarchy schemes for driven-dissipative multimode quantum systems

Wim Casteels, Stefano Finazzi, Alexandre Le Boité +2 more

Abstract We present a method to describe driven-dissipative multi-mode systems by considering a truncated hierarchy of equations for the correlation functions. We consider two hierarchy truncation schemes with a global cutoff on the correlation order, which is the sum of the exponents of the operators involved in the correlation functions: a ‘hard’ cutoff corresponding to an expansion around the vacuum, which applies to a regime where the number of excitations per site is small; a ‘soft’ cutoff which corresponds to an expansion around coherent states, which can be applied for large excitation numbers per site. This approach is applied to describe the bunching-antibunching crossover in the driven-dissipative Bose–Hubbard model for photonic systems. The results have been successfully benchmarked by comparison with calculations based on the corner-space renormalization method in 1D and 2D systems. The regime of validity and strengths of the present truncation methods are critically discussed.

•Sep 2, 2016•17 citations

Space-time first-order correlations of an open Bose Hubbard model with incoherent pump and loss

Martina Zündel, Leonardo Mazza, Léonie Canet +1 more

We investigate the correlation properties in the steady state of driven-dissipative interacting bosonic systems in the quantum regime, as for example non-linear photonic cavities. Specifically, we consider the Bose-Hubbard model on a periodic chain and with spatially homogeneous one-body loss and pump within the Markovian approximation. The steady state is non-thermal and is formally equivalent to an infinite-temperature state with finite chemical potential set by the dissipative parameters. While there is no effect of interactions on the steady state, we observe a nontrivial behaviour of the space-time two-point correlation function, obtained by exact diagonalisation. In particular, we find that the decay width of the propagator is not only renormalised at increasing interactions, as it is the case of a single non-linear resonator, but also at increasing hopping strength. Furthermore, we numerically predict at large interactions a plateau value of the decay rate which goes beyond perturbative results in the interaction strength. We then compute the full spectral function, finding that it contains both a dispersive free-particle like dispersion at low energy and a doublon branch at energy corresponding to the on-site interactions. We compare with the corresponding calculation for the ground state of a closed quantum system and show that the driven-dissipative nature – determining both the steady state and the dynamical evolution – changes the low-lying part of the spectrum, where noticeably, the dispersion is quadratic instead of linear at small wavevectors. Finally, we compare to a high temperature grand-canonical equilibrium state and show the difference with respect to the open system stemming from the additional degree of freedom of the dissipation that allows one to vary the width of the dispersion lines.

arXiv:2405.19972•SciPost Physics•Mar 17, 2025•1 citation

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 correlation hierarchy schemes for driven-dissipative multimode quantum systems

Wim Casteels, Stefano Finazzi, Alexandre Le Boité +2 more

We present a method to describe driven-dissipative multi-mode systems by considering a truncated hierarchy of equations for the correlation functions. We consider two hierarchy truncation schemes with a global cutoff on the correlation order, which is the sum of the exponents of the operators involved in the correlation functions: a 'hard' cutoff corresponding to an expansion around the vacuum, which applies to a regime where the number of excitations per site is small; a 'soft' cutoff which corresponds to an expansion around coherent states, which can be applied for large excitation numbers per site. This approach is applied to describe the bunching-antibunching crossover in the driven-dissipative Bose–Hubbard model for photonic systems. The results have been successfully benchmarked by comparison with calculations based on the corner-space renormalization method in 1D and 2D systems. The regime of validity and strengths of the present truncation methods are critically discussed.

arXiv:1605.00882•New Journal of Physics•Sep 2, 2016•18 citations
Methods for modeling large driven-dissipative quantum systems are becoming increasingly urgent due to recent experimental progress in a number of photonic platforms. We demonstrate the positive-P method to be ideal for this purpose across a wide range of parameters, focusing on the archetypal driven-dissipative Bose-Hubbard model. Notably, these parameters include intermediate regimes where interactions and dissipation are comparable, and especially cases with low occupations for which common semiclassical approximations can break down. The presence of dissipation can alleviate instabilities in the method that are known to occur for closed systems, allowing the simulation of dynamics up to and including the steady state. Throughout the parameter space of the model, we determine the magnitude of dissipation that is sufficient to make the method useful and stable, finding its region of applicability to be complementary to that of the truncated Wigner method. We then demonstrate its use in a number of examples with nontrivial quantum correlations, including a demonstration of solving the urgent open problem of large and highly nonuniform systems with tens of thousands of sites.
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.
Abstract We present a method to describe driven-dissipative multi-mode systems by considering a truncated hierarchy of equations for the correlation functions. We consider two hierarchy truncation schemes with a global cutoff on the correlation order, which is the sum of the exponents of the operators involved in the correlation functions: a ‘hard’ cutoff corresponding to an expansion around the vacuum, which applies to a regime where the number of excitations per site is small; a ‘soft’ cutoff which corresponds to an expansion around coherent states, which can be applied for large excitation numbers per site. This approach is applied to describe the bunching-antibunching crossover in the driven-dissipative Bose–Hubbard model for photonic systems. The results have been successfully benchmarked by comparison with calculations based on the corner-space renormalization method in 1D and 2D systems. The regime of validity and strengths of the present truncation methods are critically discussed.
We investigate the correlation properties in the steady state of driven-dissipative interacting bosonic systems in the quantum regime, as for example non-linear photonic cavities. Specifically, we consider the Bose-Hubbard model on a periodic chain and with spatially homogeneous one-body loss and pump within the Markovian approximation. The steady state is non-thermal and is formally equivalent to an infinite-temperature state with finite chemical potential set by the dissipative parameters. While there is no effect of interactions on the steady state, we observe a nontrivial behaviour of the space-time two-point correlation function, obtained by exact diagonalisation. In particular, we find that the decay width of the propagator is not only renormalised at increasing interactions, as it is the case of a single non-linear resonator, but also at increasing hopping strength. Furthermore, we numerically predict at large interactions a plateau value of the decay rate which goes beyond perturbative results in the interaction strength. We then compute the full spectral function, finding that it contains both a dispersive free-particle like dispersion at low energy and a doublon branch at energy corresponding to the on-site interactions. We compare with the corresponding calculation for the ground state of a closed quantum system and show that the driven-dissipative nature – determining both the steady state and the dynamical evolution – changes the low-lying part of the spectrum, where noticeably, the dispersion is quadratic instead of linear at small wavevectors. Finally, we compare to a high temperature grand-canonical equilibrium state and show the difference with respect to the open system stemming from the additional degree of freedom of the dissipation that allows one to vary the width of the dispersion lines.
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 present a method to describe driven-dissipative multi-mode systems by considering a truncated hierarchy of equations for the correlation functions. We consider two hierarchy truncation schemes with a global cutoff on the correlation order, which is the sum of the exponents of the operators involved in the correlation functions: a 'hard' cutoff corresponding to an expansion around the vacuum, which applies to a regime where the number of excitations per site is small; a 'soft' cutoff which corresponds to an expansion around coherent states, which can be applied for large excitation numbers per site. This approach is applied to describe the bunching-antibunching crossover in the driven-dissipative Bose–Hubbard model for photonic systems. The results have been successfully benchmarked by comparison with calculations based on the corner-space renormalization method in 1D and 2D systems. The regime of validity and strengths of the present truncation methods are critically discussed.