Nature Communications

Nat Commun. 10: 796

Divalent EuRh₂Si₂ as a reference for the Luttinger theorem and antiferromagnetism in trivalent heavy-fermion YbRh₂Si₂

M. Güttler¹, A. Generalov², S. I. Fujimori³, K. Kummer⁴, A. Chikina⁵, S. Seiro⁶⁷, S. Danzenbächer¹, Yu. M. Koroteev⁸⁹, E. V. Chulkov⁸¹⁰¹¹¹², M. Radovic⁵, M. Shi⁵, http://orcid.org/0000-0002-2334-8494N. C. Plumb⁵, C. Laubschat¹, J. W. Allen¹³, http://orcid.org/0000-0002-0671-7729C. Krellner¹⁴, http://orcid.org/0000-0002-9694-4174C. Geibel⁷, D. V. Vyalikh

Denis.Vyalikh@dipc.org

¹¹¹⁵

1. 0000 0001 2111 7257grid.4488.0Institut für Festkörper- und Materialphysik, Technische Universität Dresden, D-01062 Dresden, Germany

2. 0000 0001 0930 2361grid.4514.4MAX IV Laboratory, Lund University, Box 118, 22100 Lund, Sweden

3. 0000 0001 0372 1485grid.20256.33Materials Sciences Research Center, Japan Atomic Energy Agency, Sayo, Hyogo, 679-5148 Japan

4. European Synchrotron Radiation Facility, 71 Avenue des Martyrs, 38043 Grenoble, France

5. 0000 0001 1090 7501grid.5991.4Swiss Light Source and Swiss FEL, Paul Scherrer Institute, CH-5232 Villigen-PSI, Switzerland

6. 0000 0000 9972 3583grid.14841.38IFW Dresden, Helmholtzstr. 20, 01069 Dresden, Germany

7. 0000 0004 0491 351Xgrid.419507.eMax-Planck-Institut für Chemische Physik fester Stoffe, 01187 Dresden, Germany

8. 0000 0001 1088 3909grid.77602.34Tomsk State University, Lenina Av., 36, Tomsk, Russia 634050

9. 0000 0001 0094 8940grid.467103.7Institute of Strength Physics and Materials Science, RAS, Tomsk, Russia 634055

10. 0000 0004 1762 5146grid.482265.fCentro de Física de Materiales CFM-MPC and Centro Mixto CSIC-UPV/EHU, 20018 San Sebastián/Donostia, Spain

11. 0000 0004 1768 3100grid.452382.aDonostia International Physics Center (DIPC), 20080 San Sebastian, Spain

12. 0000 0001 2289 6897grid.15447.33Saint Petersburg State University, Saint Petersburg, Russia 198504

13. 0000000086837370grid.214458.eRandall Laboratory, University of Michigan, 450 Church St, Ann Arbor, MI 48109-1040 USA

14. 0000 0004 1936 9721grid.7839.5Kristall- und Materiallabor, Physikalisches Institut, Goethe-Universität Frankfurt, Max-von-Laue Strasse 1, 60438 Frankfurt am Main, Germany

15. 0000 0004 0467 2314grid.424810.bIKERBASQUE, Basque Foundation for Science, 48011 Bilbao, Spain

© The Author(s) 2019

DOI: 10.1038/s41467-019-08688-y

Published online: 15 February 2019

Abstract

Application of the Luttinger theorem to the Kondo lattice YbRh₂Si₂ suggests that its large 4f-derived Fermi surface (FS) in the paramagnetic (PM) regime should be similar in shape and volume to that of the divalent local-moment antiferromagnet (AFM) EuRh₂Si₂ in its PM regime. Here we show by angle-resolved photoemission spectroscopy that paramagnetic EuRh₂Si₂ has a large FS essentially similar to the one seen in YbRh₂Si₂ down to 1 K. In EuRh₂Si₂ the onset of AFM order below 24.5 K induces an extensive fragmentation of the FS due to Brillouin zone folding, intersection and resulting hybridization of the Fermi-surface sheets. Our results on EuRh₂Si₂ indicate that the formation of the AFM state in YbRh₂Si₂ is very likely also connected with similar changes in the FS, which have to be taken into account in the controversial analysis and discussion of anomalies observed at the quantum critical point in this system.

Abstract

Heavy-fermion materials have unusual electronic behavior due to a dual localized-itinerant character of 4f electrons. Here, by studying divalent EuRh₂Si₂, the authors gain insight into the electronic states of the trivalent heavy fermion system YbRh₂Si₂ and show that it experimentally demonstrates Luttinger’s theorem.

Introduction

The determination of Fermi surfaces (FSs) plays a central role in understanding the many-body physics of strongly correlated metals. This is because most of the exciting phenomena in these systems are governed by thermally excited quasiparticles confined to a narrow window around the Fermi energy (EF)¹–⁵. Furthermore, FSs are often the key to discriminating between rivaling theoretical models describing particular correlated metals¹,⁵. For periodic Kondo lattice (KL) materials, where local orbital electrons both carry a magnetic moment and are coupled to conduction bands, the central issue is the extent to which the moment bearing electrons are included in the FS and the way in which this might be accomplished. A fundamental guiding theoretical principle was established by the seminal realization of R. M. Martin⁶,⁷ that the Luttinger FS sum rule⁷ can be applied even to such strongly correlated systems.

