Doping dependence and electron-boson coupling in the ultrafast relaxation of hot electron populations in Ba(Fe_{1-x}Co_x)_2As_2

Using femtosecond time- and angle-resolved photoemission spectroscopy we investigate the effect of electron doping on the electron dynamics in Ba(Fe_{1-x}Co_x)_2As_2 in a range of 0<x<0.15 at temperatures slightly above the N\'eel temperature. By analyzing the time-dependent photoemission intensity of the pump laser excited population as a function of energy, we found that the relaxation times at 0<E-E_F<0.2 eV are doping dependent and about 100 fs shorter at optimal doping than for overdoped and parent compounds. Analysis of the relaxation rates also reveals the presence of a pump fluence dependent step in the relaxation time at E-E_F = 200meV which we explain by coupling of the excited electronic system to a boson of this energy. We compare our results with static ARPES and transport measurements and find disagreement and agreement concerning the doping-dependence, respectively. We discuss the effect of the electron-boson coupling on the energy-dependent relaxation and assign the origin of the boson to a magnetic excitation.


I. INTRODUCTION
The discovery of Fe-based high-T c superconductivity is one of the breakthroughs in the last two decades of condensed matter research [1] because in addition to the cuprates scientists identified a second class of compounds displaying superconductivity with high T c enhancing the interest and hope in explaining this phenomenon. However, the richness of the physics contained in these materials does not make things smoother and continuously triggers new challenges in the research [2]. Understanding the microscopic origins of electronic phases in high-T c superconductors (HTSC) is important for elucidating the mechanism of superconductivity.
At high temperatures the parent compounds are bad paramagnetic metals with a Fermi surface consisting of three hole pockets at the Brillouin zone center (Γ) and two electron pockets at the zone corner (X/M ) [3,4]. The strong nesting between the pockets leads, at temperatures below the Néel temperature T N , to a transition into an antiferromagnetic (AFM) state, where the Fe spins assume a spin density wave ordering accompanied by a simultaneous or slightly precedent structural tetragonal-to-orthorhombic transition [5,6]. A consequence of the magnetic ordering is the backfolding of the electron bands to Γ where they hybridize strongly with the non-folded states, leading to the opening of energy gaps [7]. Before the AFM order is completly established, nematic order sets in, characterized by orbital anisotropy [8] and quantum fluctuations that persist above the structural transition [9,10]. Superconductivity arises by tuning a control parameter, as doping or pressure, that gradually reduces the AFM order until it is suppressed. A mainstream view is that the strong antiferromagnetic fluctuations at this point, called a quantum critical point (QCP), could be a good candidate for the pairing mechanism that leads to superconductivity and that these fluctuations would also account for the normal state non-Fermi-liquid (NFL) behavior. A QCP has been theoretically predicted [11,12], measurements of the London penetration depth [13] support a quantum critical point of view and experimental efforts have observed anomalous normal state electronic properties in the resistivity [14,15], optical conductivity [16], thermal [17] and NMR [18] studies that would explain the NFL-behavior of the metallic state and therefore support the quantum critical scenario. On the other hand, there are other models explaining the normal state properties of HTSC [19] and our own recent ARPES study [20] showed no evidence of quantum critical features but rather of a BCS-Bose-Einstein crossover state in the superconductive phase.
In this frame, further study is required and a valuable contribution could come from other experimental techniques like, for instance, doping-dependent time-resolved (tr) studies, which are lacking in this field. Tr-optical reflectivity [21] and trARPES [22] had been important for the study of the phonons in the system under non-equilibrium conditions after fs laser excitation and the determination of the electron-phonon coupling strength [23][24][25], pointing out that a pairing mediated by phonons is unlikely and the attention had to be turned to other "non-conventional" mechanisms. In particular trARPES is a powerful and versatile tool combining the direct access to the electronic structure in momentum space with the non-equilibrium regime excited by femtosecond laser pulses and probing the ultrafast relaxation dynamics [26]. In the work presented here, we report on a systematic trARPES study on parent (T N = 140 K) and Co-substituted, i.e. electron doped BaFe 2 As 2 with optimal (x = 0.06, T c = 23 K) and overdoped (x = 0.14) concentration, to investigate the influence of doping on the dynamics of the system. By analyzing the hot electron population decay as a function of energy above E F we find a trend of the relaxation for different doping in agreement with transport measurements [27,28] and a blocking of the relaxation at energies below 200 meV. The latter indicates coupling to a boson which is most probably of magnetic origin. We discuss our results in the context of the latest developed scenarios.

