Crossover from spin-waves to diffusive spin excitations in underdoped Ba(Fe{}_{1-x}Co{}_{x}){}_{2}As{}_{2}

Crossover from spin-waves to diffusive spin excitations in underdoped Ba(Fe$_{1-x}$Co$_{x}$)$_{2}$As$_{2}$

Abstract

Using inelastic neutron scattering, we show that the onset of superconductivity in underdoped Ba(FeCo)As coincides with a crossover from well-defined spin waves to overdamped and diffusive spin excitations. This crossover occurs despite the presence of long-range stripe antiferromagnetic order for samples in a compositional range from , and is a consequence of the shrinking spin-density wave gap and a corresponding increase in the particle-hole (Landau) damping. The latter effect is captured by a simple itinerant model relating Co doping to changes in the hot spots of the Fermi surface. We argue that the overdamped spin fluctuations provide a pairing mechanism for superconductivity in these materials.

pacs:
74.70.Xa, 75.25.-j, 61.05.fg

The key to many unconventional superconductors lies in their proximity to an ordered antiferromagnetic (AFM) phase (3); (2); (1). As a ground state that competes with superconductivity (SC), the suppression of AFM order (by chemical tuning or applied pressure) is required for the SC state to appear. However, the vestiges of AFM order that remain as correlated spin fluctuations have been proposed to provide the glue that pairs electrons in the SC state (4). It is, therefore, very important to understand how the spin excitations evolve from collective spin waves in the AFM ordered state to the overdamped correlated spin fluctuations characteristic of the SC state. In the iron pnictide Ba(FeCo)As, the suppression of AFM ordering upon Co substitution of a few percent allows a SC ground state to appear (5) in the presence of substantial spin fluctuations at the AFM wavevector, . Unlike some unconventional superconductors, e.g., BaKFeAs, the competing AFM ordered and SC states actually coexist microscopically in a limited compositional range from , the so-called underdoped compositions (6). This allows one to investigate how the normal state spin fluctuations provide the conditions for SC to emerge, even in the presence of weak AFM order.

Given the important connection between superconductivity and magnetism, extensive studies of the magnetic dynamics have been performed in these compounds as a function of composition. The magnetic dynamics of electron doped compounds, Ba(Fe)As ( Co, Ni) have been studied in some detail by inelastic neutron scattering (INS) (7); (8); (9); (10); (11); (12). These investigations found that the high-energy spin dynamics ( meV) are relatively insensitive to electron doping whereas the low-energy spin dynamics show a strong dependence on electron doping. Deep within the AFM ordered state of the parent BaFeAs compound (Néel transition temperature, K), the low-energy spin dynamics are dominated by a large spin gap meV that characterizes the ordered AFM state (13). Above the spin gap, very steep spin wave excitations propagate within the Fe layer, while much lower spin wave velocities are found for modes propagating perpendicular to the layers, indicative of quasi-two-dimensional magnetism. In BaFeAs (14), as well as CaFeAs (15) and SrFeAs (16), the low-energy spin waves have very long lifetimes (no substantial energy-dependent damping). The large spin gap and small damping of the collective spin wave modes highlight the robust AFM state of the parent compounds. Within an itinerant spin-density wave picture for the AFM order in the iron pnictides, such behavior indicates that the electronic spin-density wave (SDW) gap is large, estimated to be meV via optical conductivity measurements (17), thereby gapping out particle-hole (Landau) damping mechanisms. Note that while the spin gap is related to anisotropies in spin space (such as single-ion anisotropy), the SDW gap is proportional to the magnetization and, therefore, to the energy gain in the magnetically ordered state.

At the opposite extreme, those compositions without long-range AFM order ( for Co substitutions) display low-energy spin excitations that are diffusive (overdamped) in nature, and typical of systems close to a critical point (18); (8). The low-energy spin fluctuations are still centered at , but appear gapless and are characterized by a finite spin-spin correlation length () and relaxational energy scale () related to the Landau damping. The presence of substantial magnetic spectral weight at low energies (as obtained from a gapless spectrum with , with superconducting transition temperature ) is considered an important ingredient for magnetically-mediated SC (2).

