## Isoelectronic tuning of heavy fermion systems: Proposal to synthesize Ce Sb Pd 3 4 3

Jan M. Tomczak 1, ∗

1 Institute of Solid State Physics, Vienna University of Technology, A-1040 Vienna, Austria

The study of (quantum) phase transitions in heavy-fermion compounds relies on a detailed understanding of the microscopic control parameters that induce them. While the influence of external pressure is rather straight forward, atomic substitutions are more involved. Nonetheless, replacing an elemental constituent of a compound with an isovalent atom is-effects of disorder aside-often viewed as merely affecting the lattice constant. Based on this picture of chemical pressure, the unit-cell volume is identified as an empirical proxy for the Kondo coupling. Here instead, we propose an 'orbital scenario' in which the coupling in complex systems can be tuned by isoelectronic substitutions with little or no effect onto cohesive properties. Starting with the Kondo insulator Ce Bi Pt , we consider-within band-theory-isoelectronic substitutions of the pnictogen (Bi 3 4 3 → Sb) and/or the precious metal (Pt → Pd). We show for the isovolume series Ce 3 Bi (Pt 4 1 -x Pd ) x 3 that the Kondo coupling is in fact substantially modified by the different radial extent of the 5 d (Pt) and 4 d (Pd) orbitals, while spin-orbit coupling mediated changes are minute. Combining experimental Kondo temperatures with simulated hybridization functions, we also predict effective masses m ∗ , finding excellent agreement with many-body results for Ce 3 Bi Pt . Our analysis motivates studying 4 3 the so-far unknown Kondo insulator Ce Sb Pd , for which we predict 3 4 3 m /m ∗ band = O (10).

PACS numbers: 71.27.+a, 75.30.Mb, 71.30.+h

## I. INTRODUCTION.

The physics of materials with f -electrons is determined by a delicate interplay of energy scales. The competition between the latter can be tuned by various meansincluding external fields, pressure, carrier dopingresulting in enormously rich phase diagrams 1-4 . Of particular interest are stimuli that control the (Kondo) coupling strength while leaving (nominal) valence states inert. In the case of heavy-fermion metals this allows, for instance, to directly explore Doniach's phase diagram . 5 Isoelectronic tuning can be achieved, e.g., through external or chemical pressure. In the latter case, elemental constituents are (at least partially) replaced with atoms having the same valence configuration but a different principle quantum number n .

In binary compounds, such isoelectronic substitutions directly control the unit-cell volume, as n typically correlates with the element's bonding radius. Furthermore, calculations of the hybridization of the f -states with their surrounding display-in binary compounds-an empirical one-to-one correspondence to the experimental Kondo temperature . Hence the notion that, both, external and 6 chemical pressure manipulate the Kondo coupling J K primarily through the change in lattice constants.

For more complex Kondo systems-ternary compounds and above-we here demonstrate that the unitcell volume (perhaps unsurprisingly) is not a good proxy for the size of the Kondo coupling. Consider that in isoelectronic substitutions both the shape and the radial extent of the wavefunction is different as the principle quantum number changes. The will affect the hybridization and bonding not just with the atoms hosting the 4 f -orbitals but also between all other atoms. Therefore it is not a priori clear whether isovalent substitutions

FIG. 1. Trends in the charge gap. Top: experimental gaps of our focus compounds (see also Fig. 2); the value for Ce Sb Pd 3 4 3 is unknown. Bottom: theoretical gaps within band-theory. Arrows indicate trends for pnictogen (Bi ↔ Sb) and precious metal (Pt ↔ Pd) substitution. We seek to identify the microscopic control parameters that these substitutions engage. Replacing the precious metal was shown 7-9 (is predicted, see Section IV A) to be isovolume for the Bi (Sb) compounds. Note that the gap ∆ N is measured from the Fermi level to the N-point, for reasons explained in the text.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00840v1-figure-f1b53145-2522-4151-b795-f3efbeed14c3.jpg)

modify (i) cohesive properties (largely dominated by hybridizations not involving the f -orbitals), (ii) the Kondo coupling J K (controlled by the hybridization of the f -electrons with their surrounding), or (iii) both of them.

Here, we explore how isoelectronic tuning affects

the itineracy of f -electrons in the family of Kondo insulators 4,10,11 Ce A M 3 4 3 (A=Bi,Sb, M=Pt,Pd), which is of high current interest 12-18 . We consider both possible isovalent substitutions: that of the pnictogen and that of the precious metal. We find that the difference between the Bi- and the Sb-derived materials are mainly driven by the unit-cell volume-a conventional case of isoelectronic tuning. Interchanging the precious metal element in the Bi-compounds, on the other hand, does not change the unit-cell volume 7-9,12 . Nonetheless, we recently established 4 that the Kondo couplings of Ce Bi Pt 3 4 3 and Ce Bi Pd 3 4 3 differ significantly and surmised that the main driver is the different radial extent of the 5d (Pt) and 4d (Pd) orbitals. In the current work, we confirm this proposal and investigate the detailed interplay of volume, orbital extent, hybridization, and charge gaps within band-theory. In other words, we elucidate the microscopic parameters that tune Ce A M 3 4 3 (A=Bi,Sb, M=Pt,Pd) along the arrows in Fig. 1(bottom). Finally, by combining insights from our band-structure calculations with experimental magnetic susceptibilities, we predict the magnitude of mass renormalizations (Tab. I): We find (i) excellent agreement with recent many-body simulations for Ce Bi Pt 3 4 3 14 ; (ii) the Pd-based system have significantly larger renormalizations than the Pt-ones; (iii) the Sb-members have larger Kondo scales than their Bi-counterparts. The latter highly motivates synthesizing and characterizing the so far unknown compound Ce Sb Pd 3 4 3 as well as the series Ce Sb (Pt 3 4 1 -x Pd ) , in which the isovalent substitution x 3 can be studied at more accessible temperatures than for the Bi-series.

## II. SURVEY OF EXPERIMENTAL RESULTS AND PREVIOUS INTERPRETATIONS.

We summarize pertinent experimental findings and their interpretation for the materials under consideration (for a broader perspective see Refs. 4 and 11): As seen in the collection of resistivities (top) and magnetic susceptibilities (bottom) in Fig. 2, the charge and spin gap, ∆ and ∆ , s and the crossover temperature, T max χ , indicative of the onset of Kondo screening, vary significantly in these materials: Interchanging of pnictogen and/or precious metal elements markedly tunes the energy scales. The tendencies of the charge gap ∆-on which we will focus our attention-are condensed in Fig. 1(top): The Pt-based compounds are prime examples of Kondo insulators 7,10,11,19 , with the Sb-member having a slightly larger charge gap than its Bi-counterpart. Exchanging Pt with Pd in Ce Bi Pt 3 4 3 diminishes the charge gap 7,20 , while the volume remains constant 7-9 . The isovolume nature of the precious-metal substitution was recently confirmed by studying the continuous series Ce Bi (Pt 3 4 1 -x Pd ) x 3 12 . The end member Ce Bi Pd 3 4 3 9 has been characterized as a semi-metal 12 , while more recent experiments 17 extract a gap ∆ ∼ 2meV

