arxiv:cond-mat/ v1 [cond-mat.str-el] 18 Apr 2005

Transcription

1 Determination of the magnetic omain size in the ferromagnetic superconuctor UGe by three imensional neutron epolarization S. Sakarya an N. H. van Dijk Interfaculty Reactor Institute, Delft University of Technology, Mekelweg 15, 69 JB Delft, The Netherlans arxiv:con-mat/ v1 [con-mat.str-el] 18 Apr 005 E. Brück Van er Waals Zeeman Institute, University of Amsteram, Valckenierstraat 65, 1018 XE Amsteram, The Netherlans (Date: July 4, 018) Three imensional neutron epolarization measurements have been carrie out on singlecrystalline UGe between 4 K an 80 K in orer to etermine the average ferromagnetic omain size. It is foun that below T C = 5 K uniaxial ferromagnetic omains are forme with an estimate magnetic omain size of 4 5 µm. PACS numbers: q, Cc, Ch I. INTRODUCTION Recently, the compoun UGe has attracte much attention because superconuctivity was foun to coexist with ferromagnetism. 1, Until this iscovery, only superconucting compouns exhibiting antiferromagnetic orer ha been known like DyMo 6 S 8, GMo 6 S 8, an TbMo 6 S 8. 3,4,5 Coexistenceofantiferromagnetismansuperconuctivity was also foun in the heavy fermion compouns like CeIn 3, CeP Si, an UP Al 3. 6,7,8 In these cases superconuctivity an antiferromagnetism appear simultaneously, because the Cooper pairs are insensitive to the internal fiels arising from the antiferromagnetic orering when the superconucting coherence length ξ is much larger than the perioicity of the static antiferromagnetic orere structure. However, in a ferromagnetic structure we expect that the internal fiels o not cancel out on the length scale of ξ an therefore have their influence on the Cooper pairs. That ferromagnetic orer exclues superconuctivity, is nicely emonstrate in ErRh 4 B 4 9,10 an HoMo 6 S 8 11, where stanar BCS singlet-type superconuctivity is suppresse when ferromagnetic orer sets in. Otherwise, if one woul consier unconventional spin-triplet superconuctivity meiate by ferromagnetic spin fluctuations, the pairing is relatively insensitive to a local magnetic fiel an can therefore coexist with ferromagnetic orer. On the other han, when the ferromagnetic omain size is much smaller than the superconucting coherence length ξ, one effectively has no internal magnetic fiel. The coherence length ξ for UGe is estimate 1,1 to be Å. Interestingly, Nishioka et al. 13,14, consiering jumps in the magnetization at regular intervals of magnetic fiel an at very low temperatures, estimate the ferromagneticomain sizetobe ofthe orerof40å, by attributing the jumps to fiel-tune resonant tunnelling between quantum spin states. Since woul be several times smaller than ξ, it was propose that the ferromagnetism can be cancelle out on the scale of the coherence length of the Cooper pairs. This woul imply that the pairing mechanism for superconuctivity might be of the singlet-type after all. UGe crystallizes in the orthorhombic ZrGa crystal structure (space group Cmmm) 15 with unit cell imensions a = Å, b = Å, an c = Å, containing 4 formula units. Ferromagnetic orer sets in at T C = 5 K. The saturate magnetic moment at ambient pressureis 1.4 µ B /U, irecte along the a axis. Magnetic measurements inicate a very strong magnetocrystalline anisotropy. 16 Superconuctivity is foun only in a limite pressure range between 10 an 16 kbar with a maximum transition temperature of T c 0.7 K. In this pressure range, the magnetic moment is still 1 µ B /U. Within the ferromagnetic phase, a secon transition occurs at T X 5 K at ambient pressure, below which the magnetic moment is enhance. Therefore the temperature region from T C to T X was name the weakly polarize phase, while the lower temperature region T < T X was calle the strongly polarize phase. 17 In this paper we report on three imensional neutron epolarization measurements performe on singlecrystalline UGe at ambient pressure between 4 K an 80 K. Our principal aim was to etermine the ferromagnetic omain size in UGe an compare the value to the size of 40 Å estimate by Nishioka et al. 13 on the basis of the hypothesis of fiel-tune resonant tunnelling between spin quantum states. Since the neutron is a very sensitive probe to local magnetic fiels, neutron epolarization is an excellent technique to measure the average omain size an the omain wall with. II. EXPERIMENT A. Experimental The measurements were performe on the Poly Axis Neutron Depolarization Analyzer (PANDA) at the Interfaculty Reactor Institute (IRI) of the Delft University of Technology. The use neutron wave length was.03 Å

