Temperature and time scaling of the peak-effect vortex configuration in FeTe{}_{0.7}Se{}_{0.3}

Temperature and time scaling of the peak-effect vortex configuration in FeTeSe

Marco Bonura Marco.Bonura@unige.ch Département de Physique Appliquée - Université de Genève, rue de l’École de Médecine 20, 1211 Genève, Switzerland    Enrico Giannini Département de Physique de la Matière Condensée, Université de Genève, Quai Ernest-Ansermet 24, 1211 Genève, Switzerland    Romain Viennois Département de Physique de la Matière Condensée, Université de Genève, Quai Ernest-Ansermet 24, 1211 Genève, Switzerland Institut Charles Gerhardt Montpellier UMR 5253, Université Montpellier II and CNRS, Place Eugène Bataillon, 34095 Montpellier, France    Carmine Senatore Département de Physique Appliquée - Université de Genève, rue de l’École de Médecine 20, 1211 Genève, Switzerland Département de Physique de la Matière Condensée, Université de Genève, Quai Ernest-Ansermet 24, 1211 Genève, Switzerland

An extensive study of the magnetic properties of FeTeSe crystals in the superconducting state is presented. We show that weak collective pinning, originating from spatial variations of the charge carrier mean free path ( pinning), rules in this superconductor. Our results are compatible with the nanoscale phase separation observed on this compound and indicate that in spite of the chemical inhomogeneity spatial fluctuations of the critical temperature are not important for pinning. A power law dependence of the magnetization vs time, generally interpreted as signature of single vortex creep regime, is observed in magnetic fields up to 8 T. For magnetic fields applied along the axis of the crystal the magnetization curves exhibit a clear peak effect whose position shifts when varying the temperature, following the same dependence as observed in YBaCuO. The time and temperature dependence of the peak position has been investigated. We observe that the occurrence of the peak at a given magnetic field determines a specific vortex configuration that is independent on the temperature. This result indicates that the influence of the temperature on the vortex-vortex and vortex-defect interactions leading to the peak effect in FeTeSe is negligible in the explored range of temperatures.

74.70.Xa, 74.25.Ha, 74.25.Wx

I Introduction

The study of the vortex properties in type-II superconductors is of extreme interest both for investigating the basic physics of the superconductivity and for evaluating the quality of the materials in view of practical applications. Many aspects of this subject are still topical, especially with regard to the high temperature superconductors (HTS) and the recently discovered iron-based superconductors. To date, five families of Fe-based superconductors have been discovered: REOFeAs, (”1111”, RE=rare earth),Kamihara () AFeAs (”122”, A=alkaline earth), Rotter () XFeAs (”111”, X=Li; Na),Wang (); Parker () Fe(Se,Ch) (”11”, Ch=S, Te) Hsu (); Yeh () and the most recently discovered ”21311” family of SrMOFePn (M=Sc, V, Cr and Pn=pnictogen).Ogino () Among these families, iron chalcogenides are considered of particular interest because of their simple crystal structure consisting of Fe ions tetrahedrally coordinated by Se and Te arranged in layers stacked along the c-axis, without any other interlayer cations, as occurs in the pnictides. For this reason, iron chalcogenides are generally considered an ideal candidate for understanding some open issues of high-temperature superconductivity.

One of the most intriguing phenomena observed in the study of the vortex properties in type-II superconductors is the so called peak (or fishtail) effect. A second peak in the magnetization curve () of low- superconductors has been observed as early as in the sixtiesDe Sorbo () and a first theoretical explanation has been suggested by Pippard in 1969.Pippard () The peak has been found in proximity of in NbSn, Lortz () CeRu, Banerjee () NbSe, Banerjee () VSiIsino () and MgB. Pissas (); Angst ()
The peak effect has been clearly observed also in the cuprate superconductors but its features vary from one compound to the other. In YBaCuO crystals, the peak is broad and temperature dependent; its origin has been attributed to the presence of inhomogeneities in the oxygen contentKlein (); Wen (); Tobos () and associated to different mechanisms.Krusin-Elbaum (); Civale (); Abulafia (); Kokkaliaris (); Rosenstein () In the more anisotropic BiSrCaCuO the peak in the magnetization loop has usually been associated to a 3D to 2D transition of the vortex structure. Yang1 (); Lee () However, significant differences in the second-peak features have been observed within this family: in the highly anisotropic Bi-2212 and Bi-2223 the peak is sharp and temperature independent Khaykovivh (); Piriou1 (); Piriou2 (); Tamegai (); Avraham () while in Bi-2201 crystals it exhibits a strong temperature dependence. Piriou3 () Regarding the iron-based superconductors, the peak effect has been observed in FeAs-1111, Moore (); Senatore (); Proz-NJP () in FeAs-122 Yang2 (); Prozorov (); Proz-PhisicaC (); Shen (); Salem () and in FeTeSe,Das () even if a clear understanding of the physical origin of the phenomenon in this family has not been achieved, yet. However, since it has been shown by small angle neutron scattering,Eskildsen (); Inosov (); Eskildsen_2 () Bitter decorationEskildsen () and magnetic force microscopyInosov () that an ordered vortex lattice is not present in these compounds, the origin of the peak effect due to a phase transition of the vortex lattice seems to be excluded.

Beyond the investigation of the magnetization curves, an useful approach for investigating the vortex dynamics is the study of the relaxation processes of the magnetization. The critical state of the vortex lattice, which determines the hysteresis of the magnetization in type-II superconductors, is a metastable state. It follows that vortices tend to hop out of their pinning-potential well in order to reach the configuration of absolute minimum energy. Such motion usually arises from thermal activation, but it can also arise from quantum tunneling (at low temperatures) or can be stimulated by external perturbations, such as microwave shaking of the vortex lattice. Mikitik () Magnetic relaxation processes have been observed in various low-temperature superconductors. Yeshurun () However, the subject has become of even greater interest after the discovery of high- superconductors, because of the higher operating temperatures and of the small activation energies related to the short coherence length and large anisotropy. The concept of thermally-induced hopping of the flux lines has been first treated by Anderson and Kim. Anderson (); AndersonKim () In the framework of their model a logarithmic dependence of the magnetization () on the time () is expected. This behavior has been verified in various superconductors, both low and high ; however, in many experiments, deviations from the logarithmic dependence of have been observed, indicating a failure of the approximations that bring this expectation. Yeshurun (); Zola () In the framework of the collective pinning theory a more complex expression for the dependence of the potential energy barrier height (the so-called “interpolation formula”) is proposed: where is a parameter varying as a function of the vortex-vortex and vortex-pinning center interaction. Feigelman (); Blatter () For example, in three dimensional systems when the creep is dominated by the motion of individual flux lines, in the case of collective creep of small bundles and when the bundle size is much larger than the Larkin correlated volume Blatter (). In 1991, Vinokur, Feigel’man and Geshkenbein proposed a theoretical model for the thermally activated flux creep, assuming a logarithmic dependence of the activation energy on the critical current, . Vinokur () This dependence of over is a good approximation for the creep activation barrier in the single vortex creep regime (limit ) and provides a proper fit of the experimental dependence observed in different superconductors. Zeldov (); Maley (); McHenry (); Zhang () The divergence of when can be understood in the context of the collective pinning considering that in the vortex-glass state vortex motion is possible only in the presence of a current. Yeshurun () In the case of this logarithmic dependence of the activation energy on the critical current, a power-law time dependence of the magnetization is expected. Vinokur ()

