arxiv: v2 [cond-mat.supr-con] 9 Dec 2017

Transcription

1 arxiv: v2 [cond-mat.supr-con] 9 Dec 217 Transport properties of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 single crystals in the mixed state Shanshan Sun 1, Shaohua Wang 1, Chenghe Li 1, and Hechang Lei 1, 1 Department of Physics and Beijing Key Laboratory of Opto-electronic Functional Materials & Micro-nano Devices, Renmin University of China, Beijing 1872, China hlei@ruc.edu.cn Abstract. We study the electric transport properties of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 single crystals with z = and.6 in the mixed state. Thermally-activated flux-flow, vortex glass and flux-flow Hall effect (FFHE) behaviors are observed. Experimental results show that there are rich vortex phases existing in these systems and the vortex liquid states occupy broad regions of phase diagrams. Further analysis suggests that thermal fluctuation plays an important role in the vortex phase diagrams of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2. Moreover, for Li x (NH 3 ) y Fe 2 Se 2, there is no sigh reversal of FFHE in the mixed state and a scaling behavior ρ xy (µ H) = Aρ xx (µ H) β with β 2. is observed.

2 Transport properties of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 single crystals in the mixed state 2 1. Introduction Since the discovery of LaOFeAs[1], iron-based superconductors (IBSCs) have became another family of high-t c superconductors besides cuprate SCs and attracted intense attention. On the one hand, because of high superconducting transition temperature T c and large upper critical fields µ H c2, IBSCs are important not only for basic science but also for practical application. On the other hand, the values of Ginzburg number Gi ranging from 1 4 to 1 2 in IBSCs indicate that the effect of thermal fluctuation varies in different members of IBSCs. It leads to rich vortex phase diagrams in these SCs[2]. Among IBSCs, the iron-chalcogenide SC FeSe has a relative low and nearly isotropic µ H c2 with low T c ( 8 K)[3, 4], when compared to iron-pnictide SCs[5]. After intercalating alkali metals A in between FeSe layers using high-temperature synthesis method, the T c can be dramatically enhanced up to about 31 K for A x Fe 2 y Se 2 [6]. However, the obvious mesoscopic phase separation between the superconducting phase and the inter-grown antiferromagnetic insulating phase in these compounds make it difficult to study on their intrinsic physical properties[7]. On the other hand, superconducting (Li 1 x Fe x OH)FeSe synthesized by a novel hydrothermal method exhibits the features of high T c (> 4 K), rather large µ H c2 (> 6 T for H c) and Gi ( )[8, 9, 1, 11]. It results in the existence of various vortex phases in (Li 1 x Fe x OH)FeSe single crystals[11, 12]. Besides A x Fe 2 y Se 2 and (Li 1 x Fe x OH)FeSe, AM-NH 3 cointercalated FeSe (AM = alkali, alkali-earth, and rare-earth metals) forms another class of iron-chalcogenide SCs with high T c (> 4 K)[13, 14, 15, 16, 17]. But the absent of single crystals impedes the study on transport properties and vortex dynamics of these SCs. Very recently, we have grownli x (NH 3 ) y Fe 2 Se 2 singlecrystalssuccessfully anditshowsratherhight c andµ H c2 with significant anisotropy[19], similar to (Li 1 x Fe x )OHFeSe. In this work, we present a comprehensive study on the electric transport properties of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 single crystals with z = and.6 (denoted as LiFeSe-122 and LiFeTeSe-122 for brevity) in the mixed state. The observed thermally-activated flux-flow (TAFF), vortex glass (VG) and flux-flow Hall effect (FFHE) behaviors indicate that there are diverse vortex phases and rather broad vortex liquid region in these systems. Large Gi in LiFeSe-122 suggests the thermal fluctuation is essential to vortex dynamics. Moreover, the scaling behavior of ρ xy (µ H) = Aρ xx (µ H) β with β 2 implies that the strength of pinning force is relatively weak in LiFeSe-122 when H c and temperature is close to T c. On the other hand, with Te doping, the Gi decreases accompanying with decreased T c. 2. Experiment Single crystals of LiFeSe-122 and LiFeTeSe-122 were synthesized by the low-temperature ammonothermal technique. The detailed experimental procedure and characterizations of crystals were described in previous work[18, 19]. The elemental analysis was

