signature=69a22092d66a85285836ed87476f1378,Long-range Kondo signature of a single magnetic impurity

The interaction of a single magnetic impurity with the surrounding electron gas of a non-magnetic metal leads to fascinating phenomena in the low-temperature limit, which are encompassed by the term Kondo effectTK, a correlated electronic state develops and the impurity spin is screened. The most prominent fingerprint of this many-body singlet state is a narrow resonance at the Fermi energy ɛF in the single-particle spectrum of the impurity, called Kondo or Abrikosov–Suhl resonance. The existence of this Kondo resonance has been experimentally confirmed for dense systems with high-resolution photoemission electron spectroscopy and inverse photoemission

A theoretical prediction for the local density of states (LDOS)—the key quantity measured in STS experiments—was first provided by Újsághy and colleagues2, 3) many experiments on magnetic atoms and molecules on metal surfaces have been carried out, all of which show Kondo fingerprintset al. detected one resonance on the adatom and another localized signature in the second focus of the ellipse. However, no line shape variations with distance were reported (for a review on adatom Kondo systems see ref. 13).

In this work we follow a new route and investigate single isolated magnetic impurities buried below the surface with a low-temperature STM operating at 6 K. It has recently been shown that the anisotropy of the copper Fermi surface leads to a strong directional propagation of quasi-particles, which is called electron focusingSupplementary Information).

As a first striking example of how the Kondo effect influences the energy-dependent scattering behaviour on the millivolt scale, Fig. 1a shows STM topographies of a fourth-layer Fe impurity for different bias voltages V near zero millivolts. The local minimum of the LDOS present in the centre of the interference pattern for V <0 develops into a plateau-like maximum for V >0. In Fig. 1b the differential conductance as a function of bias voltage and one spatial coordinate y across the impurity pattern is shown. The crossover observed in the topographies occurs very close to zero bias. In Fig. 1c, four spectra for different positions A–D are shown, illustrating very clearly that a single Kondo atom buried below the surface of copper induces long-range spectral signatures that depend on the distance to the impurity. A second possibility for investigating the Kondo effect versus distance is to look at impurities situated at different depths d below the surface. Figure 2 shows, as an example, STS data of Co atoms, comparing single spectra (purple curves) measured directly above the impurities (y=0). The lateral variations of dI/dV are again depicted as sections (upper part of Fig. 2). All show a constriction of the pattern for positive bias voltages. The comparison between Co and Fe data demonstrates that both impurity species show similar behaviour on completely different energy scales (for instance compare the fourth-layer Fe in Fig. 1b with a fourth-layer Co in Fig. 2).

Figure 1: Variation of the LDOS with the lateral distance of the tip.

a, STM constant current topography (1.8 nm×1.8 nm,2 nA) of a fourth-layer Fe impurity for different bias voltages. For negative bias voltage the LDOS directly above the impurity is reduced. By increasing the voltage the central ring contracts to a plateau-like maximum and the lateral extension of the focusing pattern decreases. b, Section of the differential conductance ΔdI/dV=dI/dV−dI0/dV along the [010] direction. The tip height was adjusted to give a differential conductance for the free surface dI0/dV of 13.3 pA mV−1. An energy-dependent phase shift of the interference pattern can be observed. The vertical black line would resemble the phase front for a energy-independent scattering phase and indicates that the overall phase shift caused by the resonance is less than π. c, Single spectra (purple curve) for four different lateral distances A–D (marked in b by black arrows). To illustrate the relative amplitude of the Kondo signal, the data is normalized to the differential conductance of the free surface (blue scale on the right side). Black curve: calculated spectra obtained by fitting a phenomenological expression of the Kondo resonance to the differential conductance (see Supplementary Information).

Figure 2: STS data of subsurface Co impurities in three–seven layers below the Cu(100) surface.

Spectrum sections ΔdI/dV =dI/dV −dI0/dV along the [010] direction are shown in the upper part. A strong change in the scattering behaviour can be seen for all impurity depths. The black bar corresponds to a length of 0.5 nm for every section. Single ΔdI/dV spectra (purple curve) measured with a tip position direct above the impurity (d=3–7 ML) are depicted in the lower part. The tip height was adjusted to give a differential conductance for the free surface dI0/dV of 3.2 pA mV−1 (3, 4, 6 ML), 5.3 pA mV−1 (5 ML) and 6.4 pA mV−1 (7 ML). To illustrate the relative amplitude of the Kondo signal, the data are normalized to the differential conductance of the free surface (blue scale on the right side). Black curves: calculated spectra obtained by fitting a phenomenological expression of the Kondo resonance to the differential conductance (see Supplementary Information). Note that the colour scale for the spectroscopic data is similar to that used in Fig. 1, except the range varies, and is set by the corresponding single spectra.