Here, we study the evolution of the normal state spin dynamics between these two extremes. Most of the compositions are underdoped, possessing both weak AFM ordering and superconductivity at low temperatures (i.e. small SDW and SC gaps). In the normal state of the underdoped compounds, we find clear signatures of diffusive behavior in the low-energy spin dynamics (spatial disorder and a gapless spectrum with overdamped dynamics) despite the AFM ordering. The crossover of the spin dynamics is associated both with the collapse of the spin-density wave gap, see Ref. (19), and the subsequent development of strong Landau (particle-hole) damping. This crossover coincides with the appearance of SC in underdoped samples.

The INS measurements were carried out on the HB3 spectrometer at the High Flux Isotope Reactor at the Oak Ridge National Laboratory. Samples were grown and characterized as outlined in Ref. (5) and were mounted in the scattering plane. We define = = in reciprocal lattice units as referenced to the tetragonal unit cell. Details about the instrumental configuration and resolution are given elsewhere.

Figure 1: Temperature dependence of INS data for Ba(FeCo)As with [(a, c, e)] and [(b, d, f)] plus fits to the spinwave model. (a, b) Energy scans at performed at the indicated temperatures are offset vertically. (c, d) Reduced temperature dependence of spinwave model parameters (open symbols) and (filled symbols). (e, f) Reduced temperature dependence of the ordered magnetic moment, , normalized to its low-temperature value. Light gray symbols in panels (a) and (b) represent measured intensity which was excluded from fitting due to concerns with the validity of background estimates at those points.

Figure 1 shows the spectrum of Ba(FeCo)As at several different temperatures for lightly doped and nonsuperconducting and . The spectra are dominated by a large spin gap at 10 meV for both compositions. These data can be fit to a damped spin wave form for ,

(1)

where is a spin gap, is the spin wave velocity, is a damping rate, and is the reduced momentum transfer. The full anisotropic form for the damped spin wave susceptibility is given elsewhere. For at K, meV is small in comparison to other energy scales and, in principle, can arise from a combination of different damping processes (such as Landau damping for energy scales larger than the SDW gap, or magnon-magnon interactions). The fit to the K data shows a large spin gap meV characteristic of the parent AFM ordered state.

The solid lines in Fig. 1 (a) and (b) represent independent fits to the damped spin wave model where the gap and damping rate are allowed to vary freely. The magnitude of the spin gap is determined to be nearly constant with temperature up to our closest approach of where the ordered magnetic moment . Similar to the results described for NaFeAs (25), BaFeAs (25), and LaFeAsO (26), we find that the spin gap energy, , is roughly 9 meV in the ordered state, regardless of size of the ordered moment or the concentration, , and the dynamics become overdamped as is approached.

Figure 2: Background subtracted INS intensity of Ba(FeCo)As corrected for the Bose thermal population factor and the Fe single-ion magnetic form factor plus best fit lines to the diffusive (light green lines) and the damped spin-wave (black lines) models. (a-e) Constant- energy scans at for five compositions. (f-h) Constant- scans in the -direction across at meV. (i-j) Constant- scans in the -direction across at meV. (k-m) Constant- scans in the -direction, perpendicular to the Fe layer, across at meV. (n-o) Constant- scans in the -direction across at meV. Light gray symbols represent measured intensity which was excluded from fitting due to concerns with the validity of background estimates at those points.

Figure 2 shows a series of representative low-energy INS scans taken in the AFM ordered and normal state () for each composition. Upon increased Co substitution, the spin gap appears to gradually close [Figs. 2(a)-(e)] and is completely absent at . One can also observe a gradual reciprocal space broadening of the longitudinal cut [Figs. 2(f)-(j)] with increasing Co composition. Finally, the modulations along [Figs. 2(k)-(o)] are reduced, signaling a gradual evolution to two-dimensional spin dynamics. Within the damped spin wave model of Eq. (1), the data at all compositions have been successfully fit by assuming that: in accordance with our temperature-dependent results, the spin gap remains constant; the damping increases dramatically with ; and both the in-plane and inter-plane spin wave velocities are reduced with . However, it is clear from high-energy INS investigations that the in-plane spin velocities are independent of composition (see the discussion in Ref. (11)) and constraining the in-plane velocity to this value leads to poorer and poorer agreement of the low-energy data with the damped spin wave model (as shown by the black lines in Fig. 2).

