Site specific spin dynamics in BaFe2As2: tuning the ground state by orbital differentiation.

The role of orbital differentiation on the emergence of superconductivity in the Fe-based superconductors remains an open question to the scientific community. In this investigation, we employ a suitable microscopic spin probe technique, namely Electron Spin Resonance (ESR), to investigate this issue on selected chemically substituted BaFe2As2 single crystals. As the spin-density wave (SDW) phase is suppressed, we observe a clear increase of the Fe 3d bands anisotropy along with their localization at the FeAs plane. Such an increase of the planar orbital content is interestingly independent of the chemical substitution responsible for suppressing the SDW phase. As a consequence, the magnetic fluctuations in combination with this particular symmetry of the Fe 3d bands are propitious ingredients for the emergence of superconductivity in this class of materials.

The structural parameters in low-symmetry layered systems, such as the iron pnictides/chalcogenides, cuprates and the heavy fermion (HF) compounds have played an important role in determining both symmetry and dimensionality of the magnetic fluctuations123. An immediate consequence of a structural control (or tuning) parameter is its influence on the crystalline electrical field (CEF) effects, which are often relevant in determining the ground state in strongly correlated materials. Particularly for the HF superconductors, the f-electrons strongly hybridize with conduction electron (ce) bands and the single-ion anisotropy, mostly defined by the CEF effects, influences the magnetic fluctuations at the Fermi surface (FS)45678. In fact, the superconducting transition temperature (T) scales with the lattice parameter ratio c/a for several HF systems48. Yet for the high-T cuprates, some structural parameters, such as the Cu-apical oxygen distance and the bond length between Cu and the in-plane oxygen, have been suggested as control parameters of T. Remarkably, the apical oxygen distance to the CuO2 plane also displays a reasonable scaling with T910. In the same fashion, the Fe-pnictide distance has been recently suggested as a similar control parameter tuning the CEF levels in the Fe-based compounds and, consequently, changing the orbital contributions to the bands at the Fermi level1112. Moreover, it is notorious that this orbital interplay is indeed crucial to the formation of the spin-density wave (SDW) phase due to the symmetry break (or xz/yz orbitals). Nonetheless, the ultimate role of the orbital degrees of freedom on the emergence of supercondutivity in Fe-pnictides materials remains an unsolved puzzle.

In this regard, several spin probes such as electron spin resonance (ESR), nuclear magnetic resonance (NMR), muon spin rotation (µSR), Mössbauer spectroscopy, and magnetic neutron scattering, can in principle be used to explore the CEF and orbital differentiation effects on the spin dynamics of a system and their possible influence over its superconducting order parameter. The observation of homogeneous resonances by such techniques in strongly correlated metals allows one to obtain the imaginary part of the spin susceptibility χ″ (q, ω). In particular, the ESR absorption line of a paramagnetic local probe is a measure of χ″ (q, ω), which in turn is directly related to the spin-lattice relaxation rate (1/T1) of the localized resonating spins through the ce and thus to the lattice1314. In a general approach for dilute local moments in metals, the ESR linewidth is linear in temperature and can be written as ΔH = a + bT, the so-called Korringa relaxation15. The residual linewidth, a, arises from the spin-spin interactions of local moments, lattice defects, g-factor distribution, and perhaps unresolved hyperfine and crystalline electrical field (CEF) fine structures. In addition, the thermal broadening b is defined as: Here, 〈J(q)2〉1/2 is the effective exchange interaction between the local moment and the ce in the presence of ce momentum transfer averaged over the whole FS16, η(E) is the “bare” density of states (DOS) for one spin direction at the Fermi level, k is the Boltzmann constant, µ is the Bohr magneton, g is the local moment g-value and K(α) is the Korringa exchange enhancement factor due to electron-electron exchange interaction1920. Here ‘bottleneck”, “dynamic” effects, and multiple bands effects are not present21. These are reasonable assumptions because when “dynamic” effects are present the g-values are usually strongly T-dependent and when the bottleneck effect is relevant the Korringa rate b decreases with an increasing concentration of the magnetic ions21. Concerning the case of multiple band effects, one would expect the calculated g-shift to be much smaller than the measured one222324. As an example, we consider the pure EuFe2As2 compound in the simplest scenario for the treatment of exchange interaction between Eu2+ spins and conduction electron spins, i.e., b = (π.k/g.µ)Δg. By replacing Δg by its experimental value Δg ~ 0.05 one gets a Korringa rate of b ~ 150 Oe/K, which is much larger than the measured one. This result is not only valid for the pure compound but also for the samples studied in the present work, which present the same g-shift. Thus, none of the above mentioned effects have been observed in the studied temperature range in previous works26, nor in our present ESR data. In fact, a smaller experimental Korringa rate implies that one needs to take into account q-dependent exchange interaction and perhaps conduction electron-electron correlations as the dominant effects, which were already included in Eq. (1) of the present manuscript. The g-shift Δg can thus be written as25: At the local moment site, the g-shift probes the ce polarization (q = 0). On the other hand, the Korringa rate probes the ce momentum transfer (0 ≤ q ≤ 2k) averaged over the FS16. Thus, from Eqs. 1 and 2 one obtains information about the density of states at the Fermi level as well as the strength and q-dependence of the exchange interaction.

