\section*{Introduction}

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 symmetry \cite{Kane_2005,Fu_2007,Zhang_2009}. 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 SmB$_6$\cite{Menth_1969} came immediately into consideration\cite{Dzero_2010, Takimoto_2011, Dzero_2012, Dzero_2013, Alexandrov_2013, Lu_2013} and first investigations\cite{Jiang_2013,Frantzeskakis_2013,Kim_2014,Li_2014,Xu_2014,Ishida_2015} show that the robust metallicity of the sample surface could be a promising concept for explaining its long-time mysterious low temperature residual conductivity\cite{Riseborough2000}. Increasing the interest for correlated topological insulators is the possibility to find phenomena that might be even more exotic\cite{PhysRevX.3.011016,1742-5468-2013-09-P09016,PhysRevB.88.115137,PhysRevX.3.041016} 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 behavior\cite{Dzero_2010, Takimoto_2011,Lu_2013}. One of the propositions is that rare earth Kondo or mixed-valent insulators should have all the necessary ingredients for realizing topological insulating properties\cite{Dzero_2010, Takimoto_2011, Dzero_2012, Dzero_2013, Alexandrov_2013, Lu_2013}: 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 cubic\cite{Dzero_2010,  Takimoto_2011} and tetragonal symmetry\cite{Dzero_2012, Dzero_2013} a gap always opens when the ground state is a nodeless $J_z$\,=\,$\pm 1/2$ Kramers doublet and when the $4f$ occupation is sufficiently away from integer.

Here CeRu$_4$Sn$_6$ comes into play. It is a tetragonal, non-centrosymmetric ($I\bar{4}2m$) \cite{Pottgen_1997}, potential Kondo insulator with a $J$\,=\,5/2 Hund's rule ground state. A hybridization-induced gap was reported first in 1992\,\cite{Das_1992}, supported further by several experimental studies\,\cite{Strydom_2005,Paschen_2010,Bruning_2010}, and later estimated to be of the order of 100\,K from single crystal resistivity data\,\cite{Paschen_2010,Winkler_2012}. Optical conductivity measurements point towards remnant metallicity\,\cite{Guritanu_2013} and yet the low temperature resistivity is in the $milli$Ohm-cm range, indicating clearly the depletion of the charge carriers\,\cite{Das_1992,Paschen_2010,Winkler_2012,Guritanu_2013}. 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)\,\cite{Paschen_2010,Winkler_2012}. $\mu$SR measurements did not find magnetic order down to 50\,mK, thus pointing out the importance of Kondo-type interactions in this material \,\cite{Strydom_2007} .

One of the pressing open questions is whether CeRu$_4$Sn$_6$ fulfills the conditions mentioned in the theoretical studies about strongly correlated topological insulators\,\cite{Dzero_2010, Takimoto_2011, Dzero_2012, Dzero_2013, Alexandrov_2013, Lu_2013}. To answer this question, we determine here whether the hybridization between conduction and $4f$ electrons ($cf$-hybridization) in CeRu$_4$Sn$_6$ 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_K$. 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 \textsl{Methods}.

\section*{Results}
\subsection*{Cerium $4f$ valence and Kondo temperature}
Trivalent cerium (Ce$^{3+}$) has a $4f$ occupancy $n_f$\,=\,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 $|\Psi_{\mathrm{GS}}\rangle$\,=\,$c_0$\,$|f^0\rangle$\,+\,$c_1$\,$|f^1$\uline{L}$\rangle$\,+\,$c_2$\,$|f^2$\uuline{L}$\rangle$, i.e$.$ with additional contributions of tetravalent ($f^0$) and divalent ($f^2$) states\cite{GunnarssonPRB28}. Here \uline{L} and \uuline{L} denote the number of ligand holes. For Ce the resulting occupation of the $4f$ shell is then given by $n_f$\,=\,2$c_2^2$\,+\,$c_1^2$\,$\leq$\,1, whereas $c_0^2$ quantifies the degree of delocalization. The hybridization is described in terms of the formation of a Kondo singlet state at $k_B T_K$ 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_K$ is obtained from the temperature variation of $n_f$, or equivalently of $c_0^2$ if $c_2^2$\,$\ll$\,$c_1^2$ or if $c_2^2$ does not change with $T$\cite{Bickers_1987,Tjeng_1993,Dallera_2002}.

Core level spectroscopy is a powerful tool for determining the occupations of the various $4f^n$ configurations\cite{GunnarssonPRB28}. 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 $L_{\alpha_1}$ recombination process ($3d_{5/2} \rightarrow 2p_{3/2}$) with partial fluorescence yield (PFY-XAS) at the Ce L$_3$-edge ($2p \rightarrow 5d$) have tremendously improved the accuracy of determining valencies\cite{Kotani_2001,Dallera_2002,Rueff_2010}. We have performed such a PFY-XAS experiment on CeRu$_4$Sn$_6$, determined $c_0^2$ and obtained $T_K$ from its temperature dependence.

