Magnetic field-tuned Fermi liquid in a Kondo insulator

Kondo insulators are expected to transform into metals under a sufficiently strong magnetic field. The closure of the insulating gap stems from the coupling of a magnetic field to the electron spin, yet the required strength of the magnetic field–typically of order 100 T–means that very little is known about this insulator-metal transition. Here we show that Ce\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{3}$$\end{document}3Bi\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{4}$$\end{document}4Pd\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{3}$$\end{document}3, owing to its fortuitously small gap, provides an ideal Kondo insulator for this investigation. A metallic Fermi liquid state is established above a critical magnetic field of only \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${B}_{{\rm{c}}}\approx$$\end{document}Bc≈ 11 T. A peak in the strength of electronic correlations near \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${B}_{{\rm{c}}}$$\end{document}Bc, which is evident in transport and susceptibility measurements, suggests that Ce\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{3}$$\end{document}3Bi\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{4}$$\end{document}4Pd\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{3}$$\end{document}3 may exhibit quantum criticality analogous to that reported in Kondo insulators under pressure. Metamagnetism and the breakdown of the Kondo coupling are also discussed.

K ondo insulators are a class of quantum materials in which the coupling between conduction electrons and nearly localized f -electrons may lead to properties that are distinct from those of conventional band insulators [1][2][3] . Remarkable properties of current interest include topologically protected surface states that are predicted [4][5][6] and reportedly confirmed by experiments [7][8][9][10][11][12] , and reports of magnetic quantum oscillations originating from the insulating bulk [13][14][15][16] . The very same magnetic field that produces quantum oscillations also couples to the f -electron magnetic moments, driving the Kondo insulator inexorably towards a metallic state 17,18 . The required magnetic fields in excess of 100 T 17,19 have, however, proven to be prohibitive for a complete characterization of the metallic ground state.
Thermodynamic experiments have thus far provided evidence for the presence of electronic correlations in Kondo insulators at high-magnetic fields. Quantum oscillation experiments, for example, have found moderately heavy masses within the insulating phase in strong magnetic fields 15 . Furthermore, heat capacity experiments have shown that the electronic contribution undergoes an abrupt increase with increasing magnetic field 20,21 : in one case 20 the increase occurs within the insulating phase, suggesting the presence of in-gap states 22 , whereas in another, it coincides 21 with the onset of an upturn in the magnetic susceptibility 18,19 and reports of metallic behavior 7,18 . An unambiguous signature of metallic behavior, such as electrical resistivity that increases with increasing temperature at accessible magnetic fields, has yet to be reported.
We show here that insulating Ce 3 Bi 4 Pd 3 , once polished to remove surface contamination, exhibits properties consistent with it being a reduced gap variant of its sister Kondo insulator Ce 3 Bi 4 Pt 3 22,23 , which is also a topological Kondo insulator candidate 4,24,25 . The atypically small magnetic field of B c % 11 T required to overcome the Kondo gap of insulating Ce 3 Bi 4 Pd 3 (see Fig. 1), makes this material 26 ideal for investigating the metallization of a Kondo insulator in strong magnetic fields. We find a magnetic field-induced metallic state exhibiting an electrical resistivity that varies as T 2 at low temperatures (where T is the temperature), thereby revealing a Fermi liquid ground state in the high-magnetic field metallic phase. We also identify a magnetic field-tuned collapse of the Fermi liquid temperature scale T FL near B c , indicating that Ce 3 Bi 4 Pd 3 may exhibit a magnetic fieldtuned quantum critical point analogous to that observed in SmB 6 as a function of pressure [27][28][29][30][31] . The origin of the collapsing Fermi liquid temperature scale in Ce 3 Bi 4 Pd 3 is revealed by susceptibility measurements, which find a peak as a function of magnetic field at B c that can be traced back to a broad maximum in the susceptibility as a function of temperature at T M in weak magnetic fields. This maximum implies the formation of Kondo singlets at a temperature T K ¼ cT M (where 3 < c < 4) 20 , yielding 15 < T K < 20 K. A gradual evolution of the susceptibility from a crossover along the temperature axis to a sharper peak along the magnetic field axis at B c is therefore unveiled.

