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:2412.03641 · Proceedings of the National Academy of Sciences

A solvable model for strongly interacting nonequilibrium excitons

Zhenhao Song, Tessa Cookmeyer, Leon Balents

Published Mar 14, 2025• 1 citation
arXivHTMLPDF

Abstract

We study the driven-dissipative Bose-Hubbard model with an all-to-all hopping term in the system Hamiltonian, while subject to incoherent pumping and decay from the environment. This system is naturally probed in several recent experiments on excitons in WS 2 /WSe 2 moiré systems, as well as quantum simulators. By positing a particular form of coupling to the environment, we derive the Lindblad jump operators and show that, in certain limits, the system admits a closed-form expression for the steady-state density matrix. Away from the exactly solvable regions, the steady state can be obtained numerically for 100s to 1,000s of sites. We study the nonequilibrium phase diagram and phase transitions, which qualitatively matches the equilibrium phase diagram, agreeing with the intuition that increasing the intensity of the light is equivalent to changing the bosonic chemical potential. However, the steady states are far from thermal states, and the nature of the phase transitions is changed.

DOI
10.1073/pnas.2424663122

Related papers

Symmetry between repulsive and attractive interactions in driven-dissipative Bose-Hubbard systems

Adil A. Gangat, Ian P. McCulloch, Ying-Jer Kao

The driven-dissipative Bose-Hubbard model can be experimentally realized with either negative or positive onsite detunings, inter-site hopping energies, and onsite interaction energies. Here we use one-dimensional matrix product density operators to perform a fully quantum investigation of the dependence of the non-equilibrium steady states of this model on the signs of these parameters. Due to a symmetry in the Lindblad master equation, we find that simultaneously changing the sign of the interaction energies, hopping energies, and chemical potentials leaves the local boson number distribution and inter-site number correlations invariant, and the steady-state complex conjugated. This shows that all driven-dissipative phenomena of interacting bosons described by the Lindblad master equation, such as "fermionization" and "superbunching", can equivalently occur with attractive or repulsive interactions.

arXiv:1706.07349•Scientific Reports•Feb 21, 2018•2 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

Nonequilibrium phase diagram of a driven and dissipative many-body system

Andrea Tomadin, Sebastian Diehl, P. Zoller

We study the nonequilibrium dynamics of a many-body bosonic system on a lattice, subject to driving and dissipation. The time evolution is described by a master equation, which we treat within a generalized Gutzwiller mean field approximation for density matrices. The dissipative processes are engineered such that the system, in the absence of interaction between the bosons, is driven into a homogeneous steady state with off-diagonal long-range order. We investigate how the coherent interaction affects the properties of the steady state of the system qualitatively and derive a nonequilibrium phase diagram featuring a phase transition into a steady state without long-range order. The phase diagram also exhibits an extended domain where an instability of the homogeneous steady state gives rise to a persistent density pattern with spontaneously broken translational symmetry. In the limit of low particle density, we provide a precise analytical description of the time evolution during the instability. Moreover, we investigate the transient following a quantum quench of the dissipative processes and we elucidate the prominent role played by collective topological variables in this regime.

arXiv:1011.3207•Physical Review A•Jan 18, 2011•110 citations

One-dimensional Bose-Hubbard model far from equilibrium

Iztok Pižorn

We study the nonequilibrium steady state of the one-dimensional Bose-Hubbard model coupled to Lindblad reservoirs, using the density-matrix renormalization group in operator space. We observe a transition from a flat particle density profile in the noninteracting limit to the linearly increasing profile with onset of the interparticle interaction. Analyzing the effect of coherent pumping on the nonequilibrium steady state, we find a subspace which remains unaffected by the pumping in the noninteracting limit with the protection gradually diminishing due to interparticle interaction. In the equilibrium situation with one or more symmetric reservoirs, we show analytically that the steady state of the system is a product state for any interaction strength. We also provide analytical results in the noninteracting limits, using the method of the third quantization in operator space.

arXiv:1308.3195•Physical Review A•Oct 25, 2013•34 citations