The sum rule states that the FS volume is conserved in the presence of interactions so long as the interactions do not produce a phase transition. In this context, the essence of the sum rule is that if the local orbital magnetic moments are quenched, for example by the Kondo effect, then the underlying local orbital electrons must be counted in the FS. If the moments are not quenched, for example in a magnetically ordered state, the underlying electrons are excluded from the FS. For the KL, the sum rule is satisfied by the formation of heavy quasi-particle states at the Fermi energy, the analog of the Kondo resonance for a local orbital impurity site. The behavior along the line of the Fermi energy labeled EFYb in Fig. 1a shows, for the hole situation that will be relevant for this paper, how hybridization of the local orbital quasiparticle into the conduction bands enlarges the hole-like Fermi surface size from kF to k′F to count the local orbital hole. These two situations are now known as the small and large Fermi surfaces, respectively⁸,⁹.

Fig. 1View larger versionFig. 1. Formation of the large Fermi surface in a Kondo lattice. a Schematic hybridization model for Yb case: the localized and renormalized level of the 4f hole (dashed red line) hybridizes with a light hole-like conduction band (blue) and forms a heavy-quasiparticle band (purple) right below the Fermi level, EF. The Fermi vector shifts from kF to k′F. The red sun symbol indicates the region of f-d hybridization. b ARPES band map of YbRh₂Si₂ (compare ref. ⁹). The Brillouin zone cut corresponds to the black line shown in d. The crystal-field-split Yb 4f bands (highlighted in red) hybridize with a bunch of projected bulk bands (shaded in blue) forming a heavy hole-like quasi-particle band at low temperature and shifting the Fermi vector outwards. c, d Schematic representation of a small and large Fermi surface: due to the formation of the Kondo lattice and the Luttinger theorem, the small hole pocket in c transforms into a large hole pocket in d, which incorporates now the additional 4f hole

Even though the initial theoretical realization is now many years in the past, clear experimental demonstrations of the both Luttinger theorem and its workings for KL materials are still very few. Measurements with the de Haas van Alphen (dHvA) technique have provided some observations of both the small and large FS¹⁰. But for angle resolved photoemission spectroscopy (ARPES) the large FS has been observed thus far only in YbRh₂Si₂⁸,⁹, while for another paradigm system CeCoIn₅ the situation remains problematic and controversial¹¹,¹². This paucity is particularly unfortunate because ARPES is the FS technique most evocative for this concept, in that it allows to see directly both the Fermi surface and the underlying dispersing quasi-particles that define it.

Here, we present an ARPES study which reveals both the workings of the Luttinger sum rule and also another very general FS principle pointed out by Müller–Hartmann¹³, that if the single particle self energy for a dispersing quasi-particle is independent of the crystal momentum k, as is highly plausible for the 4f electrons of a rare-earth material, not only the FS size but also its shape is conserved for symmetry preserving interactions. These two concepts together are the reason why the FSs found in density functional theory (DFT) can often be in good agreement with experiment for rare-earth materials.

As explained further below, we compare two isoelectronic and isostructural rare-earth materials, YbRh₂Si₂ (YRS) and EuRh₂Si₂ (ERS), both in the paramagnetic (PM) state, the former where the FS encloses one strongly correlated 4f hole and N weakly correlated conduction band holes, and the latter where the FS encloses N + 1 weakly correlated holes, i.e. the same total hole count for both, implying a large FS for both. The two situations are illustrated in Fig. 1a by the EF lines labeled Yb and Eu, respectively. Our very pretty result is that, both in ARPES and DFT, the FSs of both materials are essentially similar in size and shape, even though the former has the dispersion of a heavy mass and the latter has a light mass.

Results

Strongly correlated YbRh₂Si₂ and weakly correlated EuRh₂Si₂

Heavily studied YRS¹⁴–¹⁷ is a canonical heavy-fermion KL material, with antiferromagnetic (AFM) order below a Néel temperature TN of 70 mK and superconductivity below 2 mK¹⁷. The AFM transition lies very close to a quantum-critical point (QCP) induced by a tiny magnetic field of 60 mT applied along the basal plane, that suppresses the magnetic order. The Yb valence ν is near 3⁺ which gives one hole in the 4f shell. If magnetic moments are Kondo quenched in the low temperature paramagnetic phase just above TN, we would expect this hole to be included in the FS. As will be seen in more detail below, DFT shows that the generic RERh₂Si₂ (RE = rare earth) FS consists of two sheets called the Doughnut (D) and the Jungle-gym (J)⁸–¹⁰,¹⁸. The Jungle-gym forms a largely interconnected Fermi surface sheet across the whole k-space, whereas the Doughnut consists of a square-like, large hole pocket around the Z-point of 3D Brillouin zone (BZ).

