Abstract

A limiting factor in the ability to interpret isotope effect measurements in cuprates is the absence of sufficiently accurate data for the whole phase diagram; there is precise data for , but not for the antiferromagnetic transition temperature . Extreme sensitivity of to small changes in the amount of oxygen in the sample is the major problem. This problem is solved here by using the novel compound ()() for which there is a region where is independent of oxygen doping. Meticulous measurements of for samples with and find the absence of an oxygen isotope effect on with unprecedented accuracy. A possible interpretation of our finding and existing data is that isotope substitution affects the normal state charge carrier density.


Isotope substitution is a powerful experimental tool used to investigate complex systems. Ideally, the isotope substitution affects only one parameter, for example, a phonon frequency which is directly related to the nuclear mass. However, the strong coupling of many parameters in cuprates highly limit the ability to interpret isotope effect (IE) experiments. Even though the oxygen isotope substitution is known to affect the superconductivity transition temperatures , it is unclear whether it primarily impacts phonons, polarons, magnons, doping, or other physical properties [1].

The isotope effect is usually described using the isotope coefficient via the relation where is a phase transition temperature, is the isotope mass and , and for the superconducting (SC), antiferromagnetic (AFM), and spin glass critical temperatures, respectively. In many conventional superconductors was found to be very close to 0.5 [2, 3]. The explanation of this in terms of Cooper pairs glued by phonons was one of the triumphs of the BCS theory for metallic superconductors [4].

In cuprates the isotope effect is much more complicated and is not single valued and varies across the phase diagram. The consensus today for YBCO like compounds is that in the SC phase, close to optimal doping, for oxygen substitution [5]. On the SC dome increases as the doping decreases [6, 7]. In the glassy state less data is available, but it seems that the isotope effect reverses sign and becomes negative. In extremely underdoped samples, where long range AFM order prevails at low temperatures, data is scarce, controversial, and has relatively large error bars [5, 8]. The most recent measurements with show that in the parent compound [5]. There are also several theories dealing with the variation of along the phase diagram [9โ€“14], but since and are within an error bar of each other one cannot contrast these theories with experiments. In particular, it is impossible to tell whether the same glue that holds the spins together holds the cooper pair together, or not. Increasing the accuracy of the IE measurements of the Nรฉel temperature will shed light on the role of isotope substitution.

As mentioned above, a major limitation on the measurements of is the strong dependence of on doping. For example, in , decreases with increasing doping at a rate โ€‰K per . This strong temperature dependence is, of course, common to many other cuprates. As a consequence, it is very difficult to prepare two samples with exactly the same even with the same isotope. The smallest fluctuation in either or may lead to a huge fluctuation in regardless of the isotope effect. This is not the case for ()() where is constant for oxygen density , as can be seen in Figure 1 [15].

For our experiments, four sintered pellets were prepared using standard techniques [16]. Two of the pellets were enriched with isotope and two with isotope in the same procedure: the samples were placed simultaneously in two closed tubes, each with different isotope gas, and then they were heated to allow the isotope to diffuse into the sample. In order to achieve a higher percentage of gas, the enrichment was repeated several times.

The isotope content in the samples was determined based on measurements of gas composition being in equilibrium with the sample during the exchange. Balzers Prisma mass spectrometer was used to analyze in situ isotopic composition of the atmosphere. After the exchange process was performed, the weight increase of the sample was also determined as the light isotope was exchanged with the heavy . The isotope enrichment in the samples measured by both methods looked to be higher than 80%. Finally a Thermal Analysis (TA) experiment was performed for the investigated samples after all experiments described in this work were fulfilled in order to verify the isotope enrichment. The samples were heated up to 1200ยฐC in the NETSCH STA 449C Jupiter analyzer in a stream of helium. During the TA experiments, ion current signals for the , OO, and molecules were measured using a mass spectrometer (ThermoStar Pfeiffer Vacuum). The results are shown in Figure 2. The isotope content deduced from these measurements (comparing peak areas for the signals of particular oxygen molecules) was larger than 70%.

