Dome-shaped magnetic order competing with high-temperature superconductivity at high pressures in FeSe

Abstract

The coexistence and competition between superconductivity and electronic orders, such as spin or charge density waves, have been a central issue in high transition-temperature (T_c ) superconductors. Unlike other iron-based superconductors, FeSe exhibits nematic ordering without magnetism whose relationship with its superconductivity remains unclear. Moreover, a pressure-induced fourfold increase of T_c has been reported, which poses a profound mystery. Here we report high-pressure magnetotransport measurements in FeSe up to ∼15 GPa, which uncover the dome shape of magnetic phase superseding the nematic order. Above ∼6 GPa the sudden enhancement of superconductivity (T_c ≤38.3 K) accompanies a suppression of magnetic order, demonstrating their competing nature with very similar energy scales. Above the magnetic dome, we find anomalous transport properties suggesting a possible pseudogap formation, whereas linear-in-temperature resistivity is observed in the normal states of the high-T_c phase above 6 GPa. The obtained phase diagram highlights unique features of FeSe among iron-based superconductors, but bears some resemblance to that of high-T_c cuprates.

ARTICLE

Received 19 Feb 2016 | Accepted 5 Jun 2016 | Published 19 Jul 2016

J.P. Sun1,*, K. Matsuura2,*, G.Z. Ye1,3,*, Y. Mizukami2, M. Shimozawa4, K. Matsubayashi5, M. Yamashita4,T. Watashige6, S. Kasahara6, Y. Matsuda6, J.-Q. Yan7,8, B.C. Sales7, Y. Uwatoko4, J.-G. Cheng1 & T. Shibauchi2

The coexistence and competition between superconductivity and electronic orders, such as spin or charge density waves, have been a central issue in high transition-temperature (Tc)

superconductors. Unlike other iron-based superconductors, FeSe exhibits nematic ordering without magnetism whose relationship with its superconductivity remains unclear. Moreover, a pressure-induced fourfold increase of Tc has been reported, which poses a profound mystery. Here we report high-pressure magnetotransport measurements in FeSe up to

B15 GPa, which uncover the dome shape of magnetic phase superseding the nematic order. Above B6 GPa the sudden enhancement of superconductivity (Tcr38.3 K) accompanies a suppression of magnetic order, demonstrating their competing nature with very similar energy scales. Above the magnetic dome, we nd anomalous transport properties suggesting a possible pseudogap formation, whereas linear-in-temperature resistivity is observed in the normal states of the high-Tc phase above 6 GPa. The obtained phase diagram highlights unique features of FeSe among iron-based superconductors, but bears some resemblance to that of high-Tc cuprates.

1 Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China. 2 Department of Advanced Materials Science, University of Tokyo, Kashiwa, Chiba 277-8561, Japan. 3 School of Physical Science and Technology, Yunnan University, Kunming 650091, China. 4 The Institute for Solid State Physics, The University of Tokyo, Kashiwa, Chiba 277-8581, Japan. 5 Department of Engineering Science, The University of Electro-Communications, Chofu, Tokyo 182-8585, Japan. 6 Department of Physics, Kyoto University, Sakyo-ku, Kyoto 606-8502, Japan.

7 Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA. 8 Department of Materials Science and Engineering, University of Tennessee, Knoxville, Tennessee 37996, USA. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to J.-Q.Y. (email: mailto:[email protected]

Web End [email protected] ) or to J.-G.C. (email: mailto:[email protected]

Web End [email protected] ) or to T.S. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 7:12146 | DOI: 10.1038/ncomms12146 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

DOI: 10.1038/ncomms12146 OPEN

Dome-shaped magnetic order competing with high-temperature superconductivity at high pressures in FeSe

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12146

In most iron-based superconductors, superconductivity emerges on the verge of a long-range antiferromagnetically ordered state1, which is a common feature to many

unconventional superconductors2,3 including the cuprates and heavy-fermion materials. It has been shown that the antiferro-magnetic order in the iron-pnictide materials accompanies or follows the tetragonal-to-orthorhombic structural transition at Ts.

In striking contrast, the structurally simplest FeSe exhibits a high TsE90 K but no magnetic order appears at lower temperatures47, and still its ground state is an unconventional superconducting state with TcE9 K (refs 810). This material is also intriguing in that in the form of one-unit-cell-thick lms a very high Tc (up to 109 K) has been reported recently1113, which is likely associated with a carrier-doping effect14,15 from the substrate. In bulk FeSe, a signicant electronic anisotropy is found below Ts in the nonmagnetic orthorhombic phase, which is often called a nematic state9,10,16,17. In the nematic phase, very small Fermi surfaces with strong deviations from the rst-principles calculations have been observed10,18,19, and the occurrence of superconductivity with such small Fermi energies is quite unusual, implying that the system is deep in the crossover regime between the weak-coupling BardeenCooperSchrieffer and strong-coupling BoseEinsteincondensate limits10.

