1 Introduction
Narrow-gap semiconductors have emerged as a prominent research frontier in condensed matter physics and materials science due to their unique electronic structures and tunable physical properties [
1−
3]. These materials typically exhibit bandgaps on the order of several hundred millielectronvolts (meV), corresponding to excitation energies within the infrared (IR) regime, making them widely applicable in fields such as infrared photodetection, remote sensing, and low-power IR devices [
4−
10]. For instance, indium antimonide (InSb, bandgap ~ 0.17 eV) exhibits high sensitivity to infrared radiation in the 1−5 μm wavelength range, making it a key material for cryogenic short-wave infrared (SWIR) detectors [
11,
12]. On the other hand, the large carrier effective masses in narrow-gap semiconductors can yield Seebeck coefficients that greatly exceed those of metals and wide-gap semiconductors. When combined with tunable carrier densities, this facilitates the achievement of a high thermoelectric figure of merit,
ZT =
S2σT/
κ, where
S,
σ,
T, and
κ represent the Seebeck coefficient, electrical conductivity, absolute temperature, and total thermal conductivity, respectively [
1]. A paradigmatic case is Bi
2Te
3 (bulk bandgap ~ 0.15 eV) [
13], whose large Seebeck coefficient and intrinsically low lattice thermal conductivity make it one of the most outstanding thermoelectric materials near room temperature and have enabled its deployment in commercial solid-state cooling devices [
14−
17].
Although narrow-gap semiconductors have achieved significant applications in infrared technology and thermoelectric, current research actively searches for new narrow-gap compounds that exhibit unconventional electronic or thermoelectric behavior, with the goal of extending operating temperature ranges, improving efficiencies, introducing new functionalities, and reducing reliance on rare-earth elements [
18]. While most established narrow-gap materials can be optimized through band engineering and doping, their physical properties largely define by the bandgap within the single-particle picture without any electronic correlations. To move beyond this conventional framework of semiconductors, research attention is shifting toward non-trivial topology, reduced dimensionality, and strong electronic correlations. These systems may leverage many-body effects to achieve superior or novel functionalities, catering to more extreme application demands. Among these emerging candidate materials, iron diantimony (FeSb
2) has attracted considerable interest due to the series of anomalous physical properties observed at low temperatures [
1].
Initially, FeSb
2 was regarded as a conventional narrow-gap semiconductor. Subsequent studies revealed that it exhibits a record-breaking Seebeck coefficient, reaching approximately 45 mV/K at around 10 K, which far exceeds conventional theoretical predictions [
19]. Such an exceptional thermoelectric response suggests immense potential for thermoelectric applications. If the lattice thermal conductivity can be effectively suppressed, the thermoelectric figure of merit could be significantly enhanced, enabling the material to play a pivotal role in solid-state cryogenic cooling and thermoelectric power generation in the extreme low-temperature regime.
Since the discovery of its exceptionally large low-temperature thermoelectric response, FeSb
2 has drawn sustained attention from both the condensed-matter and thermoelectric communities. Early studies focused mainly on the microscopic origin of this remarkable macroscopic behavior, leading to extensive discussion of mechanisms such as carrier excitation, correlation-driven renormalization, and phonon drag. At the same time, considerable effort has been directed toward improving its thermoelectric performance through elemental substitution, defect control, and nanostructure engineering aimed at reducing lattice thermal conductivity and enhancing the thermoelectric figure of merit [
20−
28]. More recently, the scope of FeSb
2 research has expanded significantly. Between 2020 and 2025, conceptual and experimental approaches from topological quantum materials, strongly correlated electron physics, and X-ray spectroscopy have increasingly been applied to this system. Meanwhile, the emergence of ideas such as altermagnetism has renewed interest in spin-dependent phenomena in transition-metal compounds and has opened additional perspectives for understanding the electronic structure and magnetism of FeSb
2. As a result, current research extends well beyond phonon drag and thermal transport to encompass issues such as surface-bulk separation and correspondence, lattice dynamics, and their consequences for electronic structure and charge transport [
29−
36].
Beyond its thermoelectric behavior and these newly emerging directions, FeSb
2 also exhibits a range of unusual fundamental properties reminiscent of strongly correlated electron systems. Although FeSb
2 exhibits anisotropic magnetic susceptibility at low temperatures, it shows diamagnetic (
H||(001)) or nonmagnetic (
H is along other direction) at the low temperature [
37]. With the temperature ramping up, the magnetic response get profoundly enhanced to paramagnetic at temperature around 100 K. Correspondingly, its electrical resistivity exhibits a temperature-controlled metal−insulator-like transition. These behaviors are highly analogous to those of classic correlated semiconductors with Kondo interaction, such as SmB
6 and FeSi, suggesting the existence of a novel electronic state in FeSb
2 dominated by electron correlations [
38−
40]. Therefore, elucidating the microscopic origin of FeSb
2’s anomalous properties is crucial for both advancing our fundamental understanding of correlation-driven phenomena and guiding the rational design of materials with tailored functionalities.
This review provides a systematic overview of the intriguing physical phenomena in FeSb2, with the aim of developing a comprehensive understanding of this material and offering insights into other correlated semiconductors with similar properties. The article is organized as follows. Section 2 introduces the puzzling and exotic properties of FeSb2 and summarizes the proposed physical explanations. Section 3 focuses on the physics of phase transitions, examining interpretations based on the Kondo insulator (KI) model and comparing FeSb2 with representative Kondo materials such as SmB6. Section 4 reviews alternative theoretical models, mainly the spin-state excitation (SSE) hypothesis, and provides a systematic analysis of ongoing controversies. Section 5 concentrates on microscopic evidence from recent spectroscopic probes, including X-ray absorption and photoemission spectroscopy, and discusses the constraints on theoretical models. Finally, Section 6 outlines future research directions for FeSb2, highlighting potential applications in thermoelectricity, spintronics, and novel magnetic phenomena like altermagnetism. By integrating these interconnected aspects, this review aims to elucidate the rich and nontrivial physics hidden within this seemingly simple material and to offer perspectives on future research directions.
2 Unusual physical properties of FeSb2
FeSb2 exhibits unusually strong sensitivity of its transport properties to subtle variations in microscopic disorder and sample conditions, including stoichiometric deviations, residual impurities or defects, strain, and growth conditions. As a result, transport and magnetic measurements reported in the literature can differ qualitatively even among nominally similar samples. Unless otherwise stated, the following discussion of macroscopic transport behavior and microscopic modeling takes high-quality, near-intrinsic single crystals as the reference case.
2.1 Anomalous electrical transport
FeSb
2 crystallizes in the orthorhombic marcasite structure with the space group
Pnnm [Fig. 1(a)] [
41]. While early studies on polycrystalline FeSb
2 suggested it as narrow-gap semiconductor, the successful growth of high-quality single crystals in the early 2000s revealed a far more complex physical behavior. A notable feature of FeSb
2 is its profound anisotropy in electrical transport. Petrovic
et al. [
42] demonstrated that while the resistivity along the
a and
c axis exhibits typical semiconducting behavior (d
ρ/d
T < 0) over the entire temperature range, the
R−
T along
b-axis behaves in a qualitatively different manner. Along this axis, the resistivity is metallic (d
ρ/d
T > 0) at high temperatures, passes through a minimum around
Tmin ~ 40 K, and subsequently evolves into insulating-like behavior upon further cooling. However, not all reports reproduce a stark divergence among the three crystallographic axes. Some measurements [
19,
43] have instead observed nearly isotropic semiconducting behavior along all axes down to low temperatures, arguing that the large anisotropy and the metallic segment along
b-axis may depend sensitively on subtle structural distortions or sample-specific differences (e.g., strain, defects) [Fig. 1(c)]. Later on, Raman spectroscopy, identified the phonon anomaly along the
b-axis at
T ~ 100 K, indicating a possible lattice distortion along the
b-axis [
44]. Thus, one must treat the metal to insulator transition along the
b-axis with caution, acknowledging possible the sample-dependence due to possible off-stoichiometry effect.
Detailed analysis of the temperature-dependent resistivity revealed the presence of multiple energy scales. Through the fitting analysis of Arrhenius equation, several research groups independently identified two distinct thermal activation energy gaps. A larger gap (Δ
1) of ~ 30−40 meV is observed at higher temperatures (
T > 30 K) and is typically associated with the intrinsic semiconductor band gap. Below this temperature, however, a much smaller gap (Δ
2) of ~ 6−10 meV dominant the transport behaviors [
42,
45,
46]. The origin of this smaller gap remains an open question. Although it was initially suspected to be an impurity band, its reproducible observation in high-purity crystals suggests that it is possibly an intrinsic feature and may be associated with Kondo-like hybridization, a topic explored later in this review.
The most enigmatic feature is a distinct resistivity saturation plateau observed below 10 K, a feature inconsistent with simple semiconducting models. For years, the origin of this plateau was debated, as conventional four-probe measurements could not separate bulk from potential surface contributions. A breakthrough came from Corbino-disk geometry experiments [
47], which isolate bulk conduction from surface paths, demonstrated that the bulk resistivity of FeSb
2 does not saturate but continues to increase exponentially upon cooling. This compelling result proves that the bulk of FeSb
2 is truly insulating and that the observed resistivity plateau is definitively caused by highly conductive metallic surface channels. This finding has fundamentally shifted the focus of the field. The central question is no longer whether non-bulk conduction channels exist or not, but rather what the exactly nature of these conductive surface states is and whether they possess non-trivial topological character [
38,
48].
2.2 Temperature-dependent magnetic susceptibility
Complementing the complex electric transport phenomena, the magnetic properties of FeSb
2 provided the evidence for the possible role of strong electron correlations. Systematic measurements of the magnetic susceptibility [χ(
T)] on high-quality single crystals revealed a highly unconventional temperature dependence [Fig. 2(a)] [
42,
49]. At low temperatures, the system is diamagnetic with a nearly flat, negative susceptibility. As the temperature increases, χ(
T) crosses over to a paramagnetic state, passing through a broad maximum centered around 100 K before following a Curie-Weiss-like behavior at higher temperatures. This characteristic behavior has since been confirmed in numerous studies, including those under high pressure which showed the feature to be robust but tunable [Fig. 2(b)] [
50].
The paramagnetic anomaly observed in FeSb
2 is unexpected given its low-spin Fe configuration (3
d4 or 3
d6), as confirmed by X-ray absorption spectroscopy (XAS) [
51]. However, this non-monotonic behavior is strikingly reminiscent of the magnetic signature of the canonical KI candidate, FeSi [
52]. In FeSi, both the thermally activated resistivity and the crossover from a low-temperature diamagnetic state to a high-temperature paramagnetic state are well described within a hybridization-gap framework supplemented by dynamical mean-field theory (DMFT), which captures a temperature-driven evolution from a nonmagnetic insulating ground state to a correlated paramagnetic semiconductor [Fig. 2(c)] [
53]. The striking similarity between FeSi and FeSb
2 therefore indicates a possibly shared microscopic mechanism rooted in strong hybridization and electron–electron correlations that drive gap formation [
54]. From this perspective, the magnetic anomaly redefined the classification of FeSb
2 from a simple semiconductor to a candidate correlated electron system. This analogy to FeSi has strongly influenced the theoretical understanding ever since, serving as a primary piece of evidence for models based on a Kondo-like hybridization gap−a picture that remains central to the ongoing debate over the material’s microscopic origins.
2.3 The colossal Seebeck effect and the thermoelectric dilemma
In addition to its intriguing transport and magnetic behaviors, the observation of a colossal thermoelectric effect has made FeSb
2 a material of significant scientific interest. Bentien
et al. [
19] reported in a series of studies that high-quality FeSb
2 single crystals exhibit a record-breaking Seebeck coefficient at low temperatures, reaching a peak value of approximately −45 000 μV/K around 10−12 K [Fig. 3(a)]. Combined with the relatively low resistivity in the corresponding temperature range, this resulted in a peak thermoelectric power factor (PF =
S2σ) as high as 2300 μW/(cm·K
2), a value several tens of times greater than that of state-of-the-art commercial thermoelectric materials like Bi
2Te
3, indicating significant potential for low-temperature refrigeration applications. Based on this, Sun
et al. [
54,
55] showed that FeSb
2’s power factor far exceeds that of the isostructural compounds FeAs
2 and RuSb
2, arguing that the enhancement of
ZT is intrinsic to the narrow, correlated bands of FeSb
2 rather than a universal property of its crystal structure [Fig. 3(a)]. Complementarily, Takahashi
et al. [
56] demonstrated a strong dependence of the Seebeck coefficient on sample size and geometry, with smaller crystals showing enhanced peak values and narrower peak widths around 10 K. This behavior indicates that the thermoelectric response in FeSb
2 is markedly influenced by the phonon mean free path and boundary conditions [Fig. 3(b)]. Together, these compelling experimental findings spurred the development of theoretical frameworks in understanding the colossal thermoelectricity in FeSb
2. In particular, Tomczak
et al. [
57] involved many-body theory to demonstrate that conventional band pictures cannot reproduce the extreme low-
T Seebeck coefficient of FeSb
2. They further proposed that vertex corrections (e.g., electron−phonon drag) beyond simple quasiparticle renormalization in the calculation are essential, placing FeSb
2 in a distinct regime among correlated semiconductors.
However, despite the enormous efforts of experimental and theoretical work dedicated to understanding its thermoelectric properties, the practical application of FeSb
2 remains severely constrained: the thermoelectric figure of merit of FeSb
2 is extremely low, calculated to be only about 0.005, far below the practical threshold (
ZT ~ 1). The primary reason lies in the exceptionally high lattice thermal conductivity (
κ) of FeSb
2, which also peaks [around 300–550 W/(m·K)] in the same temperature range (~ 15 K) coincidentally [Fig. 3(c)] [
41]. This high thermal conductivity facilitates rapid heat flow through the material, severely compromising its ability to establish and maintain a temperature gradient, thus drastically reducing the thermoelectric conversion efficiency. A proper understanding of this thermoelectric dilemma requires an explicit disentanglement between the diffusion and phonon-drag contributions to the Seebeck coefficient. The colossal low-temperature thermopower in FeSb
2 is widely attributed to a strong phonon-drag contribution associated with long-lived phonons of long mean free path [
58−
61]. This provides a microscopic explanation for the apparent trade-off: the same phonon properties that enhance phonon drag also lead to an unusually large lattice thermal conductivity. As a result, conventional strategies intended to decouple these quantities, including nanostructuring, solid-solution alloying, and defect engineering, have generally produced only limited improvements in
ZT. The underlying issue is an inherent trade-off: phonon scattering introduced to reduce the lattice thermal conductivity also tends to weaken phonon drag, while the associated disorder reduces the exceptionally high carrier mobility that underpins the electrical conductivity of FeSb
2. Consequently, the power factor is often reduced more rapidly than the thermal conductivity. This interplay poses a central challenge for FeSb
2: how to suppress heat transport without significantly compromising the phonon-drag contribution and carrier mobility that together produce the large power factor. The aforementioned disputed experimental phenomena — particularly the characteristic resistivity transition, the magnetic susceptibility peak resembling that of SmB
6/FeSi, and the concomitantly occurring giant thermoelectric effect and high thermal conductivity at low temperatures — collectively form a complex puzzle. Any theory attempting to explain the microscopic mechanism of FeSb
2 must uniformly address a core question: what is the fundamental physical mechanism driving the profound reconstruction of its electronic state? Current discussions primarily revolve around two competing paradigms: i) The Kondo lattice picture views FeSb
2 as a correlated semiconductor or KI. ii) The spin-state excitation picture propose that the spin-state of the Fe ions may changes with temperature. These two theoretical paradigms offer distinct yet self-consistent physical pictures for understanding the exotic properties of FeSb
2. In the following sections, we will first delve into how the Kondo model explains these phenomena, drawing comparisons with typical KI like SmB
6 (Section 3). Subsequently, we will examine the theoretical predictions and experimental evidence for the spin-state excitation model, particularly direct insights provided by recent spectroscopic studies (Section 4), aiming to provide a clear perspective insight of the prolonged debate surrounding the microscopic origin of FeSb
2.
3 The physics of phase transition: An analogy to Kondo insulators
Metal−insulator-like transitions in compounds with partially filled d shells can arise from multiple microscopic mechanisms, including Slater- or Peierls-type instabilities and Mott physics. In FeSb2, however, the emergence of a small low-temperature gap together with a crossover in the magnetic response points to a correlation-driven renormalization of the low-energy electronic structure. These observations have motivated the proposal that FeSb2 may be analogous to Kondo insulators, a framework originally developed for f-electron heavy-fermion systems, where hybridization and many-body coherence produce a narrow gap. In this section, we review and critically assess the evidence for a KI-like interpretation of FeSb2, and compare its key experimental signatures with those of established KI materials to clarify both the advantage and the limitations of this analogy.
3.1 The Kondo insulator analogy: Experimental evidence
The KI model describes materials in which a lattice of localized magnetic moments (typically from rare-earth
f or transition-metal
d orbitals) interacts with itinerant conduction electrons [
62−
68]. At high temperatures, the localized moments behave incoherently, leading to paramagnetic metallic state [Fig. 4(a)]. Upon cooling below a characteristic coherence temperature, the localized moments become screened by the conduction electrons, forming a non-magnetic spin-singlet ground state [Fig. 4(b)] [
69]. This process opens a narrow hybridization gap at the Fermi level, driving a transition to an insulating state [Fig. 4(c)]. Owing to its striking experimental signatures, FeSb
2 has been widely discussed as a candidate for the extension of this theoretical framework to
d-electron-based systems [
52,
70].
Within the KI description, the resistivity evolution with temperature of FeSb
2 can be understood as a sequence of different transport regimes governed by the formation of the hybridization gap. The larger gap Δ
1, observed at higher temperatures (
T > 30 K), is generally associated with the intrinsic band gap of the electronic structure. The smaller gap Δ
2, which dominates the transport at lower temperatures (
T < 30 K), is often identified with the many-body hybridization gap itself. This energy scale reflects the strength of the Kondo hybridization and is regarded as a hallmark of KI-like behavior in FeSb
2 [
70]. Furthermore, while the Kondo model alone cannot account for the large Seebeck coefficient observed in FeSb
2, the hybridization gives rise to sharp resonant features in the density of states near the narrow energy gap. This electronic structure, when strongly coupled to lattice vibrations through a phonon-drag mechanism, is widely considered responsible for the unusually large thermopower [
56]. However, to further assess the applicability of the KI model, it is instructive to compare FeSb
2 with established KI.
3.2 Comparison with other Kondo insulators
While the Kondo model provides a coherent basis for interpreting the magnetic and electric phase transition in FeSb2, direct comparison with well-established KIs such as SmB6 reveals both striking parallels and critical distinctions that complicate a straightforward classification.
SmB
6, the prototypical
f-electron KI with non-trivial topology, manifests Kondo correlation-driven insulating behavior [
71−
77]. It is featured with the coexistence of a correlated insulating bulk states and robust topological metallic surface states that dominate transport at low temperatures, a phenomenon confirmed by Angle-resolved photoemission spectroscopy (ARPES) and various transport measurements [Figs. 5(a, b)] [
72,
73,
78,
79]. Intriguingly, ARPES measurements on FeSb
2 have revealed an electronic band structure surprisingly similar to that of SmB
6, identifying in-gap surface states that cross the Fermi level and a quasi-flat band located beneath the Fermi level through the facet of FeSb
2 (010) [Figs. 5(c, d)]. While the overall phenomenology resembles that of SmB
6, the quasi-flat band in FeSb
2 originates from the weakly dispersive
d-orbital states rather than
f-electron states [
38,
48]. Complementary transport experiments, including the previously discussed Corbino-disk geometry experiments, further support a scenario of an insulating bulk coexisting with a metallic surface channel [
47]. However, both the origin of the gap and the nature of the in-gap states on the FeSb
2 (010) surface remain unsettled, owing to significant inconsistencies among ARPES studies. For example, Xu
et al. [
38] reported pronounced surface-derived features on the (010) surface and interpreted them as evidence for a KI-like hybridization gap with topologically nontrivial surface states [Fig. 5(e)]. In contrast, Chikina
et al. [
48] observed in-gap states with totally different dispersion and Fermiology, which were successfully interpreted as the trivial unbonded surface bands by calculations [Fig. 5(f)]. Taken together, the currently available surface-sensitive ARPES data suggest the surface states are highly sample-dependent, thus do not provide conclusive support for a topological Kondo-insulator picture in FeSb
2. This uncertainty motivates the use of bulk-sensitive probes, such as XAS, to establish the intrinsic electronic and magnetic evolution of the system and to distinguish between competing microscopic scenarios.
Furthermore, the gap observed in FeSb2 (30−40 meV) is substantially larger than the hybridization gap typically reported for prototypical f-electron Kondo insulators such as SmB6 (10−15 meV). Whereas the Kondo gap in SmB6 is closely tied to the coherence temperature and becomes strongly thermally perturbed at temperatures of only a few tens of Kelvin, the gap in FeSb2 remains robust over a much broader temperature range, extending to nearly room temperature. More importantly, if FeSb2 were a conventional KI, its charge and spin responses would be expected to arise from the same underlying hybridization mechanism. However, the experimental results indicate a clear separation between these two channels: the magnetic susceptibility of FeSb2 exhibits a pronounced crossover near 100 K, whereas the electrical resistivity shows no corresponding anomaly at that temperature. This mismatch between the transport and magnetic energy scales is difficult to reconcile with a conventional KI picture.
A more rigorous distinction between
d- and
f-electron systems also requires consideration of the hierarchy between spin-orbit coupling (SOC) and the crystal electric field (CEF). In
f-electron systems, SOC is intrinsically strong because of the large atomic number of rare-earth elements and typically defines the dominant local energy scale. For example, in Ce-based heavy-fermion compounds with a 4
f1 configuration, the SOC splitting is typically on the order of 0.25−0.3 eV, whereas the CEF splitting is usually below 100 meV [
80]. Consequently, the non-magnetic ground state in these materials emerges from the collective many-body Kondo screening process, while CEF effects mainly modify the fine level structure. By contrast, in
d-electron systems such as FeSb
2, the more spatially extended 3d orbitals couple strongly to the ligand environment, making the crystal field the dominant energy scale. For Fe 3
d states in FeSb
2, the effective crystal-field splitting is on the order of 1−2 eV, far exceeding the intrinsic SOC (~ 50 meV). Such a large crystal-field splitting favors a paired, diamagnetic low-spin ground state at the single-ion level, from which thermally activated spin-state excitations can occur.
In light of these fundamental differences, the classification of FeSb
2 as a genuine Kondo insulator remains under active debate. For example, several studies have shown that an insulating ground state can be reproduced from local Coulomb interactions among Fe 3
d orbitals, without invoking Kondo hybridization [
81,
82]. These ongoing debates have prompted the exploration of alternative mechanisms, such as temperature-driven spin-state excitations, which will be discussed in the next section.
4 The spin-state excitation paradigm
As discussed in the previous section, although the KI model could provide a phenomenologically consistant analogy for certain macroscopic observations, whether genuine many-body Kondo coherence exists in FeSb2 remains controversial. These open questions have prompted the exploration of alternative microscopic mechanisms. Among them, the spin-state excitation (SSE) model has emerged as a promising framework, providing a more direct microscopic explanation for the anomalous physical properties of this d-electron system without invoking heavy-fermion hybridization.
SSE has been observed in transition metal compounds, particularly in systems based on 3
d elements with near degenerate spin configurations [Fig. 6(a)] [
83−
86]. In these materials, partially filled
d-orbitals on transition-metal ions (e.g., Fe or Co) may manifest more than one spin state — typically low-spin (LS) and high-spin (HS) configurations — which are nearly degenerate in energy. The ground states are mainly governed by the energy scale of the crystal electric field, which plays more prominent role in an insulator [Fig. 6(a)] [
87]. Transitions between these states can be induced by external stimuli such as temperature, pressure, or magnetic fields, and light, which overcome the multiplet energy gap between the LS/HS. These distinct spin states exhibit different properties: LS state is nonmagnetic (
S = 0) and has smaller effective ionic radius, while HS state possesses finite magnetic moments (
S > 0) and a larger ionic volume [Fig. 6(b)] [
88−
90]. A temperature-driven spin-state crossover thus provides a possible mechanism by which a nonmagnetic, insulating ground state evolves into a paramagnetic and typically more conductive state at higher temperatures.
A prototypical SSE system is LaCoO
3, which crystallizes in a distorted perovskite structure with Co
3+ ions in an octahedral (CoO
6) environment and a 3
d6 electronic configuration, analogous to Fe
2+ (3
d6) in octahedral FeSb
6 units in FeSb
2 [Fig. 6(c)]. It has been studied for decades to understand its SSE induced magnetic transition. LaCoO
3 displays a magnetic and metal-insulator crossover that can be understood in terms of thermally driven transitions from a LS ground state to HS state [Fig. 6(d)] [
90−
94]. Spectroscopic studies, including XAS and inelastic scatterings (RIXS) by Haverkort
et al. [
95], have provided smoking-gun evidence for such transitions. Further work by Takegami
et al. [
89] highlighted the critical role of lattice relaxation in enabling the thermal population of spin-excited states and successfully reproducing the observed temperature dependence.
A similar crystal field scenario may be applicable to FeSb
2, where Fe ions reside in a distorted octahedral environment, similar to the Co ions in LaCoO
3. In such a field, the 3
d orbitals split into lower energy
t2g and higher-energy
eg levels [
87]. When the crystal field splitting energy is comparable to the Hund’s exchange energy of the 3
d-orbitals, the energy difference between the LS and HS would be shrinked, enabling thermally activated transitions between LS and HS configurations. Unlike the KI model, which requires coherent hybridization over extended length scales, the SSE model emphasizes local thermal excitations. This allows for a correlated but localized interpretation of the temperature evolution of both resistivity and magnetic susceptibility, without invoking long-range electronic coherence. Moreover, the pronounced transport anisotropy observed in FeSb
2 — which is difficult to reconcile with the KI model — can emerge in the SSE framework: the crystal-field splitting that determines the LS-HS energy separation is inherently sensitive to structural distortions and bonding geometry [
44,
96]. Consequently, anisotropic thermal lattice distortions naturally lead to axis-dependent transport properties. Thus, the SSE model provides a self-consistent narrative that captures multiple key features of FeSb
2 and merits consideration alongside Kondo descriptions.
Nevertheless, the SSE framework remains under experimental scrutiny. For instance, inelastic neutron scattering (INS) studies by Zaliznyak
et al. [
97] did not observe spin-related excitations within 60 meV corresponding to LS-HS transitions at low temperature, but leaving signals between 15 meV to 30 meV not easy to interpret. High-field Mössbauer spectroscopy on FeSb
2 by Farhan
et al. [
98] have reported the presence of two distinct Fe sites with different hyperfine parameters, which hinted two different spin components from a single valence condition. Moreover, direct spectroscopic evidence confirming temperature dependent spin state transitions in FeSb
2 remains lacking, leaving the actual role of the SSE model in FeSb
2 unconfirmed. The next section will examine available spectroscopic studies of FeSb
2 in greater detail, assessing whether experimental data can help distinguish between competing theoretical pictures.
5 X-ray absorption evidence for spin-state excitations in FeSb2
The transport and magnetic measurements of FeSb
2 can be qualitatively captured by two fundamentally different models: the Kondo hybridization framework and the SSE picture. As outlined in Section 4, previous spectroscopic studies failed to resolve this issue, and in some cases even deepened the debate. In this context, XAS, a technique highly sensitive to both the valence state and local spin configuration of transition metal ions [
99], provides a promising route to probe this problem. A recent study conducted a systematic Fe 3
d L-edge XAS study on FeSb
2, combined with atomic multiplet simulations (AMS), to investigate the spin state configuration of Fe in high quality FeSb
2 single crystals [
100].
By comparing experimental spectra to AMS simulations, it was found that the Fe valence was confirmed to be in a 3
d6 configuration [Fig. 7(a)]. A pure low-spin (LS, S = 0) state successfully captured the overall spectral shape, but failed to account for key pre-edge peak structures observed near the
L3 and
L2 edges. These features could not be attributed to alternative valence states (e.g., 3
d4 or 3
d5), but emerged naturally when a high-spin (HS, S = 2) state component was introduced into the simulation. Similar to the XAS results reported by Sun
et al. [
51] for various iron-based compounds, the data suggest the possible existence of distinct spin states in FeSb
2 [Fig. 7(b)]. A mixed spin state simulation (composed of 85% LS and 15% HS) yielded excellent agreement with the full experimental spectrum, including the pre-edge peaks. This may provide the direct spectroscopic evidence that the ground state of FeSb
2 is not a pure nonmagnetic LS state, but a correlated mixture of LS and HS configurations.
Further evidence comes from temperature dependent XAS experiments. In the process of ramping up the temperature of FeSb2 sample from 20 K to 300 K, the overall spectral profile did not change significantly, but difference spectra (subtracting the 20 K data) reveal systematic evolution [Figs. 8(a, b)]. Specifically, the spectral feature (peak) at the energies associated with HS states (e.g., ~ 705.7 eV) gradually gains the intensity, while the loss of spectral intensity (dip) are observed at LS related features (~ 707.5 eV and ~ 709.4 eV). This peak-dip pattern in the differential spectra provides evidence for a thermal redistribution of spin populations from LS state to HS state. AMS simulations confirm that increasing the HS fraction from approximately 15% to ~20% reproduces these spectral changes quantitatively [Fig. 8(c)]. This result links the macroscopic transition of FeSb2 with a microscopic increase in thermally activated HS occupancy [Fig. 8(d)].
A potential concern is that the properties of topological correlated materials may differ between surface and bulk. Given the surface-sensitive nature of XAS in the total electron yield (TEY) mode, it is necessary to determine whether the SSE features originate from the surface states previously detected by ARPES on (010) facets [
38,
48] or from the intrinsic bulk properties. To address this ambiguity, the study employed ARPES measurement to characterize the (001) facet of FeSb
2. In contrast to the FeSb
2 (010) facet, which hosts metallic surface states, the (001) facet shows no quasiparticle intensity at
E =
EF (Fermi level) [Fig. 9(a)]. ARPES cuts along high symmetry directions reveal a band gap of ~ 30−40 meV, consistent with the activation gap Δ
1 observed in bulk transport measurements [Figs. 9(b−e)]. This established the (001) facet as an ideal, surface-states-free platform for a control experiment, on which temperature-dependent XAS measurements were repeated. The resulting XAS response and its temperature evolution were found to be nearly identical to those measured on the (010) facet. These results strongly indicating that the SSE is intrinsic to the bulk, rather than a surface dominated feature.
The presence of surface states on the (010) facet has previously been taken as evidence supporting a topological Kondo narrative of FeSb2, while the absence of surface state on the (001) facet implicates the trivial topology of FeSb2. As a bulk-surface correspondence, the non-trivial topology would give rise to surface state on any facet that breaks translational symmetry, which is not relying on specific facet. Thus, the identification of surface-state-free (001) facet also unambiguously exclude the assumption of non-trivial topology in FeSb2.
To further explore the tunability of spin states, the effects of doping were studied. RuSb2 is isostructural with FeSb2 and thus offers a way to chemically tune the electronic environment without introducing significant structural disorder. At 20 K, the 5% Ru-doped FeSb2 exhibited stronger pre-edge HS features in the XAS spectrum [Fig. 10(a)]. AMS fitting confirmed an increased HS fraction (~25%), suggesting that Ru doping reduces the LS−HS energy gap [Fig. 10(b)]. These spectroscopic trends to be parallel with macroscopic transport data, which show that Ru doping decreases the activation gap and enhances low-temperature conductivity.
First principles calculation gives an insight to the physical origin of the spin-state transition and the doping effect [Figs. 10(c, d)]. For instance, Density functional theory (DFT)-based calculations confirmed again that LS state and HS state of 3d6 in FeSb2 are energetically almost degenerate. The energy difference of the ground state between the two states is about 7 meV, which is rather comparable to the thermal excitation of ~ 100 K. Moreover, this energy scale is sensitive to the SOC energy. The substitution with a heavy element would eliminate the ground state energy difference. According to the calculations, replacing Fe with Ru can further lowers the HS energy level and promotes HS occupancy, regardless how the Hubbard U or Hund’s coupling J values are tuned, clearly illustrate why HS-state is populated through Ru-doping.
Although first-principles calculations capture the near-degeneracy of spin states and the trend induced by Ru substitution, FeSb
2 has shown intriguing and, at times, contradictory behavior in DFT-based calculations. These discrepancies in the resulting ground states are often taken as indications of strong-correlation or spin-fluctuation effects beyond static mean-field theory, but they can also arise from the extreme sensitivity of narrow-gap systems in first-principles calculations. Specifically, the metallic and magnetic ground state reported in certain studies (e.g., Refs. [
30,
101]) was obtained using the standard PBE-GGA (generalized gradient approximation (GGA) of the Perdew−Burke−Ernzerhof (PBE)) functional. It is well known that PBE systematically underestimates the band gaps of semiconductors [
102−
104]. Because the experimental band gap of FeSb
2 is exceptionally small (30−40 meV), well below the typical uncertainty of standard DFT functionals (~ 0.1−0.2 eV), PBE could artificially close the gap and thus yield a spurious metallic state. By contrast, gap-correcting approaches such as the modified Becke−Johnson (mBJ) potential can recover a small gap; however, the resulting electronic structure remains highly sensitive to small changes in the computational parameters. Because the target gap lies near the limit of computational precision, different methods, or even minor variations in computational parameters, can easily tip the balance and lead to qualitatively different outcomes (metallic versus semiconducting, or magnetic versus non-magnetic). Consequently, although these contrasting theoretical results underscore the complexity of FeSb
2, they should be interpreted with caution, and intrinsic physical mechanisms should be carefully distinguished from the methodological limitations inherent in describing extremely narrow gaps.
In summary, the XAS and ARPES measurement provide essential spectroscopic evidence for spin-state excitation in FeSb2, supporting a ground state comprising a mixture of low spin and high spin Fe2+ configurations. The temperature dependent and Ru-doping dependent evolution of the high spin fraction correlates closely with macroscopic transport trends, reinforcing the relevance of spin-state physics in this system. The experimental and theoretical analysis further suggests that crystal field splitting and Hund’s coupling alone are insufficient to determine the spin state hierarchy. SOC, which energetically favors the high spin state, plays a non-negligible role. For a given SOC strength, the delicate balance between CEF and Hund’s coupling can be thermally tuned when the LS and HS states are nearly degenerate, leading to the observed redistribution of spin populations and transition behaviors.
However, while the AMS simulations captures key spectral features and provides a useful model for interpreting XAS data, it remains a simplified toy-model that does not provide an exact quantitative description of the complex solid environment. Moreover, the small energy scale associated with these spin state transitions presents a significant challenge for theoretical treatment, requiring high precision models capable of resolving fine balance between correlation, spin−orbit interaction, and lattice dynamics. Future work combining high-resolution spectroscopy with advanced theoretical approaches will be highly desired for fully resolving the microscopic origin of FeSb2’s spin state behavior and its role in the material’s anomalous properties.
6 Summary and outlook
Through decades of investigation, FeSb2 has been known a compelling model system at the intersection of electronic correlated physics and functional materials science. This review has examined the material’s unusual transition properties and two competing theoretical interpretations: the KI explanation and the SSE narrative. The accumulated evidence, especially the most recent ARPES and XAS work, clearly substantiate the magnetic transition, as well as the transport anomaly is directly derived from the SSE mechanism. The exclusion of the KI candidate makes the resistivity saturation and 2D transport nature of the cryogenic regime still an open question. However, it is important to emphasize that there is no theorem which excludes the coexistence of the spin-state excitation and Kondo hybridization. While current bulk sensitive evidence firmly establishes SSE as the dominant driver of the anomalous temperature evolution, we cannot fully rule out a coexistence or coupled scenario. It is highly plausible that a Kondo-like coherence or hybridization at low temperatures subtly contributes to stabilizing the low-spin ground state through correlation-driven renormalization. As temperature increases, the suppression of this coherence may in turn facilitate the thermal population of higher spin states. Future investigations may still need to consider whether elements of Kondo interactions play a subtle role in these unresolved low-temperature phenomena. However, as we discussed before, the energy scales of the interactions of CEF and SOC within d-orbital and f-orbital are totally different. Whether or not the single-ion impurity model could extend to a Kondo lattice coherence in the materials of d-orbital basis is very interesting.
We notice that upon doping, FeSb
2 has recently been discussed as a candidate of altermagnet, which would manifest the spin-splitting band structure under the time-reversal symmetry preservation [
101,
105−
107]. We are convinced that understanding the spin configuration, as well as the low-energy electronic band structure is vital in advancing the both new altermagnetic physics and thermoelectricity-oriented applications. To finish this review, we provide a brief perspective on both aspects.
Thermoelectric optimization. One of the core challenges in FeSb
2 research is resolving the so-called “thermoelectric dilemma”: the coexistence of an exceptionally large power factor and a high lattice thermal conductivity, which limits the thermoelectric figure of merit. Various strategies have been attempted to address this issue. Elemental doping (e.g., Bi, Te) has been shown to suppress thermal conductivity but also reduces electrical conductivity and power factor, resulting in limited
ZT enhancement [
108,
109]. Substitutions of both Fe (e.g., Co, Cr) sites have been explored, with limited improvement in enhancing thermoelectric performance [
110]. These reports highlight a general dilemma in doping strategies: reducing lattice thermal conductivity often comes at the expense of carrier mobility and electronic performance. An alternative approach has focused on fabricating the thin films devices. For instance, magnetron sputtered FeSb
2 films with controlled Sb vacancies show reduced thermal conductivity but inferior power factor compared to single crystals [
111]; MBE-grown films allow exploration of phonon drag effects and heterostructure engineering, enhancing the thermopower to 208 μV/K [
112]. Other efforts focus on defect engineering in single crystals to enhance phonon scattering while minimizing disruption to electronic transport [
59−
61,
113−
116]. Although
ZT value of 0.071 has been achieved in FeSb
2 at 55 K [
111], the challenge for practical application remains significant. These efforts collectively involve different efforts to tune the thermoelectricity, including chemically doping, alloying, or nanostructuring (e.g., thin films, nanowires, composites) [
109,
117−
123] aiming at selectively scattering phonons without degrading carrier mobility. In parallel, theoretical studies and high throughput screening may aid in identifying optimal doping or design strategies. A deeper understanding of the phonon drag mechanism and its coupling to spin and electronic structure will be essential [
58,
124−
127].
Altermagnetism. Altermagnetism is a recently proposed class of materials characterized by momentum-dependent spin-splitting despite vanishing net magnetization, attracting attention for their potential applications in spintronics [
128−
135]. Although FeSb
2 itself is nonmagnetic, its distinctive crystal symmetry is potential to break the
PT-symmetry to host the exotic magnetic order. Recent theoretical works repeatly proposed that doping FeSb
2 (e.g., with Co) could induce AFM order and derive an altermagnetic phase in this system [
101]. Experimentally verifying this prediction presents a compelling research direction, requiring a combination of magnetization measurements and advanced spectroscopic techniques (such as spin-resolved ARPES or XMCD) to map the magnetic phase diagram and spin textures of doped FeSb
2. Confirming altermagnetism in this system would not only expand the family of known altermagnetic materials, but also position FeSb
2 as a potential platform for next-generation spintronic devices.
In conclusion, FeSb2 continues to stand out as a unique platform linking narrow-gap semiconductors, strongly correlated electron systems, and high-performance thermoelectric materials. Deepening the understanding of the interplay between electronic correlations, spin states, lattice dynamics, and potential exotic magnetic order in this material may not only resolve its underlying physical puzzles, but also push forward the frontier of research for next-generation materials with enhanced energy-conversion efficiency and spintronic functionalities.