## Evolution of the Nematic Susceptibility in LaFe 1 -x Co AsO x

Xiaochen Hong 1 , ∗ , Federico Caglieris 1 , Rhea Kappenberger , Sabine Wurmehl , Saicharan Aswartham , Francesco 1 1 1 Scaravaggi 1 2 , , Piotr Lepucki , Anja U. B. Wolter , Hans-Joachim Grafe , Bernd B¨chner 1 1 1 u 1 2 3 4 , , , and Christian Hess 1 3 , , † 1 Leibniz-Institute for Solid State and Materials Research (IFW-Dresden), 01069 Dresden, Germany 2 Institute of Solid State and Materials Physics, TU Dresden, 01069 Dresden, Germany 3 Center for Transport and Devices, TU Dresden, 01069 Dresden, Germany 4 W¨ urzburg-Dresden Cluster of Excellence ct.qmat , TU Dresden, 01062 Dresden, Germany

(Dated: August 4, 2020)

We report a systematic elastoresistivity study on LaFe 1 -x Co AsO single crystals, which have well x separated structural and magnetic transition lines. All crystals show a Curie-Weiss-like nematic susceptibility in the tetragonal phase. The extracted nematic temperature is monotonically suppressed upon cobalt doping, and changes sign around the optimal doping level, indicating a possible nematic quantum critical point beneath the superconducting dome. The amplitude of the nematic susceptibility shows a peculiar double-peak feature. This could be explained by a combined effect of different contributions to the nematic susceptibility, which are amplified at separated doping levels of LaFe 1 -x Co AsO. x

Unconventional superconductivity is intimately related to other electronic symmetry breaking states. A common feature among unconventional superconductors is that by varying a tuning parameter x , such as pressure or chemical substitution, the superconductivity exists under a dome-like region in the T -x phase diagram, and the extrapolation of some transition line hits zero temperature inside or close to that dome [1-4]. The formation of a spin density wave is the most recognized order in the normal state of unconventional superconductors, and its ubiquitousness stimulated theories to explain the origin of unconventional superconductivity mediated by fluctuations of the magnetic order [5-10]. Recently, the report of newly recognized electronic orders beside magnetism in unconventional superconductors is infectious [2, 11], including nematicity, an electronic ordered state that spontaneously breaks the rotational symmetry of its host crystal [11-13]. Nematic fluctuation by themselves are considered as possible mediator of electronic pairing [14-18]. Indeed, a nickel pnictide Ba 1 -x Sr x Ni As 2 2 was just discovered to be a nematicity-boosted superconductor [19].

ing force for other transitions [23, 24].

In the case of iron-based superconductors, the stoichiometric parent compounds are always antiferromagnetically (AFM) ordered below T N , following a structural transition at T S [3, 4]. The paramagnetic orthorhombic phase between T S and T N acquires highly anisotropic electronic properties regardless of the small lattice anisotropy [20]. In this electronic nematic phase, all ordered states are highly intertwined, making it hard to discern the leading instability among them. One can find clues by studying their pertinent fluctuations at higher temperatures [21-23]. For nematicity, its fluctuations have been probed by various approaches, including the elastoresistivity measurements [24, 25]. The divergent nematic susceptibility has been found in different families of iron-based compounds [26-28]. These facts clearly eliminated the structural instability as the driv-

One prevailing understanding of the nematicity in ironbased compounds is to treat it as the result of fluctuating AFM order [29-32]. This scenario received strong support from the discovery of a good scaling between the lattice softening and the magnetic fluctuating amplitude in Ba(Fe 1 -x Co ) As x 2 2 over a wide doping range [33]. As illustrated in Fig. 1(a), the T S and T N transition lines are very close to each other in Ba(Fe 1 -x Co ) As . x 2 2 The extension of these lines penetrate into a sign-reversed s -wave superconducting dome [3], which is the fingerprint of AFM fluctuation driven superconductivity [34, 35]. In stark contrast, as Fig. 1(b) shows, there is no static magnetic order over the whole phase diagram of FeSe 1 -x S . x Theories proposed a different kind of superconductivity of sign-preserved pairing, based on unequal orbital occupancy [36-40]. The nuclear magnetic resonance measurements saw no slowing down of spin fluctuations above T S in FeSe and much smaller changes of spin-lattice relaxation rates across T S than in other iron-based superconductors [41, 42], lent support to the orbital scenarios and alternative origin of nematicity. However, subsequent neutron scattering studies can identify spin fluctuations above T S [43]. Moreover, these spin fluctuations are found changing across T S and being anisotropic in the nematic phase [44, 45]. Thus, the importance of magnetic fluctuations is restored in FeSe. The leading instability and therefore the primary pairing mediator of the seemingly simple FeSe system remain elusive.