3 Transport properties of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 single crystals in the mixed state 3 (a) ln[ ( cm)] ln[ ( cm)] K (m cm).2 T c, onset = 44.59K z = /T (1/K) K H(T)= z = /T (1/K).8 Figure 1. (a) and The natural logarithm of ρ(t,µ H) of LiFeSe-122 as a function of inverse of temperature at various fields for H ab and H c, respectively. The orange solid and black dashed lines are fitting results from the Arrhenius relation and Eq. (5). Inset : temperature dependence of in-plane resistivity ρ(t) for LiFeSe-122 at zero field. performed using the inductively coupled plasma atomic emission spectroscopy (ICP- AES) and the atomic ratio of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 single crystals used in this study is Li : Fe : Se =.18 : 1 :.9 for z = and Li : Fe : Te : Se =.16 : 1 :.6 :.38 for z =.6, respectively. Electrical transport measurements were performed using a four-probe configuration with current flowing in the ab plane of crystals in a Quantum Design PPMS-14. The Hall resistivity was obtained from the difference of the transverse resistivity measured at the positive and negative fields in order to remove the longitudinal resistivity contribution due to voltage probe misalignment, i.e., ρ xy (µ H) = [ρ(+µ H) ρ( µ H)]/2. 3. Results and Discussions Figs. 1 and 2 present the longitudinal resistivity ρ(t,µ H) of LiFeSe-122 and LiFeTeSe- 122 single crystals near the superconducting transition region for H ab and H c, respectively. The field-induced broadenings of resistive transitions and resistive tail behaviors are obvious, especially for the field along c axis. Similar behaviors have also been observed in (Li, Fe)OHFeSe, SmFeAsO.85 and cuprates[9, 12, 2, 21]. It can be ascribed to the field-induced TAFF. Based on theoretical model describing the TAFF behavior, the resistivity in TAFF region can be expressed as[22, 23, 24], ρ = (2ν LB/J)exp( J c BVL/T)sinh(JBVL/T) (1)

4 Transport properties of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 single crystals in the mixed state 4 (a) ln[ ( cm)] ln[ ( cm)] K (m cm) 2 T c, onset = 2.78K z = /T (1/K) 2.1K H (T)= /T (1/K) 1 z =.6 Figure 2. (a) and lnρ(t,µ H) vs. 1/T for LiFeTeSe-122 at various fields with H ab and H c, respectively. The orange solid and black dashed lines are fitting results from the Arrhenius relation and Eq. (5). Inset : temperature dependence of in-plane resistivity ρ(t) for LiFeTeSe-122 at zero field. where ν is an attempt frequency for a flux bundle hopping, L is the hopping distance, B is the magnetic induction, J is the applied current density, J c is the critical current density in the absence of flux creep, V is the bundle volume and T is the temperature. If J is small enough and JBVL/T 1, the Eq. (1) can be expressed as ρ = (2ρ c U/T)exp( U/T) = ρ f exp( U/T) (2) where U = J c BVL is the thermally-activated energy (TAE) or flux pinning energy, and ρ c = ν LB/J c, which is usually considered to be temperature independent. According to the condensation model, U(T,µ H) = U (µ H)(1 t) q, where t = T/T c (T c is the superconducting transition temperature), q = 2 n/2 (n depends on the dimensionality of the vortex system with the range from to 3) and U (µ H) is the apparent activated energy[22, 24, 25]. Generally assuming n = 2 and the prefactor ρ f is a constant as for cuprate SCs[22], we can obtain and lnρ(t,µ H) = lnρ (µ H) U (µ H)/T (3) lnρ (µ H) = lnρ f +U (µ H)/T c (4) Thus, there is a linear relationship between lnρ(t,µ H) and 1/T in TAFF region (Arrhenius relation) and the slope and y-axial intercept correspond to U (µ H) and lnρ (µ H), respectively. As shown in Fig. 1(a) and, the solid lines represent the