The topology of the Doughnut can be used to distinguish the large (itinerant 4f hole) and the small (localized 4f hole) FSs⁸–¹⁰, via the size of its necks which, respectively, do or do not extend into neighboring BZs, as shown schematically in Fig. 1c, d. Our past ARPES work for YRS clearly shows the large FS of Fig. 1d down to 800 mK⁸. The YRS ARPES spectrum of Fig. 1b shows the heavy 4f quasi-particle and the itinerant band states dispersing across EF marked with colors corresponding to the schematic of Fig. 1a. As a calibration, the small FS of Fig. 1c has been conclusively verified by ARPES and DFT for the stable trivalent, local moment AFM, weakly correlated, isoelectronic reference system YbCo₂Si₂⁸,¹⁹, see also Supplementary Note 1.

ERS, the material studied in this paper, is also a weakly correlated, local moment AFM isoelectronic reference system²⁰–²², but with a valency of 2⁺ it is divalent. Thus, in accordance with the Luttinger sum rule, ERS in its paramagnetic phase should possess the large FS, and in the first part of the paper we will show that this is indeed the case. Further, not only can the electronic structure and Fermi surface for the AFM phase of ERS be calculated easily with DFT, but also the transition at TN = 24.5 K can be explored conveniently with ARPES. In the second part of the paper we present the respective DFT and ARPES results for the AFM phase and, finally, we use these data to discuss possible changes in the electronic structure of YRS in its AFM phase below 70 mK and the possible implications for its QCP behavior.

Paramagnetic phase of EuRh2Si2 from DFT and UV-ARPES

The DFT-derived FS of ERS in the PM phase is shown in Fig. 2a, b as a side and top view, respectively. It very closely resembles the one established for YRS in the Kondo PM regime and both the Doughnut and the Jungle-gym can clearly be seen. The center and inset panels of Fig. 2c compare, respectively, DFT and UV-ARPES data for PM ERS, taken with 45 eV photon energy. For this photon energy and the resulting photoelectron kinetic energy, the photoelectron has a very small relaxation time and a short mean free path. In consequence²³ our UV-ARPES data is both surface sensitive and kz broadened so that the bulk electronic structure is mostly k-summed along kz. Therefore the DFT-derived electronic structure is shown projected along kz onto the 2D surface BZ, which is depicted as a dashed-red line in Fig. 2b. High symmetry points of the surface BZ are marked in both panels with overlined symbols. The large Doughnut sheet can be well distinguished and is seen as a homogenous gray area. The projected Jungle-gym is weakly visible in the first BZ, but due to photoemission cross-section effects it is much more pronounced in the second BZ. Those parts that build the border of the projected gap around the M¯-point and are not masked by the Doughnut, are indicated by arrows. In addition to these features predicted by the bulk band-structure calculations, sharp spectral patterns around the Γ¯-point as well as the diamond-shaped feature around the M¯-point result from surface resonances and Shockley surface states. The Shockley state has been discussed in detail in earlier works²⁴–²⁶. Because of their surface origin these features are obviously missing in the calculated bulk bands in Fig. 2, but can be reproduced and identified with the help of slab band-structure calculations and thus be well separated from bulk contributions to the spectra.

Fig. 2View larger versionFig. 2. Fermi surface for the paramagnetic EuRh₂Si₂. Three-dimensional representation of the bulk Fermi surface with the Doughnut and Jungle-gym sheets in a side a and top b views. The red dotted line in b depicts the surface BZ. c ARPES-derived Fermi-surface measured with 45-eV photons combined with the projected bands from bulk band-structure calculations

The projected Fermi surface excluding the surface-related states is in perfect agreement with the calculated 3D Fermi surface and can thus be attributed to the projected bulk band structure in the PM phase. This result, together with the recently presented large Fermi surface for YRS⁸,⁹ also derived by ARPES, shows nicely and clearly the working of the Luttinger theorem for the PM phases of YRS and ERS (see also Supplementary Note 2 and Supplementary Figure 1) and leaves us with a solid starting point for exploring in more detail how the Fermi surface will be modified in the AFM phase of ERS.

For the AFM phase of ERS, an incommensurate magnetic order with a propagation vector Q = (0, 0, τ) with τ ≈ 0.79 has been established experimentally²⁴. If we approximate the propagation vector close to the sample surface with τ = 1, the unit cell in real space is doubled due to the lowered symmetry and the BZ is halved changing from a truncated octahedron to a simple tetrahedron upon entering the AFM phase, as shown in Fig. 3a. Correspondingly, this transition should be accompanied by a back-folding of electron bands into the halved BZ, i.e. the Z-point will be folded onto the Γ-point and the electronic structure at Z and Γ of the PM phase should become equivalent as two Γ points of the AFM phase. Observation of this basic change requires resolution along kz, which precludes using the kz broadened UV-ARPES spectra. Therefore we performed truly bulk sensitive soft X-ray ARPES measurements.