In this letter, we report the study of nematicity of a series of LaFe 1 -x Co AsO single crystals. x As Fig. 1(c) indicates, the T S and T N of LaFe 1 -x Co AsO are well x separated [46-48]. Superconductivity emerges after the total suppression of the AFM order [47, 49], making LaFe 1 -x Co AsO interesting to search for separated magx netic and nematic fluctuation regions. By utilizing elastoresistivity measurements, the persistence of a CurieWeiss-like nematic susceptibility deep in the tetragonal

Temperature

FIG. 1: (Color online). Schematic T -x phase diagrams of three representative iron based superconductors [3, 27, 47, 49]: (a) Ba(Fe 1 -x Co ) As , (b) FeSe x 2 2 1 -x S , (c) LaFe x 1 -x Co AsO. x The orange area is the AFM ordered state, the green dome is the superconducting state, and the solid blue line shows the structural transition (with dashed line as a guide to the eye). The nematic order exists below the structural transition line and above the AFM region.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00484v2-figure-1ca9b6ba-65e1-400f-939c-32e7271b9d9f.jpg)

phase of LaFe 1 -x Co AsO was found. x The Curie-Weiss temperature changes sign around the optimal doping level x = 0.06, with an enhanced nematic susceptibility on its top. This indicates the existence of a nematic quantum critical point (nQCP) under the superconducting dome, and supports the proposals of nematicity as the driver for superconductivity [14-18]. Besides, another enhancement of the nematic susceptibility is resolved around x = 0.04, close to the end point of the AFM order.

LaFe 1 -x Co AsO single crystals were prepared by the x solid state single crystal growth method [46]. The oriented samples were cut into a typical size of 1.2 mm × 0.5 mm in plane and cleaved to around 20 µ m thickness to ensure an efficient strain transmission to the sample [24]. The experimental setup is illustrated in the insert of Fig. 2(c). The sample was glued to a commercial piezoelectric actuator (Piezomechanik PSt 150 / 5 × × 5 7) by using an adhesive epoxy (Devcon, No. 14250), after electrically contacting voltage and current leads directly on the fresh surface with sliver paint. A strain gauge was glued on the other side of the actuator and was measured with a built-in Wheatstone bridge circuit. To eliminate a possible temperature change caused by driving the piezoelectric actuator, data were taken after waiting several seconds at each strain step.

Elastoresistivity measured along the [110] direction of the LaFeAsO single crystal is shown in Fig. 2. By changing the voltage across a piezo actuator, one can precisely tune the strain applied to the sample glued on the actuator's edge. In the nematic ordered state, the resistivity anisotropy of the in-plane axes is very large [50, 51]. This anisotropy represents the electronic nematicity. Due to the electron-lattice coupling, strain can induce a resistivity anisotropy above the nematic ordered state if nematic fluctuations exist [24]. As illustrated in Fig. 2 (b), the fractional change of sample resistivity ( η = ∆ ρ/ρ ) shows a perfect linear relationship to the strain ( ε x = ∆ L /L x x )

FIG. 2: (Color online). Elastoresistance measurements of the parent compound LaFeAsO. (a) The relative change of strain ( ε x = ∆ L /L x x ) and sample resistivity ( η = ∆ ρ/ρ ) according to the voltage applied to the piezo actuator at a fixed temperature T = 150 K. The strain was measured along the current direction, which was aligned against the [110] axis of the crystal. (b) The relationship between ε x and η at some representative temperatures. The nematic susceptibility was obtained from ˜ = n -( δη/δε x ) in the ε x = 0 limit. (c) The temperature dependence of ˜ n is shown as blue open squares in the upper panel. The structural transition temperature determined in Ref. [46] is indicated by the green arrow. The red curve represents a Curie-Weiss fitting to data in the temperature range between the vertical dotted lines. Inset is a schematic of the experimental setup. The inverse nematic susceptibility (˜ n -˜ n 0 ) -1 and the Curie constant C = (˜ n -˜ n 0 )( T -T nem ) are shown in the lower panel, indicating the fitting quality.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00484v2-figure-655f3f87-f897-44b2-a8b5-dd10ca207aec.jpg)

