Electron-phonon coupling in 122 Fe pnictides analyzed by femtosecond time-resolved photoemission

Based on results from femtosecond time-resolved photoemission, we compare three different methods for determination of the electron-phonon coupling constant {\lambda} in Eu and Ba-based 122 FeAs compounds. We find good agreement between all three methods, which reveal a small {\lambda}<0.2. This makes simple electron-phonon mediated superconductivity unlikely in these compounds.


Introduction
Albeit an enormous amount of research has been conducted in the past few years on the Fe pnictide high-temperature superconductors (HTSCs) [1,2], the search for the superconducting pairing mechanism in these materials is still ongoing. Besides other excitations such as spin fluctuations [3], the electron-phonon (e-ph) coupling responsible for Cooper pairing in conventional superconductors in the Bardeen-Cooper-Schrieffer (BCS) theory is considered a potential candidate. Therefore, the quantitative determination of the e-ph coupling constant λ in the Fe pnictides is of particular interest. Several methods to determine λ have been established in the literature and successfully used in various studies. In the energy domain, angle-resolved photoemission spectroscopy (ARPES) analyzes the renormalization of the single-particle spectral function by the electronic self energy Σ, which includes the e-ph coupling. However, such type of analysis requires very high sample and data quality, which is difficult to obtain in the Fe pnictides, and only few ARPES studies of FeAs compounds were able to determine effects of Σ [4,5,6]. In addition, in thermal equilibrium, contributions from other degrees of freedom to Σ are often difficult to disentangle from e-ph coupling.
In contrast, femtosecond (fs) time-resolved spectroscopies allow to separate e-ph scattering from other relaxation channels like e.g. electron-electron (e-e) scattering or heat diffusion due to their different intrinsic timescales out of thermal equilibrium [7,8,9]. Here, time-resolved ARPES (trARPES) is used, which combines the energy and momentum resolution of ARPES with fs time resolution into a powerful tool to directly investigate the dynamics of the electronic structure in a non-equilibrium state [10,11].
The e-ph relaxation in metals after photoexcitation has been successfully described using the two-temperature model (2TM) [7,8,9,12,13,14]. This model describes a system of two coupled heat baths for the conduction electrons and the ion lattice with temperatures T e and T l by coupled rate equations for T e and T l , respectively. One key assumption of this model is that e-e and phonon-phonon (ph-ph) scattering occurs on a much faster timescale than the e-ph scattering, leading to a thermal distribution within each subsystem. However, in correlated electron systems like the cuprate or Fe pnictide HTSCs, e-ph scattering might occur on similar timescales [15], which poses some questions on the applicability of the 2TM for these materials and care has to be taken. Therefore, we analyze our data in terms of a suitable version of the 2TM and compare the results with two complementary methods for the determination of the e-ph coupling strength. We find that all three methods to analyze the e-ph coupling strength from time-resolved photoemission data yield results for the second moment of the Eliashberg coupling function λ ω 2 which agree well within error bars and reproduce the trends for three different 122 FeAs compounds consistently.

Experiments and results
For trARPES experiments, single crystals of EuFe 2 As 2 and BaFe 2 As 2 parent compounds and optimally doped BaFe 1.85 Co 0.15 As 2 (T c = 23K) were cleaved in ultrahigh vacuum (p < 10 −10 mbar) at T = 100 K where most measurements were carried out. The experimental setup is sketched in fig. 1(a): The output of a commercial regenerative Ti:Sapphire amplifier (Coherent RegA 9050) delivering ultrashort laser pulses at hν 1 = 1.5 eV photon energy with a pulse duration of 55 fs, and operating at 300 kHz is used to optically excite the samples (pump pulse). A time-delayed frequency-quadrupled probe pulse at hν 2 = 6.0 eV photon energy with a pulse duration of 80 fs leads to the emission of photoelectrons, which are detected using an electron time-of-flight (TOF) spectrometer with an acceptance angle of ±3 • . The energy resolution of 50 meV is mainly determined by the spectral width of the probe pulses, and the overall temporal resolution is < 100 fs. For details see [9].