FIG. 2. Survey of experimental data. Shown are the resistivities from Ref. 20 and Ref. 23 (top panel) and magnetic susceptibilities from Ref. 7 and Ref. 24 (bottom panel) of Ce Bi Pt 3 4 3 and Ce Sb Pt , respectively. The susceptibility of 3 4 3 Ce Bi Pd 3 4 3 is from Ref. 12. Using phenomenological fits as indicated (in grey, see also Ref. 4), the charge and spin gaps are extracted. For χ , the fitting ignores low-temperature CurieWeiss contributions likely to derive from impurities. The gap in Ce Bi Pd 3 4 3 (no data shown) is estimated as ∼ 2meV 17 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00840v1-figure-327f2c5c-8f44-4a36-8f71-0e9b22d3db28.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00840v1-figure-c0f8b709-7592-4327-ac0e-e07cb1338669.jpg)

from the Hall resistivity and the critical magnetic field. Under the assumption of a direct correspondence between unit-cell volume and the Kondo coupling J K , it was initially asserted 12 that J K is independent of x in isovolume Ce Bi (Pt 3 4 1 -x Pd ) . x 3 It was then concluded 12 that the decrease of energy scales with increasing x is driven by the change in the spin-orbit coupling, owing to the mass differences of Pt and Pd-an exciting perspective in view of the recent proposals of topological Kondo insulators 21,22 . Here, instead, we argue in favour of a substantial evolution of the Kondo coupling in the series Ce Bi (Pt 3 4 1 -x Pd ) . The compound Ce Sb Pd x 3 3 4 3 has-to the best of our knowledge-not yet been synthesized.

## III. CRYSTAL-STRUCTURE AND METHODS.

The compounds Ce A M 3 4 3 (A=Bi,Sb, M=Pt,Pd) are a family of cubic intermetallics that crystallize in the body-centred Y Sb Au 3 4 3 structure 25 which has the non-symmorphic space-group I-43d, two formula units per unit-cell and a lattice constant a = 10 051-. 10 058 ˚ . A 7-9,12 ( a = 9 814 ˚ . A 26 ) for Ce Bi Pt 3 4 3 and Ce Bi Pd 3 4 3 (Ce Sb Pt ). 3 4 3 In our band-structure calcu-

FIG. 3. Experimental and theoretical lattice constants . Shown is the relative energy difference ∆ E for varying lattice constant a , using the PBE (symbols and solid lines) and the LDA (dashed lines) potential, where lines are parabolic fits. Arrows indicate experimental lattice constants (see text). The colours are chosen as follows: Ce Bi Pt 3 4 3 l , Ce Bi Pd 3 4 3 l , Ce Sb Pt 3 4 3 n , Ce Sb Pd 3 4 3 n . In that order, the computed equilibrium lattice constants are 10.12, 10.10, 9.90, 9.88 ˚ A with PBE and 9.89, 9.85, 9.70, 9.66 ˚ using LDA. A

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00840v1-figure-763f97d9-4b49-4fa5-84c6-ba6887a97252.jpg)

lations we used a = 10 051 ˚ ( . A a = 9 814 ˚ ) for Ce . A 3 Bi M 4 3 (Ce Sb M ), as motivated below in Section IV A. 3 4 3 In all four compounds, the Ce atoms occupy in the Wyckoff position 12a (3/8, 0, 1/4), Pt/Pd sits at site 12b (7/8, 0, 1/4), and Bi/Sb is found at (u,u,u). Here, we have followed Ref. 27 and used the ideal position u = 1 12, / which is very close to u = 0 0857 found for Ce Bi Pt . 3 4 3 8 and u = 0 084 of Ce Bi Pd . 3 4 3 9 .

In this setting, we performed electronic structure calculations using the wien2k package 28 . We employed the PBE functional, and included spin-orbit coupling (SOC) in the scalar relativistic approximation for all atoms. Volume optimizations were performed with 12 3 reducible kpoints, self-consistency for experimental lattice constants was reached with a k-mesh of 20 3 points. We use the auxiliary Kohn-Sham spectrum as proxy for the electronic structure of the studied materials. Statements about inter-atomic hybridizations are based on: (a) maximally localized Wannier functions 29 constructed in the total angular momentum J -basis for the subspace of Ce-4 f , Ap and Md orbitals with Wannier90 30 via wien2wannier 31 and using a tool of J. Fernandez Afonso 32 to rotate the local spin coordinate system; (b) a local projection formalism due to Haule et al. 33 . The investigation of the electronic structure of the Pt-members has been pioneered by Takegahara et al. 27 . The band-structures and density of states of the two Bi-compounds have been previously shown in the review article Ref. 4.

## A. Comment: Band-theory for Kondo insulators

Given that effective one-particle theories fail in most respects for systems with f -electrons, a few comments are in order. Indeed, the presence of quasi-localized states makes the use of genuine many-body methods a prerequisite for an adequate description, see, e.g., Refs. 14, 16, 33-

41. Nonetheless, conventional electronic structure theory can give meaningful insights into at least the trends of properties in families of f -electron compounds 6,42,43 .

This is particularly true for (Ce-based) Kondo insulators: In the language of the Kondo lattice model (KLM), the charge gap in Kondo insulators is interaction-driven through the Kondo coupling J K . However, this spinfermion interaction is not a fundamental force, but a quantity that emerges in the renormalization group-like flow to a low-energy description. Indeed, it is well known that the KLM can be obtained as a limiting case of the periodic Anderson model (PAM), which contains the Hubbard U interaction as a remnant of the fundamental Coulomb force: The strong U -coupling limit of the PAM yields the weak J K -coupling regime of the KLM with J ∝ V 2 /U 44 . The above illustrates that (bare) hybridizations V between conduction electrons and the f -level and the Hubbard interaction U , together, generate the Kondo coupling. Note that, while the KLM is metallic for vanishing interaction, J = 0, in the absence of Hubbard U -interactions, the (symmetric, half-filled) PAM is a band-insulator. In this sense, the Kondo insulating gap in the PAM is not interaction-driven.

Similar to the non-interacting PAM, band-theory finds Ce A M 3 4 3 to be covalent band-insulators. Band-theory of course grossly overestimates the hybridization V of f -states with their surroundings. Yet, it is sensitive to changes in the Lanthanide environment (external pressure, isoelectronic substitutions) and can therefore qualitatively capture modifications in the essential interatomic hybridizations and the Kondo coupling. Hence, trends in the hybridization gap of Ce 3 A M 4 3 within DFT can be used as a gauge for the Kondo coupling.

Let us note that this situation is facilitated by (i) the nominal 4 f 1 configuration in Ce-based compounds 45-47 , as multiplet effects are less crucial than in systems with more electrons in the f -shell 48 ; (ii) the f -electrons in Ce-materials are more delocalized than, e.g., in the holeanalogues based on trivalent Yb; (iii) contrary to Mott insulators, the self-energy of Kondo insulators is benign near the Fermi level 4 , yielding a renormalized bandstructure below the Kondo temperature.

## IV. RESULTS &amp; DISCUSSION.

## A. Volume optimization.

We assess the isovolume nature of Ce Bi (Pt 3 4 1 -x Pd ) x 3 by optimizing the volume of the end members ( x = 0 1) , within DFT for fixed fractional atomic positions. As can be seen in Fig. 3 the theoretical lattice constants of Ce Bi Pt 3 4 3 and Ce Bi Pd -identified 3 4 3 as the minima in the energy curves-are indeed very similar. We find, using PBE and the Pt-compound as reference: ∆ a/a = 0 2%. Experiments find instead that Ce Bi Pd . 3 4 3 is minutely larger than Ce Bi Pt : 3 4 3 ∆ a/a = -0 07% . 12 . We repeat the same calculations for the pair Ce Sb M 3 4 3