One major assumption of our data analysis using the spin wave model is that the spin gap is independent of composition. If the spin gap is due to single-ion anisotropy, then its magnitude should be proportional to some power of (24). Data fitting in which the spin gap was allowed to freely vary resulted in an increase of the gap with composition, and fits in which the spin gap was constrained to be proportional to gave worse results.

The increased reciprocal space broadening suggests that another length scale must be introduced for low-energy magnetic fluctuations, such as a spin-spin correlation length. Considering also the gapless form of the magnetic excitations, the data at higher compositions resemble the diffusive response that has been used to describe the optimal and overdoped samples.(8); (20) This diffusive response has the form

(2)
(3)

where is the spin-spin correlation length, is the tetragonal lattice constant, and is the Landau damping coefficient. One can also define as the spin relaxation rate. Fits to the diffusive form for are shown as light green lines in Fig. 2. While the diffusive form does a poor job at the lowest compositions where the spin gap is sharp, it works exceptionally well at the higher compositions where the spectrum appears gapless and the increased reciprocal space broadening for longitudinal scans shown in Fig. 2(f)-(j) is captured by a smaller correlation length.

Figure 3: (a) Phase diagram of Ba(FeCo)As showing regions of AFM, SC, and their coexistence; colored symbols show the locations in phase-space of the measurements performed in this study. (b-f) Select model parameters as a function of composition for the diffusive (green diamonds) and spin wave (black circles) models. All data points shown in (b-h) were determined at the lowest temperature indicated in (a), and the lightly shaded background indicates compositions which exhibit SC at low temperature. Spin-wave model: (c) spin gap (filled), damping (open), and (diamonds); (e) Inter-plane spin wave velocity. Diffusive model: (b) Landau damping and the corresponding theoretical prediction; (d) Spin relaxation characteristic energy (filled) and the effective magnetic energy (open); (f) correlation length . In (c-d) an estimate for the SDW gap – derived from from Ref. (6) and from Ref. (17) – is also shown (tan line). (g) Residual for each fit model. (h) Spectral weight of the excitation. The spectral weight is the -averaged energy integration of the trace of the imaginary component of the magnetic susceptibility tensor. The averaging-range in here is , , ; the energy integration is over the range meV. All errorbars represent the combined errors for all function parameters. The solid green and black lines in (c-h) are guides to the eye.

Figure 3 shows the locations of our measurements in a phase-space diagram, the fitting parameters for both the damped spin wave and diffusive models in Eqs. (1) and (3), a measure of the goodness-of-fit for the constant- and constant- scans presented in Fig. 2 for these two models, and the composition-dependence of the low-energy spectral weight. For , the damped spin wave model is the best and is consistent with underdamped dynamics. For intermediate composition, , both models are of comparable quality. As shown in Fig. 3(c), within the damped spin wave model and the dynamics have become overdamped causing the spin gap to disappear. In the limit where , the overdamped spin wave model also takes on a relaxational form with ; as shown in Figs. 3(c) and (d) , , and the effective magnetic energy, , decrease as the critical concentration for which the AFM order is fully suppressed is approached, as indicated by vanishing . As seen in Figs. 2(k-o) the excitation becomes increasingly two-dimensional with as captured by the damped spin-wave model parameter , Fig. 3(e). For , , and , the diffusive model becomes the better fit, as the smaller correlation length [Fig. 3(f)] is able to capture the reciprocal space broadening of the in-plane spin fluctuations. A comparison of the residual for each fit model in Fig. 3(g) clearly shows the crossover from spin-wave- to diffusive-like excitations. In Fig. 3(h), a sharp increase in the low-energy spectral weight ( meV) coincides with the appearance of SC.

