The bosonic single-impurity Anderson model (B-SIAM) is studied to understand the local dynamics of an atomic quantum dot (AQD) coupled to a Bose-Einstein condensation (BEC) state, which can be implemented to probe the entanglement and the decoherence of a macroscopic condensate. Our recent approach of the numerical renormalization-group calculation for the B-SIAM revealed a zero-temperature phase diagram, where a Mott phase with local depletion of normal particles is separated from a BEC phase with enhanced density of the condensate. As an extension of the previous work, we present the calculations of the local dynamical quantities of the B-SIAM which reinforce our understanding of the physics in the Mott and the BEC phases.
The bosonic single-impurity Anderson model (B-SIAM) is studied to understand the local dynamics of an atomic quantum dot (AQD) coupled to a Bose-Einstein condensation (BEC) state, which can be implemented to probe the entanglement and the decoherence of a macroscopic condensate. Our recent approach of the numerical renormalization-group calculation for the B-SIAM revealed a zero-temperature phase diagram, where a Mott phase with local depletion of normal particles is separated from a BEC phase with enhanced density of the condensate. As an extension of the previous work, we present the calculations of the local dynamical quantities of the B-SIAM which reinforce our understanding of the physics in the Mott and the BEC phases.
The density matrix renormalization group method is introduced in energy space to study Anderson impurity models. The method allows for calculations in the thermodynamic limit and is advantageous for studying not only the dynamical properties, but also the quantum entanglement of the ground state at the vicinity of an impurity quantum phase transition. This method is applied to obtain numerically exactly the ground-state phase diagram of the single-impurity Anderson model on the honeycomb lattice at half-filling. The calculation of local static quantities shows that the phase diagram contains two distinct phases, the local moment (LM) phase and the asymmetric strong coupling (ASC) phase, but no Kondo screening phase. These results are supported by the local spin and charge excitation spectra, which exhibit qualitatively different behavior in these two phases and also reveal the existence of the valence fluctuating point at the phase boundary. For comparison, we also study the low-energy effective pseudogap Anderson model using the method introduced here. Although the high-energy excitations are obviously different, we find that the ground-state phase diagram and the asymptotically low-energy excitations are in good quantitative agreement with those for the single-impurity Anderson model on the honeycomb lattice, thus providing a quantitative justification for the previous studies based on low-energy approximate approaches. Furthermore, we find that the lowest entanglement level is doubly degenerate for the LM phase, whereas it is singlet for the ASC phase and is accidentally threefold degenerate at the valence fluctuating point. This should be contrasted with the degeneracy of the energy spectrum because the ground state is found to be always singlet. Our results therefore clearly demonstrate that the low-lying entanglement spectrum can be used to determine with high accuracy the phase boundary of the impurity quantum 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.
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.
The dynamical mean-field theory (DMFT), which is successful in the study of strongly correlated fermions, was recently extended to boson systems [K. Byczuk and D. Vollhardt, Phys. Rev. B 77, 235106 (2008)]. In this paper, we employ the bosonic DMFT to study the Bose-Hubbard model which describes on-site interacting bosons in a lattice. Using exact diagonalization as the impurity solver, we get the DMFT solutions for the Green's function, the occupation density, as well as the condensate fraction on a Bethe lattice. Various phases are identified: the Mott insulator, the Bose-Einstein condensed (BEC) phase, and the normal phase. At finite temperatures, we obtain the crossover between the Mott-type regime and the normal phase, as well as the BEC-to-normal phase transition. Phase diagrams on the $\ensuremath{\mu}/U\text{\ensuremath{-}}\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{t}/U$ plane and on the $T/U\text{\ensuremath{-}}\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{t}/U$ plane are produced ($\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{t}$ is the scaled hopping amplitude). We compare our results with the previous ones, and discuss the implication of these results to experiments.
The single impurity Anderson model (SIAM) is one of the most recurrent model Hamiltonians in condensed matter physics, playing an important role in the description of systems where electron-electron interactions can not be neglected, such as strongly interacting quantum dots and doped metals with magnetic impurities. In order to address the variety of scenarios in which the model can be applied many impurity solvers were developed for the SIAM. It is consensus that there is no universal method that provides reliable results for the model parameter space of experimental interest. The main methods
found in the literature di ers signi cantly on the computational cost involved in the numerical implementation and in the range of the model parameters where the used method is reliable. In this undergraduate dissertation we study the reliability of some impurity solvers for the SIAM, namely, the equations of motion (EOM) method, the numerical renormalization group (NRG), the non-crossing approximation (NCA), and the one crossing approximation (OCA). We study the Green's function of a quantum dot using the aforementioned impurity solvers and then calculate the density of states
as a function of the temperature and the chemical potential of the quantum dot . Our analysis is focused on the Coulomb blockade and the Kondo physics. We report that some standard approximations based on the EOM technique display a rather unexpected poor behaviour in the Coulomb blockade regime even at high temperatures.
We formulate a bosonic dynamical mean-field theory (B-DMFT) which provides a comprehensive, thermodynamically consistent framework for the theoretical investigation of correlated lattice bosons. The B-DMFT is applicable for arbitrary values of the coupling parameters and temperature and becomes exact in the limit of high spatial dimensions d or coordination number Z of the lattice. In contrast to its fermionic counterpart, the construction of the B-DMFT requires different scalings of the hopping amplitudes with Z depending on whether the bosons are in their normal state or in the Bose-Einstein condensate (BEC). A detailed discussion of how this conceptual problem can be overcome by performing the scaling in the action rather than in the Hamiltonian itself is presented. The B-DMFT treats normal and condensed bosons on equal footing and thus includes the effects caused by their dynamic coupling. It reproduces all previously investigated limits in parameter space, such as the Beliaev-Popov and Hartree-Fock-Bogoliubov approximations, and generalizes the existing mean-field theories of interacting bosons. The self-consistency equations of the B-DMFT are those of a bosonic single-impurity coupled to two reservoirs corresponding to bosons in the condensate and in the normal state, respectively. We employ the B-DMFT to solve a model of itinerant and localized, interacting bosons analytically. The local correlations are found to enhance the condensate density and the BEC transition temperature TBEC. This effect may be used experimentally to increase TBEC of bosonic atoms in optical lattices.