In addition to these distinct electronic characteristics of FeSe, remarkable properties have been reported under high pressure2027. First of all, the initial study on powder samples has shown that the relatively low TcE9 K at ambient pressure can be enhanced by more than fourfold to B37 K under B8 GPa, pushing it into the class of high-Tc superconductors21. More recent studies under better hydrostatic pressure conditions revealed a complex temperaturepressure (TP) phase diagram featured by a suppression of Ts around 2 GPa, a sudden development of static magnetic order above B1 GPa (ref. 23), and an enhancement of Tc in a three-plateau process24, that is, TcB10(2) K for 02 GPa, TcB20(5) K for 35 GPa, and TcB35(5) K for 68 GPa. The rst jump of Tc from B10 to B20 K seems to coincide with the suppression of the nonmagnetic nematic state and the development of the long-range magnetic order at Tm evidenced by mSR measurements22. The observation that both Tc and Tm increase with pressure in the pressure range 12.5 GPa has been taken as evidence for the cooperative promotion of superconductivity by the static magnetic order22. Such a scenario, however, does not t into the general scope of iron-based superconductors, in which the optimal superconductivity is realized when the long-range magnetic order is close to collapse1,28. This issue remains unclear unless the fate of magnetic order at Tm under higher pressures is sorted out. Due to the technical limitations of probing small-moment magnetic order above 3 GPa, this task only becomes possible very recently when a clear signature at Tm is visible in the resistivity23,25,26 of high-quality FeSe single crystals29. We also note that more recently the pressure-induced magnetic order in these single crystals has been conrmed below Tm by Mssbauer30 and nuclear magnetic resonance (NMR) measurements31.

Here by performing the high-pressure resistivity r(T) measurements up to B15 GPa on high-quality single crystals, we construct for bulk FeSe the most comprehensive TP phase diagram mapping out the explicit evolutions with pressure of Ts,

Tc and Tm. We uncover a previously unknown dome-shaped Tm(P), having two end points situated on the boundaries separating the three plateaus of Tc(P). Our results thus provide compelling evidence linking intimately the sudden enhancement of Tc to 38 K to the suppression of long-rang magnetic order. This highlights a competing nature between magnetic order and high-Tc superconductivity in the phase diagram of FeSe, which is a key material among the iron-based superconductors.

ResultsLow-pressure region. The tetragonal-orthorhombic structure transition at TsE90 K for bulk FeSe (blue square in Fig. 1) is manifested as a slight upturn in resistivity, which can be taken as a signature to track down the evolution of Ts with pressure. Our resistivity r(T) data measured with a self-clamped piston cylinder cell (PCC) up to B1.9 GPa are shown in Fig. 2a. As can be seen, Ts is suppressed progressively to below 50 K at B1.5 GPa, above which the anomaly at Ts becomes poorly dened. Meanwhile, a second anomaly manifested as a more profound upturn in r(T) emerges at TmB20 K and moves up steadily with pressure. In light of the recent high-pressure mSR, Mssbauer, and NMR studies22,30,31, this anomaly at Tm corresponds to the development of long-range magnetic order. We also note that in this magetically ordered state below Tm, the

orthorhombic structure similar to the one (space group Cmma) in the nematic phase has been reported recently30,31. Ts and Tm seem to cross around B2 GPa. In this pressure range, the superconducting transition temperature Tc (dened as the zero-resistivity temperature) rst increases and then decreases slightly before rising again. This features a small dome-shaped Tc(P)

peaked at B1.2 GPa (Fig. 1), which roughly coincides with the pressure where the long-range magnetic order at Tm starts to emerge. These results in this relatively low-pressure range are in general consistent with those reported previously23,25,26.

High-pressure region. To further track down the evolution of Tm, we turn to r(T) measurements in cubic anvil cells (CACs)

that can maintain a quite good hydrostaticity up to B15 GPa (refs 3234). Figure 2bd displays the r(T) data measured in two self-clamped CACs and one constant-loading CAC (see Methods for experimental details). In line with the results of PCC in Fig. 2a, the sudden upturn is clearly visible at Tm in both measurements in the pressure range up to B2.5 GPa (Figs 2b,d and 3a,b), above which the upturn anomaly disappears and instead a kink appears in r(T) followed by a gradual drop before reaching the superconducting transition (Figs 2b,d and 3c,d). Our

100

80

60

40

20

0 0 2 4 6 8 10P (GPa)

T(K)

Piston cylinder cell (PCC)

Clamp-type cubic anvil cell (CAC)

Tc

Constant loading CAC

Ts

Tm

Nematic

SDW

SC

