Main

The magneto-structural phase transition is among the most thoroughly studied phenomena in iron-based materials. When Fe is substituted by Co in BaFe2As2 the structural transformation at Ts precedes the magnetic ordering at TSDW < Ts (ref. 22). The nematic phase between Ts and TSDW is characterized by broken C4 symmetry but preserved O(3) spin rotational symmetry (no magnetic order). Nematic fluctuations are present even above Ts in the tetragonal phase, as has been demonstrated by both elastoresistance measurements4,6 and studies of the elastic constants23,24. In strained samples, one observes orbital ordering in the photoemission spectra25 and electronic nematicity by transport4,26. However, it is rather difficult to derive the dynamics and critical momentum typical for the underlying fluctuations and to identify which of the ordering phenomena drives the instabilities.

Raman scattering provides experimental access to all types of dynamic nematicity although only the charge sector has been studied in more detail9,14,27,28. However, also in the case of spin-driven nematic order the technique can play a prominent role for coupling to a two-spin operator whereas a four-spin correlation function is the lowest order contribution to the neutron cross-section5. We exploit this advantage here and study the low-energy Raman response of Ba(Fe1−xCox)2As2 experimentally and interpret the results in terms of a microscopic model for a spin-driven nematic phase. In addition to earlier work9,14,20 we study the temperature dependence of the response in the crucial range between Ts and TSDW in great detail and address the spectral shape and the selection rules enabling us to explain the structural and magnetic transitions in a unified microscopic picture.

We study Ba(Fe1−xCox)2As2 single crystals, having x = 0, x = 0.025 and x = 0.051 (see Methods), as a function of photon polarization in the temperature range 4.2 < T ≤ 300 K. For the symmetry assignment we use the 1 Fe unit cell in which the fluctuations appear in B1g symmetry. We use the appearance of twin boundaries and of the As A1g (x2 + y2) phonon line as internal thermometers for the structural and the magnetic phase transitions, respectively. In this way, Ts and TSDW can be determined with a precision of typically ±0.2 and ±1 K, respectively.

Figure 1 shows the Raman response Rχ′′(Ω, T) for Ba(Fe0.975Co0.025)2As2 for various temperatures in A1g and B1g (1 Fe per unit cell) symmetry. B2g spectra were measured only at a few temperatures and found to be nearly temperature independent in agreement with previous data14. Results for other doping levels x are shown in Supplementary Figs 2 and 3. The spectra comprise a superposition of several types of excitation including narrow phonon lines and slowly varying continua arising from electron–hole (e–h) pairs; hence, the continuum reflects the dynamical two-particle behaviour. The A1g and B1g spectra predominantly weigh out contributions from the central hole bands and the electron bands, respectively29. The symmetry-dependent initial slope τ0, μ(T) (μ = A1g, B1g, B2g) (see Fig. 1a, c) can be compared to transport data. [τ0, μ(T)]−1 corresponds to the static transport relaxation rate Γ0, μ(T) of the conduction electrons. The memory function method facilitates the quantitative determination of the dynamic relaxation Γ(Ω, T) in absolute energy units. The static limit can be obtained by extrapolation, Γ0, μ(T) = Γμ(Ω → 0, T) (see Methods and Supplementary Information III). In Fig. 1d we show the result for x = 0.025 corresponding to the spectra of Fig. 1a–c. The results for all doping levels studied are compiled in Supplementary Fig. 4 and compared to the scattering rates derived from the resistivities22.

Figure 1: Raman results of Ba(Fe0.975Co0.025)2As2.
figure 1

a–c, Raman response Rχ′′(Ω, T) (raw data after division by the Bose–Einstein factor) in B1g above (a) and below (b) Ts, and A1g symmetry (c) at temperatures as indicated. The initial slopes shown in a and c as grey arrows are proportional to the static two-particle lifetime in symmetry μ = A1g, B1g. d, Raman relaxation rates Γ0, μ(T) in A1g (blue circles) and B1g (red diamonds) symmetry as a function of temperature (see Methods). The fluctuation range Ts < T < Tf and the nematic phase TSDW < T < Ts are indicated in green and magenta, respectively.

Figure 1d shows one of the central results: above approximately 200 K Γ0, μ(T) varies slowly and similarly in both symmetries. The more rapid decrease of Γ0, B1g(T) below 200 K is accompanied by a strong intensity gain in the range 20–200 cm−1 (see Fig. 1a) as observed before in similar samples9,14,28. The intensity gain indicates that there is an additional contribution superposed on the e–h continuum that, as will be shown below, arises from fluctuations. Therefore, the kink in Γ0, B1g(T) is labelled Tf and marks the crossover temperature below which nematic fluctuations can be observed by Raman scattering. At least for low doping, Tf is relatively well defined. The kink allows us to separate the two regimes of the low-energy response above and below Tf as being dominated by carrier excitations and fluctuations, respectively.