Previous studies on the series Ba1−EuFe2As2 have shown that decreases with decreasing Eu content26. The scaling of with T suggests that the electron bands with appreciable overlap with the Eu2+ 4f states become more anisotropic (less s-like) as the SDW phase is suppressed. Moreover, these electron bands also move, on average, further away from the Eu2+ sites in real space26. Therefore, in order to generalize such a scenario and further investigate the role of symmetry changes in highly occupied 3d bands in superconducting members of this family, we perform in the present report a systematic ESR study on Ba1−EuFe2−MAs2 single crystals (with M standing for transition metals Co, Cu, Mn, Ni, and Ru). It is worth pointing out that the Ba and FeAs planes are probed by the ESR active Eu2+ and Mn2+/Cu2+ paramagnetic spins, respectively. The (x = 0.2; y = 0) compound was taken as reference for the evolution of the spin dynamics as a function of M substitution (y content) because it is the lowest Eu concentration presenting Korringa relaxation.

Starting off with the macroscopic physical properties, Fig. 1a displays the temperature dependence of the normalized electrical resistivity for the studied single crystals. Room-T values of ρ(T) vary within the range 0.2–0.8 mΩ.cm. Except for M = Ru, which displays the smallest substitution content (y = 0.01), the electrical resistivity increases at T, as typically found for doped BaFe2As2 samples27. Besides, for the chosen reference compound (x = 0.2; y = 0), a metallic behaviour is observed down to T where a sudden drop at 137 K is identified in its dataset. A suppression of the SDW phase with lower transition temperatures of 135 K, 102 K, 98 K, and 86 K is clearly observed as Fe is replaced by transition metals M = Ru, Ni, Cu, and Co, respectively. Superconductivity also emerges for Co and Ni substitutions with T = 22 K and 6 K, respectively (inset to Fig. 1a). In good agreement, the SDW phase is also suppressed for BaFe2−MAs2 (i.e., x = 0) in comparison with the parent BaFe2As2 with .This can be noted for the y = 0.1 compound with T = 78 K. In particular, the SDW phase is fully suppressed for y = 0.1 giving way to superconductivity at T = 3.8 K (inset to Fig. 1a).

Figure 1

Thermodynamic properties of the Ba1−EuFe2−MAs2 single crystals.

The insets show a) the zero electrical resistivity at the superconducting transition and, b) zero-field cooling/field cooling magnetic susceptibilities.

Fig. 1b shows the in plane magnetic susceptibility as a function of temperature for H = 1 kOe. For all samples, χ(T) is well fit to the Curie-Weiss law plus a T-independent Pauli term, χ(T) = χ0 + C/(T − θ) (solid lines). We obtain an effective moment µ ≈ 8µ from the Eu2+ ions and an electronic spin susceptibility of χ0 = 2(1) × 10−3 emu/mol-Oe for all compounds. The superconducting transition at low-field (H = 20 Oe) is defined by the arrows in the inset of Fig. 1b. It is worth noting that, for M = Co and Ni, the diamagnetic response below T is overcome by the Eu2+ paramagnetic one. Fig. 1c shows the T-dependence of the specific heat for the selected single crystals. For the reference compound with x = 0.2 and y = 0, the SDW transition is defined by a jump at 137 K (green arrow). In agreement with the ρ data of Figs. 1a–1b, we observed a suppression of the SDW phase as Fe is substituted by transition metals M.

Now we turn our attention to the ESR results on the selected single crystals in the paramagnetic (PM) state (T > TSDW). A single asymmetric Dysonian ESR resonance is observed for all compounds, meaning that the skin depth is smaller than the sample dimensions28. Nonetheless, as presented in Fig. 2, the X-Band () ESR spectra at room-T become Lorentzian-like after gently crushing the single crystals in order to improve the signal to noise ratio. From fits to the resonances using the appropriate admixture of absorption and dispersion (solid lines), we obtained the T-dependence of both the linewidth ΔH and g-values presented in Fig. 3.

