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:2302.02977 · Physical Review A

Nonequilibrium steady states and critical slowing down in the dissipative Bose-Hubbard model

Robbe Ceulemans, Michiel Wouters

Published Jul 19, 2023• 6 citations
arXivHTMLPDF

Abstract

Motivated by recent experiments, we study the properties of large Bose-Hubbard chains with single-particle losses at one site using classical field methods. We construct and validate a compact effective model that reduces computations to only a few sites. We show that in the mean-field approach the description captures the stationary states of the dissipative mode very well. Not only is there good quantitative agreement in the hysteresis loop, but the dark soliton state can be reproduced as well due to the preservation of the $\text{U}(1)$ symmetry. Bimodality of the steady states, observed on longer timescales, is studied using the truncated Wigner method. We compare the switching statistics and derive the effective Liouvillian gap as a function of the tunneling, showing that the effective description underestimates fluctuations.

DOI
10.1103/physreva.108.013314

Related papers

Three-site Bose-Hubbard model subject to atom losses: Boson-pair dissipation channel and failure of the mean-field approach

V. S. Shchesnovich, Dmitri S. Mogilevtsev

We employ the perturbation series expansion for derivation of the reduced master equations for the three-site Bose-Hubbard model subject to strong atom losses from the central site. The model describes a condensate trapped in a triple-well potential subject to externally controlled removal of atoms. We find that the $\ensuremath{\pi}$-phase state of the coherent superposition between the side wells decays via two dissipation channels, the single-boson channel (similar to the externally applied dissipation) and the boson-pair channel. The quantum derivation is compared to the classical adiabatic elimination within the mean-field approximation. We find that the boson-pair dissipation channel is not captured by the mean-field model, whereas the single-boson channel is described by it. Moreover, there is a matching condition between the zero-point energy bias of the side wells and the nonlinear interaction parameter which separates the regions where either the single-boson or the boson-pair dissipation channel dominate. Our results indicate that the $M$-site Bose-Hubbard models, for $M>2$, subject to atom losses may require an analysis which goes beyond the usual mean-field approximation for correct description of their dissipative features. This is an important result in view of the recent experimental works on the single-site addressability of condensates trapped in optical lattices.

arXiv:1007.2544•Physical Review A•Oct 25, 2010•35 citations

Dynamical hysteresis properties of the driven-dissipative Bose-Hubbard model with a Gutzwiller Monte Carlo approach

Dolf Huybrechts, Michiel Wouters

We study the dynamical properties of a driven-dissipative Bose-Hubbard model in the strongly interacting regime through a quantum trajectory approach with a cluster Gutzwiller Ansatz for the wave function. This allows us to take classical and quantum correlations into account. By studying the dynamical hysteresis surface that arises by sweeping through the coherent driving strength we show that the phase diagram for this system is in qualitative correspondence with the Gutzwiller mean-field result. However, quantitative differences are present and the inclusion of classical and quantum correlations causes a significant shift of the critical parameters. Additionally, we show that approximation techniques relying on a unimodal distribution such as the mean field and $1/z$ expansion drastically underestimate the particle number fluctuations. Finally, we show that a proposed mapping of the driven-dissipative many-body Bose-Hubbard model onto a single driven-dissipative Kerr model is not accurate for parameters in the hysteresis regime.

arXiv:2008.07316•Physical Review A•Nov 5, 2020•21 citations

Dark-state phase transition to a dark soliton in a dissipative Bose-Hubbard chain

Robbe Ceulemans, Samuel E. Begg, Matthew J. Davis +1 more

We show that a Bose-Hubbard model with local loss exhibits a first-order dark-state phase transition between a dark soliton and a uniform superfluid. Semiclassical simulations and stability analysis reveal that the soliton acts as a many-body dark state, whose suppressed fluctuations modify the transition relative to the comparable phenomenon of optical bistability in Kerr resonators. This mechanism quantitatively captures the bistable phase boundary observed in the experiment of Labouvie [], resolving long-standing discrepancies in previous theoretical modeling. Our results identify a general route by which spatial coherence and engineered loss can generate nonequilibrium phase transitions in ultracold gases.

arXiv:2509.25707•Physical Review Research•Jul 2, 2026•0 citations

Bloch oscillations in a Bose–Hubbard chain with single-particle losses

Bradley Longstaff, Eva-Maria Graefe