Quantum dynamics of disordered bosons in an optical lattice

Chien-Hung Lin, Rajdeep Sensarma, K. Sengupta +1 more

We study the equilibrium and nonequilibrium properties of strongly interacting bosons on a lattice in the presence of a random bounded disorder potential. Using a Gutzwiller projected variational technique, we study the equilibrium phase diagram of the disordered Bose-Hubbard model and obtain the Mott insulator, Bose glass, and superfluid phases. We also study the nonequilibrium response of the system under a periodic temporal drive where, starting from the superfluid phase, the hopping parameter is ramped down linearly in time, and back to its initial value. We study the density of excitations created, the change in the superfluid order parameter, and the energy pumped into the system in this process as a function of the inverse ramp rate $\ensuremath{\tau}$. For the clean case the density of excitations goes to a constant, while the order parameter and energy relax as $1/\ensuremath{\tau}$ and $1/{\ensuremath{\tau}}^{2}$ respectively. With disorder, the excitation density decays exponentially with $\ensuremath{\tau}$, with the decay rate increasing with the disorder, to an asymptotic value independent of the disorder. The energy and change in order parameter also decrease as $\ensuremath{\tau}$ is increased.

arXiv:1208.4852•Physical Review B•Dec 26, 2012•12 citations

Emergent equilibrium in many-body optical bistability

Michael Foss‐Feig, Pradeep Niroula, Jeremy T. Young +4 more

Many-body systems constructed of quantum-optical building blocks can now be realized in experimental platforms ranging from exciton-polariton fluids to ultracold Rydberg gases, establishing a fascinating interface between traditional many-body physics and the driven-dissipative, nonequilibrium setting of cavity QED. At this interface, the standard techniques and intuitions of both fields are called into question, obscuring issues as fundamental as the role of fluctuations, dimensionality, and symmetry on the nature of collective behavior and phase transitions. Here, we study the driven-dissipative Bose-Hubbard model, a minimal description of numerous atomic, optical, and solid-state systems in which particle loss is countered by coherent driving. Despite being a lattice version of optical bistability, a foundational and patently nonequilibrium model of cavity QED, the steady state possesses an emergent equilibrium description in terms of a classical Ising model. We establish this picture by making new connections between traditional techniques from many-body physics (functional integrals) and quantum optics (the system-size expansion). To lowest order in a controlled expansion---organized around the experimentally relevant limit of weak interactions---the full quantum dynamics reduces to nonequilibrium Langevin equations, which support a phase transition described by model A of the Hohenberg-Halperin classification. Numerical simulations of the Langevin equations corroborate this picture, revealing that canonical behavior associated with the Ising model manifests readily in simple experimental observables.

