CeRu4Sn6: a strongly correlated material with nontrivial topology.

Topological insulators form a novel state of matter that provides new opportunities to create unique quantum phenomena. While the materials used so far are based on semiconductors, recent theoretical studies predict that also strongly correlated systems can show non-trivial topological properties, thereby allowing even the emergence of surface phenomena that are not possible with topological band insulators. From a practical point of view, it is also expected that strong correlations will reduce the disturbing impact of defects or impurities, and at the same increase the Fermi velocities of the topological surface states. The challenge is now to discover such correlated materials. Here, using advanced x-ray spectroscopies in combination with band structure calculations, we infer that CeRu4Sn6 is a strongly correlated material with non-trivial topology.

A metallic surface is one of the necessary characteristics of a 3-dimensional topological insulator. The surface charge carriers have a massless Dirac dispersion and are protected from backscattering by time reversal symmetry123. So far the materials used in the research field of topological insulators are based exclusively on simple semiconductors and only very recently it was realized that also strongly correlated systems can be topologically non-trivial. Here the Kondo or intermediate-valent insulator SmB6 4 came immediately into consideration567891011 and first investigations121314151617 show that the robust metallicity of the sample surface could be a promising concept for explaining its long-time mysterious low temperature residual conductivity18. Increasing the interest for correlated topological insulators is the possibility to find phenomena that might be even more exotic19202122 than the ones of topological semiconductors.

The important issue is now to find other strongly correlated materials with non-trivial topology. This is not a trivial task since quantitative numerical tools often applied to characterize the topological nature of semiconductors are lacking for correlated materials, due to the fact that the reliability and accuracy of mean field solutions of the electronic structure in the presence of strong correlations are far from sufficient. It is therefore important that recent theoretical studies have drafted effective models for strongly correlated systems that should produce topological non-trivial behavior5610. One of the propositions is that rare earth Kondo or mixed-valent insulators should have all the necessary ingredients for realizing topological insulating properties5678910: a fully spin-compensated insulating ground state formed by orbitals having opposite parity, e.g. 4f vs. 5d, and inverted in energy in particular sections in the Brillouin zone. It is even specifically mentioned that for Ce systems with cubic56 and tetragonal symmetry78 a gap always opens when the ground state is a nodeless J = ±1/2 Kramers doublet and when the 4f occupation is sufficiently away from integer.

Here CeRu4Sn6 comes into play. It is a tetragonal, non-centrosymmetric 23, potential Kondo insulator with a J = 5/2 Hund’s rule ground state. A hybridization-induced gap was reported first in 199224, supported further by several experimental studies252627, and later estimated to be of the order of 100 K from single crystal resistivity data2628. Optical conductivity measurements point towards remnant metallicity29 and yet the low temperature resistivity is in the milliOhm-cm range, indicating clearly the depletion of the charge carriers24262829. Magnetization and resistivity are highly anisotropic strongly suggesting the J = 5/2 Hund’s rule ground state is split by the tetragonal crystal-electric field (CEF)2628. μSR measurements did not find magnetic order down to 50 mK, thus pointing out the importance of Kondo-type interactions in this material30.

One of the pressing open questions is whether CeRu4Sn6 fulfills the conditions mentioned in the theoretical studies about strongly correlated topological insulators5678910. To answer this question, we determine here whether the hybridization between conduction and 4f electrons (cf-hybridization) in CeRu4Sn6 is strong enough and whether the CEF ground state has a favorable symmetry to develop a gap of sizable magnitude. The degree of hybridization has been investigated by performing a temperature dependent x-ray absorption spectroscopy (XAS) experiment at the Ce L-edge to obtain the Ce valence and the Kondo temperature T. The CEF ground state symmetry was studied by measuring the Ce M-edge in an XAS experiment and the Ce N-edge in a non-resonant inelastic x-ray scattering (NIXS) experiment. We have also performed band structure calculations to understand better our experimental findings and to test the robustness of our conclusions. Details about the experimental set-ups and the single crystal preparation and quality can be found in the Section Methods.

Results

Cerium 4f valence and Kondo temperature

Trivalent cerium (Ce3+) has a 4f occupancy n = 1 (f 1). The presence of strong cf-hybridization leads to a partial delocalization of the 4f electrons into the conduction electron sea and the resulting non-integer valent 4f ground state can be written in a local picture as a mixed state , i.e. with additional contributions of tetravalent (f 0) and divalent (f 2) states31. Here and denote the number of ligand holes. For Ce the resulting occupation of the 4f shell is then given by , whereas quantifies the degree of delocalization. The hybridization is described in terms of the formation of a Kondo singlet state at kT below the energy of the 4f 1 configuration, which breaks down as thermal fluctuations match the same energy. Hence, information about the Kondo temperature T is obtained from the temperature variation of n, or equivalently of if or if does not change with T 323334.