The additional B1g signal below Tf has to be treated in a way different from that in A1g symmetry and in B1g above Tf. As it is rather strong it can be separated out with little uncertainty by subtracting the e–h continuum. We approximate the continuum at Tf by an analytic function that is then determined for each temperature according to the variation of the resistivity and the A1g spectra and subtracted from all spectra at lower temperatures. The details are explained in Supplementary Information VII. The results of the subtraction procedure are shown in Fig. 2. The response increases rapidly towards Ts without however diverging, and the maximum moves to lower energies.

Figure 2: Fluctuation contribution to the Raman response of Ba(Fe0.975Co0.025)2As2.
figure 2

a,b, Results above Ts = 102.8 K (a) and between Ts and TSDW = 98 K (b). The red lines are theoretical predictions on the basis of Aslamazov–Larkin diagrams17 describing the exchange of a pair of fluctuations (for details see Supplementary Information IV and V).

As a surprise, the fluctuations do not disappear directly below Ts (Fig. 2b) as one would expect if long-ranged order would be established. Rather, the intensity decreases continuously and the maximum stays approximately pinned implying that the correlation length does not change substantially between the two transitions at Ts = 102.8 ± 0.2 K and TSDW = 98 ± 1 K. The persistence of the fluctuations down to TSDW strongly favours their magnetic origin.

We first compare the data to the theoretical model for thermally driven spin fluctuations associated with the striped magnetic phase ordering along Qx = (π, 0) or Qy = (0, π). In leading order two non-interacting fluctuations carrying momenta Q and −Q are exchanged. Electronic loops (see Fig. 3 and Supplementary Fig. 5) connect the photons and the fluctuations and entail Q-dependent selection rules, which were derived along with the spectral response R0, μ(Ω) in ref. 17 and are summarized in Supplementary Information IV. In brief, because the response results from a sum over all electronic momenta close to the Fermi surface cancellation effects may occur if Q connects parts on different Fermi surface sheets having form factors γμ(k) with opposite sign. For the ordering vectors (π, 0) and (0, π) the resulting selection rules explain the enhancement of the signal in B1g symmetry and its absence in the A1g and B2g channels.

Figure 3: Higher-order Aslamazov–Larkin diagrams for interacting fluctuations.
figure 3

The paramagnetic coupling between fluctuations mediated by fermions is obtained by inserting quaternions into the leading-order diagram. The re-summed box is shaded grey. The first index of the matrix B denotes the type α = X/Y of entering spin fluctuations, and the second index the type of exiting spin fluctuations.

However, the lowest-order diagrams alone can account only for the spectral shape whereas the variation of the intensity around Ts remains unexplained. To describe this aspect, we consider the interaction of fluctuations among themselves and with the lattice, all of which becomes crucial in the vicinity of the nematic transition19,20.

The interactions between spin fluctuations can be represented by a series of quaternion paramagnetic couplings mediated by fermions inserted into the leading-order Aslamazov–Larkin diagrams as shown in Fig. 3. The inserted fermionic boxes effectively resemble the dynamic nematic coupling constant g of the theory.

We have analysed the problem by extending SU(2) → SU(N) and taking the large N limit. For small frequencies Ω and in the large-N limit, after re-summing an infinite number of such box-like Aslamazov–Larkin diagrams, the Raman response function  reads:

Equation (1) states that the Raman response is proportional to the electronic contribution to the susceptibility of the nematic order parameter,

χmag(q) represents the magnetic susceptibility, which diverges at TSDW. For g ≠ 0 χnemel(0) has a Curie-like |T − T∗|−1 divergence at T∗ ≥ TSDW.

If the spins (or charges) couple to the lattice the susceptibility of the nematic order parameter is given by4,16,20

where λsl denotes the magneto-elastic coupling, and c0s is the bare elastic constant. Obviously, χnem(0) diverges at higher temperature than χnemel(0). We identify Ts ≥ T∗ with the structural transition and conclude that the Raman response (equation (1)) develops only a maximum rather than a divergence at Ts in agreement with the experiment here and recent theoretical work19,20.

Close to Ts, we expect equation (1) to hold qualitatively also inside the nematic phase, TSDW < T < Ts. We argue13 that χnemel(0) and, according to equation (1), the Raman amplitude are smaller than in the disordered (tetragonal) state but different from zero. This explains the continuous reduction of the Raman response of spin fluctuations on entering the nematic state. One can also show that the A1g response gets even further suppressed if one includes collisions between the fluctuations20.