FIG. 4. DFT band-structures. Ce Bi Pt 3 4 3 &amp; Ce Bi Pd 3 4 3 (a) and Ce Sb Pt 3 4 3 &amp; Ce Sb Pd 3 4 3 (b). On the left are indicated the fundamental indirect gap ∆ and the gap ∆ N used as a proxy for the apparent gap in the density of states in Fig. 5(a).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00840v1-figure-1f4f009c-0790-4303-a236-aedbb78630cb.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00840v1-figure-ec7ab575-766e-453d-8ff0-090d212e7b04.jpg)

(M=Pt,Pd). We find that exchanging Bi with Sb yields significantly smaller compounds. Crucially, the exchange Pt ↔ Pd leads again only to very small changes: ∆ a/a = 0 2%. . 49 We interpret the overall small changes in lattice constant, incurred upon precious-metal substitution, to warrant the label 'isovolume'. For Ce Bi M 3 4 3 this supports previous experimental findings 7-9,12 , while the isovolume nature of Ce Sb (Pt 3 4 1 -x Pd ) x 3 is a prediction.

Let us remark that replacing Pt with Pd is not affecting the lattice constant is not a generic behaviour. Indeed, the difference in lattice constant between CePtSn 50 and CePdSn 51 , amounting to ∆ a/a = 0 8%, is experimentally . an order of magnitude larger than in Ce 3 Bi (Pt 4 1 -x Pd ) . x 3 We can speculate that in the more complex unit-cell of Ce A M 3 4 3 (in which, e.g., every 'M' has 4 nearest 'A' neighbours) small atomic differences can be better compensated than in CeMB (where each 'M' has only 3 neighbouring 'B').

We further note that-in line with common experience-PBE(LDA) over(under)estimates the lattice constant. In the following electronic structure calculations we have therefore settled to employing experimental values: we use the lattice constant of Ce Bi Pt 3 4 3 for Ce Bi M , and the one of Ce Sb Pt 3 4 3 3 4 3 for Ce Sb M , where M=Pt,Pd. 3 4 3

## B. Band-structures.

Fig. 4 displays the band-structures of all Ce A M 3 4 3 compounds considered here, while Fig. 5(a) shows the corresponding density of states (DOS). For a given pnictogen 'A', the valence states are rather similar, irrespective of the precious metal 'M'. Changing the pnictogen, results in perceptible changes of the valence states, in particular near the Brillouin zone centre. Overall, the valence states are-congruent with photoemission spectroscopy 52,53 -comparatively dispersive as they are dominantly of pd -character. Indeed, as expected with only one nominal 4 f -electron, quasi-atomic narrow bands are prevalent only above the Fermi level, where also the DOS exhibits pronounced peaks.

In this work, we are mostly interested in the trends of the charge gap under isovalent substitutions. The gaps

∆ extracted from the band-structures are reported in Fig. 5(b). All gaps are indirect and spread, for most compounds, from around the H-point to the zone centre Γ. An exception is Ce Bi Pt , where the top of the va3 4 3 lence bands emerges between the P and the Γ-point. Another speciality of Ce 3 Bi Pt 4 3 is the fact that the Γ-point conduction-band minimum is particularly pronounced. While the DOS-dominated by the f -spectral weightsuggests that the charge gap increases when substituting Pd for Pt, the fundamental gap of Ce 3 Bi Pd -controlled 4 3 by the dispersive band at Γ-is actually larger than in Ce Bi Pt 3 4 3 within DFT. Since we believe the physics of these compounds to be dominated by the f -electrons, we introduce-as a proxy for the trends in the f -band position-the quantity ∆ N which measures the gap between the Fermi level and onset of the conduction states at the N -point. To disentangle putative control parameters of the charge gap, we report in Fig. 5(b) also the values of ∆ and ∆ N when turning off the SOC on selected atoms. Furthermore, so as to demonstrate the effect of the unit-cell volume onto the charge gap, we also show results for a fictitious variant of Ce 3 Bi Pt 4 3 that has the lattice constant of Ce 3 Sb Pt . 4 3

## C. Analysis of electronic structure

In this discussion section, we asses the influence of the spin-orbit coupling, the unit-cell volume, and the radial extent of orbitals onto the charge gap.

## 1. Spin-orbit coupling

As seen in Fig. 5(b), the SOC is-within DFTirrelevant for the qualitative trends of the charge gap in this family of compounds: Turning selectively off the SOC for Pt in Ce Bi Pt -to simulate 3 4 3 the change in atomic mass of the precious metal atoms-does not result in a band gap comparable to isovolume Ce Bi Pd . 3 4 3 In fact, this procedure yields a ∆ and ∆ N that surpass the values obtained with SOC on all atoms. Hence, the difference in masses of Pt and Pd and the resulting change

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00840v1-figure-729fc0a8-3332-4654-9aea-719714d477f8.jpg)

(b)

FIG. 5. Density of states and charge gaps. Panel (a) displays the density of states (DOS) for all four materials. Panel (b) collects the fundamental indirect gap ∆ occurring between the valence band maximum near the H-point and the conduction band minimum at Γ. Also shown are the larger values ∆ N that are measured between the valence maximum and onset of conduction states at the N-point. Owing to the small weight of conduction states at Γ, ∆ N is more representative for the trends of the gap as apparent in the DOS, see also Fig. 4. As illustrated by selectively turning off the spin-orbit coupling (SOC) for the pnictogen and precious metal atoms, it is seen that within DFT the SOC does not control the overall trends within this family of compounds. Additionally shown are gap values of fictitious Ce 3 Bi Pt 4 3 that uses the lattice constant/volume of Ce 3 Sb Pt . 4 3

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00840v1-figure-a2639dfe-a358-4f69-83dc-ff579a297e4a.jpg)

in the magnitude of relativistic effects is likely not the driver of the Kondo insulator to semi-metal transition in Ce Bi (Pt 3 4 1 -x Pd ) x 3 of Ref. 12 (see also the discussion in Ref. 4). Also, eliminating the SOC of Bi in Ce 3 Bi Pt 4 3 to simulate a smaller pnictogen mass when going to Ce Sb Pt -cannot bridge the disparity of band-gaps in 3 4 3 these two materials. The spin-orbit coupling propels the properties neither along the vertical (Bi ↔ Sb) nor the horizontal (Pt ↔ Pd) arrow in Fig. 1.

## 2. Hybridization nature of the gap

Afar more common tuning parameter in heavy fermion systems is the unit-cell volume. While the precious metal substitution series Ce 3 Bi (Pt 4 1 -x Pd ) x 3 is isovolume 7-9,12 , changing the pnictogen atom from Bi to Sb decreases the lattice constant by 2% 7,8,26 . To isolate the effect of the lattice, we report, in Fig. 5(b), the charge gap of fictitious Ce 3 Bi Pt 4 3 where the lattice constant has been artificially shrunken to the value realized in Ce 3 Sb Pt . 4 3 With respect to true Ce Bi Pt , we find a significantly 3 4 3 increased gap. In fact, the change in the unit-cell volume accounts for almost 80% of the gap difference between Ce Sb Pt 3 4 3 and Ce Bi Pt . 3 4 3 This finding identifies the unit-cell volume as the leading protagonist for the trends between the Bi and the Sb pairs of materials.