Figure 1 | Temperaturepressure phase diagram of bulk FeSe. The structural (Ts, blue), magnetic (Tm, green), and superconducting transition temperatures (Tc, red and black) as a function of hydrostatic pressure in high-quality single crystals determined by anomalies in resistivity r(T)

measured in the PCC (open circles), clamp-type CAC (closed circles), and constant-loading type CAC (closed squares). Tc values determined from the ac magnetic susceptibility (w(T)) measurements in the clamp-type CAC are also shown (solid triangles). The magnetic phase is most likely a spin density wave (SDW) phase. Colour shades for the nematic, SDW, and superconducting (SC) states are guides to the eyes.

2 NATURE COMMUNICATIONS | 7:12146 | DOI: 10.1038/ncomms12146 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12146 ARTICLE

300

200

150

100

50

200

150 100

100

200

300

400

500

600

a b c d

150

100

100

200

300

400

500

T

T

0 GPa 0 GPa

250

0.2 GPa

0.43 GPa

0.67 GPa0.9 GPa1.2 GPa1.38 GPa1.51 GPa

1.63 GPa1.72 GPa1.83 GPa

T

0 GPa

1.8 GPa

2.3 GPa

2.8 GPa

3.3 GPa

3.8 GPa4.8 GPa5.8 GPa6.3 GPa

6.8 GPa

7.8 GPa8.8 GPa

15 GPa

0 0

300

250

200

150

100

50

0

13.5 GPa

12 GPa

8 GPa

10 GPa

T

T

2 GPa

(cm)

T

T

4 GPa

6 GPa

8 GPa

0

50

100

1.89 GPa

150

T

T

T

0

0 50 100 0 50 100 100 200 300 0 50 100 150 200 T (K) T (K) T (K) T (K)

Figure 2 | Temperature dependence of resistivity in FeSe single crystals under high pressure. (a) r(T) curves below 100 K at different pressures up to B1.9 GPa measured in the self-clamped PCC. (b) Data up to B8.8 GPa measured in the self-clamped CAC. (c) Data up to B15 GPa measured in the smaller self-clamped CAC. (d) r(T) curves below 200 K up to 8 GPa measured in the constant-loading CAC. Except for (c) the data are vertically shifted for clarity. The resistive anomalies at transition temperatures Ts, Tm, and Tc are indicated by the arrows.

self-clamped CAC enables us to measure r(T) under different magnetic elds32,33, which clearly resolves that this anomaly is of the same origin as the magnetic order at lower pressures. As seen in Fig. 4b, on the application of magnetic elds the resistivity kink at 2.8 GPa restores to the upturn anomaly and this anomaly shows little eld dependence as that observed at 1.8 GPa (Fig. 4a). On further increasing pressures, Tm moves up gradually and reaches about 45 K at 4.8 GPa, where the resistivity kink becomes sharper (Fig. 2b). In striking contrast with the progressive enhancement of Tm in this pressure range of 25 GPa, Tc remains nearly unchanged at the second plateau of B20 K.

The situation changes markedly when further increasing pressure above B5 GPa. As shown in Figs 2b and 3d, r(T) at5.8 GPa exhibits an abrupt drop at 38.5 K before reaching zero resistivity at Tc 27.5 K, and then at 6.3 GPa the abrupt drop

develops into a very sharp superconducting transition at Tc 38.3

K with a transition width of only 0.4 K (Figs 2b and 3e). Thus, Tc is nearly doubled from B20 K at 4.8 GPa to 38.3 K at 6.3 GPa, and the abrupt drop at 38.5 K under 5.8 GPa corresponds to the onset of superconductivity rather than the magnetic order, which is supported by the large change under magnetic elds (Fig. 4d). Then, what is the fate of Tm? A closer inspection of the r(T) data at 5.8 GPa reveals an inection point at the temperature slightly above the abrupt drop. This feature can be well resolved from the temperature derivative of resistivity dr/dT in Fig. 3d. Here the small peak centred at 41 K corresponds to Tm, which is conrmed by the eld independence (Fig. 4d). At 6.3 GPa, such a magnetic anomaly is absent in the dr/dT(T) curve at zero eld, but becomes visible under 9-T eld when Tc is suppressed to lower temperatures, as shown in Figs 3e and 4i. On further increasing pressure, the superconducting transition remains very sharp and Tc moves down slowly to 33.2 K at 8.8 GPa. Due to the large upper critical eld, Tm cannot be dened at P46.3 GPa (Fig. 4jl). These above ndings are further conrmed by separate high-pressure resistivity measurements with a constant-loading CAC, as shown in Fig. 2d.

At higher pressures, the resistivity curves measured by using a small CAC (Fig. 2c) show a sudden change from metallic and superconducting behaviour to semiconducting and

