Arxiv:1908.05448V1 [Cond-Mat.Supr-Con] 15 Aug 2019 ﬀc.Svrlognccnutr,Sc As Such Screening Conductors, Weak Organic the Several Thus and Density, Eﬀect

arXiv:1908.05448v1 [cond-mat.supr-con] 15 Aug 2019 ffc.Svrlognccnutr,sc as such screening conductors, weak organic the Several thus and density, effect. carrier because low suited their best of the are superconductors lay- organic ered superconductivity, on interactions mecha- electron-electron pairing SC is of states possibility SC the nism. and expand CO the to coexisting between crucial even relationship and the neighboring, for identification exper- Further imental fluc- superconductivity. the between and charge interplay instability the the from charge suggesting superconductor in [6] study, organic increase optical the an the in observed fact, is In the tuations tem- when [5] low very formation to pair peratures. suppressed is SC transition the transi- ordering to charge ordering contribute charge can the that tion near shows fluctuations model theoretical charge Hubbard a extended the though, the Even on state. based SC study strong the itin- of the the nature with with erant conflicts severely coincides which ordering localization, difficult carrier charge is the realization because experimental in- charge However, the features near SC stability. superconductivity intriguing for critical found that magnetic be expect would the naturally near we super- found points, unconventional been many have Since conductors [4] mag- contributes instability. the interaction netic near on-site superconductivity unconventional the form as to manner in interaction same pairing the the (SC) superconducting be that the address can of to interactions source us electron-electron leading the long-range [1–3] state, in the CO appears insulat- the (CO) sometimes of vicinity ordered Superconductivity charge interac- state. the for ing Coulomb responsible long-range is the tion in band, While quarter-filled state. insulating a Mott-Hubbard the moving stabilize from on- to carriers strong the the prevents band, repulsion half-filled Coulomb a site of case in- the electron-electron In strong the teractions. by caused is band duction oe he-odC tt,[0 hc a ananthe maintain can which [10] state, CO three-fold a posed in.[,79 hoeia tde for studies Theoretical transi- 7–9] [3, ordering charge tions. apparent show trathiafulvalen) θ tp ETTFsls(ETTF bisethyldithiote- (BEDT-TTF: salts BEDT-TTF -type oivsiaeeprmnal h ffc flong-range of effect the experimentally investigate To con- filled partially a in transition metal-insulator The irsoi oxsec fsprodciiyadcag o charge and superconductivity of coexistence Microscopic eateto hsc,Gaut colo cec,Hokkaid Science, of School Graduate Physics, of Department TTF) hscnutv hreodrdsaei investigated. is Rela state molecules. ordered mol other charge which on conductive in s this carriers state, ordered significan conducting ordered charge charge not with exotic is the coexists the conductivity In by in-plane explained the is material. while one gapped, in is state ordered charge h lcrnprmgei eoac td o nogncsu organic an for study resonance paramagnetic electron The superconductor 4 [(H 3 O)Ga(C 2 .Iaa .Mrb,S uuk,adA Kawamoto A. and Fukuoka, S. Moribe, K. Ihara, Y. O 4 ) 3 ] · C 6 β H θ ′′ tp at pro- salts -type 5 -(BEDT-TTF) NO α 2 tp and -type Dtd uut1,2019) 16, August (Dated: eel htsprodciiyceit nfrl ihth with uniformly coexists superconductivity that reveals odC tt a o enfudin found been not will has state state three- CO liquid such fold however, pinball Experimentally, [11–13] the unconventional as realized. be charges, to carri- referred localized state itinerant the metallic the with Further when interact sites. that molecular ers show three carri- of studies one the theoretical at of localized part are some ers when even conductivity metallic h o-rva Osaevr ls oteS transition SC the to show SC close which very the temperature. materials state to CO impact study non-trivial its should the 11] and we [10, state formation, CO phase. pair novel horizontal vertical the the the reveal with of To limited competes stability a phase the at stripe only where appears space, state parameter CO three-fold the cause TTF) ntblt tatmeauevr ls oteSC the to close very and temperature order temperature charge a transition the at charge explore shows instability to compound this material because superconductivity, ideal an is eouin ewr o bet xld h possibility the exclude to able not microscopic were a spectrum we insufficient from the resolution, of state because However, electronic The viewpoint. technique the powerful investigate most coexists the not. to state of one CO or is the spectroscopy superconductivity the NMR if from clarify with separated should is we microscopically sample [19] the of segregation, part, resistivity part SC phase the CO of the as is, However, possibility that a realization state. suggest experimental liquid experiments the pinball for the candidates of best the of hs rniina h hreodrn temperature ordering charge the magnetic at to T due transition anomaly no phase detected experiment EPR a ensgetd 1]A h hrelclzto in localization below and conductivity, charge superconductivity metallic with the the even with As the coexists state [18] pattern CO between charge the suggested. three-fold ratio the been intensity [16] sites, has poor NMR transition. and rich SC the charge the from near fluctuations Besides, charge the in fteNRsetu pitn.[8 lo h NMR the dynamics Also, spin [18] at low-energy increases the splitting. that showed spectrum experiments NMR the of 4 CO [(H aee rai superconductor organic layered A 8 = 4 3 inhpbtensprodciiyand superconductivity between tionship [(H l oie.Ti nstoi behavior anisotropic This modified. tly O)Ga(C nvriy apr 6-80 Japan 060-0810, Sapporo University, o clrst eetv are localization carrier selective ecular-site . 5 3 K O)Ga(C ae h nepaesi exchange spin interplane the tate, 1,1]wihwsdtrie steonset the as determined was which 17] [16, , T CO hc ssgetv fteincrease the of suggestive is which , 2 perconductor 2 O O 4 4 ) ) 3 3 T ] ] c · · C C dri organic in rder .Peiu M and NMR Previous K. 7 = 6 H 6 H 5 NO 5 β NO ′′ 2 -(BEDT- T θ c tp at,[4 be- [14] salts, -type ( , 2 β β ′′ ′′ G)sl [15] salt -Ga) G ati one is salt -Ga e ∗ β ′′ -(BEDT- 2

(a) (b) β''-Ga salt b axis 2.012 A layer c* axis H || 45 deg. from b

factor 2.008 g b axis 20 K B layer 2.004

c* axis 0 60 120 180 10 K Angle from b axis (deg.)

EPR Intensity (c) ET molecules 8 K A layer

Insulating layer

c* 7 K B layer 341 342 343 344 345 H (mT ) b 0 6 K EPR intensity EPR (arb .units)

3.6 K FIG. 1: (a) X-band EPR spectra of β′′-Ga salt at 3.6 K in various external field directions. The two-peak spectrum observed in arbitral field direction merges into a single peak in fields along the symmetric axes. (b) Field direction de- 343 344 345 346 pendence of the g factors determined from the EPR peak positions. Two branches with sinusoidal angle dependence µ0H (mT) originate from the A and B layers. (c) The double-layer ′′ crystal structure of β -Ga salt. The conducting BEDT-TTF FIG. 2: The EPR peak splitting observed below TCO = 8.5 molecules layers are separated by the insulating anion layers. ◦ The long axis of the BEDT-TTF molecules in layer A (B) is K in the 45 field, where the peak separation at the low- ◦ ◦ inclined by 76 (−76 ) with respect to the b axis, which coin- est temperature becomes the maximum. The spectrum split- cides with the maximum in the angle dependence of g factor ting without the increase in the spectrum width indicates the (b). anomaly in the interplane cross relaxation (see text) with neg- ligible modification in the in-plane spin relaxation. of phase segregation An alternative probe with higher resolution was desired. tuned by the single axis rotator, with which the sam- In this study, we show that the X-band electron para- ple can be rotated in situ around the crystalline a axis. magnetic resonance (EPR) experiment is one such probe. We also performed the in-plane (b axis) and interplane The EPR signal in β′′-Ga salt originates from the π (c∗ axis) resistivity measurements by the conventional 4 electrons in the highest occupied molecular orbital of probe method. [20] the BEDT-TTF molecules. We succeed in detecting the Figure 1 (a) shows the EPR spectra at the lowest tem- charge anomaly on the clearly resolved EPR spectrum by perature of 3.6 K in the bc∗ plane fields. The field di- taking advantage of the anisotropy of g factors and the rection is determined by the angle from the b axis. The bilayer crystal structure of β′′-Ga salt. Thus, the EPR angular dependence of the EPR spectrum originates from experiment allows us to observe the phase segregation, if the anisotropic g tensor, for which the principal axes any, as the additional component of the EPR spectrum. coincide with the symmetric axes of the BEDT-TTF The present results, which are explained by a single EPR molecule, and the principal values were determined for contribution at any temperatures, clearly evidence the β-(BEDT-TTF)2I3 salt as g = (2.011, 2.008, 2.002). [21] uniform coexistence between SC and CO states. We also The g factors for each field direction were determined conducted the resistivity measurement to confirm that by fitting the spectra with two-peaks Lorenzian func- the CO state actually involves the conducting carriers, tion. As the result, two sinusoidal branches were ob- because the EPR experiment is in principle a spin sensi- tained as shown in Fig.1(b). These branches are assigned tive probe. The results of resistivity and EPR measure- to the EPR signals from the A and B layers displayed in ments are explained consistently, which unambiguously Fig. 1 (c). Assuming that the principal axes of g tensor suggest a SC state coexisting with the conductive CO are fixed to the BEDT-TTF molecules, we can deter- state. mine the principal values of g tensor for β′′-Ga salt as The single crystalline samples were grown by the stan- g = (2.016, 2.010, 2.001). Good agreement with previ- dard electrochemical reaction. [15] The X-band EPR ex- ous study [21] suggests that the entire EPR spectra can periments were performed with a commercial spectrom- be explained by the anisotropic g tensor and the double- eter (Bruker EMX Plus). One single crystal with a di- layer crystal structure of β′′-Ga salt, meaning that our mension of 2 × 0.5 × 0.2 mm3 was used for the experi- crystal is single phase and the electronic state is uniform. ment. The orientation of the external magnetic field was Next, we measured the temperature dependence of the 3

change because of the extremely high upper critical field 344.7 (a) β''-Ga salt (Bc2 > 30 T). The effect of SC transition was observed T = 8.5 K H || 45 deg. from b CO in the EPR spectrum as the reduction of the integrated 344.4 intensity below 7 K (Fig. 3 (b)). Such decrease in in- ◦ tensity was not observed in the 45 field, because Tc is suppressed below 3 K. This result confirms that the elec- 344.1 tronic spins that would show superconductivity in zero field contribute to the EPR intensity when superconduc- Peak positions positions Peak (mT) tivity is suppressed by 45◦ field. Thus, if SC part of 343.8 the sample did not show the charge ordering transition, T = 7 K which is the case for the macroscopic phase segregation, (b) c 1.0 an additional EPR peak originating from the electrons in a normal metallic state should be observed at the center of the two-peak spectrum. However, such extra contribution was not observed at the lowest temperature of 3.6 0.5 K, as shown at the bottom of Fig. 2. The clear two-peak ◦ H || b spectrum in the 45 field allows us to conclude that the H || 45 deg. from b SC state coexists uniformly with the CO state. We note that the EPR intensity decreases gradually below Tc in EPR intensity intensity EPR units) (arb. 0.0 the b-axis field, and finite intensity remains even at 3.6 0 5 10 15 20 K. This behavior is contrasting to the conventional be- T (K) havior expected for a homogeneous SC state, in which EPR signal should disappear. The EPR in the SC state FIG. 3: (a) The temperature dependence of the EPR peak may originate from the nearly localized electrons in the fields. Abrupt peak splitting was observed below TCO = 8.5 CO state, for instance, the pin site in the pinball liquid K. (b) The temperature dependence of the EPR intensity in state. We exclude the possibility that the entire EPR ◦ the b-axis field, and 45 field. The intensity is normalized intensity originates from the free spins induced by impu- at 12 K to correct the intensity difference due to the sample rity or defects, because the EPR intensity increases from shape. The SC transition is detected as the reduction of EPR 20 K up to room temperature, which is contrast to the intensity below Tc in the b axis field. Curie behavior expected for the free spins. [20] In fact, a typical Curie-Weiss type temperature dependence was observed above Tc in the C and Al-doped MgB2, [22– EPR spectrum in the field applied to the direction 45◦ 24] in which the EPR signal originates from the defects rotated from the b axis (45◦ field). In this field direction, and/or impurities. Tc is suppressed below 3 K by a field of approximately Now, to study in detail the origin of the EPR peak 300 mT. With increasing temperatures, the peak separa- splitting in the CO state, we fit the EPR spectra using tion becomes small, as shown in Fig. 2, and a single peak the Bloch model with cross spin relaxation between the was observed at the temperatures higher than 8.5 K. We neighboring layers A and B (Fig. 1(c)). [25] The time found a trace of the two-peak structure at 8 K as the evolution of the electron magnetization on the weakly wiggle around the center of the spectrum at 344.4 mT. coupled layers A and B, M A, M B are written as Therefore, the peak positions were determined by the A A AB two-component Lorentzian fit for the spectra below 8 K dM g µB ∆M = M A × B − RA + , (1) (filled symbols in Fig.3 (a)), and by the single compo- dt h TX nent Lorentzian fit above 8.5 K (open symbols in Fig.3 M B B M BA d g µB B B ∆ (a)). The abrupt increase in the peak separation below = M × B − R + , (2) dt h T 8 K clearly evidences a phase transition. This anomaly X agrees with the charge ordering transition at T = 8.5 Rα α α α − CO where = (Mx /T2,My /T2, (Mz M0)/T1) with α = K previously observed from the 13C NMR study. [18] AB B A α A, B layers, ∆M = M −M , and g , µB, h are the g The origin of the EPR peak splitting will be discussed factor for the α layer, the Bohr magneton and the Plank’s later. constant. When the intrinsic spin relaxation time T2 be- As the charge ordering anomaly is successfully ob- comes short, the whole spectrum is broadened in pro- served in the EPR spectrum, then we compare the EPR portion to 1/T2, and the broad spectrum smears out the spectra in b-axis and 45◦ fields to unveil the relationship two-peak structure. Contrastingly, when the cross relax- between the CO and SC states. Typical EPR spectra ation time TX becomes faster than the timescale equiva- for b-axis fields are presented in the supplimental ma- lent to the spectrum separation, the two-peak structure terial. [20] In the b-axis field of 345 mT, Tc does not is lost due to the dynamical narrowing effect, and a single 4

The anisotropy developing in the CO state is also ob- (a) b served from the electrical transport experiment, as shown TCO = 8.5 K ''-Ga salt 10 in Fig. 4 (b). [20] The in-plane resistance (Rperp) and interplane resistance (Rpara) were measured in magnetic (MHz)

X ﬁelds of 0.3 T and 2.5 T applied perpendicular to the T

1/ conducting plane to suppress Tc. The increase in the 1 ratio R⊥/R|| below TCO means that the interplane resistivity increases due to the gap opening in the inter- 1200 plane transfer channel. A kink was observed at a lower (b) temperature of 4.5 K because of the precursor to the 900 SC transition. The anomaly at TCO was not clearly para observed in the previously reported in-plane resistivity R / 600 measurements, [15, 16] because the highly conducting

perp in-plane channel is not significantly modified in the CO R 300 state, which is consistent with the invariant T2. The highly conducting in-plane channel is also suggested from 0 the quantum oscillation study, in which clear Shubnikov- 0 5 10 15 20 de Haas oscillation was observed even in the CO state, T (K) where resistivity shows a semiconducting temperature dependence. [19, 26] Observation of the Fulde-Ferrell- FIG. 4: The temperature dependence of (a) the cross relax- Larkin-Ovchinnikov superconducting state in high mag- ation rate 1/TX, and (b) the ratio between the in-plane and netic fields is anothor evidence for a clean electronic state interplane resistance. The gray area in (a) represents the fre- in the conducting plane. [27] quency that corresponds to the EPR peak separation at the Finally, a question is whether the CO state supports or lowest temperature, above which 1/TX cannot be determined suppress superconductivity. In a series of β′′ type BEDT- because of the motional narrowing. At low temperatures, TTF salts with various guest molecules and metallic ions, 1/TX follows an exponential function with a gap energy of 16 K as shown by the dotted line. Tc higher than 5 K is found only in the salts with resistivity upturn at low temperatures. [26] Whereas in the ′′ β -Rh salt with Tc =2.5 K, almost no resistivity upturn was observed. [28] As the charge ordering transition in- sharp peak will be observed. The sharp EPR spectra at creases the SC transition temperatures, we suggest the high temperatures above 8.5 K in the present study is ex- importance of the exotic electronic state with the spatial plained by this dynamical narrowing effect. Figure 4 (a) charge modulation [12, 13] to induce this unconventional shows the temperature dependence of TX determined by type of superconductivity. fitting the experimental spectra with the coupled Bloch To conclude, we performed the EPR experiment for the model. [20] The gray region in Fig. 4 (a) represents the organic superconductor β′′-Ga salt, and detected clearly timescale shorter than the spectrum separation of ap- the charge ordering transition at T =8.5 K. In the CO proximately 20 MHz. CO state, the interplane spin exchange channel is gapped as In the CO state, the temperature dependence of 1/TX the result of the partial charge localization. In the same shows an exponential behavior with a gap energy of sample, we observed the SC transition as the reduction of ∆/kB = 16 K. This gap size is consistent with that ex- EPR intensity below Tc. We found only a single spectrum pected for the weak coupling charge density wave state, component both in the SC and CO states, which sug- 2∆/kBTCO =3.8. The charge gap was detected in the in- gests a uniform coexistence between superconductivity terplane spin exchange channel because the already small and charge order. To understand this exotic coexisting interplane transfer integral in the metallic state is com- state, we suggest a three-fold CO state, in which partial pletely lost in the CO state due to the partial charge charge localization coexists with the high conductivity. localization. Whereas for the in-plane spin relaxation We would like to acknowledge Y. Oshima, L. Martin T2, we found that the spectrum width is invariant above for fruitful discussions, and M. Fujiwara for the support and below TCO, which leads us to conclude that T2 is in carrying out the EPR experiment in the Institute for not significantly modified by the charge ordering tran- Molecular Science. This study was partly supported by sition. The invariant T2 allows us to exclude the pos- the Suhara memorial foundation. sibility of magnetic phase transition. This anisotropic behavior is consistently explained by the three-fold CO state, in which partial carrier localization immediately extinguishes the weak interplane transfer channel, while ∗ maintaining the coherent in-plane transfer channel. Electronic address: [email protected] 5

[1] H. Kobayashi, R. Kato, A. Kobayashi, Y. Nishio, K. Ka- [17] Y. Ihara, Y. Futami, and A. Kawamoto, Journal of the jita, and W. Sasaki, Chemistry Letters 15, 789-792 Physical Society of Japan 85, 014601 (2016). (1986). [18] Y. Ihara, M. Jeong, H. Mayaffre, C. Berthier, M. Hor- [2] M. Kurmoo, M. J. Rosseinsky, P. Day, P. Auban, vatić, H. Seki, and A. Kawamoto, Physical Review B 90, W. Kang, D. Jérôm, and P. Batail, Synthetic Metals 27, 121106(R) (2014). A425 (1988). [19] A. I. Coldea, A. F. Bangura, J. Singleton, A. Ardavan, [3] H. Mori, S. Tanaka, and T. Mori, Physical Review B 57, A. Akutsu-Sato, H. Akutsu, S. S. Turner, and P. Day, 12023 (1998). Physical Review B 69, 085112 (2004). [4] N. D. Mathur, F. M. Groshe, S. R. Julian, I. R. Walker, [20] See Supplemental Material [ulr], which includes Refs. [29, D. M. Freye, R. K. Haselwimmer, and G. G. Lonzarich, 30]. Nature (London) 394 39 (1998). [21] N. Kinoshita, M. Tokumoto, H. Anzai, and G. Saito, [5] J. Merino, and R. H. McKenzie, Physical Review Letters Journal of the Physical Society of Japan 54, 4498 (1985). 87, 237002 (2001). [22] A. Bateni, S. Repp, R. Thomann, S. Acar, E. Erdem, and [6] S. Kaiser, M. Dressel, Y. Sun, A. Greco, J. A. Schlueter, M. Somer, Applied Physics Letters 105, 202605 (2014). G. L. Gard, and N. Drichko, Physical Review Letters [23] A. Bateni, E. Erdem, S. Repp, S. Acar, I. Kokal, 105, 206402 (2010). W. Häßler, S. Weber, and M. Somer, Journal of Applied [7] T. Kakiuchi, Y. Wakabayashi, H. Sawa, T. Takahashi, Physics 117, 153905 (2015). and T. Nakamura, Journal of the Physical Society of [24] A. Bateni, E. Erdem, S. Repp, S. Weber, and M. Somer, Japan 76, 113702 (2007). Applied Physics Letters 108. 202601 (2016). [8] Y. Nogami, J.-P. Pouget, M. Watanaba, K. Oshima, [25] A.´ Antal, T. Fehér, E. Tátrai-Szekeres, F. Fülöp, H. Mori, S. Tanaka, and T. Mori, Synthetic Metals 103, B. Náfrádi, L. Forró, and A. Jánossy, Physical Review 1911 (1999). B 84, 075124 (2011). [9] K. Miyagawa, A. Kawamoto, and K. Kanoda, Physical [26] A. F. Bangura, A. I. Coldea, J. Singleton, A. Ardavan, 62 Review B , R7679 (2000). A. Akutsu-Sato, H. Akutsu, S. S. Turner, P. Day, T. Ya- 72 [10] T. Mori, Journal of the Physical Society of Japan , mamoto, and K. Yakushi, Physical Review B 72, 014543 1269 (2003). (2005). [11] M. Kaneko, and M. Ogata, Journal of the Physical Soci- [27] S. Uji, Y. Iida, S. Sugiura, T. Isono, K. Sugii, 75 ety of Japan , 014710 (2006). N. Kikugawa, T. Terashima, S. Yasuzuka, H. Akutsu, [12] C. Hotta, N. Furukawa, A. Nakagawa, and K. Kubo, Y. Nakazawa, D. Graf, and P. Day, Physical Review B 75 Journal of the Physical Society of Japan , 123704 97, 114505 (2018). (2006). [28] L. Martin, A. L. Morritt, J. R. Lopez, H. Akutsu, [13] J. Merino, A. Greco, N. Drichko, and M. Dressel, Physi- Y. Nakazawa, S. Imajo, and Y. Ihara, Inorganic Chem- 96 cal Review Letters , 216402 (2006). istry 56, 717 (2017). [14] M. Watanabe, Y. Noda, Y. Nogami, and H. Mori, [29] Y. Goto, A. Miura, C. Moriyoshi, Y. Kuroiwa, T. D. Mat- 73 Jouranl of the Physical Society of Japan , 116 (2004). suda, Y. Aoki, and Y. Mizuguchi, Scientific Reports 8, [15] H. Akutsu, A. Akutsu-Sato, S. S. Turner, D. Le Peve- 12852 (2018). len, P. Day, V. Laukhin, A.-K. Klehe, J. Singleton, [30] A. J. S. Machado, S. T. Renosto, C. A. M. dos San- D. A. Tocher, M. R. Probert, and J. A. K. Howard, Jour- tos, L. M. S. Alvas, and Z. Fisk, “Superconductors – 124 nal of the American Chemical Society , 12430 (2002). Materials, Properties and Applications–” ch. 3. Inte- [16] Y. Ihara, H. Seki, and A. Kawamoto, Journal of the Phys- chOpen (2012). ical Society of Japan 82, 083701 (2013).