arXiv:1611.02284•Physical Review A•Apr 17, 2017•139 citations
We study the driven-dissipative Bose-Hubbard model with an all-to-all hopping term in the system Hamiltonian, while subject to incoherent pumping and decay from the environment. This system is naturally probed in several recent experiments on excitons in WS 2 /WSe 2 moiré systems, as well as quantum simulators. By positing a particular form of coupling to the environment, we derive the Lindblad jump operators and show that, in certain limits, the system admits a closed-form expression for the steady-state density matrix. Away from the exactly solvable regions, the steady state can be obtained numerically for 100s to 1,000s of sites. We study the nonequilibrium phase diagram and phase transitions, which qualitatively matches the equilibrium phase diagram, agreeing with the intuition that increasing the intensity of the light is equivalent to changing the bosonic chemical potential. However, the steady states are far from thermal states, and the nature of the phase transitions is changed.
The driven-dissipative Bose-Hubbard model can be experimentally realized with either negative or positive onsite detunings, inter-site hopping energies, and onsite interaction energies. Here we use one-dimensional matrix product density operators to perform a fully quantum investigation of the dependence of the non-equilibrium steady states of this model on the signs of these parameters. Due to a symmetry in the Lindblad master equation, we find that simultaneously changing the sign of the interaction energies, hopping energies, and chemical potentials leaves the local boson number distribution and inter-site number correlations invariant, and the steady-state complex conjugated. This shows that all driven-dissipative phenomena of interacting bosons described by the Lindblad master equation, such as "fermionization" and "superbunching", can equivalently occur with attractive or repulsive interactions.
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.
We study the nonequilibrium dynamics of a many-body bosonic system on a lattice, subject to driving and dissipation. The time evolution is described by a master equation, which we treat within a generalized Gutzwiller mean field approximation for density matrices. The dissipative processes are engineered such that the system, in the absence of interaction between the bosons, is driven into a homogeneous steady state with off-diagonal long-range order. We investigate how the coherent interaction affects the properties of the steady state of the system qualitatively and derive a nonequilibrium phase diagram featuring a phase transition into a steady state without long-range order. The phase diagram also exhibits an extended domain where an instability of the homogeneous steady state gives rise to a persistent density pattern with spontaneously broken translational symmetry. In the limit of low particle density, we provide a precise analytical description of the time evolution during the instability. Moreover, we investigate the transient following a quantum quench of the dissipative processes and we elucidate the prominent role played by collective topological variables in this regime.
We study the nonequilibrium steady state of the one-dimensional Bose-Hubbard model coupled to Lindblad reservoirs, using the density-matrix renormalization group in operator space. We observe a transition from a flat particle density profile in the noninteracting limit to the linearly increasing profile with onset of the interparticle interaction. Analyzing the effect of coherent pumping on the nonequilibrium steady state, we find a subspace which remains unaffected by the pumping in the noninteracting limit with the protection gradually diminishing due to interparticle interaction. In the equilibrium situation with one or more symmetric reservoirs, we show analytically that the steady state of the system is a product state for any interaction strength. We also provide analytical results in the noninteracting limits, using the method of the third quantization in operator space.
We study the equilibrium and nonequilibrium properties of strongly interacting bosons on a lattice in the presence of a random bounded disorder potential. Using a Gutzwiller projected variational technique, we study the equilibrium phase diagram of the disordered Bose-Hubbard model and obtain the Mott insulator, Bose glass, and superfluid phases. We also study the nonequilibrium response of the system under a periodic temporal drive where, starting from the superfluid phase, the hopping parameter is ramped down linearly in time, and back to its initial value. We study the density of excitations created, the change in the superfluid order parameter, and the energy pumped into the system in this process as a function of the inverse ramp rate $\ensuremath{\tau}$. For the clean case the density of excitations goes to a constant, while the order parameter and energy relax as $1/\ensuremath{\tau}$ and $1/{\ensuremath{\tau}}^{2}$ respectively. With disorder, the excitation density decays exponentially with $\ensuremath{\tau}$, with the decay rate increasing with the disorder, to an asymptotic value independent of the disorder. The energy and change in order parameter also decrease as $\ensuremath{\tau}$ is increased.
Many-body systems constructed of quantum-optical building blocks can now be realized in experimental platforms ranging from exciton-polariton fluids to ultracold Rydberg gases, establishing a fascinating interface between traditional many-body physics and the driven-dissipative, nonequilibrium setting of cavity QED. At this interface, the standard techniques and intuitions of both fields are called into question, obscuring issues as fundamental as the role of fluctuations, dimensionality, and symmetry on the nature of collective behavior and phase transitions. Here, we study the driven-dissipative Bose-Hubbard model, a minimal description of numerous atomic, optical, and solid-state systems in which particle loss is countered by coherent driving. Despite being a lattice version of optical bistability, a foundational and patently nonequilibrium model of cavity QED, the steady state possesses an emergent equilibrium description in terms of a classical Ising model. We establish this picture by making new connections between traditional techniques from many-body physics (functional integrals) and quantum optics (the system-size expansion). To lowest order in a controlled expansion---organized around the experimentally relevant limit of weak interactions---the full quantum dynamics reduces to nonequilibrium Langevin equations, which support a phase transition described by model A of the Hohenberg-Halperin classification. Numerical simulations of the Langevin equations corroborate this picture, revealing that canonical behavior associated with the Ising model manifests readily in simple experimental observables.