Results
Evidence for a small gap Kondo insulator. Typical signatures of a Kondo insulating state include an electrical resistivity that exhibits a thermally activated insulating ∂ρ xx =∂T < 0 behavior (where ρ xx is the electrical resistivity) over a wide range of temperatures and a Curie-Weiss behavior in the magnetic susceptibility at high temperatures that crosses over via a maximum in the susceptibility into a Kondo gapped state at low temperatures 1,3,22,[32][33][34][35] . Zero magnetic field electrical resistivity and magnetic susceptibility measurements reveal Ce 3 Bi 4 Pd 3 to exhibit such signatures (see Figs. 1 and 2 and Supplementary Fig. 1), but with a Kondo temperature scale that is significantly reduced relative to those found in the archetypal Kondo insulators SmB 6 and Ce 3 Bi 4 Pt 3 . Typical Kondo insulators have Kondo gaps of several millielectronvolts, accompanied by inverse residual resistivity ratios in the range 10 2 < R À1 RRR < 10 5 and maxima in the magnetic susceptibility and heat capacity 20,36 in the range of temperatures 10 K ( T M ≲ 100 K. Ce 3 Bi 4 Pd 3 , however, appears to have a significantly smaller Kondo gap, accompanied by an inverse residual resistivity ratio of only R À1 RRR $ 5 (see Fig. 2a) and maxima in the susceptibility (see Fig. 2b) and heat capacity (see Fig. 3) at T M % 5 K. A possible explanation for the smaller Kondo gap in Ce 3 Bi 4 Pd 3 compared to its sister compound Ce 3 Bi 4 Pt 3 is the reduction in the strength of the hybridization identified in electronic structure calculations 37 .
The small size of the Kondo gap in Ce 3 Bi 4 Pd 3 in Fig. 2a gives rise to a shallow slope of the resistivity Arrhenius plot in Fig. 2a (see Supplementary Fig. 2 for the magnetic field B-dependence), which therefore has a strong likelihood of being affected by the presence of in-gap states and changes in the transport scattering rate with temperature. To obtain an estimate of the Kondo gap that is less dependent on the scattering rate, we turn to the approximately activated behavior observed in Hall effect measurements 38 (see Fig. 4a), whereupon we obtain Δ ¼ 1.8 ± 0.5 meV (% 21 K). This value is quantitatively consistent with estimates obtained from modeling the heat capacity and susceptibility (see below). An important factor in being able to extract the carrier density is that ρ xy (shown in Fig. 4a) is linear in B at low B. We therefore find that Δ % k B T K , which is quantitatively consistent with the hybridized many-body band    39 -the only significant difference being that Ce 3 Bi 4 Pd 3 has a several times smaller gap and a higher concentration of in-gap states, n Ã 0 (corresponding to 0.2% holes per Ce), than Ce 3 Bi 4 Pt 3 . Consistent with Ce 3 Bi 4 Pd 3 having a several times smaller gap, the insulating state is found to be suppressed by magnetic fields that are several times weaker in strength than those required in Ce 3 Bi 4 Pt 3 (see Fig. 4b).
The atypically small Kondo gap and finite concentration of ingap states in Ce 3 Bi 4 Pd 3 imply that this system lies closer to the insulator-metal threshold 40 than archetypal Kondo insulators, leaving open the possibility that its Kondo insulating behavior is tuned by impurities and defects, as has been found to be the case in CeNiSn 41,42 . Although subsequent samples of CeNiSn prepared with fewer impurities were found to exhibit increased metallic behavior, the opposite has thus far been found in Ce 3 Bi 4 Pd 3 40,43 . Further, in nominally pure CeNiSn, a residual carrier density quantitatively similar to n Ã 0 in the inset to Fig. 4a yields a coherent metallic state with an a-axis resistivity of between 20 and 50 μΩcm or less 42,44 , whereas this same value of residual carrier density in Ce 3 Bi 4 Pd 3 yields an insulating behavior with an electrical resistivity that is 20 to 50 times higher. Indications are therefore that the residual carriers in Ce 3 Bi 4 Pd 3 are incoherent or extrinsic in contrast to CeNiSn (see Methods). Here, we consistently prepare samples of Ce 3 Bi 4 Pd 3 with insulating behavior at B ¼ 0 (two examples of which are shown in Fig. 2) by polishing away surface oxidation (see Methods) and by avoiding superconducting binaries inherent to the flux technique used to grow these materials (see Supplementary  Fig. 3). All samples prepared in such a manner have a thermally activated carrier density at low-magnetic fields (see Fig. 4a) of similar magnitude to that in Ce 3 Bi 4 Pt 3 , and have an electrical resistivity that increases with decreasing T for all T ≲ 250 K (see Fig. 2a), which includes the absence of saturating behavior at low temperatures (see Methods and Supplementary Fig. 4) 45 . The above observations are consistent with the finding made by way of electronic structure calculations performed in the noninteracting limit (in which the f -electrons are regarded as itinerant) that Ce 3 Bi 4 Pt 3 and Ce 3 Bi 4 Pd 3 both similarly exhibit a charge gap 37,46 , with the robustness of the Kondo gap being largely contingent upon the strength of T K . e ÀΔ=k B T ) behavior for 5 ≲ T ≲ 30 K, Ce 3 Bi 4 Pd 3 appears to exhibit such behavior over a range 6 ≲ T ≲ 35 K (with a smaller slope) that lies outside the fitting range, suggesting that Δ cannot be estimated reliably from ρ xx in the case of Ce 3 Bi 4 Pd 3 . The shallow slope and departure from activated behavior at low T is suggestive of a temperature-dependent scattering rate or residual in-gap states associated either with bulk defects, inclusions of other phases or the sample surface (see Methods). Also shown is a schematic of the crystal structure (see Methods). b Susceptibility χ versus T for Ce 3 Bi 4 Pd 3 (at B ¼ 0.1 and 5 T in gray and blue, respectively) compared with that for Ce 3 Bi 4 Pt 3 (at B ¼ 0.1 T) from ref. 35 , with T M indicating the maxima. The inset shows 1=χ versus T for Ce 3 Bi 4 Pd 3 , evidencing Curie-Weiss behavior (indicated by red dashed fitted line) over a broad range of T and a Curie constant C % 0.8 mol Ce per emu K. NATURE COMMUNICATIONS | https://doi.org/10.1038/s41467-019-13421-w ARTICLE from insulating behavior for B < B c to metallic behavior for B > B c . We identify metallic behavior as an electrical resistivity having a value that decreases with decreasing temperature in the limit T ! 0 (i.e., ∂ρ xx =∂T > 0). Further evidence for the transition or crossover from insulating to metallic behavior with the closure of the Kondo insulating gap in Ce 3 Bi 4 Pd 3 is provided by Hall effect measurements at high-magnetic fields. The Hall resistivity ρ xy begins to drop on the approach to B c , which indicates an increase in carrier density in advance of the closing of the Kondo gap, likely caused by thermal excitations of carriers across the reduced gap. The nonlinearity of ρ xy above B c indicates the presence of multiple bands and a possible anomalous Hall contribution from skew scattering 47 (see Methods). In the absence of a dominant anomalous Hall component from skew scattering, the positive Hall coefficient implies either that the hole carriers have numerical supremacy over electron carriers, as would be expected on the basis of band structure calculations 37 , or that the hole carriers are of higher mobility 38 .
Fermi liquid behavior. On entering the high-magnetic field regime, we find the low-temperature electrical resistivity to display the Fermi liquid form ρ xx ¼ ρ 0 þ AT 2 , where ρ 0 is a constant and A is a coefficient proportional to the square of the quasiparticle effective mass. The quadratic-in-temperature form is demonstrated by plotting the electrical resistivity versus T 2 in Fig. 5, yielding linear behavior at the lowest temperatures. The Fermi liquid temperature scale T FL , which we define as the point beyond which the resistivity is observed to exhibit a downward curvature from T 2 behavior on increasing the temperature, is found to collapse on approaching B c in Fig. 1. At the same time, A is observed to undergo a steep upturn on approaching B c from high-magnetic fields in Fig. 1. Both the upturn in A and collapse of T FL are indicative of a collapsing Fermi liquid energy scale in the vicinity of B c , which is one of the primary signatures of quantum criticality in Fermi liquid systems 48,49 . However, rather than peaking precisely at B c , A is observed to peak at B % 15 T. This is a likely consequence of the metallic contribution to the conductivity vanishing due to the collapse of the sizes of the Fermi surface pockets in the limit B ! B c . Other possible factors include sample inhomogeneities causing B c to be non-uniform throughout the sample or an anisotropic gap causing the gap closure field to be momentum-space-dependent.
Magnetic susceptibility measurements. Thermodynamic evidence linking B c to the breakdown of the Kondo coupling accompanying the closing of the gap is provided by measurements of magnetization and magnetic susceptibility in Fig. 6. A strong argument for the reduction in Kondo screening is provided by the approach of the magnetization towards saturation in the inset to Fig. 6a, which is indicative of a state in which the f -electron moments are to a large extent polarized and subject to reduced fluctuations. More directly, we find the insulator-tometal transition or crossover near B c to be accompanied by a lowtemperature peak in the magnetic susceptibility χ in Fig. 6a, which we trace back to the peak in χ at T M along the temperature axis. On performing magnetic susceptibility measurements over a range of temperatures and magnetic fields, we construct a pseudo phase boundary from the locus of the peak that continuously connects B c with T M within the B; T plane (see Fig. 1), revealing that B c ¼ BðT M ! 0Þ. In the low-magnetic field limit, the broadness of the peak is indicative of a crossover, yet it becomes increasingly sharp on increasing the magnetic field and reducing the temperature.
Specific heat measurements. Further thermodynamic evidence for the magnetic field-induced closing of the Kondo gap is provided by the finding of a Schotte-Schotte anomaly in the specific heat in Fig. 3 (see also Supplementary Fig. 5), which undergoes a continuous change in shape in an applied magnetic field. A Schotte-Schotte 50,51 anomaly occurs under circumstances in which a peak in the electronic density-of-states in the shape of a Lorentzian is offset from the chemical potential by a gap (see Supplementary Fig. 6 and Methods). Fits (see Fig. 3a) to the heat capacity measured at B ¼ 0 yield Δ ¼ 1.7 ± 0.2 meV, which is in excellent agreement with Δ estimated from Hall effect  Fig. 4a. c ρ xx versus T for the same sample measured at various steady static magnetic fields (as indicated). d The same ρ xx as in c versus T 2 . e and f Same as for c and d but measured on a different sample (as indicated). Dashed lines indicate the linearity in T 2 while arrows indicate T FL (the temperature above which ρ xx appears to deviate from T 2 behavior). The slopes of the lines versus T 2 provide an estimate of A (plotted in Fig. 1). A small departure from T 2 behavior in the limit T 2 ! 0 is sample-dependent, suggesting it to be of extrinsic origin.  . 4a). On considering B ≠ 0 and fixing all parameters except g eff , we find that the evolution of the Schotte-Schotte anomaly under a magnetic field can be explained solely in terms of Zeeman splitting, where g eff % 2.9 ± 0.3 is an effective g-factor for pseudospins of ± 1 2 . The critical magnetic field B c therefore roughly corresponds to g eff μ B B % Δ % k B T K . As a further consistency check, a similar value of the gap (Δ ¼ 1.5 ± 0.2 meV) qualitatively accounts for the maximum in the magnetic susceptibility as a function of temperature (see Methods and Supplementary Fig. 7) and its shifting to lower temperatures under a magnetic field.