Core level spectroscopy is a powerful tool for determining the occupations of the various 4f configurations31. This is based on the fact that the impact of the core-hole potential depends strongly on the configuration. Moreover, recent XAS experiments measuring the recombination process with partial fluorescence yield (PFY-XAS) at the Ce L3-edge have tremendously improved the accuracy of determining valencies343536. We have performed such a PFY-XAS experiment on CeRu4Sn6, determined and obtained T from its temperature dependence.

Figure 1(a) shows the PFY-XAS data of CeRu4Sn6 at 5 K (open circles). The strongest absorption is due to the f 1 contribution in the ground state, the small but well resolved hump at about 5740 eV is due to the f 0 contribution, and the small shoulder on the low energy tail of the f 1 peak (just above 5720 eV) is due to some f 2 spectral weight. It is not straightforward to describe the L3-edge spectral line shape since it is determined by the unoccupied 5d density of states in the presence of a 2p core hole. Yet, for extracting the f weights, it is sufficient to fit the data with a set of Gaussian or similar functions as a mean to obtain their integrated intensities373839. We have applied such a fit by describing each spectral weight (I(f 0), I(f 1), and I(f 2)) with identical line shapes, consisting of three Gaussian profiles each. Only the total intensities and central energy positions were varied. Each profile includes a tanh-type step function (black lines) with the maximum slope at the peak maximum to fit the background. The background was adjusted so that the modulations due to the extended absorption fine structure (EXAFS) at ≈5749, ≈5754, and ≈5759 eV are averaged out. The red line in Fig. 1(a) is the result of such a fit.

Figure 1

(a) PFY-XAS spectra of the Ce L-edge at T = 5 K (open circles). The size of the circles indicates the error bars. The red line is the total fit, the regions of different gray refer to the I(f 2), I(f 1) and I(f 0) spectral weights, taking into account a tanh-type background (black solid lines). (b) Temperature evolution of the Ce L3-edge spectra and the I(f 0) peak region on an enlarged scale. (c) Inset with as function of log(T). The relative error bars amount to ±0.25%. The black, dashed line is a guide to the eye, the red arrow indicates its inflection point which is taken as Kondo temperature T.

For the spectral weights we obtain I(f 0) = 0.059 ± 0.015, I(f 1) = 0.931 ± 0.015 and I(f 2) = 0.011 ± 0.015 which corresponds to and under the assumption . This is an approximation since it neglects hybridization effects in the final states39. As in core level photoemission40, I(f 2) overestimates and I(f 0) underestimates so that the above values for n and should be taken as upper and lower limits, respectively. It is therefore safe to conclude that the cf-hybridization is very strong in CeRu4Sn6.

Figure 1(b) shows the temperature dependence of the PFY-XAS data, normalized to the integrated intensity. The inset in Fig. 1(b) shows the f 0 region on an enlarged scale after smoothing the data, and it is clearly seen that the main changes occur within the temperature interval of 150 to 250 K. Applying the same fit by varying only the total intensity of the f 1 and f 0 contributions and adjusting the height of the background according to the f1 absorption line, yields the temperature variation [upper inset, Fig. 1(c)]. The dashed line is a broadened step function describing the temperature evolution. We determine T empirically as the inflection point of 323334, yielding T = (170 ± 25) K supporting the analysis of the temperature dependence of the single crystal resistivity in ref. 29.

CEF ground state symmetry

The local point symmetry at the Ce site of D2 23 in CeRu4Sn6 gives rise to a tetragonal CEF that splits the J = 5/2 Hund’s rule ground state into three Kramers’ doublets: one of the type , and two mixed states of type  and , in the J-representation . We have investigated the CEF ground state wave function with two spectroscopic methods: soft XAS at the Ce M-edge 414243 and hard x-ray NIXS at the Ce N-edge 44. We did not use inelastic neutron scattering since Kondo insulators exhibit spin gaps and strongly broadened magnetic excitations45464748, complicating the extraction of magnetic spectral weights which are required for the determination of wave functions. There are some neutron data on CeRu4Sn6 which also exhibit a spin gap-type response that peaks at 30 meV and broadened magnetic excitations above it in the high energy data plus some quasielastic intensity in the low energy data49. In contrast, XAS and NIXS target the ground state symmetry directly: In the soft XAS experiment, the dipole selection rules for linearly polarized light give access to the initial state symmetry. Excited CEF states only contribute to the polarization as they get increasingly populated with rising temperature. A NIXS experiment with hard x-rays is truly bulk sensitive, and when performed at large momentum transfers the signal is dominated by higher multipole orders of the scattering cross section5051. In the NIXS experiment the -dependence (-direction with respect to the sample orientation) replaces the polarization dependence in XAS and, accordingly, -dependent multipole selection rules give access to the ground state symmetry44.