In this paper we study the magnetic properties of FeTeSe crystals in the superconducting state by investigating the dependence of the magnetic moment () on the applied magnetic field () the temperature () and the relaxation time (). In Sec. II we report details on the examined samples and on the experimental techniques used. The experimental results are shown and discussed in two different sections. In Sec. III.1, dedicated to the vortex pinning properties, we analyze the curves and deduce the pinning mechanism in the investigated crystals. The magnetic relaxation processes are studied in Sec. III.2. Conclusions are reported in Sec. IV.

Ii Sample and Experimental

The FeTeSe crystals investigated were grown by a modified Bridgman-Stockbarger method, starting from a nominal Fe:(Te,Se) ratio of 0.9:1. Details of the procedure used for preparing the samples are reported in Ref. Viennois (). The refined composition of the samples is FeTeSe, indicating a very low iron excess.Viennois () It has been shown that reducing the Fe excess in FeTeSe favors the occurrence of superconductivity and weakens the antiferromagnetic order. Viennois (); Liu () However, according to the phase diagram of the compound, Okamoto () a little excess of Fe is needed for stabilizing the structure. It is also worth noting that, due to Se-doping for Te, less Fe is allowed to occupy the additional site, since both the effects of reducing and increasing result in shrinking and re-shaping the FeTe tetrahedra; so, the lower is the Se content, the more difficult is to limit the Fe excess.Viennois ()

Two different crystals, weighting mg (Sample A) and mg (Sample B), have been investigated in the present study. The susceptibility curve of the two samples has been deduced from a zero field cooled measurement performed at 10 Oe by a Superconducting Quantum Interference Device (SQUID) with the field applied along the axis. Both curves have been corrected for the demagnetization effects according to the formula reported in Ref.Stoner (). The results are shown in Fig.1: the curve of Sample A exhibits a superconducting transition with  K,  K whilst for Sample B it results  K,  K. Both the samples exhibit bulk superconductivity with .

Figure 1: Temperature dependence of the DC volume susceptibility (demagnetization corrected) for the investigated samples.

The samples have been characterized also by x-ray diffraction (XRD) both in a powder diffractometer (on manually ground crystals) and in a 4-circle diffractometer. The rocking curves indicate the presence of a distribution of the c-axis crystallites with a full width at half maximum (FWHM) of about for the Sample A and about for Sample B. Our results indicate that, whereas the sample B is of high crystalline quality, the sample A is more likely to be composed of few single-crystalline domains well aligned along the crystallographic c-axis.

The magnetization curves () of the FeTeSe samples were acquired by means of a SQUID magnetometer and a VSM (Vibrating Sample Magnetometer). The maximum applied magnetic field was 4.5 T for the SQUID and 8.8 T for the VSM. Measurements were performed in the range of temperatures from 4.25 K to , with the magnetic field applied both parallel and perpendicularly to the axis. The VSM measurements were executed with different values of the magnetic field sweep rate () namely 0.05 T/min, 1 T/min and 2 T/min. Furthermore, by means of the VSM, we studied the relaxation of magnetization over periods of time up to 6000 sec, for fixed values of temperature and magnetic field. In order to perform measurements in the trapped flux configuration, the relaxation measurements were performed with the following procedure: i) the sample is zero-field cooled; ii) the magnetic field is increased from zero to 8.8 T with = 2T/min; iii) the field is decreased from 8.8 T to 8 T with = 0.1 T/min; iv) the magnetic relaxation at 8 T is recorded for 6000 sec; iv) the field is decreased by 1 T with = 0.1 T/min; v) the magnetic relaxation is recorded for 6000 sec. Steps iv) and v) are repeated for registering the relaxation down to 0 T, in steps of 1 T.

Iii Results and Discussion

iii.1 Vortex pinning properties

Fig. 2 shows the magnetic-moment curves of the Sample A acquired with the VSM at different temperatures, for (a) and (b), for a sweep rate of 2 T/min. Increasing the temperature, the pinning becomes weaker and consequently the width of the hysteresis loop decreases. A second peak in the () curve is clearly observed for , up to temperatures close to . A second peak is also present for but less pronounced and detectable only for temperatures higher than 5 K. In the inset of Fig. 2 the results obtained for the Sample B at  K and  K, are reported for the two different orientations of the field. In this case, a second peak in the () curve is not observed for , whereas it is clearly evident for , in analogy to what observed in YBaCuO crystalsKlein () and in other Fe-based superconductors.Pramanik (); Prozorov () Since Sample A presents a wider c-axis distribution of single-crystalline domains than Sample B, we argue that the weak peak effect revealed for in this sample is due to the slight misalignment of the crystallites forming the sample. We would also like to specify that the experimental curves of Sample B for are qualitatively similar with those of Sample A at comparable temperatures. This suggests that a misalignment of the layers forming the crystal is crucial only for magnetic fields applied perpendicularly to the axis of the sample. The curves obtained for exhibit a paramagnetic-like background arising from the experimental setup. This background is not observed in the curves measured for in the SQUID. In any case, the background does not alter the results obtained in the following, which rely only on the separation between the positive and negative branches of the magnetic-moment loop.

Figure 2: Magnetic moment vs. external field curves obtained at different temperatures in Sample A for (a) and (b). The insets show the results obtained in Sample B at  K and  K.