non-superconducting one at around 12 GPa. Indeed, the crystal structure studies of FeSe polycrystals at high pressures35,36 have shown that high pressures above B10 GPa will stabilize the collapsed three-dimensional orthorhombic Pbnm structure, which is non-superconducting. Therefore, we will focus our attention to the phase diagram at Po12 GPa, where the crystal structure symmetries should be the same as the ones in the ambient pressure.

The superconductivity is also checked by the ac susceptibility measurements under pressure using the self-clamped CAC (Fig. 5). The onset temperature Twc of the diamagnetic signal is quantitatively consistent with the zero-resistivity Tc at each pressure (Fig. 1). Moreover, the magnitude of the diamagnetic signals at lowest temperatures does not show signicant pressure dependence, suggesting the bulk superconducting nature in a wide pressure range covering the nematic, magnetic and nonmagnetic phases.

DiscussionWhen the obtained Tm(P) and Tc(P) data are plotted in Fig. 1, we uncover a previously unknown dome-shaped magnetic order with two ends situated near the boundaries separating the three-plateau Tc(P). This observation suggests that the superconductivity in FeSe is intimately correlated with the magnetism. Despite the absence of static magnetic order within the nematic phase at ambient pressure57, strong spin uctuations with an in-plane wave vector q (p, 0) have been identied recently37. Recent

theoretical studies38 have attributed the absence of static magnetic order to unusual magnetic frustration among multiple competing magnetic orders with q (p, x) with 0rxrp/2. According to the

calculations, this magnetic frustration can be removed under pressure, and gives way to a long-range magnetic order with q (p, 0). Indeed, the stabilization of static magnetic order at Tm

is concomitant with the destruction of nematic order as seen in Fig. 1, thus highlighting the competing nature between nematic and magnetic order. The qa0 nature of the pressure-induced magnetic phase is supported by the Fermi surface reconstruction recently reported from quantum oscillations25, and also

NATURE COMMUNICATIONS | 7:12146 | DOI: 10.1038/ncomms12146 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12146

a b

80

60

40

20

100

80

60

40

20

010 20 30 40 50 30 40 50

60

20

10

5

5

10

0

15

10

5

5

10

0

1.8 GPa

2.3 GPa

15

Tcpeak

60

(cm)

(cm)

(cm)

40

Tm

Tc

0

0 10 20 30 40 50 60 10 20 30 40 50 60

15

10

5

0 0

20

d /dT (cm K1 )d /dT (cm K1 ) d /d T ( cm K 1 )

20

c d

20

15

10

5

0

80

3.3 GPa

5.8 GPa

60

40

e f

60

40

20

0

20

15

10

5

0

60

6.3 GPa

9 T 0 T

40

20

20

0

8.8 GPa

30 40 50

10

8

6

4

2

0

40

60 80 100

T (K) T (K)

Figure 3 | Determination of magnetic and superconducting transition temperatures from the resistivity curves under high pressure. (af) Temperature dependence of resistivity r(T) (left axis) and temperature derivative of resistivity dr/dT(T) (right axis) at low temperatures. The temperature at the maximum slope of superconducting transition is dened as Tpeakc and the zero resistivity temperature as Tc. The magnetic transition temperature Tm is determined by the dip or peak in dr/dT(T). At 6.3 GPa, the data at 9 T is also shown where Tm anomaly appears when the superconducting transition is suppressed to a lower temperature (e). At 8.8 GPa, the normal-state resistivity can be tted to a T-linear dependence (dashed line in (f)).

consistent with the eld independence of Tm found in the present study (Fig. 4ad). Thus the pressure-induced order that supersedes the nematic order above B2 GPa is most likely a spin density wave (SDW), although the previous neutron scattering experiment was unsuccessful due to the small magnetic moment22. The most recent NMR study as well as the orthorhombicity found below Tm also support the stripe-type

SDW order with q (p, 0)30,31. Most importantly, our results

reveal unambiguously that the sudden jump of Tc around 6 GPa comes concomitantly with the suppression of magnetic order, demonstrating another competition between the SDW order and high-Tc superconductivity at high pressures. The peak temperature (E45 K) of the Tm(P) dome is quite close to the maximum Tc( 38.3 K), indicating that these competing orders

have very similar energy scales. This is different from other iron-based superconductors where the antiferromagnetic transition temperatures in the mother materials are signicantly higher than the maximum Tc values of the derived superconducting

phases1,28. We note that the non-monotonic dependence of Tc(P) is somewhat akin to the two-dome shape of Tc(x) recently found in LaFeAsO1 xHx (refs 1,39). However, the complete

dome shape of Tm(P) with two ends in a single-phase diagram of FeSe without changing of carrier balance is distinctly different from the case of LaFeAsO1 xHx, in which two