Three-Temperature model
Exemplary trARPES data of EuFe 2 As 2 are shown in fig. 1(b) for various pump-probe delays on a logarithmic intensity scale. After excitation, a pronounced distribution of excited charge carriers is formed at E − E F > 0.1 eV, which deviates from the Fermi-Dirac distribution of the electronic system before excitation and which originates from hot, non-thermalized electrons. Subsequently, e-e scattering quickly leads to thermalization of these non-thermal electrons and to the formation of a hot thermalized electron distribution [9]. At t > 250 fs, the non-thermal contribution has decayed to less than 1% of the electron population. The temperature of the thermalized part of the electronic system T e can be extracted from the trARPES data by fitting a Fermi-Dirac distribution function to the high-energy cutoff of the transient spectra, multiplied by a phenomenological densityof-states function and convoluted with an instrumental resolution function [9,14], as shown in fig. 1(b) (solid lines). T e determined by the fitting is shown in fig. 1(c) as a function of pump-probe delay. After the steep rise of the electronic temperature at zero pump-probe delay, we find a relaxation on two distinct timescales τ α and τ β . Such a behavior indicates the selective coupling of hot electrons to a subset of strongly coupled phonon modes on the faster timescale τ α and the subsequent energy transfer to the rest of the phonons by anharmonic ph-ph scattering, determined by τ β . Thus, the system is described by an extended version of the 2TM, which had been developed for the cuprate HTSCs [14], where a similar relaxation dynamics was found, and which has been also used recently for the Fe pnictides [16,17]. This three-temperature model (3TM) consists of a system of three coupled differential equations describing the temperature of the electrons, T e , of a hot phonon distribution T p and of the rest of the lattice modes, T l : The source term S describes the optical excitation, C e , C p and C l are the specific heat capacities of electrons, strongly and weakly coupled phonon modes, respectively. The anharmonic decay of phonons is described by τ β . For the energy transfer from the electrons to the more strongly coupled phonons, the formula derived by Allen [7] is used: where k B is the Boltzmann constant. This relation allows determination of the second moment of the Eliashberg e-ph coupling function λ ω 2 . A fit of a numerical solution of the 3TM to T e is shown in fig. 1(c) and yields a good agreement to the experimental data. Due to an ambiguity of the model parameters such as the electron and lattice specific heat capacities reported in the literature, a Pump-probe delay (fs)  range of values for λ ω 2 is retrieved, where the fits show similar good agreement to the data. Details of this analysis can be found in [17]; here, only the main results are discussed. From the 3TM, we find values of λ ω 2 = 56 − 65 meV 2 for EuFe 2 As 2 , while in BaFe 1.85 Co 0.15 As 2 it is slightly smaller (λ ω 2 = 46 − 55 meV 2 ) and even smaller in undoped BaFe 2 As 2 (λ ω 2 = 30 − 46 meV 2 ).

