Nematic quantum critical point without magnetism in FeSe{}_{1-x}S{}_{x} superconductors

Nematic quantum critical point without magnetism in FeSeS superconductors

S. Hosoi    K. Matsuura    K. Ishida    Hao Wang    Y. Mizukami Department of Advanced Materials Science, University of Tokyo, Kashiwa, Chiba 277-8561, Japan    T. Watashige    S. Kasahara    Y. Matsuda Department of Physics, Kyoto University, Sakyo-ku, Kyoto 606-8502, Japan    T. Shibauchi Department of Advanced Materials Science, University of Tokyo, Kashiwa, Chiba 277-8561, Japan
July 12, 2019

The importance of antiferromagnetic fluctuations are widely acknowledged in most unconventional superconductors Moriya2003 (); Scalapino2012 (); Hirschfeld2011 (). In addition, cuprates and iron pnictides often exhibit unidirectional (nematic) electronic correlations Keimer2015 (); Fradkin2015 (); Fernandes2014 (), including stripe and orbital orders, whose fluctuations may also play a key role for electron pairing Maier2014 (); Metlitski2015 (); Lederer2015 (); Kontani2010 (). However, these nematic correlations are intertwined with antiferromagnetic or charge orders, preventing us to identify the essential role of nematic fluctuations. This calls for new materials having only nematicity without competing or coexisting orders. Here we report systematic elastoresistance measurements in FeSeS superconductors, which, unlike other iron-based families, exhibit an electronic nematic order without accompanying antiferromagnetic order. We find that the nematic transition temperature decreases with sulphur content , whereas the nematic fluctuations are strongly enhanced. Near , the nematic susceptibility diverges towards absolute zero, revealing a nematic quantum critical point. This highlights FeSeS as a unique nonmagnetic system suitable for studying the impact of nematicity on superconductivity.

In underdoped cuprate superconductors, unidirectional electronic correlations (stripe correlations) appear in the pseudogap state, whose relation with superconductivity is a centre of debate. It has become more complicated after the charge density wave (CDW) order has been observed in a portion of this pseudogap region of the phase diagram Keimer2015 (). In iron pnictides, the tetragonal-to-orthorhombic structural transition always precedes or coincides with the antiferromagnetic (AFM) transition Hirschfeld2011 (). Below the structural transition temperature , electronic nematicity that represents a large electronic anisotropy breaking the rotational symmetry, is observed Fernandes2014 (), which may have a similar aspect with the stripe correlations in underdoped cuprates. In both cases, however, the nematicity is largely coexisting and intertwined with other CDW and AFM orders. Large nematic fluctuations have been experimentally observed in BaFeAs systems above , and these nematic fluctuations are strikingly enhanced with approaching the end point of the structural transition Yoshizawa2012 (); Kuo2015 (). This suggests the presence of a nematic quantum critical point (QCP), but in this case the nematic QCP coincides with, or locates very close to the AFM QCP, where antiferromagnetic fluctuations are also enhanced Shibauchi2014 (). This raises a fundamental question as to which fluctuations are the main driving force of the Cooper pairing in this system.

Recently there is growing evidence that the iron-chalcogenide FeSe displays remarkable properties. The superconducting transition temperature of  K at ambient pressure is largely enhanced up to 38 K under pressure Sun2015 (). Further enhancement of has been reported in monolayer FeSe Wang2012 (). What is important in this system is that despite its high nematic transition temperature of  K Boehmer2013 (), no magnetic order occurs down to  K. Above , large nematic fluctuations are reported but, unlike iron pnictides, no sizable AFM fluctuations are observed Baek2015 (). The enhanced nematic fluctuations most likely have an orbital origin, because a momentum-dependent orbital polarization has been found in angle-resolved photoemission spectroscopy in the nematic phase below Suzuki2015 (). More recently, it has been shown that the high-energy spin excitations measured by the inelastic neutron scattering do not show significant temperature dependence above Shamoto2015 (), which indicates no direct correlations between dynamical spin susceptibility and the nematic susceptibility, inconsistent with the spin-nematic scenario Fernandes2014 (). Although these results point to that the nematic transition in FeSe is orbital driven, the most important issue, namely the presence of a nonmagnetic nematic QCP at which the nematic fluctuation diverges with  K, remains open. The presence or absence of such a QCP should provide pivotal information on the superconductivity of iron-based high- materials. It has been shown very recently that the applying pressure leads to the suppression of , but near the critical pressure at which vanishes, the magnetic ordering is induced Sun2015 (), preventing our understanding of the interplay between nematicity and superconductivity.