From the () curves measured for at different , we determined the temperature dependence of the magnetic field value corresponding to the second peak in magnetization (). The results for samples A and B are shown in Fig. 3. It follows that for the FeTeSe superconductor the position of the second-peak shifts towards lower fields monotonically on increasing the temperature, in analogy to what observed in the YBaCuO superconductor. Zhukov () For both samples investigated, the peak becomes undetectable at  K. The continuous lines in Fig. 3 are the best fit curves obtained supposing the same -dependence of as observed in YBaCuO:Klein () , where is a constant and is the temperature at which the peak is undetectable. The best fit parameters (with the relative statistical errors) are:  T,  K for Sample A and  T,  K for Sample B. This remarkable similarity between the results obtained in the YBaCuO and in the FeTeSe superconductors, might suggest that the peak effect in the two systems has an analogous origin.

Figure 3: Temperature dependence of the magnetic field value at which the peak effect occurs, for the two measured samples. Continuous lines are the best fit curves obtained supposing the same dependence observed in YBaCuO.

The magnetic field dependence of the critical current density () has been extracted from the () curves, for different values of the temperature, using the Bean critical state formulae.Bean1 (); Bean2 () For slab samples in a perpendicular magnetic field: , where () is the separation between the two branches of the magnetic-moment loop, and are the length and the width of the sample () and is the thickness. The () curves obtained for the Sample A at different temperatures are presented in Fig. 4; in the inset of the same figure the () curve at  K of the Sample B is shown. For both samples investigated the values at 4.25 K are of the order of  A/cm. This value is smaller than what measured in FeTeSe crystals with higher values than our samples.Taen (); Liu2 () However, the values become comparable if related to the same reduced temperature .Liu2 () The fact that the values obtained for Sample A are very similar to those obtained for the high-quality single crystal Sample B, where no weak links are present, indicates that weak links effects are not important for Sample A either.

The critical current density is always limited by the depairing current density , where and are the thermodynamic critical field and the GL penetration depth, respectively. Griessen () Important information can be deduced from the ratio . For low- superconductors pinning is usually strong () resulting from the interaction of vortices with extended defects such as, e.g., precipitates or grain boundaries.Blatter () On the contrary, for the cuprate family pinning is usually weak () normally arising from point defects, e.g., oxygen vacancies.Blatter () For the FeTeSe we estimate indicating that also in the case of the Fe-based 11 family pinning is weak. The microscopic origin of pinning in FeTeSe will be discussed in the following.

Figure 4: Field dependence of the critical current density at different temperatures and for , in Sample A. The inset shows the results obtained at  K and for for the Sample B.

As already mentioned, the presence of the peak effect in the magnetization curve has been widely documented in the literature, both in low and high superconductors. Its origin has been associated to different processes, depending on the particular system investigated. Pippard (); Krusin-Elbaum (); Rosenstein (); Avraham () In order to shed light on the mechanisms that rule pinning in the FeTeSe superconductor, we have investigated the magnetic-field dependence of the pinning force density, . It has been shown for a large variety of low- and high- superconductors that the curves of vs. obtained at different temperatures may be scaled into a unique curve if they are plotted as a function of the reduced field, .Fietz (); Dew-Hughes (); Koblischka1 () The empirical formula that accounts for the scaling is:


where is a proportionality constant, and are two parameters whose values depend on the origin of the pinning mechanism. Dew-Hughes () The different contributions to flux-pinning are usually catalogued into two main categories: i) (or normal) pinning, arising from spatial variations in the charge carrier mean free path (); ii) (or ) pinning, associated with spatial variations of the Ginzburg parameter () due to fluctuations in the transition temperature . Dew-Hughes (); Blatter () In the last category fall also non-superconducting metallic particles whose dimension is smaller than the coherence length () since, in this case, the superconductivity is induced by the proximity effect. A classification is also made for pinning centers, as a function of the number of dimensions that are large with respect to the inter-vortex distance . Following the definition given by Dew-Hughes in Ref.Dew-Hughes (), in this paragraph we refer to point pins as regions whose dimensions in all directions are less than , line pins have one dimension larger than , grain- and twin- boundaries, which have two dimensions greater than , act as surface pins, volume pins have all dimensions large with respect to .Dew-Hughes ()

In the framework of the Dew-Hughes model,Dew-Hughes () different values for and in Eq. 1 are expected, as a function of the specific pinning mechanism involved. Correspondingly, the theoretical vs. curves present a maximum at different values. In the case of pinning, the maximum is expected at (, ) for point pins and at (, ) for surface pins, such as grain boundaries (the same dependence has been predicted for shear-breaking in the case of a set of planar pins Kramer ()); no maximum is expected in the case of volume pinning, being and . The maximum of the curve is expected at higher values in the case of pinning; in particular it occurs at (, ) for point pins, at (, ) for surface pins and at (, ) for volume pins. Therefore, important information on the physical origin of the pinning mechanisms can be achieved by analyzing the scaled curves.

The comparison of the predicted pinning functions with the experiments requires an estimation of , defined as the -value at which . Starting from the experimental curves, is usually determined as the extrapolated zero value in the so-called Kramer plot,Kramer () where is plotted as a function of . This procedure has been successfully used in the case of wires, polycrystalline samples and crystals.Flukiger (); Ekin (); Ekin2 (); Putti () However, if pinning is ruled only by defects whose dimensions are smaller than the intervortex distance, it is expected .Dew-Hughes () In this case, may be determined by doing a linear extrapolation down to zero of the vs. curve.Dew-Hughes (); Montgomery (); Shapira () In Fig. 5 (a) and (b) we report the Kramer plots along with the vs. curves, obtained at two different temperatures, for Sample A and Sample B, respectively. We observe that the Kramer plot presents a wide linear behavior for both samples. On the other hand, the curve exhibits a linear behavior in a wide range of magnetic fields, only for Sample B. For this reason, we have extracted from the Kramer plot for Sample A, whereas for Sample B we cannot discriminate a-priori which is the most correct procedure. However, a more conclusive result can be achieved if one considers simultaneously both the vs. dependence and the scaled curves. If one hypothesizes that pinning is ruled by point pins and consequently extracts from the vs. curve then the corresponding curve should exhibit a maximum at .

Figure 5: Kramer plot (right axis) and vs. (left axis) curves obtained at  K and  K for sample A (a) and sample B (b), as described in the text.

Fig. 6 (a) and (b) present the normalized pinning-force density for samples A and B at different temperatures as a function of the reduced field based on values deduced from the Kramer plots, whereas the inset of Fig. 6 (b) is based on values determined from a linear vs. dependence. For the sake of clearness, the superscript of and in the figure indicate respectively the values extracted from the Kramer plot or hypothesizing .