That the lattice has such a large effect onto the electronic structure finds its origin in the hybridization nature of the gap in these systems. Indeed, as has been demonstrated in Ref. 4 it is the hybridization of the Ce4 f orbitals to their surrounding that accounts for the gap opening: Treating the Cef (J=5/2) as core-states results in metallic ground-states, in analogy to a periodic Anderson model with vanishing hybridization between localized f and conduction states.

Here, we want to further distinguish the individual contributions of atoms/orbitals to the hybridization gap ∆. To this end we constructed maximally localized Wan- nier functions (see method section) and expressed the projected DFT Hamiltonian in that basis. Fig. 6 shows what happens to the relative magnitude of the gap in (a) Ce Bi M 3 4 3 and (b) Ce Sb M 3 4 3 when the hybridization (off-diagonal elements in the Hamiltonian) between the Ce-4 f orbitals and the M -d (M=Pt,Pd) or the Bi/Sb -p orbitals is scaled with a factor α ≤ 1. 54 We find the gap to be most sensitive to changes in the hybridization between the 4 f -orbitals of Ce and the p -orbitals of Bi and Sb. For the Sb compounds, a smaller scaling factor α is needed to collapse the gap, congruent with Ce Sb M 3 4 3 having larger hybridizations than their Bi-counterparts. The gap is much less sensitive to hybridizations between the Cef and the precious metal M- -orbitals. d We note that for both A=Bi and A=Sb the gap is more quickly diminished with α if the precious metal is M=Pt. The relative suppression of ∆( α ) when scaling the precious metal-hybridization is similar irrespective of the type of the pnictogen 'A'. This finding supports the assertion that the anticipated diminishing of the gap in Ce Sb (Pt 3 4 1 -x Pd ) x 3 with x is not volume driven-akin to what happens in Ce Bi (Pt 3 4 1 -x Pd ) . x 3

## 3. Radial extent of orbitals

Now we will disentangle changes in the inter-atomic hybridization mediated by (i) the unit-cell volume and (ii) the radial extent of pertinent orbitals. To this end we and analyse the orbitally ( -)resolved l spread Ω l = 〈 | l R 2 | l 〉 - 〈 | l R l | 〉 2 of the Wannier functions, the sum of which is minimized in the maximally localization procedure 29 . We define the mean extent of a type L of orbitals as the average of the square-root of the l -resolved Wannier spread: R L = 1 /N L ∑ l ∈ L √ Ω . l For example, in case of the Ce-4 f orbitals, i.e. L = 4 f , l runs over the N L = 14 spin-orbitals per Ce-atom. Fig. 7 visualizes the thus measured size of orbitals to scale with nearest-neighbour atomic distances for (a) Ce Bi M , 3 4 3

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00840v1-figure-df457c5d-0000-4098-9b9d-7257c9d1c83d.jpg)

(b)

α

FIG. 6. Hybridization nature of the charge gap. Shown is the relative magnitude of the gap ∆ when scaling the hybridization between the Ce-4 f and the Md (M=Pt, Pd), or Ap (A=Bi,Sb) orbitals in a maximally localized Wannier basis with a factor α ≤ 1 for Ce 3 Bi M 4 3 (a) and Ce Sb M 3 4 3 (b).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00840v1-figure-a241590f-5902-43e1-b7d7-3869a7c647ab.jpg)

(b) Ce Sb M 3 4 3 where M=Pt,Pd. Panel (c) displays the trends of R L for precious metal and pnictogen substitutions relative to the materials' lattice constant a .

substitution, the intuitive correspondence between the size of ( p -)orbitals that mediate bonding and the resulting lattice constant holds to a large extent.

First, we discuss the trends upon interchanging Pt ↔ Pd. Congruent with the precious metals' van der Waals 55 and metallic 56 radii indicated in Fig. 7(a), we find the Pd-4 d orbitals to be smaller than the Pt-5 d ones, irrespective of the pnictogen. The Bi and Sbp orbitals instead are slightly larger when the precious metal is palladium, which in part compensates for the chemical pressure of the smaller Pd. It is thus suggestive to ascribe this expansion of the pnictogen to play a role in the Ce 3 Bi (Pt 4 1 -x Pd ) x 3 series being isovolume. Nonetheless, computing the volume of all 'Wannier spheres' (of radii R L ) in the unit-cell, V Wannier = 4 / 3 π × ( 6 × 14 R f 3 +4 × 6 R p 3 +3 × 10 R p 3 ) , we find for Ce Bi Pt 3 4 3 ( V Wannier = 4340˚ ) A 3 a value still larger by 3% than for Ce Bi Pd 3 4 3 ( V Wannier = 4203˚ ), A 3 see Fig. 7(d). The relative difference is similar for the Sbbased compounds Ce Sb Pt 3 4 3 ( V Wannier = 3843 ˚ ) and A 3 Ce Sb Pd 3 4 3 ( V Wannier = 3732˚ ). A 3 This suggests that the crystal-structure is held together dominantly by the (much larger) pnictogen p -orbitals, while the contribution to bonding from the precious metal d -orbitals is subleading. The small dependence of the hybridization gap on the Ce-4 f to Pt/Pdd hopping (see Fig. 6) supports this statement.

The above is further substantiated when looking at the exchange of the pnictogen: Bi ↔ Sb. There, we find the smaller orbitals of the Sb-based compounds to be in line with the smaller unit-cell volume as compared to their Bicounterparts (lattice constants: 10 051 ˚ . A for Ce Bi M 3 4 3 and 9 814 ˚ . A for Ce Sb Pt , see Section IV A). 3 4 3 Indeed, the Wannier volume (as defined above) decreases by more than 10% when replacing Bi with Sb, see Fig. 7(d). Note, that this trend with (pnictogen) chemical pressure is opposite to the generic behaviour of Wannier functions under external pressure, where the spread increases when the volume shrinks 57,58 . When rescaling V Wannier with the unit-cell volume V 0 = a 3 , see Fig. 7(e), we find the relative differences between the pairs with fixed precious metal 'M' to more than halve. Hence, for the Bi ↔ Sb

Finally, we discuss the extent of the Ce-4 f orbitals. The differences in R f between the four compounds may seem minute on the scale of Fig. 7(a,b,c), yet their trends are revealing, see Fig. 7(f): Indeed, while for a given precious metal, a swapping of the pnictogen does not change the size of the 4 f -shell by much, replacing the precious metal Pt with Pd significantly decreases R f despite leaving the unit-cell volume inert. This suggest that-for the chosen compounds-the impact of (i) the unit-cell volume and (ii) the relative extent of the 4 f -orbitals can be approximately disentangled, as discussed in the next section.

## D. Isoelectronic tuning: Summary &amp; Perspective

Fig. 8 visualizes the central result of this paper: It characterizes the charge gap of our four compounds as a function of their lattice constant a and the relative extent R /a f of their f -orbitals: ∆ N = ∆ N ( a, R f /a ). Quite intuitively, ∆ N increases with shrinking volume and larger relative orbital extent, as both lead to larger orbital overlaps/hybridizations. Most importantly, from this representation it is apparent, that the isoelectronic pnictogen replacement primarily acts on the lattice constant, while the relative extent of the f -shell is well-neigh untouched. The isoelectronic tuning mediated by an exchange of the precious metal, on the other hand, leaves the volume fixed, but expands or shrinks the f -orbitals' reach relative to the surrounding atoms.