We theoretically investigate Bloch oscillations in a one-dimensional Bose–Hubbard chain, with single-particle losses from the odd lattice sites described by a Lindblad equation. For a single particle the time evolution of the state is completely determined by a non-Hermitian effective Hamiltonian. We analyse the spectral properties of this Hamiltonian for an infinite lattice and link features of the spectrum to observable dynamical effects, such as frequency doubling in breathing modes. We further consider the case of many particles in the mean-field limit leading to complex nonlinear Schrödinger dynamics. Analytic expressions are derived for the generalised nonlinear stationary states and the nonlinear Bloch bands. The interplay of nonlinearity and particle losses leads to peculiar features in the nonlinear Bloch bands, such as the vanishing of solutions and the formation of additional exceptional points. The stability of the stationary states is determined via the Bogoliubov–de Gennes equation and is shown to strongly influence the mean-field dynamics. Remarkably, even far from the mean-field limit, the stability of the nonlinear Bloch bands appears to affect the quantum dynamics. This is demonstrated numerically for a two-particle system.

arXiv:2002.12727•Journal of Physics B Atomic Molecular and Optical Physics•Jul 17, 2020•3 citations

Dynamic structure factor of a driven-dissipative Bose–Hubbard model

Subhanka Mal, Anushree Dey, Kingshuk Adhikary +1 more

Abstract Dynamic structure factor (DSF) is important for understanding excitations in many-body physics; it reveals information about the spectral and spatial correlations of fluctuations in quantum systems. Collective phenomena like quantum phase transitions of ultracold atoms are addressed by harnessing density fluctuations. Here, we calculate the DSF of a nonequilibrium spinless Bose–Hubbard model from the perspective of dissipative phase transition (DPT) in a steady state. Our methodology uses a homogeneous mean-field approximation to make the single-site hierarchy simpler and applies the Lindbladian perturbation method (LPM) to go beyond the single site, limited by the ratio of the inter-site hopping term to the Liouvillian gap as a small parameter. Our results show that the DSF near a DPT point is characteristically different from that away from the transition point, providing a clear density spectral signature of the DPT. In addition to comparing the two numerical frameworks, the mean-field results serve as a benchmark for proof-of-principle robustness of LPM. Despite the numerical difficulty, our methodology provides a computationally accessible route for studying density fluctuations in an open lattice quantum system without requiring large-scale computation.

arXiv:2506.09615•New Journal of Physics•Oct 27, 2025•1 citation

Critical slowing down in driven-dissipative Bose-Hubbard lattices

Filippo Vicentini, Fabrizio Minganti, R. Rota +2 more

We explore theoretically the dynamical properties of a first-order dissipative phase transition in coherently driven Bose-Hubbard systems, describing, e.g., lattices of coupled nonlinear optical cavities. Via stochastic trajectory calculations based on the truncated Wigner approximation, we investigate the dynamical behavior as a function of system size for one-dimensional (1D) and 2D square lattices in the regime where mean-field theory predicts nonlinear bistability. We show that a critical slowing down emerges for increasing number of sites in 2D square lattices, while it is absent in 1D arrays. We characterize the peculiar properties of the collective phases in the critical region.

