per Cu

4
Anomalous slowing down of Cu-spin fluctuations observed by muon spin relaxation in the Zn-substituted Bi 2 Sr 2 Ca 1 2x Y x Cu 1 2y Zn y 2 O 8 1d system around the hole concentration of 1 8 per Cu I. Watanabe Muon Science Laboratory, The Institute of Physical and Chemical Research (RIKEN), 2-1 Hirosawa, Wako, Saitama 351-0198, Japan M. Akoshima and Y. Koike Department of Applied Physics, Graduate School of Engineering, Tohoku University, Aramaki Aoba-yama 08, Aoba-ku, Sendai 980-8579, Japan K. Nagamine Muon Science Laboratory, The Institute of Physical and Chemical Research (RIKEN), 2-1 Hirosawa, Wako, Saitama 351-0198, Japan and Meson Science Laboratory, Institute of Materials Structure Science, High Energy Accelerator Research Organization (KEK-MSL), 1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan ~Received 1 March 1999! Muon spin relaxation ( m SR) measurements were applied to a partially Zn-substituted Bi-2212 system, Bi 2 Sr 2 Ca 12x Y x (Cu 12y Zn y ) 2 O 81d , to study a possibility of the existence of the so-called ‘‘ 1 8 effect’’ which has been established in La systems. The muon spin depolarization rate has been found to increase with decreasing temperature below 10 K in Zn-substituted samples around the hole concentration p 5 1 8 per Cu where the high-T c superconductivity is anomalously suppressed. It has been revealed that the enhancement of the muon spin depolarization rate is attributed to the slowing down of the Cu-spin fluctuations and that both Zn substitution and p 5 1 8 per Cu are essential for the slowing down of the Cu-spin fluctuations. Although no three-dimensional long-range coherent ordering of the Cu spins appeared down to 0.30 K, the present m SR results support the existence of the 1 8 effect in a Bi-2212 system also. @S0163-1829~99!50838-2# Anomalous suppression of the high-T c superconductivity ~SC! within a narrow region of the hole concentration p ; 1 8 per Cu was first observed in La 2 2x Ba x CuO 4 ~LBCO! around x 50.125. 1,2 This phenomenon is well known as the so-called 1 8 effect. A structural phase transition from orthorhombic ~space group: Bmab! to tetragonal ~space group: P 4 2 / ncm ) was discovered in this hole-concentration range, indicating a relationship between the anomalous suppression of the SC and the structural change. 3 A similar phenomenon was also found in La 2 2x Sr x CuO 4 ~LSCO!, although no structural tran- sition was observed except for instability of a phonon mode. 4 Muon spin relaxation ( m SR) measurements using a posi- tive muon beam were performed for both LBCO and LSCO and revealed the existence of an antiferromagnetically ~AF! ordered state where the SC was anomalously suppressed. 5–7 The magnetic transition temperature showed a maximum at x 50.125 for LBCO and at x 50.115 for LSCO ~Ref. 8!. Subsequently, from neutron-diffraction measurements on Nd-substituted LSCO ~LNSCO! undergoing a structural transition, 9 Tranquada et al. suggested the existence of a stripe structure of a spin-density wave ~SDW! state and a charge-density wave ~CDW! state. 10 However, the detailed mechanism of the 1 8 effect and the ordered state of spins and charges are still in controversy. At the moment, the 1 8 effect has been established only in La systems. Thus, it is important to investigate it in other high-T c cuprates which have different crystal structures from those of the La systems in order to clarify the relationship between the crystal structure and the 1 8 effect. Recently, Akoshima et al. showed one possibility. That is, they found anomalous suppression of the high-T c SC in the partially Zn-substituted Bi-2212 system, Bi 2 Sr 2 Ca 1 2x Y x (Cu 1 2y Zn y ) 2 O 8 1d , around p 5 1 8 per Cu ( x 50.3125), though a structural transition similar to those ob- served in LBCO ~Ref. 3! and LNSCO ~Ref. 9! did not take place. 11 They suggested that the substituted Zn atoms acted as charge and/or spin pinning centers, and that the SC was suppressed by the similar mechanism to the 1 8 effect. The anomalous suppression of the SC has been confirmed by transport measurements on the Zn-substituted Bi-2212 sys- tem, while the existence of a magnetically ordered state, which is expected to be derived from the 1 8 effect, has not yet been confirmed. In this paper, we report on the results of zero field ~ZF! and longitudinal field ~LF! m SR measurements in a partially Zn-substituted Bi-2212 system. We have observed enhance- ment of the muon spin depolarization rate at low tempera- tures below 10 K in Zn-substituted samples around p 5 1 8 per Cu, indicating the slowing down of the Cu-spin fluctuations. It has been found that both Zn-substitution and p 5 1 8 per Cu are essential for the slowing down of the Cu-spin fluctua- tions. Our results support the existence of the 1 8 effect in the Bi-2212 system. Samples of Bi 2 Sr 2 Ca 1 2x Y x (Cu 1 2y Zn y ) 2 O 8 1d were pre- pared by the conventional solid-state reaction method. 11 We prepared three sorts of polycrystalline samples with different hole concentrations: x 50.20 ( p . 1 8 per Cu!, x 50.3125 ( p ; 1 8 per Cu!, and x 50.40 ( p , 1 8 per Cu!. Zn concentrations RAPID COMMUNICATIONS PHYSICAL REVIEW B 1 OCTOBER 1999-II VOLUME 60, NUMBER 14 PRB 60 0163-1829/99/60~14!/9955~4!/$15.00 R9955 ©1999 The American Physical Society