Fig. 3View larger versionFig. 3. PM-to-AFM transition in EuRh₂Si₂ measured by soft X-ray ARPES. a Brillouin zones of body-centered (PM) and simple tetragonal (AFM) EuRh₂Si₂. b Navigation through the 3D k-space. The solid (dashed) lines indicate the boarders of the Brillouin zones along with corresponding high symmetry points in the PM (AFM) phase. The photoemission measurement arcs, probing the crystal momenta close to the Γ- and Z- points are shown in red (hν = 670 eV) and green (hν = 730 eV), respectively. c Fermi surface sections derived from soft X-ray ARPES performed at 40 K (for the PM phase) and at 10 K (for the AFM phase)

Soft X-ray ARPES on EuRh₂Si₂

At such high photon energies, the probing depth is quite high compared to that of UV-ARPES, the probability to excite surface-related states is significantly reduced, and we can study the bulk kz dependence (a complete kx - kz map obtained by soft X-ray ARPES is presented in Supplementary Figure 2 and discussed in Supplementary Note 3). The energy resolution is reduced but nonetheless adequate for the purpose at hand. We used photon energies of 670 eV and 730 eV, which probe the electronic structure along the measurement arcs shown in Fig. 3b. These provide planes of approximately constant kz near the Γ- and Z- point of the paramagnetic BZ. In Fig. 3c, we compare ARPES-derived Fermi surface maps obtained in the PM phase at T = 40 K and the AFM phase at T = 10 K, see also Supplementary Figure 3 and Supplementary Note 4. As expected, in the PM phase we see distinctly different spectra in the regions of the inequivalent Γ- and Z- points. In particular, a hole-like pocket can be seen around Z, whereas no such pocket exists around Γ. Upon cooling the system down to 10 K below the Néel temperature, a similar hole-like band as that around Z now appears around the Γ-point, where the intensity was previously depleted. At the same time, the hole pocket around Z seems to be unaffected by the AFM transition. The appearance of spectral weight around Γ in the AFM phase, which in the PM phase could be seen exclusively around Z, is fully consistent with the change of the BZ and the resulting back-folding of electronic bands.

Antiferromagnetic phase of EuRh₂Si₂ from DFT and UV-ARPES

Having experimentally established the back folding of electron bands along the Γ - Z direction across the AFM transition, we turn now to the results of our band structure calculations to analyze in more detail the changes of the electronic structure upon entering the AFM phase. The DFT-derived FS is shown in Fig. 4a, b as a side and top view, respectively. The Doughnut and parts of the Jungle-gym are folded into the new BZ. The Doughnut is folded onto the Γ-point along with the previous Z-point of the PM phase. The new Z-plane, which lies halfway between the Γ- and Z-points of the PM structure, cuts the folded Doughnut slightly, resulting in small oval hole pockets labeled (3) in Fig. 4b and marked in green (see Supplementary Figure 4). These pockets are oriented towards the corners of the simple-tetragonal BZ and are a further important indicator for the folded FS in the AFM phase. The large Jungle-gym then hybridizes with the Doughnut, splitting up into two distinct sheets. The resulting splitting of the FS is labeled (2) in Fig. 4b and marked in cyan. The outermost sheet enclosing the BZ corners is labeled (1) in Fig. 4b and marked in red. This banana-shaped feature results from the folding of the Jungle-gym onto itself, splitting off due to the exchange interaction to form long tubes with almost parallel surfaces along kz. Similarly, the outer of the two sheets of feature (2) also forms long, but warped, tubes along kz. As a general conclusion, the bulk DFT calculation predicts a strong fragmentation of the PM FS upon the AFM transition.

Fig. 4View larger versionFig. 4. Fermi surface of AFM-ordered EuRh₂Si₂ derived by UV-ARPES and DFT calculations. a 3D view and b top view of the calculated FS along kz. c Comparison of the projected bulk FS (top left panel) and the ARPES-derived FS (bottom right panel) taken at 10 K. The yellow dashed line depicts the border of the Brillouin zone. d Additionally to the projected bulk FS, the surface states and surface resonances have been extracted from slab calculations for Eu-termination (cyan, top left panel) and Si-termination (green, bottom left). The spin polarization of the termination-specific surface-related states is indicated in red and blue on the right side

We now compare the band structure calculation to the FS experimentally obtained by UV-ARPES at a temperature of 10 K, sufficiently below the Néel temperature. As with the PM phase we must account for the kz broadening of bulk features and the sensitivity to surface features. We consider the bulk features first and show in the upper left part of Fig. 4c the calculated bulk FS projected along kz. As labeled, it still shows nicely the same three key features already called out in Fig. 4b. The gap labeled (2) lies within the projected Doughnut close to the borders. The ARPES-derived FS in the lower right half of Fig. 4c shows these same three bulk features, as indicated by dashed lines with the same color coding. Although the ARPES intensity of the PM phase Jungle-gym is quite low in the first BZ, the borders of the outer FS sheet, which coincide with the banana-shaped feature (1), are clearly visible. The ARPES-derived FS also shows the gap labeled (2) and additionally there are four oval gap-like features which appear exactly where the oval hole pockets (3) are expected.

