We present a method to study the time evolution of the single impurity Anderson model which exploits a mean field decoupling of the interacting impurity and the non-interacting bath (in form of a chain). This is achieved by the introduction of a pair of auxiliary Majorana fermions between the impurity and the chain. After decoupling, we obtain a self-consistent set of equations for the impurity and chain. First, we study the behavior of the system in equilibrium at zero temperature. We obtain a phase transition as a function of the interaction at the impurity and the coupling between the impurity and the chain between the Kondo regime, where the mean field parameters are zero and, hence, we have a well-defined spin at the impurity, to a phase where mean field parameters acquire finite values leading to a screening of the impurity spin by conduction bath electrons. In the latter case, we observe charge and spin fluctuations at the impurity site. Starting from this equilibrium ground state at zero temperature we quench in the interaction strength at the impurity and/or the hybridization strength between the impurity and the chain and study the time evolution of the system. We find that for quenches to weak to intermediate coupling the system converges to the equilibrium state defined by the final set of parameters after the quench. We analyze the oscillation frequency and as well as the thermalization rate during this quench. A quench to a strong interaction value results in persistent oscillations and a trapping of the system in a non-thermal state. We speculate that these two regimes of different long-time behavior are separated by a dynamical phase transition.
We present a method to study the time evolution of the single impurity Anderson model which exploits a mean field decoupling of the interacting impurity and the non-interacting bath (in form of a chain). This is achieved by the introduction of a pair of auxiliary Majorana fermions between the impurity and the chain. After decoupling, we obtain a self-consistent set of equations for the impurity and chain. First, we study the behavior of the system in equilibrium at zero temperature. We obtain a phase transition as a function of the interaction at the impurity and the coupling between the impurity and the chain between the Kondo regime, where the mean field parameters are zero and, hence, we have a well-defined spin at the impurity, to a phase where mean field parameters acquire finite values leading to a screening of the impurity spin by conduction bath electrons. In the latter case, we observe charge and spin fluctuations at the impurity site. Starting from this equilibrium ground state at zero temperature we quench in the interaction strength at the impurity and/or the hybridization strength between the impurity and the chain and study the time evolution of the system. We find that for quenches to weak to intermediate coupling the system converges to the equilibrium state defined by the final set of parameters after the quench. We analyze the oscillation frequency and as well as the thermalization rate during this quench. A quench to a strong interaction value results in persistent oscillations and a trapping of the system in a non-thermal state. We speculate that these two regimes of different long-time behavior are separated by a dynamical phase transition.
An essential role for the description of correlated quantum many-particle systems is played by impurity problems. These consist of a small number of localized orbitals, where the electrons are subject to a Coulomb-interaction, coupled to noninteracting bath degrees of freedom. Impurity models are used for the description of magnetic impurities in metals, of quantum points in nanostructures, and in the context of the dynamical mean field theory (DMFT). For their solution a large number of numerical methods exists. We mention here quantum Monte Carlo methods, the numerical renormalization group, or the exact diagonalization method. In this thesis we introduce a new approach to impurity problems, which is based on the functional renormalization group (fRG). In contrast to the methods mentioned above this approach is approximative, but it can be applied with comparatively lower numerical effort. The aim of this thesis is to examine the advantages and limits of this method. An important role in this connection is played by the single impurity Anderson model (SIAM), which is already understood very well with other methods and is seen as a prototype for models with strong local correlations. The spin fluctuations of this model are governed by the so-called Kondo energy scale and a central question of this thesis is whether this scale can be reproduced with our approach. We begin with a short introduction to the Green's function formalism for correlated quantum many particle systems. Then the fRG flow equations for the one-particle irreducible vertex functions are deduced, and several approximation schemes are discussed. In chapter 4 we introduce the fRG scheme used in this thesis. At first the SIAM is mapped to a semi-infinite chain, in which the interacting orbital is given by the first site of this chain. The system is then subdivided into two parts. The first part, which is called 'core', consists of the interacting orbital and the first L bath sites. The remaining bath sites form the second part. At the beginning both parts are decoupled from each other such that the core can be solved exactly. Starting with this exact solution the coupling between the core and the remaining bath sites is switched on slowly in a renormalization group flow. This way it is possible to calculate the one-particle and two-particle correlation functions on the correlated site. We call this flow scheme 'hybridization flow'. The flow equations are formulated in an effective theory on the first bath site outside the core, which turns out to be advantageous compared to other implementations of the hybridization flow. In chapter 5 this hybridization flow method is applied to the SIAM with semi-elliptic bath density of states. We discuss the differences which arise for different core-sizes (L = 0,1,2,3) and for different truncations of the fRG flow equations. In the cases L=0 and L=2 the local Fermi liquid properties of the SIAM and the associated Kondo scale can not be reproduced in both approximation schemes. For the core sizes L=1 and L=3 the ground state of the core is already a spin singlet state in the beginning of the flow and the Fermi liquid properties are reproduced. The Kondo energy scale, which determines the width of the central quasi-particle peak and the size of the local spin-susceptibility, is however not accurately resolved. This is partly due to the finite temperature used in our implementation such that for larger interaction strength the temperature is above the Kondo energy scale. The application of the hybridization flow method to quantum impurity problems in the context of the DMFT is discussed in chapter 6. In this approach lattice models with a local interaction, like the Hubbard model, are mapped to an effective impurity problem for a single lattice site embedded in a dynamical mean field representing the influence of the other electrons. This scheme is exact in infinitely many dimensions and a nonperturbative approximation method for finite-dimensional systems. In particular the Mott metal-insulator transition in the Hubbard model is successfully described by this method. One class of quantities, which is easily accessible from the hybridization flow scheme, is given by the local two-particle correlation functions. These play an essential role in non-local extensions of the DMFT. Therefore we mainly focus on these quantities. At first we apply our flow scheme to the Hubbard model on the Bethe lattice in infinite dimensions and calculate the local one-particle irreducible vertex functions in the insulating and the metallic phase in good agreement with previous calculations, which use the exact diagonalization method as impurity solver. After this we extend our flow scheme to a cluster DMFT method, which includes short-range antiferromagnetic correlations and calculate again the local and non-local vertex functions for the two-dimensional Hubbard model. We compare these to the results of the single-site DMFT.
Quantum impurity models with strong electron correlations, such as the paradigmatic Anderson Impurity Model (AIM), are central to our understanding of a range of physical phenomena including local moment formation, Coulomb blockade and Kondo screening. They describe magnetic atoms and molecules on surfaces, quantum dot circuits, and correlated materials through dynamical mean field theory. The physics of such systems in strongly non-equilibrium conditions is particularly complex and challenging to capture, whereas Gaussian models of free fermions can be easily solved. Here we show that the time-evolving dynamics of the AIM after a quench can be described by a completely non-interacting version of the model, at the expense of coupling to additional static auxiliary degrees of freedom. Starting from the full solution of the quenched AIM using ED and DMRG, we study the properties of this mapping using numerical optimization, and uncover intriguing structure in the auxiliary system. The method allows us to understand interacting non-equilibrium dynamics through the simpler lens of an effective non-interacting system of larger dimension.
Using the non-perturbative numerical renormalisation group (NRG) we analyse purely fermionic as well as more complex Bose-Fermi quantum impurity systems on their equilibrium properties and beyond: We present an extensive exploration of real-time dynamics for various quenches using the time-dependent NRG. For the purely fermionic system, we use the Single-impurity Anderson model (SIAM) in which the impurity is a localised orbital hybridising with a surrounding conduction band. The model enables charge scattering with a pseudo-gap scattering rate and features a quantum critical point separating the local moment (LM) phase with an unscreened spin moment from a symmetric strong coupling (SSC) phase with a fractionally screened spin moment. The real-time dynamics of the local double occupancy nicely thermalises for quenches regarding the SSC phase. Whereas for quenches within or into the LM phase the double occupancy systematically deviates from thermal equilibrium for long times. A comparison with results of a recently published time-dependent Gutwiller approach reveals clear discrepancies since this approach restricts the formation of the local moment only on the impurity site. We extend the system by an additional bosonic bath leading to the Bose-Fermi Anderson model (BFAM) and study the influence of the additional bosonic degrees of freedom on the real-time dynamics. In equilibrium the bosonic bath features an additional quantum critical point leading to a localised phase once the coupling to a sub-ohmic bosonic bath exceeds a critical value. In general the coupled bosonic bath leads to a renormalisation of the Coulomb interaction. We show that the dynamics is influenced by the specific choice of the bath exponents. For a F-type combination the non-equilibrium dynamics can be completely reproduced by an effective purely fermionic SIAM while for a B-type combination the nature of the bosonic bath and the coupling strength has a profound impact on the electron dynamics which can only be partially understood by an effective SIAM. Furthermore, we investigate the spin dynamics of an initially polarised impurity spin which is coupled to a fermionic and a bosonic bath leading to the Bose-Fermi Kondo model (BFKM). Through the coupling to the fermionic bath the polarisation vanishes with time. We show that in analogy to the spin-boson model (SBM) an initially coupling to the bosonic bath affects the real-time dynamics fundamentally leading to a remaining polarisation for long times.
This work is motivated by the fact that the investigation of non-equilibrium phenomena in strongly correlated electron systems has developed into one of the most active and exciting branches of condensed matter physics as it provides rich new insights that could not be obtained from the study of equilibrium situations. However, a theoretical description of those phenomena is missing. Therefore, in this thesis, we develop a numerical method that can be used to study two minimal models -- the Hubbard model and the Anderson impurity model with general parameter range and time dependence. We begin by introducing the theoretical framework and the general features of the Hubbard model. We then describe the dynamical mean field theory (DMFT), which was first invented by Georges in 1992. It provides a feasible way to approach strongly correlated electron systems and reduces the complexity of the calculations via a mapping of lattice models onto quantum impurity models subject to a self-consistency condition. We employ the non-equilibrium extension of DMFT and map the Hubbard model to the single impurity Anderson model (SIAM). Since the fundamental component of the DMFT method is a solver of the single impurity Anderson model, we continue with a description of the formalism to study the real-time dynamics of the impurity model staring at its thermal equilibrium state. We utilize the non-equilibrium strong-coupling perturbation theory and derive semi-analytical approximation methods such as the non-crossing approximation (NCA) and the one-crossing approximation (OCA). We then use the Quantum Monte-Carlo method (QMC) as a numerically exact method and present proper measurements of local observables, current and Green's functions. We perform simulations of the current after a quantum quench from equilibrium by rapidly applying a bias voltage in a wide range of initial temperatures. The current exhibits short equilibrium times and saturates upon the decrease of temperature at all times, indicating Kondo behavior both in the transient regime and in the steady state. However, this bare QMC solver suffers from a dynamical sign problem for long time propagations. To overcome the limitations of this bare treatment, we introduce the ``Inchworm algorithm'', based on iteratively reusing the information obtained in previous steps to extend the propagation to longer times and stabilize the calculations. We show that this algorithm greatly reduces the required order for each simulation and re-scales the exponential challenge to quadratic in time. We introduce a method to compute Green's functions, spectral functions, and currents for inchworm Monte Carlo and show how systematic error assessments in real time can be obtained. We illustrate the capabilities of the algorithm with a study of the behavior of quantum impurities after an instantaneous voltage quench from a thermal equilibrium state. We conclude with the applications of the unbiased inchworm impurity solver to DMFT calculations. We employ the methods for a study of the one-band paramagnetic Hubbard model on the Bethe lattice in equilibrium, where the DMFT approximation becomes exact. We begin with a brief introduction of the Mott metal insulator phase diagram. We present the results of both real time Green's functions and spectral functions from our nonequilibrium calculations. We observe the metal-insulator crossover as the on-site interaction is increased and the formation of a quasi-particle peak as the temperature is lowered. We also illustrate the convergence of our algorithms in different aspects.
We provide a detailed formulation of the recently proposed variational approach [Y. Ashida, T. Shi, M.-C. Ba\~nuls, J. I. Cirac, and E. Demler, Phys. Rev. Lett. 121, 026805 (2018)] to study ground-state properties and out-of-equilibrium dynamics for generic quantum spin-impurity systems. Motivated by the original ideas of Tomonaga, Lee, Low, and Pines, we construct a canonical transformation that completely decouples the impurity from the bath degrees of freedom. By combining this transformation with a Gaussian ansatz for the fermionic bath, we obtain a family of variational many-body states that can efficiently encode the strong entanglement between the impurity and fermions of the bath. We give a detailed derivation of equations of motions in the imaginary- and real-time evolutions on the variational manifold. We benchmark our approach by applying it to investigate ground-state and dynamical properties of the anisotropic Kondo model and compare results with those obtained using the matrix-product state (MPS) ansatz. We show that our approach can achieve an accuracy comparable to MPS-based methods with several orders of magnitude fewer variational parameters than the corresponding MPS ansatz. Comparisons to the Yosida ansatz and the exact solution from the Bethe ansatz are also discussed. We use our approach to investigate the two-lead Kondo model and analyze its long-time spatiotemporal behavior and the conductance behavior at finite bias and magnetic fields. The obtained results are consistent with the previous findings in the Anderson model and the exact solutions at the Toulouse point.
We present a numerical method for the study of correlated quantum impurity problems out of equilibrium, which is particularly suited to address steady-state properties within dynamical mean field theory. The approach, recently introduced by Arrigoni et al. [Phys. Rev. Lett. 110, 086403 (2013)], is based upon a mapping of the original impurity problem onto an auxiliary open quantum system, consisting of the interacting impurity coupled to bath sites as well as to a Markovian environment. The dynamics of the auxiliary system is governed by a Lindblad master equation whose parameters are used to optimize the mapping. The accuracy of the results can be readily estimated and systematically improved by increasing the number of auxiliary bath sites, or by introducing a linear correction. Here, we focus on a detailed discussion of the proposed approach including technical remarks. To solve for the Green's functions of the auxiliary impurity problem, a non-Hermitian Lanczos diagonalization is applied. As a benchmark, results for the steady-state current-voltage characteristics of the single-impurity Anderson model are presented. Furthermore, the bias dependence of the single-particle spectral function and the splitting of the Kondo resonance are discussed. In its present form, the method is fast, efficient, and features a controlled accuracy.