Transcript of per Cu

Page 1: per Cu

Japan

Japan

RAPID COMMUNICATIONS

PHYSICAL REVIEW B 1 OCTOBER 1999-IIVOLUME 60, NUMBER 14

Anomalous slowing down of Cu-spin fluctuations observed by muon spin relaxationin the Zn-substituted Bi2Sr2Ca12xYx„Cu12yZny…2O81d system

around the hole concentration of 18 per Cu

I. WatanabeMuon Science Laboratory, The Institute of Physical and Chemical Research (RIKEN), 2-1 Hirosawa, Wako, Saitama 351-0198,

M. Akoshima and Y. KoikeDepartment of Applied Physics, Graduate School of Engineering, Tohoku University, Aramaki Aoba-yama 08, Aoba-ku,

Sendai 980-8579, Japan

K. NagamineMuon Science Laboratory, The Institute of Physical and Chemical Research (RIKEN), 2-1 Hirosawa, Wako, Saitama 351-0198,

and Meson Science Laboratory, Institute of Materials Structure Science, High Energy Accelerator Research Organization(KEK-MSL), 1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan

~Received 1 March 1999!

Muon spin relaxation (mSR) measurements were applied to a partially Zn-substituted Bi-2212 system,Bi2Sr2Ca12xYx(Cu12yZny)2O81d , to study a possibility of the existence of the so-called ‘‘1

8 effect’’ whichhas been established in La systems. The muon spin depolarization rate has been found to increase withdecreasing temperature below 10 K in Zn-substituted samples around the hole concentrationp5

18 per Cu

where the high-Tc superconductivity is anomalously suppressed. It has been revealed that the enhancement ofthe muon spin depolarization rate is attributed to the slowing down of the Cu-spin fluctuations and that both Znsubstitution andp5

18 per Cu are essential for the slowing down of the Cu-spin fluctuations. Although no

three-dimensional long-range coherent ordering of the Cu spins appeared down to 0.30 K, the presentmSRresults support the existence of the1

8 effect in a Bi-2212 system also.@S0163-1829~99!50838-2#

edic

gS

o-e

si-O

d

ola

nd

inem

hi

t-

,

b-

tedwas

bysys-te,

lynce-ra-

s.

a-

ent

Anomalous suppression of the high-Tc superconductivity~SC! within a narrow region of the hole concentrationp; 1

8

per Cu was first observed in La22xBaxCuO4 ~LBCO! aroundx50.125.1,2 This phenomenon is well known as the so-call18 effect. A structural phase transition from orthorhomb~space group:Bmab! to tetragonal~space group:P42 /ncm)was discovered in this hole-concentration range, indicatinrelationship between the anomalous suppression of theand the structural change.3 A similar phenomenon was alsfound in La22xSrxCuO4 ~LSCO!, although no structural transition was observed except for instability of a phonon mod4

Muon spin relaxation (mSR) measurements using a potive muon beam were performed for both LBCO and LSCand revealed the existence of an antiferromagnetically~AF!ordered state where the SC was anomalously suppresse5–7