Discussion
The signatures of strongly enhanced electronic correlations in the vicinity of the critical magnetic field B c % 11 T in Figs. 1, 4, 5 and 6, and their association with an insulator-to-metal transition, are strikingly similar to those recently identified in SmB 6 as a function of pressure 28,31 . In SmB 6 , the entry into the metallic phase with increasing pressure can be qualitatively understood in terms of increased tendency for magnetism accompanying a continuous increase in the valence towards a trivalent value [27][28][29][30][31] . In Ce 3 Bi 4 Pd 3 , the polarization of the f -electron moments by a magnetic field causes the magnetic field to behave like an effective negative pressure; the lattice volume invariably increases under a magnetic field in Ce compounds 52,53 , causing a reduction in the valence towards a trivalent value. The relative ease by which f -electron moments are polarized by a magnetic field in Ce 3 Bi 4 Pd 3 (see inset to Fig. 6a) compared to Ce 3 Bi 4 Pt 3 54 suggests that the smaller Kondo temperature and gap in the former is a likely consequence of it lying closer to integer valence than the latter 37 .
While a peak in the magnetic susceptibility at B c has been predicted to occur on closing the gap in a Kondo insulator 55 , along with signatures of quantum criticality, other scenarios are possible at B c , including the case found in itinerant electron metamagnets [56][57][58] . Here, quantum criticality may occur in the absence of broken symmetry due to the tuning of the highest temperature end point of a notional line of first order phase transitions to absolute zero 56 . Our finding of a peak in the susceptibility that becomes increasingly sharp on lowering the temperature suggests that Ce 3 Bi 4 Pd 3 remains in the crossover regime, however, with the end point remaining inaccessible in the limit T ! 0. It has yet to be determined whether this is an intrinsic property of the metal-insulator transition in Ce 3 Bi 4 Pd 3 or whether it is caused by the broadening of the quasiparticle bands by defects and impurities.
An alternative possibility, given the suppression of the Kondo coupling in strong magnetic fields, is that the Kondo coupling vanishes completely at B c 59,60 -in effect causing the f -electrons to behave in a more localized fashion for B > B c . Indeed, one of the predicted consequences of a localization of the f -electrons in Ce 3 Bi 4 Pd 3 (and Ce 3 Bi 4 Pt 3 ) is a metallic state 37 . The metallic state induced by the localization of the f -electrons has also recently been identified by way of dynamical mean field-theory as a strong candidate for a line-node semimetal 46 , in which case B c would correspond to a insulator-to-line-node semimetal transition or, potentially, a topological quantum critical point 61 . The collapse of the Hall resistivity to a small value for B ≳ B c is indeed consistent with the Fermi liquid regime also being semimetallic.
The contiguity between a vanishing carrier density for B < B c and suppressed Kondo coupling for B > B c has the potential to make the nature of the magnetic field-tuned transition somewhat unique compared to that in other d-and f -electron systems. Whether the strong correlations are caused by proximity to a metamagnetic end point or a magnetic field-induced localization of the f -electrons, in both cases the reduction in electronic correlations beyond B c is generally expected to result in a gradual reduction in the effective mass m Ã57,62 , or equivalently the A coefficient 56,63 , with increasing magnetic field. The reduction in A with increasing field in Ce 3 Bi 4 Pd 3 is indeed qualitatively similar to that previously observed in d-and f -electron metamagnets (for example, Sr 3 Ru 2 O 7 and CeRu 2 Si 2 ) 56,62,63 .
Finally, there is the question of whether strong electronic correlations in proximity to a magnetic field-tuned quantum critical point contribute to the growth of the electronic contribution to the heat capacity in the vicinity of the field-induced transitions into the metallic states of the Kondo insulators Ce 3 Bi 4 Pt 3 and YbB 12 20,21 , or to the temperature dependence of the quantum oscillation amplitudes within the insulating phases of SmB 6 and YbB 12 13,15 . Our findings suggest, at the very least, that a growth of electronic interactions at the expense of a declining hybridization needs to be taken into consideration in future more complete models of the insulating state in strong magnetic fields. A further possibility is that the Kondo insulator gradually transforms into an excitonic insulator 64 in advance of metallization. , with vertical bars indicating the maxima in the temperature dependence. Note that at B ¼ 2 T, a shoulder rather than a peak is observed due to a low T contribution from magnetic impurities (or defects) at low B and T (see Methods). A similar low T contribution is observed in Ce 3 Bi 4 Pt 3 in Fig. 2b, although, in that case, it is farther from T M .

