No Arabic abstract
The density matrix renormalization group method is applied to obtain 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. 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. 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 the first 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 three fold degenerate at the valence fluctuating point. 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.
We introduce a block Lanczos (BL) recursive technique to construct quasi-one-dimensional models, suitable for density-matrix renormalization group (DMRG) calculations, from single- as well as multiple-impurity Anderson models in any spatial dimensions. This new scheme, named BL-DMRG method, allows us to calculate not only local but also spatially dependent static and dynamical quantities of the ground state for general Anderson impurity models without losing elaborate geometrical information of the lattice. We show that the BL-DMRG method can be easily extended to treat a multi-orbital Anderson impurity model. We also show that the symmetry adapted BL bases can be utilized, when it is appropriate, to reduce the computational cost. As a demonstration, we apply the BL-DMRG method to three different models for graphene: (i) a single adatom on the honeycomb lattice, (ii) a substitutional impurity in the honeycomb lattice, and (iii) an effective model for a single carbon vacancy in graphene. Our analysis reveals that, for the particle-hole symmetric case at half filling of electron density, the ground state of model (i) behaves as an isolated magnetic impurity with no Kondo screening while the ground state of the other two models forms a spin singlet state. We also calculate the real-space dependence of the spin-spin correlation functions between the impurity site and the conduction sites for these three models. Our results clearly show that, reflecting the presence of absence of unscreened magnetic moment at the impurity site, the spin-spin correlation functions decay as $r^{-3}$, differently from the non-interacting limit ($r^{-2}$), for model (i) and as $ r^{-4}$, exactly the same as the non-interacting limit, for models (ii) and (iii) in the asymptotic $r$, where $r$ is the distance between the impurity site and the conduction site.
We have developed a new efficient and accurate impurity solver for the single impurity Anderson model (SIAM), which is based on a non-perturbative recursion technique in a space of operators and involves expanding the self-energy as a continued fraction. The method has no special occupation number or temperature restrictions; the only approximation is the number of levels of the continued fraction retained in the expansion. We also show how this approach can be used as a new approach to Dynamical Mean Field Theory (DMTF) and illustrate this with the Hubbard model. The three lowest orders of recursion give the Hartree-Fock, Hubbard I, and Hubbard III approximations. A higher level of recursion is able to reproduce the expected 3-peak structure in the spectral function and Fermi liquid behavior.
We analyze the ground-state energy, magnetization, magnetic susceptibility, and Kondo screening cloud of the symmetric single-impurity Anderson model (SIAM) that is characterized by the band width $W$, the impurity interaction strength $U$, and the local hybridization $V$. We compare Gutzwiller variational and magnetic Hartree-Fock results in the thermodynamic limit with numerically exact data from the Density-Matrix Renormalization Group (DMRG) method on large rings. To improve the DMRG performance, we use a canonical transformation to map the SIAM onto a chain with half the system size and open boundary conditions. We compare to Bethe-Ansatz results for the ground-state energy, magnetization, and spin susceptibility that become exact in the wide-band limit. Our detailed comparison shows that the field-theoretical description is applicable to the SIAM on a ring for a broad parameter range. Hartree-Fock theory gives an excellent ground-state energy and local moment for intermediate and strong interactions. However, it lacks spin fluctuations and thus cannot screen the impurity spin. The Gutzwiller variational energy bound becomes very poor for large interactions because it does not describe properly the charge fluctuations. Nevertheless, the Gutzwiller approach provides a qualitatively correct description of the zero-field susceptibility and the Kondo screening cloud. The DMRG provides excellent data for the ground-state energy and the magnetization for finite external fields. At strong interactions, finite-size effects make it extremely difficult to recover the exponentially large zero-field susceptibility and the mesoscopically large Kondo screening cloud.
We introduce a method to obtain the specific heat of quantum impurity models via a direct calculation of the impurity internal energy requiring only the evaluation of local quantities within a single numerical renormalization group (NRG) calculation for the total system. For the Anderson impurity model, we show that the impurity internal energy can be expressed as a sum of purely local static correlation functions and a term that involves also the impurity Green function. The temperature dependence of the latter can be neglected in many cases, thereby allowing the impurity specific heat, $C_{rm imp}$, to be calculated accurately from local static correlation functions; specifically via $C_{rm imp}=frac{partial E_{rm ionic}}{partial T} + 1/2frac{partial E_{rm hyb}}{partial T}$, where $E_{rm ionic}$ and $E_{rm hyb}$ are the energies of the (embedded) impurity and the hybridization energy, respectively. The term involving the Green function can also be evaluated in cases where its temperature dependence is non-negligible, adding an extra term to $C_{rm imp}$. For the non-degenerate Anderson impurity model, we show by comparison with exact Bethe ansatz calculations that the results recover accurately both the Kondo induced peak in the specific heat at low temperatures as well as the high temperature peak due to the resonant level. The approach applies to multiorbital and multichannel Anderson impurity models with arbitrary local Coulomb interactions. An application to the Ohmic two state system and the anisotropic Kondo model is also given, with comparisons to Bethe ansatz calculations. The new approach could also be of interest within other impurity solvers, e.g., within quantum Monte Carlo techniques.
We investigate the effect of the Coulomb interaction, $U_{cf}$, between the conduction and f electrons in the periodic Anderson model using the density-matrix renormalization-group algorithm. We calculate the excitation spectrum of the half-filled symmetric model with an emphasis on the spin and charge excitations. In the one-dimensional version of the model it is found that the spin gap is smaller than the charge gap below a certain value of $U_{cf}$ and the reversed inequality is valid for stronger $U_{cf}$. This behavior is also verified by the behavior of the spin and density correlation functions. We also perform a quantum information analysis of the model and determine the entanglement map of the f and conduction electrons. It is revealed that for a certain $U_{cf}$ the ground state is dominated by the configuration in which the conduction and f electrons are strongly entangled, and the ground state is almost a product state. For larger $U_{cf}$ the sites are occupied alternatingly dominantly by two f electrons or by two conduction electrons.