Magnetism in Na-filled Fe-based skutterudites.

The interplay of superconductivity and magnetism is a subject of ongoing interest, stimulated most recently by the discovery of Fe-based superconductivity and the recognition that spin-fluctuations near a magnetic quantum critical point may provide an explanation for the superconductivity and the order parameter. Here we investigate magnetism in the Na filled Fe-based skutterudites using first principles calculations. NaFe4Sb12 is a known ferromagnet near a quantum critical point. We find a ferromagnetic metallic state for this compound driven by a Stoner type instability, consistent with prior work. In accord with prior work, the magnetization is overestimated, as expected for a material near an itinerant ferromagnetic quantum critical point. NaFe4P12 also shows a ferromagnetic instability at the density functional level, but this instability is much weaker than that of NaFe4Sb12, possibly placing it on the paramagnetic side of the quantum critical point. NaFe4As12 shows intermediate behavior. We also present results for skutterudite FeSb3, which is a metastable phase that has been reported in thin film form.

An important feature of the Fe-based superconductors and the related antiferromagnetic phases is that the energy scale for the spin fluctuations is quite high, for example with spin wave dispersions extending well over 100 meV, hardly indicative of materials in which the magnetism is readily suppressed1234567891011. Nonetheless, this is the case and in particular the magnetic orderings can generally be readily suppressed in favor of superconductivity by doping, isovalent alloying and pressure.

One scenario is that the ordering is suppressed by quantum fluctuations. In this case, at the mean field level the compounds would be strongly magnetic over the entire range of the phase diagrams, similar to what has been discussed near magnetic quantum critical points1213. For example, standard density functional theory (DFT) calculations (done with the PBE generalized gradient approximation and the experimental crystal structure) predict ordered magnetism with moments of ~2 μ over the entire composition range of LaFeAs(O,F)14, while in fact the ordered moment is ~0.5 μ for LaFeAsO and is rapidly suppressed by fluorine doping89. This type of overestimate of magnetic tendencies within density functional calculations is unusual. It typically arises when spin fluctuations associated with a nearby quantum critical point are strong enough to renormalize the mean-field like magnetic state predicted by standard approximate density functional calculations1516. Thus these results suggest the presence of strong renormalization by quantum fluctuations in the Fe-based superconductors. These would presumably include longitudinal fluctuations17, recently observed in BaFe2As218. This is supported by experimental measurements for example the observation by xray-absorption spectroscopy of a large exchange splitting of the Fe 3s state in superconducting CeFeAs(O,F) as if the compound were magnetically ordered, even though it is not magnetic at the composition studied19. This is qualitatively similar to but stronger than what is seen in NbFe2, which is known to be close to a quantum critical point202122.

Turning to the chemistry of the Fe-based superconductors, all the phases have layers of Fe in a formal divalent or close to divalent state and tetrahedral coordination by the pnictogens, P or As, or Se or alloys with Te and/or S. There has been theoretical work suggesting stronger magnetism and possibly superconductivity in hypothetical Fe antimonide compounds analogous to the superconducting arsenides23, but to date these compounds have not been synthesized. While there are many other families of Fe compounds that show magnetism, and in some cases the magnetism can be suppressed by alloying or pressure, in general these other compounds do not display the interesting magnetic properties of the Fe-based superconductors. Nonetheless, it may be that the strong renormalizations present in the Fe-based superconductors may not be unique to those compounds and it is of considerable interest to find other iron compounds displaying such behavior.

Recently, calculations suggest that such a suppression of magnetic order by quantum fluctuations with qualitatively similar size to the Fe-based superconductors can occur in the Ge compound, YFe2Ge22425. This is a compound that shows other indications of nearness to a magnetic quantum critical point, and also possibly low temperature superconductivity2627. The electronic structure of this compound is considerably more three dimensional than any of the Fe-based superconductors. It is therefore of interest to search for other compounds that might show related physics.

Filled skutterudites are a large family of transition metal pnictide compounds that have been extensively studied in the context of thermoelectrics, but also show other interesting behaviors including Fe-based magnetism in NaFe4Sb12282930, and both conventional (e.g. LaRu4P12)31 and unconventional superconductivity (e.g. PrOs4Sb12)3233.