Considering next the surface, by comparing the ARPES data to slab DFT calculations we find that all remaining features can be understood as a superposition of spectra from both Si and Eu terminations present together on our measured surfaces. In our slab DFT calculations, a supercell consisting of a thick slab of four conventional unit cells (32 atomic layers) surrounded by empty space mimics the two surfaces of the material with several bulk layers in between and enables us to study the surface modifications of the band structure. In Fig. 4d, we overlaid the projected bulk bands (shaded in gray) with the surface states and surface resonances of the FS obtained from our slab calculation (colored). To identify the surface character, only states with a predominant contribution in the topmost four atomic layers are extracted from the overall slab band structure. The two segments at the top colored cyan and red/blue are calculated for the Eu-terminated half of the slab, the two segments at the bottom colored green and red/blue are obtained from the Si-terminated part. The cyan and green states indicate surface states and resonances for the respective terminations. The red and blue coloring on the right side additionally points out the spin polarization of these states.

We have already identified in the PM phase UV-ARPES spectra a pronounced well known Shockley state around the M¯-point. We see that this strictly two-dimensional Shockley state splits off from the border of the projected bulk band gap, occurs only for Si-termination, and is strongly spin split due to exchange interactions with the ferromagnetically polarized Eu 4f moments, both in the calculated AFM FS and in the AFM ARPES data. This characteristic fingerprint thus clearly indicates the presence of Si-terminations in UV-ARPES.

For Eu-termination there are a bunch of surface resonances that arise from bulk states, are pinned to the topmost Eu surface layer, and are further spin split. These account nicely for what is perhaps the most striking and instructive feature in the UV-ARPES derived FS, a dense bunch of bright and rather sharp states, which appear within the Doughnut following roughly its shape. Multiply repeated experiments convincingly indicate reproducibility of this spectral pattern. These states obviously cannot be observed in the PM phase (compare Fig. 2c) and are not present in the calculated bulk FS of the AFM phase.

Discussion

What are the relevance and the implications of our observations on the AFM state of EuRh₂Si₂ for the system YbRh₂Si₂? In the debate on the nature of the field induced quantum critical point in YbRh₂Si₂, the evolution of the Fermi surface from the magnetically non-ordered state to the AFM state is a pivotal issue¹⁵,¹⁶,²⁷. Most of the discussion has focused on the possibility of a change in the nature of the quasiparticle, while the effect of the magnetic ordering on its own has not been considered. However, the present observations on EuRh₂Si₂ show that magnetic ordering results in large changes in the Fermi surface through band folding and band splitting. These results cannot be taken over directly to YbRh₂Si₂ because the AFM structure of YbRh₂Si₂ is yet unknown. However, results on both static magnetic properties (e.g. magnetization) as well as dynamic magnetic properties (e.g. relaxation behavior) (see Supplementary Note 5) indicate that magnetic correlations in YbRh₂Si₂ are very similar to those in EuRh₂Si₂. In particular both display a strong competition between FM and AFM order. With the Fermi surfaces being so similar in the PM states of EuRh₂Si₂ and YbRh₂Si₂, and also a strong similarity in their magnetic correlations, EuRh₂Si₂ becomes a relevant model system for gauging the impact of AFM order on the Fermi surface of YbRh₂Si₂.

One concern with this approach might be that EuRh₂Si₂ is an RKKY mediated large local moment system whereas the extremely small moment of YbRh₂Si₂ strongly suggests that it should be regarded as a spin-density-wave (SDW) system (here SDW refers to its original sense as the ordering of delocalized electrons). However, for the relevance and resulting change of the Fermi surface, both cases are actually very similar. Because in RKKY exchange the transfer of information on the polarization from one 4f site to the next 4f site is exclusively due to conduction electrons, the generalized susceptibility χ(q) of the conduction electrons, and therefore e.g. nesting properties, is as much relevant for local moment ordering as for the SDW scenario. Accordingly, modifications of the band structure, at least on a qualitative level, should be very similar in both kinds of systems. Both our experimental ARPES results and our calculations show that for this Fermi surface, AFM order leads to significant changes, because of both the large amount of band folding and the large band splittings. Such large changes cannot be ignored in discussing the anomalies seen in a number of transport properties at the field induced quantum critical point (QCP) in YbRh₂Si₂²⁸, of special significance because these anomalies are at the heart of a debate on the nature of this QCP¹⁶,²⁷. For example, the Hall effect has received considerable attention¹⁶,²⁷ and it surely must react to large Fermi surface changes. Thus, our results on EuRh₂Si₂ indicate that the formation of the AFM state in YbRh₂Si₂ very likely results in strong changes in the large Fermi surface. Within this picture, the Fermi surface changes would be due to the magnetic order and not to a “large to small” transition arising from Kondo breakdown. These changes could well be sufficient for understanding the anomalous properties observed at the QCP on entering the AFM state. Folding and splitting of the bands is a priori expected to affect the thermodynamic and transport properties only within the AFM state, but not the critical behavior in the paramagnetic regime above TN. However, within this scenario the unconventional critical behavior observed in the paramagnetic regime can e.g. be accounted for by interaction between fluctuations of the antiferromagnetic order parameter and critically renormalized electron quasiparticles²⁷. More experimental studies of the properties of YbRh₂Si₂ in the ordered phase and corresponding microscopic theory making use of the down-folded Fermi surface structure are needed to further explore this electronic structure model. But it is a testable hypothesis.