The magnetic transition temperature showed a maximumx50.125 for LBCO and atx50.115 for LSCO~Ref. 8!.Subsequently, from neutron-diffraction measurementsNd-substituted LSCO~LNSCO! undergoing a structuratransition,9 Tranquadaet al. suggested the existence ofstripe structure of a spin-density wave~SDW! state and acharge-density wave~CDW! state.10 However, the detailedmechanism of the18 effect and the ordered state of spins acharges are still in controversy.

At the moment, the18 effect has been established only

La systems. Thus, it is important to investigate it in othhigh-Tc cuprates which have different crystal structures frothose of the La systems in order to clarify the relationsbetween the crystal structure and the1

8 effect.

PRB 600163-1829/99/60~14!/9955~4!/$15.00

aC

.

.at

n

r

p

Recently, Akoshimaet al. showed one possibility. Thais, they found anomalous suppression of the highTc

SC in the partially Zn-substituted Bi-2212 systemBi2Sr2Ca12xYx(Cu12yZny)2O81d , aroundp5 1

8 per Cu (x50.3125), though a structural transition similar to those oserved in LBCO~Ref. 3! and LNSCO~Ref. 9! did not takeplace.11 They suggested that the substituted Zn atoms acas charge and/or spin pinning centers, and that the SCsuppressed by the similar mechanism to the1

8 effect. Theanomalous suppression of the SC has been confirmedtransport measurements on the Zn-substituted Bi-2212tem, while the existence of a magnetically ordered stawhich is expected to be derived from the1

8 effect, has not yetbeen confirmed.

In this paper, we report on the results of zero field~ZF!and longitudinal field~LF! mSR measurements in a partialZn-substituted Bi-2212 system. We have observed enhament of the muon spin depolarization rate at low tempetures below 10 K in Zn-substituted samples aroundp5 1

8 perCu, indicating the slowing down of the Cu-spin fluctuationIt has been found that both Zn-substitution andp5 1

8 per Cuare essential for the slowing down of the Cu-spin fluctutions. Our results support the existence of the1

8 effect in theBi-2212 system.