Recent advances in the crystal growth of iron chalcogenides by using the chemical vapour transport technique Boehmer2013 () allow us to grow high-quality single crystals of FeSeS. The high quality of the obtained crystals is evidenced by the observation of quantum oscillations in a wide range of (Kasahara, S. et al., and Coldea, A. I. et al., unpublished results). Here we report that the partial substitution of isoelectric S for the Se site Mizuguchi2008 (); Guo2014 (); Watson2015 (); Moore2015 () leads to the complete suppression of this nematic transition without inducing the AFM order. Moreover shows a maximum at in FeSeS, suggesting the presence of superconducting dome. Therefore FeSeS appears to be a key material to investigate a possible correlation between the nonmagnetic nematic fluctuations and the superconductivity.

The temperature dependence of the in-plane resistivity in these high-quality crystals shows a kink at the nematic transition at , which is gradually decreased to lower temperatures (Fig. 1a). At high concentrations , the kink anomaly of disappears, which is more clearly seen in the temperature derivative (Fig. 1b). In contrast to the physical pressure case where another clear anomaly at a temperature different from is reported indicating the appearance of a magnetic transition, we find no such anomaly in the entire range of this study, which implies that the ground states of this system are nonmagnetic up to at least .

An elegant way to evaluate experimentally the nematic fluctuations has been developed by the Stanford group, that is based on the elastoresistance measurements by using a piezoelectric device Chu2012 (); Kuo2013 (). We utilize this technique in our FeSeS to extract the temperature dependence of the so-called nematic susceptibility . Here the change in the resistance induced by the strain that can be controlled by the voltage applied to the piezoelectric device is measured as a function of strain (Fig. 2a). This quantity is proportional to the nematic order parameter, and thus its fluctuations can be quantified by the nematic susceptibility , which is given by the slope of the versus curve Chu2012 (). To cancel the effect of anisotropic strain induced in the piezoelectric device (the Poisson effect) Kuo2013 (), we measured the resistance changes and for two current directions orthogonal each other in FeSe () crystals (Figs. 2a, b). The nematic susceptibility data obtained from and are essentially identical to those from alone (Fig. 3a). Therefore, we measure only for substituted samples to avoid the effect of small difference of the composition in different crystals from the same batch. The sign of the nematic susceptibility above is positive, which corresponds to the positive value of where and represent the orthorhombic crystal axes (). This positive sign is consistent with the previous results in FeSe Watson2015 (); Tanatar2015 () and FeTeSe Kuo2015 (), but is opposite to the Co-doped BaFeAs case Chu2012 (). This sign difference may be related to the difference of Fermi surface structure, which determines the anisotropic scattering of quasiparticles in the nematic state Blomberg2013 ().

The temperature dependence of the obtained nematic susceptibility shows a systematic trend as a function of the sulphur content (Figs. 3a-e). In the tetragonal state above , can be fitted to the Curie-Weiss type temperature dependence given by


where and are constants, and is the Weiss temperature, which represents the electronic nematic transition temperature without the lattice coupling. The observed changes in the resistance is given by , the last two terms on the right hand side represent geometric factors related to changes in the length and cross-sectional area Kuo2013 (). These terms contribute to the temperature-independent term in Eq. (1), which is dominant in simple metals. In the present FeSeS system, however, has strong temperature dependence, indicating that nematic fluctuations are very large. With a proper choice of , follows nearly -linear dependence in a wide temperature range of the tetragonal phase, consistent with the Curie-Weiss law (Eq. (1)).

The obtained phase diagram based on our single-crystal study is given by Fig. 4. The nematic transition temperature monotonically decreases and vanishes at . The Weiss temperature crosses the zero line and changes sign at , where the magnitude of is strongly enhanced. These results provide clear evidence that the nematic QCP is present near this composition , where nematic fluctuations are diverging toward  K. We stress that the present measurements are performed under identical conditions in a series of samples with similar dimensions and thus the obtained systematic results represent intrinsic properties of this system. We also note that the recent quantum oscillations and low-field magnetoresistance experiments indicate that the electronic structure of high-concentration samples () is quite different from that of low- crystals (Kasahara, S. et al., and Coldea, A. I. et al., unpublished results), which is indicative of distinctly different ground states in the two regimes, consistent with the quantum phase transition found in the present study.