\begin{SCfigure}
	\centering
	\includegraphics[width=0.5\linewidth]{Fig01.eps}
	\caption{(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 $L_3$-edge spectra and the $I$($f^0$) peak region on an enlarged scale. (c) Inset with $I$($f^0$)\,$\simeq$\,$c_0^2$ as function of log($T$). The relative error bars amount to $\pm$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_K$.}
	\label{Fig01}
\end{SCfigure}

Figure~\ref{Fig01}(a) shows the PFY-XAS data of CeRu$_4$Sn$_6$ 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 $L_3$-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^n$ weights, it is sufficient to fit the data with a set of Gaussian or similar functions as a mean to obtain their integrated intensities\cite{RueffPRL106,YamaokaPRL107,Kotani_2014}. 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 red line in Fig.~\ref{Fig01}(a) is the result of such a fit.

For the spectral weights we obtain $I$$(f^0)\!=\!0.059 \pm 0.015$, $I$$(f^1)\!=\!0.931 \pm 0.015$ and $I$$(f^2)\!=\!0.011 \pm 0.015$ which corresponds to $n_f \simeq 0.953$ and $c_0^2 \simeq 0.059 \pm 0.015$ under the assumption I$(f^n) \simeq c_n^2$. This is an approximation since it neglects hybridization effects in the final states\cite{Kotani_2014}. As in core level  photoemission\cite{Strigari_2015}, $I$($f^2$) overestimates $c_2^2$ and $I$($f^0$) underestimates $c_0^2$ so that the above values for $n_f$ and $c_0^2$ should be taken as upper and lower limits, respectively. It is therefore safe to conclude that the $cf$-hybridization is very strong in CeRu$_4$Sn$_6$.

Figure~\ref{Fig01}(b) shows the temperature dependence of the PFY-XAS data, normalized to the integrated intensity. The inset in Fig.~\ref{Fig01}(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 $f^1$ absorption line, yields the temperature variation $c_0^2(T)$ [upper inset, Fig.~\ref{Fig01}(c)]. The dashed line is a broadened step function describing the temperature evolution. We determine $T_K$ empirically as the inflection point of $c_0^2(T)$\cite{Bickers_1987,Tjeng_1993,Dallera_2002}, yielding $T_K$\,=\,$(170 \pm 25)$\,K supporting the analysis of the temperature dependence of the single crystal resistivity in Ref.\,\citenum{Guritanu_2013}.

\begin{SCfigure}
\centering
\includegraphics[width=0.6\linewidth]{Fig02.eps}
	\caption{Temperature dependence of the linearly polarized $M$-edge XAS spectra of CeRu$_4$Sn$_6$. The arrows mark the spectral weight due to the $f^0$ contribution in the ground state.}
	\label{Fig02}
\end{SCfigure}

\begin{SCfigure}
	\centering
	\includegraphics[width=0.5\linewidth]{Fig03.eps}
	\caption{(a) Linearly polarized XAS data (dots) at the Ce $M$-edges for $T$\,=\,4\,K and simulation with a $\Gamma_6$\,=\,$|5/2,\pm1/2\rangle$ ground state (lines). Blue corresponds to the polarization $E$\,$\parallel$\,[001] and red to $E$\,$\perp$[001]. (b) Experimental (dots) and simulated (colored lines) linear dichroism (LD), I$_{E\perp[001]}-\mathrm{I}_{E\parallel[001]}$, of the pure $J_z$ 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 $\Gamma_6$ symmetry.}
	\label{Fig03}
\end{SCfigure}

\subsection*{CEF ground state symmetry}
The local point symmetry at the Ce site of $D_{2d}$\cite{Pottgen_1997} in CeRu$_4$Sn$_6$ 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 $\Gamma_6\!=\!|5/2,\pm 1/2\rangle$, and two mixed states of type $\Gamma_7^1\!=\!\sqrt{1-\alpha^2}$\,$|5/2,\pm3/2\rangle+\alpha$\,$|5/2,\mp5/2\rangle$ and $\Gamma_7^2\!=\!\sqrt{1-\alpha^2}$\,$|5/2,\pm5/2\rangle-\alpha$\,$|5/2,\mp3/2\rangle$, in the $J_z$-representation $|J,J_z\rangle$. We have investigated the CEF ground state wave function with two spectroscopic methods: soft XAS at the Ce $M$-edge ($3d \rightarrow 4f$)\cite{HansmannPRL100,WillersPRB80,Strigari_2012} and hard x-ray NIXS at the Ce N-edge ($4d \rightarrow 4f$)\cite{WillersPRL109}. 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 $|\vec{q}|$ the signal is dominated by higher multipole orders of the scattering cross section\cite{Gordon2008,GordonJElecSpec184}. In the NIXS experiment the $\vec{q}$-dependence ($\vec{q}$-direction with respect to the sample orientation) replaces the polarization dependence in XAS and, accordingly, $\vec{q}$-dependent multipole selection rules give access to the ground state symmetry\cite{WillersPRL109}.

Figure~\ref{Fig02} shows the soft XAS data at the Ce M$_{4,5}$ edge for the two polarizations $E$\,$\parallel$\,[001] and $E$\,$\perp$[001] 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 $\geq$\,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~\ref{Fig03}(a) shows the linearly polarized low-temperature data (red and blue dots). For the simulation of the spectra a full multiplet code is used (see \textsl{Methods, Analysis}). Figure~\ref{Fig03}(a) shows a simulation assuming a pure $|5/2,\pm1/2\rangle$\,=\,$\Gamma_6$ ground state (red and blue lines), which is the state with rotational symmetry [see inset in Fig.~\ref{Fig03}(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.

In Fig.~\ref{Fig03}(b) the experimental linear dichroism (LD) ($I$$_{E\perp[001]}$-$I$$_{E\|001}$, black dots) is compared with the simulated LD of the three pure $|5/2,J_z\rangle$ states (colored lines). The LD of the pure $|5/2,\pm1/2\rangle$\,=\,$\Gamma_6$ (orange) and also of the $|5/2,\pm3/2\rangle$ (blue), corresponding to a $\Gamma_7$-symmetry with $\alpha$\,=\,0, have the correct sign, but the LD of the latter is too small and any $|5/2,\mp5/2\rangle$-admixture (green) in a $\Gamma_7$-type ground state would make the discrepancy even larger. Even if we adopt a $\Gamma_7$-type ground state of predominantly $|5/2,\pm3/2\rangle$, the measured LD will be larger than the simulated one. We therefore conclude that although the ground state is not a pure $\Gamma_6$\,=\,$|5/2,\pm1/2\rangle$, 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.~\ref{Fig04} (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 $\vec{q}$\,$\parallel$\,[001] and blue for $\vec{q}$\,$\parallel$\,[100]. Again there is a clear dichroic effect. We compare the measured $N_{4,5}$ edge with simulations for the pure $|5/2,J_z\rangle$ states (red and blue lines), which here are shown again (cf$.$ Ref.\,\citenum{WillersPRL109}) for reasons of convenience. Details about the theoretical simulations containing a full multiplet ansatz can be found in Section \textsl{Methods, Analysis}.

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

\begin{SCfigure}
 \centering
 \includegraphics[width=0.5\linewidth]{Fig04.eps}
	\caption{Dots: Experimental NIXS intensity at the Ce $N_{4,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 S($\vec{q},\omega$) for the pure $J_z$ states of the $J$\,=\,5/2 multiplet. Data and simulation are shown for $\vec{q}$\,$\parallel$\,[100] (blue) and $\vec{q}$\,$\parallel$\,[001] (red).}
	\label{Fig04}
\end{SCfigure}

\section*{Discussion}
What could be the reason for the deviation between the experiment and the simulation starting from a pure $\Gamma_6$\,=\;$|5/2,\pm1/2\rangle$? Already in 1986 van der Laan \textsl{et al.} pointed out that in the presence of strong hybridization higher lying states contribute to the $M$-edge signal\,\cite{Laan_1986}. In CeRu$_4$Sn$_6$ we find that the Kondo temperature (170 K) is not much smaller than the CEF excitations ($\geq$\,25\,meV\,$\hat{=}$\,290\,K) 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.~\ref{Fig03}(b) and (c), it is evident that the LD of this composite ground state will be reduced from that of the pure $|5/2,\pm1/2\rangle$  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 CeRu$_4$Sn$_6$ is predominantly of $|5/2,\pm1/2\rangle$ symmetry. This is in good agreement with the anisotropy of the low-temperature magnetization\,\cite{Winkler_2012} and susceptibility \cite{Sidorenko-unp}, thereby supporting  also the results of LDA+DMFT calculations at 290\,K\,\cite{Wissgott_Cond_Mat}.

We now will use density functional theory based band structure calculations to obtain a deeper understanding of the experimental findings (details in Section \textsl{Materials and Methods, Band structure calculations}). In Fig.~\ref{Fig05}\,(a) we display our result of a full-relativistic non-magnetic calculation employing the local density approximation (LDA) for CeRu$_4$Sn$_6$. We can clearly observe that the band structure has a true band gap in agreement with an earlier study in \,\cite{Guritanu_2013}, 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.~\ref{Fig05}\,(a) that the lowest unoccupied Ce $4f$ state is surprisingly given by the $J_z$\,=\,$|5/2,\pm3/2\rangle$ (blue dots) and not by the $|5/2,\pm1/2\rangle$ (red dots). Yet, looking at the Ce $4f$ partial density of states depicted in the right panel of Fig.~\ref{Fig05}\,(a) we find that it is the $|5/2,\pm1/2\rangle$ which has the largest occupation, as found in our experiment. We therefore can infer from these calculations that the empty $|5/2,\pm1/2\rangle$ 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, acquire a larger occupation.

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_z$ state so that it becomes more occupied, i.e. as to mimic the Ce $4f^1$ valence state thereby making the inversion with the Ru $4d$. Figure~\ref{Fig05}\,(b) shows the result for moving down the $|5/2,\pm1/2\rangle$ orbital and Fig.~\ref{Fig05}\,(c) the $J_z$\,=\,$|5/2,\pm3/2\rangle$. We can observe that in the $J_z$\,=\,$|5/2,\pm3/2\rangle$ case the system becomes metallic with quite a large number of bands crossing the Fermi level\,\cite{footnote1}. For the the $|5/2,\pm1/2\rangle$ 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 Galitzki \cite{Dzero_2012,Dzero_2013}, 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 CeRu$_4$Sn$_6$. 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_z$\,=\,$\pm 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).

\begin{figure}
	\centering
	\includegraphics[width=0.52\linewidth]{Fig05a.eps}
	\centering
	\includegraphics[width=0.495\linewidth]{Fig05b.eps}
	\centering
	\includegraphics[width=0.495\linewidth]{Fig05c.eps}
	\caption{(a) Full relativistic (LDA+spin-orbit) band structure of CeRu$_4$Sn$_6$, wherein the Ce 4$f$ weight is mostly above the Fermi level. From the $J_{z}$ projected DOS, we note the largest contribution to the occupied valence sector is from $J$ = 5/2 and $J_{z}$ = $\pm$ 1/2, (b) Applying a constant potential to the $J_z$\,=\,1/2 state. (c) Applying a constant potential to the $J_z$\,=\,3/2 state.}
	\label{Fig05}
\end{figure}

To summarize, we thus have found that CeRu$_4$Sn$_6$ fulfills the conditions mentioned by Dzero and Galitzki\cite{Dzero_2012,Dzero_2013} to be a strong topological insulator. It has a ground state symmetry given primarily by the $\Gamma_6$\,=\,$|5/2,\pm1/2\rangle$ 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 CeRu$_4$Sn$_6$ is a strongly correlated material with non-trivial topology.

\section*{Methods}
\subsection*{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 furnace\cite{Prokofiev_2012}. 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\,\AA\ and $c$\,=\,9.7520\,\AA\ as given in Ref.\,\citenum{Pottgen_1997}. The single crystals were pre-aligned with a Laue camera. A peculiarity of the present tetragonal structure is that $\sqrt{2}$$a$\,$\approx$\,$c$ within $\approx$\,0.2\%. As a consequence the structure appears pseudo-cubic in Laue images and thus the final alignment of the single crystals had to be done via magnetization measurements\cite{Paschen_2010}.

\subsection*{Set-ups of x-ray techniques} The PFY-XAS experiment was performed at the GALAXIES inelastic scattering endstation at the SOLEIL synchrotron in France\cite{GalaxiesPaper}. The syn\-chro\-tron radiation was monochromatized using a Si(111) nitrogen-cooled fixed-exit double-crystal monochromator ($\Delta E/E\!\approx\!1.4 \times 10^{-4}$), followed by a Pd-coated spherical collimating mirror. The x-rays were then focused to a spot size of 30\,$\mu$m\,x\,90\,$\mu$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$\times$10$^{13}$ photons/second on the sample. The incident x-ray energy was scanned across the Ce L$_3$ absorption edge ($\approx$\,5723\,eV) while monitoring the intensity of the Ce L$_{\alpha_1}$ de-excitation ($\approx$\,4839\,eV) using a Ge\,331 spherical bent crystal analyzer (Bragg angle 80.8$^{\circ}$, 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 CeRu$_4$Sn$_6$ single crystals were cleaved \textsl{in situ} to guarantee clean sample surfaces. The TEY signal was normalized to the incoming photon flux $I_0$ as measured at the refocusing mirror. The energy resolution at the cerium M$_{4,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 $E$\,$\parallel$\,[001] and $E$\,$\perp$\,[001], [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 \textsl{in-plane} scattering angles of $2\,\Theta = 140^\circ, 146^\circ$, and $153^\circ$, corresponding to an averaged momentum transfer of $|\vec{q}|\!\approx\!(9.4 \pm 0.2)$\,\AA$^{-1}$. This configuration yields an instrumental energy resolution of FWHM\,=\,1.3\,eV and made the Ce N$_{4,5}$ edge measurable beyond the dipole limit\cite{Gordon2008,GordonJElecSpec184,WillersPRL109}. Two crystals with polished (100) and (001) surfaces were mounted in a cryostat and measured in specular geometry with $\vec{q}$\,$\|$\,[001] and $\vec{q}$\,$\bot$\,[100], i.e$.$ parallel and perpendicular to the $c$ axis.

\subsection*{Analysis} The soft XAS $M$-edge data were analyzed with the full multiplet code \textsc{XTLS}~8.3 as described in Ref.\,\citenum{HansmannPRL100}. The atomic parameters are given by Hartree-Fock values and the reduction factors for the Coulomb and spin orbit interaction interaction are determined by fitting the isotropic spectrum I$_\mathrm{iso}$\,=\,1/3\,(I$_{E\parallel[001]}$\,+\,2I$_{E\perp[001]}$) in Fig.~\ref{Fig06}. 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 compounds\cite{HansmannPRL100,WillersPRB80,Strigari_2012}. For the reproduction of the polarized data, the CEF has to be taken into account. The calculations are purely ionic.

\begin{SCfigure}
\centering
\includegraphics[width=0.5\linewidth]{Fig06.eps}
	\caption{Dots: isotropic spectrum as constructed from the experimental data $I_{iso}$\,=\,1/3\,($I_{E\parallel[001]}$\,+\,2$I_{E\perp[001]}$). Solid line: fit with full muliplet routine after adjusting reduction factors (see text).}
	\label{Fig06}
\end{SCfigure}

The NIXS data were simulated using the full multiplet code \textsl{Quanty}\cite{MauritsQuanty}. 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 4$f$-4$f$ and 4$d$-4$f$ Coulomb interactions, respectively. The simulation procedure is explained in detail in Ref.\,\citenum{WillersPRL109}.

\subsection*{Band structure calculations}
Density functional theory based calculations were performed using the full-potential non-orthogonal local orbital (\texttt{FPLO})\cite{Koepernik_1999} code 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.\,\citenum{Wissgott_Cond_Mat}.
To simulate the enhancement in the occupation of the various $J_{z}$ states, we have applied a constant potential shift to the on-site energies (1 eV) of the relevant orbitals.

\section*{Acknowledgment}
We thank K.~Held, P.~Thalmeier and D.T.~Adroja for fruitful discussions. M.S., F.S, T.W. and A.S. benefited from support of the German funding agency DFG (Projects 600575 and 615811) and H.W., A.P. and S.P. from the Austrian Science Fund (Project I623-N16) and the U.S. Army Research Office (Grant W911NF-14-1-0497). We further acknowledge SOLEIL for provision of synchrotron radiation facilities (proposal 20120919).

\section*{Author contributions statement}
A.S. and L.H.T designed the research,
M.S., F.S., T.W., J.M.A., J.-P.R., D.S., E.W, M.M.S., A.A-Z., L.H.T, and A.S. performed the experiments,
H.W., A.P., and S.P. made high quality singly crystals,
A.T. and M.W.H. contributed codes for data analysis,
M.S., F.S., D.K., and A.S. analyzed data,
M.S., L.H.T., S.P. and A.S. wrote the paper.
All authors reviewed the manuscript.


\section{Introduction}

Sixty eight percent of our universe\cite{Planck2013cosmological} consists of a still mysterious component called ``dark energy'' (DE), which is believed to be responsible for the  present acceleration of the universe\cite{reiss1998, perlmutter1999}. Among a wide range of alternatives for the dark energy  (see Ref.~\refcite{copeland2006dynamics} for review),  which includes the cosmological constant, scalar or vector fields, modifications of gravity and different kinds of cosmological fluids, the usage of a canonical scalar field, called ``quintessence'', is a viable and natural candidate\cite{peebles1988,ratra1988,Frieman1992,Frieman1995,Caldwell:1997ii}. Another striking attempt  to explain the acceleration comes from holography. The holographic principle states that the degrees of freedom of a physical system scales with its boundary area  rather than its volume\cite{'tHooft:1993gx,Susskind:1994vu}. Cohen and collaborators\cite{Cohen:1998zx} suggested that the dark energy should obey this principle, thus its energy density  has an upper limit and the fine-tuning problem for the cosmological constant is eliminated. Refs.~\refcite{Hsu:2004ri,Li:2004rb} followed the previous ideas regarding holography and argued that the holographic dark energy (HDE) has an energy density given by  $\rho_D=3c^2/L^2$ ($M_{p}^{-2}\equiv 8\pi G=1$), where $c$ is a constant and $L$ is an infrared (IR) cutoff. The first choice for $L$ was  the Hubble radius $H^{-1}$, however it led to an equation of state that describes dust\cite{Hsu:2004ri}. The correct equation of state for dark energy was obtained by \cite{Li:2004rb}, when he chose the future event horizon as the IR cutoff.
 The problem with the Hubble radius as the IR cutoff could be avoided assuming the interaction between dark energy and dark matter (DM)\cite{Pavon:2005yx}. Such interaction was first proposed in the context of quintessence\cite{Wetterich:1994bg,Amendola:1999er}, and since the energy  densities of the DE and DM are comparable, the interaction can alleviate the coincidence problem\cite{Zimdahl:2001ar,Chimento:2003iea}. Taking the coupling into account, the HDE with Hubble radius as IR cutoff could lead to an accelerated expansion of the universe and also solve the coincidence problem\cite{Pavon:2005yx}. Li's proposal could also be generalized assuming the interaction between the two components of the dark sector\cite{Wang:2005jx,Wang:2005pk,Wang:2005ph,Wang:2007ak}.

   From the point of view of theoretical physics, it would be interesting find out a model of dark energy from first principles. Since supergravity is  the low-energy limit of the superstring theory, it is natural to investigate if it can provide a model that describes the accelerated expansion of the universe. The simplest supergravity case is with one supersymmetry ($\mathcal{N}=1$)\footnote{ Extended supergravities can be applied to cosmology as well, see for instance \cite{kallosh2002supergravity,kallosh2002gauged}. They are closer to string theory, however minimal supergravity can work as an effective theory or an approximation.} and in this framework  Refs.~\refcite{Brax1999,Copeland2000,Landim:2015upa} presented some models that try to describe dark energy through quintessence.  Due to the prominent role of the AdS/CFT correspondence\cite{Maldacena:1997re} to relate both supergravity and holography concepts, it is natural to ask if there is any connection between supergravity (even for $\mathcal{N}=1$) and HDE. As showed in Ref.~\refcite{Sheykhi:2011cn}, it is possible to establish a connection between HDE and different kind of scalar fields, such as quintessence, tachyon and K-essence, through reconstructed scalar potentials.

    In this paper we show that the HDE in interaction with DM can be embedded in minimal supergravity plus matter with a single chiral superfield. The interaction between HDE and DM is taken into account because they are generalizations of some uncoupled cases, namely, the models presented in Ref.~\refcite{Hsu:2004ri,Li:2004rb}. We use Planck units ($\hbar=c=8\pi G=1$) throughout the text.

  We assume that the dark energy is coupled with dark matter, in such a way that total energy-momentum is still conserved. In the flat Friedmann-Robertson-Walker background with a scale factor $a$, the continuity equations for both components are

\begin{equation}\label{contide}
\dot{\rho_D}+3H(1+w_D)\rho_D=-Q, \quad \dot{\rho_M}+3H\rho_M=Q,
\end{equation}

\noindent  where $H=\dot{a}/a$ is the Hubble rate,  $Q$ is the coupling and the dot is a derivative with respect to the cosmic time $t$. The index $M$ stands for dark matter. The case $Q>0$ corresponds to dark energy transformation into dark matter, while $Q<0$ is the transformation in the opposite direction. In principle, the coupling  can depend on several variables $Q=Q(\rho_m,\rho_\phi, H,\dots)$, but assuming it is a small quantity, the interaction term can be written as a Taylor expansion, giving rise to three main kernels: (i) $Q\propto H\rho_M$, (ii)   $Q\propto H\rho_D$ and (iii)  $Q\propto H(\rho_M+\rho_D)$. In the context of HDE, the second kernel was used in Ref.~\refcite{Pavon:2005yx}), while the third one is found in Ref.~\refcite{Wang:2005jx,Wang:2007ak}. Both kernels generalize  previous results\cite{Hsu:2004ri,Li:2004rb}, respectively, thus the correspondent models are good choices to analyze whether they can be embedded in minimal supergravity or not.

\section{HDE with Hubble radius as IR cutoff}

 First we will consider the size $L$  as the Hubble radius $H^{-1}$, thus the energy density for the dark energy becomes

  \begin{equation}\rho_D=3c^2H^2, \label{eqH}
   \end{equation}

  \noindent which describes a pressureless fluid\cite{Hsu:2004ri} in the absence of interaction with dark matter. When one considers such interaction\cite{Pavon:2005yx}, the fluid  has an equation of state that can describe the dark energy. We consider the interaction $Q=3b^2H\rho_D$, where $b$ is a constant which measures the strength of the interaction. The equation of state for dark energy with this interaction is written in terms of the ratio of the energy densities $r\equiv \rho_M/\rho_D=(1-c^2)/c^2$

   \begin{equation}w_D=-\left(1+\frac{1}{r}\right)b^2=-\frac{b^2}{1-c^2}.
   \label{eqPavon}
   \end{equation}

In order to reconstruct a holographic scalar field model we should relate (\ref{eqPavon}) with the scalar field $\phi$. Using the energy density and the pressure for the scalar field $\rho_\phi=\dot{\phi}^2/2+V(\phi)$ and $p_\phi=\dot{\phi}^2/2-V(\phi)$, we have

   \begin{equation}\dot{\phi}^2=(1+w_\phi)\rho_\phi, \quad V(\phi)=(1-w_\phi)\frac{\rho_\phi}{2}.
   \label{scalarV}
   \end{equation}

Using (\ref{eqH}) and (\ref{eqPavon}) into (\ref{scalarV}) we have the potential\cite{Sheykhi:2011cn}

 \begin{equation}
 V(\phi)=B^{-2}e^{-\sqrt{2}B\phi},
   \label{V1}
   \end{equation}

   \noindent where $B=\frac{3k}{2\sqrt{2}c}\left(3-\frac{3b^2}{1-c^2}\right)^{-1/2}$ and $k=1-b^2c^2/(1-c^2)$.

On the other hand, the scalar potential in minimal supergravity with no D-terms is given by

\begin{equation}V= e^{K}\left( K^{\Phi\bar{\Phi}}\left|W_\Phi+K_\Phi W\right|^2-3|W|^2\right),\label{potential}
   \end{equation}

    \noindent for a single complex scalar field. The potential above  depends on a real function $K\equiv K(\Phi,\bar{\Phi})$, called K\"ahler potential, and a holomorphic function $W\equiv W(\Phi)$, the superpotential. $K_{\Phi\bar{\Phi}}\equiv \frac{\partial^2 K}{\partial \Phi\partial\bar{\Phi}}$ is the K\"ahler metric, $W_\Phi\equiv\frac{\partial W}{\partial \Phi}$ and $K_\Phi\equiv\frac{\partial K}{\partial \Phi}$. We will use the same choices of Ref.~\refcite{Copeland2000} for $K$ and $W$, namely, the string-inspired K\"ahler potential $K=-\ln(\Phi+\bar{\Phi})$, which is present at the tree level for axion-dilaton field in string theory, and the superpotential $W=\Lambda^2\Phi^{-\alpha}$, where $\Lambda$ is a constant. With these choices, with the imaginary part of the scalar field stabilized at zero $\langle$Im $\Phi\rangle=0$ and the field redefinition Re $\Phi=e^{\sqrt{2}\phi}$, we get the scalar potential expressed in terms of the canonical normalized field $\phi$

 \begin{equation}V(\phi)= \frac{\Lambda^4}{2}(\beta^2-3)e^{-\sqrt{2}\beta\phi},\label{potentialsugra}
   \end{equation}

   \noindent where  $\beta=2\alpha+1$. Comparing (\ref{V1}) and (\ref{potentialsugra}) we have $B=\beta$ and $\Lambda^2=\sqrt{2}/(\beta\sqrt{\beta^2-3})$, with $\beta^2\geq 3$. Thus, the two parameters $b$ and $c$ determine the exponent $\alpha$ of the superpotential and the parameter $\Lambda$. We notice in this case that $\Lambda$ is not an independent parameter, but it depends on $b$ and $c$ as well. Even if we would have a way to know what is the value of $\alpha$ or rather, $\beta$, we could not know the specific values of $b$ and $c$. The opposite direction is favored, because once one finds out observationally $b$ and $c$, the HDE can be embedded in a specific supergravity model. If there is no interaction ($b=0$) $\beta=\sqrt{3}/(2c\sqrt{2})$, but $w_D=0$.

   The kernel used so far simplifies $w_D$ in such a way that (\ref{V1}) was written with no need of any approximation. When one considers the other possibilities for the interaction $Q$ [(i) or (iii)], the scalar potential cannot be written as easy as it was in (\ref{V1}). We illustrate this possibility now, but with other IR cutoff.

  \section{HDE with future event horizon as IR cutoff}
   We will analyze another possibility of HDE with the kernel $Q=3b^2H(\rho_M+\rho_D)$ and with the future event horizon $R_E=a\int_t^\infty{dt/a}=c\sqrt{1+r}/H$ as the IR cutoff.  The choice of $L$ is done because when $b=0$ the original Li's model of HDE\cite{Li:2004rb} is recovered. The energy density for the HDE is

    \begin{equation}\label{rho2}
   \rho_D=\frac{3H^2}{1+r}
   \end{equation}

\noindent and the equation of state for HDE becomes

     \begin{equation}
     w_D=-\frac{1}{3}\left(1+\frac{2\sqrt{\Omega_D}}{c}+\frac{3b^2}{\Omega_D}\right),
   \label{eqabdalla}
   \end{equation}

  \noindent where $\Omega_D\equiv\rho_D/(\rho_D+\rho_M)=(1+r)^{-1}$ is the density parameter for DE. Using (\ref{rho2}) and (\ref{eqabdalla}) in (\ref{scalarV}) we have

        \begin{equation}
    \dot{\phi}^2=2H^2\Omega_D\left(1-\frac{\sqrt{\Omega_D}}{c}-\frac{3b^2}{2\Omega_D}\right),
   \label{phidot}
   \end{equation}

      \begin{equation}
    V(\phi)=H^2\Omega_D\left(2+\frac{3\sqrt{\Omega_D}}{c}+\frac{9b^2}{2\Omega_D}\right).
   \label{V}
   \end{equation}

 \noindent From the equations above we see that we have to know the evolution of $\Omega_D$ in order to fully determine  $V(\phi)$ as a function of $\phi$. The evolution equation  for $\Omega_D$ was found in Ref.~\refcite{Wang:2005jx} and it is

      \begin{equation}
   \frac{\Omega_D'}{\Omega_D^2}=(1-\Omega_D)\left(\frac{1}{\Omega_D}+\frac{2}{c\sqrt{\Omega_D}}-\frac{3b^2}{\Omega_D(1-\Omega_D)}\right),
   \label{omegaevolut}
   \end{equation}

   \noindent where the prime is the derivative with respect to $\ln a$. To solve this differential equation we let $y=1/\sqrt{\Omega_D}$, thus the equation (\ref{omegaevolut}) becomes

    \begin{equation}
   y^2y'=(1-y^2)\left(\frac{1}{c}+\frac{y}{2}+\frac{3b^2y^3}{2(1-y^2)}\right).
   \label{yevolut}
   \end{equation}

 In order to have an analytic solution of (\ref{yevolut}) we will investigate the asymptotic behavior of $\Omega_D$, for small and large $a$. For very small $a$ we have $\Omega_D\rightarrow 0$ and $y\rightarrow \infty$, therefore (\ref{yevolut}) is approximately

    \begin{equation}
   y'\approx(1-y^2)\left(\frac{3b^2}{2}-\frac{1}{2}\right)y,
   \label{yevolut2}
   \end{equation}

  \noindent which leads to the solution $\Omega_D=\Omega_0 a^{-(3b^2-1)}\approx \Omega_0 a $, provided that $b$ should be small\cite{Wang:2005jx}. Since  $\Omega_D$ scales with $a$, the contribution of the dark energy in the early universe is negligible. We also see from (\ref{phidot}) that the term between brackets should be positive and small, then $ \dot{\phi}\propto H\sqrt{2\Omega_D} $ is small. Therefore, we have $\phi(a)\sim a^{1/2}$. Due to the small contribution of  $\Omega_D$, we neglect it at this limit and we will focus on large $a$, where HDE is dominant.

  For large $a$ we have $\Omega_D\rightarrow 1$, and $y\rightarrow 1$. Thus, the last term in the equation (\ref{yevolut}) is the dominant one, so $y'\approx\frac{3b^2}{2}y$, and $\Omega_D=\Omega_0a^{-3b^2}$, where $\Omega_0$ is the value of $\Omega_D$ at $a_0=1$. Since $\Omega_D=\Omega_0a^{-3b^2}=\Omega_0a^{-3(1+w_D)}$, we have $1+w_D=b^2$ which is very small if we consider the present value of $w_D$ given by Planck\cite{Planck2013cosmological}. Thus we neglect the last term in (\ref{phidot}) and we assume the dark energy dominance. We have  $\Omega_D=\Omega_0\approx 1$ at large $a$ and (\ref{phidot}) yields

   \begin{equation}
    \dot{\phi}^2\approx 2H^2\Omega_D\left(1-\frac{1}{c}\right).
   \label{phidot2}
   \end{equation}

  \noindent The equation above implies that $w_D\approx -\frac{1}{3}(1+2/c)$, which is the Li's proposal\cite{Li:2004rb}. Equation (\ref{phidot2}) has the solution

     \begin{equation}
   \phi(a)=\left(1-\frac{1}{c}\right)^{1/2}\sqrt{2\Omega_D}\ln a.
   \label{phidota}
   \end{equation}

   The second Friedmann equation at large $a$ leads to  $H=\frac{3c}{(c-1)t}$ and $a=t^{3/(c-1)}$. The scalar potential (\ref{V}) becomes

    \begin{equation}
V(\phi)\approx\left(2+\frac{3}{c}\right)H^2=\left(2+\frac{1}{c}\right)\frac{9c^2}{(c-1)^2}e^{-\frac{1}{3}\left(\frac{c-1}{c}\right)^{1/2}\phi}.\label{V2}
   \end{equation}

   Similarly to before, comparing (\ref{V2}) with (\ref{potentialsugra}) we see that both $\beta$ and $\Lambda$ are determined by the constant $c$. Equation (\ref{V2}) was deduced for a dark-energy-dominated universe, in such a way that the interaction with dark matter is absent, as it should be in this limit.

   A more realistic scenario could not be found analytically, in opposite to the previous case, where the scalar potential was written with no approximations.  Different kernels for the coupling $Q$ may be tried, although the main features of the method were expressed in these two cases. The other alternatives are similar.

  \section{Conclusions}
 In this paper we embedded two models of HDE in the minimal supergravity with one single chiral superfield. In the first model we used the Hubble radius as the IR cutoff and an interaction proportional to $\rho_D$, while in the second one we used the future event horizon as IR cutoff and the kernel proportional to $\rho_D+\rho_M$. In both cases the free parameters  of our superpotential could be expressed in terms of the constants $c$ and $b$, depending on the model. The second case is embedded only for a dark-energy-dominated universe, while the first one is more general. There are other ways to embed HDE in minimal supergravity, as for instance the choices of K\"ahler potential and superpotential made in Ref.~\refcite{Nastase:2015pua} for the inflationary scenario. However, the main results are the same, that is, a way to relate HDE and supergravity. Due to the nature of the holographic principle and recalling that extended supergravity is the low-energy limit of string theory, the relation presented here has perhaps a deeper meaning, when one takes a quantum gravity theory into account.

\section*{Acknowledgments}
I thank Elcio Abdalla for comments. This work is supported by FAPESP Grant No. 2013/10242-1.

\bibliographystyle{ws-ijmpd}