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.
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 g2(\ensuremathτ) which, as a function of time delay \ensuremathτ, 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.
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.
Dissipative phase transitions (DPTs) are traditionally characterized through the spectrum of a time-independent Liouvillian superoperator. However, this definition does not apply to time-periodic (Floquet) systems that cannot be exactly recast as time-independent problems. Here, we develop a general framework to characterize DPTs in time-periodic open quantum systems through the spectrum of the Floquet propagator. We first study driven-dissipative Kerr resonators, known to display a DPT, showing that counter-rotating terms in the drive shift the critical point and significantly change the time scales associated with the transition. We then investigate DPTs in the driven quantum Rabi model and its time-independent approximation, the driven Jaynes-Cummings model, finding that the Rabi model exhibits distinct critical features as the ultrastrong coupling regime is approached. Moreover, our Floquet analysis unveils the disappearance of the DPT in the deep strong coupling regime, due to light-matter decoupling. Our approach sets the stage for the study of dissipative criticality in a broad class of time-dependent open quantum systems.
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 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.
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 study the self-trapping transition in dissipative double wells using mean-field theory and exact diagonalization. We find that, similar to the conventional interaction-driven transition from oscillatory dynamics to self-trapping, dissipation alone can also induce self-trapping-like localization in one well. The resulting state, although transient and decaying, persists for a substantial time. At the mean-field level, the interaction- and dissipation-driven transitions are discontinuous and continuous, respectively. When interaction and dissipation both are present, their interplay produces three dynamical regimes: (i) at weak interaction, increasing dissipation drives a transition from oscillations to self-trapping; (ii) at moderate interaction, self-trapping at weak dissipation, oscillatory dynamics at intermediate dissipation, and re-entrant self-trapping at strong dissipation; and (iii) at strong interaction, self-trapping persists but is weakened by dissipation. We connect the transition with the Liouvillian spectrum. It is observed that purely dissipation-driven transition coincides with an exceptional point in the spectra, and a peak in the Liouvillian gap, marking the onset of the quantum-Zeno regime. In the presence of interactions, however, the onset of the Zeno regime is shifted to dissipation strengths larger than that required for the self-trapping transition.