We have thus identified a viable choice of essentially independent variables that connect our four materials in parameter space: Fig. 8 mirrors Fig. 1 in a one-to-one fashion (if rotated by 40 ◦ to the right).

Let us put these findings into perspective: It is a wellknown fact that the Kondo coupling in f -electron materials can be dramatically affected by (i) pressure or by (ii) isoelectronic substitutions with atoms of varying size 2 . Here, we demonstrated that these two routes can be linearly independent : Besides the effect of exercising

FIG. 7. The radial extent of orbitals. Panels (a-b) visualize (to scale) the average radial extent (square-root of the Wannier spread) of the Cef , Md (M=Pt,Pd), Ap (A=Bi,Sb) orbitals in relation to inter-atomic nearest-neighbour distances. Also shown in (a) is the Van der Waals radius R V dW of Pt (1.75 ˚ ) and Pd (1.63 ˚ ) A A 55 in solid (dashed) gray, as well as the metallic radii R M of Pt (1.39 ˚ ), Pd (1.37 ˚ ), and Ni (1.24 ˚ ) A A A 56 . Panel (c) display the average radial extent R L relative to the respective lattice constant a . Panel (d) shows the summed volume V Wannier of all Wannier orbital-spheres of radius R L in the unit-cell. Panel (e) shows the same as (c) but divided by the unit-cell volume V 0 . Panel (f) displays the average radial extent R f of the Ce-4 f orbitals; the relative extent R /a f is shown in Fig. 8.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00840v1-figure-ae305c59-4fa1-40dc-bd9e-9d80204e1f79.jpg)

chemical pressure, the change in the principle quantum number n of substituted atoms-here Pt: n = 5, Pd: n = 4-may largely change the charge gap while having little or no effect on cohesive properties. The isoelectronic tuning scenario that we propose extents the empirical volume picture (valid in binary compounds) to include orbital effects at constant volume. In this regard, the Ce A M 3 4 3 system is an ideal case, in that the changes in the lattice constant induced by different 'A' atoms and the relative orbital extent of 'M'atoms form an essentially orthogonal coordinate system. In general, volume and R /a f will provide linear independent yet skewed coordinates. Near orthogonality of the control parameters R /a f and the volume are likely realized for replacements of atoms-like Pt ↔ Pd-that have comparable electronegativities and (empirical, van der Waals, or metallic) atomic radii (see below for a 'non-orthogonal' example). It also stands to reason that an isovolume tuning of inter-atomic hybridizations, and hence the Kondo coupling, is more prevalent in complex materials-ternary compounds and above-and for substitutions of atoms that contribute only subleadingly to the charge gap/bonding (as is the case for the precious metal 'M' in Ce A M , cf. Fig. 6 and the discussion 3 4 3 above).

The proposed orbital scenario will be relevant in other heavy fermion systems. A putative example is the Ce(Ni,Pd)Sn system: CePdSn is an antiferromagnetic Kondo metal with a N´ eel temperature of about 7K 51,59,60 , while CeNiSn is a paramagnetic anisotropic Kondo insulator 61,62 . In the standard Doniach phase diagram , 5 CePdSn falls into the regime dominated by the RudermanKittelKasuyaYosida interaction, whereas CeNiSn realizes a larger Kondo coupling J K that impedes long-range magnetic order. This interpretation is congruent with the 'volume scenario': CeNiSn has a unit-cell volume (263.6 ˚ A ) smaller by 6% than CePdSn 3 (281.7 ˚ A ) 3 60 arguably leading to larger hybridizations and, hence, bigger J K . However, it has been highlighted 2 that there is a significant difference between the isoelectronic substitution series, Ce(Pd 1 -x Ni )Sn, that conx nects the two compounds, and CePdSn under hydrostatic pressure. Indeed, while the N´el temperature is fully supe pressed for x ≥ 0 75 . 63 , the pressure needed to achieve the

same is at least 6GPa 60 , which-using the bulk modulus 55.6GPa of CeNiSn 64 -corresponds to a volume reduction of more than 10%: an external compression of almost twice as much is required to equal the effect of chemical pressure. With the 'orbital scenario' proposed above, this large quantitative discrepancy is readily explained: As suggested by the metallic radii R M of Pd and Ni 56 (both are indicated in Fig. 7(a)), the 3 d -orbitals of Ni will be smaller than the 4 d of Pd. We therefore hypothesize that for growing x in CeNi x Pd 1 -x Sn the shrinking of the unit-cell volume and the decrease of the orbital extent conspire to jointly reduce the Kondo couplingdriving the system faster towards smaller J K -coupling in the Doniach phase diagram. The scenario might also be of relevance in the comparison of CeRu Sn 4 6 65 under external 66,67 and chemical (Sn → Ge) 67 pressure. In fact, it will be interesting to gain systematic insights into larger groups of materials to assess the described physics. Indeed, we believe that R /a f could be a useful quantity in databases for heavy-fermion material properties 6,43 . Its approximation R M /a , using tabulated metallic radii R M , could even be used as descriptor in high-throughput surveys. 68

## E. Anticipated many-body renormalizations

Finally, we transition to a more many-body oriented language. In the spirit of the Anderson impurity model, the embedding of the correlated Ce-4 f orbitals into the solid is often described by a so-called hybridization function ∆( ω ) = ω + µ -H loc ff -( G -1 loc ) ff . Here, µ is the chemical potential, and H loc and G loc are the local projections 33 of the Hamiltonian and the Greensfunction, respectively, while the subscripts f denotes the evaluation within the f -orbital subspace. Indeed, ∆( ω ) is the back-bone of dynamical mean-field theory (DMFT) 69 , where it is promoted to a self-consistent Weiss field. However, already on the level of DFT, the hybridization function can divulge important physical insight, in particular regarding trends among material families 6,43 . Moreover, ∆( ω ) can be used to classify insulators 14 , e.g., distinguishing between Mott and Kondo insulators.

For our four materials, Fig. 9(a) displays the imaginary parts of the DFT hybridization function, while panel (b) shows the corresponding positions ω max and peak amplitudes -Im∆( ω max ). In both cases, the origin of energy corresponds to the top of the valence bands, and we have limited the presentation to the m J = 1 2 and / m J = 5 2 / components of the J = 5 2 multiplet (the / m J = 3 2 / components are the smallest).

Overall, the Bi-compounds exhibit peaks that are sharper and closer to the Fermi level, indicating that these materials are closer to the flat-band limit. In the latter, an f -level hybridizes with a dispersionless conduction band, yielding -Im∆( ω ) = πV 2 δ ω ( + µ -ω max ), where V k = V is the amplitude of the hybridization between the conduction band and the Ce-4 f orbitals

FIG. 8. Control parameters. The figure shows a (linearly interpolated) intensity map of the charge gap ∆ N of Ce Bi Pt 3 4 3 l , Ce Bi Pd 3 4 3 l , Ce Sb Pt 3 4 3 n , and Ce Sb Pd 3 4 3 n as a function of the lattice constant a and the relative extent R /a f of the Ce-4 f orbitals. The substitution of the precious metal Pt ↔ Pd (horizontal arrow) only changes the size of the orbitals, while the pnictogen replacement Bi ↔ Sb (vertical arrow) dominantly engages the unit-cell volume via the lattice constant a .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00840v1-figure-3be24c30-2595-46b9-b366-a0c3372a6e35.jpg)

