We develop a method to study quantum impurity models, small interacting quantum systems bilinearly coupled to an environment, in the presence of an additional Markovian quantum bath, with a generic nonlinear coupling to the impurity. We aim at computing the evolution operator of the reduced density matrix of the impurity, obtained after tracing out all the environmental degrees of freedom. First, we derive an exact real-time hybridization expansion for this quantity, which generalizes the result obtained in the absence of the additional Markovian dissipation and which could be amenable to stochastic sampling through diagrammatic Monte Carlo. Then, we obtain a Dyson equation for this quantity and we evaluate its self-energy with a resummation technique known as the noncrossing approximation. We apply this novel approach to a simple fermionic impurity coupled to a zero temperature fermionic bath and in the presence of Markovian pump, losses, and dephasing.
We develop a method to study quantum impurity models, small interacting quantum systems bilinearly coupled to an environment, in the presence of an additional Markovian quantum bath, with a generic nonlinear coupling to the impurity. We aim at computing the evolution operator of the reduced density matrix of the impurity, obtained after tracing out all the environmental degrees of freedom. First, we derive an exact real-time hybridization expansion for this quantity, which generalizes the result obtained in the absence of the additional Markovian dissipation and which could be amenable to stochastic sampling through diagrammatic Monte Carlo. Then, we obtain a Dyson equation for this quantity and we evaluate its self-energy with a resummation technique known as the noncrossing approximation. We apply this novel approach to a simple fermionic impurity coupled to a zero temperature fermionic bath and in the presence of Markovian pump, losses, and dephasing.
We develop a diagrammatic Monte Carlo method for the real-time dynamics of dissipative quantum impurity models. These are small open quantum systems with interaction and local Markovian dissipation, coupled to a large quantum bath. Our algorithm samples the hybridization expansion formulated on a single real-time contour, rather than on the double Keldysh one, as it naturally arises in the thermofield/vectorized representation of the Lindblad dynamics. We show that local Markovian dissipation generally helps the convergence of the diagrammatic Monte Carlo sampling by reducing the sign problem, thus allowing to reach longer timescales as compared to the conventional unitary case. We apply our method to an Anderson impurity model in the presence of local dephasing and discuss its effect on the charge and spin dynamics of the impurity.
This thesis focuses on the study of dissipative quantum impurity models. Quantum impurity models describe a localized quantum system, known as the impurity, that interacts with a large quantum bath, enabling the exploration of essential phenomena such as the Kondo effect in magnetic impurities and the transport properties in quantum dots. An intriguing extension of these models is the exploration of the interplay between strong correlations and dissipation. Dissipation, often described by the Lindblad master equation in its Markovian limit, introduces challenges and opportunities as it interacts with coherent many-body effects. In this thesis, we want to investigate the non-equilibrium dynamics of quantum impurity systems influenced by both strong correlations and dissipation across four main parts. We begin by developing numerical methods, specifically a diagrammatic Monte Carlo technique, to analyze the real-time dynamics of dissipative quantum impurity models. Our findings demonstrate that local Markovian dissipation enhances the convergence of Monte Carlo sampling by constraining the process to a reduced subspace of diagrams. This effectively mitigates the sign problem and facilitates the exploration of longer time scales compared to conventional unitary approaches. In the second part, we derive a self-consistent dynamical map using the Non-Crossing Approximation to study the long-term behavior of dissipative quantum impurity models. This method effectively resums diagrams into a self-consistent Dyson equation, capturing essential non-Markovian environmental effects, including the Kondo effect. In a third part, we examine the relationship between Kondo physics and Markovian dissipation by investigating an interacting quantum dot coupled to a metallic bath under continuous charge monitoring. We identify a crossover from Kondo screening to the Quantum Zeno effect, effectively described by a non-Hermitian Kondo model in the weak dissipation regime. Finally, the last part explores the effects of particle losses on the Kondo effect and the Kondo-Zeno crossover. Preliminary results indicate that while single-particle losses reduce many-body effects, two-body losses introduce intriguing phenomena, such as spin decay and hints of the Kondo-Zeno crossover, along with an unusual spectral function exhibiting a low-frequency lineshape.
Theoretical studies of electron and energy transport in nanostructures usually involve nonequilibrium quantum impurity models, which feature a locally interacting quantum system (impurity) coupled to noninteracting bosonic and/or electronic baths each at its own equilibrium. In these models, the interplay of strong correlation, nonequilibrium effects and dissipation leads to rich and complex phenomena. As a nonequilibrium strongly correlated many-body problem, these models, except at some rare cases, have no exact analytical solution, and hence their treatment heavily relies upon numerical techniques. The goal of this thesis is to develop numerically accurate and relatively low-cost techniques for treatment of strongly correlated open nonequilibirum quantum systems. Such impurity solvers are in high demand both as standalone methods for simulation of response in nanoscale junctions and as part of divide-and-conquer schemes such as dynamical mean field theory (DMFT).For nonequilibirum open quantum systems, two characteristic energy scales governing the physics are intra-system interactions and strength of the system-baths couplings. Thus, systems strongly coupled to their baths with weak intra-system interactions or weakly coupled systems with strong local interactions can be treated within developed perturbation theory approaches. My research focuses on development of nonequilibirum dual techniques, which are capable to accurately treat situations where perturbative techniques are inapplicable due to absence of a small parameter in the system and which reduce to standard perturbative expansions in the limits of weak coupling or weak intra-system interaction. \tDual techniques rely on ability to solve exactly a simplified reference system. The latter includes exactly all intra-system interactions and accounts approximately for system-bath couplings. Dual techniques (described below) allow to account for deviation from the correct hybridization function in a organized and controlled manner of superperturbative expansion. Here, I use Markovian Lindbald quantum master equation (QME) as a solver for the reference system. In the thesis I discuss possibility of mapping between true non-Markov dynamics and simpler Markov Lindblad QME (mapping between physical and reference systems), numerical methods for the mapping and for solution of the resulting Lindblad QME, and application of the mapping in novel dual formulations: auxiliary-QME dual fermion and dual boson approaches. \tFirst I analytically prove that the continuous fermionic baths that the impurity is coupled to can be mapped onto a set of dissipative auxiliary modes under Lindblad QME. Based on this mapping, complicated non-Markovian dynamics of the open fermionic system can be treated by a much simpler Markovian Lindblad QME. However, nice mapping generally requires many auxiliary modes, and numerically only a small number of modes can be exactly diagonalized due to exponential increase of the Hilbert space. To perturbatively treat the difference between the true bath and the small-size auxiliary one, we introduce dual-fermion method. It can directly target the steady state and avoid long-time propagation. Besides current-voltage relation, it can also provide spectral functions, indicating itself as a potential impurity solver for DMFT. Furthermore, we generalize dual-fermion method to consider additionally the effects of bosonic environment; hence it is called dual-boson method. It can be used for studying electron-phonon interaction and light-matter interaction in nanojunctions. By comparing with numerically exact results we argue that both methods are very accurate and relatively cheap. \tFinally, matrix product state representation of the Lindblad QME enables us to treat a large number of auxiliary modes, which provide an exponentially improved mapping. Adding counting field to the auxiliary QME, we can compute the generating function of time-dependent full-counting statistics of electron transport, which yields not only electronic current but also higher-order cumulants (e.g., zero-frequency noise).
We introduce and apply a numerically exact method for investigating the real-time dissipative dynamics of quantum impurities embedded in a macroscopic environment beyond the weak-coupling limit. We focus on the spin-boson Hamiltonian that describes a two-level system interacting with a bosonic bath of harmonic oscillators. This model is archetypal for investigating dissipation in quantum systems and tunable experimental realizations exist in mesoscopic and cold-atom systems. It finds abundant applications in physics ranging from the study of decoherence in quantum computing and quantum optics to extended dynamical mean-field theory. Starting from the real-time Feynman-Vernon path integral, we derive an exact stochastic Schr\odinger equation that allows to compute the full spin density matrix and spin-spin correlation functions beyond weak coupling. We greatly extend our earlier work (P. P. Orth, A. Imambekov, K. Le Hur, Phys. Rev. A {\bf 82}, 032118 (2010)) by fleshing out the core concepts of the method and by presenting a number of interesting applications. Methodologically, we present an analogy between the dissipative dynamics of a quantum spin and that of a classical spin in a random magnetic field. This analogy is used to recover the well-known Non-Interacting-Blip-Approximation (NIBA) in the weak-coupling limit. We explain in detail how to compute spin-spin autocorrelation functions. As interesting applications of our method, we explore the non-Markovian effects of the initial spin-bath preparation on the dynamics of the coherence σx(t) and of σz(t) under a Landau-Zener sweep of the bias field. We also compute to a high precision the asymptotic long-time dynamics of σz(t) without bias and demonstrate the wide applicability of our approach by calculating the spin dynamics at non-zero bias and different temperatures.
In this thesis, we describe methods for solving quantum impurity models out of equilibrium, and consider several applications of these methods. Quantum impurity models describe systems in which a small subsystem (the impurity) is coupled to a large external environment (the bath). These models are useful in describing a wide variety of physical systems, including quantum dots, absorbed adatoms, and molecular junctions. They are also important for the study of lattice models, such as the Hubbard model, where they appear as auxiliary problems within the framework of dynamical mean field theory (DMFT). Many interesting phenomena in these systems, such as transport through a quantum dot, or optical pumping of a Mott insulator, occur away from equilibrium. Theoretical description of these phenomena requires methods that are capable of handling both real time evolution and strong correlation effects presenting a formidable challenge. In this thesis, we describe several methods for the solution of non-equilibrium quantum impurity models which meet this challenge. These methods are based on an expansion in the coupling between the impurity and the bath known as the hybridization expansion. Chapter 1 gives an introduction to the problems addressed in this thesis. Chapter 2 provides a brief overview of Keldysh formalism for nonequilibrium problems. Chapter 3 describes the bare hybridization expansion for the Anderson impurity model as well as its bold extensions. In particular we derive the non-crossing approximation (NCA), the one-crossing approximation (OCA), the inchworm expansion, and a bold expansion for the impurity Green's function. Chapter 4 describes the implementation of a quantum Monte Carlo algorithm based on the inchworm expansion. The rest of the thesis describes applications of these algorithms to problems of physical interest. Chapters 5 and 6 apply the methods previously described to understanding the Kondo effect in quantum dot systems. The first chapter examines the behavior of the Kondo effect under an applied voltage bias, a classic non-equilibrium strongly-correlated problem. The second chapter examines different ways of measuring the cloud of electrons which screen the impurity spin in the Kondo regime. Chapters 7 and 8 apply the same methods to the Hubbard model using the Dynamical mean-field theory (DMFT) mapping. The first chapter studies the Hubbard model in equilibrium using the numerically exact inchworm quantum Monte Carlo method. The second chapter studies the possibility of dynamic switching between the metal and Mott-insulator states of the Hubbard model using the one-crossing approximation (OCA). Finally chapter 9 gives some overall conclusions and suggests future directions for the development of hybridization expansion methods.
The non-Markovian master equation for open quantum systems is obtained by generalization of the ordinary Zwanzig-Nakajima (ZN) projection technique. To this end, a coupled chain of equations for the reduced density matrices of the bath ϱB(t) and the system ϱS(t) are written down. A formal solution of the equation for ϱB(t), having been inserted in the equation for the reduced density matrix of the system, in the 2-nd approximation in interaction yields a very specific extra term in the generalized master equation. This term, being nonlinear in ϱS(t), is related to the intrinsic bath dynamics and vanishes in the Markovian limit. To verify the consistence and robustness of our approach, we applied the generalized ZN projection scheme to a simple dephasing model. It is shown that consideration of the lowest order in interaction is insufficient to describe time evolution of the system coherence adequately. We explain this fact by analyzing the exact and approximate forms of ϱB(t) and give some hints how to take the dynamic correlations (which originate from the spin-bath coupling) into account.