II. EXPERIMENTAL DETAILS
Samples are grown by a self-flux method [29] and cleaved in ultrahigh vacuum (p < 10 −10 mbar) at a temperature equal or lower than 80 K. Pump-probe measurements were performed after having increased the temperature to 160 K. A commercial regenerative   Fig. 1 (a). The overall time resolution was < 90 fs and the overall energy resolution, determined by the TOF spectrometer and laser bandwidth, was about 50 meV. For details about the experimental setup see Refs. [30,31]. Photoelectron detection was carried out in normal emission in the vicinity of the Γ point at T = 160 K, so that all the compounds have been initially, before optical excitation, in the paramagnetic phase and the only tuned parameter was the electron doping concentration.

III. EXPERIMENTAL RESULTS
Exemplary trARPES spectra are shown in Fig. 1 (b) for the parent compound at an incident laser pump fluence F = 540 µJ/cm 2 . Upon photoexcitation we observe an increase of spectral weight above E F up to 1 eV generated by electrons excited from outside the hole-pockets and subsequent relaxation dominated by cooling of hot electrons to the lattice, mediated by electron-phonon coupling. Coherent oscillations superimposed to the incoherent exponential decay are caused by the dominant excitation of the A 1g phonon mode [24,32]. In this work we focus on the analysis of the relaxation mechanism of the excited electrons as function of doping concentration and start here by investigating pump fluences of 50 − 540 µJ/cm 2 . To analyze the relaxation of the excited electrons we integrated the photoelectron intensity in an energy range from E F to about 1.5 eV. We describe the obtained time-dependent data with a single exponential decay plus a damped sine function (to account for the coherent oscillations) convoluted with the temporal pump-probe envelope, see Fig.1 (c). The resulting exponential decay fits are in the sub-picosecond timescale, in agreement with [22,23] and are summarized in Fig. 1 (d) for all samples and for different incident pump fluences. At low fluences, up to ∼ 400 µJ/cm 2 , the parent compound relaxes clearly slower than the Co-substituted compounds, showing also a dependence on fluence where, except for the first data point, we see a decrease from 570 fs to 490 fs for more than a tripling of F . Such an F -dependence is not observed in this analysis in the case of the optimally doped (OP) and overdoped (OD) compounds, where the relaxation times remain constant but even up to 200 fs shorter at 140 µJ/cm 2 . At higher fluences different relaxation times for different doping merge to joint relaxation times a bit below 500 fs.
At this point we emphasize that the previous analysis accounts for the decay of the energy integrated hot electron population. Contribution of electron-electron (e − e) scattering at higher energies and the respective secondary electrons closer to E F contribute in a far from trivial manner in addition to electron-phonon scattering which mediates close to E F energy transfer out of the electron system into phonons [25,[33][34][35]. The relaxation times reported in Fig.1(d) are considered as an effective characteristic quantity that cannot be compared with energy and momentum-dependent single particle scattering rates determined from linewidth analysis in static ARPES. The reason is that ARPES line widths are determined by i) elastic and inelastic scattering processes of ii) individual quasiparticle excitations, while 5 trARPES analyzes the population decay to which secondary electrons contribute in a nontrivial manner here after energy integrating the transient population [36,37]. We will come back to this point below and remark that nevertheless these effective relaxation times (Fig   1d) already exhibit a clear doping dependence. Further insight can be gained, as we show in the following, from an energy resolved analysis, though the difference between static and time-resolved ARPES will persist. We analyze time-dependent photoemission intensities in different energy windows starting from 20 meV to about 1 eV above E F , to obtain energydependent relaxation times (τ ) using a single exponential decay fit function I = A exp(t/τ )+ We observe a doping dependence until E − E F = 250 meV. As for the energy integrated intensities, the parent compound relaxes on a longer timescale than the Co-substituted compounds but it becomes now clear that the x = 0.06 (optimally doped) relaxes faster than the x = 0.14 (overdoped) compound. If we compare, for example, τ at E −E F = 0.1 eV we get 280 fs for the x = 0.06, while we obtain for x = 0.14 and parent 370 fs and 470 fs, respectively. These differences in τ in the first 250 meV are clearly visible at 140 µJ/cm 2 and 240 µJ/cm 2 but they vanish at the highest F of 540 µJ/cm 2 .
Two effects are noteworthy: the dependence of the relaxation times on x does not monotonously follow the carrier concentration and the fastest relaxation is found at optimal doping. In the limit of Fermi liquid theory for a three dimensional metal, the scattering probability due to e − e scattering is proportional to n −5/6 , with n being the number of carriers [38]. Therefore, at T > T N , where the system is metallic, one would expect a monotonic dependence of the relaxation on x. Interestingly similar anomalies were also found in transport measurements in the resistivity [27]. The derived scattering times [28] do not scale monotonically with Co-substitution and are shorter at optimal doping than at overdoping.   Another interesting point to consider is that the difference in the energy-dependent relaxation times for different Co-substitution is only present until E − E F = 250 meV. At this energy a change in the curves of τ vs E − E F appears as a discontinuity, which is more pronounced at the lowest F as marked by the arrow in Fig. 2(d) and absent at the maximal F . A closer look at data in that energy region, is instructive. In Fig. 3 time dependent intensities of the x = 0.06 compound for energy windows centered around 150 meV (orange) and 250 meV (light blue) are compared for 140 mJ/cm 2 (Fig. 3(a)) and for 240 µJ/cm 2 ( Fig.   3(b)). At the lowest F the difference in relaxation time between the two subsequent energy windows is obvious, while it is weaker at the intermediate F . Such behavior becomes also evident in Fig. 4 where relaxation rates Γ = 1/τ as function of binding energy are shown.
At the lowest F ( Fig. 4(a)) a step-like increase in the energy-dependent relaxation rates from E − E F < 200 meV to the next values above this energy is observed for all x. Interestingly, the visibility of the step is reduced with increasing F (Fig. 4 (b)) until it completely vanishes (Fig. 4 (c)).