and glyph[epsilon1] k = ω max the position of the conduction band. For conduction-states glyph[epsilon1] k and couplings V k with finite momentum-dependence, the hybridization strength is spread over larger energies. The latter is the case for the Sb-members of the Ce A M 3 4 3 family. For a given pnictogen 'A', the hybridization function is much more pronounced for Pt than for Pd, as could be anticipated from the size of orbitals (5 d vs. 4 d ) as discussed in the previous section.

Using the DFT hybridization function in conjunction with the experimental magnetic susceptibility, we can estimate the magnitude of many-body renormalizations, see Tab. I. We first assess the Kondo temperature from the spin degrees of freedom using the Kondo lattice expression, T K = 4 T max χ / (2 J + 1) 70 with J = 5 2 / and T max χ from Fig. 2. Note that the estimate for T K of Ce Bi Pd 3 4 3 is of the same order of magnitude of what has been deduced from the experimental specific heat 12 . The Kondo temperature can also be determined from the charge degrees of freedom: k B T K = -( π/ 4) Im∆( ω max ) × Z/ (2 J +1) 70,71 . With the T K from above, we then extract the effective mass enhancement 1 /Z . For Ce Bi Pt , 3 4 3 we find Z ≈ 1 / 10, in excellent agreement with recent realistic DMFT calculations 14 . While for the other compounds, a quantitative theoretical assessment has to await the application of many-body electronic structure methods 69,72 73 , , we believe the fol-

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00840v1-figure-94d5575d-c94b-47ee-a117-a9bf20df2c2a.jpg)

(b)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00840v1-figure-c8680a65-5ed7-4860-8b42-b11cdc6a21f0.jpg)

FIG. 9. Hybridization function ∆( ω ) . Panel (a) shows the imaginary part of the hybridization function ∆( ω ) of all compounds, resolved into their m J components. Panel (b) indicates the trends in the maximal amplitude of the m J = 1 2 / component of ∆( ω ) and the frequency ω max where the latter occurs. The origin of energy corresponds to the valence-band maximum.

|                | exp         |         | theory                        | prediction         | prediction   |
|----------------|-------------|---------|-------------------------------|--------------------|--------------|
| material       | T max χ [K] | T K [K] | 〈- Im∆( ω max ) 〉 m J [meV] | m ∗ /m band = 1 /Z | ∆ [meV]      |
| Ce 3 Bi 4 Pt 3 | 80          | 53      | 400                           | 11                 | 17           |
| Ce 3 Bi 4 Pd 3 | 10          | 7       | 290                           | 66                 | 2.6          |
| Ce 3 Sb 4 Pt 3 | 400         | 265     | 480                           | 2.7                | 90           |
| Ce 3 Sb 4 Pd 3 | ?           | ?       | 300                           | ?                  | ?            |

## V. CONCLUSIONS.

TABLE I. Collection of experimental T max χ (see Fig. 2), estimated Kondo temperatures T K = 4 T max χ / (2 J + 1) with J = 5 2, / m J -averaged peak magnitude of the hybridization function -Im∆( ω max ), and anticipated mass renormalization 1 /Z estimated via k B T K = -( π/ 4) Im∆( ω max ) × Z/ (2 J +1); the gap ∆ is predicted via Z × ∆ DFT .

lowing simple estimates to reflect at least the hierarchy of effective masses in the Ce 3 A M 4 3 family: The largely reduced Kondo temperature when substituting Pd for Pt accounts for Ce Bi Pd 3 4 3 having the largest effective mass (1 /Z = 66), as well as for Ce Sb Pt 3 4 3 being less correlated than Ce Bi Pt . 3 4 3 Naively applying the mass renormalization of Ce 3 Bi Pd 4 3 to its DFT band-gap yields 2.6meV, in good agreement with the experimental estimate of 1 8 . ± 0 5meV . 17 . Doing the same for Ce Sb Pt , 3 4 3 however, significantly overestimates the gap. We can surmise that the Sb-compounds are further from a description in terms of fully localized f -states, which renders the Kondo lattice formula used in the determination of T K less reliable.

Given the trends in the band gaps, hybridization functions, and the Kondo temperature (that itself depends on the two former), we nonetheless gauge the Ce Sb Pd 3 4 3 compound-that we propose experimentalists to synthesize-to exhibit mass enhancements of the order of 10.

We have shown the suppression of coherence and energy scales in isoelectronic and isovolume Ce Bi (Pt 3 4 1 -x Pd ) -previously x 3 attributed to changes of the spin-orbit coupling 12 -to derive from a reduced Kondo coupling. The latter is tunable-even at constant volume-through the principle quantum number of orbitals that host conduction electrons. We have identified the relative extent, R /a f , of the lanthanide's f -orbitals as (indirect) control parameter of the Kondo coupling and demonstrated it to be linearly independent from changes mediated by the unit-cell volume V 0 : In the Ce A M 3 4 3 system, the precious metal Pt ↔ Pd (pnictide Bi ↔ Sb) substitution only changes the f -orbitals' hybridizations through R /a f ( V 0 ). Our 'orbital scenario' generalizes the empirical 'volume scenario' (applicable to binary compounds) to isoelectronic tuning in complex heavyfermion materials. In particular, the two parameter dependency J K = J K ( a, R f /a ) explains why there can be notable differences between external and chemical pressure, cf. the Ce(Ni,Pd)Sn system.

Since Ce Bi Pd 3 4 3 has a very small coherence scale, its ground-state remains controversial 12,17 . We therefore propose to synthesise the new material Ce Sb Pd . 3 4 3 We predict it to have a unit-cell volume very close to that of Ce Sb Pt , and a mass enhancement of the order of 3 4 3 10. Anticipating larger gaps in Ce 3 Sb (Pt 4 1 -x Pd ) , the x 3 series provides an attractive testing ground for the proposed 'orbital scenario' of isoelectronic (and isovolume) tuning of energy and coherence scales.

## VI. ACKNOWLEDGEMENTS

The author is indebted to Juan Fernandez Afonso. Exchanges with S. Dzsaber, G. Eguchi, K. Held, J. Kuneˇ s, S. Paschen, and Q. Si are acknowledged. This work has been supported by the Austrian Science Fund (FWF) through project 'LinReTraCe' P 30213-N36. Some calculations were performed on the Vienna Scientific Cluster (VSC).