2 which correspons to a velocity of 1949 m/s. The neutron epolarization measurements on UGe were performe on a single-crystalline sample with imensions a b c = mm 3. The b axis was oriente along the transmitte neutron beam (x) with a transmission length L an the easy axis for magnetization a along the vertical axis (z) within the plate of the sample. The crystal has been grown from a polycrystalline ingot using a Czochralski tri-arc technique. No subsequent heat treatment was given to the crystal. The illuminate area was a rectangle with imensions y z = 1 mm centere at the mile of the sample. The measurements in zero fiel were performe uring a temperature sweep from K up to 80 K an own to Kwith a lowsweeprate of10k/hr. The measurements in non-zero fiel (4 an 8 mt) were one uring a similar temperature sweep with a sweep rate of 5 K/hr. The sample was first zero-fiel coole, whereafter the fiel was switche on at the start of the measurements. The subsequent measurements were performe uring heating an cooling in a constant fiel. B. Neutron epolarization The neutron epolarization(nd) technique is base on the loss of polarization of a polarize neutron beam after transmission through a (ferro)magnetic sample. Each neutron unergoes only a series of consecutive rotations on its passage through the (ferro)magnetic omains in the sample. It is important to note that the beam cross section covers a huge number of omains, which results in an averaging over the magnetic structure of the whole illuminate sample volume. This averaging causes a loss of polarization, which is etermine by the mean omain size an the mean irection cosines of the omains. The rotation of the polarization uring transmission probes the average magnetization. The 3 3 epolarization matrix D in a ND experiment expresses the relation between the polarization vector before( P 0 ) an after( P 1 ) transmissionthroughthe sample ( P 1 = D P 0 ). The polarization of the neutrons is create an analyze by magnetic multilayer polarization mirrors. In orer to obtain the complete matrix D, one polarization rotator is place before the sample an another one right after the sample. Each rotator provies the possibility to turn the polarization vector parallel or antiparallel to the coorinate axes x, y, an z. The resultant neutron intensity is finally etecte by a 3 He etector. The polarization rotators enable us to measure any matrix element D ij with the ai of the intensity of the unpolarize beam I S : I S = I ij +I ij. (1) where I ij is the intensity for P 0 along i an P 1 along j. The matrix element D ij is then calculate accoring to D ij = 1 P 0 I S I ij I S () where P 0 is the egree of polarization in the absence of a sample. In our case we have P 0 = We now introuce the correlation matrix α ij : L α ij = x B i (x,y,z) B j (x,y,z), (3) 0 where B( r) = B( r) B is the variation of the magnetic inuction,... enotes the spatial average over the sample volume an where the integral is over the neutron transmission length L through the sample. Assuming α ij 0 for i j we efine the correlation function ξ as ξ = i α ii. (4) With these two quantities it can be shown that if there is no macroscopic magnetization ( B = 0) the epolarization matrix is iagonal an uner the assumption of α ij 0 for i j given by 18,19,0 D ii = e γ v L{ξ αii} i = x,y,z. (5) where γ = s 1 T 1 is the gyromagnetic ratio of the neutron an v its velocity. Intrinsic anisotropy is the epolarization phenomenon that for magnetically isotropic meia the epolarization epens on the orientation of the polarization vector with respect to the propagation irection of the neutron beam. The origin of this intrinsic anisotropy is the emagnetization fiels aroun magnetize volumes in the sample. In the following we will assume that the emagnetization fiels are negligible for neele shape magnetic omains. We now iscuss the case B 0. When the sample shows a net magnetization, the polarization vector will rotate in a plane perpenicular to the magnetization irection. If the sample shape gives rise to stray fiels, the rotation angle φ is relate to the net magnetization M by φ = η γ v Lµ 0 M = η γ v Lµ 0M S m (6) where η is a geometrically factor given in Eq. A14 for a rectangular shape sample an m = M/M S the reuce sample magnetization in terms of the saturation magnetization M S = M S (T). If the mean magnetic inuction B inthesampleisorientealongthez-axis,the epolarization matrix is, for φ (γ/v) α xx α yy L/ (the weak amping limit), given by 18,19,0 } { D xx = D yy = e γ v L { D xy = D yx = e γ v L ξ αxx+αyy ξ αxx+αyy D zz = e γ v L{ξ αzz}, D xz = D zx = D zy = D yz = 0. cosφ, } sinφ, (7)

3 3 D xx UGe et D 1. et D 1. D yy D zz mt FIG. 1: The iagonal elements of the measure epolarization matrix D for increasing an ecreasing temperature for UGe. All other elements of the epolarization matrix are zero within the experimental uncertainty. For D xx an D yy the experimental uncertainty is within the symbol size. With the net magnetization along the z-axis, the rotation angle φ of the beam polarization is obtaine from the measurements by ( ) Dxy D yx φ = arctan (8) D xx +D yy an ξ is calculate with ξ = v ln{etd}/γ L. (9) As mentione earlier, ND also provies information about the mean-square irection cosines of the magnetic inuction vector in the (ferro)magnetic omains. These are irectly given by the quantities γ i = α ii /ξ, where i = x,y,z, an can be estimate from the measurements by γ i = 1 lnd ii /ln{etd}. (10) This equation is only vali for those irections which show no net rotation of the beam polarization. III. RESULTS A. Measurements in zero fiel In Fig. 1 we show the iagonal elements of the epolarization matrix for UGe, measure in zero magnetic fiel. All off-iagonal elements are zero within the experimental uncertainty in the stuie temperature range. The measurements for increasing temperature are qualitatively the same as those for ecreasing temperature, as expecte. The Curie temperature of T C = 5 K is clearly inicate in Fig. 1 by the kink in D xx an D yy. Note that (ra) UGe UGe 0. 4 mt mt 3.0 (ra) FIG. : The eterminant of the measure epolarization matrix etd an the rotation angle φ of the beam polarization after passage through the sample of UGe in 4 mt an 8 mt for increasing an ecreasing temperature. D xx D yy below T C inicates that there is no intrinsic anisotropy an hence that the magnetic omains prouce virtually no stray fiels. Furthermore, D zz 1 inicates that all moments are oriente along the a axis. B. Measurements in small fiel In Fig. we show the eterminant of the epolarization matrix etd an the rotation angle φ after passage through the sample for measurements in fiels of respectively 4 an 8 mt (after zero-fiel cooling). The ata of φ have been correcte by subtracting the mean value above T C, since this rotation is merely ue to the applie fiel. At low temperatures the magnetic fiels (applie after zero fiel cooling) are too small to fully align the magnetic omains. Therefore, the measurements for increasing an ecreasing temperature o not yiel the same results. Whereas for increasing temperature the rotation shows an increase, for ecreasing temperature the ata represent a monotonous magnetization curve, as expecte for a fiel-coole ferromagnet. Close to T C there is no ifference between fiel cooling or fiel warming. Fig. shows that for 4 mt the epolarization is at the same level as for 0 mt. Above T X, however, extra epolarization occurs. This means the system gets more inhomogeneous, i.e. the omains grow an the magnetic correlation length increases (ξ in Eq. 7), leaing to extra epolarization. Close to T C the epolarization isappears, because the magnetic moment ecreases sharply. For ecreasing temperature the eterminant has the same shape as in the case of 0 mt. At 8 mt the eterminant is alreay reuce below T X, inicating larger omains. Again, the Curie temperature of T C = 5 K is clearly inicate by the kink in etd an φ. Also note the abrupt increase in φ aroun T X 0 K. Eviently

4 4 the system passes, with increasing temperature, from a strongly polarize phase to a weakly polarize phase, as reporte earlier. 17 IV. DISCUSSION A. Moel an results The measurements confirm that UGe is a highly anisotropic uniaxial ferromagnet. Further, the magnetic omains are long compare to their (average) with, because D xx D yy inicates relatively weak stray fiels prouce by the magnetic omains. This allows us to assume B( r) = µ 0 M( r) insie the omains. In orer to analyze the ata we consier a moel where the sample is split into N long neeles along the a axis with a fixe with an a magnetic inuction B S = µ 0 M S along the a axis. With N (N ) the number of omains with a magnetic inuction pointing upwars (ownwars), we can efine the reuce macroscopic magnetization of the sample, pointing along the z-irection, as m z = N N N +N = B z B S (11) Each neele will have magnetic inuction or with probability p = (1 + m z )/ an p = (1 m z )/, respectively. The polarize neutron beam traversing the sample will therefore see a binomial istribution of an, which results in a epolarization matrix D with elements B L D xx = D yy = e γ S v D xy = D yx = e γ B S L v ( ) (1 m z ) γbs L cos m z, v ) ( (1 m z ) γbs L sin m z v D zz = 1, (1) an all other elements equal to 0. (Note that, since we have not taken into account the macroscopic stray fiels, the angle φ shoul be correcte by the factor of η (Eq. 6) before calculating m z in Eq. 1.) Within this binomial istribution moel it is easy to show that for the case m z = 0 the average ferromagnetic omain size is equal to. Given a omain wall (i.e. two ajacent neeles with an opposite magnetic inuction), the probability of forming a omain of n neeles is ( 1 n ) an the average is calculate by n=1 n( 1 n ) =. When a fiel is applie, we have to istinguish between a omain (with size ) in which the magnetic inuction is parallel to the fiel an a omain (with size ) with opposite inuction. The probability of forming a omain of size n is p n 1 p = p n 1 (1 p ) which leas to / = 1/(1 p ) = /(1 m z ). Similarly, / = /(1+ m z ). In orer to estimate the omain wall thickness δ we assume that m z changes sinusoially from +1 to -1 over, <m z>, ( m) UGe 4 mt <m z>, ( m) UGe mt FIG. 3: Calculate values of the reuce macroscopic magnetization m z an the average ferromagnetic omain sizes with magnetic inuction parallel ( ) or antiparallel ( ) to the applie magnetic fiel of 4 mt an 8 mt in UGe for increasing an ecreasing temperature. aistance δ in the formofablochwall. The consequence isthat D zz is slightlyless than1in the orerestate. For such a omain wall it is straightforwar to show that the omain wall thickness δ can be estimate by γ z = m z = 1 1 ( ) δ, (13) which can be measure irectly by Eq. 10. For the values of B S neee in Eq. 1, we use the experimental magnetic moment of Ref. 1 which we convert to magnetic inuction, remembering there are 4 formula units per unit cell. For the value of η in Eq. 6 we take η =. From Fig. 1 it is clear that the ata for increasing an ecreasing temperature give slightly ifferent results for the ferromagnetic omain size in zero magnetic fiel. The values foun for = are 5.1() µm when cooling own slowly an 4.4(1) µm when heating up after fast cooling. Both values are inepenent of temperature. These values inicate the omain size perpenicular to the a axis (along the b axis). Along the a axis which we assume the omain size is much larger. The magnetic omain wall thickness δ, ivie by the magnetic omain size, is calculate with Eq. 13 from the experimental ata in Fig. 1, an amounts to δ/ = 47(3), inepenent of temperature. This gives δ = 0. µm. The size of the omain wall thickness is thus foun to be only a minor fraction of the omain size. Analysis of the ata in a small magnetic fiel (Fig. ) with Eq. 1 gives the results shown in Fig. 3 an table I. For 4 mt an increasing temperature (after zero fiel cooling), the reuce magnetization m z remains equal to 0 up to T X 0 K. As a consequence is equal to an of the same orer of the zero fiel values. Above T X however the system gets magnetically soft an m z starts to increase linearly towars 0.7. Domain walls

5 5 TABLE I: Ferromagnetic omain sizes in UGe for increasing temperature after zero fiel cooling (ZFC) an ecreasing temperature in fiel (FC). The sizes of the omains with magnetization parallel to the applie magnetic fiel is enote by an the omains with antiparallel magnetization by. Below T X the omain sizes are temperature inepenent. Above T X the omains grow. The values shown are at a few Kelvin below T C. µ 0H Temp. (µm) (µm) (µm) (µm) (mt) incr./ecr. T < T X T < T X T T C T T C 0 ZFC, incr. 4.4(1) 4.4(1) 4.4(1) 4.4(1) 0 FC, ecr. 5.1() 5.1() 5.1() 5.1() 4 ZFC, incr. 3.9(1) 3.8(1) 150(0) 5(5) 4 FC, ecr. 46.4(8) 9.5() 60(10) 13() 8 ZFC, incr. 17.9() 11.4(1) 85(0) 10() 8 FC, ecr. 45(5) 8.(1) 85(0) 10(5) are expelle above T X, since increases much faster than. (Note the logarithmic vertical scale.) While getsofthe orerof100µm, onlyreaches5µm. When the omains grow in with, at a certain moment it is no longer allowe to assume B( r) = µ 0 M( r), because stray fiels prouce by the omains have to be taken into account. The moel therefore is no longer appropriate close to T C. For fiel cooling in 4 mt, the system has m z = 68(1) for the whole temperature range below T C. The values of the omain size are shown in table I. When after zero fiel cooling a fiel of 8 mt is turne on, the sample oes get a macroscopic magnetization, in contrast to the case of 4 mt. Up to T X 0 K the reucemagnetization m z = 0.1()isinepenent of temperature. Then m z startsto increaseup to 0.718(3) aroun 30 K an is constant afterwars up to T C. When cooling own in 8 mt, m z = 0.708(1) over the whole temperature range below T C. For fiel warming after zero fiel cooling, the calculation of the omain sizes yiels unexpecte temperature epenencies of omain sizes above T X. As can be seen in Fig. 3, accoring to the moel the omain sizes grow abovet X to ecreasein sizeagainat highertemperature. Clearly there is another source of epolarization, not accounteforbythe moel. Sincethe fielisstrongenough to penetrate the sample, aitional epolarization arises from an inhomogeneous magnetic omain structure. B. Discussion If the omain with becomes relatively large compare to its length, stray fiels become important an the simple moel assuming B( r) = µ 0 M( r) is no longer vali. Calculation of the mean-square irection cosine along the z irection, γ z, witheq.10, ineeshowsaecreasefrom unity above T X, inicating that the magnetic inuction B is not along the a axis throughout the sample. The moel can of course be improve if we no longer assume a length/with ratio of infinity (no stray fiels) for the omains. The simple moel together with our measurements, however, o show that the magnetic omain sizes in zero fiel are a few µm an that by applying small fiels the omains grow. Our measurements therefore inicate that the omain sizes in UGe at ambient pressure an own to K are certainly larger than the 40 Å preicte by Nishioka et al. 13,14 In Fig. 1 it is shown that D zz is less than unity below T C. This can be cause by the omain walls, but can also be accounte for by a misalignment. A simple calculation shows that a misalignment of 8 woul fully account for the values of D zz below T C. The state value of δ = 0. µm (or δ/ = 47(3)) shoul therefore be regare as an upper limit. From the above consierations we conclue that the omain structure of UGe behaves like in a conventional ferromagnet. The magnetic omain size largely excees the superconucting correlation length of the Cooper pair. The magnetic omain bounaries can therefore only give seconary effects on the superconucting orer. V. CONCLUSION The ferromagnetic omain sizes of UGe was stuie by means of three imensional neutron epolarization at ambient pressure. We conclue that the existence of fiel-tune resonant tunneling between spin quantum states 13,14 is highly unlikely. The requirement of this moel is a ferromagnetic omain size of 40 Å while our measurements inicate a size a factor of 1000 larger. The observe jumps in the magnetization shoul be attribute to a Barkhausen effect as iscusse by Lhotel et al.. The superconuctivity therefore exists within a single ferromagnetic omain. The omain walls are not expecte to strongly affect the bulk Cooper pair wave function, as suggeste by Nishioka et al. 13,14, since the omain wall is less than a few percent of the average omain size. APPENDIX A: EFFECT OF STRAY FIELDS INDUCED BY A HOMOGENEOUSLY MAGNETIZED SAMPLE In this appenix we will calculate the magnetic inuction B generate by a uniformly magnetize sample with length l, with w, an thickness t (Fig. 4). Moreover, analytical expressions will be given for the line integrals of B along the path of a neutron. The center of the sample is taken as the origin of the reference frame.

6 6 z y x l FIG. 4: Schematic layout of a homogeneously magnetize rectangular sample. Our starting point is the Biot-Savart law: B(x,y,z) = µ 0 M n r 4π S r 3 S + µ 0 M r 4π τ r 3 τ (A1) where µ 0 = 4π 10 7 H/m, M is the magnetization, n the unit vector perpenicular to the surface S of the sample, r the vector pointing from the surface S to the point (x,y,z), an τ the volume of the sample. Since the sample has a homogeneous magnetization, the secon term vanishes an, with m = M/M, we have B = µ 0M 4π S w t M m n r r 3 S. (A) A straightforwar but teious calculation shows that B x = ǫ 1 ǫ ǫ 3 ln(s (y +ǫ w)) B y = B z = where S = ǫ 1,ǫ,ǫ 3=±1 ǫ 1,ǫ,ǫ 3=±1 ǫ 1,ǫ,ǫ 3=±1 ǫ 1 ǫ ǫ 3 ln(s (x+ǫ 1 t)) (A3) ( ) z +ǫ 3 l ǫ 1 ǫ ǫ 3 arctan (x+ǫ 1 t)(y +ǫ w) S (x+ǫ 1 t) +(y+ǫ w) +(z +ǫ 3 l). (A4) The rotation of the polarization of a neutron beam epens on the line integral of the magnetic fiel along the neutron path. From the Larmor equation tp(t) = γp(t) B(t), orequivalently xp(x) = (γ/v) P(x) B(x) where γ = s 1 T 1 the gyromagnetic ratio of the neutron an v its velocity, we get the generalsolution { ( γ x )} P(x,y,z) = exp B(x,y,z)x P(,y,z) v (A5) where we have efine the magnetic fiel tensor B as B(x,y,z) = 0 B z B y B z 0 B x B y B x 0 (x,y,z) (A6) Thus, in orer to calculate the rotation of the neutron beam polarization ue to a homogeneously magnetize sample, the following line integrals are require: X(y,z) = Y(y,z) = Z(y,z) = B x (x,y,z)x = 0 B y (x,y,z)x = tµ 0M ǫ ǫ 3 ln ((y+ǫ w) +(z +ǫ 3 l) ) 4π ǫ,ǫ 3 tµ 0 M π B z (x,y,z)x = ( ) z +ǫ3 l ǫ ǫ 3 arctan y+ǫ ǫ w,ǫ 3 (A7) For completeness we also give the line integrals in the case the neutron beam is along the z-irection: X (x,y) = Y (x,y) = Z (x,y) = B x (x,y,z )z = 0 B y (x,y,z )z = 0 B z (x,y,z )z = µ 0 M l (A8) From Eq. A5 an the above line integrals we get for the final polarization P(,y,z) = D(y,z) P(,y,z) with D(y,z) equal to Σcosa Σ Z Σsina Σ Y Σsina Σ 1 Z Σsina Σ Z cosa Σ+Y YZ ( 1 cosa Σ ) Σ Y Σsina Σ YZ ( 1 cosa Σ ) Y cosa Σ+Z (A9) where Σ(y,z) = Y (y,z)+z (y,z) an a = γ/v. Now Eq. A9 relates the initial polarization to the final polarization for a beam passing through the sample at (y,z). For a neutron beam with finite cross section the matrix D shoul be integrate over the beam cross section. If the integration is symmetric relative to the origin, then we can make use of the fact that B y (ǫ 1 x,ǫ y,ǫ 3 z) = ǫ ǫ 3 B y (x,y,z), (A10) where ǫ 1,ǫ,ǫ 3 = ±1. This means that B y (x,y,z) an thereforey(y,z)areantisymmetric with respectto y an z. Therefore, from Eq. A9, we only have to integrate the matrix 1 Σ Σcosa Σ Z Σsina Σ 0 Z Σsina Σ Z cosa Σ+Y Y cosa Σ+Z (A11) over the cross section of the neutron beam. An infinitely narrow neutron beam passing exactly through the mile of the sample will only have its

7 7 polarization rotate by B z (x,0,0) since B x (x,0,0) an B y (x,0,0) vanish. As long as Y (y,z) is small compare to Z (y,z), which is vali if (y,z) is sufficiently far from the eges, Eq. A11 is a pure rotation matrix. It is now possible to calculate the magnetization of the sample from the measure rotation angle (Eq. 8). If no stray fiels were present, the rotation angle woul be given by tµ 0 Mγ/v. However, the stray fiels reuce the rotationangleto γ/vz(y,z)with Z(y,z)giveninEq. A7. We can therefore efine the geometrical factor η as φ = ηµ 0 Mtγ/v (A1) Since η(0,0) is a sale point (η(0,z) has a local maximum an η(y,0) a local minimum) an average over the cross section of the neutron beam, centere aroun the mile of the sample, will yiel a result very close to the value of η(0,0), which is given by η = π arctan ( l w ). (A14) where η(y,z) is given by η(y,z) = 1 ( ) z +ǫ3 l ǫ ǫ 3 arctan. (A13) π y +ǫ ǫ w,ǫ 3 Electronic aress: S.Sakarya@iri.tuelft.nl 1 S. S. Saxena, P. Agarwal, K. Ahilan, F. M. Grosche, R. K. W. Haselwimmer, M. J. Steiner, E. Pugh, I. R. Walker, S. R. Julian, P. Monthoux, et al., Nature (Lonon) 406, 587 (000). D. Aoki, A. Huxley, E. Ressouche, D. Braithwaite, J. Flouquet, J.-P. Brison, E. Lhotel, an C. Paulsen, Nature (Lonon) 413, 613 (001). 3 D. E. Moncton, G. Shirane, W. Thomlinson, M. Ishikawa, an O. Fischer, Phys. Rev. Lett. 41, 1133 (1978). 4 C. F. Majkrzak, G. Shirane, W. Thomlinson, M. Ishikawa, O. Fischer, an D. E. Moncton, Soli State Commun. 31, 773 (1979). 5 W. Thomlinson, G. Shirane, D. E. Moncton, M. Ishikawa, an O. Fischer, Phys. Rev. B 3, 4455 (1981). 6 I. R. Walker, F. M. Grosche, D. M. Freye, an G. G. Lonzarich, Physica C 8, 303 (1997). 7 F. M. Grosche, S. J. S. Lister, F. V. Carter, S. S. Saxena, R. K. W. Haselwimmer, N. D. Mathur, S. R. Julian, an G. G. Lonzarich, Physica B 39, 6 (1997). 8 N. K. Sato, N. Aso, K. Miyake, R. Shiina, P. Thalmeier, G. Varelogiannis, C. Geibel, F. Steglich, P. Fule, an T. Komatsubara, Nature (Lonon) 410, 340 (001). 9 W. A. Fertig, D. C. Johnston, L. E. DeLong, R. W. McCallum, M. B. Maple, an B. T. Matthias, Phys. Rev. Lett. 38, 987 (1977). 10 D. E. Moncton, D. B. McWhan, P. H.Schmit, G. Shirane, W. Thomlinson, M. B. Maple, H. B. MacKay, L. D. Woolf, Z. Fisk, an D. C. Johnston, Phys. Rev. Lett. 45, 060 (1980). 11 M. Ishikawa an O. Fischer, Soli State Commun. 3, 37 (1977). 1 N. Tateiwa, T. C. Kobayashi, K. Hanazono, K. Amaya, Y. Haga, R. Settai, an Y. Onuki, J. Phys. Conens. Matter 13, L17 (001). 13 T. Nishioka, G. Motoyama, S. Nakamura, H. Kaoya, an N. K. Sato, Phys. Rev. Lett. 88, 3703 (00). 14 T. Nishioka, N. K. Sato, an G. Motoyama, Phys. Rev. Lett. 91, 0970 (003). 15 P. Boulet, A. Daoui, M. Potel, H. Noel, G. M. Gross, G. Anre, anf. Bouree, J. Alloys Comp. 47, 104(1997). 16 Y. Onuki, I. Ukon, S. W. Yun, I. Umehara, K. Satoh, T. Fukuhara, H. Sato, S. Takayanagi, M. Shikama, an A. Ochiai, J. Phys. Soc. Jpn. 61, 93 (199). 17 A. Huxley, I. Sheikin, an D. Braithwaite, Physica B 84-88, 177 (000). 18 M. T. Rekvelt, Z. Phys. 59, 391 (1973). 19 R. Rosman an M. T. Rekvelt, Z. Phys. B 79, 61 (1990). 0 R. Rosman an M. T. Rekvelt, Phys. Rev. B 43, 8437 (1991). 1 C. Pfleierer an A. D. Huxley, Phys. Rev. Lett. 89, (00). E. Lhotel, C. Paulsen, an A. D. Huxley, Phys. Rev. Lett. 91, (003).