5 Transport properties of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 single crystals in the mixed state 5 results of linear fitting in TAFF region using the Arrhenius relation. All the linear fittings cross at T cross approximately, which is about K and K for H ab and H c, respectively. Ideally, all the lines at different fields should be crossed into one same point, T cross, which should equal to T c [26]. Obviously, the values of T cross are consistent with the T c,onset in the ρ(t) curves for both field directions (44.59 K and K). Although the lnρ(t,µ H) vs. 1/T can be fitted linearly, there are relative large fitting errors for H c because the Arrhenius relation can only be satisfied in a narrow region. It suggests that the assumptions of linear temperature dependence of U(T,µ H) and temperature-independent ρ f may not be valid, leading to lnρ(t,µ H) vs. 1/T deviating from Arrhenius relation[27]. Using the relation U(T,µ H) = U (µ H)(1 t) q, the Eq. (2) can be rewritten as lnρ = ln(2ρ c U )+qln(1 t) lnt U (1 t) q /T (5) where ρ c and U are temperature independent and value of T c is derived from Arrhenius relation. All fits (black dashed lines) using Eq. (5) are also plotted in Fig. 1(a) and. It can be seen that this more general method is in better agreement with experimental data than Arrhenius relation, especially for H c. Same analysis procedure can be applied to LiFeTeSe-122. As shown in Fig. 2 (a) and, T cross obtained from the linear fitting is about 2.28 K and 2.1 K for H ab and H c, respectively. The values of T cross are close to the T c,onset in the ρ(t) curves (2.78 K and 2.79 K). In contrast to LiFeSe-122, the lnρ(t,µ H) of LiFeTeSe-122 for both field directions are not linearly proportional to 1/T, i.e., the Arrhenius relation can not capture the bending trends of curves. However, the data can be fitted well by using Eq. (5). ThefittedU andq atvariousfieldsforbothsamplesareshowninfig. 3. According totheanderson-kimmodel, U (µ H)isrelatedtotheeffectivepinning energy. U (µ H) for H ab are much larger than that for H c in both crystals (Fig. 3(a) and ), confirming that the flux pinning is much stronger for H ab. It is consistent with the layered structure of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2. When compared with U (µ H) 25 K at µ H =.5 T for FeSe[28], the U (µ H) for H c at same field enhances to about 1 4 K in LiFeSe-122, indicating a much stronger vortex pinning energy accompanying with an increased T c after intercalating the Li-NH 3 layers in between FeSe layers. In contrast, the value of U (µ H) remains in the same order of magnitude ( K) after intercalating the Li-NH 3 layers into FeTe 1 x Se x [29, 3]. That is to say, their effective pinning energy are comparable. According to Kramer s law[31], the field dependence of U (µ H) can be expressed as U (µ H) = a(µ H) γ (1 µ H/µ H irr ) δ (6) where µ H irr is an irreversible field and γ, δ are scaling parameters. In LiFeSe-122, when assuming δ = 2, the field is parallel to ab plane, the value of γ is.57(3) and.67(5) for H ab and H c, respectively (Fig. 3(a)). On the other hand, the value of γ