Samples of Bi2Sr2Ca12xYx(Cu12yZny)2O81d were pre-pared by the conventional solid-state reaction method.11 Weprepared three sorts of polycrystalline samples with differhole concentrations:x50.20 (p. 1

8 per Cu!, x50.3125 (p; 1

8 per Cu!, andx50.40 (p, 18 per Cu!. Zn concentrations

R9955 ©1999 The American Physical Society

Page 2: per Cu

rie

Ae

on

aThym

searte

d

etinitehe. Aizt.ceti

y

thha

oxyethsa

se

theftheon

edys-ntme

Omons

ds

-

inso-ysis

epo-.30lesea-

thnce.per-

hat

d

p

po-

ow

RAPID COMMUNICATIONS

R9956 PRB 60I. WATANABE, M. AKOSHIMA, Y. KOIKE, AND K. NAGAMINE

werey50 andy50.025 for eachx value. All of the sampleswere checked to be of the single phase, and the charactetion of the samples was performed by resistivity and susctibility measurements.

mSR measurements were carried out at the RIKEN-RMuon Facility at the Rutherford-Appleton Laboratory in thUnited Kingdom. A spin-polarized pulsed surface mubeam with a momentum of 27 MeV/c was used. Forwardand backward counters were located at the upstreamdownstream sides of the sample in a beam direction.muon spin polarization was parallel to the beam. The asmetry parameter was defined asA(t)5@F(t)2B(t)#/@F(t)1B(t)#, whereF(t) and B(t) were the total muon eventcounted by the forward and backward counters at timt,respectively. A LF was applied along the muon spin polization. The time evolution of the asymmetry parame(mSR time spectrum! was measured.

Figure 1 shows ZF-mSR time spectra of a Zn-substitutesample withx50.3125 andy50.025 in which the high-TcSC is anomalously suppressed. The time spectrum abovK shows a Gaussian-type depolarization behavior, indicathat no significant magnetic effect exists at the muon sexcept for randomly distributed nuclear dipolar fields. Ttime spectrum deviates from a Gaussian type below 10 Klow temperatures, it shows an exponential-type depolartion behavior, as reported from a preliminary experimen12

This change in the time spectrum suggests the appearanan additional effect due to local fields. Since the magneproperties of this system are ruled by the Cu spins, itinferred that this additional effect reflects their static or dnamical properties.

In the analysis of the time spectra, the difference ofmuon sites should be taken into account. It is known tmuons injected in a high-Tc oxide stop near oxygen ions.13

From the comparison with previousmSR studies in Y~Ref.14! and Tl systems,15 there are expected to be roughly twkinds of muon sites in a Bi-2212 system; sites near the ogen ions constituting the CuO5 pyramid and sites near thoxygen ions in the Bi-O plane. Muons stopping nearCuO5 pyramid feel strong internal fields from the Cu spinand muons stopping near the Bi-O plane which are far awfrom the CuO5 pyramid feel weaker internal fields, becau

FIG. 1. Zero-field mSR time spectra of Zn-substituteBi2Sr2Ca12xYx(Cu12yZny)2O81d with x50.3125 andy50.025 at19.2, 1.74, and 0.30 K. The solid lines are the best fit of Eq.~1!.The high-Tc superconductivity of this sample is anomalously supressed~Ref. 11!.

za-p-

L

nde-

-r

10gs

ta-

ofics-

et

-

e,y

the internal fields at the muon sites are dominated bydipolar fields of the Cu spins.16 Because of this difference ointernal fields, the depolarization of the muon spins nearCuO5 pyramid is expected to be faster than that of the muspins near the Bi-O plane. In fact, two components~fast andslow! of the depolarization behavior were already observin a magnetically ordered state of an insulating Bi-2212 stem by Nishidaet al.16 Therefore, we take a two-componefunction to describe the observed exponential-type tispectrum, as given by

A1e2l1t1A2e2l2t3GZ~D,HLF ,t !. ~1!

The first term describes muons which are close to the Cu5pyramid and depolarize mainly by the additional fields frothe Cu spins. The second term describes the other muwhich are far away from the CuO5 pyramid and depolarizemainly by the randomly distributed nuclear dipolar fielwith a weak effect of the Cu spins. GZ(D,HLF ,t) indicatesthe static Kubo-Toyabe function, whereD is the distributionwidth of the nuclear dipolar fields at the muon sites.17 HLF isthe LF. A1 andA2 are the initial asymmetries of each component att50. In the case of ZF-mSR, we put zero in HLF .

The solid lines in Fig. 1 show the best fit of Eq.~1! to thetime spectra.l1 is close tol2 at high temperatures, whichmeans that the effect of the internal fields from the Cu spbecomes small. Therefore, it is difficult to apply the twcomponent fitting. In this case, a single-component analusing only the second term is adopted for convenience.

Figure 2 shows the temperature dependence of the dlarization rate of all of the measured samples down to 0K. The depolarization rate of the non-Zn-substituted sampshows a very weak temperature dependence within the msured temperature range. The Zn-substituted sample wix50.20 also shows a similar weak temperature dependeThese results indicate that the static and dynamical proties of the Cu spins do not change in the case ofp. 1

8 per Cueven in a Zn-substituted sample. However, it is found t

-

FIG. 2. Temperature dependence of zero-field muon spin delarization rate of Bi2Sr2Ca12xYx(Cu12yZny)2O81d with x50.20,0.3125, 0.40 andy50, 0.025. The anomalous enhancement bel4 K of Zn-substituted samples withx50.3125 andx50.40 is due toa temperature dependence of the fast depolarization rate,l1. Inthese cases, the slow depolarization rates,l2, show similar smallvalues with those obtained at higher temperatures.

Page 3: per Cu

hur

heeasiavo

raashthn

avhepiontebeutc

e-tp

iontllaa

df

n

b-

toigys-

theedhat

ioningis

p-r

tice of

t.hisis

fsisis

tateandbu-

isldsis

escethisera-ces-the

e

eof

RAPID COMMUNICATIONS

PRB 60 R9957ANOMALOUS SLOWING DOWN OF Cu-SPIN . . .

the depolarization rate of the Zn-substituted samples witx50.3125 andx50.40 increases with decreasing temperatbelow about 10 K.

If the muon spins depolarize by dynamical effects of tCu spins, the enhancement of the depolarization rate mthat the Cu-spin fluctuations are suppressed with decreatemperature, exhibiting the so-called slowing-down behior. On the other hand, if they depolarize by static effectsthe Cu spins, the enhancement of the depolarizationmeans that some Cu spins freeze below 10 K and makeditional static internal fields at the muon sites which cauthis exponential type depolarization of the muon spins. Tstatic depolarization process can be realized in the casethe distribution of the static internal fields is of an exponetial type. Actually, such an exponential depolarization behior by the static internal fields was already observed in otmaterials.17,18In any case, the enhancement of the muon sdepolarization rate suggests that the Cu-spin fluctuatslow down with decreasing temperature in Zn-substitusamples. This fact implies that the magnetic correlationtween Cu spins is expected to be enhanced in Zn-substitsamples resulting in the slowing down of the Cu-spin flutuations.

Figure 3 shows thex dependence, namely, the holconcentration dependence of the depolarization rate ofZn-substituted samples at 0.30, 0.50, and 1.8 K. The delarization rate exhibits a maximum atx50.3125 wherep; 1

8 per Cu and the high-Tc SC is anomalouslysuppressed.11,19 Since no enhancement in the depolarizatrate is observed in the non-Zn-substituted samples, ifound that both Zn substitution andp; 1

8 per Cu are essentiafor the slowing down of the Cu-spin fluctuations. This retionship between the slowing down of the Cu-spin fluctutions and the SC, revealed in the presentmSR studies, isanalogous to that of the18 effect which has been establishein La systems.6,8 Therefore, in conjunction with the results othe transport measurements,11 we can conclude that the18effect is likely to exist in the Bi-2212 system also. This coclusion suggests that the18 effect is characteristic of theCuO2 plane of high-Tc cuprates.

A three-dimensional~3D! long-range ordered state, as o

FIG. 3. x dependence, namely, the hole-concentration depdence of the zero-field muon spin depolarization rate,l1, of Zn-substituted Bi2Sr2Ca12xYx(Cu12yZny)2O81d with y50.025 at0.30, 0.50, and 1.8 K. The depolarization rate is anomalouslyhanced atx50.3125 (p; 1

8 per Cu!. Dotted lines show thex de-pendence ofTc determined by Akoshimaet al. from transport mea-surements withy50 and 0.025~Ref. 11!.

e

nsng-fted-eisat

--r

ins

d-

ed-

heo-

is

--

-

served in La systems,6,8,10 has not been observed down0.30 K in a Zn-substituted Bi-2212 system. This is a bdifference in the state of the Cu spins between the La stems and the Bi-2212 system aroundp5 1

8 per Cu. This re-sult shows that the coherent long-range ordered state ofCu-spins is difficult to be stabilized in the Zn-substitutBi-2212 system. In the case of LNSCO, it is suggested tthe coherent long-range ordered state~the stripe structure ofSDW and CDW states! is derived from the pinning of thedynamical stripe correlations of spins and holes.10 Assumingthat the1

8 effect is due to the stripe structure, no observatof the long-range ordered state indicates that the pinnforce due to Zn atoms in a Zn-substituted Bi-2212 systemmuch weaker than that in LNSCO. A 3D ordering may apear at lower temperatures.mSR measurements at lowetemperatures using a3He-4He dilution refrigerator are nowplanned to search for the long-range ordered state.

LF was applied to the Zn-substituted sample withx50.3125 andy50.025 at 0.30 K to investigate the magnestate of the Cu spins. Figure 4 shows the LF dependencthe time spectrum at 0.30 K. The LF-mSR time spectra areanalyzed with Eq.~1!. Solid lines in the figure are the best fiThe depolarization is found to be suppressed by LF. Tmeans that the additional effect from the Cu spinsquenched by LF.

It is found from this analysis that the distribution width othe static internal fields,D, increases linearly up to a few tenGauss with increasing LF. When a polycrystalline sampleused, the distribution of the applied fields in the sampleusually inhomogeneous and complicate in the bulk SC sbecause of intricate giometrical conditions among grainslarge anisotropy of the electronic state. Since the contrition of the nuclear dipolar fields is at most 1 G in this sys-tem, the wider field distribution seen at the muon sitesexpected to be due to the inhomogeneity of the applied fiein the sample. This result may suggest that the bulk SCstill left in the sample.

The decoupling behavior in LF shown in Fig. 4 indicata possibility that the additional internal fields which induan exponential-type depolarization behavior are static. Incase, the Cu spins are expected to be frozen at low temptures and strongly disordered, because no coherent presion of the muon spins is observed. On the other hand,

n-

n-

FIG. 4. Longitudinal-field dependence ofmSR time spectra ofthe Zn-substituted Bi2Sr2Ca12xYx(Cu12yZny)2O81d with x50.3125 andy50.025 at 0.30 K. The solid lines are the best fitEq. ~1!.

Page 4: per Cu

gpila

ypemviCu

oysheguee

enoptith

iponx-

d

entedratewn0 Kneipnd-aport

ur

lp-heby af

o bych-ion

RAPID COMMUNICATIONS

R9958 PRB 60I. WATANABE, M. AKOSHIMA, Y. KOIKE, AND K. NAGAMINE

effect of additional fields can be also understood assumindynamical reason, namely, the slowing down of the Cu-sfluctuations. This is because a similar change in the depoization behavior from a Gaussian type to an exponential twas observed in La systems just above the 3D ordering tperature, and this exponential-type depolarization behawas concluded to be due to dynamically fluctuatingspins.6

Although more detailed measurements are needed totain results on the depolarization process in LF, if the dnamical depolarization process is assumed, the fitting reof the LF-mSR shows an interesting feature. That is, tdynamical depolarization rate,l1, decreases with increasinLF with roughly a square-root field dependence. Such spression of the depolarization behavior by LF has also bobserved in one-dimensional materials.20 In those cases, thedepolarization of the muon spins is dominated by a magncally excited state moving along one direction. Althoughclear evidence has been obtained yet, there may be twosibilities for a one-dimensional excitation. One is a magneexcitation of the in-plane short-range ordered state alongdirection perpendicular to the CuO2 plane through the weakinterplane AF interaction. The other is based on the strstructure. That is, the stripe structure is periodic alongdirection in the CuO2 plane. Therefore, one-dimensional e

anr-e-

or

b--ult

p-n

ti-

os-ce

ee

citation in the CuO2 plane may be possible. More detaileLF-mSR measurements are necessary to be conclusive.

In summary, ZF- and LF-mSR measurements have becarried out to study magnetic properties of a Zn-substituBi-2212 system. The enhancement of the depolarizationof the muon spins, which corresponds to the slowing doof the Cu-spin fluctuations, has been observed below 1aroundp5 1

8 per Cu in the Zn-substituted samples. Both Zsubstitution andp5 1

8 per Cu are essential to induce thslowing down of the Cu-spin fluctuations. The relationshbetween the slowing down of the Cu-spin fluctuations athe SC, revealed by the presentmSR measurements, is analogous to that of the1

8 effect which has been established in Lsystems. Thus, taking into account the results of the transmeasurements,11 it is concluded that the18 effect is likely toexist in the Bi-2212 system as well as in La systems. OmSR results suggest that the1

8 effect is characteristic of theCuO2 plane of high-Tc cuprates.

The authors would like to acknowledge S. Ohira for heful technical support. A part of this study, especially tsample preparation and characterization, was supportedGrant-in-Aid for Scientific Research from the Ministry oEducation, Science, Sports and Culture, Japan and alsCREST ~Core Research for Evolutional Science and Tenology! of Japan Science and Technology Corporat~JST!.

1A.R. Moodenbaughet al., Phys. Rev. B38, 4596~1988!.2K. Kumagaiet al., J. Magn. Magn. Mater.76&77, 601 ~1988!.3J.D. Axeet al., Phys. Rev. Lett.62, 2751~1989!.4T. Fukaseet al., Physica B165-166, 1289~1990!.5G.M. Luke et al., Physica C185-189, 1175~1991!.6I. Watanabeet al., J. Phys. Soc. Jpn.61, 3058~1992!.7K. Kumagaiet al., Physica C235-240, 1715~1994!.8I. Watanabeet al., Hyperfine Interact.86, 603 ~1994!.9M.K. Crawford et al., Phys. Rev. B44, 7749~1991!.

10J.M. Tranquadaet al., Nature~London! 375, 561 ~1995!; Phys.Rev. B54, 7489~1996!.

11M. Akoshimaet al., Phys. Rev. B57, 7491~1998!.12I. Watanabeet al., Physica B259-261, 557 ~1999!.13S.B. Sulaimanet al., Hyperfine Interact.79, 901 ~1993!.14N. Nishida and H. Miyatake, Hyperfine Interact.63, 183

~1990!.15I. Watanabeet al., Hyperfine Interact.105, 181 ~1997!.16N. Nishidaet al., Physica C156, 625 ~1988!.17Y.J. Uemuraet al., Phys. Rev. B31, 546 ~1985!.18A. Kratzeret al., Hyperfine Interact.64, 449 ~1990!.19S.D. Obertelliet al., Phys. Rev. B46, 14 928~1992!.20K. Nagamineet al., Phys. Rev. Lett.53, 1763~1984!.