Electronic excess energy
The analysis within the 3TM was based on the energy relaxation within the thermalized part of the transient electronic distribution function and thus neglected the non-thermal electrons present at early delay times (see fig. 1(b)). However, the analysis of the energy relaxation of these non-thermal, excited electrons can also provide information about the strength of e-ph coupling. The scattering rate of electrons excited at energy = E−E F above the Fermi level is in the self-energy formulation of many-body theory determined by Γ = τ −1 = 2ImΣ( ), where ImΣ is the imaginary part of the electronic self energy [18]. Important energydependent contributions to Γ arise from e-e and e-ph scattering, where e-e scattering is considered to follow the quadratic energy scaling of the Fermi liquid theory, Γ e−e = 2β [(πk B T ) 2 + 2 ] ( fig. 2(a)). As the proportionality coefficient β is rather small in typical metals (in the order of 10 −2 eV −1 to 10 −1 eV −1 ) [19,20], this contribution is negligible compared to the e-ph scattering at low enough excitation energies. The latter increases up to the maximal phonon energy ω max and is constant above ( fig. 2(a)). For > ω max and T = 0 K, Γ e−ph results for an Einstein mode ω 0 to [21] Γ e−ph = π λω 0 .
Within the energy window between ω max and the crossover regime, where Γ e−e becomes dominant (shaded area in fig. 2(a)), the rate of energy dissipation of an electron due to the emission of a phonon with energy ω 0 is given by which leads to a linear relaxation of the electron energy [22]. The rate of energy relaxation of excited electrons can be extracted from the experimental trARPES intensity I( , t) by analyzing the mean excess energy, However, this integral represents the mean excess energy within a selected partition of a transient distribution function and not the energy relaxation of single individual electrons. Thus, special care has to be taken with the determination of the integration boundaries, 0 and 1 . A careful investigation of the influence of the integration boundaries and the excitation fluence reveals that a reasonable choice of 0 = 70 meV and 1 = 300 meV allows to give a good estimate of the e-ph coupling strength, as detailed in [23]. Experimental data of EuFe 2 As 2 , BaFe 1.85 Co 0.15 As 2 and BaFe 2 As 2 near k F and at T = 100 K are shown in fig. 2(b). To minimize lattice heating and the influence of the hot thermalized electron distribution, low excitation fluences of F ∼ 50 µJcm −2 have been used. The linear fits to equation (6) within the first 100 fs yield values of λ ω 2 = 65(5) meV 2 for EuFe 2 As 2 , 50(3) meV 2 for BaFe 1.85 Co 0.15 As 2 and 34(6) meV 2 for BaFe 2 As 2 . Here, error bars represent the numerical uncertainties of the fits, while the overall errors have to be considered larger due to the uncertainty of the integration boundaries discussed above and of the elevated lattice temperature. Nevertheless, we find a good agreement with the values obtained in the 3TM simulations and a considerably higher e-ph coupling in EuFe 2 As 2 than in BaFe 1.85 Co 0.15 As 2 and even lower λ ω 2 in BaFe 2 As 2 .

Temperature-dependent hole relaxation rates
Finally, another estimate of the e-ph coupling strength can be gained from the temperature dependence of quasiparticle (QP) relaxation rates [7,15,24,25]. A recent theoretical investigation of the Boltzmann equation for e-e and e-ph interaction found an analytic solution for the temperature dependence of the QP relaxation rate τ [24]. In the limit of a bad metal at elevated temperatures, τ depends linearly on the lattice temperature T [15,24,25]: This relation allows for an independent determination of the e-ph coupling strength λ ω 2 . Experimentally, the temperature-dependent relaxation of holes in the hole pocket near Γ in EuFe 2 As 2 are used, which are found to be independent from the antiferromagnetic transition [26]. The respective hole relaxation time τ holes is determined from the temperature-dependent trARPES intensity by fitting exponential decay curves ( fig. 3(a)) and is shown in fig. 3(b) as a function of temperature. The fit to equation (8) for T > 100 K reveals a good agreement to the linear behavior and yields a value of λ ω 2 = 90(30) meV 2 , where the error represents a confidence interval of 80%. This value, albeit considerably higher than the values found in the 3TM and in the evaluation of the excess energy, is in agreement with the other evaluations within error bars.