different antiferromagnetic orders exist near the opposite doping ends of the superconducting phase. Altogether with the presence of competing nematic order at low pressures, this phase diagram of FeSe thus exhibits unique features among iron-based superconductors. It should be noted that a phase diagram bearing some resemblance has been recently found in hole-doped cuprate superconductors, where the high-Tc superconducting dome is partially suppressed by the presence of dome-shaped charge density wave (CDW) order3.

It is also intriguing to point out that above the SDW dome, the resistivity curves in the pressure range of B36 GPa show anomalous sub-linear (convex) temperature dependence

4 NATURE COMMUNICATIONS | 7:12146 | DOI: 10.1038/ncomms12146 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12146 ARTICLE

a b c d

80 1.8 GPa 2.8 GPa 3.3 GPa 5.8 GPa

8.8 GPa

(cm)

(cm) 0H (T)

0H (T)

0 e f g h

Tc Tm

Tcpeak

0 10 20 30 40 0 10 20 30 40 0 10 20 30 40 0 10 20 30 40

0 10 20 30 40 0 10 20 30 40 0 10 20 30 40 0 10 20 30 40T (K) T (K) T (K) T (K)

T (K)

T (K) T (K) T (K)

i j k l

6.3 GPa

9 T7 T5 T3 T 1 T 0 T

6.8 GPa

7.8 GPa

0 m n o p

Figure 4 | Effects of magnetic elds on the magnetic and superconducting transitions under high pressure. (ad, il) r(T) curves at different magnetic elds applied parallel to the c axis. The magnetic transition temperature Tm is eld independent and marked by the green arrow (ad). At 6.3 GPa, the Tm anomaly is only visible at high elds (i). (eh, mp) H-T phase diagrams. The zero resistivity is attained below Tc and the Tpeakc line (determined by the peak in dr/dT(T)) is a lower bound of upper critical eld Hc2. The sharp superconducting transitions under magnetic elds for P46.3 GPa indicate narrowed regimes of the vortex liquid state.

0.0

0.2

0.4

0.6

0.80 10 20 30 40 50

T (K)

(Fig. 2b,d), which mimics the suppression of the quasiparticle scattering by the pseudogap formation above the CDW phase of underdoped cuprates. At higher pressures, this convex anomaly in the temperature dependence of resistivity becomes less pronounced. At 8.8 GPa, r(T) displays a perfect linear-in-T dependence in a wide temperature range of the normal state (Figs 2b and 3f), which is a hallmark of non-Fermi-liquid behaviour. The observation of the non-Fermi-liquid behaviour at high pressures where the magnetic order is vanishing is an indication of the presence of strong critical uctuations. Similar T-linear r(T) above Tc has been observed near a quantum critical point (QCP) at xE0.3 in the BaFe2(As1 xPx)2 system28. In that

system, the critical elds near the QCP show anomalous features pointing to an enhancement of the energy of superconducting vortices possibly due to a microscopic mixing of anti-ferromagnetism and superconductivity40. The analysis of the upper critical eld Hc2, which is related to the effective mass m* in a simple picture by m 2 / 1=Tc

dH=dT

TTc

Tc

1.8 GPa2.3 GPa2.8 GPa3.3 GPa3.8 GPa4.8 GPa5.8 GPa6.3 GPa6.8 GPa7.8 GPa8.8 GPa

4

Figure 5 | Temperature dependence of ac susceptibility v under high pressure. The data are taken by using the clamp-type CAC. The onset temperatures Twc of diamagnetic signals are marked by arrows.

j (ref. 41), has

been shown to be much less sensitive to the mass enhancement observed in other techniques near the QCP. In our similar

NATURE COMMUNICATIONS | 7:12146 | DOI: 10.1038/ncomms12146 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12146

Methods

Sample preparation and characterization. High-quality FeSe single crystals used in the present study have been grown by two different methods: the chemical vapour transport technique29 (Kyoto University) and the ux method43 (Oak Ridge National Laboratory). These two methods yield FeSe single crystals with similar quality in terms of the phase transition temperatures (Ts 8790 K and

Tc 8.59 K) and the residual resistivity ratio RRR B40. However, the ux method

produces much larger (up to B10 mm) and thicker (B1 mm) crystals and thus we cleaved and cut these samples to t in the pressure cells. The crystals are well characterized by means of X-ray diffraction, energy dispersive spectroscopy, magnetic and transport properties under ambient pressure.