Methods
Sample preparation. Ce 3 Bi 4 Pd 3 crystallizes in a cubic structure of space group I 43d with a noncentrosymmetric unit cell that contains four formula units and a lattice parameter of a ¼ 10.05 Å 26 . The crystals are grown by the Bi-flux technique with starting composition Ce:Pd:Bi = 1:1:1.8. The reagents are put in an alumina crucible placed inside an evacuated quartz tube. The quartz tube is then heated to 1050°C at 100°C/h and is kept there for 8 h. The solution is then cooled down to 450°C at 2°C/h. The excess of Bi is removed by spinning the tube in a centrifuge. The crystallographic structure is verified by single-crystal diffraction at room temperature using Mo radiation in a Bruker D8 Venture diffractometer. The upturn in the magnetic susceptibility at low temperature indicates the presence of uncompensated Ce moments and magnetic impurities (such as other rare earths) 1,3,32-34 , which are frequently seen in Kondo insulators. However, inclusions of other crystalline phases such as CePdBi 2 , which has a sharp antiferromagnetic transition at T N % 6 K 65 , and superconducting Bi-Pd binary compounds, which have a variety of different transition temperatures T c , can also be present. The latter include tetragonal PdBi 2 (T c ¼ 4.25 K), BiPd (T c ¼ 3.7 K) and monoclinic PdBi 2 (T c ¼ 1.73 K). The absence of a sharp feature in the heat capacity at % 6 K rules out a significant presence of CePdBi 2 . On fitting the magnetic susceptibility to a Kondo gap plus a low-temperature Curie-Weiss term to model uncompensated Ce moments and isolated magnetic impurities (see simulations of the magnetic susceptibility below), we estimate their concentration to be equivalent to $2% per Ce. In the case of Tb impurities, which are determined from a chemical analysis of Ce to be the leading impurity, the larger moment would imply a significantly lower concentration of 0.1%. In samples that have been polished, no significant negative contribution to the susceptibility is found that would indicate a significant presence of superconducting impurity phases within the bulk, suggesting that superconductivity seen in electrical transport measurements occurs mostly on the surface (see electrical transport experimental details below).
At this early stage of sample growth, the electrical resistivity becomes higher with improved sample quality (see e.g., the difference between early work on these samples 40 , this work and more recent work by another group 43 ).
Electrical transport experimental details. Samples of Ce 3 Bi 4 Pd 3 are cut into rectangular shapes and polished for electrical transport measurements. All resistivity measurements, including those in static and pulsed magnetic fields are performed using the four wire technique. Whereas freshly polished samples (green line in Supplementary Fig. 4) exhibit insulating behavior down to % 2 K, thermally cycled samples exhibit additional superconductivity (red line in Supplementary  Fig. 3). After re-polishing, most of the additional superconductivity is removed (blue line in Supplementary Fig. 3).
Despite the opening of a gap in the energy spectrum, Kondo insulators are prone to in-gap states that can reside either at the surface or the bulk and generally cause the electrical resistivity to acquire a finite value in the limit T ! 0. In a similar manner to SmB 6 45 , Ce 3 Bi 4 Pd 3 is found to have an inverse resistivity ρ À1 xx that varies linearly in T at low temperatures (see Supplementary Fig. 4), suggesting that scattering from defects may be responsible for the departure from activated behavior at the lowest temperatures (at least down to % 2 K) 45 .
To delineate regions of the phase diagram where insulating ∂ρ xx =∂T < 0 and metallic ∂ρ xx =∂T < 0 are observed, the data in Fig. 5 are differentiated (see Supplementary Fig. 1) and the points where ∂ρ xx =∂T ¼ 0 are plotted in Fig. 1.
Magnetic susceptibility experimental details. Magnetic susceptibility measurements in pulsed magnetic fields are performed using an extraction magnetometer 66 . Magnetic susceptibility measurements in static magnetic fields are performed using a Quantum Design vibrating sample magnetometer. Owing to isolated impurities giving rise to an upturn in the susceptibility at low temperatures 1,3,[32][33][34] , and the small energy scale associated with the Kondo gap in Ce 3 Bi 4 Pd 3 , a steady magnetic field of B ≳ 4 T is required to uncover the maximum in the susceptibility. Only a shoulder is observed at lower magnetic fields (see simulations of the magnetic susceptibility below). In a similar manner to Ce 3 Bi 4 Pt 3 54 , the magnetic susceptibility of Ce 3 Bi 4 Pd 3 continues to remain large as the magnetic field is increased, which is suggestive of a band Van Vleck paramagnetism effect 23 .
Heat capacity experimental and modeling details. Heat capacity measurements on Ce 3 Bi 4 Pd 3 , made using a Quantum Design Physical Properties Measuring System, reveal a peak in the overall heat capacity centered at %5 K (see Supplementary Fig. 5). In Ce 3 Bi 4 Pt 3 , by contrast, a peak is observed in the heat capacity centered at %50 K 20,23 . The heat capacity is found to be reproducible in different samples for T ≳ 1.5 K. Features that appear below %1.5 K are sample-dependent, suggesting them to be of extrinsic origin.
A Schotte-Schotte anomaly occurs when a peak in the electronic density-ofstates in the shape of a Lorentzian is offset by an energy gap Δ relative to the chemical potential (see Supplementary Fig. 6a). In the single impurity limit, the width 2W of the Lorentzian is compared to the Kondo temperature. However, it will also include the effect of quasiparticle scattering from nonmagnetic defects and impurities. While earlier it had been assumed that the Kondo gap is symmetric with respect to the chemical potential (at ε ¼ 0) 20 , band structure calculations of Ce 3 Bi 4 Pt 3 and Ce 3 Bi 4 Pd 3 have indicated an asymmetric gap with flat bands situated just below the chemical potential 37 . Flat bands are also situated above the chemical potential, but are considerably further away in energy. Given what is known about the electronic structure, we model the Schotte-Schotte anomaly using an electronic density of states of the form where σg eff μ B B represents Zeeman splitting. Here, σ represents pseudospins of ± 1 2 while g eff is an effective g-factor. The free energy given by F ¼ Àk B T R 1 À1 DðεÞln½1 þ expðÀε=k B TÞdε while the heat capacity given by C p % C v ¼ T∂ 2 F=∂T 2 j v , which we calculate numerically 67 and fit to the experimental data over the temperature range 1.5 ≲ T ≲ 10 K. During fitting, we assume a phonon contribution C ph =T ¼ βT 2 where β ¼ 1 mJ mol -1 Ce -1 K -4 . It is important to note, that on fitting to the Schotte-Schotte model, one cannot distinguish whether Δ is positive (corresponding to a peak below the chemical potential) or whether it is negative (corresponding to a peak above the chemical potential), or whether there are peaks located at both ±Δ as has sometimes been assumed 20 .
We initially fit the data only at B ¼ 0, from which we obtain Δ ¼ 1.7 ± 0.2 meV, which is comparable to the value estimated from Hall effect measurements (see Fig. 4a), and W ¼ 1.1 ± 0.2 meV. For B ≠ 0, we hold D 0 , W, and Δ constant and allow only g eff to vary. Zeeman splitting alone, with g eff ¼ 2.9 ± 0.3, is able to explain the evolution of the Schotte-Schotte peak under a magnetic field. Δ À W % 0.6 meV (%7 K) is of comparable order to the reduced gap estimated from the slope of the Arrhenius plot in the inset to Fig. 2a.
A distinguishing feature of the Schotte-Schotte anomaly that appears to be present in the experimental data (see Fig. 3 and ref. 40 ) is an electronic contribution to the heat capacity of the form (see Supplementary Fig. 5b) in the limit T ! 0, where ðβ 0 À βÞT 2 exceeds the contribution βT 2 from phonons and where γ 0 is the Sommerfeld coefficient characterizing in-gap states (that also depends on the strength of the magnetic field) 67 . Supplementary Fig. 5b shows that Ce 3 Bi 4 Pd 3 exhibits the form given by Eq. (2) in the limit T ! 0, with γ 0 increasing by %30 mJmol À1 K À2 on increasing the magnetic field from 0 to 4 T. A comparable increase was reported in Ce 3 Bi 4 Pt 3 on increasing the magnetic field from 0 to 40 T 20 . The observed consistency with the Schotte-Schotte model therefore includes (i) a peak at %5 T, (ii) a residual electronic contribution γ 0 and (iii) a magnetic fielddependent low-temperature behavior of the form C el =T % γ 0 þ ðβ 0 À βÞT 2 (see Supplementary Fig. 5 and Methods). It remains to be determined to what extent the Lorentzian line-shape (of halfwidth W ¼ 1.1 ± 0.2 meV) is intrinsic (e.g., related to the Kondo coupling and electronic bandwidths) or affected by quasiparticles scattering from defects and impurities.
Simulations of the magnetic susceptibility. As a further consistency check, we show that similar parameters as used in the modeling of the heat capacity also provide a qualitative explanation for the approximate location of the peak in the susceptibility as a function of temperature and its suppression to lower temperatures as the magnetic field is increased. We calculate the susceptibility using (valid for low-magnetic fields) where f 0 FD ðεÞ is the derivative of the Fermi-Dirac distribution function and C is a constant renormalized to the experimental susceptibility at 300 K, the results of which are plotted in Supplementary Fig. 7. Equation (3) does not include the low-temperature Curie-Weiss-like contribution from impurities and defects. Partly for this reason, it overestimates the the initial growth in the susceptibility with increasing field. The model performs poorly once B $ B c , however, which is likely a consequence of electronic correlations not being included. Electrical transport measurements indicate electronic correlations to be important in the vicinity of B c . On comparing the calculated peaks in the susceptibility (B < B c ) with those obtained in Fig. 6b, we estimate Δ ¼ 1.5 ± 0.2 meV and g eff ¼ 2.7 ± 0.2 from the susceptibility.
To obtain a rough estimate of the concentration of uncompensated Ce moments and magnetic impurities, we fit the low-temperature portion of the magnetic susceptibility at 0.1 T (inset to Supplementary Fig. 7) to χ fit ¼ ð1 À f Þχ þ f χ imp , where f is the volume fraction of impurities, χ is the susceptibility given by Eq. (3) and χ imp corresponds to the impurity concentration. We model this by assuming a modified version of χ in which DðεÞ is replaced by a delta function (for simplicity, we assume the total number of states R 1 À1 DðεÞdε and effective g-factor are the same as used in modeling the Kondo gap). We estimate f $ 2%. Anomalous Hall effect estimate. We estimate the anomalous Hall effect contribution from skew scattering within the high-magnetic field metallic phase using ρ A xy ¼ γχρ xx 47 , where γ $ 0.08 KT -1 is a typical value of the coefficient for Ce compounds 47 . At B ¼ 0, the reduced susceptibility isχ % χ peak 5T =C % 0.03 K -1 , but this falls to %0.014 K -1 at 60 T. Using ρ xx % 160 μΩcm at 60 T, we obtain ρ A xy $ 1. 5 10 -8 Ωm at 60 T, which is $10% of ρ xy .

Data availability
The data presented in this manuscript is available from the corresponding author on reasonable request.