applied along the current direction. As demonstrated in the references [24-27], the slope of η ε ( x ) curve probes the nematic susceptibility of the sample. The negative slope means that the in-plane resistivity of LaFeAsO is higher along the shorter orthorhombic direction. For simplicity, we use -( δη/δε x ) in the small strain limit as the definition of ˜, which is a gauge of the nematic susceptibility. n The temperature dependence of ˜ in LaFeAsO is shown n in Fig. 2(c). A clear kink is resolved at the structural transition T S , above which the ˜( n T ) curve can be well fitted to the Curie-Weiss type temperature dependence

$$\tilde { n } = \tilde { n } _ { 0 } + \frac { C } { T - T ^ { n e m } },$$

in which ˜ n 0 is the intrinsic piezoresistivity effect unrelated to electronic nematicity, T nem is the nematic transition temperature in the mean field theory, and C is the Curie constant which indicates the magnitude of the nematic susceptibility.

Similar measurements and analyses were conducted for

FIG. 3: (Color online). Divergent nematic susceptibility in LaFe 1 -x Co AsO single crystals. Similar data as Fig. 2(c) for eleven x samples of different doping level ranging from underdoped to overdoped is shown in (a - k). Cobalt content x is the nominal composition, which agrees well with the actual concentration determined by energy dispersive X-ray spectroscopy on crystals from the same batches. The temperatures between which Curie-Weiss fitting was performed were determined by minimizing the systematic deviation of the Curie constant (see Ref. [26]). The structural transition temperature is indicated by the green arrow and is extracted from the references [48, 49]. There is no arrow in Fig. 3(e) because the Co 4.2% sample was not studied in the references [48, 49].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00484v2-figure-22bede90-96fe-4045-8a4f-f95d909ac478.jpg)

LaFe 1 -x Co AsO single crystals across the substitution x series ( x up to 0.075). As shown in Fig. 3, for all the samples we have studied, their ˜( n T ) curves diverge in a clear Curie-Weiss-like form. The kink at T S becomes less pronounced with doping and is unresolvable in our data set for x &gt; 0 03, . while T S could be traced up to x = 0 045 in thermal expansion measurements (indicated . by green arrows) [49]. One may attribute the smearing of the T S feature to the increasing impurity scattering effect. In such a case, according to Ref. [52], ˜ should n remain practically unaffected by the defects and impurity concentrations in the tetragonal phase. However, the blurred feature of T S could also be an intrinsic property of LaFe 1 -x Co AsO, since the structural transition x anomaly observed in the thermal expansion experiments becomes broader with increasing x [49]. Indeed, a nuclear quadrupole resonance study on LaFe 1 -x Co AsO showed x that competing local charge environments exist at the nanoscale in the underdoped region [47], similar to the sister compound LaFeAsO 1 -x F x [53, 54]. In that sense, it is natural that the features of T S from averaged properties are muddled by doping.

fitting range of the ˜( n T ) curves, especially for those without (a clear) T S , we followed the procedure described in Ref. [26]. Additionally, the fitting range is allowed to alter by 20 K for each sample in order to estimate the uncertainty caused by improper fitting ranges [55], set as error bars of the fitted parameters T nem and C in the following discussion.

We analyze the ˜( n T ) data for the finite doping levels in the same Curie-Weiss manner as for the parent compound (solid red lines). In order to choose the proper

The information extracted from Fig. 3 is summarized in Fig. 4, plotted together with the phase diagram of LaFe 1 -x Co AsO derived from the doping dependence of x T S , T N , and T c [47, 49]. First of all, one can see that the nematic transition temperature T nem is linearly suppressed by doping, and finally changes its sign around the optimal doping level x ≈ 0 06. . It is known that superconductivity could get boosted via fluctuations related to electronic instabilities. The most studied case is that of AFM fluctuations with certain wave vectors around a magnetic critical point. As shown in Fig. 1(c), the magnetic transition is totally suppressed before the superconductivity emerges in LaFe 1 -x Co AsO [47, 49], suggesting x that the magnetic fluctuations near the critical point may not contribute to the superconducting transition temperature T c directly. Theoretically, a nQCP and its pertinent fluctuation can also promote the superconductivity via the attractive forces in all paring channels [14-17],

0

50

100

150

200

250

300

(K)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00484v2-figure-93ef03cf-9720-46e8-a51f-5dd455d4ca64.jpg)

T

FIG. 4: (Color online). Phase diagram and nematic susceptibility of LaFe 1 -x Co AsO. The white triangles, x light blue diamonds, and green circles indicate T S , T N , and T c respectively, according to the references [47-49]. The orange area encloses the AFM ordered state, while the superconducting region is represented by the green area. The nematic critical temperature T nem (red star) and the divergency coefficient (Curie constant) C (magenta square, right axis, in log scale) are obtained by the Curie-Weiss fits of the nematic susceptibility data shown in Fig. 2 and Fig. 3. The error bars of them come from altering the fit range (see text). The dashed lines are guides to the eye. Evolution of the magnitude of ˜ is n shown by the color plot, which shows a double-peak feature.