Discussion
The values of λ ω 2 determined by the three methods discussed above are compared in table 1. Despite the limitations and drawbacks of the various methods, we find a perfect  (6) agreement of all three methods within error bars. Especially, the trend for larger e-ph coupling in EuFe 2 As 2 than in BaFe 1.85 Co 0.15 As 2 and in BaFe 2 As 2 is nicely represented.
Based on our results of λ ω 2 we can estimate the value of the e-ph coupling constant λ for a particular value of ω. Considering the Raman active A 1g mode at 23 meV, which can be coherently excited [17,28,29] and therefore shows enhanced e-ph coupling in the system, we find λ < 0.2 for all compounds. This estimate is in agreement with calculations of various Fe-pnictide compounds, which report values of λ < 0.35 [30,31]. Taking the mean of the phonon spectrum as reference, λ gets even smaller, in agreement with other publications [16,25]. Even if we consider the lowest coupled modes around 12 meV to be most important for e-ph coupling, λ does not exceed a value of 0.5. Similarly small values for λ have been found in the cuprate HTSCs [14,15], which suggests limited importance of e-ph coupling for the pairing mechanism in both classes of materials.
These boundaries for the e-ph coupling constant λ allow to estimate the superconducting critical temperature T c , assuming a conventional BCS-type pairing based on e-ph coupling. Here, T c in isotropic systems and in a strong coupling regime is given by McMillan's formula [32], modified by Allen and Dynes [33], where µ * is the effective Coulomb repulsion and ω log is the logarithmic average of the phonon spectrum. Taking µ * = 0 and ω log = 205 K [30], which provides an upper limit for T c , we find T c = 0.33K and T c = 7.5K for λ = 0.2 and λ = 0.5, respectively. Comparing these results with the experimentally found T c = 24 K in optimally doped BaFe 1.85 Co 0.15 As 2 [34] and even up to T c = 38 K in K doped BaFe 2 As 2 [35] demonstrates that e-ph coupling in a conventional BCS pairing scenario cannot explain the high critical temperatures found in the pnictide HTSCs. In addition, albeit the stronger e-ph coupling, the critical temperatures found in EuFe 2 As 2 upon Co [36] and K [37] substitution are considerably smaller compared to BaFe 2 As 2 .
However, e-ph coupling might still play a significant role in the Cooper pair formation in the Fe pnictides. For instance, a very strong sensitivity of the Fe magnetic moment on the As height in the FeAs tetrahedra was found in DFT band structure calculations [38], with a rate of 6.8 µ B /Å. Inelastic x-ray scattering experiments on CaFe 2 As 2 concluded on strong coupling of phonons to magnetic excitations even in the high-temperature paramagnetic phase [39]. Furthermore, a strong sensitivity of the maximum critical temperature on the pnictogen height was found [1]. On the basis of these findings a strong magneto-phonon coupling in these compounds was proposed and, e.g. the strongly coupled A 1g mode perpendicular to the Fe layers, which modulates the pnictogen height, could mediate superconductivity in the Fe-pnictides [40]. Such a scenario is supported by the strong coupling of this particular mode to electronic states directly at the Fermi level [17] evidenced by the coherent excitation of this phonon mode [17,28,29].

Summary
In summary, we demonstrated the quantitative analysis of the e-ph coupling strength from femtosecond time-and angle-resolved photoemission spectroscopy. In detail, we compared three different methods to determine the e-ph coupling strength in three 122 FeAs compounds from time-and angle-resolved photoemission experiments. The transient temperature of the thermalized electronic distribution is analyzed in EuFe 2 As 2 , BaFe 1.85 Co 0.15 As 2 and BaFe 2 As 2 by a three-temperature model, while the rate of energy relaxation of the non-thermal electrons yields direct information on the second moment of the Eliashberg coupling function λ ω 2 . Finally, the temperature-dependence of hole relaxation rates in EuFe 2 As 2 also allows to determine λ ω 2 . All three methods consistently yield a small λ ω 2 < 100 meV 2 , which results in λ < 0.2 using a reasonable choice of the phonon spectrum. This value is discussed to be too small to explain superconductivity in the Fe pnictides by a conventional BCS-type pairing.

Acknowledgements
Experiments were carried out at the Freie Universität Berlin. We acknowledge discussions with Martin Wolf and Hermann Dürr in the early stage of the project and experimental support by Setti Thirupathaiah. This work was supported by the Mercator Research Center Ruhr and the Deutsche Forschungsgemeinschaft through BO 1823/2 and SPP 1458. R.C. acknowledges the Alexander von Humboldt Foundation.