6 Transport properties of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 single crystals in the mixed state 6 (a) 1 U ( 1 3 K) 1.1 (c) 2. q z = z =.67(5).57(3) 1 H (T) H (T) U ( 1 3 K) (d) 2.2 q z =.6 z =.6.17(5).38(5) 1 H (T) H (T) Figure 3. Field dependence of (a, b) U and (c, d) q obtained from the fits of the resistivity in the TAFF region using Eq. (5) for LiFeSe-122 and LiFeTeSe-122, respectively. The solid lines in (a) and are fits using Eq. (6) and the values beside the lines are fitted γ. is.17(5) for H ab and.38(5) for H c in LiFeTeSe-122 (Fig. 3). The weak field dependence of U implies that single-vortex pinning should be dominant in LiFeTeSe In contrast, the faster decreases of U (µ H) in LiFeSe-122 imply that the collective flux pinning maybe dominate[32]. The better fits using Eq. (5) than Arrhenius relation can be partially ascribed to considering the temperature dependence of prefactor ρ f. When the assumption of U T with q = 1 is satisfied, the Arrhenius relation is valid and it can well describe the TAFF behavior, as in LiFeSe-122 for H ab. But because the obtained U of LiFeSe-122 for H c and LiFeTeSe-122 for both field directions are much smaller than that of LiFeSe-122 for H ab, the assumption of temperature-independent ρ f, i. e. U T, might be improper and the Arrhenius relation will become invalid. As shown in Fig. 3(c) and (d), the q for both crystals are almost independent of the field strength for both directions. When H c, the values of q for both samples are about 2, but it changes from 1 in LiFeSe-122 to 2 in LiFeTeSe-122 when H ab. The q = 1 for H ab in LiFeSe-122 is consistent with the good linear behavior of lnρ(t,µ H) vs. 1/T shown in Fig. 1(a). In contrast, the value of q = 2 should be another reason leading to the deviation of experimental curves from Arrhenius relation. According to the VG model[33, 34], the linear resistivity disappears as a power law ρ = ρ T T g 1 s (7) close to the glass transition temperature, T g. In Eq. (7), ρ is identified as a characteristic resistivity and should in some way be related to the normal state resistivity

7 Transport properties of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 single crystals in the mixed state 7 (a) 1 z = (c) 1 z =.6 (dln /dt) -1 T * 5 5 T * T g H (T)= T g z = (d) z =.6 (dln /dt) -1 1 T g T * (dln /dt) -1 (dln /dt) -1 1 T g T * Figure 4. (dlnρ/dt) 1 vs. T at µ H = 8 and 14 T of (a, b) LiFeSe-122 and (c, d) LiFeTeSe-122 for H ab and H c, respectively. (a) 15 H c2 H * H g (c) 15 H (T) (d) H (T) 1 5 z = pinned vortex liquid vortex glass z = unpinned vortex liquid normal state H (T) H (T) z = z = Figure 5. The vortex phase diagrams (µ H T) of (a, b) LiFeSe-122 and (c, d) LiFeTeSe-122 for H ab and H c, respectively. µ H c2 is the upper critical field defined from the 9% ρ n ; µ H represents the transition from unpinned to strongly pinned vortex region and µ H g represents the transition from vortex liquid to glass region. and s is the glass critical exponent. T g can be extracted by applying the relation (dlnρ/dt) 1 (T T g ) to the resistive tail, and T is defined as the temperature deviating from the straight line, as shown in Fig. 4. Based on the values of T g and T extracted from VG model and µ H c2 defined from the 9% ρ n, we plot the vortex phase diagrams of LiFeSe-122 and LiFeTeSe-122 (Fig. 5). The different regions indicate different flux-pinning mechanism. The region below the µ H g line is a vortex glass phase and the region above µ H c2 is the normal

