Driven-dissipative Bose-Hubbard Model Research Articles

Steady-state quantum Zeno effect of driven-dissipative bosons with dynamical mean-field theory

We study a driven-dissipative Bose-Hubbard model in presence of two-particle losses and an incoherent single-particle drive on each lattice site, leading to a finite-density stationary state. Using dynamical mean-field theory (DMFT) and an impurity solver based on exact diagonalization of the associated Lindbladian, we investigate the regime of strong two-particle losses. Here, a stationary-state quantum Zeno effect emerges, as can be seen in the on-site occupation and spectral function. We show that DMFT captures this effect through its self-consistent bath. We show that, in the deep Zeno regime, the bath structure simplifies, with the occupation of all bath sites except one becoming exponentially suppressed. As a result, an effective dissipative hard-core Bose-Hubbard dimer model emerges, where the auxiliary bath site has single-particle dissipation controlled by the Zeno dissipative scale.

Discrete time crystal in a driven-dissipative Bose-Hubbard model with two-photon processes

We propose a scheme for discrete time crystals (DTCs) based on the metastability of the driven-dissipative Bose-Hubbard model, linking the two metastable states as such using rotation operations and two-photon processes. Numerical analyses in finite size, through the Gutzwiller Monte Carlo approach, validate our scheme, where the couplings between different sites are taken into account. Requiring only local dissipations and uniform rotations, our scheme is feasible on the state-of-the-art noisy intermediate-scale quantum platforms, e.g., superconducting circuits and ultracold atoms. The robustness of the proposed DTC is guaranteed by construction due to the metastability. Our scheme provides a simple test bed for studying the less explored field of nonequilibrium phases of matter in physics and exploring more possibilities of quantum simulations.

Photon Transport in a Bose-Hubbard Chain of Superconducting Artificial Atoms.

We demonstrate nonequilibrium steady-state photon transport through a chain of five coupled artificial atoms simulating the driven-dissipative Bose-Hubbard model. Using transmission spectroscopy, we show that the system retains many-particle coherence despite being coupled strongly to two open spaces. We find that cross-Kerr interaction between system states allows high-contrast spectroscopic visualization of the emergent energy bands. For vanishing disorder, we observe the transition of the system from the linear to nonlinear regime of photon blockade in excellent agreement with the input-output theory. Finally, we show how controllable disorder introduced to the system suppresses nonlocal photon transmission. We argue that proposed architecture may be applied to analog simulation of many-body Floquet dynamics with even larger arrays of artificial atoms paving an alternative way towards quantum supremacy.

Time- and frequency-domain two-particle correlations of a driven dissipative Bose-Hubbard model

We investigate theoretically the time- and frequency-domain two-particle correlations of a driven-dissipative Bose-Hubbard model at and near a dissipative phase transition (DPT). We compute the Hanbury Brown--Twiss (HBT)-type two-particle temporal correlation function ${g}^{2}(\ensuremath{\tau})$ which, as a function of time delay $\ensuremath{\tau}$, exhibits oscillations with frequencies determined by the imaginary part of the Liouvillian gap. As the gap closes near a transition point, the oscillations at that point die down. For parameters slightly away from the transition point, the HBT correlations show oscillations from superbunching to antibunching regimes. We show that the Fourier transform of HBT correlations into the frequency domain provide information about DPT and Liouvillian dynamics. We numerically solve the many-body Lindblad master equation and calculate the Wigner distribution of the system in the steady state to ascertain the DPT. Below a certain drive strength, the Fourier transform shows a two-peak structure, while above that strength it exhibits either a Lorentzian-like single-peak structure or a structure with two dips. The width of the single-peak structure is minimal at the phase-transition point and the peak of this structure always lies at zero frequency. The positions of the two symmetrical peaks in case of a two-peak structure are given by the imaginary parts of the Liouvillian gap while their half width at half maximum is given by the real part of the gap. The positions and widths of the two dips are also related to low-lying eigenvalues of the Liouvillian operator. We discuss quantum statistical properties of the model in terms of the HBT correlation function and its Fourier transform.

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

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.

Pattern Formation and Exotic Order in Driven-Dissipative Bose-Hubbard Systems.

Modern experimental platforms such as superconducting circuit arrays call for the exploration of bosonic tight-binding models in unconventional situations with no counterpart in real materials. Here we investigate one such situation in which excitations are driven and damped by pairs, leading to pattern formation and exotic bosonic states emerging from a nonequilibrium quantum many-body system. Focusing on a two-dimensional driven-dissipative Bose-Hubbard model, we find that its steady states are characterized by the condensation of bosons around momenta lying on a "Bose surface," a bosonic analog of the Fermi surface in solid-state systems. The interplay between instabilities generated by the driving, the nonlinear dissipative mode coupling, and the underlying lattice effect allows the system to equilibrate into an exotic superfluid state of bosons condensed on a closed ring in momentum space instead of discrete points. Such an unconventional state with a spatially uniform density distribution goes beyond the traditional scope of pattern formation and thus has no counterpart in the classical literature. In addition, it is a state connected to several open problems in modern condensed-matter physics. Here we provide the means to stabilize it, opening the way to its experimental study. Moreover, we also provide a concrete experimental implementation of our model in currently available superconducting circuit arrays. We also investigate the relaxation spectrum around the condensate, which shows a characteristic purely diffusive behavior.

Excitations and correlations in the driven-dissipative Bose-Hubbard model

Using a field-theoretic approach within the Schwinger-Keldysh formalism, we study a Bose-Hubbard model in the presence of a driving field and dissipation due to one-body losses. We recover the bistability diagram from the Gross-Pitaevski equation and analyze the different phases with respect to their elementary excitations and correlations. We find the low-density solution to be subdivided into a dynamically instable, a gapped, and a gapless regime. The correlations decay exponentially, but a substantial increase of correlation length marks the regime of gapless excitations.

