Capturing non-Markovian dynamics of open quantum systems is generally a challenging problem, especially for strongly interacting many-body systems. In this Letter, we combine recently developed non-Markovian quantum state diffusion techniques with tensor network methods to address this challenge. As a first example, we explore a Hubbard-Holstein model with dissipative phonon modes, where this new approach allows us to quantitatively assess how correlations spread in the presence of non-Markovian dissipation in a 1D many-body system. We find regimes where correlation growth can be enhanced by these effects, offering new routes for dissipatively enhancing transport and correlation spreading, relevant for both solid state and cold atom experiments.
Capturing non-Markovian dynamics of open quantum systems is generally a challenging problem, especially for strongly interacting many-body systems. In this Letter, we combine recently developed non-Markovian quantum state diffusion techniques with tensor network methods to address this challenge. As a first example, we explore a Hubbard-Holstein model with dissipative phonon modes, where this new approach allows us to quantitatively assess how correlations spread in the presence of non-Markovian dissipation in a 1D many-body system. We find regimes where correlation growth can be enhanced by these effects, offering new routes for dissipatively enhancing transport and correlation spreading, relevant for both solid state and cold atom experiments.
We introduce a general numerical method to compute the dynamics and multitime correlations of chains of quantum systems, where each system may couple strongly to a structured environment. The method combines the process tensor formalism for general (possibly non-Markovian) open quantum systems with time-evolving block decimation for one-dimensional chains. It systematically reduces the numerical complexity originating from system-environment correlations before integrating them into the full many-body problem, making a wide range of applications numerically feasible. We illustrate the power of this method by studying two examples. First, we study the thermalization of individual spins of a short XYZ Heisenberg chain with strongly coupled thermal leads. Our results confirm the complete thermalization of the chain when coupled to a single bath, and they reveal distinct effective temperatures in low-, mid-, and high-frequency regimes when the chain is placed between a hot and a cold bath. Second, we study the dynamics of diffusion in a longer XY chain, when each site couples to its own bath.
Recent advances in numerical methods significantly pushed forward the understanding of electrons coupled to quantized lattice vibrations. At this stage, it becomes increasingly important to also account for the effects of physically inevitable environments. Here, we combine state-of-the-art tensor-network and quantum trajectories methods in order to study the impact of dissipation on realistic condensed matter models including highly excited phononic modes. In particular, we study the transport properties of the Hubbard-Holstein Hamiltonian that models a large class of materials characterized by strong electron-phonon coupling, in contact with a dissipative environment. We combine the non-Markovian hierarchy of pure states method and the Markovian quantum jumps method with the newly introduced projected purified density-matrix renormalization group, creating powerful tensor-network methods for dissipative quantum many-body systems. Investigating their numerical properties, we find a significant speedup up to a factor $\ensuremath{\approx}30$ compared to conventional tensor-network techniques. We apply these methods to study dissipative quenches, aiming for an in-depth understanding of the formation, stability, and quasiparticle properties of bipolarons. Surprisingly, our results show that in the metallic phase dissipation localizes the bipolarons, which is reminiscent of an indirect quantum Zeno effect. However, the bipolaronic binding energy remains mainly unaffected, even in the presence of strong dissipation, exhibiting remarkable bipolaron stability. These findings shed light on the problem of designing real materials exhibiting phonon-mediated high-TC superconductivity.
The paradigm of considering open quantum systems—i.e., focusing only on the system of interest, and treating the rest of the world as an effective environment—has proven to be a highly effective way to understand a range of quantum systems, across areas of study such as quantum optics, cold atoms, superconducting qubits, and impurities in solid-state systems. A common approach in many of these contexts has been to consider simplified approaches based on the Born and Markov approximations. While these approximations are indeed often appropriate in contexts such as quantum optics, the widespread application of these approximations has been driven more by simplicity than by accuracy. In particular, these Markovian treatments will fail in many cases, such as when coupling to the environment is not weak, when the environment is structured and has resonances, when the system couples to low-frequency modes of the environment, or when the questions of interest involve the propagation of information through the environment. Despite the fact that many real problems are non-Markovian, the Markov approximation is still widely used, as it is often assumed that a fully non-Markovian treatment is too complex to be practical. In this Perspective we discuss a recently developed set of techniques that address this challenge. Centering our discussion around the notion of the process tensor, we demonstrate that the generality of the process tensor concept, coupled with efficient tensor-network methods, opens the door to the description of a wide range of observable non-Markovian processes in a wide range of open quantum systems.
The nonequilibrium effects of dissipation and drive play a key role in an immense variety of nanoscale and mesoscale quantum systems. To understand the behaviour of open quantum systems, we need accurate methods that capture the influence of the environment on the system, while managing the exponentially large Hilbert space required to describe the system. Tensor network algorithms offer an efficient way to approach this challenge. In this thesis, we develop and apply tensor network techniques to study the dynamics and steady states of various open quantum systems. The first part of the thesis focuses on the driven dissipative many body physics in coupled cavity arrays described by Born-Markov master equations. We extend transfer matrix product operator methods to Liouvillian dynamics, and utilize them to compute dynamical correlation functions and fluorescence spectrum of an infinite coupled cavity array in 1D. We also investigate thermalization, and observe the emergence of a quasi-thermal steady state with a negative effective temperature. In another study, we use infinite projected entangled pair state (iPEPS) methods to compute steady states of coupled cavity lattices in 2D. We find that a straightforward adaptation of iPEPS to Liouvillian dynamics is unstable, contradicting a recent publication in the field. The second part investigates more general systems involving strong couplings and structured environments that induce non-Markovian dynamics. We develop a powerful time-evolving matrix product operator (TEMPO) algorithm that builds on Feynman-Vernon influence functional formalism, and uses matrix product states (MPS) to represent the temporal non-Markovian correlations efficiently. We apply TEMPO to study the localization phase transition of the spin-boson model and the dynamics of two spatially separated two-level systems coupled to a common environment. Finally, we propose the Toblerone TEMPO algorithm, which extends TEMPO to many-body systems interacting with general bosonic environments.
Open quantum many body systems describe several quantum simulators, ranging from arrays of superconducting circuits hosting correlated states of light, to ultracold atoms in optical lattices in presence of controlled dissipative processes. Their theoretical understanding is hampered by the exponential scaling of their Hilbert space and by their intrinsic nonequilibrium nature, limiting the applicability of many traditional approaches. In this work we extend Dynamical Mean Field Theory (DMFT), a powerful nonperturbative technique developed for strongly correlated quantum many-body systems, to Markovian open quantum systems. Within Open-DMFT, a master equation describing a lattice of bosonic particles is mapped, at large but finite connectivity, onto an impurity problem describing a single site coupled to Markovian and non-Markovian environments, where the latter accounts for quantum fluctuations beyond Gutzwiller mean-field theory. We introduce a new non-perturbative approach to solve the impurity problem, which dresses the impurity, featuring Markovian dissipative channels, with the non-Markovian bath, in a self-consistent scheme based on a resummation of non-crossing diagrams. Finally, we apply our Open-DMFT to solve a driven-dissipative Bose-Hubbard model, which is relevant to current experiments with dissipative ultracold bosons. We show that this model undergoes a dissipative phase transition towards a superfluid phase, which can be equally interpreted as a quantum many-body synchronization transition of an array of quantum van der Pol oscillators. We show how Open-DMFT captures crucial effects due to finite lattice connectivity, such as hopping-induced dissipation and heating, which are completely missed by mean field theories and which lead to a drastic reduction of the broken symmetry phase that we interpret as a desynchronization driven by quantum fluctuations.
The study of many-body quantum physics has been recently challenged by the appearance of an increasing number of experimental platforms where genuine quantum phenomena take place notwithstanding the presence of an environment and of dissipation. Exciton-polaritons, lossy atomic and molecular gases, cavity QED arrays, arrays of trapped ions and Rydberg atoms, are only few prominent examples of a long list. The present work is ideally divided into two parts, the former is devoted to methodological results; the latter to the theoretical characterization of experimental data. Whereas the notion of equilibrium has been a fruitful guide to the development of standard many-body physics, these setups are inherently out of equilibrium and their modelisation requires the introduction of new theoretical tools. Several methods for addressing this dissipative out-of-equilibrium dynamics have been proposed, based for instance on quantum trajectories, truncated Wigner expansions, tensor networks, machine learning, but the solution of many-body physics for open quantum systems remains a formidable task. Clearly, techniques developed in the framework of Hamiltonian closed systems are a continuous source of inspiration for novel developments, and in this work we present the generalization of one such technique, the so-called flow equations, to the dissipative framework where the generator of the dynamics is a linear non-Hermitian operator; notable examples include non-Hermitian Hamiltonians and Lindblad master equations. In the second part of this work we focus on two-body lossy dynamics in strongly-correlated cold atomic gases. In particular, we derive an effective description for the long-time dynamics of a spin ½ fermionic system by extending a pre-existing theory for the bosonic case.