Figure 2

X-Band (~ 9.5 GHz) ESR spectra of Ba1−EuFe2−MAs2 single crystals at room temperature.

Figure 3

T-dependence of a-b) the ESR linewidth ΔH and c) the g-values measured for Ba1−EuFe2−MAs2 compounds.

A linear (Korringa) increase of ΔH with increasing-T is observed for the Eu2+ ESR signal in the PM state, as shown in Fig. 3a. From linear fits to ΔH(T) above TSDW we extracted the values of the Korringa rate b ≡ ΔH/ΔT. It is evident that b decreases systematically for any of the transition metal substitution. On the other hand, Fig. 3b shows that, when x = 0 and M = Mn2+ or Cu2+, the ESR signal from these magnetic probes broadens as the temperature is lowered indicating that a competition between the spin-spin and spin-ce interactions must be occurring. However, in the high-T range, the Korringa-like relaxation dominates for the y = 0.1 samples. We will discuss the behavior of the dilute samples (y = 0.01) below. Thus, linear fits to ΔH(T) allowed us to obtain Oe/K and Oe/K, respectively.

The T-dependence of the Eu2+ ESR g-value in the normal state is presented in Fig. 3c. Before the long range SDW order sets in, a roughly T-independent is commonly observed for the Eu2+ ESR signal, as previously suggested by reports on EuFe2As2, Eu0.5K0.5Fe2As2 and EuFe2−CoAs2323334. In fact, the high-T Eu2+ g-values of our Ba1−EuFe2−MAs2 single crystals are, on average, equal to 2.04(4) and independent of the Eu2+ content x. For the narrowest ESR linewidths, we are able to determine the g-value with higher precision, g = 2.05(3). As it can be verified in Fig. 3c, such a value is indeed a reasonable mean g-value for all x with T > TSDW.

Based on the framework discussed in Ref. 26, we now show that possible changes in the DOS as well as any evolution of the electron-electron interaction cannot account for the Korringa suppression observed as a function of M substitution. In order to evaluate the possible contribution of e-e correlations in our ESR data, we note that the electronic spin susceptibility obtained experimentally is χ0 = 2(1) × 10−3 emu/mol-Oe for all studied compounds. This value is about one order of magnitude larger than the magnetic Pauli susceptibility estimated from the heat capacity coefficient γ, indicating that an e-e exchange enhancement contributes to the ce spin susceptibility in the series. Thus, in the presence of such an enhancement, the host metal ce spin susceptibility can be approximated by: where α accounts for the e-e interaction, (1 - α)−1 is the Stoner enhancement factor and η (E) is the “bare” DOS for one spin direction at E3031. Now, assuming that the enhancement in χ0 is only due to e-e interaction, an α ≈ 0.85(5) can be estimated. Consequently, based on Ref. 20 one can determine the corresponding K(α) = 0.2(1). Then, using the relation (πk/gµ) = 2.34 × 104 Oe/K and replacing , η(E) = 3.34 states/eV mol-spin and b values from Table 1 into Eqs. 1 and 2, we are were able to extract J(0) = 2(1) meV and for all compounds. These values are also presented in Table 1.

Table 1

Experimental ESR Korringa rate b, the g-values, the g-shift Δg, the exchange interaction J(0) and calculated for all samples studied in this work. As defined in Eq. 2, Δg ≡ g − g. For Eu2+, g is well known to be 1.993(1)13. For Mn2+, an average of g value in similar insulating compounds leads to 2.00(5)13. For Cu2+, we assume that the magnetic moment is entirely due to the spin and the g ~ 2. Although this assumption leads to a significant error in J(0), it does not affect the Korringa-like analysis (Eq. 1)161718

Sample	b (Oe/K)	g-value	Δg	Jf s(0) meV
Ba0.8Eu0.2Fe2As2 (ref. compound)	4.3(2)	2.04(4)	0.05(3)	2(1)	1.4(8)
BaFe1.9Cu0.1As2	3.5(2)	2.08(4)	0.08(3)	2(1)	1.2(6)
Ba0.78Eu0.22Fe1.91Cu0.09As2	2.8(2)	2.05(4)	0.06(4)	2(1)	1.1(8)
Ba0.8Eu0.2Fe1.89Ni0.11As2	2.5(2)	2.04(4)	0.05(3)	2(1)	1.0(8)
Ba0.8Eu0.2Fe1.88Co0.12As2	1.7(2)	2.05(3)	0.06(3)	2(1)	0.9(7)
Ba0.75Eu0.25Fe1.99Ru0.01As2	1.5(2)	2.04(3)	0.05(3)	2(1)	0.8(7)
BaFe1.88Mn0.12As2	1.1(2)	2.04(3)	0.05(3)	2(1)	0.7(6)