Figure 2 shows the soft XAS data at the Ce M4,5 edge for the two polarizations and for several temperatures between 4 and 300 K. At all temperatures there is a strong polarization effect which decreases only slightly between 200 and 300 K, thus indicating the first excited CEF state is situated at ≥25 meV. The data at 4 K will therefore yield ground state information. In addition to the main absorption line, there are smaller humps at slightly higher photon energies (see arrows) due to some f 0 weight in the initial state. Its spectral weight decreases with rising temperature in agreement with the findings from the PFY-XAS data.

Figure 2

Temperature dependence of the linearly polarized M-edge XAS spectra of CeRu4Sn6.

The arrows mark the spectral weight due to the f 0 contribution in the ground state.

Figure 3(a) shows the linearly polarized low-temperature data (red and blue dots). The spectra are simulated ionically with a full muliplet code (see Methods, Analysis) to account for the f 1 ground state symmetry. Transitions due to f 0 in the ground state or transitions to continuum states are not accounted for. Figure 3(a) shows a simulation assuming a pure ground state (red and blue lines), which is the state with rotational symmetry [see inset in Fig. 3(a)]. The simulation reproduces the main features of the spectra and their polarization dependence, but the match is by no means perfect. This issue is investigated further in the following.

Figure 3

(a) Linearly polarized XAS data (dots) at the Ce M-edges for T = 4 K and simulation with a ground state (lines). Blue corresponds to the polarization and red to . (b) Experimental (dots) and simulated (colored lines) linear dichroism (LD), , of the pure J states of the J = 5/2 multiplet. (c) Same as (b) for the J = 7/2 multiplet. The inset in (a) shows the spatial distribution of the 4f electrons for Γ6 symmetry.

In Fig. 3(b) the experimental (LD) (, black dots) is compared with the simulated LD of the three pure states (colored lines). Note here the f 0 and continuum contributions have no or negligible polarization dependence. The LD of the pure (orange) and also of the (blue), corresponding to a Γ7-symmetry with α = 0, have the correct sign, but the LD of the latter is too small and any -admixture (green) in a Γ7-type ground state would make the discrepancy even larger. Even if we adopt a Γ7-type ground state of predominantly , the measured LD will be larger than the simulated one. We therefore conclude that although the ground state is not a pure , it is still mainly of that character.

A NIXS experiment at the Ce N-edge was performed in order to rule out surface effects being responsible for the unsatisfactory agreement between the soft XAS data and the simulation. The data in Fig. 4 (red and blue dots) were obtained after 240 s counting per point and after subtracting a rising, linear background (Compton). The colors are chosen as before, i.e. red for and blue for . Again there is a clear dichroic effect. We compare the measured N4,5 edge with simulations for the pure states (red and blue lines), which here are shown again (cf. ref. 44) for reasons of convenience. Details about the theoretical simulations containing a full multiplet ansatz can be found in Section Methods, Analysis.

Figure 4

Dots: Experimental NIXS intensity at the Ce N4,5 edge after subtracting a linear background.

The point size reflects the statistical error, the counting time was 240 s per point. Lines: Calculated scattering function for the pure J states of the J = 5/2 multiplet. Data and simulation are shown for (blue) and (red).

The strongest dichroism in the data appears at the shoulder structures at about 105 and 113 eV. Qualitatively the simulation for the state reproduces the data quite well. In contrast, the state suggests a crossing of both directions on both shoulders, which is not observed experimentally, and the dichroism of the state has the wrong sign. We therefore conclude also from this experiment that the ground state must be governed by a symmetry, but at the same time is not purely given by the Γ6.

Discussion