8 Transport properties of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 single crystals in the mixed state 8 (a).2 (m cm) xx.1 xy ( cm) K -6 2 K 23 K 25 K 3 K K 4 25 K 5 K H (T) Figure 6. Field dependence of (a) longitudinal resistivity ρ xx (µ H) and Hall resistivity ρ xy (µ H) of LiFeSe-122 at various temperatures for H c. state, while the region between µ H g and µ H c2 is the vortex liquid phase. Besides, µ H represents the divider of unpinned and strongly pinned regions of vortex liquid. Obviously, the vortex liquid region is broad in these crystals, especially in LiFeSe-122 with H c, suggesting a weak vortex-pinning ability even at relative low field in this system when T approaching T c. Because the resistive transition curve ρ(t,µ H) not only shifts but also broadens when field increases, it suggests that thermal fluctuation is essential to vortex dynamics in these materials. The strength of thermal fluctuations can be quantified by the Gi = (γ 2 /2)[(µ k B T c )/(4πB 2 c ()ξ3 ab ())]2 [37], where γ is the electronic mass anisotropy, B c () = Φ /(2 2πλ ab ()ξ ab ()) is the thermodynamic critical field, Φ is the flux quantum, ξ ab () and λ ab () is the coherence length and penetration depth at K for H ab. With γ = ξ ab /ξ c, we obtain Gi = [(2πµ K B T c λ 2 ab())/(φ ξ c ())] 2 /2, where ξ c () is the coherence length at K for H c. Assumed λ ab ()= 2 nm same as (Li, Fe)OHFeSe[11], the value of Gi is and for LiFeSe-122 and LiFeTeSe- 122 with ξ c () =.27 and.9 nm derived from µ H c2 (). The value of Gi for LiFeSe-122 is comparable with (Li, Fe)OHFeSe ( ) and YBCO (1 2 ), but much larger than that of LiFeTeSe-122[11, 37]. It is consistent with the more obvious broadening of ρ(t,µ H) in LiFeSe-122, especially for H c. The remarkable thermally-activated vortex dynamics in LiFeSe-122 also lead to the wider region of vortex glass/liquid states in vortex phase diagrams when compared to LiFeTeSe-122. The Hall effect in the mixed state is another measurement to study vortex dynamics. Fig. 6(a) shows the field dependence of longitudinal resistivity ρ xx (µ H) of LiFeSe- 122 in the temperature range of 17 K - 5 K when H c. With increasing field,

9 Transport properties of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 single crystals in the mixed state 9 xy ( cm) X 2.5X (m cm) xx 1.2X 17K 2K 23K 25K 3K 35K 4K Figure 7. ρ xy vs. ρatvarioustemperaturesforlifese-122. Thesolidlinesarefitting results using the scaling behavior ρ xy (µ H) = Aρ xx (µ H) β. Inset: temperature dependence of β(t). the superconductivity is suppressed gradually and the transition of ρ(µ H) shifts to lower magnetic fields when temperature increases. ρ xx (µ H) shows a weak positive magnetoresistance at the normal state (T = 5 K), consistent with previous results[19]. As shown in Fig. 6, at the mixed state, the Hall resistivity ρ xy (µ H) at low field is zero and becomes negative with increased absolute values at the high fields. The highfield values of ρ xy (µ H) gradually reach that in the normal state when temperature is slightly higher than T c. The sign of the Hall resistivity is negative, indicating that the electron type carriers dominate in the mixed state as well as in the normal state, consistent with the electron doping in LiFeSe-122. Moreover, there is no sign reversal of ρ xy (µ H) in the mixed state, which is often observed in cuprates[38, 39]. Furthermore, there is a scaling behavior of ρ xy (µ H) = Aρ xx (µ H) β in LiFeSe-122 (Fig. 7) and the fitted values of β are close to 2 in the whole measuring temperature range. Different values of β have been observed in IBSCs and cuprate SCs, such as Fe(Te, S) (β =.9-1.)[26], Ba(Fe.9 Co.1 )As 2 (β = 2.(2))[4], and Bi 2 Sr 2 CaCu 2 O y (β = 2. ±.1)[41]. A number of theories have been proposed to explain the value of β. For example, when considering the effect of pinning on the Hall resistivity, Vinokur et al proposed a phenomenological model where β = 2. in the TAFF region[42]. And a unified theory for the Hall effect including both the pinning effect and thermal fluctuations was developed by Wang, Dong, and Ting[43, 44]. They explained scaling behavior by taking into account the backflow current effect on flux motion due to pinning and pointed out β changing from 2 to 1.5 for increasing pinning force. The exponent β 2. and the absence of sign change of ρ xy (µ H) for LiFeSe-122 single crystals imply that the strength of pinning force is relatively weak when H c and T T c, consistent with the low TAE at same field direction shown above.