From Fig. 3, regardless of the model used to fit our data, it is clear that upon approaching the optimally-doped composition, damping becomes stronger, the spin fluctuations acquire a more two-dimensional character, and the energy scale associated with these fluctuations ( or ) become smaller. These features, as well as the crossover from spin-wave to diffusive excitations, are consistent with a suppression of the SDW gap upon doping. In Figs. 3(c-d), we show the experimentally determined suppression of obtained by combining the doping evolution of the zero-temperature ordered magnetic moment, , from Ref. (6) with the optical conductivity derived value of from Ref. (17), using the fact that (6); (22); (21); (23).

Based on this information, we can conclude that the presence of sub-gap spectral weight which appears with either an increase in temperature or Co composition is driven entirely by damping. For the temperature-driven transition, we find an increase of damping close to . Given the similarities between the spin fluctuations above and below near optimal doping and the smallness of the spin-wave gap in this regime [see Fig. 3(c)], we compare the fitted damping rate with the calculated Landau damping due to the decay of spin excitations into particle-hole pairs near the Fermi level in Fig. 3(b). Using a simplified two-band model, which was previously shown to successfully capture the coexistence of SC and AFM (22); (21); (23), the Landau damping is given by (27), where and are respectively the Fermi velocities of the electron and hole pockets at the hot spots (i.e. points connected by the AFM ordering vector ). Upon Co substitution, electrons are introduced into the system, making the hole pocket shrink and the electron pocket expand. As revealed by ARPES (28), this moves the hot spots, making their Fermi velocities become nearly parallel around optimal doping. As a result, , as seen experimentally. Note that describes well the data only in compositions near optimal doping, indicating that in slightly-doped compositions the damping comes from another mechanism, such as magnon-magnon interactions.

In summary, we have shown that in Ba(FeCo)As the low-energy spin dynamics, which are most strongly tied to excitations in close proximity to the Fermi surface, display a crossover from gapped spin waves to a regime of strong damping and short correlation length, even though weak AFM order persists. The appearance of strong Landau damping near coincides with the appearance of superconductivity, suggesting that the corresponding increase of low energy spectral weight below the spin gap is a key ingredient for superconductivity to develop. In theories where pairing is mediated by spin fluctuations, their energy scale () is usually positively correlated to (3); (2); (1). We instead observe that decreases with increasing (and ), see Fig. 3(d). To avoid the apparent contradiction one must also consider that AFM and SC order compete (29) thereby effectively decreasing for underdoped samples (and eliminating SC for the parent compound). Indeed, when the long-range AFM order is suppresed by pressure, for lower samples is enhanced beyond that for optimal doping(30). Therefore, and are infact positively correlated. Iron pnictide compositions on either side of the SC region — such as Ba(FeCo)As with or (31), or Ba(FeNi)As (32) — lack overdamped spin fluctuations, in contrast to the underdoped SC compositions presented here; this provides further evidence that overdamped spin fluctuations are a necessary component in the paring mechanism for superconductivity in the iron pnictides.

Acknowledgements.
The authors would like to thank D. C. Johnston and A. Kaminski for useful discussions. Work at the Ames Laboratory was supported by the Department of Energy – Basic Energy Sciences, Division of Materials Sciences and Engineering under Contract No. DE-AC02-07CH11358. Part of the research conducted at ORNL’s High Flux Isotope Reactor was sponsored by the Scientific User Facilities Division, Office of Basic Energy Sciences, U.S. Department of Energy. Some work at Oak Ridge (BCS, ASS) was supported by the Department of Energy, Basic Energy Sciences, Materials Sciences and Engineering Division.