What could be the reason for the deviation between the experiment and the simulation starting from a pure ? Already in 1986 van der Laan et al. pointed out that in the presence of strong hybridization higher lying states contribute to the M-edge signal52. In CeRu4Sn6 we find that the Kondo temperature (170 K) is not much smaller than the CEF excitations so that a contribution of higher lying CEF states to the ground state is indeed to be expected. Given the fact that the LD of these higher lying CEF states is smaller or even of opposite sign, see Fig. 3(b,c), it is evident that the LD of this composite ground state will be reduced from that of the pure state. We could also reverse the argumentation: the observation of a reduced LD can be taken as evidence for the presence of strong hybridization. Nevertheless, we can conclude that the ground state of CeRu4Sn6 is predominantly of symmetry. This is in good agreement with the anisotropy of the low-temperature magnetization28 and susceptibility53, thereby supporting also the results of LDA + DMFT calculations at 290 K54.

We now will use density functional theory based band structure calculations to obtain a deeper understanding of the experimental findings (details in Section Materials and Methods, Band structure calculations). In Fig. 5(a) we display our result of a full-relativistic non-magnetic calculation employing the local density approximation (LDA) for CeRu4Sn6. We can clearly observe that the band structure has a true band gap in agreement with an earlier study in ref. 29, comprising of unoccupied Ce 4f bands above the Fermi level and occupied Ru 4d bands below the Fermi level. The presence of the gap is important since it confirms that the system has the even number of electrons per unit cell necessary to allow for the formation of a fully spin-compensated, non-magnetic state. We also notice in the left panel of Fig. 5(a) that the lowest unoccupied Ce 4f state is surprisingly given by the (blue dots) and not by the (red dots). Yet, looking at the Ce 4f partial density of states depicted in the right panel of Fig. 5(a) we find that it is the which has the largest occupation, as found in our experiment. We therefore can infer from these calculations that the empty has indeed the stronger hybridization with the occupied Ru 4d bands so that it gets pushed up higher in energy and, at the same time, acquires a larger occupation.

Figure 5

(a) Full relativistic (LDA + spin-orbit) band structure of CeRu4Sn6, wherein the Ce 4f weight is mostly above the Fermi level. From the J projected DOS, we note the largest contribution to the occupied valence sector is from J = 5/2 and J = ±1/2. (b) Applying a constant potential to the J = 1/2 state. (c) Applying a constant potential to the J = 3/2 state.

Next, we will make an attempt to estimate the effect of correlations on the relevant symmetries of the system. We do this by introducing artificially an extra potential for a particular J state so that it becomes more occupied, i.e. as to mimic the Ce 4f1 valence state thereby making the inversion with the Ru 4d. Figure 5(b) shows the result for moving down the orbital and Fig. 5(c) the . We can observe that in the case the system becomes metallic with quite a large number of bands crossing the Fermi level55. For the the case, we also see that the system is metallic, but only with a few crossings and small Fermi surface areas, suggesting that the number of charge carriers will be low. Important for the topology is that the band structure is that of a semi-metal, i.e. that there is a warped gap, so that a unique assignment can be made about the character of the bands. This aspect did not necessarily follow from the conditions set by Dzero and Galitzki78, since there could always have been some other non-bonding bands crossing the Fermi level, which is not considered in their effective model. Yet, our calculations indicate specifically that such is not the case for CeRu4Sn6. It is also important to notice that the warped gap is substantial, i.e. 20 meV or larger, which corroborates nicely with the fact that the J = ±1/2 is nodeless and therefore guarantees the presence of hybridization with the Ru 4d bands at all k-points of the Brillouin Zone (except at the time reversal inversion points).

To summarize, we thus have found that CeRu4Sn6 fulfills the conditions mentioned by Dzero and Galitzki78 to be a strong topological insulator. It has a ground state symmetry given primarily by the and has a 4f occupation sufficiently far away from integer with a very high Kondo temperature (170 K) so that the effective degeneracy is larger than two. We have in addition verified that there are no non-bonding bands that otherwise could spoil the formation of a warped gap. We therefore conclude that CeRu4Sn6 is a strongly correlated material with non-trivial topology.

Methods

Samples

All experiments were performed on single crystals. Single crystals were grown with the self-flux floating-zone melting method using optical heating in a four-mirror furnace56. The quality of the crystals was checked with Laue and x-ray diffraction as well as scanning electron microscopy and energy dispersive x-ray spectroscopy, which show that the samples are single grain and single phase materials. The lattice parameters are in agreement with a = 6.8810 Å and c = 9.7520 Å as given in ref. 23. The single crystals were pre-aligned with a Laue camera. A peculiarity of the present tetragonal structure is that within ≈0.2%. As a consequence the structure appears pseudo-cubic in Laue images and the final alignment of the single crystals had to be done via magnetization measurements26.

Set-ups of x-ray techniques