Optimal stochastic unraveling of disordered open quantum systems: Application to driven-dissipative photonic lattices

We propose an efficient numerical method to compute configuration averages of observables in disordered open quantum systems whose dynamics can be unraveled via stochastic trajectories. We prove that the optimal sampling of trajectories and disorder configurations is simply achieved by considering one random disorder configuration for each individual trajectory. As a first application, we exploit the present method to the study the role of disorder on the physics of the driven-dissipative Bose-Hubbard model in two different regimes: (i) for strong interactions, we explore the dissipative physics of fermionized bosons in disordered one-dimensional chains; (ii) for weak interactions, we investigate the role of on-site inhomogeneities on a first-order dissipative phase transition in a two-dimensional square lattice.

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

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.

Spatial correlations in driven-dissipative photonic lattices

We study the nonequilibrium steady-state of interacting photons in cavity arrays as described by the driven-dissipative Bose–Hubbard and spin-1/2 XY model. For this purpose, we develop a self-consistent expansion in the inverse coordination number of the array () to solve the Lindblad master equation of these systems beyond the mean-field approximation. Our formalism is compared and benchmarked with exact numerical methods for small systems based on an exact diagonalization of the Liouvillian and a recently developed corner-space renormalization technique. We then apply this method to obtain insights beyond mean-field in two particular settings: (i) we show that the gas–liquid transition in the driven-dissipative Bose–Hubbard model is characterized by large density fluctuations and bunched photon statistics. (ii) We study the antibunching–bunching transition of the nearest-neighbor correlator in the driven-dissipative spin-1/2 XY model and provide a simple explanation of this phenomenon.

Spontaneous Beliaev-Landau scattering out of equilibrium

We investigate Beliaev-Landau scattering in a gas of interacting photons in a coherently driven array of nonlinear dissipative resonators, as described by the 1D driven-dissipative Bose-Hubbard model. Due to the absence of detailed balance in such an out-of-equilibrium setup, steady-state properties can be much more sensitive to the underlying microscopic dynamics. As the popular truncated Wigner approximation dramatically fails in capturing this physics, we present an alternative approach, based on a systematic expansion beyond the Bogoliubov approximation, which includes the third-order correlation functions in the dynamics. As experimentally accessible signatures of Beliaev-Landau processes, we report a small but nonnegligible correction to the Bogoliubov prediction for the steady-state momentum distribution showing a characteristic series of peaks and dips, as well as non-Gaussian features in the statistics of the cavity output field.

Nonequilibrium gas-liquid transition in the driven-dissipative photonic lattice

We study the nonequilibrium steady state of the driven-dissipative Bose-Hubbard model with Kerr nonlinearity. Employing a mean-field decoupling for the intercavity hopping $J$, we find that the steep crossover between low and high photon-density states inherited from the single cavity transforms into a gas$-$liquid bistability at large cavity-coupling $J$. We formulate a van der Waals like gas$-$liquid phenomenology for this nonequilibrium situation and determine the relevant phase diagrams, including a new type of diagram where a lobe-shaped boundary separates smooth crossovers from sharp, hysteretic transitions. Calculating quantum trajectories for a one-dimensional system, we provide insights into the microscopic origin of the bistability.

On the robustness of strongly correlated multi-photon states in frustrated driven-dissipative cavity lattices

We present a theoretical study on the robustness of multi-photon states in a frustrated lattice of coupled nonlinear optical cavities, which are described by a driven-dissipative Bose-Hubbard model. In particular, we focus here on a Lieb lattice with two elementary cells and periodic boundary conditions. Due to the geometric frustration of the lattice, the non-equilibrium steady state can exhibit dark sites with low photon density and strong correlations, ascribable to the population of multi-photon modes. We explore the sensitivity of such strong correlations on the random inhomogeneity of the lattice parameters. We show that the correlations are more sensitive to the inhomogeneity of the cavity frequencies than to the random fluctuations of the hopping strength.

Emergent equilibrium in many-body optical bistability.

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.

Truncated correlation hierarchy schemes for driven-dissipative multimode quantum systems

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.

Bose-Hubbard model: Relation between driven-dissipative steady states and equilibrium quantum phases

We present analytical solutions for the mean-field master equation of the driven-dissipative Bose-Hubbard model for cavity photons, in the limit of both weak pumping and weak dissipation. Instead of pure Mott insulator states, we find statistical mixtures with the same second-order coherence as a Fock state with n photons, but a mean photon number of n/2. These mixed states occur when n pump photons have the same energy as n interacting photons inside the nonlinear cavity and survive up to a critical tunneling coupling strength, above which a crossover to classical coherent state takes place. We also explain the origin of both antibunching and superbunching predicted by P-representation mean-field theory at higher pumping and dissipation. In particular, we show that the strongly correlated region of the associated phase diagram cannot be described within the semiclassical Gross-Pitaevski approach.

Steady-state phases and tunneling-induced instabilities in the driven dissipative Bose-Hubbard model.

We determine the steady-state phases of a driven-dissipative Bose-Hubbard model, describing, e.g., an array of coherently pumped nonlinear cavities with a finite photon lifetime. Within a mean-field master equation approach using exact quantum solutions for the one-site problem, we show that the system exhibits a tunneling-induced transition between monostable and bistable phases. We characterize the corresponding quantum correlations, highlighting the essential differences with respect to the equilibrium case. We also find collective excitations with a flat energy-momentum dispersion over the entire Brillouin zone that trigger modulational instabilities at specific wave vectors.