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 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.
This thesis investigates the physics of strongly interacting driven–dissipative bosonic systems in one dimension, combining field-theoretical approaches with various numerical simulation techniques. While weakly interacting condensates are accurately described by a generalized Gross–Pitaevskii equation, the strongly correlated regime of the same microscopic models remains far less understood, particularly out of equilibrium. This thesis aims to bridge this gap by developing controlled quantum-many-body simulations for open systems governed by Lindblad dynamics, as well as investigating the corresponding action under the functional renormalization group.In the first part, we formulate a field-theoretical description of driven–dissipative one-dimen-sional bosons using the Schwinger–Keldysh formalism and the functional renormalization group (fRG).Two models are analyzed: (i) a driven-dissipative condensate whose phase dynamics is expected to exhibit emergent Kardar–Parisi–Zhang (KPZ) behavior, for which we search the emergent KPZ fixed point starting from the bosonic field theory; and (ii) a two-photon–driven system with two-body losses realizing dissipative cat qubits. For the latter, we determine the phase diagram and identify a continuous phase transition between symmetry-broken and symmetric phases, which motivates the construction of an fRG Ansatz capturing this critical behavior.The second part of the thesis focuses on the open system dynamics using exact diagonalization and tensor-network methods. For the Bose–Hubbard model with incoherent pump and loss, we characterize the steady state, and study spectral functions across the full parameter space. We show that key ground-state features are absent under dissipation, while the high-temperature unitary spectrum captures aspects of the open-system dynamics. In an extended model with additional two-body and non-local losses, we establish the semi-classical emergence of KPZ scaling and extract the KPZ dynamical exponent from the spectral function. A finite-size analysis exposes the current limitations of the tensor-network methods for reaching the high occupancy regime in the full quantum evolution.Finally, a last project develops a two-loop perturbative RG analysis of a spin-1 chain at its continuous symmetry breaking quantum critical point, revealing critical exponents distinct from the Ising universality class.Together, these results advance the understanding of strongly interacting open quantum systems, elucidating their symmetries, effective field theories, and dynamical properties, and outlining key directions for future theoretical and numerical work.
We consider a one-dimensional driven-dissipative Bose-Hubbard model subjected to incoherent pump and one- and two-body losses, and analyze it by studying its two-point space-time correlations. By employing a combination of numerical methods, such as stochastic semi-classical simulations and tensor network methods, as well as perturbation theory within the Keldysh formalism, we characterize the system at varying filling and interaction strength. We present results for two complementary regimes: i) at large filling and weak interactions, where we show that Kardar-Parisi-Zhang scaling is visible in the linewidth of the spectral function and ii) at weak filling and strong interactions, where we analyze what remains of the mean field transition and identify a change of nature of the excitations. Our study covers two important regions of the phase diagram of such a system.
We study the non-equilibrium behavior of optically driven dissipative coupled resonator arrays. Assuming each resonator is coupled with a two-level system via a Jaynes–Cummings interaction, we calculate the many-body steady state behavior of the system under coherent pumping and dissipation. We propose and analyze the many-body phases using experimentally accessible quantities such as the total excitation number, the emitted photon spectra and photon coherence functions for different parameter regimes. In parallel, we also compare and contrast the expected behavior of this system assuming the local nonlinearity in the cavities is generated by a generic Kerr effect as described by the Bose–Hubbard (BH) model rather than a Jaynes–Cummings interaction. We find that the behavior of the experimentally accessible observables produced by the two models differs for realistic regimes of interactions even when the corresponding nonlinearities are of similar strength. We analyze in detail the extra features available in the Jaynes–Cummings–Hubbard (JCH) model originating from the mixed nature of the excitations and investigate the regimes where the BH approximation would faithfully match the JCH physics. We find that the latter is true for values of the light–matter coupling and losses beyond the reach of current technology. Throughout the study we operate in the weak pumping, fully quantum mechanical regime where approaches such as mean field theory fail, and instead use a combination of quantum trajectories and the time evolving block decimation algorithm to compute the relevant steady state observables. In our study we have assumed small to medium size arrays (from 3 up to 16 sites) and values of the ratio of coupling to dissipation rate g / γ ∼ 20, which makes our results implementable with current designs in circuit QED and with near future photonic crystal set ups.
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.
Abstract We investigate signatures of a self-trapping transition in the driven-dissipative Bose Hubbard dimer, in presence of incoherent pump and single-particle losses. For fully symmetric couplings the stationary state density matrix is independent of any Hamiltonian parameter, and cannot therefore capture the competition between hopping-induced delocalization and the interaction-dominated self-trapping regime. We focus instead on the exact quantum dynamics of the particle imbalance after the system is prepared in a variety of initial states, and on the frequency-resolved spectral properties of the steady state, as encoded in the single-particle Green’s functions. We find clear signatures of a localization-delocalization crossover as a function of hopping to interaction ratio. We further show that a finite a pump-loss asymmetry restores a delocalization crossover in the steady-state imbalance and leads to a finite intra-dimer dissipation.
Driven-dissipative light-matter systems can exhibit collective nonequilibrium phenomena due to loss and gain processes on the one hand and effective photon-photon interactions on the other hand. As a generic example we study a bosonic lattice system implemented via an array of driven-dissipative coupled nonlinear resonators with linearly increasing resonance frequencies across the lattice. The model also describes a driven-dissipative Bose-Hubbard model in a tilted potential without a particle-conservation constraint. We numerically predict a diverse range of stationary and nonstationary states resulting from the interplay of the tilt, tunneling, on-site interactions, and loss and gain processes. Our key finding is that, under weak on-site interactions, the bosons mostly condense into a selected, single-particle Wannier-Stark state without exhibiting the expected Bloch oscillations. As the strength of the on-site interactions increase, a nonstationary regime emerges which, surprisingly, exhibits periodic Bloch-type oscillations. As a direct consequence of the driven-dissipative nature of the system we predict a highly nontrivial phase diagram including regular oscillating as well as chaotic dynamical regimes. While a straightforward photonic implementation using microwave or optical modes is possible, such dynamics might also be observable for an ultracold gas in a vertical lattice with gravity or a tilted external potential.