- ∗ tomczak.jm@gmail.com
- 1 P. Gegenwart, Q. Si, and F. Steglich, Nat Phys 4 , 186 (2008).
- 2 Q. Si and S. Paschen, physica status solidi (b) 250 , 425 (2013).
- 3 S. Wirth and F. Steglich, Nat. Rev. Materials 1 , 16051 EP (2016).
- 4 J. M. Tomczak, J. Phys.: Condens. Matter (Topical Review) 30 , 183001 (2018).
- 5 S. Doniach, Physica B+C 91 , 231 (1977).
- 6 H. C. Herper, T. Ahmed, J. M. Wills, I. Di Marco, T. Bj¨rkman, D. Iu¸ san, A. V. Balatsky, o and O. Eriksson, Phys. Rev. Materials 1 , 033802 (2017).
- 7 M. F. Hundley, P. C. Canfield, J. D. Thompson, Z. Fisk, and J. M. Lawrence, Phys. Rev. B 42 , 6842 (1990).
- 8 G. H. Kwei, J. M. Lawrence, P. C. Canfield, W. P. Beyermann, J. D. Thompson, Z. Fisk, A. C. Lawson, and J. A. Goldstone, Phys. Rev. B 46 , 8067 (1992).
- 9 W. Hermes, S. Linsinger, R. Mishra, and R. P¨ ottgen, Monatshefte f¨r Chemie - Chemical Monthly u 139 , 1143 (2008).
- 10 Z. Fisk, J. Sarrao, J. Thompson, D. Mandrus, M. Hundley, A. Miglori, B. Bucher, Z. Schlesinger, G. Aeppli, E. Bucher, J. DiTusa, C. Oglesby, H.-R. Ott, P. Canfield, and S. Brown, Physica B: 206 , 798 (1995).
- 11 P. S. Riseborough, Advances in Physics 49 , 257 (2000).
- 12 S. Dzsaber, L. Prochaska, A. Sidorenko, G. Eguchi, R. Svagera, M. Waas, A. Prokofiev, Q. Si, and S. Paschen, Phys. Rev. Lett. 118 , 246601 (2017).
- 13 S. Dzsaber, X. Yan, G. Eguchi, A. Prokofiev, T. Shiroka, P. Blaha, O. Rubel, S. E. Grefe, H.-H. Lai, Q. Si, and S. Paschen, arXiv e-prints , arXiv:1811.02819 (2018), arXiv:1811.02819 [cond-mat.str-el].
- 14 J. M. Tomczak, arXiv e-prints , arXiv:1904.01346 (2019), arXiv:1904.01346 [cond-mat.str-el].
- 15 S. Dzsaber, D. A. Zocco, A. McCollam, F. Weickert, R. McDonald, M. Taupin, X. Yan, A. Prokofiev, L. M. K. Tang, B. Vlaar, L. Stritzinger, M. Jaime, Q. Si, and S. Paschen, arXiv e-prints , arXiv:1906.01182 (2019), arXiv:1906.01182 [cond-mat.str-el].
- 16 C. Cao, G.-X. Zhi, and J.-X. Zhu, arXiv e-prints , arXiv:1904.00675 (2019), arXiv:1904.00675 [cond-mat.strel].
- 17 S. K. Kushwaha, M. K. Chan, J. Park, S. M. Thomas, E. D. Bauer, J. D. Thompson, F. Ronning, P. F. S. Rosa, and N. Harrison, arXiv e-prints , arXiv:1906.01740 (2019), arXiv:1906.01740 [cond-mat.str-el].
- 18 D. J. Campbell, Z. E. Brubaker, C. Roncaioli, P. Saraf, Y. Xiao, P. Chow, C. Kenney-Benson, D. Popov, R. J. Zieve, J. R. Jeffries, and J. Paglione, arXiv e-prints , arXiv:1907.09017 (2019), arXiv:1907.09017 [cond-mat.strel].
- 19 M. Kasaya, K. Katoh, and K. Takegahara, Solid State Communications 78 , 797 (1991).
- 20 K. Katoh and T. Takabatake, Journal of Alloys and Compounds 268 , 22 (1998).
- 21 M. Dzero, J. Xia, V. Galitski, and P. Coleman, Annual Review of Condensed Matter Physics 7 , 249 (2016).
- 22 P.-Y. Chang, O. Erten, and P. Coleman, Nature Physics 13 , 794 EP (2017).
- 23 C. D. W. Jones, K. A. Regan, and F. J. DiSalvo, Phys. Rev. B 58 , 16057 (1998).
- 24 M. Kasaya, K. Katoh, M. Kohgi, T. Osakabe, and N. Sato, Physica B: Condensed Matter 199-200 , 534 (1994).
- 25 A. E. Dwight, Acta Crystallographica Section B 33 , 1579 (1977).
- 26 C. D. W. Jones, K. A. Regan, and F. J. DiSalvo, Phys. Rev. B 60 , 5282 (1999).
- 27 K. Takegahara, H. Harima, Y. Kaneta, and A. Yanase, Journal of the Physical Society of Japan 62 , 2103 (1993).
- 28 P. Blaha, K. Schwarz, G.-K.-H. Madsen, D. Kvasnicka, and J. Luitz, Vienna University of Technology, Austria (2001), iSBN 3-9501031-1-2.
- 29 N. Marzari, A. A. Mostofi, J. R. Yates, I. Souza, and D. Vanderbilt, Rev. Mod. Phys. 84 , 1419 (2012).
- 30 A. A. Mostofi, J. R. Yates, Y.-S. Lee, I. Souza, D. Vanderbilt, and N. Marzari, Computer Physics Communications 178 , 685 (2008).
- 31 J. Kunes, R. Arita, P. Wissgott, A. Toschi, H. Ikeda, and K. Held, Computer Physics Communications 181 , 1888 (2010).
- 32 J. Fernandez Afonso, Excitonic condensation in strongly correlated materials (PhD thesis, TU Wien, 2019).
- 33 K. Haule, C.-H. Yee, and K. Kim, Phys. Rev. B 81 , 195107 (2010).
- 34 K. Held, A. K. McMahan, and R. T. Scalettar, Phys. Rev. Lett. 87 , 276404 (2001).
- 35 B. Amadon, S. Biermann, A. Georges, and F. Aryasetiawan, Phys. Rev. Lett. 96 , 066402 (2006).
- 36 J. M. Tomczak, L. V. Pourovskii, L. Vaugier, A. Georges, and S. Biermann, Proc. Natl. Acad. Sci. USA 110 , 904 (2013).
- 37 L. V. Pourovskii, P. Hansmann, M. Ferrero, and A. Georges, Phys. Rev. Lett. 112 , 106407 (2014).
- 38 A. B. Shick, L. Havela, A. I. Lichtenstein, and M. I. Katsnelson, Scientific Reports 5 , 15429 EP (2015).
- 39 P. Wissgott and K. Held, Eur. Phys. J. B 89 , 5 (2016).
- 40 E. Plekhanov, P. Hasnip, V. Sacksteder, M. Probert, S. J. Clark, K. Refson, and C. Weber, Phys. Rev. B 98 , 075129 (2018).
- 41 S. Seong, K. Kim, E. Lee, C.-J. Kang, T. Nam, B. I. Min, T. Yoshino, T. Takabatake, J. D. Denlinger, and J.-S. Kang, Phys. Rev. B 100 , 035121 (2019).
- 42 G. Zwicknagl, Rep. Prog. Phys. 79 , 124501 (2016).
- 43 H. Hafiz, A. I. Khair, H. Choi, A. Mueen, A. Bansil, S. Eidenbenz, J. Wills, J.-X. Zhu, A. V. Balatsky, and T. Ahmed, npj Computational Materials 4 , 63 (2018).
- 44 Actually, there exists an exact mapping of the PAM to the KLM beyond the Kondo regime 74 .
- 45 J. C. Fuggle, F. U. Hillebrecht, Z. Zoglyph[suppress]nierek, l R. L¨ asser, C. Freiburg, O. Gunnarsson, and K. Sch¨nhammer, Phys. o Rev. B 27 , 7330 (1983).
- 46 F. Strigari, T. Willers, Y. Muro, K. Yutani, T. Takabatake, Z. Hu, Y.-Y. Chin, S. Agrestini, H.-J. Lin, C. T. Chen, A. Tanaka, M. W. Haverkort, L. H. Tjeng, and A. Severing, Phys. Rev. B 86 , 081105 (2012).
- 47 F. Strigari, M. Sundermann, Y. Muro, K. Yutani, T. Takabatake, K.-D. Tsuei, Y. Liao, A. Tanaka, P. Thalmeier, M. Haverkort, L. Tjeng, and A. Severing, Journal of Electron Spectroscopy and Related Phenomena 199 , 56 (2015).