Figure 6: Normalized pinning force density curves as a function of the reduced field, , for sample A (a) and B (b). indicates that the irreversibilty field has been deduced by a linear extrapolation in the Kramer Plot. The inset in Fig. 6 (b) shows the results obtained in sample B, scaled by using the irreversibility field, , deduced by the vs. curves

The pinning force curves of Sample A, obtained at different temperatures, scale well and present a maximum at while, for Sample B, the maximum occurs at if one considers or at if . These results indicate that, for the single crystal Sample B, the pinning mechanisms is ruled by defects whose dimensions are smaller that the intervortex distance in the investigated field range and as a consequence, the most appropriate procedure for determining is from the vs. curve. On the contrary, for Sample A the position of the maximum in the curve suggests that a contribution to pinning coming from surface pinning cannot be excluded. From the scaled pinning-force density curves we can also deduce information on the type of pinning centers. In particular, the fact that the maximum in the curve occurs for values of the reduced field is an indication that the pinning centers in FeTeSe are of type. It is worth noting that this conclusion does not depend on the specific procedure used for determining H. In fact, in the case of pinning contribution coming from -type pins, the maximum is expected to occur at , Koblischka () which is far away from what observed in our samples.

The analysis of the scaled pinning force curves in the frame of the Dew-Hughes model revealed to give important information about pinning in low-, cuprates as well as Fe-based superconductors.Koblischka1 (); Koblischka (); Tobos (); Yang2 (); Tarantini (). However since in our samples pinning is weak () and thermal fluctuations are large (Ginzburg number ), more conclusive results about pinning origin in FeTeSe are expected to be achieved by analyzing the experimental results in the framework of the collective pinning theory.Blatter () Following the theoretical approach proposed by Griessen et al, Griessen () in the case of -type weak pinning in the single vortex regime it is expected that the critical-current density variation with respect to the reduced temperature () is described by the following expression:Griessen ()


while, for pinning, it is:


In Fig. 7, we plot the normalized data obtained at  T and  T for , along with the theoretical curves expected within the scenario of and pinning. The values have been extracted from the () curves obtained at various temperatures. Analogous results have been obtained for sample B. A remarkably good agreement between the experimental results and the pinning theoretical curve is obtained, confirming that pinning in the FeTeSe samples originates form spatial variation of the mean free path.

The weak collective pinning theory considers that pinning originates from fluctuation in the density and force of defects whose dimension is smaller than the coherence length. In this case small values for the critical current density are expected ().Blatter () On the contrary, if pinning was mainly ruled by extended defects such as, e.g., precipitates, twinning or grain boundaries, or columnar defects such as dislocation lines, higher values for the ratio would be expected ().Blatter () A phase separation in FeTeSe single crystals with different and values has been experimentally observed by scanning tunneling microscopy and electron energy-loss spectroscopy. Hu (); He () In particular, fluctuations of the local Te concentration (even by 20 from the average local composition) in nanometric regions have been observed in FeTeSe crystals. Hu () These fluctuations could be the origin of the observed weak pinning.

In spite of the inhomogeneous chemical distribution, our results indicate that local fluctuations of the critical temperature are not important for pinning. This interesting result is in agreement with tunneling spectroscopy experiments performed at  K in optimal doped FeTeSe crystals, indicating that the compound is chemically inhomogeneous but electronically homogeneous.He ()

Figure 7: Normalized critical current density data, as a function of the reduced temperature (), obtained at  T and  T for . The continuous lines are the theoretical curves expected in the case of and pinning.

It is worth mentioning that a coexistence of superconducting and magnetic orders has been proposed for the FeTeSe system. In particular such coexistence has been inferred by Khasanov et al. based on muon-spin rotation experiments.Kasanov () On the other hand, Li et al. have ascribed the apparent coexistence of superconductivity and magnetism in the FeTeSe system to a phase separation in the real space.Li () As already mentioned, a phase separation has also been observed in FeTeSe crystals by STM.He () Other authors have assumed a spin-glass state to exist between the long range antiferromagnetic order (at low Se content) and the bulk superconducting state, at high Se content.Katayama () The low-field DC susceptibility measurements of our samples, reported in Fig.1, do not show any anomaly or feature attributable to the coexistence of superconductivity and magnetic order. In a previous article we have deeply investigated the effect of the actual chemical composition and its effect on the crystal chemistry and on the magnetic and superconducting phase diagram.Viennois () We have proved the combined effect of Fe-excess () and Se-substitutions () on the transition from a superconducting to an antiferromagnetic state and highlighted the importance of controlling the real composition in a three dimensional phase diagram. By keeping under control a low Fe-excess, we have obtained bulk superconductivity at Se-content even lower than .Viennois () In 2010, Bendele et al.Bendele2 () has confirmed our study on the excess Fe and merged data from previous publications in a three dimensional phase diagram, still claiming a coexistence of superconductivity and magnetism. On the basis of our experimental results and the contradicting results reported in the literature, it is not possible to discriminate whether such coexistence occurs intrinsically at the atomic level in FeTeSe or a phase separation and local composition fluctuations are responsible for the magnetic behavior observed by some authors. Further studies on very clean and chemically homogeneous samples, carried out with different experimental techniques, are needed in order to clarify this point.

At the end of our analysis on the pinning properties we would like to remind that both the Dew-Hughes and the Griessen models have been developed in the single-vortex pinning regime, i.e. neglecting inter-vortex interactions. The validity of this approximation in our case is confirmed by the magnetic relaxation study reported in the next section.

iii.2 Vortex dynamics properties

In this section we report a study of the dynamical properties of vortices in the FeTeSe superconductor performed by magnetic relaxation measurements. The analysis has been limited to Sample A, since the measured magnetic moment of Sample B was too small for achieving a good resolution in a wide range of times and temperatures. Fig. 8 shows the relaxation of the magnetization () normalized to its maximum value, obtained at  K for , plotted in a log-log graph. At all the fields investigated and for times greater than about 50 sec, a linear dependence of has been observed, indicating a power-law dependence of vs. . An analogous behavior has also been detected at higher temperatures and for . It has been shown that deviations from the expected curves at short times could be due to the magnetic field overshoot occurring when the external field ramp is stopped.Zola () This overshoot produces a shielded flux zone in proximity of the surface of the sample, which affects the initial relaxation process. The inset of Fig. 8 shows the values obtained at different values; large values of have already been observed in the FeTeSe system as well as in other Fe-based superconductors.Yadav (); Prozorov (); Yang-Sm () Since the value is related to the pinning potential energy barrier height, Blatter () the minimum observed in the () curve is associated to the presence of the peak effect in the () curves.