Near the QCP, the data exhibit some deviations from the Curie-Weiss temperature dependence (grey regions in Figs. 3c-e, h-j). Such deviations have also been frequently found in iron-based superconductors when becomes low Kuo2015 (). The origin of these deviations from the Curie-Weiss law has not been fully understood, but one should consider possible effects of impurity scattering that may be modified by quantum criticality. In the theory of magnetic quantum critical point, various physics quantities show non-Fermi liquid behaviour, and magnetic susceptibility may deviate from the Curie-Weiss law Moriya2003 (). Magnetic susceptibility directly couples to the correlation length, which depends on dimension and dynamical exponent of the system according to the Hertz-Millis theory Lohneysen2007 (). In analogy to this, the nematic correlation length in the presence of disorder would have some characteristic temperature dependence different from the Fermi-liquid behaviour, and then the deviations from Curie-Weiss law may be expected in the nematic susceptibility. When the effects of disorder is strong, it has been discussed that the so-called Griffith singularity may occur, which modifies the critical behaviour Kuo2015 (). In any case, these disorder effects should weaken the divergent behaviour especially near the QCP, which is consistent with our results.

The present results provide evidence that a nonmagnetic nematic QCP lies in the superconducting phase, which has the same topology as the phase diagrams of iron pnictides Shibauchi2014 () and heavy fermions Lohneysen2007 () in which AFM QCP locates inside the superconducting dome. However, this new kind of QCP has several different aspects from the well-studied AFM QCP. First, the resistivity data (Fig. 1a) shows non- dependence for all the sample we measured, and we find no significant enhancement in the residual resistivity near the nematic QCP, which appears to be different from the expectations in the magnetic QCP case Lohneysen2007 (). Second, in iron pnictides shows the maximum at around the AFM QCP Shibauchi2014 (), wheres in the present system shows a maximum deep inside the nematic ordered phase, not at the nematic QCP. This suggests that the critical nematic fluctuations do not have simple correlations with the enhancement, although some theoretical calculations propose a peak at a nematic QCP Lederer2015 (). Possible origins of this unexpected result are related to the fact that the nematic order is most likely a ferro-type order with the wave vector , which is different from AFM order with a finite . In a ferro-type order, the Fermi surface change is not as dramatic as in the antiferro-type order in which band-folding occurs, and thus the competition between superconductivity and ferro-nematic order may be much more modest. Indeed, the unusual lack of coupling between superconductivity and orthorhombic distortion in nonmagnetic FeSe has been reported by the thermal expansion measurements Boehmer2013 (). Another factor which may be relevant here is the enhanced quasiparticle damping effect near a ferro-type QCP, that can suppress . In fact, such a scenario has been discussed in ferromagnetic superconductors, such as UGe, in which the peak of the dome is not located at a ferromagnetic end point Lohneysen2007 (). However, the resistivity curves do not show a dramatic change near the present nematic QCP, and thus further studies are necessary to clarify this point.

It should be noted that in cuprate superconductors the importance of nematic electronic correlations has been suggested as well Lederer2015 (); Maier2014 (), and the quantum criticality of the CDW or the pseudogap phase has been discussed Keimer2015 (). In YBaCuO Ramshaw2015 () and YBaCuO Putzke2016 (), however, opposite correlations between and the quasiparticle mass enhancement associated with the putative quantum phase transition have been reported, causing controversy on the relationship between quantum fluctuations and superconductivity in cuprates. To solve this issue it would be important to separate the effects of CDW order from the nematic correlations in the pseudogap state.

The present discovery of the nematic QCP and enhanced ferro-type nematic fluctuations without coexisting or competing magnetic order in FeSeS may provide a unique avenue to study the unconventional superconductivity mediated by exotic mechanisms other than spin fluctuations.


Single crystals of FeSeS were grown by the chemical vapour transport technique. A mixture of Fe, Se and S powders together with KCl and AlCl powders was sealed in an evacuated SiO tube. The ampoule was heated to C on one end while the other end was kept at C. In these conditions, we were able to cover a much wider range than previous single crystal studies Watson2015 (); Moore2015 (), which is important to completely suppress the nematic transition. The actual sulphur composition of is determined by the energy dispersive X-ray spectroscopy, and is found to be about 80% of the nominal S content (Fig. 1, inset). The experimental setup for the measurement of nematic susceptibility is shown in the inset of Fig. 2 Kuo2013 (). The samples are cut into bar-shape along the orthogonal crystal axis (tetragonal ) and are glued on a piezoelectric stack. After mounted on the piezoelectric stack, samples are cleaved so that the thickness becomes about m and the strain is sufficiently transmitted to them. The strain is induced by applying the voltage to piezoelectric stack and measured by strain gauges glued to the surface of piezoelectric stack. We simultaneously measure the changes in the resistance when varying the strain and this induced anisotropy directly couples to the electronic nematic order parameter.