High-pressure measurements. High-pressure resistivity r(T, P) measurements have been performed under hydrostatic pressures up to B2 GPa with a PCC and up to B15 GPa with CACs, respectively. For the crystals grown with the ux method, r(T, P, H) data under high pressures and magnetic elds were measured in the Institute of Physics, Chinese Academy of Sciences, by using the self-clamped type CACs and PCC. The data up to 9 GPa have been taken by a 4-mm CAC, in which we also measured ac susceptibility by a conventional mutual inductance technique. Higher-pressure data up to 15 GPa were taken by a 2.5-mm CAC. The high-pressure resistivity of crystals grown with the chemical vapour transfer method was measured in the Institute for Solid State Physics, University of Tokyo with a constant-loading type CAC, which can maintain a nearly constant pressure over the whole temperature range from 300 to 2 K. The pressure value inside the PCC was determined by monitoring the shift of superconducting transition of lead (Pb), while those of CACs were calibrated at room temperature by observing the characteristic transitions of bismuth (Bi). (It should be noted that for the self-clamped pressure cells, both PCC and CAC, the pressure value at room temperature is slightly different from that at low temperature due to the solidication of liquid pressure transmitting medium and the different thermal contraction of cell components. Therefore, some corrections have been made in comparison with the data obtained from the constant-loading CAC). For all these high-pressure resistivity measurements, we employed glycerol as the pressure transmitting medium. All resistivity measurements were performed with the conventional four-terminal method with current applied within the ab plane and magnetic eld perpendicular to the ab plane.

Data availability. The data that support the ndings of this study are available on request from the corresponding authors (J.-Q.Y., J.-G.C., or T.S.).

References

1. Hosono, H. & Kuroki, K. Iron-based superconductors: Current status of materials and pairing mechanism. Phys. C 514, 399422 (2015).

2. Davis, J. C. S. & Lee, D.-H. Concepts relating magnetic interactions, intertwined electronic orders, and strongly correlated superconductivity. Proc. Natl. Acad. Sci. USA 110, 1762317630 (2013).

3. Keimer, B., Kivelson, S. A., Norman, M. R., Uchida, S. & Zaanen, J. From quantum matter to high-temperature superconductivity in copper oxides. Nature 518, 179186 (2015).

4. McQueen, T. M. et al. Tetragonal-to-orthorhombic structural phase transition at 90 K in the superconductor Fe1.01Se. Phys. Rev. Lett. 103, 057002 (2009).

5. Imai, T., Ahilan, K., Ning, F. L., McQueen, T. M. & Cava, R. J. Why does undoped FeSe become a high-Tc superconductor under pressure? Phys. Rev.

Lett. 102, 177005 (2009).6. Baek, S.-H. et al. Orbital-driven nematicity in FeSe. Nat. Mater. 14, 210214 (2015).

7. Bhmer, A. E. et al. Origin of the tetragonal-to-orthorhombic phase transition in FeSe: A combined thermodynamic and NMR study of nematicity. Phys. Rev. Lett. 114, 027001 (2015).

8. Hsu, F. C. et al. Superconductivity in the PbO-type structure a-FeSe. Proc. Natl. Acad. Sci. USA 105, 1426214264 (2008).

9. Song, C.-L. et al. Direct observation of nodes and twofold symmetry in FeSe superconductor. Science 332, 14101413 (2011).

10. Kasahara, S. et al. Field-induced superconducting phase of FeSe in the BCS-BEC cross-over. Proc. Natl. Acad. Sci. USA 111, 1630916313 (2014).

11. He, S. L. et al. Phase diagram and electronic indication of high-temperature superconductivity at 65 K in single-layer FeSe lms. Nat. Mater. 12, 605610 (2013).

12. Tan, S. Y. et al. Interface-induced superconductivity and strain-dependent spin density wave in FeSe/SrTiO3 thin lms. Nat. Mater. 12, 634640 (2013).

13. Ge, J.-F. et al. Superconductivity above 100 K in single-layer FeSe lms on doped SrTiO3. Nat. Mater. 14, 285289 (2015).

14. Miyata, Y., Nakayama, N., Sugawara, K., Sato, T. & Takahashi, T. High-temperature superconductivity in potassium-coated multilayer FeSe thin lms. Nat. Mater. 14, 775779 (2015).

15. Shiogai, J., Ito, Y., Mitsuhashi, T., Nojima, T. & Tsukazaki, A. Electric-eld-induced superconductivity in electrochemically etched ultrathin FeSe lms on SrTiO3 and MgO. Nat. Phys. 12, 4246 (2016).

peak|0.5

|( 0 dH/dT| T = T c peak )/T c

(T0.5 K1 )

0.5

0.4

0.3

0.2

0.1

0.0

Kaluarachchi et al. This work

0 2 4 6 8 10P (GPa)

Figure 6 | Mass analysis from HT curves under high pressure. The pressure dependence of the square root of the initial slope of the Tpeakc line

divided by Tpeakc, which corresponds to the effective mass in a simple picture. For comparison, the low-pressure data from ref. 26 are also indicated (blue squares).