Depending on the filling atom, the filled Fe-based skutterudites form with all of P, As and importantly Sb (e.g. LaFe4P12, LaFe4As12, LaFe4Sb12). This is unlike the Fe-based superconducting structures, for which antimonides do not exist. This provides an opportunity for studying Fe-pnictides near magnetism in compounds that allow comparison of antimonides with arsenides and phosphides. However, also different from the Fe-based superconductors, which are based on nominally divalent Fe tetrahedrally coordinated by pnictogens, filled skutterudites have a more covalent Fe with octahedral coordination. Further, the filled skutterudites are cubic, while the Fe-based superconductors have layered structures. This is an important difference from the point of view of quantum fluctuations, since dimensionality it thought to play an important role in the behavior of quantum critical systems34. Nonetheless, these filled skutterudites do show a composition dependent cross-over between a paramagnetic metal to a magnetically ordered state with signatures of itinerant character and perhaps sizable effects due to quantum fluctuations in the cross-over region353637. These include transport signatures, highly enhanced specific heats, strong orderings induced by rare-earth moments, and magnetotransport signatures.

The purpose of this paper is to elucidate the chemical range of the magnetism in these compounds and specifically the extension from the well studied antimonides to the arsenides and phosphides. We also report results for the binary skutterudite, FeSb3 motivated by recent experimental results38, in relation to reported DFT calculations39.

Results

The band structures (Fig. 1) show rough similarity between the La-filled and Na-filled compounds within a family with the exception of the placement of the Fermi energy. The filled skutterudite band structures show a light pnictogen derived band. This light band comes from the band that crosses the gap and forms the top valence band in the unfilled compounds4041. The lower bands are heavier and have substantial Fe 3d contributions. It is these heavy Fe 3d derived bands that lead to magnetism, as was discussed previously283037. The different families (P, As and Sb) differ in the shape and placement of the light band with respect to the heavier bands and in the details of the hybridization of the Fe and pnictogen states. This is also seen in the density of states (Fig. 2). However, the details of the band structures of the Na and La compounds are not simply related by rigid band behavior. This is because the position of the light band is influenced by the charge on the filling atom42, and in particular the light band is lower in energy, closer to the heavy bands, in the La filled compounds, corresponding to the higher charge on trivalent La. This then affects the shape of the density of states in a way that can work against the effect of the change in band filling as regards magnetism (see below).

Figure 1

Band structures of the various compounds. The Fermi energies are at 0 eV. Note that in all cases the Na compounds have Fermi levels cutting heavy transition element derived bands.

Figure 2

Total (top) and Fe-d projections (bottom) of the density of states of NaFe4P12, NaFe4As12 and NaFe4Sb12. The Fe-d projections are onto LAPW sphere radii of 2.1 Bohr for Fe. The Fermi energy is at 0 eV. The dotted line in the bottom panel marks 2 eV−1 (see text).

The relative placement of the light and heavy bands plays a central role in the thermoelectric performance of p-type filled skutterudites41, which has been one of the main interests in the study of filled skutterudites. In this regard, the light band extends higher above the heavy bands in the phosphide skutterudites, which makes them less favorable as thermoelectrics43. Perhaps for this reason the phosphide and arsenide filled skutterudites have been less studied as thermoelectrics than the antimonides. This also has consequences for magnetism.

Within Stoner theory444546, an instability towards ferromagnetism occurs in the mean field when the Stoner criterion N(E) ≥ I−1 is met. Here N(E) is the electronic density of states at the Fermi energy, E, and I is the Stoner parameter, which depends on the chemical details, but for an Fe pnictide compound is governed by the Fe contribution because the pnictogen p orbitals are very spatially extended, and therefore have a small Hund’s coupling. Thus the criterion is roughly , where N(E)(d) is the Fe 3d contribution to N(E) on a per spin and per atom basis, and I, which measures the Hunds coupling on Fe, is ~1 eV4748.

Therefore a magnetic instability at mean field can be expected when the Fe d projected density of states at E significantly exceeds 2 eV−1 (i.e. 1 eV−1 per spin). This threshold is indicated by the dotted line in the density of states projection (lower panel of Fig. 2). The size of the magnetic instability (i.e. the moments) is controlled by the width of the energy region in the density of states, N(E) over which N(E) is large4546. Turning to the shape of the density of states, there are two competing effects when the charge on the filling atom is reduced going from La to Na. First of all, E is lowered which raises N(E) due to the heavy bands. The second is that the light band rises relative to the heavy bands, which works against this.

