Numerical Renormalization Group for Impurity Quantum Phase Transitions: Structure of Critical Fixed Points

The numerical renormalization group method is used to investigate zero temperature phase transitions in quantum impurity systems, in particular in the particle-hole symmetric soft-gap Anderson model. The model displays two stable phases whose fixed points can be built up of non-interacting single-particle states. In contrast, the quantum phase transitions turn out to be described by interacting fixed points, and their excitations cannot be described in terms of free particles. We show that the structure of the many-body spectrum of these critical fixed points can be understood using renormalized perturbation theory close to certain values of the bath exponents which play the role of critical dimensions. Contact is made with perturbative renormalization group calculations for the soft-gap Anderson and Kondo models. A complete description of the quantum critical many-particle spectra is achieved using suitable marginal operators; technically this can be understood as epsilon-expansion for full many-body spectra.Comment: 14 pages, 12 figure

Spectral Densities of Response Functions for the O(3) Symmetric Anderson and Two Channel Kondo Models

The O(3) symmetric Anderson model is an example of a system which has a stable low energy marginal Fermi liquid fixed point for a certain choice of parameters. It is also exactly equivalent, in the large U limit, to a localized model which describes the spin degrees of freedom of the linear dispersion two channel Kondo model. We first use an argument based on conformal field theory to establish this precise equivalence with the two channel model. We then use the numerical renormalization group (NRG) approach to calculate both one-electron and two-electron response functions for a range of values of the interaction strength U. We compare the behaviours about the marginal Fermi liquid and Fermi liquid fixed points and interpret the results in terms of a renormalized Majorana fermion picture of the elementary excitations. In the marginal Fermi liquid case the spectral densities of all the Majorana fermion modes display a |omega| dependence on the lowest energy scale, and in addition the zero Majorana mode has a delta function contribution. The weight of this delta function is studied as a function of the interaction U and is found to decrease exponentially with U for large U. Using the equivalence with the two channel Kondo model in the large U limit, we deduce the dynamical spin susceptibility of the two channel Kondo model over the full frequency range. We use renormalized perturbation theory to interpret the results and to calculate the coefficient of the ln omega divergence found in the low frequency behaviour of the T=0 dynamic susceptibility.Comment: 26 pages, 18 figures, to be published in Eur. Phys. J.

Gap formation and soft phonon mode in the Holstein model

We investigate electron-phonon coupling in many-electron systems using dynamical mean-field theory in combination with the numerical renormalization group. This non-perturbative method reveals significant precursor effects to the gap formation at intermediate coupling strengths. The emergence of a soft phonon mode and very strong lattice fluctuations can be understood in terms of Kondo-like physics due to the development of a double-well structure in the effective potential for the ions

Numerical Renormalization Group Calculations for the Self-energy of the impurity Anderson model

We present a new method to calculate directly the one-particle self-energy of an impurity Anderson model with Wilson's numerical Renormalization Group method by writing this quantity as the ratio of two correlation functions. This way of calculating Sigma(z) turns out to be considerably more reliable and accurate than via the impurity Green's function alone. We show results for the self-energy for the case of a constant coupling between impurity and conduction band (ImDelta = const) and the effective Delta(z) arising in the Dynamical Mean Field Theory of the Hubbard model. Implications to the problem of the metal-insulator transition in the Hubbard model are also discussed.Comment: 18 pages, 9 figures, submitted to J. Phys.: Condens. Matte

On X-ray-singularities in the f-electron spectral function of the Falicov-Kimball model

The f-electron spectral function of the Falicov-Kimball model is calculated within the dynamical mean-field theory using the numerical renormalization group method as the impurity solver. Both the Bethe lattice and the hypercubic lattice are considered at half filling. For small U we obtain a single-peaked f-electron spectral function, which --for zero temperature-- exhibits an algebraic (X-ray) singularity ($|\omega|^{-\alpha}$) for $\omega \to 0$. The characteristic exponent $\alpha$ depends on the Coulomb (Hubbard) correlation U. This X-ray singularity cannot be observed when using alternative (Keldysh-based) many-body approaches. With increasing U, $\alpha$ decreases and vanishes for sufficiently large U when the f-electron spectral function develops a gap and a two-peak structure (metal-insulator transition).Comment: 8 pages, 8 figures, revte

Numerical Renormalization Group for Bosonic Systems and Application to the Subohmic Spin-Boson Model

We describe the generalization of Wilson's Numerical Renormalization Group method to quantum impurity models with a bosonic bath, providing a general non-perturbative approach to bosonic impurity models which can access exponentially small energies and temperatures. As an application, we consider the spin-boson model, describing a two-level system coupled to a bosonic bath with power-law spectral density, J(omega) ~ omega^s. We find clear evidence for a line of continuous quantum phase transitions for subohmic bath exponents 0<s<1; the line terminates in the well-known Kosterlitz-Thouless transition at s=1. Contact is made with results from perturbative renormalization group, and various other applications are outlined.Comment: 4 pages, 5 figs, (v2) final version as publishe

Anderson impurity in pseudo-gap Fermi systems

We use the numerical renormalization group method to study an Anderson impurity in a conduction band with the density of states varying as rho(omega) \propto |omega|^r with r>0. We find two different fixed points: a local-moment fixed point with the impurity effectively decoupled from the band and a strong-coupling fixed point with a partially screened impurity spin. The specific heat and the spin-susceptibility show powerlaw behaviour with different exponents in strong-coupling and local-moment regime. We also calculate the impurity spectral function which diverges (vanishes) with |omega|^{-r} (|\omega|^r) in the strong-coupling (local moment) regime.Comment: 8 pages, LaTeX, 4 figures includes as eps-file

Numerical renormalization group study of the symmetric Anderson-Holstein model: phonon and electron spectral functions

We study the symmetric Anderson-Holstein (AH) model at zero temperature with Wilson's numerical renormalization group (NRG) technique to study the interplay between the electron-electron and electron-phonon interactions. An improved method for calculating the phonon propagator using the NRG technique is presented, which turns out to be more accurate and reliable than the previous works in that it calculates the phonon renormalization explicitly and satisfies the boson sum rule better. The method is applied to calculate the renormalized phonon propagators along with the electron propagators as the onsite Coulomb repulsion $U$ and electron-phonon coupling constant $g$ are varied. As $g$ is increased, the phonon mode is successively renormalized, and for $g \gtrsim g_{co}$ crosses over to the regime where the mode splits into two components, one of which approaches back to the bare frequency and the other develops into a soft mode. The initial renormalization of the phonon mode, as $g$ is increased from 0, depends on $U$ and the hybridization $\Delta$; it gets softened (hardened) for $U \gtrsim (\lesssim) U_s (\Delta)$. Correlated with the emergence of the soft mode is the central peak of the electron spectral function severely suppressed. These NRG calculations will be compared with the standard Green's function results for the weak coupling regime to understand the phonon renormalization and soft mode.Comment: 18 pages, 4 figures. Submitted to Phys. Rev.

Spectral properties of locally correlated electrons in a BCS superconductor

We present a detailed study of the spectral properties of a locally correlated site embedded in a BCS superconducting medium. To this end the Anderson impurity model with superconducting bath is analysed by numerical renormalisation group (NRG) calculations. We calculate one and two-particle dynamic response function to elucidate the spectral excitation and the nature of the ground state for different parameter regimes with and without particle-hole symmetry. The position and weight of the Andreev bound states is given for all relevant parameters. We also present phase diagrams for the different ground state parameter regimes. This work is also relevant for dynamical mean field theory extensions with superconducting symmetry breaking.Comment: 22 pages, 12 figure