arXiv:1709.04238•Physical Review A•Jan 31, 2018•120 citations
Motivated by recent experiments, we study the properties of large Bose-Hubbard chains with single-particle losses at one site using classical field methods. We construct and validate a compact effective model that reduces computations to only a few sites. We show that in the mean-field approach the description captures the stationary states of the dissipative mode very well. Not only is there good quantitative agreement in the hysteresis loop, but the dark soliton state can be reproduced as well due to the preservation of the U(1)\text{U}(1)U(1) symmetry. Bimodality of the steady states, observed on longer timescales, is studied using the truncated Wigner method. We compare the switching statistics and derive the effective Liouvillian gap as a function of the tunneling, showing that the effective description underestimates fluctuations.
We employ the perturbation series expansion for derivation of the reduced master equations for the three-site Bose-Hubbard model subject to strong atom losses from the central site. The model describes a condensate trapped in a triple-well potential subject to externally controlled removal of atoms. We find that the $\ensuremath{\pi}$-phase state of the coherent superposition between the side wells decays via two dissipation channels, the single-boson channel (similar to the externally applied dissipation) and the boson-pair channel. The quantum derivation is compared to the classical adiabatic elimination within the mean-field approximation. We find that the boson-pair dissipation channel is not captured by the mean-field model, whereas the single-boson channel is described by it. Moreover, there is a matching condition between the zero-point energy bias of the side wells and the nonlinear interaction parameter which separates the regions where either the single-boson or the boson-pair dissipation channel dominate. Our results indicate that the MMM-site Bose-Hubbard models, for M>2M>2M>2, subject to atom losses may require an analysis which goes beyond the usual mean-field approximation for correct description of their dissipative features. This is an important result in view of the recent experimental works on the single-site addressability of condensates trapped in optical lattices.
We study the dynamical properties of a driven-dissipative Bose-Hubbard model in the strongly interacting regime through a quantum trajectory approach with a cluster Gutzwiller Ansatz for the wave function. This allows us to take classical and quantum correlations into account. By studying the dynamical hysteresis surface that arises by sweeping through the coherent driving strength we show that the phase diagram for this system is in qualitative correspondence with the Gutzwiller mean-field result. However, quantitative differences are present and the inclusion of classical and quantum correlations causes a significant shift of the critical parameters. Additionally, we show that approximation techniques relying on a unimodal distribution such as the mean field and 1/z1/z1/z expansion drastically underestimate the particle number fluctuations. Finally, we show that a proposed mapping of the driven-dissipative many-body Bose-Hubbard model onto a single driven-dissipative Kerr model is not accurate for parameters in the hysteresis regime.
We show that a Bose-Hubbard model with local loss exhibits a first-order dark-state phase transition between a dark soliton and a uniform superfluid. Semiclassical simulations and stability analysis reveal that the soliton acts as a many-body dark state, whose suppressed fluctuations modify the transition relative to the comparable phenomenon of optical bistability in Kerr resonators. This mechanism quantitatively captures the bistable phase boundary observed in the experiment of Labouvie [], resolving long-standing discrepancies in previous theoretical modeling. Our results identify a general route by which spatial coherence and engineered loss can generate nonequilibrium phase transitions in ultracold gases.
We theoretically investigate Bloch oscillations in a one-dimensional Bose–Hubbard chain, with single-particle losses from the odd lattice sites described by a Lindblad equation. For a single particle the time evolution of the state is completely determined by a non-Hermitian effective Hamiltonian. We analyse the spectral properties of this Hamiltonian for an infinite lattice and link features of the spectrum to observable dynamical effects, such as frequency doubling in breathing modes. We further consider the case of many particles in the mean-field limit leading to complex nonlinear Schrödinger dynamics. Analytic expressions are derived for the generalised nonlinear stationary states and the nonlinear Bloch bands. The interplay of nonlinearity and particle losses leads to peculiar features in the nonlinear Bloch bands, such as the vanishing of solutions and the formation of additional exceptional points. The stability of the stationary states is determined via the Bogoliubov–de Gennes equation and is shown to strongly influence the mean-field dynamics. Remarkably, even far from the mean-field limit, the stability of the nonlinear Bloch bands appears to affect the quantum dynamics. This is demonstrated numerically for a two-particle system.
Abstract Dynamic structure factor (DSF) is important for understanding excitations in many-body physics; it reveals information about the spectral and spatial correlations of fluctuations in quantum systems. Collective phenomena like quantum phase transitions of ultracold atoms are addressed by harnessing density fluctuations. Here, we calculate the DSF of a nonequilibrium spinless Bose–Hubbard model from the perspective of dissipative phase transition (DPT) in a steady state. Our methodology uses a homogeneous mean-field approximation to make the single-site hierarchy simpler and applies the Lindbladian perturbation method (LPM) to go beyond the single site, limited by the ratio of the inter-site hopping term to the Liouvillian gap as a small parameter. Our results show that the DSF near a DPT point is characteristically different from that away from the transition point, providing a clear density spectral signature of the DPT. In addition to comparing the two numerical frameworks, the mean-field results serve as a benchmark for proof-of-principle robustness of LPM. Despite the numerical difficulty, our methodology provides a computationally accessible route for studying density fluctuations in an open lattice quantum system without requiring large-scale computation.
We explore theoretically the dynamical properties of a first-order dissipative phase transition in coherently driven Bose-Hubbard systems, describing, e.g., lattices of coupled nonlinear optical cavities. Via stochastic trajectory calculations based on the truncated Wigner approximation, we investigate the dynamical behavior as a function of system size for one-dimensional (1D) and 2D square lattices in the regime where mean-field theory predicts nonlinear bistability. We show that a critical slowing down emerges for increasing number of sites in 2D square lattices, while it is absent in 1D arrays. We characterize the peculiar properties of the collective phases in the critical region.