10 Transport properties of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 single crystals in the mixed state 1 4. Conclusion In summary, the field-induced resistive broadenings of superconducting transitions and resistive tail behaviors have been observed in Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 single crystals. The detailed analysis shows that these results are associated with the TAFF and VG behaviors. Moreover, the thermal fluctuation plays an important role in the vortex dynamics ofli x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 singlecrystals, especially forlifese-122, resulting in broad vortex liquid regions in the vortex phase diagrams. The FFHE with β 2 indicates the relatively weak vortex pinning force in LiFeSe-122 when magnetic field is along the c-axis and T approaches T c. Acknowledgments This work was supported by the Ministry of Science and Technology of China (216YFA354), the National Natural Science Foundation of China (No ), the Fundamental Research Funds for the Central Universities, and the Research Funds of Renmin University of China (RUC) (15XNLF6, 15XNLQ7). References [1] Kamihara Y, Watanabe T, Hirano M and Hosono H 28 J. Am. Chem. Soc [2] Gurevich A 211 Rep. Prog. Phys [3] Hsu F C, Luo J Y, Yeh K W, Chen T K, Huang T W, Wu P M, Lee Y C, Huang Y L, Chu Y Y, Yan D C and Wu M K 28 Proc. Natl. Acad. Sci. USA [4] Vedeneev S I, Piot B A, Maude D K and Sadakov A V 213 Phys. Rev. B [5] Hunte F, Jaroszynski J, Gurevich A, Larbalestier D C, Jin R, Sefat A S, McGuire M A, Sales B C, Christen D K and Mandrus D 28 Nature (London) [6] Guo J G, Jin S F, Wang G, Wang S C, Zhu K X, Zhou T T, He M and Chen X L 21 Phys. Rev. B (R) [7] Chen F, Xu M, Ge Q Q, Zhang Y, Ye Z R, Yang L X, Jiang J, Xie B P, Che R C, Zhang M, Wang A F, Chen X H, Shen D W, Hu J P and Feng D L 211 Phys. Rev. X [8] Lu X F, Wang N Z, Wu H, Wu Y P, Zhao D, Zeng X Z, Luo X G, Wu T, Bao W, Zhang G H, Huang F Q, Huang Q Z and Chen X H 215 Nat. Mater [9] Dong X, Jin K, Yuan D, Zhou H, Yuan J, Huang Y, Hua W, Sun J, Zheng P, Hu W, Mao Y, Ma M, Zhang G, Zhou F and Zhao Z 215 Phys. Rev. B [1] Wang Z, Yuan J, Wosnitza J, Zhou H, Huang Y, Jin K, Zhou F, Dong X and Zhao Z 217 J. Phys. : Condens. Matter [11] Yi X, Wang C, Tang Q, Peng T, Qiu Y, Xu J, Sun H, Luo Y and Yu B 216 Supercond. Sci. Technol [12] Lin H, Xing J, Zhu X Y, Yang H and Wen H H 216 Sci. Chi. PMA [13] Ying T P, Chen X L, Wang G, Jin S F, Zhou T T, Lai X F, Zhang H and Wang W Y 212 Sci. Rep [14] Scheidt E-W, Hathwar V R, Schmitz D, Dunbar A, Scherer W, Mayr F, Tsurkan V, Deisenhofer J and Loidl A 212 Eur. Phys. J. B [15] Ying T P, Chen X L, Wang G, Jin S F, Lai X F, Zhou T T, Zhang H, Shen S J and Wang W Y 213 J. Am. Chem. Soc [16] Burrard-Lucas M, Free D G, Sedlmaier S J, Wright J D, Cassidy S J, Hara Y, Corkett A J, Lancaster T, Baker P J, Blundell S J and Clarke S J 213 Nat. Mater