- 48 See, however, e.g., Refs. 37 and 75 for multiplet effects in Ce-compounds under external pressure.
- 49 All relative changes are slightly larger when using LDA instead of PBE.
- 50 T. Takabatake, H. Iwasaki, G. Nakamoto, H. Fujii, H. Nakotte, F. de Boer, and V. Sechovsk, Physica B: Condensed Matter 183 , 108 (1993).
- 51 S. K. Malik, D. T. Adroja, S. K. Dhar, R. Vijayaraghavan, and B. D. Padalia, Phys. Rev. B 40 , 2414 (1989).
- 52 K. Breuer, S. Messerli, D. Purdie, M. Garnier, M. Hengsberger, G. Panaccione, Y. Baer, T. Takahashi, S. Yoshii, M. Kasaya, K. Katoh, and T. Takabatake, EPL (Europhysics Letters) 41 , 565 (1998).
- 53 Y. Takeda, M. Arita, H. Sato, K. Shimada, H. Namatame, M. Taniguchi, K. Katoh, F. Iga, and T. Takabatake, J Electron Spectrosc 101103 , 721 (1999).
- 54 See also Ref. 4 where such a construction helped distinguishing the gap nature in covalent FeSi and ionic LaCoO 3 .
- 55 A. Bondi, The Journal of Physical Chemistry 68 , 441 (1964), https://doi.org/10.1021/j100785a001.
- 56 in Chemistry of the Elements (Second Edition) , edited by N. Greenwood and A. Earnshaw (Butterworth-Heinemann, Oxford, 1997) second edition ed., pp. 1144 - 1172.
- 57 J. M. Tomczak, T. Miyake, R. Sakuma, and F. Aryasetiawan, Phys. Rev. B 79 , 235133 (2009), preprint arXiv:0906.4398.
- 58 J. M. Tomczak, T. Miyake, and F. Aryasetiawan, Phys. Rev. B 81 , 115116 (2010), preprint arXiv:1006.0565.
- 59 M. Kasaya, T. Tani, F. Iga, and T. Kasuya, Journal of Magnetism and Magnetic Materials 76-77 , 278 (1988).
- 60 F. Iga, M. Kasaya, H. Suzuki, Y. Okayama, H. Takahashi, and N. Mori, Physica B: 186 , 419 (1993).
- 61 T. Takabatake, F. Teshima, H. Fujii, S. Nishigori, T. Suzuki, T. Fujita, Y. Yamaguchi, J. Sakurai, and D. Jaccard, Phys. Rev. B 41 , 9607 (1990).
- 62 T. E. Mason, G. Aeppli, A. P. Ramirez, K. N. Clausen, C. Broholm, N. St¨cheli, E. Bucher, u and T. T. M. Palstra, Phys. Rev. Lett. 69 , 490 (1992).
- 63 M. Kasaya, T. Tani, H. Suzuki, K. Ohoyama, and M. Kohgi, Journal of the Physical Society of Japan 60 , 2542 (1991), http://dx.doi.org/10.1143/JPSJ.60.2542.
- 64 K. Ikushima, H. Yasuoka, Y. Uwatoko, and Y. Isikawa, Phys. Rev. B 60 , 14537 (1999).
- 65 I. Das and E. V. Sampathkumaran, Phys. Rev. B 46 , 4250 (1992).
- 66 K. Sengupta, K. K. Iyer, R. Ranganathan, and E. V. Sampathkumaran, Journal of Physics: Conference Series 377 , 012029 (2012).
- 67 J. Zhang, S. Zhang, Z. Chen, M. Lv, H. Zhao, Y. feng Yang, G. Chen, and P. Sun, Chinese Physics B 27 , 097103 (2018).
- 68 For non-cubic compounds, such as CeNiSn and CeRu Sn , 4 6 the lattice parameter a might have to be replaced with relevant nearest neighbour distances.
- 69 A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, Rev. Mod. Phys. 68 , 13 (1996).
- 70 P. Coleman, Phys. Rev. B 28 , 5255 (1983).
- 71 A. C. Hewson, The Kondo Problem to Heavy Fermions (Cambridge University Press, 1993).
- 72 J. M. Tomczak, P. Liu, A. Toschi, G. Kresse, and K. Held, The European Physical Journal Special Topics 226 , 2565 (2017).
- 73 For the temperature dependence of ∆( ω ) in Ce Bi Pt , 3 4 3 see Ref. 14; for spectral properties of Ce Bi Pt 3 4 3 and Ce Bi Pd , see Refs. 14 and 16. 3 4 3
- 74 P. Sinjukow and W. Nolting, Phys. Rev. B 65 , 212303 (2002).
- 75 W. M. Temmerman, A. Svane, L. Petit, M. Lders, P. Strange, and Z. Szotek, Phase Transitions 80 , 415 (2007), http://dx.doi.org/10.1080/01411590701228703.
Jan M. Tomczak
Isoelectronic tuning of heavy fermion systems: Proposal to synthesize Ce3Sb4Pd3
The study of (quantum) phase transitions in heavy-fermion compounds relies on
a detailed understanding of the microscopic control parameters that induce
them. While the influence of external pressure is rather straight forward,
atomic substitutions are more involved. Nonetheless, replacing an elemental
constituent of a compound with an isovalent atom is---effects of disorder
aside---often viewed as merely affecting the lattice constant. Based on this
picture of chemical pressure, the unit-cell volume is identified as an
empirical proxy for the Kondo coupling. Here instead, we propose an "orbital
scenario" in which the coupling in complex systems can be tuned by
isoelectronic substitutions with little or no effect onto cohesive properties.
Starting with the Kondo insulator Ce$_3$Bi$_4$Pt$_3$, we consider---within
band-theory---isoelectronic substitutions of the pnictogen (Bi$\rightarrow$Sb)
and/or the precious metal (Pt$\rightarrow$Pd). We show for the isovolume series
Ce$_3$Bi$_4$(Pt$_{1-x}$Pd$_x$)$_3$ that the Kondo coupling is in fact
substantially modified by the different radial extent of the $5d$ (Pt) and $4d$
(Pd) orbitals, while spin-orbit coupling mediated changes are minute. Combining
experimental Kondo temperatures with simulated hybridization functions, we also
predict effective masses $m^*$, finding excellent agreement with many-body
results for Ce$_3$Bi$_4$Pt$_3$. Our analysis motivates studying the so-far
unknown Kondo insulator Ce$_3$Sb$_4$Pd$_3$, for which we predict
$m^*/m_{band}=\mathcal{O}(10)$.