As shown in equation (1) the full Raman response is proportional to the bare response R0, μ(Ω) and to the electronic nematic susceptibility χnemel(0). Hence, the spectral shape is essentially given by R0, μ(Ω), which is therefore used in Fig. 2 to fit the data, whereas the intensity is dominated by the prefactor |T − T∗|−1. As the theoretical model is valid only in the limit of small frequencies we argue that the initial slope reflects the temperature dependence of the intensity and is proportional to χnemel(0), at least close to the transition. For generally reflecting the spectral shape above TSDW (equation (1)), R0, B1g(Ω, T) enables us to directly extract the initial slope of the experimental spectra by plotting R0, B1g(Ω, T)/Ω for all temperatures (Supplementary Information VI). These results are compiled in Fig. 4 along with the variation of χnemel(0, T) expected from mean-field theory. For low doping, we find qualitative agreement in the ranges T > TSDW. For higher doping the interactions between fluctuations become dominant and the mean-field prediction breaks down (Fig. 4c).

Figure 4: Temperature dependence of the initial slope of the fluctuation response.
figure 4

a–c, The initial slope is given in arbitrary units because only the temperature dependence matters. The error bars in b represent the differences originating in the subtraction of the e–h continuum. They are large below Ts because of surface contaminations and smaller than the size of the points above. The grey curves in a,b represent |T − T∗|−1 with T∗ < Ts (see text).

The Raman response was also studied at various other doping levels in the range 0 ≤ x ≤ 0.085. Up to 6.1% Co substitution fluctuations were observed. In contrast to other publications14 we were not able to clearly identify and isolate the response of fluctuations at 8.5% although the kink in the B1g relaxation rate used to define Tf is clearly observed (Supplementary Fig. 4). The results for χnemel(0) up to 5.1% are unambiguous and are represented as a colour scale on the phase diagram in Fig. 5. Our phase diagram compares rather well to that derived from the elastic constant m66 (ref. 6). In addition, we show Tf up to x = 0.085. The fluctuations can be observed over a temperature range of approximately 70–100 K. This is more than in most of the other experiments on unstrained samples and comparable to what is found in the cuprates18,30.

Figure 5: Phase diagram of Ba(Fe1−xCox)2As2.
figure 5

The full lines limiting the nematic phase (magenta) and the blue squares representing the transition temperature Tc of superconducting samples were derived in ref. 22. Grey diamonds represent doping and temperature positions of the current Raman data. The coloured field between Ts and Tf represents the initial slope of the spectra according to the colour scale on the right. Tf is a crossover and corresponds to χnemel(0) → 0.

Methods

Samples.

The single crystals of undoped31 and Co-substituted Ba(Fe1−xCox)2As2 were grown using a self-flux technique and have been characterized elsewhere22. The cobalt concentration was determined by microprobe analysis. Ts and TSDW are close to 134 K in the undoped sample and cannot be distinguished. At nominally x = 0.025 we find Ts = 102.8 ± 0.1 K and TSDW = 98 ± 1 K by directly observing the appearance of twin boundaries and a symmetry-forbidden phonon line, respectively (see Supplementary Information I for details). The extremely sharp transition at Ts having ΔTs ≍ 0.2 K indicates that the sample is very homogeneous in the area of the laser spot.

Raman experiment.

The experiments were performed with standard light scattering equipment. For excitation either a solid-state laser (Coherent, Sapphire SF 532-155 CW) or an Ar ion laser (Coherent, Innova 300) was used emitting at 532 or 514.5 nm, respectively. The samples were mounted on the cold finger of a He-flow cryostat in a cryogenically pumped vacuum. The laser-induced heating was determined experimentally (see Supplementary Information I) to be close to 1 K per milliwatt of absorbed power. The spectra represent the response Rχμ′′(Ω, T) (μ = A1g, B1g, A2g and B2g), which is obtained by dividing the measured (symmetry resolved) spectra by the Bose–Einstein factor {1 + n(T, Ω)} = [1 − exp(−ℏΩ/kBT)]−1. χμ′′(Ω, T) is the imaginary part of the response function, and R is an experimental constant that connects the observed photon count rates with the cross-section and the van Hove function and accounts for units. For simplicity, the symmetry index μ is dropped in most of the cases. The symmetry selection rules refer to the 1 Fe unit cell (see Supplementary Fig. 1), which is more appropriate for electronic and spin excitations.

Static relaxation rates.

Static relaxation rates Γ0, μ(T), with μ denoting the symmetry, are used in various places (see, for example, Fig. 1b),

Γ0, μ(T) can be considered a symmetry-resolved ‘Raman resistivity’ that can be compared to, for example, transport data32. As the scattering intensity and consequently R are not known in absolute units Γ0, μ(T) cannot directly be derived. The problem was solved a while ago by adopting the memory function method33,34 for Raman scattering35. Then Γ0, μ(T) can be derived by extrapolating the dynamic Raman relaxation rates Γμ(Ω, T) to zero energy. More details can be found in Supplementary Information III.