We thank J.-H. Chu, A. I. Coldea, R. M. Fernandes, I. R. Fisher, C. W. Hicks, H. Kontani, S. Lederer, A. H. Nevidomskyy, T. Takimoto, M. D. Watson, and T. Wolf for helpful comments and discussions. This work was supported by Grants-in-Aid for Scientific Research (KAKENHI) from Japan Society for the Promotion of Science (JSPS), and by the ‘Topological Material Science’ Grant-in-Aid for Scientific Research on Innovative Areas from the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan.


  • (1) Moriya, T. & Ueda, K. Antiferromagnetic spin fluctuation and superconductivity. Rep. Prog. Phys. 66 1299-1341 (2003).
  • (2) Scalapino, D. J. A common thread: the pairing interaction for the unconventional superconductors. Rev. Mod. Phys. 84, 1383 (2012).
  • (3) Hirschfeld, P. J., Korshunov, M. M., & Mazin, I. I. Gap symmetry and structure of Fe-based superconductors. Rep. Prog. Phys. 74, 124508 (2011).
  • (4) Keimer, B., Kivelson, S. A., Norman, M. R., Uchida, S., & Zaanen, J. From quantum matter to high-temperature superconductivity in copper oxides. Nature 518, 179-186 (2015).
  • (5) Fradkin, E., Kivelson, S. A., & Tranquada, J. M. Theory of intertwined orders in high temperature superconductors. Rev. Mod. Phys. 87, 457-482 (2015).
  • (6) Fernandes, R. M., Chubukov, A. V., & Schmalian, J. What drives nematic order in iron-based superconductors? Nat. Phys. 10, 97-104 (2014).
  • (7) Maier, T. A. & Scalapino, D. J. Pairing interaction near a nematic quantum critical point of a three-band CuO model. Phys. Rev. B 90, 174510 (2014).
  • (8) Metlitski, M. A., Mross, D. F., Sachdev, S., & Senthil, T. Cooper pairing in non-Fermi liquids. Phys. Rev. B 91, 115111 (2015).
  • (9) Lederer, S., Schattner, Y. Berg, E., & Kivelson, S. A. Enhancement of superconductivity near a nematic quantum critical point. Phys. Rev. Lett. 114, 097001 (2015).
  • (10) Kontani, H. & Onari, S. Orbital-fluctuation-mediated superconductivity in iron pnictides: analysis of the five-orbital Hubbard-Holstein model. Phys. Rev. Lett. 104, 157001 (2012).
  • (11) Kuo, H.-H., Chu, J.-H., Kivelson, S. A., & Fisher, I. R. Ubiquitous signatures of nematic quantum criticality in optimally doped Fe-based superconductors. arXiv:1503.00402 (2015).
  • (12) Yoshizawa, M. et al. Structural quantum criticality and superconductivity in iron-based superconductor Ba(FeCo)As. J. Phys. Soc. Jpn. 81, 024604 (2012).
  • (13) Shibauchi, T., Carrington, A., & Matsuda, Y. A quantum critical point lying beneath the superconducting dome in iron-pnictides. Annu. Rev. Condens. Matter Phys. 5, 113-135 (2014).
  • (14) Sun, J. P. et al. Dome-shaped magnetic order competing with high-temperature superconductivity at high pressures in FeSe. arXiv:1512.06951 (2015).
  • (15) Wang, Q.-Y. et al. Interface-induced high-temperature superconductivity in single unit-cell FeSe films on SrTiO. Chin. Phys. Lett. 29, 037402 (2012).
  • (16) Böhmer, A. E. et al. Lack of coupling between superconductivity and orthorhombic distortion in stoichiometric single-crystalline FeSe. Phys. Rev. B 87, 180505(R) (2013).
  • (17) Baek, S.-H. et al. Orbital-driven nematicity in FeSe. Nat. Mater. 14, 210-214 (2015).
  • (18) Suzuki, Y. et al. Momentum-dependent sign inversion of orbital order in superconducting FeSe. Phys. Rev. B 92, 205117 (2015).
  • (19) Shamoto, S. et al. Spin nematic susceptibility studied by inelastic neutron scattering in FeSe. arXiv:1511.04267 (2015).
  • (20) Mizuguchi, Y., Tomioka, F., Tsuda, S., Yamaguchi, T., & Takano, Y. Substitution effects on FeSe superconductor. J. Phys. Soc. Jpn. 78, 074712 (2009).
  • (21) Guo, J., Lei, H., Hayashi, F., & Hosono, H. Superconductivity and phase instability of NH-free Na-intercalated FeSeS. Nat. Commun. 5, 4756 (2014).
  • (22) Watson, M. D. et al. Suppression of orbital ordering by chemical pressure in FeSeS. Phys. Rev. B 92, 121108(R) (2015).
  • (23) Moore, S. A. et al. Evolution of the superconducting properties in FeSeS. Phys. Rev. B 92, 235113 (2015).
  • (24) Chu, J.-H., Kuo, H.-H., Analytis, J. G., & Fisher, I. R. Divergent nematic susceptibility in an iron arsenide superconductor. Science 337, 710-712 (2012).
  • (25) Kuo, H.-H., Shapiro, M. C., Riggs, S. C. & Fisher, I. R. Measurement of the elastoresistivity coefficients of the underdoped iron arsenide Ba(FeCo)As. Phys. Rev. B 88, 085113 (2013).
  • (26) Tanatar, M. A. et al. Origin of the resistivity anisotropy in the nematic phase of FeSe. arXiv:1511.04757 (2015).
  • (27) Blomberg, E. C. et al. Sign-reversal of the in-plane resistivity anisotropy in hole-doped iron pnictides. Nat. Commun. 4, 1914 (2013).
  • (28) v. Löhneysen, H., Rosch, A., Vojta, M., & Wölfle, P. Fermi-liquid instabilities at magnetic quantum phase transitions. Rev. Mod. Phys. 79, 1015-1075 (2007).
  • (29) Ramshaw, B. J. et al. Quasiparticle mass enhancement approaching optimal doping in a high- superconductor. Science 348, 317-320 (2015).
  • (30) Putzke, C. et al. Inverse correlation between quasiparticle mass and in a cuprate high- superconductor. Sci. Adv. 2, e1501657 (2016).