analysis (Fig. 6) we also nd a rather at behaviour at high pressure. This suggests that the vortex state of FeSe at high pressures is also nontrivial.

We also note that the unusually strong vortex pinning has been found in clean lms of BaFe2(As1 xPx)2 near the QCP42. Our

H-T phase diagrams at different pressures indicate that the vortex-liquid region between the width of transition (Tc and

Tpeakc) shows a dramatic narrowing when the SDW is suppressed below the superconducting transition (Fig. 4eh, mn), which can also be accounted for by the enhanced pinning near the QCP. However, a peculiar feature here is that the optimal Tc takes place at the point where Tm just crosses Tc, rather than the possible QCP at \8 GPa where Tm completely vanishes. This is again rather similar to the case of hole-doped cuprate superconductors in a sense that a putative QCP is located at the slightly overdoped side3.

To summarize, by using high pressure as a clean tuning knob, we have claried in bulk FeSe the interplay of three competing orders, that is, nematic, magnetic and superconductivity, and elucidated how the high-Tc superconducting phase is achieved.

The relatively low Tc at ambient pressure should be attributed to the presence of competing nematic order; the application of high pressure initially suppresses gradually the nematic order so as to enhance superconductivity. At the same time, the suppression of nematic order also relieves the magnetic frustration so as to stabilize the magnetic order, which then competes again with superconductivity, leading to a local maximum of Tc(P) around 1 GPa. Then, a complete elimination of nematic order around 2 GPa gives rise to a more profound enhancement of Tc to the second plateau of B20 K; the opposite effects on Tc imposed by the strong uctuations of nematic order and the long-range SDW order produce the local minimum of Tc(P) around 1.5 GPa. Above B2 GPa, only the static SDW order is present and competes with superconductivity so that Tc keeps nearly unchanged while Tm rises until

B5 GPa, where the magnetic order at Tm eventually becomes destabilized by pressure. Then, the suppression of magnetic order leads to a great enhancement of Tc to 38 K. In such a way,

FeSe exhibits a unique phase diagram with three competing orders, nematic, SDW and high-Tc superconductivity; the latter two of which have very close energy scales. This newly constructed diagram, which shares some similarities to that of cuprates, may offer important clues for discussing the unconventional origins of the high-Tc superconductivity in this class of materials.

6 NATURE COMMUNICATIONS | 7:12146 | DOI: 10.1038/ncomms12146 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12146 ARTICLE

16. Watson, M. D. et al. Emergence of the nematic electronic state in FeSe. Phys. Rev. B 91, 155106 (2015).

17. Suzuki, Y. et al. Momentum-dependent sign inversion of orbital order in superconducting FeSe. Phys. Rev. B 92, 205117 (2015).

18. Terashima, T. et al. Anomalous Fermi surface in FeSe seen by Shubnikov-de Haas oscillation measurements. Phys. Rev. B 90, 144517 (2014).

19. Watson, M. D. et al. Dichotomy between the hole and electron behaviour in multiband superconductor FeSe probed by ultrahigh magnetic elds. Phys. Rev. Lett. 115, 027006 (2015).

20. Mizuguchi, Y., Tomioka, F., Tsuda, S., Yamaguchi, T. & Takano, Y. Superconductivity at 27K in tetragonal FeSe under high pressure. Appl. Phys. Lett. 93, 152505 (2008).

21. Medvedev, S. et al. Electronic and magnetic phase diagram of b-Fe1.01Se with

superconductivity at 36.7 K under pressure. Nat. Mater. 8, 630633 (2009).22. Bendele, M. et al. Coexistence of superconductivity and magnetism in FeSe1 x under pressure. Phys. Rev. B 85, 064517 (2012).

23. Terashima, T. et al. Pressure-induced antiferromagnetic transition and phase diagram in FeSe. J. Phys. Soc. Jpn. 84, 063701 (2015).

24. Miyoshi, K. et al. Enhanced superconductivity on the tetragonal lattice in FeSe under hydrostatic pressure. J. Phys. Soc. Jpn. 83, 013702 (2014).

25. Terashima, T. et al. Fermi surface reconstruction in FeSe under high pressure. Phys. Rev. B 93, 094505 (2016).

26. Kaluarachchi, U. S. et al. Non-monotonic pressure evolution of the upper critical eld in superconducting FeSe. Preprint at http://arxiv.org/abs/1511.08766

Web End =http://arxiv.org/abs/ http://arxiv.org/abs/1511.08766

Web End =1511.08766 (2015).

27. Jung, S.-G. et al. Enhanced critical current density in the pressure-induced magnetic state of the high-temperature superconductor FeSe. Sci. Rep. 5, 16385 (2015).

