We investigate the influence of a Markovian environment on the dynamics of interacting spinful fermionic atoms in a lattice. To explore the physical phenomena occurring at short times, we develop a method based on a slave-spin representation of fermions that is amenable to the investigation of the dynamics of dissipative systems. We apply this approach to two different dissipative couplings that can occur in current experiments: a coupling via the local density and a coupling via the local double occupancy. We complement our study based on this method, with results obtained using the adiabatic elimination technique and with an exact study of a two-site model. We uncover that the decoherence is slowed down by increasing either the interaction strength or the dissipative coupling (the Zeno effect). We also find, for the coupling to the local double occupancy, that the final steady state can sustain single-particle coherence.
We investigate the influence of a Markovian environment on the dynamics of interacting spinful fermionic atoms in a lattice. To explore the physical phenomena occurring at short times, we develop a method based on a slave-spin representation of fermions that is amenable to the investigation of the dynamics of dissipative systems. We apply this approach to two different dissipative couplings that can occur in current experiments: a coupling via the local density and a coupling via the local double occupancy. We complement our study based on this method, with results obtained using the adiabatic elimination technique and with an exact study of a two-site model. We uncover that the decoherence is slowed down by increasing either the interaction strength or the dissipative coupling (the Zeno effect). We also find, for the coupling to the local double occupancy, that the final steady state can sustain single-particle coherence.
Remarkable experimental progress and theoretical advancements in the study of strongly corre- lated systems have opened up new doors towards a better understanding of interesting states of matter. Such states could possibly be used to develop novel technological applications based on the principles of quantum mechanics. In recent years, the investigation of both complex solid state and ultracold atom systems has lead to a number of fascinating discoveries with regards to their non-equilibrium properties. One route to probe non-equilibrium physics is to consider the evolution of open systems defined by the coupling to an external environment introducing dissipation. Engineering the environment according to predefined requirements can then be used to dynamically prepare states which are difficult to access otherwise. This thesis explores the dynamics of many-body open quantum systems aiming at under- standing the effects resulting from the interplay of kinetic motion, interactions and dissipation in different physical settings. First, we build upon recent developments concerning the coupling of ultracold atoms to optical resonators which are subject to photon losses. We show how the dissipation can trigger an attractor dynamics towards a topologically non-trivial steady state that is robust against perturbations. The numerically obtained full open system evolution en- ables us to not only characterize in detail the steady state but also to analyze the features of the approaching dynamics. Second, we study the propagation of two-time correlation functions in interacting quantum spin-1/2 chains in contact with an environment causing dephasing. Pro- vided with the numerically exact full time-evolution, we identify an algebraic scaling regime where time-translation invariance breaks down and uncover the emergence of aging dynam- ics. Finally, we study a system of interacting fermionic atoms in a one-dimensional lattice with particle loss at the central site. There we report on the dissipation-induced deceleration of the non-equilibrium dynamics known as the quantum Zeno effect. As in the majority of cases quantum many-body systems are not solvable analytically, another central aspect of this thesis deals with adapting numerical methods based on matrix product state algorithms to investigate open quantum systems dynamics.
We present two approaches capable of describing the dynamics of an interacting many-body system on a lattice coupled globally to a dissipative bosonic mode.Physical realizations are, for example, ultracold atom gases in optical lattice coupled to a photonic mode of an optical cavity or electronic gases in solids coupled to THz cavity fields.The first approach, applicable for large dissipation strengths and any system size, is a variant of the many-body adiabatic elimination method for investigating the long-time dynamics of the system.The second method extends the time-dependent matrix product techniques to capture the global coupling of the interacting particles to the bosonic mode and its open nature.It gives numerically exact results for small to intermediate system sizes.As a benchmark for our methods we perform the full quantum evolution of a Bose-Hubbard chain coupled to a cavity mode.We show that important deviations from the mean-field behavior occur when considering the full atoms cavity coupling [C.
We develop a discrete truncated Wigner method to analyze the real-time evolution of dissipative SU(N) spin systems coupled with a Markovian environment. This semiclassical approach is not only numerically efficient but also particularly capable of accurately capturing local loss processes due to its local linearity in the dynamical equations. We apply the method to a state-of-the-art experiment involving an analog quantum simulator of a three-dimensional dissipative Bose-Hubbard model in a strongly interacting regime. Our numerical results show good agreement with experimental data, specifically capturing the continuous quantum Zeno effect in the dynamics subjected to a gradual change of the ratio between the hopping amplitude and the onsite interaction across the superfluid Mott insulator crossover. Furthermore, we present comparative analyses with the continuous truncated Wigner method, derived as an effective Fokker-Planck equation for SU(N) classical spin variables, showing that the discrete method outperforms the continuous one in simulating the long-time dynamics of SU(2) and SU(3) spin models. The discrete phase-space framework offers a versatile and powerful tool for exploring a wide range of open quantum many-body systems in dimensions higher than one dimension, where numerically exact methods are impractical due to the exponential growth of the Hilbert-space dimension.
This thesis is about the theory of open quantum many-body physics with a particular focus on driven-dissipative photonic and spin lattices. After a review of the main physical platforms and theoretical concepts, we describe our original results. In Chapter 2, we present a theory of the single-mode Kerr model with a time-dependent pump. This model describes a single-mode optical cavity with a third-order optical nonlinearity. In the regime of parameters where the semiclassical analysis shows bistability, we find that a dynamic hysteresis loop appears in the exact solution. The hysteresis area as a function of the sweep time shows a double power-law decay where the second exponent is independent of the system parameters. We show how these effects are related to the emergence of a dissipative phase transition. We also describe concisely the experimental results which have recently confirmed such theoretical predictions. In Chapter 3, we present the corner-space renormalization method. In order to obtain the steady-state density matrix of the lattice, we solve the Lindblad master equation in a subspace of the Hilbert space (the "corner"). The states spanning the corner space are selected iteratively using eigenvectors of the density matrix of smaller lattice systems, merging in real space two lattices at each iteration and selecting M pairs of states by maximizing their joint probability. The accuracy of the results is then improved by increasing the dimension M of the corner space until convergence is reached. The method has been benchmarked on a two-dimensional Bose-Hubbard model with coherent driving. The strength and limitations of the method are critically discussed.In Chapter 4, we investigate a dissipative phase transition in the two-dimensional anisotropic Heisenberg XYZ model. Using the corner-space renormalization method we present a finite-size analysis of steady-state observables. In particular, we show the critical behaviour of the magnetic susceptibility, the entropy growth and the entanglement witnesses, providing a first evaluation of the critical exponents characterizing the transition. A study of the dynamics of finite-size systems is also consistent with a critical slowing down. For comparison, we present the corresponding analysis for one-dimensional arrays, showing the absence of criticality due to the reduced dimension.In Chapter 5, we explore the non-equilibrium photonic phases of a dissipative Bose-Hubbard model with incoherent pumping of coupled two-level systems. Within a Gutzwiller mean-field approach, we determine the steady-state phase diagram of the system. We predict a second-order phase transition between an incompressible Mott-like phase and a coherent delocalized phase.
We study the out-of-equilibrium dynamics of a dilute, lattice-confined Bose-Fermi mixture initialized in a highly excited state consisting of boson-fermion pairs (doublons) occupying single lattice sites. This system represents a paradigmatic case for studying relaxation dynamics in strongly correlated systems, and provides a versatile platform for studying thermalization and localization phenomena. We provide analytical expressions for the short-time decay of isolated doublons and small doublon clusters due to the competition between tunneling and interparticle interactions. We also discuss a mechanism for long-time decay that crucially depends on the quantum statistics of the particles constituting the doublon, namely, the conversion of pairs of neighboring doublons into an unpaired fermion and a site with a fermion and two bosons. Building on these insights, we develop a cluster expansion method to describe the dynamics in extended systems and compare it to numerically exact matrix product state simulations in one dimension. Finally, we discuss how our predictions can be observed in experiments with ultracold heteronuclear molecules.
All systems are open to an essentially uncontrollable environment that acts as a source of decoherence and dissipation. In some cases the environment's only effect is to add a weak relaxation mechanism and thus can be ignored for short timescales. In others, however, the presence of the environment can fundamentally alter the behavior of the system. Such is the case in mesoscopic superconductors where the environment can stabilize superconductivity and in spin-boson systems where the environment induces a localization transition. Likewise, in technological applications we are often interested in systems operating far from equilibrium. Here the environment might act as a particle reservoir or strong driving force. In all these examples, we need accurate methods to describe the influence of the environment on the system and to solve for the resulting dynamics or equilibrium states. In this thesis, we develop computational and conceptual approaches to efficiently simulate quantum systems in contact with an environment. Our starting point is the use of numerical renormalization techniques. Thus, we restrict our attention to one-dimensional lattices or small quantum systems coupled to an environment. We have developed several complementary algorithms: a superoperator renormalization algorithm for simulating real-time Markovian dynamics and for calculating states in thermal equilibrium; a blocking algorithm for simulating integro-differential equations with long-time memory; and a tensor network algorithm for branched lattices, which can be used to simulate strongly dissipative systems. Further, we provide support for an idea that to generically and accurately simulate the real-time dynamics of strongly dissipative systems, one has to include all or part of the environment within the simulation. In addition, we discuss applications and open questions.