In summary, we have studied the Fermi surfaces for both the PM and AFM phases of EuRh₂Si₂ using soft X-ray and UV-ARPES combined with ab initio full relativistic DFT band structure calculations. In the PM state we identify the bulk Fermi surface in EuRh₂Si₂ (with divalent Eu) as being very similar to that of the Kondo system YbRh₂Si₂ (with nearly trivalent Yb). This is a beautiful new experimental visualization of the working of the Luttinger sum rule. For PM Kondo systems it predicts that the f degree of freedom, in the case of trivalent Yb the 4f hole, despite its localized character, must be counted in the Fermi volume, leading to the same Fermi volume as for a homolog compound with a divalent rare earth. Furthermore, our ARPES data and our calculations demonstrate that the onset of AFM order leads to pronounced changes in the Fermi surfaces due to band folding and band splitting. The strong changes between the PM and AFM states calculated and observed for EuRh₂Si₂ indicate that the formation of the AFM state in YbRh₂Si₂ very likely also results in strong changes in its large Fermi surface. Therefore, such kinds of changes, as an alternative to changes in the quasi-particle due to Kondo breakdown, have to be considered in the controversial analysis and discussion of the origin of the transport properties observed on crossing the QCP into the ordered state in this system.

As a final note, an important instructive result of our work is the clear demonstration of different sets of spin-polarized states persisting on two distinct surface terminations of EuRh₂Si₂. The Si terminated surface reveals mainly itinerant magnetism confined to the surface due to the presence of Shockley surface states, which couple via exchange interaction with the Eu 4f moments in the fourth layer below the Si surface. However, the Eu-terminated surface exhibits mainly bulk electron states, which show strong resonant behavior at the Eu surface. Due to the exchange coupling with the topmost polarized Eu 4f moments these states become themselves spin-polarized in the vicinity of the surface. These results clearly show that a magnetically ordered hetero system may reveal curious and profoundly different magnetism at its surface depending on the constituting sorts of surface atoms which is challenging to predict and can be interesting for technological applications.

Methods

Experiment

UV-ARPES and Soft X-ray ARPES experiments were carried out at the Swiss Light Source (SIS X09LA instrument) and SPring-8 synchrotron radiation facility (Beamline BL23SU), respectively. The UV-ARPES experiments were discussed in details in²⁴, while the Soft X-ray ARPES experiments were comprehensively presented in ref. ²⁹. The single crystals of EuRh₂Si₂ were grown in indium flux as described in²⁰,²¹.

Theory

The full-relativistic band structure calculations have been performed with the FPLO code³⁰, a density functional theory code based on the full-potential local orbital method, within the generalized gradient approximation (GGA)³¹. The Eu 4f basis states have been fixed to the experimentally found occupation (n4f ≈ 7.0³², frozen core approximation) using either an unpolarized configuration in the paramagnetic phase or a local-moment scenario of 7 μB per Eu atom with ferromagnetic coupling within the Eu ab-planes and antiferromagnetic interaction along the c-axis in the AFM phase. Surface effects were simulated by an asymmetric slab with 32 atomic layers terminated on one side by Si and on the other side by Eu atoms. Further details on the computations can be found in²⁴,³³.

Supplementary Material

Open In Web Browser

Footnotes

Journal peer review information: Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work.

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Supplementary Information accompanies this paper at 10.1038/s41467-019-08688-y.

Acknowledgements

This work was supported by the German Research Foundation (DFG) (through grant Nos. GE 602/4-1, KR3831/5-1 and Fermi-NEst; GRK1621 and SFB1143), Saint-Petersburg State University under Research Grant 15.61.202.2015, Spanish Ministry of Science and Innovation (Grant No. FIS2016-75862-P), and Tomsk State University Academic D.I. Mendeleev Fund Program (Grant No.8.1.01.2017). The SX-ARPES experiment was performed under Proposal Numbers 2016A3810 and 2018B3811 at SPring-8 BL23SU. The present work was financially supported by JSPS KAKENHI Grant Numbers 16H01084 (J-Physics) and 18K03553.

Author contributions

C.G. and D.V.V. designed research; S.S., C.K., and C.G. grew and characterized the samples; M.G., A.G., A.C., S.D., and D.V.V. performed the UV-ARPES experiments; S.I.F. performed Soft X-ray ARPES experiment; Theoretical studies were performed by M.G., Yu.M.K. and E.V.C.; Operation of the UV-ARPES facility was carried out by M.R., M.S., and N.C.P; The obtained results were discussed together with K.K., C.L., and J.W.A.; The manuscript was written by M.G., C.G., J.W.A. and D.V.V.