or by reducing the intra-pocket repulsion [18, 22]. Our observation of the sign change of T nem around x = 0.06 strongly corroborates those theoretical expectations.

Interestingly, upon investigating the doping dependence of the magnitude of ˜, n one can clearly observe a significant enhancement of this quantity at the optimal doping level. This can be inferred from the data in the upper panels of Fig. 3(a-k) and also the contour plot presented in Fig. 4, revealing that the enhanced ˜ n mimics the superconducting dome. The latter trend is also clearly visible in the evolution of the Curie constant displayed in Fig. 4, too, which represents the divergency coefficient of the ˜( n T ) curve. This finding of an enhanced ˜ on top of the critical point of n T nem and maximum T c strongly underpins the above notion of a putative nQCP near the optimal doping level.

The further inspection of ˜( n T ) in Figs. 3 and 4 at lower doping clearly reveals a second peak-like enhancement of it, very close to the magnetic critical point where the AFM order disappears but still at doping levels with an AFM ground state. An alternative nematic critical point as the origin of the additional peak in the under- doped region at x ≈ 0.04 can be ruled out because the fitting of ˜( n T ) gives unambiguously a finite T nem ≈ 35 K at x ≈ 0.04 [55]. One may speculate that the enhanced ˜ around n x = 0.04 is caused by magnetic fluctuations, which are critically enhanced at the doping levels just below the disappearance of the AFM order (like in many other iron-based superconductors [4]), and remain small up to the optimal doping level x = 0.06 [47]. This notion is supported by the established T nem ≈ 35 K for x ≈ 0.04 which roughly matches the observed magnetic ordering temperature T N around this doping level. Thus, altogether our data seem to reveal that in the underdoped region of the phase diagram for x ≤ 0.05, ˜ foln lows the trend of the magnetic fluctuations which are connected to the incipient order. At higher doping, ˜ is n enhanced and seems decoupled from the magnetic fluctuations, while on top of the presumptive nQCP, ˜ and n T c domes exist despite a reduced magnetic susceptibility. It is worthwhile to point out that the seeming appearance of two separated criticalities suggests nematicity itself can exist as a primary fluctuation and may be responsible for the superconductivity. Furthermore, the comparable Curie constant C around the two critical doping levels indicates a similar strength of the elasto-electronic and elasto-magnetic couplings in this material. This is an interesting point to look at in further investigations.

The dichotomic-origin of the ˜ in LaFe n 1 -x Co AsO prox vides fresh input for rationalizing the importance of nematic fluctuations in other members of the iron-based superconductor family. In electron-doped BaFe As , the 2 2 peak of ˜ n has a broad appearance with the tendency to maximize at a slightly underdoped level [24, 26, 56]. Since the magnetic and structural transition lines of those materials are close to each other, the measured ˜ nR , a superposition of the two profiles which maximize at slightly separated end points, could appear like a broadened peak in between. We note that FeSe 1 -x S x has recently been proven to be a remarkable system which possesses an nQCP under its superconducting region [27, 57]. Although no AFM order exists in its ambient pressure phase diagram, magnetic fluctuations are still interpreted as the booster of superconductivity in FeSe 1 -x S x [58]. As we have mentioned above, however, magnetic fluctuations are suppressed at the superconducting doping levels in LaFe 1 -x Co AsO, suggesting it is a more suitable protox type system to search for the relationship between superconductivity and nematicity. Finally, we point out that recent reports also claimed separated QCPs masked below the wide-spread superconducting dome in cobalt- and nickel-doped NaFeAs, extracted from the nuclear magnetic resonance and neutron scattering spectra respectively [59, 60]. We therefore suspect a similar doublepeak ˜ profile would also be detected in the NaFeAs famn ily.

To summarize, we found that the Curie-Weiss-like ˜ n also exists in LaFe 1 -x Co AsO, in line x with the other

iron-based superconductors. The sign-change of T nem around x = 0.06 and the divergent amplitude of ˜ n by approaching this doping level are consistent with a possible nQCP at the optimal doping level. An additional peak of ˜ in the underdoped region around n x ≈ 0.04 is present close to the end point of the AFM transition line in LaFe 1 -x Co AsO, which we attribute to a coupling of x the nematic fluctuation to critical magnetic fluctuations. Altogether, the double-peak feature of ˜ n suggests another origin of nematicity in iron-based compounds, besides the well accepted vestigial magnetism explanation [29, 30]. A detailed map out of the magnetic fluctuations in LaFe 1 -x Co AsO is important to further understand x nematicity in this material.