Figure 1: Temperature dependence of resistivity in FeSeS single crystals. a, In-plane resistivity as a function of temperature in several crystals with different . Each curve is shifted vertically for clarity. Sulphur content determined by energy dispersive X-ray spectroscopy is linearly proportional to the nominal composition (inset). b, Temperature dependence of . Each curve is shifted vertically so that all the data at 120 K merge. The arrows indicate the nematic (structural) transition temperature estimated from the inflection point, which is more accurately determined by the peak position in the temperature dependence of the nematic susceptibility (Figs. 3a-c).
Figure 2: Elastoresistance measurements by using the piezoelectric device. a, Relative change in the resistance of an FeSe crystal at several temperatures above as a function of strain induced by the attached piezoelectric stack. The current direction in the crystal is parallel to the strain direction. b, Similar data for another crystal in which the current direction in perpendicular to the strain direction. Inset is a schematic experimental set-up. The crystals (black bars) are attached on top of the piezoelectric stack (green) which expands (shrinks) along the violet arrows when positive (negative) bias is applied. The strain is measured by a strain gauge (red) attached underneath the device.
Figure 3: Divergent nematic susceptibility in FeSeS single crystals. a-e, Temperature dependence of the nematic susceptibility for (a), 0,08 (b), 0.12 (c), 0.17(d), and 0.21 (e). In FeSe (), the nematic susceptibility data extracted by using both (red filled circles) and (violet open circles) show consistent results. For other only data has been used (b-e). f-j, Temperature dependence of inverse of for (f), 0,08 (g), 0.12 (h), 0.17(i), and 0.21 (j). Here is a temperature independent constant, coming from the geometrical change that leads to a resistance change under strain. Solid lines are the fits of to the Curie-Weiss law , and the obtained values of the Weiss temperature are indicated in the lower panels. The nematic transition is defined by the peak in , which is quantitatively consistent with the kink anomaly in curves (Fig. 1). In the temperature range highlighted in grey, the nematic susceptibility data above shows some deviations from the Curie-Weiss law.
Figure 4: Phase diagram and quantum criticality in FeSeS. Temperature dependence of the nematic transition (, green diamonds) and the superconducting transition temperature (, orange circles) determined by the zero resistivity criteria. The Weiss temperature obtained by the Curie-Weiss analysis of the nematic susceptibility data (Fig. 3) is also plotted (, red hexagons). The magnitude of in the tetragonal phase is superimposed in the phase diagram by a colour contour (see the colour bar for the scale). The lines are the guides for the eyes.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

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 description