We adopt a diffusive-Gutzwiller approach to investigate a phase transition in a quadratically driven-dissipative Bose-Hubbard lattice. Diffusive trajectories may lead to lower average entanglement as compared to jump-like trajectories and have a natural tendency to approach coherent states, therefore the method can be less prone to the bias induced by the fully uncorrelated Gutzwiller ansatz. Averaging over trajectories does lead to classical correlations and this allows us to address the correlation length of such 2D lattices of open quantum systems which is the main goal of this work. Under this approximations, we find negligible correlation length in the low density phase and apparently unbounded length grow in the high density phase. Additionally, we show that the effective relaxation times associated to the times scales for synchronisation in the high density phase may also diverge suggesting the vanishing of the Lindbladian gap.
We adopt a diffusive-Gutzwiller approach to investigate a phase transition in a quadratically driven-dissipative Bose-Hubbard lattice. Diffusive trajectories may lead to lower average entanglement as compared to jump-like trajectories and have a natural tendency to approach coherent states, therefore the method can be less prone to the bias induced by the fully uncorrelated Gutzwiller ansatz. Averaging over trajectories does lead to classical correlations and this allows us to address the correlation length of such 2D lattices of open quantum systems which is the main goal of this work. Under this approximations, we find negligible correlation length in the low density phase and apparently unbounded length grow in the high density phase. Additionally, we show that the effective relaxation times associated to the times scales for synchronisation in the high density phase may also diverge suggesting the vanishing of the Lindbladian gap.
We examine theoretically the driven-dissipative Bose-Hubbard dimer in the optical bistable regime. Various approximation schemes based on a Gutzwiller mean-field decoupling are applied and compared. Depending on the system parameters, we show that a decoupling with respect to the real space or to the reciprocal space can be more accurate. The Gutzwiller decoupling is applied both at the level of the density matrix and for the wave function during a quantum trajectory simulation. The latter is shown to be a more accurate approximation. A Gaussian approximation for the nonhomogeneous antibonding mode is also explored. We also show that entanglement in this system is witnessed by squeezing in reciprocal space.
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.
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 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.
The present Thesis aims to deepen the understanding of the physics of strongly-correlated systems along diverse directions with a common denominator, namely a focus on the physical role played by the collective modes shaping quantum many-body correlations. This investigation has been mostly led through the heuristic application of the Quantum Gutzwiller Approach, a recently developed theory of the quantised excitations above the Gutzwiller approximation for strongly-interacting systems on a lattice. In the first place, the intertwining between quantum correlations and collective effects is explored in the context of Bose-Hubbard models, intended both as simple yet highly non-trivial realisations of complex many-body scenarios, and as controllable environments for the quantum simulation of impurity problems. In particular, we are able to provide a comprehensive and semi-analytical description of quantum fluctuations across the critical regimes induced by the proximity of a Mott transition (Mottness), ranging from accurate predictions for the superfluid density and pair correlations across quantum critical regimes to the Andreev-Bashkin effect in a two-component bosonic mixture. As case studies for the effectiveness and flexibility of our approach, we give a detailed account of the rich behaviours of a quantum impurity, either fixed or itinerant, immersed in a Bose-Hubbard environment, showing that both the decoherence dynamics of spinful particle and the resulting Bose polaron in the second case are extremely sensitive probes of the type of critical fluctuations experienced by the surrounding system. Additionally, we extend our analysis of quantum fluctuations to the Fermi-Hubbard model, where we uncover the interplay between low-energy quasiparticles and bosonic elementary excitations in determining two-particle correlations. A more complex physical scenario addressed in this Thesis is represented by driven-dissipative interacting systems, which find a state-of-the-art technological implementation in circuit QED platforms, among others. Inspired by the experimental interest in these systems as a potential toolbox for the generation of exotic quantum states of light, we consider a strongly-correlated fluid of photons on a lattice stabilised via incoherent drive and dissipation, and examine its collective fluctuations across the out-of-equilibrium insulator-to-superfluid transition of the system. Specifically, establishing a conceptual parallelism with our approach to fluctuations at equilibrium, we develop a linear-response theory of the excitations on top of the Gutzwiller stationary state of the Lindblad dynamics, and characterise their physical role in quantum observables and response functions. Taking advantage of our formalism, we describe also dynamical instabilities and novel physical regimes arising in the ultrastrong coupling limit of the superfluid phase. Among the most remarkable results, we bring to light the paradoxical nature of the Mott-like insulating phase of the system, whose peculiar one-body dynamical fluctuations are shown to echo the behaviour of a lasing state. Encompassing diverse modern domains of many-body physics, our investigation goes in the direction of solidifying the interpretation of collective modes as fundamental actors in correlation-induced phenomena under a unified physical picture.
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.
In recent years, the control of photon-photon interactions in optical nonlinear media has led to the realization of quantum fluids of light. One of the current challenges is to increase the strength of these interactions and enter the so-called strongly correlated regime. To achieve this goal, arrays of nonlinear cavities are a very promising candidate. In this thesis, theoretical results on arrays of nonlinear cavities described by a driven-dissipative Bose-Hubbard model are presented. In particular, a general method to compute the mean-field phase diagram of this model is described. Due to the finite life time of photons, the system is intrinsically dissipative: cavity losses must be compensated by an external driving field. This nonequilibrium nature gives rise to interesting features, such as a transition between monostable and bistable phases induced by tunneling. In the limit of weak dissipation and weak driving, analytical results describing generalized Mott insulating phases are derived. These states survive up to a critical tunneling strength, above which a crossover to a classical coherent state takes place. Finally, the issue of how to go beyond the mean-field approximation is addressed by performing exact numerical simulations. Large arrays of cavities were simulated by implementing a new method specifically tailored for driven-dissipative systems.