The PFY-XAS experiment was performed at the GALAXIES inelastic scattering endstation at the SOLEIL synchrotron in France57. The synchrotron radiation was monochromatized using a Si(111) nitrogen-cooled fixed-exit double-crystal monochromator (ΔE/E ≈ 1.4 × 10−4), followed by a Pd-coated spherical collimating mirror. The x-rays were then focused to a spot size of 30 μm × 90 μm FWHM [vertical x horizontal] at the sample position by a 3:1 toroid Pd-coated mirror. The overall energy resolution for the experiment was measured to be 700 meV FWHM with about 1 × 1013 photons/second on the sample. The incident x-ray energy was scanned across the Ce L3 absorption edge (≈5723 eV) while monitoring the intensity of the Ce de-excitation (≈4839 eV) using a Ge 331 spherical bent crystal analyzer (Bragg angle 80.8°, 1 m radius of curvature) in the Rowland circle geometry which focused the emitted x-rays onto a silicon avalanche photodiode. All measurements took place in a backscattering geometry and normalization was performed using a monitor placed just before the sample that recorded the scattered radiation from a Kapton foil placed in the beam path. The flight path was filled with helium gas to reduce the air scattering and the samples were cooled in a helium cryostat.

The XAS spectra were recorded at the synchrotron light source BESSY II using the UE46 PGM-1 undulator beamline. The total electron yield (TEY) was measured under ultra high vacuum of the order of 10−10 mbar. The CeRu4Sn6 single crystals were cleaved in situ to guarantee clean sample surfaces. The TEY signal was normalized to the incoming photon flux I0 as measured at the refocusing mirror. The energy resolution at the cerium M4,5 edges was set to 0.15 eV. The undulator combined with a normal incident measurement geometry allow a change in polarization without changing the probing spot on the sample surface. The crystals were oriented such that the two polarizations and , [001] being the long tetragonal axis, could be measured. Several crystals were measured or re-cleaved in order to assure the reproducibility of the data.

The NIXS measurements were performed at the beamline ID20 at ESRF. A single Si(111) monochromator set the incident energy to 9690 eV and the scattered intensity was analyzed by three columns of three Si(660) analyzers (in total nine) at in-plane scattering angles of 2Θ = 140°, 146°, and 153°, corresponding to an averaged momentum transfer of . This configuration yields an instrumental energy resolution of FWHM = 1.3 eV and made the Ce N4,5 edge measurable beyond the dipole limit445051. Two crystals with polished (100) and (001) surfaces were mounted in a cryostat and measured in specular geometry with and , i.e. parallel and perpendicular to the c axis.

Analysis

The soft XAS M-edge data were analyzed with the full multiplet code XTLS 8.3 as described in ref. 41. The atomic parameters are given by Hartree-Fock values and the reduction factors for the Coulomb and spin orbit interaction are determined by fitting the isotropic spectrum Iiso = 1/3 in Fig. 6. We obtain reduction factors of 41% for the 4f-4f and 20% for the 3d-4f Coulomb interaction while the spin orbit values remain unchanged within 1 or 2%. These are typical values for Ce intermetallic compounds414243. For the reproduction of the polarized data, the CEF has to be taken into account. The calculations are purely ionic so that the spectral weights due to f 0 in the ground states and transitions to continuum states are not captured in the simulation.

Figure 6

Dots: isotropic spectrum as constructed from the experimental data .

Solid line: fit with full muliplet routine after adjusting reduction factors (see text).

The NIXS data were simulated using the full multiplet code Quanty58. A Gaussian and a Lorentzian broadening of FWHM = 1.3 eV and 0.2 eV, respectively, are assumed to account for the instrumental resolution and life time effects. The atomic Hartree-Fock values were adjusted via the peak positions, resulting in reductions of 30% and 20% for the 4f-4f and 4d-4f Coulomb interactions, respectively. The simulation procedure is explained in detail in ref. 44.

Band structure calculations

Density functional theory based calculations were performed using the full-potential non-orthogonal local orbital code (FPLO)59 employing the local density approximation (LDA) and including spin-orbit coupling. The band structure and density of states (DOS) are obtained from a carefully converged calculation with 1000 k-points in the irreducible wedge of the Brillouin zone. For easy comparison with literature, we have plotted the band structure along the k-path as described in ref. 54. To simulate the enhancement in the occupation of the various J states, we have applied a constant potential shift to the on-site energies (1 eV) of the relevant orbitals.

Additional Information

How to cite this article: Sundermann, M. et al. CeRu4Sn6: a strongly correlated material with nontrivial topology.. Sci. Rep. 5, 17937; doi: 10.1038/srep17937 (2015).