The calculated values of N(E) are 26.3 eV−1, 43.8 eV−1 and 29.9 eV−1 on a per formula unit basis for the phosphide, arsenide and antimonide, respectively. The corresponding contributions from Fe d states as measured by the d projection onto an Fe LAPW sphere are 4.1 eV−1, 6.7 eV−1 and 3.5 eV−1 per atom. These values are clearly above the Stoner criterion, and therefore all these Na-filled skutterudites are expected to have ferromagnetic instabilities at the mean field-like DFT level. This is reflected in our self consistent calculations, which show ferromagnetic ground states at the DFT level for all compounds.

We emphasize that in general, while the presence of a mean field ferromagnetic instability is governed by the value of N(E) in Stoner theory, once there is an instability, its size is not determined by the value of N(E) but on the shape of the density of states in the energy range around E, as described by extended Stoner theory4546. In particular, a larger instability often results from a broader density of states peak, even though the peak N(E) is lower.

Fig. 3 shows the results of fixed spin moment calculations. Fixed spin moment calculations for NaFe4Sb12 were reported previously by Leithe-Jaspers and co-workers28. As shown, the Na-filled skutterudites are all ferromagnetic at mean field. However, the Sb compound is much more ferromagnetic in terms of both the magnetic energy and the size of the predicted moment (3.0 μ per formula unit for NaFe4Sb12, close to 3.0 μ for NaFe4As12, and 1.6 μ for NaFe4P12). The As compound is intermediate, but much more similar to the antimonide than to the phosphide. We also did unconstrained calculations including spin-orbit for NaFe4Sb12 and NaFe4P12, and obtained very similar results, specifically spin moments of 2.94 μ and 1.62 μ, respectively.

Figure 3

Fixed spin moment total energy as a function of constrained spin magnetization as obtained using the PBE GGA for NaFe4P12, NaFe4As12 and NaFe4Sb12 on a per formula unit basis. The energy zero is at the energy of the non-spin-polarized case.

Although hardly unusual, it is worth noting that the magnetic energies and magnetizations do not correlate with the size of N(E). The reason for the differences is seen in the band structures and corresponding density of states, N(E). As mentioned, the light band is higher in the phosphides than in the arsenide and antimonide filled skutterudites. This leads to a lower N(E) in the phosphides as the energy is raised above E along with a wider region of low N(E) near the band edge (3 electrons per formula unit above E for the Na-filled compounds). Based on this analysis, NaFe4P12 is expected to be ferromagnetic, but with weaker ferromagnetism than NaFe4As12 and NaFe4Sb12, which is what is found in the full calculations. We also did local spin density approximation (LSDA) calculations for NaFe4P12. We find again a ferromagnetic ground state (spin moment of 1.3 μ per formula unit), supporting the conclusion from the PBE GGA calculations that this compound has a ferromagnetic ground state at the mean field like DFT level.

Mochel and co-workers38, have reported experimental investigation of the transport and susceptibility of skutterudite structure FeSb3 films. The authors extracted a band gap of 16.3 meV from a fit of the resistivity. The measured resistivity was found to decrease by ~20% with temperature from ~30 K to 300 K. The temperature dependent resistivity data does not, however, properly fit the activated form that would be expected for a semiconductor with a defined gap over any temperature range and furthermore has a maximum at ~30 K and decreases as temperature is lowered below this. The susceptibility does not show any magnetic transition as a function of temperature, but instead nearly diverges as the temperature is lowered, showing Curie-Weiss behavior with a Weiss temperature close to 0 K. These behaviors could be consistent with a material very close to a ferromagnetic instability.

Rasander and co-workers39, reported detailed lattice dynamics and electronic structure calculations for this phase, finding that the compound has a ferromagnetic ground state with a high estimated Curie temperature T = 175 K, and that it is only dynamically stable when ferromagnetism is included. The necessity of including magnetism in order to accurately describe the crystal structure even in the non-magnetic state is a characteristic also found in the Fe-based superconductors14.