The oxygen IE was measured using zero-field muon spin rotation/relaxation (SR). We particularly used the ISIS facility, which allows low muon relaxation rate and rotation frequency measurements. This is ideal for measurements near magnetic phase transitions where the muon signal varies on a long-time scale. SR data at temperatures close to the phase transition are shown in Figure 3. As the temperature is lowered from 383.4โ€‰K, the relaxation rate increases. At 378.5โ€‰K oscillations appear in the data indicating the presence of long range magnetic order. The frequency of oscillations and the relaxation rate increase as the temperature is further lowered. The formula was fitted to the muon polarization, where , , and are the polarization, relaxation rates, and frequency of muons spin in the fraction of the samples which is magnetic, and is the weighting factor between the muons experiencing transverse field and longitudinal field. This factor is close to 2/3 and temperature independent. and are the polarization and relaxation of the spin of muons that stopped in the nonmagnetic volume of the sample. The solid lines in Figure 3(a) represent the fits.

The AFM order parameter is determined by the frequency , where is the average magnetic field at muon site and is the muon gyromagnetic ratio. This frequency can be easily extracted from the SR data well below the transition but is very difficult to determined near the transition. A second approach is to treat the magnetic volume fraction as the order parameter [5, 17]. can be followed more closely to , although this parameter also has large error bars when is not well defined. We used an alternative approach similar to [18]; we define an order parameter, which does not require a fit, via the relation where is the average polarization at temperature , and is the average polarization above the transition. The denominator normalizes to 1 at zero temperature. All three order parameters for one sample are shown in Figure 3(b). The transition temperatures determined using the different order parameters are in good agreement. The advantages of are clear: it is a model-free and has very small uncertainties.

In Figure 4, we present the for two samples with and one with around 380โ€‰K. A wide temperature range from 50 to 410โ€‰K is shown in the inset for an experiment done on separate occasion, in which two samples of and one with were examined. We determine by fitting a straight line to the data in the main panel of Figure 4 in the temperature range 378 to 382โ€‰K, for each sample, and taking the crossing with the temperature axis. We find that โ€‰K and โ€‰K. For 100% isotope substitution the isotope exponent is determined by When taking into account the isotopic fraction in the samples we obtain This result indicates that (see (2)) beyond the error bars and is consistent with no isotope effect on .

One possible interpretation of these results is that magnetic excitations are not relevant for superconductivity since the isotopes affect without affecting . This approach was presented, for example, by Zhao et al. [9]. They found that samples enriched with have longer penetration depth than samples enriched with with the same amount of oxygen per unit cell. is related to the SC carrier density and effective mass by so a priory both and can be affected by isotope substitution [19]. They ruled out the possibility that the number of carriers concentration varies by demonstrating that the thermal expansion coefficient of samples with different isotopes are the same. The authors therefore concluded that the IE changes the mass of the cooper pairs, which could be explained by polaronic supercarriers.

An alternative interpretation is that the isotopes affect the efficiency of doping, as suggested in [13]. To demonstrate this interpretation we present in Figure 5(a) the critical temperatures in for the two different isotopes taken from [5]. seems to be a bit higher than , but is a bit lower than . However, if we define an efficiency parameter which relates the number of holes to the number of oxygen ions in the unit cell via , where stands for the isotope type, we can generate a unified phase diagram. This is demonstrated in Figures 5(b) and 5(c) for and , respectively. In these graphs three different values of are used to generate while keeping . When using , both curves of and versus for the two different isotopes collapse to the same curve. Similar scaling of the doping axis was applied to the () () system [15] and was explained by NQR [20].

Next we discuss the IE on the stiffness in the above scenario. We assume that the effective mass of the SC charge carrier , the critical where superconductivity starts , and where is optimal , are not affected by the isotope substitution. The stiffness can be measured by the muon transverse relaxation rate , and it is expected that [19] where is a constant. Dividing the differential of from (7) by the relaxation at optimal doping yields

The expected change in the stiffness due to isotope substitution can be calculated from (8) using , (which are extracted from Figure 5(c)), and โ€‰ฮผs-1 [21]. For , and we get . This value is consistent with the experimental value of โ€‰ฮผs-1 reported in [21]. In other words a 2% difference in the doping efficiency between the two isotopes can explain both the variations in the phase diagram and the variation in the stiffness.

Our experiment shows that oxygen isotope substitution does not affect the Nรฉel temperature and therefore does not play a role in magnetic excitations. However, the isotope effect of does not necessarily imply that phonons play a role in cuprate superconductivity. We show that an isotope-dependent doping efficiency can explain the variation in and in the magnetic penetration depth between samples rich in or .

Acknowledgments

The authors would like to thank the ISIS pulsed muon facility at Rutherford Appleton Laboratory, UK for excellent muon beam conditions. This work was funded by the Israeli Science Foundation, the joint German-Israeli DIP project, and the Posnansky research fund in high temperature superconductivity.