IV. DISCUSSION
For the investigated compounds we can reasonably assume a rigid band model [39], that is, the effect of electron doping is to add carriers to the system shifting the Fermi level in agreement with [37] where the same analysis was applied to a cuprate HTSC and is a confirmation of the fact that the lifetimes measured by static ARPES and the relaxation times obtained here close to E F represent two different quantities (see also III). In the case of trARPES, relaxation of the laser-excited non-equilibrium electron population is analyzed, while in static ARPES the scattering probability of a single particle is determined. In addition, one should recall that static ARPES probes occupied states, while trARPES in the present case probes unoccupied states. Since it is known that the self energy above and below E F is different for the case of hole-doped compounds [42], we can expect some differences in the respective relaxation times also in the present case.
On the other hand, in the doping range under investigation, qualitative agreement is found with transport measurements, where the same trend with doping found in this work is observed in resistivity measurements, and the determined e − e scattering times [27,28] are larger near optimal doping than at over doping.
A widely shared interpretation [13-16, 27, 43], sees antiferromagnetic fluctuations as being mainly responsible for these anomalies, because they occur near the QCP. The deviation of the T -dependent term of the resistivity from a quadratic to a linear scaling, accounting for the non-Fermi liquid behavior in the normal state, would qualitatively alter the Tdependence of the relaxation rate in the quantum critical regime, leading to the maximum scattering rate at optimal doping discussed above. Alternative explanation for the strange normal state properties of these systems near a QCP involves the influence of a Lifshitz transition, occurring at optimal doping [44,45], and correlation effects on the scattering rates of the charge carriers [20,41].
Given these considerations, the deviation from the behavior of a conventional metal which we observe in our trARPES study is not surprising, considering that FL theory does not hold. In discussing the relaxation process, the 3D character of the Fermi surface in these materials has to be taken into account. The electronic structure for these compounds is strongly doping dependent [46] and, in particular, the k z dispersion increases with increasing Co concentration. That means that the observation of a Lifshitz transition, where the top of the hole pocket just touches the Fermi level is detectable only at k z = 0, in the case of optimal doping. The phase space for relaxation is complex and a complete discussion is beyond the scope of this paper, nevertheless we have to consider that, due to our probe pulse energy, we work at a finite k z value corresponding to a region between Γ and Z where we detect some remnant of the hole pocket in the case of the optimally doped compound, providing relaxation phase-space for e − e scattering. Moreover, although we work at T > T N , interband scattering (coupling between electron and hole pockets) has to be considered as a possible additional contribution to the relaxation, which is instead absent in the case of the overdoped compound, in agreement with our observation of a faster relaxation at optimal doping. Although from our results we cannot quantify the intra-and interband conditions discussed above, we remark that despite ARPES and trARPES are closely related, qualitative agreement is found with results obtained by transport, another dynamical method.
Our second important observation comes from the analysis of the energy dependent relaxation rates at binding energies higher than those accessed by static ARPES. A step in the relaxation rates is recognized at low F ( Fig. 4 (a)) around E − E F = 200 meV for all x. This is an indication that at these energies an additional relaxation channel contribute, causing the system to suddenly relax faster. Band structure calculations predict the pres-ence of an unoccupied 4p As band, at Γ, to be 0.1 to 0.3 eV above the Fermi level [3,47].
This could provide phase space for additional inelastic interband scattering due to a peak in the density of states around that energy. However, if this was the origin of the step visible in our data, we would expect the central energy at which the step occurs to change with doping, due to the shift of the chemical potential. In addition to that, no signature of such a feature is visible in our trARPES intensities spectra (analogous to Fig. 1 (b)) for any of the dopings and fluences. A second possibility, which we assume as the most plausible, is that an additional relaxation channel is provided by a bosonic excitation of a certain energy Ω that would couple to the excited electronic system. Theoretical modeling of the effects of electron-boson coupling (e − b) in trARPES [48] has shown to lead to a characteristic change in the decay rates, in the form of a step. The fluence dependence of the step is a further confirmation of the presence of e − b coupling, as we explain below. Electrons excited at E − E F > Ω and the respective holes injected at E − E F < − Ω would relax faster than those excited within the energy region − Ω < E − E F < Ω, as shown originally from [49] in the case of a coupled electron-phonon system, and subsequently by Kemper et al. [50] for electron-boson coupling in general. According to this description, at equilibrium, scattering from inside the bosonic energy region to states below the Fermi level is suppressed, as the electron occupation does not provide the necessary phase space for the excited electrons to relax back to energies below E F . Therefore the relaxation times inside the boson window are larger compared to those at energies outside, at E − E F > Ω. In the time-domain analysis the infrared pump redistributes the electronic population, which modifies the phase space. In particular, relaxation channels which were blocked in equilibrium become now available, leading to shorter relaxation times under non-equilibrium conditions. At small F one expects a situation similar to equilibrium where the signatures of the boson window in the relaxation times and the electron population are only weakly perturbed [48,50]. Indeed our data in Fig. 4(a) show that at the smallest excitation density the step is larger showing higher (lower) scattering probability above (below) a boson energy of Ω = 200 meV. With increasing F (Fig. 4 (b)) this step reduces as the pump-induced phase space redistribution provides more relaxation channels until the signature of the boson window completely disappears at the highest F (Fig. 4 (c)). Such an effect was also experimentally observed in the case of cuprates [51].
To assign the origin of the bosonic excitation, we exclude phonons as a possible coupling candidate because the energy of 200 meV at which the effect is found exceeds the high energy cutoff of the phononic spectrum, around 35 meV [52,53]. A possible candidate is then a boson of magnetic origin. Spin fluctuations in this high energy region were reported in literature for electron doped compounds [54,55] and for the parent compound [10], where they are shown to persist also at T > T N . Electron-magnon coupling was also observed in ARPES measurements on ferromagnetic Fe [56], showing a saturation of the increase of the imaginary part of the self-energy at 160 meV, consistent with the energy of spin waves in Fe.

V. CONCLUSION
We have performed a femtosecond trARPES study on Co-substituted BaFe 2 As 2 with different doping concentration. From the analysis of the energy dependent relaxation times we found that the optimally compound relaxes faster than the overdoped and parent compound.
Our findings are therefore in qualitative agreement with the doping-dependent results of transport measurements and interband coupling between electron and hole pockets has to be considered as a possible contribution in addition to intraband relaxation, due to the finite relaxation phase-space at optimal doping, because we probe at a finite k z . However, it should be noted that no k z dependence of the scattering rates was observed in static ARPES. From the analysis on the energy dependent relaxation rates we could clearly identify the presence of a step around the energy E − E F = 200 meV, which is best observed for small laser pump fluences. We conclude that this indicates the coupling of the excited electronic system to a bosonic excitation of magnetic origin.
We have demonstrated that trARPES studies are a valid contribution in elucidating the underlying mechanisms present in these complex materials, complementary to other experimental technique discussed here, although further study will be required.