We would like to thank R¨diger Klingeler, Sven Saueru land, Liran Wang, Seung-Ho Baek, Christoph Wuttke, and Steffen Sykora for helpful discussions. This work has been supported by the Deutsche Forschungsgemeinschaft (DFG) through the SFB 1143 (project-id 247310070), through the Priority Programme SPP1458 (Grant No.BU887/15-1), under grant DFG-GRK 1621, through DFG Research Grants AS 523/4-1 (S.A.) CA1931/1-1 (F.C.), and through the Emmy Noether Programme WU595/3-3 (S.W.). X.C.H. acknowledges the support of the Leibniz-DAAD Research Fellowship. This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 647276-MARS-ERC-2014-CoG) and the W¨ urzburg-Dresden Cluster of Excellence on Complexity and Topology in Quantum Matter ct.qmat (EXC 2147, project-id 390858490).

- ∗ E-mail: x.c.hong@ifw-dresden.de
- † E-mail: c.hess@ifw-dresden.de
- [1] C. Pfleiderer, Superconducting phases of f-electron compounds , Rev. Mod. Phys. 81 , 1551 (2009).
- [2] B. Keimer, S. A. Kivelson, M. R. Norman, S. Uchida, and J. Zaanen, From quantum matter to high-temperature superconductivity in copper oxides , Nature 518 , 179 (2015).
- [3] G. R. Stewart, Superconductivity in iron compounds , Rev. Mod. Phys. 83 , 1589 (2011).
- [4] P. C. Dai, Antiferromagnetic order and spin dynamics in iron-based superconductors , Rev. Mod. Phys. 87 , 855 (2015).
- [5] Ar. Abanov, A. V. Chubukov, and J. Schmalian, Quantum-critical theory of the spin-fermion model and its application to cuprates: Normal state analysis , Adv. Phys. 52 , 119 (2003).
- [6] P. Monthoux, D. Pines, and G. G. Lonzarich, Superconductivity without phonons , Nature 450 , 1177 (2007).
- [7] H. v. L¨ ohneysen, A. Rosch, M. Vojta, and P. W¨ olfle, Fermi-liquid instabilities at magnetic quantum phase transitions , Rev. Mod. Phys. 79 , 1015 (2007).
- [8] L. Taillefer, Scattering and pairing in cuprate superconductors , Annu. Rev. Condens. Matter Phys. 1 , 51 (2010).
- [9] D. J. Scalapino, A common thread: The pairing interaction for unconventional superconductors , Rev. Mod. Phys. 84 , 1383 (2012).
- [10] S. Sachdev, M. A. Metlitski, and M. Punk, Antiferromagnetism in metals: from the cuprate superconductors to the heavy fermion materials , J. Phys.: Condens. Matter 24 , 294205 (2012).
- [11] E. Fradkin, S. A. Kivelson, M. J. Lawler, J. P. Eisenstein, and A. P. Mackenzie, Nematic Fermi fluids in condensed matter physics , Annu. Rev. Condens. Matter Phys. 1 , 151 (2010).
- [12] S. A. Kivelson, E. Fradkin, and V. J. Emery, Electronic liquid-crystal phases of a doped Mott insulator , Nature 393 , 550 (1998).
- [13] M. P. Lilly, K. B. Cooper, J. P. Eisenstein, L. N. Pfeiffer, and K. W. West, Anisotropic states of two-dimensional electron systems in high Landau levels: effect of an inplane magnetic field , Phys. Rev. Lett. 83 , 824 (1999).
- [14] T. A. Maier and D. J. Scalapino, Pairing interaction near a nematic quantum critical point of a three-band CuO 2 model , Phys. Rev. B 90 , 174510 (2014).
- [15] S. Lederer, Y. Schattner, E. Berg, and S. A. Kivelson, Enhancement of superconductivity near a nematic quantum critical point , Phys. Rev. Lett. 114 , 097001 (2015).
- [16] M. A. Metlitski, D. F. Mross, S. Sachdev, and T. Senthil, Cooper pairing in non-Fermi liquids , Phys. Rev. B 91 , 115111 (2015).
- [17] D. Labat and I. Paul, Pairing instability near a latticeinfluenced nematic quantum critical point , Phys. Rev. B 96 , 195146 (2017).
- [18] D. L. Maslov and A. V. Chubukov, Fermi liquid near Pomeranchuk quantum criticality , Phys. Rev. B 81 , 045110 (2010).
- [19] C. Eckberg, D. J. Campbell, T. Metz, J. Collini, H. Hodovanets, T. Drye, P. Zavalij, M. H. Christensen, R. M. Fernandes, S. Lee, P. Abbamonte, J. W. Lynn , and J. Paglione, Sixfold enhancement of superconductivity in a tunable electronic nematic system , Nat. Phys. 16 , 346 (2020).
- [20] I. R. Fisher, L. Degiorgi, and Z. X. Shen, In-plane electronic anisotropy of underdoped '122' Fe-arsenide superconductors revealed by measurements of detwinned single crystals , Rep. Prog. Phys. 74 , 124506 (2011).
- [21] A. Cano, M. Civelli, I. Eremin, and I. Paul, Interplay of magnetic and structural transitions in iron-based pnictide superconductors , Phys. Rev. B 82 , 020408(R) (2010).
- [22] R. M. Fernandes and J. Schmalian, Manifestations of nematic degrees of freedom in the magnetic, elastic, and superconducting properties of the iron pnictides , Supercond. Sci. Technol. 25 , 084005 (2012).
- [23] R. M. Fernandes, A. V. Chubukov, and J. Schmalian, What drives nematic order in iron-based superconductors? Nat. Phys. 10 , 97 (2014).
- [24] J.-H. Chu, H.-H. Kuo, J. G. Analytis, and I. R. Fisher, Divergent nematic susceptibility in an iron arsenide superconductor , Science 337 , 710 (2012).
- [25] M. C. Shapiro, A. T. Hristov, J. C. Palmstrom, J.-H. Chu, and I. R. Fisher, Measurement of the B 1 g and B 2 g components of the elastoresistivity tensor for tetragonal materials via transverse resistivity configurations , Rev. Sci. Instrum. 87 , 063902 (2016).
- [26] H.-H. Kuo, J.-H. Chu, J. C. Palmstrom, S. A. Kivelson, and I. R. Fisher, Ubiquitous signatures of nematic quantum criticality in optimally doped Fe-based superconduc-