Figure 8: Plot of the magnetization vs. time (on a log-log scale) at  K for different values of the magnetic field. The inset shows the field dependence of the relaxation rate .

As already described in the introduction, in the framework of the Anderson and Kim theory, a linear dependence of on log() is expected; this result stems from two basic assumption: i) the pinning potential energy barrier height decreases linearly with the current density: ; ii) , which allows to hypothesize that the thermal-induced hopping rate is proportional to the Arrhenius factor . Supposing a linear dependence of on is only a first-order approximation whose validity has been demonstrated to fail many times. Yeshurun () As reported by Vinokur, Feigel’man and Geshkenbein, a power-law dependence of vs. is expected if a logarithmic dependence of the activation energy on the current density is supposed: . Vinokur () In this case, the following expressions for the time dependence of are predicted:


where is the thickness of the sample, is the position of the flux front and is the time at which the sample is fully penetrated. Vinokur (); Zhang ()

In order to verify that a logarithmic dependence of () is actually present for the FeTeSe crystal investigated, we used the method proposed by Maley et al.,Maley () which allows the determination of the () curve from the experimental data obtained at different temperatures. In their paper, the authors show that choosing a proper value for a time-independent constant () it is possible to deduce the dependence of (or equivalently the () dependence) by plotting dd vs. , where M is the magnetization at the equilibrium. This procedure is valid under the assumption that the temperature dependence of is weak and that the principal effect of increasing the temperature is to produce monotonically decreasing initial values of , which is usually valid for . Maley ()

Fig. 9 shows the curves obtained at  T, scaled considering A=28, along with the best fit curve deduced supposing the following dependence of on : ; the best fit parameters (with the relative statistical errors) are: K, emu/cm and emu/cm. Similar values for have been obtained in other Fe-based superconductors, indicating that the weakness of pinning is a general characteristic of these compounds.Shen () The negative value of is due to the diamagnetic contribution of the sample holder. A very good scaling is obtained for K, while at K a discrepancy between the experimental data and the theoretical curve starts to be present. As already mentioned in the introduction, the logarithmic dependence of on is a good approximation for the creep activation barrier in the single vortex creep regime. As a consequence, our results indicate that in the FeTeSe superconductor the vortex motion develops in the single-vortex pinning limit even in magnetic fields up to 8 T. This is in general unexpected. At high fields, the intervortex distance becomes small compared to the magnetic field penetration depth and thus one would expect that vortex-vortex interactions become important and pinning involves vortex bundles. However, single vortex pinning up to 9 T has already been observed by Inosov et al. in BaFeCoAs single crystals by magnetization measurements, small-angle neutron scattering, and magnetic force microscopy.Inosov () Our results suggests that also in the FeTeSe system the vortex-defect interaction dominates up to high fields. This experimental evidence further confirms the validity of our study on the pinning properties as reported in Sec. III.1, carried out in the framework of the Dew-Hughes model where flux-lattice elasticity effects are neglected.

Figure 9: Magnetization dependence of the pinning potential energy barrier height calculated in the frame of the Maley model, scaling the data at different temperatures as described in the text.

In Sec III.1 we have shown that the curves obtained for are characterized by a second peak whose position changes with the temperature. In order to investigate the magnetic relaxation effects on the peak effect, we measured the field dependence of the magnetic moment, by means of the VSM, for different sweep rates () of the magnetic field. It has been demonstrated that the dependence of the hysteresis amplitude on the sweep rate contains basically the same information of the time dependence of the magnetization during relaxation, in particular, the higher is the sweep rate, the shorter is the effective “observing time” in the () curve.Schnack (); Jirsa () The results are shown in Fig. 10, along with the curve acquired by the SQUID and the values (full points) extracted from the curves, obtained at different fields, considering sec. Regarding the measurement performed with the SQUID, one can assume since any variation of during the acquisition can be neglected. Furthermore, since the acquisition time is of some seconds, one has also to consider that the measured value has already relaxed. A remarkable variation of the curves associated to different characteristic times is present, as a consequence of the relaxation in time of the magnetic moment. Furthermore, it is very interesting to note that the position of the second-peak moves towards low values during the relaxation, which indicates that the magnetic field value at which the second peak in the curve is obtained relaxes in time.

Figure 10: Time evolution of the second peak position investigated by magnetization measurements performed by the VSM at different sweep rates (), by the SQUID as well as by magnetic relaxation measurements.

In order to shed light on the mechanism that determines the temperature and time dependence of the peak position, , we collected the couples of values (, ) which identify its location in the () curves obtained at different temperatures by the SQUID and by the VSM and for different sweep rates, the results being shown in Fig. 11. The inset at the bottom right of the same figure presents the calculated induction-field values () corresponding to . It is very interesting to note that the pairs of values (, ) - or equivalently (, ) - obtained at different temperatures or sweep rates lie on the same curve and that, in some cases, the same pair of values (, ) has been obtained at two different temperatures but for a different “observing time” during relaxation. The continuous line shown in the figure is the best fit curve obtained by assuming . The best fit parameters, along with the corresponding statistical errors, are:  emu/T; .

Figure 11: Pairs of values (,) that identify the second peak position in the curve at different temperatures and characteristic times. Labels 1, 2, 3, indicate the measurements performed by the VSM with = 2T/min, 1T/min and 0.05 T/min, respectively; label 4 indicates the measurements performed by the SQUID. The continuous line is the best fit curve obtained as described in the text. The inset at the bottom right shows the couples of values (,) corresponding to the peak effect. The inset at the top left is an imaginary drawing of the vortex energy landscape, the line indicating the relaxation path of the vortex configuration associated to the peak effect.

Since the induction field is proportional to the number of vortices present in the SC and the magnetic moment identifies the critical current (and as a consequence the flux profile inside the sample), the pair of values (, ) identify a specific configuration of the vortex structure. The fact that all the experimental data shown in Fig. 11 fall in a unique curve suggests that there is a unique vortex configuration which determines the occurrence of the peak effect at a given value, whatever the temperature is in the range 4.2 K 8 K, corresponding to 0.4 0.75 . In particular, the same (,) pair can be found at a temperature and observing time or equivalently at a temperature at a time . If one defines the vortex energy landscape as the mapping of all the possible configurations of the vortices in the sample, determined by and , and the corresponding energy level (), the observed behavior could be justified by supposing that there is a unique path in the vortex energy landscape that describes the relaxation of the vortex configuration associated to the peak effect and that the temperature does not particularly affect the energy landscape in the range explored, at least in proximity of the path. The only effect of increasing the temperature is to allow the relaxation to start from a point in the path closer to the final equilibrium state. This result also indicates that the mechanism behind the second peak is related to a thermally driven process in the examined range of temperatures and fields. We are extending the present analysis of the peak effect properties to other superconductors, such as YBaCuO, and the results will be discussed in a forthcoming paper.