Based on the data of Figs. 1–3, the diagrams shown in Fig. 4 summarize the behaviour of TSDW and Korringa relaxation rate b against M substitution. When the ESR spin probe is out of the Fe-As plane (Fig. 4a), we observe that both TSDW and b values are suppressed with respect to their values for Ba0.8Eu0.2Fe2As2 as Fe is substituted by M. Nevertheless, when the spin probe is in-plane (Fig. 4b), the opposite behaviour takes place and the Korringa rate b increases as the SDW phase is suppressed. Interestingly, for smaller Mn and Cu concentrations of the order of ~ 0.5%, the ESR linewidth broadens with decreasing temperature in the entire temperature range. In particular, Fig. 3b includes the thermal broadening for the y = 0.01 sample. In this case, the Korringa rate, if it exists, would be at least smaller than the thermal broadening observed, otherwise it would be detectable. Thus, one can infer that the Korringa rate for the Cu sample is indeed decreasing for smaller concentrations. On the other hand, it is not straightforward to infer such conclusion from ESR data on lower concentrations of Mn since the Korringa rate for higher concentrations is already very small (~1 Oe/K). Nevertheless, for both in-plane probes, the Korringa behavior is overcome by the spin-spin interactions associated with the proximity to a magnetic SDW phase. Furthermore, it is worth noting that 3d ions in low symmetry metallic hosts are most likely to display quench of the angular momentum and thus spin-only wave functions. In fact, our experimental data show that the g-factor is approximately 2 for both Mn and Cu ESR lines. Therefore, it is reasonable to assume that the overlap with the conduction electrons will not differ dramatically between Mn and Cu and in turn their Korringa rates for 5% concentrations could be compared.

Figure 4

Spin-density wave transition temperature TSDW and Korringa rate b as a function of transition metal substitution M when the ESR spin probe is a) out of the Fe-As plane and, b) in the Fe-As plane.

Therefore, these results combined evidence that, when the probe is out of (in) the Fe-As plane, the relative value of clearly decreases (increases) with M-substitution. Given that J(q) is the Fourier transform of the spatially varying exchange interaction, by suppressing the SDW phase by M substitution, the electron bands with appreciable overlap with the Eu2+ 4f states become more anisotropic (less s-like). Moreover, these bands, on average, move further away from the Eu2+ sites in real space, i.e., they assume a more planar/xy-orbital character. In fact, in Fig. 4a one can clearly observe a scaling between the magnitude of the SDW phase suppression and that of the Korringa rate, except for the Ru substitution. This is because the particularities in the distortion of the electronic bands may behave differently for each transition metal. For instance, if the chemical substitution is not coherent, interference can occur among the 3d bands and consequently the suppression of the Korringa rate b might be more effective than that of TSDW29.

It is imperative to notice that our findings are in complete agreement with the data on the concentrated Eu regime. A slower Korringa rate b has also been observed when K or Co were introduced into EuFe2As23334. Furthermore, DFT+DMFT calculations have recently shown that orbital differentiation is crucial when magnetic correlations become relevant in BaFe2As212. By decreasing the iron-pnictogen distances, the xy orbital occupancy increases. This in turn weakens the Fe magnetic moment.

Furthermore, this change in occupancy also affects the Fermi surface since a decrease in the occupancy of the xy orbital results in the increase of the hole pocket size, as seen by ARPES studies in several AFe2As2 families3637. In particular, a comparison between ARPES measurements for BaFe2As2 and EuFe2As2 in the paramagnetic regime reveals that the xy hole-pocket is twice as large in EuFe2As2, suggesting a decrease of planar xy occupation.

In order to give this reasoning a more graphical view, we have explored how the exchange interaction J between the 3d conduction electrons and the local spin ESR probes may be affected by the orbital differentiation of the 3d bands. As such, we assume a scenario where J is proportional to the overlap between the atomic orbitals. Although simplified, this procedure can qualitatively capture the main features of the orbital differentiation through minute changes in J. Therefore, when the paramagnetic probe is out of the Fe-As plane, one can calculate the squared overlap between Eu 4f and Fe 3d orbital wavefunctions, 〈Ψ4|Ψ3〉2, as a function of the distance between the Eu and Fe ions, namely z0. More specifically, the Fe 3d wave-function is written as Ψ3 = R3(x, y, z)[B(Φ(x, y, z) + Φ(x, y, z)) + B Φ], where R3(x, y, z) is the 3d radial wave- function and Φ (x, y, z) are the cubic harmonics for each orbital. Ultimately, as shown in Eq. 1, the crucial result within this scenario is that the Korringa relaxation rate b is directly proportional to 〈Ψ4|Ψ3〉2.