Data availability

All relevant data are available from the corresponding author upon request.

Notes

Competing interests

The authors declare no competing interests.

Articles from Nature Communications are provided here courtesy of Nature Publishing Group

PMC Copyright Notice

The articles available from the PMC site are protected by copyright, even though access is free. Copyright is held by the respective authors or publishers who provide these articles to PMC. Users of PMC are responsible for complying with the terms and conditions defined by the copyright holder.

Users should assume that standard copyright protection applies to articles in PMC, unless an article contains an explicit license statement that gives a user additional reuse or redistribution rights. PMC does not allow automated/bulk downloading of articles that have standard copyright protection.

See the copyright notice on the PMC site, https://www.ncbi.nlm.nih.gov/pmc/about/copyright/, for further details and specific exceptions.

References

Lu D, et al., authors. Angle-resolved photoemission studies of quantum materials. Annu. Rev. Condens. Matter Phys. 2012;3:129–167. [CrossRef]

Damascelli A, Hussain Z, Shen ZX, authors. Angle-resolved photoemission studies of the cuprate superconductors. Rev. Mod. Phys. 2003;75:473–541. [CrossRef]

Stewart GR, author. Non-Fermi-liquid behavior in d- and f -electron metals. Rev. Mod. Phys. 2001;73:797–855. [CrossRef]

Stewart GR, author. Addendum: non-Fermi-liquid behavior in d- and f- electron metals. Rev. Mod. Phys. 2006;78:743–753. [CrossRef]

Kordyuk AA, author. ARPES experiment in fermiology of quasi-2D metals. Low Temp. Phys. 2014;40:286–296. [CrossRef]

Martin RM, author. The Fermi surface and Fermi liquid properties of periodic Kondo and mixed valence systems. J. Appl. Phys. 1982;53:2134–2136. [CrossRef]

Martin RM, author. Fermi-surface sum rule and its consequences for periodic kondo and mixed-valence systems. Phys. Rev. Lett. 1982;48:362–365. [CrossRef]

Kummer K, et al., authors. Temperature-independent fermi surface in the kondo lattice YbRh₂Si₂. Phys. Rev. X. 2015;5:011028

Danzenbächer S, et al., authors. Insight into the f- derived fermi surface of the heavy-fermion compound YbRh₂Si₂. Phys. Rev. Lett. 2011;107:267601 [CrossRef] [PubMed]

Rourke PMC, et al., authors. Magnetic-field dependence of the YbRh₂Si₂ Fermi surface. Phys. Rev. Lett. 2008;101:237205 [CrossRef] [PubMed]

Jang, S. et al. Evolution of the Kondo lattice electronic structure above the transport coherence temperature. Preprint at https://arxiv.org/abs/1704.08247 (2017).

Chen QY, et al., authors. Direct observation of how the heavy-fermion state develops in CeCoIn₅. Phys. Rev. B. 2017;96:045107 [CrossRef]

Müller-Hartmann E, author. The Hubbard model at high dimensions: some exact results and weak coupling theory. Z. für Phys. B Condens. Matter. 1989;76:211–217. [CrossRef]

Trovarelli O, et al., authors. YbRh₂Si₂: pronounced non-fermi-liquid effects above a low-lying magnetic phase transition. Phys. Rev. Lett. 2000;85:626–629. [CrossRef] [PubMed]

Custers J, et al., authors. The break-up of heavy electrons at a quantum critical point. Nature. 2003;424:524–527. [CrossRef] [PubMed]

Gegenwart P, Si Q, Steglich F, authors. Quantum criticality in heavy-fermion metals. Nat. Phys. 2008;4:186–197. [CrossRef]

Schuberth E, et al., authors. Emergence of superconductivity in the canonical heavy-electron metal YbRh₂Si₂. Science. 2016;351:485–488. [CrossRef] [PubMed]

Zwicknagl G, author. Field-induced suppression of the heavy-fermion state in YbRh₂Si₂. J. Phys. Condens. Matter. 2011;23:094215 [CrossRef] [PubMed]

Güttler M, et al., authors. Tracing the localization of 4f electrons: angle-resolved photoemis¬sion on YbCo2Si2, the stable trivalent counterpart of the heavy-fermion YbRh₂Si₂. Phys. Rev. B. 2014;90:195138 [CrossRef]

Seiro S, Geibel C, authors. From stable divalent to valence-fluctuating behaviour in Eu(Rh(1-x) Irx)₂Si₂ single crystals. J. Phys. Condens. Matter. 2011;23:375601 [CrossRef] [PubMed]

Seiro S, Geibel C, authors. Complex and strongly anisotropic magnetism in the pure spin system EuRh₂Si₂. J. Phys. Condens. Matter. 2014;26:046002 [CrossRef] [PubMed]