tors , Science 352 , 958 (2016).

- [27] S. Hosoi, K. Matsuura, K. Ishida, H. Wang, Y. Mizukami, T. Watashige, S. Kasahara, Y. Matsuda, and T. Shibauchi, Nematic quantum critical point without magnetism in FeSe 1 -x S x superconductors , Proc. Natl. Acad. Sci. USA 113 , 8139 (2016).
- [28] Y. H. Gu, Z. Y. Liu, T. Xie, W. L. Zhang, D. L. Gong, D. Hu, X. Y. Ma, C. H. Li, L. X. Zhao, L. F. Lin, Z. Xu, G. T. Tan, G. F. Chen, Z. Y. Meng, Y. F. Yang, H. Q. Luo, and S. L. Li, Unified phase diagram for iron-based superconductors , Phys. Rev. Lett. 119 , 157001 (2017).
- [29] C. Fang, H. Yao, W.-F. Tsai, J. P. Hu, and S. A. Kivelson, Theory of electron nematic order in LaFeAsO, Phys. Rev. B 77 , 224509 (2008).
- [30] A. V. Chubukov, R. M. Fernandes, and J. Schmalian, Origin of nematic order in FeSe, Phys. Rev. B 91 , 201105(R) (2015).
- [31] F. Kretzschmar, T. B¨ ohm, U. Karahasanovi´ c, B. Muschler, A. Baum, D. Jost, J. Schmalian, S. Caprara, M. Grilli, C. Di Castro, J. G. Analytis, J.-H. Chu, I. R. Fisher, and R. Hackl, Critical spin fluctuations and the origin of nematic order in Ba(Fe 1 -x Co ) As , Nat. x 2 2 Phys. 12 , 560 (2016).
- [32] X. Y. Lu, J. T. Park, R. Zhang, H. Q. Luo, A. H. Nevidomskyy, Q. M. Si, and P. C. Dai, Nematic spin correlations in the tetragonal state of uniaxial-strained BaFe 2 -x Ni As , Science x 2 345 , 657 (2014).
- [33] R. M. Fernandes, A. E. B¨ ohmer, C. Meingast, and J. Schmalian, Scaling between magnetic and lattice fluctuations in iron pnictide superconductors , Phys. Rev. Lett. 111 , 137001 (2013).
- [34] S. L. Li, Y. Chen, S. Chang, J. W. Lynn, L. J. Li, Y. K Luo, G. H. Cao, Z. A. Xu, and P. C. Dai, Spin gap and magnetic resonance in superconducting BaFe 1 9 . Ni 0 1 . As , 2 Phys. Rev. B 79 , 174527 (2009).
- [35] 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, Normal-state spin dynamics and temperature-dependent spin-resonance energy in optimally doped BaFe 1 85 . Co 0 15 . As , 2 Nat. Phys. 6 , 178 (2010).
- [36] F. Kr¨ger, S. Kumar, J. Zaanen, and J. van den Brink, u Spin-orbital frustrations and anomalous metallic state in iron-pnictide superconductors , Phys. Rev. B 79 , 054504 (2009).
- [37] W. C. Lv, J. S. Wu, and P. Phillips, Orbital ordering induces structural phase transition and the resistivity anomaly in iron pnictides , Phys. Rev. B 80 , 224506 (2009).
- [38] C.-C. Lee, W. G. Yin, and W. Ku, Ferro-orbital order and strong magnetic anisotropy in the parent compounds of iron-pnictide superconductors , Phys. Rev. Lett. 103 , 267001 (2009).
- [39] C.-C. Chen, J. Maciejko, A. P. Sorini, B. Moritz, R. R. P. Singh, and T. P. Devereaux, Orbital order and spontaneous orthorhombicity in iron pnictides , Phys. Rev. B 82 , 100504(R) (2010).
- [40] H. Kontani and S. Onari, Orbital-fluctuation-mediated superconductivity in iron pnictides: analysis of the fiveorbital Hubbard-Holstein model , Phys. Rev. Lett. 104 , 157001 (2010).
- [41] S.-H. Baek, D. V. Efremov, J. M. Ok, J. S. Kim, J. van den Brink, and B. B¨chner, u Orbital-driven nematicity in FeSe, Nat. Mater. 14 , 210 (2015).
- [42] A. E. B¨ ohmer, T. Arai, F. Hardy, T. Hattori, T. Iye, T. Wolf, H. v. L¨hneysen, K. Ishida, and C. Meingast, o Origin of the tetragonal-to-orthorhombic phase transition in FeSe: a combined thermodynamic and NMR study of nematicity , Phys. Rev. Lett. 114 , 027001 (2015).
- [43] Q. S. Wang, Y. Shen, B. Y. Pan, Y. Q. Hao, M. W. Ma, F. Zhou, P. Steffens, K. Schmalzl, T. R. Forrest, M. Abdel-Hafiez, X. J. Chen, D. A. Chareev, A. N. Vasiliev, P. Bourges, Y. Sidis, H. B. Cao, and J. Zhao, Strong interplay between stripe spin fluctuations, nematicity and superconductivity in FeSe, Nat. Mater. 15 , 159 (2016).
- [44] Q. S. Wang, Y. Shen, B. Y. Pan, X. W. Zhang, K. Ikeuchi, K. Iida, A. D. Christianson, H. C. Walker, D. T. Adroja, M. Abdel-Hafiez, X. J. Chen, D. A. Chareev, A. N. Vasiliev, and J. Zhao, Magnetic ground state of FeSe, Nat. Commun. 7 , 12182 (2016).
- [45] T. Chen, Y. Z. Chen, A. Kreisel, X. Y. Lu, A. Schneidewind, Y. M. Qiu, J. T. Park, T. G. Perring, J R. Stewart, H. B. Cao, R. Zhang, Y. Li, Y. Rong, Y. Wei, B. M. Andersen, P. J. Hirschfeld, C. Broholm, and P. C. Dai, Anisotropic spin fluctuations in detwinned FeSe, Nat. Mater. 18 , 709 (2019).
- [46] R. Kappenberger, S. Aswartham, F. Scaravaggi, C. G. F. Blum, M. I. Sturza, A. U. B. Wolter, S. Wurmehl, and B. B¨chner, u Solid state single crystal growth of threedimensional faceted LaFeAsO crystals , J. Cryst. Growth 483 , 9 (2018).
- [47] P. Lepucki, H.-J. Grafe et al. , NMR and NQR measurements of LaFe 1 -x Co AsO, unpublished. x
- [48] L. Wang, S. Sauerland, F. Scaravaggi, R. Kappenberger, S. Aswartham, S. Wurmehl, A. U. B. Wolter, B. B¨ uchner, and R. Klingeler, Nematicity and structure in LaFe 1 -x Co AsO, J. Magn. Magn. Mater. x 482 , 50 (2019).
- [49] F. Scaravaggi, A. U. B. Wolter et al. , Revisiting the phase diagram of LaFe 1 -x Co AsO x single crystals by thermodynamic methods , unpublished.
- [50] J.-H. Chu, J. G. Analytis, K. D. Greve, P. L. McMahon, Z. Islam, Y. Yamamoto, and I. R. Fisher, In-plane resistivity anisotropy in an underdoped iron arsenide superconductor , Science 329 , 824 (2010).
- [51] M. A. Tanatar, E. C. Blomberg, A. Kreyssig, M. G. Kim, N. Ni, A. Thaler, S. L. Bud ko, ′ P. C. Canfield, A. I. Goldman, I. I. Mazin, and R. Prozorov, Uniaxial-strain mechanical detwinning of CaFe As 2 2 and BaFe As 2 2 crystals: optical and transport study , Phys. Rev. B 81 , 184508 (2010).
- [52] H.-H. Kuo and I. R. Fisher, Effect of disorder on the resistivity anisotropy near the electronic nematic phase transition in pure and electron-doped BaFe As , 2 2 Phys. Rev. Lett. 112 , 227001 (2014).
- [53] G. Lang, H.-J. Grafe, D. Paar, F. Hammerath, K. Manthey, G. Behr, J. Werner, and B. B¨ uchner, Nanoscale electronic order in iron pnictides , Phys. Rev. Lett. 104 , 097001 (2010).
- [54] G. Lang, L. Veyrat, U. Gr¨fe, F. Hammerath, D. Paar, G. a Behr, S. Wurmehl, and H.-J. Grafe, Spatial competition of the ground states in 1111 iron pnictides , Phys. Rev. B 94 , 014514 (2016).
- [55] see Supplemental Materials.
- [56] C. Wuttke, C. Hess et al. , Nernst and elastoresistivity measurements of BaFe 2 -x Rh As , unpublished. x 2
- [57] S. Licciardello, J. Buhot, J. Lu, J. Ayres, S. Kasahara, Y. Matsuda, T. Shibauchi, and N. E. Hussey, Electrical resistivity across a nematic quantum critical point , Nature