We repeated these calculations to confirm this potentially important result in relation to quantum criticality in skutterudites. We did calculations with both the PBE GGA and with the LSDA, which sometimes gives weaker magnetic tendencies than the GGA. We find ferromagnetic instabilities in both cases, which supports the findings of Rasander and co-workers39. Specifically, for the optimized PBE structures we obtain a magnetization of 4 μ per unit cell (four formula units of FeSb3) with a magnetic energy of 280 meV. In the LSDA, relaxing the lattice parameter according to the LSDA, we obtain a somewhat smaller magnetization of 3.55 μ per unit cell, with an energy gain of 50 meV. This is a severe test because the LSDA volume is 6.5% smaller than the PBE volume and this compression works against magnetism (the calculated lattice parameters for the ferromagnetic state are 9.175 Å and 8.971 Å for the PBE GGA and the LSDA, respectively, while the corresponding non-spin-polarized lattice parameters are 9.151 Å and 8.960 Å). Thus both the LSDA and PBE functionals predict a ferromagnetic ground state in contrast to the experiment, again suggesting further investigation of this material for further signatures of quantum critical behavior.

Discussion

To summarize, we find that at the standard density functional level (PBE GGA), NaFe4Sb12 is ferromagnetic in accord with experiment and previously reported calculations. As mentioned, based on experimental data, this compound shows significant effects of quantum fluctuations. In general, such fluctuations weaken the magnetic instability. This effect has been observed in NaFe4Sb12 and KFe4Sb12, where ordered moments of ~1 μ per formula unit are found experimentally instead of 3 μ predicted by density functional calculations2829. This is reminiscent of other itinerant metals near a ferromagnetic quantum critical point.

The resulting renormalization is controlled by the magnetic energy surface (the bare energy as a function of constrained moment) and the size of the quantum fluctuations in the order parameter, which can vary between different materials124950. Following the arguments of Aguayo and co-workers for Ni3Al in relation to Ni3Ga, the lighter bands in NaFe4P12 relative to NaFe4Sb12 might favor weaker quantum fluctuations in the phosphide. In any case, the present results imply that NaFe4P12 is either a ferromagnet or a renormalized paramagnet in which ferromagnetism has been destroyed by spin fluctuations. This is also the case for the hypothetical compound NaFe4As12, which shows stronger magnetism at the PBE GGA level, closer to the behavior of NaFe4Sb12.

It will be of interest to experimentally study the magnetic properties of NaFe4P12 to determine if the compound has ordered ferromagnetism. In either case, transport, specific heat and magnetotransport measurements looking for signatures of quantum fluctuations will be useful. This also applies to the arsenide, NaFe4As12 if it can be made, and to FeSb3, which exists in thin film form.

Methods

The present calculations for the filled skutterudites were performed within density functional theory using the generalized gradient approximation of Perdew and coworkers (PBE GGA)51. We used the experimental lattice parameters of the filled skutterudites, NaFe4P12 and NaFe4Sb12 and relaxed the atomic position of the pnictogen, which has two free structural parameters. Skutterudite structure NaFe4As12 has not been reported in literature. For this hypothetical compound we relaxed both the lattice parameter and the internal As position using the PBE GGA. The calculated lattice parameter for this phase is 8.345 Å. The structural data is given in Table 1.

Table 1

Structures and calculated N(E ) per formula unit for the various compounds. u and v are the free parameters in the pnictogen position (u,v,0).

Compound	a(Å)	u	v	N(EF)(eV−1)
LaFe4P12	7.8316	0.3531	0.1522	14.8
LaFe4As12	8.2152	0.3443	0.1570	14.7
LaFe4Sb12	9.1395	0.3360	0.1633	19.4
NaFe4P12	7.782	0.3530	0.1484	26.3
NaFe4As12	8.345	0.3445	0.1546	43.8
NaFe4Sb12	9.1767	0.3328	0.1620	29.9

We used the general potential linearized augmented planewave (LAPW) method52 as implemented in the WIEN2k code for most of the calculations reported here53. We employed well converged basis sets including local orbitals, with the parameter Rk = 9 (R is the minimum LAPW sphere radius and k is the cutoff for the planewave sector of the basis). We used the precise standard LAPW plus local orbital approach to treat semicore states, as opposed to the so called APW+lo method5455. Dense zone samplings were used, generally 24 × 24 × 24 meshes, and we tested convergence by varying these. We also did calculations with different sphere radii as tests and also repeated some calculations using the projector augmented wave (PAW) method5657. as implemented in the VASP code with standard settings5859. The calculations for skutterudite FeSb3 were done using the VASP code.

Additional Information

How to cite this article: Xing, G. et al. Magnetism in Na-filled Fe-based skutterudites. Sci. Rep. 5, 10782; doi: 10.1038/srep10782 (2015).