Hossain Z, Trovarelli O, Geibel C, Steglich F, authors. Complex magnetic order in EuRh₂Si₂. J. Alloy Compd. 2001;323-324:396–399. [CrossRef]

Speer NJ, et al., authors. Surface vs. bulk electronic structure of silver determined by photoemission. Europhys. Lett. 2009;88:67004 [CrossRef]

Chikina A, et al., authors. Strong ferromagnetism at the surface of an antiferromagnet caused by buried magnetic moments. Nat. Commun. 2014;5:3171 [CrossRef] [PubMed]

Güttler M, et al., authors. Robust and tunable itinerant ferromagnetism at the silicon surface of the antiferromagnet GdRh₂Si₂. Sci. Rep. 2016;6:24254 [CrossRef] [PubMed]

Generalov A, et al., authors. Spin orientation of two-dimensional electrons driven by temperature-tunable competition of spin-orbit and exchange-magnetic interactions. Nano Lett. 2017;17:811–820. [CrossRef] [PubMed]

Wolfle P, Abrahams E, authors. Spin-flip scattering of critical quasiparticles and the phase diagram of YbRh₂Si₂. Phys. Rev. B. 2015;92:155111 [CrossRef]

Paschen S, et al., authors. Hall-effect evolution across a heavy-fermion quantum critical point. Nature. 2004;432:881–885. [CrossRef] [PubMed]

Saitoh Y, et al., authors. First results from the actinide science beamline BL23SU at SPring-8. Nucl. Instrum. Methods Phys. Res. A. 2001;474:253–258. [CrossRef]

Köpernik K, Eschrig H, authors. Full-potential nonorthogonal local-orbital minimum-basis band-structure scheme. Phys. Rev. B. 1999;59:1743–1757. [CrossRef]

Perdew JP, Burke K, Ernzerhof M, authors. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996;77:3865–3868. [CrossRef] [PubMed]

Chevalier B, Coey JMD, Lloret B, Etourneau J, authors. EuIr₂Si₂: a new intermediate valence compound. J. Phys. C Solid State Phys. 1986;19:4521–4528. [CrossRef]

Höppner M, et al., authors. Interplay of Dirac fermions and heavy quasiparticles in solids. Nat. Commun. 2013;4:1646 [CrossRef] [PubMed]

[Back]

Fig. 1.

Formation of the large Fermi surface in a Kondo lattice. a Schematic hybridization model for Yb case: the localized and renormalized level of the 4f hole (dashed red line) hybridizes with a light hole-like conduction band (blue) and forms a heavy-quasiparticle band (purple) right below the Fermi level, EF. The Fermi vector shifts from kF to k′F. The red sun symbol indicates the region of f-d hybridization. b ARPES band map of YbRh₂Si₂ (compare ref. ⁹). The Brillouin zone cut corresponds to the black line shown in d. The crystal-field-split Yb 4f bands (highlighted in red) hybridize with a bunch of projected bulk bands (shaded in blue) forming a heavy hole-like quasi-particle band at low temperature and shifting the Fermi vector outwards. c, d Schematic representation of a small and large Fermi surface: due to the formation of the Kondo lattice and the Luttinger theorem, the small hole pocket in c transforms into a large hole pocket in d, which incorporates now the additional 4f hole

[Back]

[Back]

Fig. 2.

Fermi surface for the paramagnetic EuRh₂Si₂. Three-dimensional representation of the bulk Fermi surface with the Doughnut and Jungle-gym sheets in a side a and top b views. The red dotted line in b depicts the surface BZ. c ARPES-derived Fermi-surface measured with 45-eV photons combined with the projected bands from bulk band-structure calculations

[Back]

[Back]

Fig. 3.

PM-to-AFM transition in EuRh₂Si₂ measured by soft X-ray ARPES. a Brillouin zones of body-centered (PM) and simple tetragonal (AFM) EuRh₂Si₂. b Navigation through the 3D k-space. The solid (dashed) lines indicate the boarders of the Brillouin zones along with corresponding high symmetry points in the PM (AFM) phase. The photoemission measurement arcs, probing the crystal momenta close to the Γ- and Z- points are shown in red (hν = 670 eV) and green (hν = 730 eV), respectively. c Fermi surface sections derived from soft X-ray ARPES performed at 40 K (for the PM phase) and at 10 K (for the AFM phase)

[Back]

[Back]

Fig. 4.

Fermi surface of AFM-ordered EuRh₂Si₂ derived by UV-ARPES and DFT calculations. a 3D view and b top view of the calculated FS along kz. c Comparison of the projected bulk FS (top left panel) and the ARPES-derived FS (bottom right panel) taken at 10 K. The yellow dashed line depicts the border of the Brillouin zone. d Additionally to the projected bulk FS, the surface states and surface resonances have been extracted from slab calculations for Eu-termination (cyan, top left panel) and Si-termination (green, bottom left). The spin polarization of the termination-specific surface-related states is indicated in red and blue on the right side

[Back]