Iv Conclusions

We have investigated the magnetic properties of FeTeSe crystals in the superconducting state by magnetization and magnetic relaxation measurements. We have shown that pinning in FeTeSe originates form spatial variation of the mean free path and is most likely related to the nanoscale chemical phase separation observed in FeTeSe on a scale of  nm.Hu (); He () Very interestingly, our results confirm that even if chemically inhomogenous FeTeSe is electronically homogeneous since pinning is not ruled by spatial fluctuations of the critical temperature. From magnetic relaxation measurements we have obtained indications that vortex motion develops in the single vortex limit even in magnetic field as high as 8 T, in agreement to what observed in the 122 Fe-based compounds.Inosov ()On applying the magnetic field along the axis, a clear peak effect in the curves has been observed, up to temperatures near . The second-peak position varies with the temperature, following the same dependence as observed in YBaCuO.Klein () The relaxation of the vortex configuration that determines the peak effect has also been studied. We have found that the pairs of values (, ) that identify the second-peak position in the curves obtained at different temperatures and different relaxation times reconstruct a unique curve. This suggests that the vortex configuration that determines the peak effect at a particular value does not depend on the temperature in the range 0.4 0.75 . It follows that the temperature influence on the vortex-vortex and vortex-defect interactions in proximity of the peak effect is negligible.

V Acknowledgements