Let us first consider the overlap between an isotropic Eu 4f configuration and the Fe t2 orbitals: xy and xz/yz. The triplet orbitals were chosen here due to their predominant role on the FS12. The squared overlap between the Eu 4f and the Fe 3d wavefunctions as a function of z0, normalized by the atomic radius of Eu2+ (z0 = 3.5 atomic units), is presented in Fig. 5a. We note that, when the planar xy orbital component is the dominant one (solid red line), its squared overlap with the 4f wavefunction decreases slightly faster than that with only xz/yz orbital components (dotted blue lines). One would expect such behaviour simply because the overlap between an 4f isotropic wave function and a 2D wave function decreases faster as the probe moves away from the Fe-As plane when compared to a more 3D wave function with z components.

Figure 5

Squared overlap between the Eu 4f isotropic wavefunction and the t2 Fe 3d wavefunctions xy (B = 1.0) and xz/yz (B = 0.0) when a) the probe is at the Fe-As plane; b) the probe is out of the Fe-As plane.

Moreover, we can now roughly simulate each Fe t2 orbital contribution to the FS using this phenomenological model. Starting with an isotropic Fe 3d configuration, i.e., a normalized population of for each considered 3d orbital, we are able to tune the xy orbital population - represented by the B coefficient - from up to 1. The inset of Fig. 5a shows at a particular Eu-Fe distance z0 = 4.0, chosen to better illustrate the behaviour mentioned above. It is clear that, by increasing the xy orbital contribution, the diminishes drastically, in agreement with the behaviour presented by the Korringa rate b in Fig. 4a for the transition metal M substitution into Ba0.8Eu0.2Fe2−MAs2.

We now turn our attention to the particular case that considered the local probe in the Fe-As plane (M = Cu). In this case, we are able to calculate as a function of the radius r from a Fe atom. Fig. 5b displays the calculated squared overlap between the Cu 3d and the Fe 3d wavefunctions with the same t2 orbital configurations of Fig. 5a. However, when the planar xy orbital component is the dominant contribution (solid red line) in this case, the overlap decreases slower than that for the xz/yz configuration. The inset of Fig. 5b shows an exponential increase of this squared overlap at r = 3.0 atomic units as a function of the xy weigh, indicating that higher planar contributions leads to higher Korringa rate b, in good agreement with our experimental data shown in Fig. 4b.

In summary, we have employed a sensitive microscopic spin probe to study the suppression of the SDW magnetic phase in BaFe2As2 via transition metal substitution. We demonstrated that the 3d electrons tend to localize at the Fe-As plane as the SDW phase is suppressed due to changes in the structural parameters, independent of the chosen transition metal substitution. When the Fe-As distances are lowered either by chemical substitution or hydrostatic pressure, the crystal field splittings change and, in turn, the planar (xy/x2−y2) orbital character of the Fe 3d bands is strongly enhanced at the Fermi level. These effects and consequent orbital differentiation corroborate to the suppression of the itinerant SDW magnetic order and the subsequent emergence of the magnetic-mediated superconductivity in these materials.

This work was supported by FAPESP, AFOSR MURI, CNPq, FINEP-Brazil and NSF-USA.

Methods

Single crystals of Ba1−EuFe2−MAs2 (M = Mn, Co, Cu, Ru, and Ni) were grown using In-flux as described in Ref. 35. The single crystals were checked by X-ray powder diffraction and submitted to elemental analysis using a commercial Energy Dispersive Spectroscopy (EDS) microprobe. No In-incorporation was detected in the crystals. From EDS analysis, we have extracted the actual x and y concentrations used throughout the text. The in-plane resistivity was measured using a four-probe method. Specific heat data were taken in a commercial small-mass calorimeter and the magnetization data was collected using a superconducting quantum interference device (SQUID) magnetometer. X-Band (ν = 9.34 GHz) ESR measurements were carried out in a commercial Bruker spectrometer with a continuous He gas-flow cryostat.

Author Contributions

P.F.S.R., C.A., T.M.G. and T. Grant grew the single crystals and performed pressure dependent transport measurements. P.F.S.R. performed ESR measurements. T.G. performed EDS measurements. P.F.S.R., C.A., T.M.G., T.G., Z.F., R.R.U., P.G.P. discussed the data and reviewed the manuscript.