References

  1. D. J. Scalapino, Rev. Mod. Phys. 84, 1383 (2012).
  2. A. Abanov, A. V. Chubukov, and J. Schmalian, Adv. Phys. 52, 119 (2003).
  3. P. Monthoux, A. V. Balatsky, and D. Pines, Phys. Rev. B 46, 14803–14817 (1992).
  4. I. I. Mazin, D. J. Singh, M. D. Johannes, and M. H. Du, Phys. Rev. Lett. 101, 057003 (2008); A. V. Chubukov, D. V. Efremov and I Eremin, Phys. Rev. B 78, 134512 (2008); K. Kuroki, S. Onari, R. Arita, H. Usui, Y. Tanaka, H. Kontani, and H. Aoki, Phys. Rev. Lett. 101, 087004 (2008); V. Cvetkovic and Z. Tesanovic, Phys. Rev. B 80, 024512 (2009); J. Zhang, R. Sknepnek, R. M. Fernandes, and J. Schmalian, Phys. Rev. B 79, 220502(R) (2009); A.F. Kemper, T.A. Maier, S. Graser, H-P. Cheng, P.J. Hirschfeld and D.J. Scalapino, New J. Phys. 12, 073030 (2010).
  5. N. Ni, M. E. Tillman, J.-Q. Yan, A. Kracher, S. T. Hannahs, S. L. Bud’ko, and P. C. Canfield, Phys. Rev. B 78, 214515 (2008).
  6. R. M. Fernandes, D. K. Pratt, W. Tian, J. Zarestky, A. Kreyssig, S. Nandi, M. G. Kim, A. Thaler, N. Ni, P. C. Canfield, R. J. McQueeney, J. Schmalian, and A. I. Goldman, Phys. Rev. B 81, 140501(R) (2010).
  7. C. Lester, Jiun-Haw Chu, J. G. Analytis, T. G. Perring, I. R. Fisher, and S. M. Hayden, Phys. Rev. 81, 064505 (2010).
  8. H. F. Li, C. Broholm, D. Vaknin, R. M. Fernandes, D. L. Abernathy, M. B. Stone, D. K. Pratt, W. Tian, Y. Qiu, N. Ni, S. O. Diallo, J. L. Zarestky, S. L. Bud’ko, P. C. Canfield, and R. J. McQueeney, Phys. Rev. B 82, 140503 (2010).
  9. L. W. Harriger, H. Q. Luo, M. S. Liu, C. Frost, J. P. Hu, M. R. Norman, and Pengcheng Dai, Phys. Rev. B 84, 054544 (2011).
  10. Pengcheng Dai, Jiangping Hu, and Elbio Dagotto, Nat. Phys. 8 709 (2012).
  11. G. S. Tucker, R. M. Fernandes, H. F. Li, V. Thampy, N. Ni, D. L. Abernathy, S. L. Bud’ko, P. C. Canfield, D. Vaknin, J. Schmalian, and R. J. McQueeney, Phys. Rev. B 86, 024505 (2012).
  12. Mengshu Liu, Leland W. Harriger, Huiqian Luo, Meng Wang, R. A. Ewings, T. Guidi, Hyowon Park, Kristjan Haule, Gabriel Kotliar, S. M. Hayden, and Pengcheng Dai, Nat. Phys. 8, 376 (2012).
  13. K. Matan, R. Morinaga, K. Iida, and T. J. Sato, Phys. Rev. B 79, 054526 (2009).
  14. R. A. Ewings, T. G. Perring, R. I. Bewley, T. Guidi, M. J. Pitcher, D. R. Parker, S. J. Clarke, and A. T. Boothroyd, Phys. Rev. B 78, 220501(R) (2008).
  15. R. J. McQueeney, S. O. Diallo, V. P. Antropov, G. D. Samolyuk, C. Broholm, N. Ni, S. Nandi, M. Yethiraj, J. L. Zarestky, J. J. Pulikkotil, A. Kreyssig, M. D. Lumsden, B. N. Harmon, P. C. Canfield, and A. I. Goldman, Phys. Rev. Lett. 101, 227205 (2008).
  16. Jun Zhao, Dao-Xin Yao, Shiliang Li, Tao Hong, Y. Chen, S. Chang, W. Ratcliff, J. W. Lynn, H. A. Mook, G. F. Chen, J. L. Luo, N. L. Wang, E. W. Carlson, Jiangping Hu, and Pengcheng Dai, Phys. Rev. Lett. 101, 167203 (2008).
  17. W. Z. Hu, J. Dong, G. Li, Z. Li, P. Zheng, G. F. Chen, J. L. Luo, and N. L. Wang, Phys. Rev. Lett. 101, 257005 (2008).
  18. D. S. Inosov, J. T. Park, P. Bourges, D. L. Sun, Y. Sidis, A. Schneidewind, K. Hradil, D. Haug, C. T. Lin, B. Keimer, and V. Hinkov, Nat. Phys. 6, 178 (2010).
  19. Subir Sachdev, Andrey V. Chubukov, and Alexander Sokol, Phys. Rev. B 51, 14874 (1995).
  20. K Matan, S. Ibuka, R. Morinaga, Songxue Chi, J. W. Lynn, A. D. Christianson, M. D. Lumsden, and T. J. Sato, Phys. Rev. B 82, 054515 (2010).
  21. A. B. Vorontsov, M. G. Vavilov, and A. V. Chubukov, Phys. Rev. B 81, 174538 (2010).
  22. R. M. Fernandes and J. Schmalian, Phys. Rev. B 82, 014521 (2010).
  23. I. Eremin and A. V. Chubukov, Phys. Rev. B 81, 024511 (2010).
  24. R. S. Fishman and S. H. Liu, Phys. Rev. B 58, R5912 (1998).
  25. J. T. Park, G. Friemel, T. Loew, V. Hinkov, Yuan Li, B. H. Min, D. L. Sun, A. Ivanov, A. Piovano, C. T. Lin, B. Keimer, Y. S. Kwon, and D. S. Inosov, Phys. Rev. B 86, 024437 (2012).
  26. M. Ramazanoglu, J. Lamsal, G. S. Tucker, J. Q. Yan, S. Calder, T. Guidi, T. Perring, R. W. McCallum, T. A. Lograsso, A. Kreyssig, A. I. Goldman, and R. J. McQueeney, Phys. Rev. B 87, 140509 (2013).
  27. S. Sachdev, Quantum phase transitions, Cambridge University Press (1999).
  28. C. Liu, T. Kondo, R. M. Fernandes, A. D. Palczewski, E. D. Mun, N. Ni, A. N. Thaler, A. Bostwick, E. Rotenberg, J. Schmalian, S. L. Bud’ko, P. C. Canfield and A. Kaminski, Nat. Phys. 6, 419 (2010).
  29. D. K. Pratt, W. Tian, A. Kreyssig, J. L. Zarestky, S. Nandi, N. Ni, S. L. Bud’ko, P. C. Canfield, A. I. Goldman and R. J. McQueeney, Phys. Rev. Lett. 103, 087001 (2009); A. D. Christianson, M. D. Lumsden, S. E. Nagler, G. J. MacDougall, M. A. McGuire, A. S. Sefat, R. Jin, B. C. Sales and D. Mandrus, Phys. Rev. Lett. 103, 087002 (2009).
  30. E. Colombier, M. S. Torikachvili, N. Ni, S. L. Bud’ko and P. C. Canfield, Supercond. Sci. Technol. 23, 054003 (2010); N. Ni, A. Thaler, J. Q. Yan, A. Kracher, E. Colombier, S. L. Bud’ko, P. C. Canfield, and S. T. Hannahs, Phys. Rev. B 82, 024519 (2010).
  31. T. J. Sato, K. Matan, S. Ibuka, R. Morinaga, Songxue Chi, J. W. Lynn, A. D. Christianson, and M. D. Lumsden, Phys. Rev. B 83, 059901(E) (2011).
  32. Meng Wang, Chenglin Zhang, Xingye Lu, Guotai Tan, Huiqian Luo, Yu Song, Miaoyin Wang, Xiaotian Zhang, E. A. Goremychkin, T. G. Perring, T. A. Maier, Zhiping Yin, Kristjan Haule, Gabriel Kotliar, and Pengcheng Dai, Nat. Commun. 4, 2874 (2013).
123139
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
Edit
-  
Unpublish
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel
Comments 0
Request comment
""
The feedback must be of minumum 40 characters
Add comment
Cancel
Loading ...

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description