We acknowledge Prof. Thierry Giamarchi and Dr. Christophe Berthod for useful discussions, Dr. Florin Buta for critical reading of the paper, Damien Zurmuehle and Alexandre Ferreira for technical support.


  • (1) Y. Kamihara, T. Watanabe, M. Hirano, H. Hosono, J. Am. Chem. Soc. 130 (2008) 3296-3297.
  • (2) M. Rotter, M. Tegel, D. Johrendt, I. Schellenberg, W. Hermes, R. Pöttgen, Phys. Rev. B 78, 020503 (2008).
  • (3) X. C. Wang, Q. Q. Liu, Y. X. Lv, W. B. Gao, L. X. Yang, R. C. Yu, F. Y. Li, C. Q. Jin, Solid State Communications 148, 538 (2008).
  • (4) D. R. Parker, M. J. Pitcher, P. J. Baker, I. Franke, T. Lancaster, S. J. Blundell, S. J. Clarke, Chem. Commun. 16, 2189 (2009).
  • (5) Fong-Chi Hsu, Jiu-Yong Luo, Kuo-Wei Yeh, Ta-Kun Chen, Tzu-Wen Huang, Phillip M. Wu, Yong-Chi Lee, Yi-Lin Huang, Yan-Yi Chu, Der-Chung Yan, and Maw-Kuen Wu, Proc. Natl. Acad. Sci. USA 105, 14262 (2008).
  • (6) Kuo-Wei Yeh, Tzu-Wen Huang, Yi-lin Huang, Ta-Kun Chen, Fong-Chi Hsu, Phillip M. Wu, Yong-Chi Lee, Yan-Yi Chu, Chi-Lian Chen, Jiu-Yong Luo, Der-Chung Yan and Maw-Kuen Wu, Eur. Phys. Lett. 84, 37002 (2008).
  • (7) H. Ogino, Y. Matsumura, Y. Katsura, K. Ushiyama, S. Horii, K. Kishio, J. Shimoyama, Supercond. Sci. Technol. 22, 075008 (2009).
  • (8) W. De Sorbo, Rev. Mod. Phys. 36, 90 (1964).
  • (9) A. B. Pippard, Phil. Mag. 19, 217 (1969).
  • (10) R. Lortz, N. Musolino, Y. Wang, A. Junod, N. Toyota, Phys. Rev. B 75, 094503 (2007).
  • (11) S. S. Banerjee, N. G. Patil, S. Saha, S. Ramakrishnan, A. K. Grover, S. Bhattacharya, G. Ravikumar, P. K. Mishra, T. V. Chandrasekhar Rao, V. C. Sahni, M. J. Higgins, E. Yamamoto, Y. Haga, M. Hedo, Y. Inada, Y. Onuki, Phys. Rev. B 58, 995 (1998).
  • (12) M. Isino, T. Kobayashi, N. Toyota, T. Fukase, and Y. Muto, Phys. Rev. B 38, 4457 (1988).
  • (13) M. Pissas, S. Lee, A. Yamamoto, and S. Tajima, Phys. Rev. Lett. 89, 097002 (2002).
  • (14) M. Angst, R. Puzniak, A. Wisniewski, J. Jun, S. M. Kazakov, and J. Karpinski, Phys. Rev. B 67, 012502 (2003).
  • (15) L. Klein, E. R. Yacoby, Y. Yeshurun, A. Erb, G. Müller-Vogt, V. Breit, H. Wühl, Phys. Rev. B 49, 4403 (R) (1994).
  • (16) Hai-hu Wen and Zhong-xian Zhao, Appl. Phys. Lett. 68, 856 (1996).
  • (17) V. Tobos, L. M. Paulius, A. M. Petrean, S. Ferguson, J. W. Snyder, R. J. Olsson, W.-K. Kwok, and G. W. Crabtree , Appl. Phys. Lett. 78, 3097 (2001).
  • (18) L. Krusin-Elbaum, L. Civale, V. M. Vinokur, F. Holtzberg, Phys. Rev. Lett. 69, 2280 (1992).
  • (19) L. Civale, L. Krusin-Elbaum, J. R. Thompson, F. Holtzberg, F., Phys. Rev. B 50, 7188 (1994).
  • (20) Y. Abulafia, A. Shaulov, Y. Wolfus, R. Prozorov, L. Burlachkov, Y. Yeshurun, D. Majer, E. Zeldov, H. Wühl, V. B. Geshkenbein, V. M. Vinokur, Phys. Rev. Lett. 77, 1596 (1996).
  • (21) S. Kokkaliaris, P. A. J. de Groot, S. N. Gordeev, A. A. Zhukov, R. Gagnon, L. Taillefer, Phys. Rev. Lett. 82, 5116 (1999).
  • (22) B. Rosenstein, A. Knigavko, Phys. Rev. Lett. 83, 844 (1999).
  • (23) G. Yang, P. Shang, S. D. Sutton, I. P. Jones, J. S. Abell, C. E. Gough, Phys. Rev. B 48, 4054 (1993).
  • (24) S. L. Lee, P. Zimmermann, H. Keller, M. Warden, I. M. Savić, R. Schauwecker, D. Zech, R. Cubitt, E. M. Forgan, P. H. Kes, T. W. Li, A. A. Menovsky, Z. Tarnawski, Phys. Rev. Lett. 71, 3862 (1993).
  • (25) B. Khaykovich, E. Zeldov, D. Majer, T. W. Li, P. H. Kes, M. Konczykowski, Phys. Rev. Lett. 76, 2555 (1996).
  • (26) A. Piriou and Y. Fasano and E. Giannini and Ø. Fischer, Physica C 460-462, 408 (2007).
  • (27) A. Piriou, Y. Fasano, E. Giannini, Ø. Fischer, Phys. Rev. B 77, 184508 (2008).
  • (28) T. Tamegai and Y. Iye and I. Oguro and K. Kishio, Physica C 213, 33 (1993).
  • (29) N. Avraham, B. Khaykovich, Y. Myasoedov, M. Rappaport, H. Shtrikman, D. E. Feldman, T. Tamegai, P. H. Kes, Ming Li, M. Konczykowski, K. van der Beek, E. Zeldov, Nature 411, 451 (2001).
  • (30) A. Piriou, E. Giannini, Y. Fasano, C. Senatore, Ø. Fischer, Phys. Rev. B 81, 144517 (2010).
  • (31) J. D. Moore, K. Morrison, K. A. Yates, A. D. Caplin, Y. Yeshurun, L. F. Cohen, J. M. Perkins, C. M. McGilvery, D. W. McComb, Z. A. Ren, J. Yang, W. Lu, X. L. Dong, and Z. X. Zhao, Supercond. Sci. Technol. 21, 092004 (2008).
  • (32) C. Senatore, R. Flükiger, M. Cantoni, G. Wu, R. H. Liu, and X. H. Chen, Phys. Rev. B 78, 054514 (2008).
  • (33) R. Prozorov, M. E. Tillman, E. D. Mun and P. C. Canfield, New J. Phys. 11, 035004 (2009).
  • (34) Huan Yang, Huiqian Luo, Zhaosheng Wang, and Hai-Hu Wen, Appl. Phys. Lett. 93, 142506 (2008).
  • (35) R. Prozorov, N. Ni, M. A. Tanatar, V. G. Kogan, R. T. Gordon, C. Martin, E. C. Blomberg, P. Prommapan, J. Q. Yan, S. L. Bud’ko, P. C. Canfield, Phys. Rev. B 78, 224506 (2008).
  • (36) R. Prozorov, M. A. Tanatar, P. Prommapan, R. T. Gordon, N. Ni, S. L. Bud’ko, P. C. Canfield, Physica C 469, 667-673 (2009).
  • (37) Bing Shen, Peng Cheng, Zhaosheng Wang, Lei Fang, Cong Ren, Lei Shan, and Hai-Hu Wen, Phys. Rev. B 81, 014503 (2010).
  • (38) S. Salem-Sugui, Jr., L. Ghivelder, A. D. Alvarenga, L. F. Cohen, K. A. Yates, K. Morrison, J. L. Pimentel, Jr., Huiqian Luo, Zhaosheng Wang, and Hai-Hu Wen, Phys. Rev. B 82, 054513 (2010).
  • (39) P. Das, Ajay D. Thakur, Anil K. Yadav, C. V. Tomy, M. R. Lees, G. Balakrishnan, S. Ramakrishnan, A. K. Grover, Phys. Rev. B 84, 214526 (2011).
  • (40) M. R. Eskildsen, L. Ya. Vinnikov, T. D. Blasius, I. S. Veshchunov, T. M. Artemova, J. M. Densmore, C. D. Dewhurst, N. and Ni, A. Kreyssig, S. L. Bud’ko, P. C. Canfield, A. I. Goldman, Phys. Rev. B 79, 100501(R) (2009).
  • (41) D. S. Inosov, T. Shapoval, V. Neu, U. Wolff, J. S. White, S. Haindl, J. T. Park, D. L. Sun, C. T. Lin, E. M. Forgan, M. S. Viazovska, J. H. Kim, M. Laver, K. Nenkov, O. Khvostikova, S. Kühnemann, V. Hinkov, Phys. Rev. B 81, 014513 (2010).
  • (42) M. R. Eskildsen, E. M. Forgan, H. Kawano-Furukawa, Rep. Prog. Phys. 74, 124504 (2011).
  • (43) G. P. Mikitik and E. H. Brandt, Journal of Low Temperature Physics 139, 221 (2005).
  • (44) Y. Yeshurun, A. P. Malozemoff, A. Shaulov, Rev. Mod. Physics 68, 911 (1996).
  • (45) P. W. Anderson, Phys. Rev. Lett. 9, 309 (1962).
  • (46) P. W. Anderson, Y. B. Kim, Rev. Mod. Phys. 36, 39 (1964).
  • (47) D. Zola, M. Polichetti, C. Senatore, and S. Pace, Phys. Rev. B 70, 224504 (2004).
  • (48) M. V. Feigel’man, V. B. Geshkenbein, A. I. Larkin, and V. M. Vinokur, Phys. Rev. Lett. 63, 2303 (1989).
  • (49) G. Blatter, M. V. Feigel’man, V. B. Geshkenbein, A. I. Larkin, and V. M. Vinokur, Rev. Mod. Physics 66, 1125 (1994).
  • (50) V. M. Vinokur, M. V. Feigel’man, V. B. Geshkenbein, Phys. Rev. Lett. 67, 915 (1991).
  • (51) E. Zeldov, N. M. Amer, G. Koren, A. Gupta, M. W. McElfresh, and R. J. Gambino, Appl. Phys. Lett. 56, 680 (1990).
  • (52) M. P. Maley, J. O. Willis, H. Lessure, M. E. McHenry, Phys. Rev. B 42, 2639 (1990).
  • (53) M. E. McHenry, S. Simizu, H. Lessure, M. P. Maley, J. Y. Coulter, I. Tanaka, H. Kojima, Phys. Rev. B 44, 7614 (1991).
  • (54) Lu Zhang, J. Z. Liu, R. N. Shelton, M. D. Lan, Phys. Rev. B 50, 7092 (1994).
  • (55) R. Viennois, E. Giannini, D. van der Marel, R. Cerny, Journal of Solid State Chemistry 183, 769-775 (2010).
  • (56) T. J. Liu, X. Ke, B. Qian, J. Hu, D. Fobes, E. K. Vehstedt, H. Pham, J. H. Yang, M. H. Fang, L. Spinu, P. Schiffer, Y. Liu, Z. Q. Mao, Phys. Rev. B 80, 174509 (2009).
  • (57) H. Okamoto, L. E. Tanner, Bull. Alloy Phase Diagrams 11 (1990).
  • (58) E. C. Stoner, Philosophical Magazine 36, 803 (1945).
  • (59) A. K. Pramanik, L. Harnagea, C. Nacke, A. U. B. Wolter, S. Wurmehl, V. Kataev, B. Büchner, Phys. Rev. B 83, 094502 (2011).
  • (60) A. A. Zhukov, H. Küpfer, G. Perkins, L. F. Cohen, A. D. Caplin, S. A. Klestov, H. Claus, V. I. Voronkova, T. Wolf, H. Wühl, Phys. Rev. B 51, 12704 (1995).
  • (61) C. P. Bean, Phys. Rev. Lett. 8, 250 (1962).
  • (62) C. P. Bean, Rev. Mod. Phys. 36, 31 (1964).
  • (63) T. Taen, Y. Tsuchiya, Y. Nakajima, T. Tamegai, Phys. Rev. B 80, 092502 (2009).
  • (64) Y. Liu, R. K. Kremer, C. T. Lin, Eur. Phys. Lett. 92, 57004 (2010).
  • (65) R. Griessen, Wen Hai-hu, A. J. J. van Dalen, B. Dam, J. Rector, H. G. Schnack, S. Libbrecht, E. Osquiguil, Y. Bruynseraede, Phys. Rev. Lett. 72, 1910 (1994).
  • (66) W. A. Fietz, W. W. Webb, Phys. Rev. 178, 657 (1969).
  • (67) D. Dew-Hughes, Phil. Mag. 30, 293 (1974).
  • (68) M. R. Koblischka, A. J. J. van Dalen, T. Higuchi, S. I. Yoo, M. Murakami, Phys. Rev. B 58, 2863 (1998).
  • (69) E. J. Kramer, J. Appl. Phys. 44, 1360(R) (1973).
  • (70) J. W. Ekin, Experimental techniques for low-temperature measurements, Oxford University Press (New York) 2006 pp 412-467.
  • (71) J. W. Ekin, Cryogenics 20, 611 (1980).
  • (72) R. Flükiger, C. Senatore, M. Cesaretti, F. Buta, D. Uglietti, B. Seeber, Supercond. Sci. Technol. 21, 054015 (2008).
  • (73) M. Putti, I. Pallecchi, E. Bellingeri, M. R. Cimberle, M. Tropeano, C. Ferdeghini, A. Palenzona, C. Tarantini, A. Yamamoto, J. Jiang, J. Jaroszynski, F. Kametani, D. Abraimov, A. Polyanskii, J. D. Weiss, E. E. Hellstrom, A. Gurevich, D. C. Larbalestier, R. Jin, B. C. Sales, A. S. Sefat, M. A. McGuire, D. Mandrus, P. Cheng, Y. Jia, H. H. Wen, S. Lee and C. B. Eom, Supercond. Sci. Technol. 23, 034003 (2010).
  • (74) D. B. Montgomery, W. Sampson, Appl. Phys. Lett. 6, 108 (1965).
  • (75) Y. Shapira, L. J. Neuringer, Phys. Rev. A 140, 163 (1965).
  • (76) M. R. Koblischka, M. Murakami, Supercond. Sci. Technol. 13, 738-744 (2000).
  • (77) C. Tarantini, S. Lee, Y. Zhang, J. Jiang1, C. W. Bark, J. D. Weiss, A. Polyanskii, C. T. Nelson, H. W. Jang, C. M. Folkman, S. H. Baek, X. Q. Pan, A. Gurevich, E. E. Hellstrom, C. B. Eom, and D. C. Larbalestier, Appl. Phys. Lett. 96, 142510 (2010).
  • (78) Hefei Hu, Jian-Min Zuo, Jinsheng Wen, Zhijun Xu, Zhiwei Lin, Qiang Li, Genda Gu, Wan Kyu Park, Laura H Greene, New Journal of Physics 13, 053031 (2011).
  • (79) Xiaobo He, Guorong Li, Jiandi Zhang, A. B. Karki, Rongying Jin, B. C. Sales, A. S. Sefat, M. A. McGuire, D. Mandrus, E. W. Plummer, Phys. Rev. B 83, 220502(R) (2011).
  • (80) R. Khasanov, M. Bendele, A. Amato, P. Babkevich, A. T. Boothroyd, A. Cervellino, K. Conder, S. N. Gvasaliya, H. Keller, H.-H. Klauss, H. Luetkens, V. Pomjakushin, E. Pomjakushina, B. Roessli, Phys. Rev. B 80, 140511(R) (2009).
  • (81) L. Li, Z. R. Yang, Z. T. Zhang, L. Pi, S. Tan, Y. H. Zhang, New Journal of Physics 12, 063019 (2010).
  • (82) N. Katayama, S. Ji, D. Louca, S. Lee, M. Fujita, T. J. Sato, J. Wen, Z. Xu, G. Gu, G. Xu, Z. Lin, M. Enoki, S. Chang, K. Yamada, J. M. Tranquada, J. Phys. Soc. Jpn. 79, 113702 (2010).
  • (83) M. Bendele, P. Babkevich, S. Katrych, S. N. Gvasaliya, E. Pomjakushina, K. Conder, B. Roessli, A. T. Boothroyd, R. Khasanov, H. Keller, Phys. Rev. B 82, 212504 (2010).
  • (84) C. S. Yadav, P. L. Paulose, Journal of Physics: Conference Series 273, 012108 (2011).
  • (85) Huan Yang, Cong Ren, Lei Shan, and Hai-Hu Wen, Phys. Rev. B 78, 092504 (2008).
  • (86) H. G. Schnack, R. Griessen, J. G. Lensink, C. J. van der Beek, P. H. Kes, Physica C 197, 337-361 (1992).
  • (87) M. Jirsa, L. Pust, H. G. Schnack, R. Griessen, Physica C 207, 85-96 (1993).
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