11 Transport properties of Li x (NH 3 ) y Fe 2 (Te z Se 1 z ) 2 single crystals in the mixed state 11 [17] Sedlmaier S J, Cassidy S J, Morris R G, Drakopoulos M, Reinhard C, Moorhouse S J, O Hare D, Manuel P, Khalyyavin D and Clarke S J 214 J. Am. Chem. Soc [18] Wang S H, Sun S S and Lei H C 217 arxiv: [19] Sun S S, Wang S H, Yu R and Lei H C 217 arxiv: [2] Lee H-S, Bartkowiak M, Park J-H, Lee J-Y, Kim J-Y, Sung N-H, Cho B K, Jung C-U, Kim J S and Lee H-J 29 Phys. Rev. B [21] Fendrich J A, Kwok W K, Giapintzakis J, van der Beck C J, Vinokur V M, Fleshier S, Welp U, Viswanathan H K and Crabtree G W 1995 Phys. Rev. Lett [22] Palstra T T M, Batlogg B, Schneemeyer L F and Waszczak J V 1988 Phys. Rev. Lett [23] Blatter G, Feigel man M V, Geshkenbein V B, Larkin A I and Vinokur V M 1994 Rev. Mod. Phys [24] Palstra T T M, Batlogg B, van Dover R B, Schneemeyer I F and Waszczak J V 199 Phys. Rev. B [25] Brandt E H 1995 Rep. Prog. Phys [26] Lei H C, Hu R W, Choi E S and Petrovic C 21 Phys. Rev. B [27] Zhang Y Z, Wen H H and Wang Z 26 Phys. Rev. B [28] Lei H C, Hu R W and Petrovic C 211 Phys. Rev. B [29] Shahbazi M, Wang X L, Ghorbani S R, Dou S X and Lin C T 215 Physica C [3] Wu Z F, Tao J, Xu X B, Qiu L, Yang S G and Wang Z H 216 Physica C [31] Thompson J R, Sorge K D, Cantoni C, Kerchner H R, Christen D K and Paranthaman M 25 Supercond. Sci. Technol [32] Yeshurun Y and Malozemoff A P 1988 Phys. Rev. Lett [33] Espinosa-Arronte B and Andersson M 25 Phys. Rev. B [34] Andersson M, Rydh A and Rapp Ö 21 Phys. Rev. B [35] Ghorbani S R, Wang X L, Shabazi M, Dou S X, Choi K Y and Lin C T 212 Appl. Phys. Lett [36] Su T S, Yin Y W, Teng M L, Zhang M J and Li X G 214 Mater. Res. Express [37] Eley S, Miura M, Maiorov B and Civale L 217 Nat. Mater [38] Cagigal M, Fontcuberta J, Crusellas M A, Vicent J L and Piñol S 1995 Physica C [39] HagenSJ,SmithAW,RajeswariM,PengJL,LiZY,GreeneRL,MaoSN,XiXX,Bhattacharya S, Li Q and Lobb C J 1993 Phys. Rev. B [4] Wang L M, Sou U-C, Yang H C, Chang L J, Cheng C-M, Tsuei K-D, Su Y, Wolf Th and Adelmann P 211 Phys. Rev. B [41] Samoilov A V 1993 Phys. Rev. Lett [42] Vinokur V M, Geshkenbein V B, Feigel man M V and Blatter G 1993 Phys. Rev. Lett [43] Wang Z D and Ting C S 1991 Phys. Rev. Lett [44] Wang Z D, Dong J and Ting C S 1994 Phys. Rev. Lett