For a quantitative analysis, we have to use advanced tools of quantum many-particle theory, as the Kondo problem is a genuine many-body effect. The well-known universal behaviour of the Kondo effect and its associated low-energy fingerprints allows one to apply the single-impurity Anderson model. Within this model the localized orbital of the impurity is described by a single level that couples to non-interacting conduction-band electrons and a Coulomb interaction U between electrons on that level.

The effect of the Kondo resonance in the spectral function can be understood in analogy to other fields of physics, for example scattering of electrons at a potential well: the resonance causes an enhanced scattering amplitude and an energy-dependent phase shift of (quasi-) particles near the resonance energy.

In many-body formalism the measured LDOS is the imaginary part of the single-electron Green’s function G(x,x,ɛ). To calculate this quantity we use the Dyson equation, which connects the Green’s function and hence the LDOS of the perturbed system to the Green’s function G0(x,x′,ɛ) of the unperturbed conduction band electrons by a T-matrix T(ximp,ɛ). Approximating the impurity to be a point scatterer at position ximp, the LDOS change is given by

Calculating the band structure of copper using a linear combination of the atomic orbitals (LCAO) approach (Fig. 3a) and treating the surface as a potential step by fixing the continuity conditions for the wavefunctions we obtain the free Green’s functions G0(x,x′,ɛ) for copper (Fig. 3b and Supplementary Information). For the evaluation of the T-matrix we apply the numerical renormalization groupG0(ximp,ximp,ɛ) at the impurity position is needed to describe the states of the conduction band that couple to the impurity. Note that the presence of the surface leads to oscillations in this quantity (see Fig. 3c), which can be understood as a standing-wave pattern of electrons being reflected at the surface. The NRG calculations result in a Kondo resonance in the T-matrix near ɛF (Fig. 3d), with a resonance height and width showing similar oscillatory behaviour as seen in the unperturbed Green’s functions G0(x,x′,ɛ).

Figure 3: Road map of the theoretical model.

a, Band structure of copper calculated using LCAO and b, the related free propagator. c, Unperturbed LDOS as function of depth, including effects of the surface. d, The Kondo resonance for impurities located at different positions. e, Calculated ΔLDOS section along the [010] direction for Co impurities situated 3–7 ML below the surface. Single spectra taken directly above the impurity position are shown below. Note that the colour scale for the spectroscopic data is similar to that used in Fig. 1, except the range varies, and is set by the corresponding single spectra.

With the T-matrix obtained by the above procedure we calculate the LDOS change using equation (1). Values for energy and hybridization of the localized orbital were taken from ab initio resultsU, was adjusted by comparison to the experiment. The calculated LDOS sections (Fig. 3e) are in excellent agreement with the measured data (compare Fig. 2). The experimentally observed periodicity (compare a third- and seventh-layer impurity) is also found in the simulation.

The rich spatial and spectroscopic information of the measured interference patterns allow further investigation of the Kondo resonance. The resonance width is proportional to the Kondo temperature (see Supplementary Information). We fit the experimental data to a phenomenological form found by Frota (see discussion further down) for all lateral positions and a constant impurity depth. Averaging over all depths, we get a Kondo temperature TK for Fe impurities of 32(16) K and for Co of 655(155) K. Results from macroscopic bulk measurements

Investigating the Kondo temperature as a function of the lateral position for a constant impurity depth we found an unexpected variation of the resonance width. As an example, the fourth-layer Fe impurity (Fig. 1c) shows a variation of TK(A)=29(9) K, TK(B)=35(24) K, TK(C)=42(18) K and TK(D)=44(26) K from position A to D. This behaviour is not included in our theory so far. One possible explanation employs the actual orbital structure of the d-level in Fe or Co, which leads to a dependence of the T-matrix on the propagation direction.

Going further, we investigate not only the long-range signature of the scattering amplitude but also its phase. During our analysis it turned out that the widely-used (complex) Lorentzian approximation of the Kondo resonance does not give the best description to both our experimental and theoretical data (see Supplementary Information). A Lorentzian for the T-matrix results in the often-used Fano line shapesFig. 1b, the phase shift measured in the experiment is also smaller than π.

The mapping of the scattering phase as a function of energy opens a new way to discriminate between single-particle resonances and Kondo scattering. Our approach allows further investigations of the properties of magnetic impurities, such as the splitting of the resonance by an applied magnetic field. This was recently shown on adsorbed atoms