567 , 213 (2019).

- [58] K. Matsuura, Y. Mizukami, Y. Arai, Y. Sugimura, N. Maejima, A. Machida, T. Watanuki, T. Fukuda, T. Yajima, Z. Hiroi, K. Y. Yip, Y. C. Chan, Q. Niu, S. Hosoi, K. Ishida, K. Mukasa, S. Kasahara, J.-G. Cheng, S. K. Goh, Y. Matsuda, Y. Uwatoko, and T. Shibauchi, Maximizing T c by tuning nematicity and magnetism in FeSe 1 -x S x superconductors , Nat. Commun. 8 , 1143 (2017).
- [59] C. G. Wang, Z. Li, J. Yang, L. Y. Xing, G. Y. Dai, X.
- C. Wang, C. Q. Jin, R. Zhou, and G.-Q. Zheng, Electron mass enhancement near a nematic quantum critical point in NaFe 1 -x Co As, Phys. Rev. Lett. x 121 , 167004 (2018).
- [60] W. Y. Wang, Y. Song, C. D. Cao, K.-F. Tseng, T. Keller, Y. Li, L. W. Harriger, W. Tian, S. X. Chi, R. Yu, A. H. Nevidomskyy, and P. C. Dai, Local orthorhombic lattice distortions in the paramagnetic tetragonal phase of superconducting NaFe 1 -x Ni As, Nat. Commun. x 9 , 3128 (2018).
Xiaochen Hong Federico Caglieris Rhea Kappenberger Sabine Wurmehl Saicharan Aswartham Francesco Scaravaggi Piotr Lepucki Anja U. B. Wolter Hans-Joachim Grafe Bernd Büchner Christian Hess
Evolution of the nematic susceptibility in LaFe$_{1-x}$Co$_x$AsO
The identification of electronic nematicity across series of iron-based
superconductors raises the question of its relationship with superconductivity
and other ordered states. Here, we report a systematic elastoresistivity study
on LaFe$_{1-x}$Co$_x$AsO single crystals, which have well separated structural
and magnetic transition lines. All crystals show Curie-Weiss-like nematic
susceptibility in the tetragonal phase. The extracted nematic temperature is
monotonically suppressed upon cobalt doping, and changes sign around the
optimal doping level, indicating a possible nematic quantum critical point
beneath the superconducting dome. The amplitude of nematic susceptibility shows
a peculiar double-peak feature. This could be explained by a combined effect of
different contributions to the nematic susceptibility, which are amplified at
separated doping levels of LaFe$_{1-x}$Co$_x$AsO.