28. Shibauchi, T., Carrington, A. & Matsuda, Y. A quantum critical point lying beneath the superconducting dome in iron-pnictides. Annu. Rev. Condens. Matter Phys. 5, 113135 (2014).

29. Bhmer, A. E. et al. Lack of coupling between superconductivity and orthorhombic distortion in stoichiometric single-crystalline FeSe. Phys. Rev. B 87, 180505(R) (2013).

30. Kothapalli, K. et al. Strong cooperative coupling of pressure-induced magnetic order and nematicity in FeSe. Preprint at http://arxiv.org/abs/1603.04135

Web End =http://arxiv.org/abs/1603.04135 (2016).

31. Wang, P. et al. Unlocked stripe-order antiferromagnetism in FeSe under pressure. Preprint at http://arxiv.org/abs/1603.04589

Web End =http://arxiv.org/abs/1603.04589 (2016).

32. Matsubayashi, K. et al. Pressure-induced heavy fermion superconductivity in the nonmagnetic quadrupolar system PrTi2Al20. Phys. Rev. Lett. 109, 187004 (2014).

33. Cheng, J.-G. et al. Integrated-n gasket for palm cubic-anvil high pressure apparatus. Rev. Sci. Instrum. 85, 093907 (2014).

34. Cheng, J.-G. et al. Pressure induced superconductivity on the border of magnetic order in MnP. Phys. Rev. Lett. 114, 117001 (2015).

35. Kumar, R. S. et al. Crystal and electronic structure of FeSe at high pressure and low temperature. J. Phys. Chem. B 114, 1259712606 (2010).

36. Uhoya, W. et al. Simultaneous measurement of pressure evolution of crystal structure and superconductivity in FeSe0.92 using designer diamonds. Eur. Phys.

Lett. 99, 26002 (2012).37. Wang, Q. et al. Strong interplay between stripe spin uctuations, nematicity and superconductivity in FeSe. Nat. Mater. 15, 159163 (2015).

38. Glasbrenner, J. K. et al. Effect of magnetic frustration on nematicity and superconductivity in iron chalcogenides. Nat. Phys. 11, 953958 (2014).

39. Hiraishi, M. et al. Bipartite magnetic parent phases in the iron oxypnictide superconductor. Nat. Phys. 10, 300303 (2014).

40. Putzke, C. et al. Anomalous critical elds in quantum critical superconductors. Nat. Commun. 5, 5679 (2014).

41. Kogan, V. G. & Prozorov, R. Orbital upper critical eld and its anisotropy of clean one- and two-band superconductors. Rep. Prog. Phys 75, 114502 (2012).

42. Kurth, F. et al. Unusually high critical current of clean P-doped BaFe2As2 single crystalline thin lm. Appl. Phys. Lett. 106, 072602 (2015).

43. Chareev, D. et al. Single crystal growth and characterization of tetragonal FeSe1 x superconductors. Cryst. Eng. Comm. 15, 19891993 (2013).

Acknowledgements

We thank T. Terashima, X.J. Zhou, T. Xiang, R. Yu, Q.M. Zhang, G. Chen, andJ.-S. Zhou for very helpful discussions. We also thank Bosen Wang for his technical help. Work at IOP CAS was supported by the National Basic Research Program of China (Grant No. 2014CB921500), National Science Foundation of China (Grant No. 11574377, 11304371), and the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDB07020100) as well as the Opening Project of Wuhan National High Magnetic Field Center (Grant No. 2015KF22), Huazhong University of Science and Technology. Work in Japan was supported by Grant-in-Aids for Scientic Research (A), (B), (S), and on Innovative Areas Topological Materials Science. Work at ORNL was supported by the US Department of Energy, Ofce of Science, Basic Energy Sciences, Materials Sciences and Engineering Division.

Author contributions

J.-Q.Y., J.-G.C. and T.S. conceived this project. J.P.S., G.Z.Y., J.-G.C. performed the high-pressure magneto-transport measurements with the self-clamped PCC and CAC;K.Matsuura, Y.Mizukami, M.S., K.Matsubayashi, M.Y., Y.U. measured the high-pressure resistivity with the constant-loading CAC; J.-Q.Y. and B.C.S. synthesized the FeSe single crystals with the ux method; K.Matsuura, T.W., S.K., Y.Matsuda grew the FeSe single crystals with the vapor transport technique. All authors discussed the results. J.-G.C. andT. S. wrote the paper with inputs from all authors. J.-G.C., Y.Matsuda, Y.U. and T. S. supervised the project.

Additional information

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Sun, J. P. et al. Dome-shaped magnetic order competing with high-temperature superconductivity at high pressures in FeSe. Nat. Commun. 7:12146 doi: 10.1038/ncomms12146 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2016

NATURE COMMUNICATIONS | 7:12146 | DOI: 10.1038/ncomms12146 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7