NUCLEAR ORIENTATION STUDIES AT LOW TEMPERATURES by PATRICK WILLIAM DALY B.Sc. Bishop's University, 1968 A THESIS SUBMITTED IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY in the Department . - of Physics We accept this thesis as conforming to the required standard THE UNIVERSITY OF BRITISH COLUMBIA March, 1973 ii In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of British Columbia, I agree that the Library shall make i t freely available for reference and study . I further agree thatxpermission for extensive copying of this thesis for scholarly purposes may be granted by the Head of my Department or by his representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. Department of Physics The University of British Columbia Vancouver 8, Canada Date TTUX^A 3o } / ?73 ABSTRACT The interaction between the nuclear magnetic moment and the hyperfine f i e l d in ferromagnetic and antiferromagnetic materials has been used to orient radioactive nuclei at low temperatures. The measurement of the resulting angular distribution of the gamma radiation has yielded a variety of information in both nuclear and solid state physics. 59 The Fe nucleus in an iron sample has been studied in an attempt to verify a report by Tschanz and Sapp (1970) who measure an unexpectedly large effect in the gamma ray distribution from that nucleus oriented in double nitrate salts. Since in their experiment assumptions have to be made regarding the fraction of iron atoms"in lattice sites, the well-defined environment of the iron lattice was chosen for the present experiment in order to permit an unambiguous interpretation of the results. In contradiction to predictions based on the work of Tschanz and Sapp, a null result was obtained at a temperature of 15 mK, 59 indicating an upper limit for the magnetic moment of I I - Fe of | u | < 0.9 u ^ . The anisotropy of gamma radiation emitted in the decay of 103 Ru oriented in ferromagnetic iron was measured to determine the magnetic moment of the isotope and to resolve mixing ratios in the . decay scheme. The small effect observed at 11 mK permitted a lower limit to be placed on the magnetic moment: |p| > 0.15 u . The sign of the E2/M1 mixing ratio for the 497 keV gamma transition was determined for the f i r s t time and found to be negative. iv The techniques of nuclear orientation have also been used to observe the spin-flop transition in antiferromagnetic MnCl^^H^O. By 54 observing the gamma ray anisotropy in the decay of radioactive Mn we find that at 60 mK the transition occurs for a c r i t i c a l field of H = 6.50 ± 0.05 kG applied along the crystallographic c-axis. £ In the spin-flop phase, the electronic spins point at an angle of 71° to this axis. These data allow values of the molecular exchange and anisotropy fields to be determined yielding respectively H = 11.7 ± 0.6 kG and = 2.0 ± 0.1 kG in agreement with the values reported by Miedema et_ a l . (1965). This experiment demonstrates that nuclear orientation is a viable technique for studying spin-flop transitions and moreover allows the measurement of the directions of electronic spin magnetization in the spin-flop phase. 207 The nuclear orientation of determine the hyperfine field. was observed. Bi in cobalt was studied to An unexpectedly small gamma ray anisotropy However, saturation of the effect occurred at a relatively high temperature. These results indicate that only a fraction (^12%) of the Bi nuclei felt a large hyperfine field (^600 kG). this was the fraction of atoms in lattice sites. Presumably Comparison of the anisotropics of the 1.064 MeV and the 1.761 MeV gamma rays shows that the former effect is enhanced, possibly because of the long h a l f - l i f e of the 1.633 MeV intermediate state. Tschanz, J. F. and Sapp, R. C. 1970. Phys. Rev. C, 2_, 2168. 5 Miedema, A. R., Wielinga, R. F. ,and Huiskamp, W. J. 1965. Physica, 31_, 835. V TABLE OF CONTENTS ABSTRACT i i i TABLE OF CONTENTS v LIST OF TABLES vii LIST OF FIGURES viii ACKNOWLEDGEMENTS CHAPTER 1 x THEORY OF NUCLEAR ORIENTATION 1.1 Introduction 1 1.2 Simple Angular Distribution Scheme 2 1.3 Details of the Angular Distribution 8 1.4 Solid Angle Correction 1.5 Methods of Orienting Nuclei . 25 1.6 Low Temperature Production 30 1.7 Other Experimental Methods 33 CHAPTER 2 22 EXPERIMENTAL APPARATUS AND PROCEDURES 2.1 The Cryostat 36 2.2 Paramagnetic Salts and Thermal Contact 39 2.3 Low Temperature Procedure 41 2.4 Gamma Ray Detection 42 2.5 Paper Tape Analysis 49 2.6 Thermometry: ^Cobalt Experiments 54 2.7 The Term "Anisotropy" 63 59 CHAPTER 3 3.1 ORIENTATION OF Fe NUCLEI W IRON Introduction 65 59 3.2 Decay Scheme of Fe 65 59 3.3 Expected Anisotropy of Fe 68 vi 3.4 The Experiment 69 3.5 Analysis 71 3.6 Results and Discussion 77 CHAPTER 4 ORIENTATION OF RUTHENIUM NUCLEI IN IRON 103 4.1 Introduction 4.2 Decay Scheme of 4.3 Expected Anisotropy of the 497 keV Gamma Ray 83 4.4 Experiments 86 4.5 Analysis 88 4.6 Results 95 CHAPTER 5 80 103 Ru 81 SPIN-FLOP TRANSITION OF MnCl -4^0-OBSERVED BY 2 NUCLEAR ORIENTATION 5.1 Introduction 100 5.2 Molecular Field Model 100 5.3 Decay Scheme and Anisotropy 104 5.4 The Experiment 108 5.5 Analysis 111 CHAPTER 6 ORIENTATION OF 2 0 7 B i NUCLEI IN COBALT 6.1 Introduction 117 6.2 Decay Scheme 118 6.3 Expected Anisotropies 118 6.4 Sample 122 6h5 The Experiment 124 6.6 Analysis 128 6.7 Discussion 132 REFERENCES 134 vii LIST OF TABLES TABLE 1 Sample output of the spectrum analysis . 59 , 60 for Fe and Co. c TABLE 2 program 74 r 60, Co intensities during Cobalt temperatures from 75 Fe runs. TABLE 3 Average 59 for Fe for each temperature range. 76 TABLE 4 U and F coefficients for possible transition modes 103 in the 497 keV and 610 keV sequences in Ru. 84 TABLE 5 Coefficients of B and B from °Co background run 93 TABLE 6 Effect in the 497 keV gamma ray as measured by 96 6 2 4 the a x i ? l detectors. TABLE 7 Measured and fitted angles of the spin orientation 113 above the spin-flop transition. TABLE 8 Fitted values of attenuation and hyperfine field 207 from the gamma rays of Bi. 130 viii LIST OF FIGURES FIGURE .•1 Typical decay scheme 3 FIGURE 2 Solid angle correction 23 FIGURE 3 Level splittings and populations during 32 adiabatic demagnetization FIGURE ^4 Outer cryostat 37 FIGURE 5 Inner cryostat 38 FIGURE 6 Flow diagram of detector output 45 FIGURE 7 Background calculation 51 FIGURE 8 FIGURE 9 Co decay scheme 60 55 W(0) and W(ir/2) for Co versus inverse 58 60 temperature FIGURE 10 Sample °Co spectra 60 FIGURE 11 ^Co anisotropy versus polarizing solenoid 6 61 current FIGURE 12 ^Co temperature versus time elapsed from ^. . t _ '62 demagnetization FIGURE 13 59 F e decay scheme 66 59 FIGURE 14 FIGURE 15 versus inverse temperature for Sample Fe/ °Co spectra 59 FIGURE 16 FIGURE 17 Fe ' 6 Ru decay scheme 70 72 82 F versus mixing ratio for 497 and 610 keV 2 85 gamma rays 103 FIGURE 18 Sample spectra from Ru runs FIGURE 19 Gamma ray scattered into detector 89 90 ix FIGURE 20 Background ratio of 497 keV window versus '60r 92 . . Co intensity FIGURE J21 Mn decay scheme 105 54 FIGURE 22 Gamma ray intensity versus cos0 FIGURE 23 Sample Mn/ Co spectra FIGURE 24 Cos Versus field along c-axis 112 FIGURE 25 Cos 03 versus field along c-axis above 115 54 for Mn 60 107 110 spin-flop FIGURE 27 207 Bi decay scheme 207 Saturation of Bi/Co sample FIGURE 28 Sample FIGURE 29 207 Bi gamma ray effects versus inverse temperature FIGURE 226 207 Bi/ °Co spectra 6 119 125 126 131 ACKNOWLEDGEMENTS The author wishes to express his sincere gratitude to The National Research Council of Canada and the Killam Foundation of Vancouver, British Columbia for financial assistance in the forms of an NRC Post-Graduate Scholarship and a Killam Pre-Doctoral Fellowship, Dr. P.W. Martin and Dr. B.G. Turrell for their supervision and guidance and especially for their encouragement during the final weeks, R.L.A. Gorling and R. Kisser for their assistance during the experimental runs, J. Lees, glassblower, for his help in the preparation of the Bi samples, R. Weissbach, low temperature technician, for maintaining a steady supply of liquid helium, ELndLDlyDaly whose innumerable services made the whole affair run smoothly. 1 CHAPTER ONE THEORY OF NUCLEAR ORIENTATION 1.1 Introduction AOfeafaSrfehofmipadi'o'actikeod'e^ ei i s that the radiation i s emitted with equal intensity in a l l directions. However, the spherical symmetry of the radiation field i s not a property of the decay process at a l l but a result of the random arrangement of the nuclei. For many nuclei are definitely not spherical, possessing various vector quantities, such as magnetic moments and angular momentum. In the study of oriented nuclei, a preferred direction for these vectors is introduced and the angular dependence of the radiation is observed, revealing secrets of the decay scheme or of the orienting interactions. The mathematical derivation of the angular distribution as a function of the orienting and radiation parameters is straight-forward, but does not lend i t s e l f to a feeling for what i s happening. Conversely, the physical picture does not lend i t s e l f to a complete and useful formulation. However, we shall f i r s t present the physical picture to offer insight; afterwards the general result will be stated, and then derived. Much use will be made in this chapter of the quantum theory of angular momentum, references to which can be found in Rose (1967). Other general references in this field are Blin-Stoyle and Grace (1957), DeGroot, Tolhoek, and Huiskamp (1965). 2 1.2 Simple Angular Distribution When a n u c l e u s o f of s p i n I^. s p i n 1^ and z component momentum L , t h e ,z component c o n s e r v e a n g u l a r momentum. priate parity, are and t h u s radiation, either z component where are its of These determine or magnetic W(e) M = +1: W(6) = l / 4 ( l - 3 c o s 2 6 M=i:0: W(6) = 3 / 2 ( c o s 2 9 of W(6) = 1, and c a l c u l a t e t h e these t r a n s i t i o n s . tions of the to decay to initial beta-decays total to spin 1 to of example, for are f i v e states the quadrupole possible the - 1.1a + 4cos46) 1.1b cos49) z axis. the 1.1c When t h e s e f i v e fields angular d i s t r i b u t i o n for d i s t r i b u t i o n from the and on t h e p r o b a b i l i t y a spin 2 state allow this) appro- distributions: for each probabilities popula- each i n i t i a l state state. 1 shows a s i t u a t i o n i n w h i c h a p a r e n t different state of These p r o b a b i l i t i e s depend on t h e r e l a t i v e (the z components nucleus of spin 3 b e t a p a r t i c l e and n e u t r i n o o m u s t which promptly Shown a r e t h e p o s s i b l e t r a n s i t i o n s states w i t h the to an i s o t r o p i c d i s t r i b u t i o n . overall a particular final Figure i n order together state angular = l/4(l-cos46) I n p r i n c i p l e , one c a n d e t e r m i n e transition For another total - (M2), t h e r e following angular M = +2: of completely the distribution. e l e c t r i c (E2) in strength, a gamma r a y quantum numbers 6 i s the d i r e c t i o n r e l a t i v e to equal decays to t h e gamma r a y m u s t b e angular values with the and z component by emdssion o f s u f f i c i e n t to gamma ray,., Scheme gamma-decays from the parent differ only Mq to have a spin 0 = 1 state. state. Since i n o r i e n t a t i o n and not in the any Figure 1 Typical decay scheme 4 i n t r i n s i c fashion, i t i s possible to find the r e l a t i v e t r a n s i t i o n p r o b a b i l i t i e s through geometric considerations only. The t r a n s i t i o n p r o b a b i l i t i e s are proportional to the square of the Clebsch-Gordan c o e f f i c i e n t i n the angular momentum coupling of the beta p a r t i c l e and f i n a l state nucleus to 2 the i n i t i a l state: that i s to |CCI^ LQ 1^; - m, m)| . For the t r a n s i - tions i n the diagram, that i s , for the transitions from MQ = 1 to the p r o b a b i l i t i e s are 1/15, 8/15 and 2/5 respectively. = 2,1,0, Each subsequent gamma decay has a p r o b a b i l i t y of u n i t y , since only one f i n a l state e x i s t s . How- ever, i n general, one must allow for several transitions from each i n i t i a l state i n the same way as was done above for the beta decay. If a l l the parent nuclei were i n the = 1 s t a t e , the gamma rays would have the angular d i s t r i b u t i o n 4 2 4 2 4 w(e) = ^5 j ( i - c o s e ) + | j j ( i - 3 c o s e + 4cos e) + |- j ( c o s e - cos e) "I 0 A ^ = £Q(9 + 12cos 0 - 5bos?9).. .. ~ i - Sc.s 1.2 By calculating and averaging W(Q) over each of the parent states, one would be able to determine the angular d i s t r i b u t i o n for any set of parent populations. Consider the case i n which a parent state of spin I , with r e l a t i v e populations a ( M O ) , decays to a state of spin 1^ by emitting radiation of angular momentum L q . I f this process i s followed by prompt emission of radiation with angular momentum to a state of spin I^, then this l a t t e r radiation has an angular d i s t r i b u t i o n 5 W ( 8 ) = 2 m^n ^ V ^ W V o~ > M m M ) ' 2 | C C I 2 , L 1 , I 1 M -m-n, n) | F ^ e ) ; o 1.3 6 where F ^ C ) i s the angular distribution of radiation with angular momentum L and z component n. It should be emphasized that the above formula requires a prompt decay from the intermediate state; that i s , the lifetime of the intermediate state must be short enough that no reorientations can occur before the _9 decay. This requirement usually demands lifetimes less than 10 sec, the typical Larmor precession time; such short lifetimes are quite normal in gamma decays. However, exceptions do exist. The angular distribution F j v for an individual gamma ray can be calculated from the electromagnetic multipole fields. Following the phase convention used by Biedenharn and Rose (1953), we shall use the following expressions for the vector potential of an electromagnetic wave in free space with angular momentum quantum number L, z component quantum number M, and parity i r : > A L M ( ma S netic )= 1 A (electric) - -1 i LM L T \ 11 = LLM L + 1 L J 2L+1 1 1 / 2 h T VI'L.L+I.M 1 + i C _ ) L+l L + 1 L 1 A 1/2 M ^J- l h T 2L+1 V^L.L-l,! IT = (-) where h^ is a spherical Hankel function of the f i r s t kind (an outgoing wave at large distances) and T T L A M = ZC(X,1,L; Vi M i s a vector spherical harmonic defined as M-y,u) Y | 1.5 r The Y,„ are spherical harmonic functions, and the E. are spherical basis AM u r vectors, related to the cartesian unit vectors thus: t - h " " 72 5 ^o =z x+iy x r 1.6 x - - y 72" In the gauge chosen for equations 1.4, the scalar potentials are zero. Hence the radiation intensity at any place is proportional to A*A*. In fact, this density is independent of the parity, so that one may use the magnetic multipole to calculate. FLM(6) -J 2 |C(L,1,L; M-y,y)| |Y L)M _ (e)| y 2 . 1.7 The proportionality constant is adjusted so that W(8) = 1 for a uniform distribution of i n i t i a l states. This outline of the principles behind calculating W(6) as a function of the i n i t i a l state populations does not produce a useful formula. The many sums involved in equations 1.3 and 1.7 can be simplified in terms of Racah coefficients and Legendre polynomials to give a single sum: W(9) = Z B U F P v v v v v (cose) . 1.8 The sum is carried over only even values of v and the maximum value of the index is 2 I Q , 2 1^, 2 L^, whichever is smallest. the summand are as follows: The factors i n 7 The term describes the i n i t i a l state populations and is a function of the i n i t i a l state alone. Bv = (2v+l) The factor U 1/2 Z It is given by C(I 0 ,v,I 0 ; M0,0)a(M0) . describes the reorientation created by the unobserved intermediate transition. It i s a function only of the transition and of the spins of the states which i t connects. Uv = (-) 1.9 Its expression is I1+In-L +v l/ £ ° [C2I 0 +lK2I 1 +l)]- ' W ( I o , I 1 > I ( ) , I 1 ; L 0 v ) where W is a Racah coefficient. 1.10 If several transitions precede the observed one, each i s represented in equation 1.8 by another U coefficient. If the transition has mixed multipolarity, the combined coefficient i s 2 U U (L) + 6 U = — 1 + 6 (L') 1.11 where 6 is the mixing ratio of the multipblarities; that i s , S is the amplitude for L' divided by the amplitude for L. The observed gamma ray transition is described by the factor F , which is a function of the multipolarity and spins of the states involved. Its formula i s I I F, I .(-:). 1" 2 +1 (2L1+1)(2I1+1) 1/2 C(L ,L ,v;l-l) W(I ,L , I ,L ; I ,v) . 1.12 8 For m i x e d m u l t i p o l e s o f a n g u l a r momentum L a n d L ' , F F ( L ) + 26F ( L , L > ) + &2 F (L') = v 1.13 where F £L,L») = f-jW 1 [(2L+1)(2L,+1)(2I + 1 ) ] 1 / 2 WCI^L.I^L'jI^v). x The angular f i n a l term P^fcosO) C(L,L',v;l-l)::r i s a Legendre polynomial, containing the dependence. The U a n d F c o e f f i c i e n t s c a n b e e v a l u a t e d Gordan 1.14 from t a b l e s o f C l e b s c h - a n d R a c a h c o e f f i c i e n t s o r t h e y may b e f o u n d f r o m t a b l e s s u c h a s t h o s e o f F e r e n t z and Rosenzwelg(1955) (F c o e f f i c i e n t s o n l y ) or Yamazaki (1967). The B c o e f f i c i e n t s a r e t h e o n e s most and S a t c h l e r populations but d i f f e r 1.3 (1955). exist following Gray Other systems f o r s p e c i f y i n g t h e s t a t e o f t h e i n i t i a l (see DeGroot, Tolhoek, and Huiskamp (1965) f o r s y n o p s i s ) only i n normalization. Details of the Angular In t h i s s e c t i o n w i l l equation 1 . 8 . The e f f e c t done i n e q u a t i o n 1 . 3 . the r a d i a t i o n , Extension commonly u s e d , other to other Distribution be d e r i v e d t h e a n g u l a r d i s t r i b u t i o n f o r m u l a o f o f multipole mixings w i l l No a s s u m p t i o n s w i l l n o t be i g n o r e d , a s was b e made c o n c e r n i n g t h e n a t u r e t h a n i t s a n g u l a r momentum p r o p e r t i e s , forms o f r a d i a t i o n i s t h e r e f o r e of u n t i l the end. e a s i l y done. 9 We will use the notation implied in figure 1; that i s , a state of spin Ij;is to decay to a state of spin emitting radiation of angular momentum L, L+l, L+2,.... The i n i t i a l state can be described by a (21^+1) x (21^+1) density matrix p . In many situations, the system possesses rotational symmetry about some axis, which we shall c a l l the z axis. In this case, the eigenstates of the Hamiltonian are also eigen- states of the z component of angular momentum, and the density matrix i s diagonal when written with these eigenstates as basis. Therefore <I1M |p|l1M1-> = aCM^d 1.15 1 where a(M^) i s the probability (or relative population) of the state After this state decays, the system contains a nucleus of spin and radiation of angular momentum L,L+l,L+2... . We shall take as basis vectors of this final state the direct product of the radiation and nuclear states. That i s , the final state basis vectors are of the form I1 I„M„>1 IAT..> where1 |AT.,> i s the wave function of the radiation. 2 2 LM LM However, the i n i t i a l state |I^M^> cannot decay to any arbitrary combination of final state vectors; i t may decay only to those combinations which have total angular momentum 1^ and z component M^, i f angular momentum i s to be conserved. Thus |ljM^> decays to I M M I M A W + * \ t 2' ' l' 2' r 2)l 2 2 l L,M1-M9 C L I M > > 1 A 6 10 where a. i s the amplitude f o r the emission of radiation of angular momentum L. This amplitude depends on the nuclear decay, and as f a r as t h i s d i s - cussion i s concerned, i s completely a r b i t r a r y , except that the sum of t h e . squares of the amplitudes must be unity (to maintain the normalization) and that with the phases chosen i n equation 1.4 the amplitudes must be real.''" The f i n a l state wave functions of equation 1.16 w i l l have the same density matrix as the o r i g i n a l s t a t e . To f i n d the p r o b a b i l i t y that the radiation w i l l be at the point with spherical co-ordinates r , 6 , <(>, one m u l t i p l i e s the density matrix by the projection operator f o r r a d i a t i o n at r , 0, <f> and takes the trace of the product. W(r,e,ij>) = Tr[p|r,e,4><r,6,4>|] = E , M I <I 1 M 1 |p|l 1 M 1 'xI 1 M 1 |x,9,«(»<r,e,*|l Mx> v J 5 a(M,)a?.,a A* ( r , 0,<j>) «AT (r,0,<Ji) l L' LT L'jM^-I^ L,M^-M2crj , j M 1 2 x C(I 2 ,L',I 1 ;M 2 ,M 1 -M 2 )C(I 2 ,L,I 1 ;M 2 ,M 1 -M 2 ) 1.17 where use has been made of equation 1.16, to substitute f o r |l.|M >, and of the orthonormality of the f i n a l nuclear s t a t e s . The product <r,6 ,<j> |ALM> ->- has been written as A^(r,0,(j)). The angular co-ordinates are contained e n t i r e l y within the product ^L'M'^LM" T ^ i s function can be expanded i n terms of spherical harmonics, "*"See Biedenharn and Rose (1953) regarding time reversal invariance. 11 the coefficients of which being partially determined by the Wigner-Eckart theorem. Let A A L'M' LM = where Z f (L',L,M;vp) Y yv f(L',L,M;vy) = /A*,M (6,*) 1.18 A Y . LM ^ ^ ' 1 A 9 The dependence of the right-hand side of equation 1.19 on the projection quantum number M can be determined by the Wigner-Eckart theorem, yielding f(L',L,M;vy) = 6 y Q C(L,v,L';M0)(L•||v||L) 1.20 where (L'||v||L) i s called a reduced matrix element and i s independent of any projection quantum numbers. It will be shown that for the electromagnetic wave functions used in this work, the reduced matrix element i s real. However, this does not mean that i t i s symmetric i n L' and L, for i t i s only f that i s symmetric. L (-) '(2L + 1) 1/2 Instead, L (L]|v||L') = (-) (2L'+l)(L'{|v||L). 1.21 The Wigner-Eckart theorem restricts y to zero, which means that 2 in the expansion in equation 1.18, only Yvo(6<)>) = ^[ ^*]P^(cos6) contributes. Tucking the 4TT into the normalization of the reduced matrix element, one can write the angular distribution function of equation 1.17 as W(9) = E a(M )a*,a 1 CU^L',^; M ^ - M ^ L LLi' M 1 2 v C (L, v.,L! ;M-M,0) 1 1/2 x CCI^L,!^ M ,M -M ) (2v+l) ' M 2 1 2 2 (L«||v||L) P^cose). 1.22 The sum over M can be eliminated by using an identity involving 2 Racah coefficients and a sum over three ClebschGGordanccoefficients. The identity i s [(2e+l) ( 2 f + l ) ] W(a,b,c,d;e,f) C(a,f ,c, 1//2 ;a,B+S) = E C(a,b,e;ct,B) C(e,d,c;a+B,6) C (b,d,f;B,6) 1.23 B where B+6 remains constant throughout the sum. By making appropriate sub- stitutions and transformations among the indices of the Clebsc-hrGofdan coefficients, one can show E M = C(I ,L,I ;M ,M -M ) C(I L» . I ^ M ^ - M ^ C (L, v,L • ;M M ,0) = 2 1 2 1 2 2> r 2 2 I (-) 1 + L _ I 2 1/2 [(2L' l)(2I l)] Wa^L.I^L-jI^v) CCI^v.I^M^O) 1.24 i / Z + 1 + Equation 1.24 allows us to remove the sum over in equation 22 and to replace i t by a Racah and a ClebschSGordanccoefficient. The remaining sum over now contains factors which are functions only of the i n i t i a l state parameters. Let B = (2v+l) v 1/2 E a p O C (I, , v, I -M. ,0) . „. 1 1 1 1 l M 1.25 13 With equations 1.24 and 1.25, equation 1.22 may be written as W(9) = E B F P V V V V (cosO) 1.26 where the term B^ contains only the i n i t i a l state parameters and the term F^ contains the radiation and nuclear spin terms. No projection quantum numbers are contained in F . We have v I F v = E a*,a (-) LL ' L Li 1 " 2 +L I L 1/2x / [ (2L ' +1) (21^1) ] Wf^ , L, ^ ,L •; 1^ v) 1 x CL'||v||L). 1.27 Normally, only two multipoles contribute, and one of these is the dominant one. The lesser multipole is measured by giving the mixing ratio 6 = a ,/a. , where the L multipole is the dominant one. T After recalling L i J-i that the amplitudes are real, and that, by equation 1.21, the summand i s symmetric in L and L', one can write equation 1.27 as F = a F (L,L) + 2a a , F (L,L') + af, F (L',L') L v LL v L'v 2 v T T T 1 F (L,L) + 26F (L,L ') + 6 F (L',L ' ) — ^ 1 + <T 2 - 1.28 It remains toe evaluate the reduced matrix element to determine completely the F coefficients. However, the reduced matrix element contains many factors which are of no importance to this discussion. Such factors are radial terms, and 14 electromagnetic normalizations. We desire to normalize the electromagnetic wave functions in such a way that W(0) = 1 when a(M^) = When this uniform distribution occurs, a l l the B (21, + l ) " ^ . 1 except for B which goes to unity. q F q to be unity as well. v vanish Since P is also unity, we require q An explicit expression for the Racah coefficient with a zero argument i s W(a,b,c,d;e,o) = (-) - " e C 6^ [ (201) (2d+l)] ~ d 1 / 2 . 1.29 Inspecting equation 1.27 i n the light of this identity, one sees that the reduced matrix element-with v=0 andL=L' must, be unity. Before calculating the matrix element between electromagnetic states, i t is necessary f i r s t to calculate i t between vector spherical harmonics, as defined in equation 1.5. This in turn will reduce to the matrix element of one spherical harmonic between two other spherical harmonics. The value of this matrix element is given by a standard theorem for these functions: / A'M. v p XM Y Y Y d f l = t 4^2A' ( 2 + ^ M',M y " 6 + 1.30 The matrix element of a spherical harmonic between two vector spherical harmonics can then be written: Y ^ L'A-•M' VOI LAM T T > = Z y C(A',l,L';M-y,y) C(A,1,L';M-y,y) C (A,v,A';M-y,0) 15 where equations 1.5 and 1.30 have been used together with the orthonorraali of the spherical basis vectors. The factor 4ir has been dropped. In equation 1.31, the sura over y contains three Clebsch-Gordan coefficients. The Racah identity of equation 1.23 can be used to replace the sum by one ClebschTGordanand one Racah coefficient; then by combining the remaining .Clebsch-Gordan coefficient with the Racah coefficient, a new sum, independent of M, can be produced. The transitions involved are: E C(A',l,L';M-y,y) C(A,l,L;M-y,y) C(A,v,A ';M-y,0) = ( _ ) - - ' A V L +1 y [(2L+1)(2A'+1)] 1/2 W(L,A,L',A';l,v) C(L,v,L';M,0) 1.32 followed by a Racah permutation WX'L,A;L: ,,A';1,V) , = W ( A , L , A ' , L ' ;1 ,v) 1.33 and transition back to a sum W(A,L,A',L';l,v) C (A,v,A';0,0) = ( - ) x E C (L ,1 ,A ' ;y-y) 1 1+L A ' V [(2A+1)(2L'+ 1)] 1 / 2 , : C(L,l,A;y-y) C (L,v',l< ;pV03 " • 1-34 y Together, equations 1.32, 1.33, 1.34 applied to equation 1.31 can give the reduced matrix element for vector spherical harmonics: x E C(L',l,A';y, -y) C(L,l,A;y,-y) C(L,v,L';y,0). y 1.35 16 At this point, a comment can be made about parity restrictions. The parity of the vector spherical harmonic is given by the integer A as (-O^^. The electromagnetic wave functions are constructed from those vector spherical harmonics of the same parity; thus any discussion of the parity of the vector spherical harmonics applies also to the electromagnetic functions. In the sum in equation 1.35, changing the sign of the index \i changes the summand by, at most, a sign. The change is (-)^ * + + V . Thus, i f A+A'+v is an odd integer, the sum, and of course the whole matrix element vanishes. gamma decay. Let us investigate the available parities in nuclear When the i n i t i a l and final nuclear states have definite parity, as they do in most cases, a l l the allowed multipoles of the gamma ray must have the same parity, since electromagnetism conserves this quantum number. This means that A and A' have the same parity and that v must be even. The restriction to even values of v means that W(6) con- tains only even ordered Legendre polynomials; that i s , only even powers of cose. As a consequence, W(ir-0) = W(0); the gamma rays exhibit reflec- tion symmetry. However, there exist a few nuclear states of mixed parity. In these situations, odd values of v can arise from the interference of multipolarities of opposite parity. See, for example, Krane et^ a l . (1971). To complete the derivation of the reduced matrix element, we note that i t is zero when the parity condition is not met, (that i s , when A+A'+v is an odd integer) and we shall limit our further discussion to those cases in which the parity condition holds. The summand in equation 1.35 i s thus symmetric i n y , so that only the u=0 and u=l terms need be considered. 17 We are now ready to look at the matrix element of a spherical harmonic between two electromagnetic wave functions. First we shall write equations 1.4 in a more compact form using the following algebraic expressions "for certain Oehs#iQ3orctan cco;e£:'£rd.rents: C(L,0,A;0,0) = 6 1.36a. X L (A=L) C(L,1,A;0,0) = 0 = [l^-] 172 = - [2LTT] 1.36b 1.36c (A=L 1) + 1/2 ( X=L - LJ • 1 - 3 6 D With these expressions we may summarize equations 1.4 as L+o A ..(p) = Z i C(L,p,L-s;0,0) h T ._ IM ^ ._ , ^ >r> L+s L L+s M s=-l + K 1 p J r T T 1.37 p = 0 for magnetic multipoles = +1 for electric multipoles. The spherical Hankel function, of the f i r s t kind, used in the above equation tends to an outgoing wave at distances large compared to the wavelength. The asymptotic expression i s iKr \CKr) - i " w . 1.38 By including the radial dependence in the normalization, one may replace by i ^ q t i o n 1.37. ne U a The reduced matrix element ( L ' | | v | | L ) can now be found from equations 1.35 and 1.37. It contains two sums of the form in equation 1.37, and one of the form in equation 1.35. 18 Let us now examine one of the sums over s for u=0 in order to show that i t vanishes. Z Consider (-) i S L + P i L S ~ C(L,p,L-s;0,0) C (L, 1, L+s ; 0,0) . 1 This sum is over s=l,0,-l. If p=0, only s=0 contributes (by equations 1.36a) and the second coefficient vanishes, by equation 1.36b. If p=l, only s=+l and -1 contribute; the summand is antisymmetric between these two values because of the factor i The u=l term may now be examined. , causing the sum to vanish. The sums over s have this form: 2 (_) i P i'L-ss s L + 1 c(L,p,L-s;0,0) C (L,l,L+s"; 1,-1) = Z i s p + S _ 1 C(L,p,L-s>;0,0) C(L,l,L+s;l,-l) . 1.39 To evaluate this, one needs the Clebs;chsGoManeexpEessipns in equations 1.36, as well as the following expressions: X=L C(L,1,X;1,-1) = 1//2 = t2C2lTlT] [2 j TT T In 1.40a 1/2 A = L + 1 ]1/2 A = L L+ " 2 1 • 1 ' A 4 0 b O c equation 1.39, i f p=0, then s=0, and the sum, by equations 1.36a and-1.40a, i s -i//2. If p=l, then s=+l and -1, and equations 1.36c, d and 1.40b, c, show that the sum i s s t i l l -i//2. Because the matrix element contains the product of two such sums, one a complex conjugate, and because 19 the u = -1 term equals the u = 1 term, the result of a l l sums i s unity. The reduced matrix element can now be written. a'lMlD - t (^i)(2ui) /2 ] 1 c ( l f V f l .. l i 0 , 1.41 (-) L + 1 (2L 1) 1 / C(L,L',v;l,-l) 2 + That equation 1 . 4 1 has the desired normalization is seen by recalling that C ( L , 0 , L ; 1 , 0 ) = 1 . The final form of the F coefficient i s found by putting equation 1 . 4 1 into equation 1 . 3 7 . FV(L,L') = * x 21 [(2L 1)C2L' 1)(2I + W C I J . L . I J . L ' ^ . V ) + 1 + 1)] 1 / 2 C(L,L',v; 1,-1) . 1.42 Equation 1 . 4 2 is valid for both parity conserving and parity violating situations. When the electromagnetic multipoles of angular momenta L and L' have the same parity, v must be even; i f opposite parity, v must be odd. We haye now derived the angular distribution formula, equation 1.8,in terms of the populations of the i n i t i a l state of spin 1^. But i t is quite normal that one does not know these populations directly, but rather those of an earlier state in the decay. In figure 1 , i t is the state of spin I Q that has been oriented; the populations of the state of spin 1^ must be deduced from the populations of the f i r s t state. Let B' be the B coefficient for the spin I, state. v coefficient which appears in equation 1 . 2 6 . R 1 It i s this We wish to find B^ in terms of B^ the orientation parameters for the spin IQ state. 20 Let aCM^) be the relative populations of the i n i t i a l state, and a'(M ) those of the intermediate one. <I M'lpll M > 1 o o1 o o = a(M o The i n i t i a l density matrix is now ) 6., . . . MM' oo . 1.43 The i n i t i a l nuclear state 11 M > decays to a combination of 1 oo intermediate nuclear and radiation states with the same angular momentum properties as the i n i t i a l state has. 11 M > o o 1 Z aT C(I,,L ,1 ;M. ,M-M1) 1| I.M. >\$ M M > L ^ 1' o o' 1' o 1 1 1 ' L ,M -M, M1 o' o 1 T E . 1.44 o So far, there has been no difference between what has just been stated and what was stated at the beginning of the angular distribution derivation. At this point, we ask, not for the probability that the radiation be at a certain position, but for the probability that the nucleus is in a particular state |l^M^>, regardless of the radiation. This projection operator is LM The probability of finding the nucleus in the intermediate state I1M1> is 21 a'CMj) = Tr[p ^(I^)] = »„E <IoMo Ipllo M' ><oIoM' I I.M. > | <LM > j ..><1<)1' )... I <I M 11 M > o ' l l LM o o | K | T |Y 1 M M' o o L M = E M L aCMplaJ^CU^L , I^Mj ,M -M ) | Q 2 1 1 . 1.46 o Since equation 1.46 is a sum over pure multipoles, we may investigate each one separately and do the sum afterwards. For the multipole L , Q the orientation parameter^, B^, calculated from the intermediate populations is B' = ( 2 v + l ) E 1/2 M = (2v+l) CCI^v^jM^O) l E 1/2 a'CMp CCI^v.I ;M ,0)|C.(I ,L , I ; M Q M..M 1 o ,M -M ^ | 2 a(M ) Q . 1.47 The Racah coefficient identity in equation 1.23 may be applied to the sum over M^. E M „ I..+I -L +v C(I ,v,I ;M 0)|C(I ,L,I ;M ,M -M )| = (-) ° [ (2^+1) (2I +1) ] 1 Z 1 1 1J 1 o 1 o 0 1 Q l x C[I ,v,I ;M ,0D o o W C I o 0 , I L F I 0 , I I ; L o , v ) . With the help of equation 1.48, B^ can be written B* = ( 2 v + l ) v ^ J = U B V V 1/2 K J U E C(I ,v,I ;M ,0)J a(M ) v o' ' o' o o M o 1.48 A / 22 where U = (-) I +1 -L ° 0 -v ,„ [(2I +1)(2I +1)] 1 O 1 7 W(I L V O 1 ,1 -L , v ) . 1 0 1.49 O If several multipoles contribute to the transition, the resultant U is an incoherent sum over their individual coefficients, by virtue of equation 1.46. This completes the derivation of the primary terms in the angular distribution formula. Of practical interest is a further term to correct for the f i n i t e solid angle subtended by the radiation detector. This geometric factor is given in the next section. 1.4 Solid Angle Correction Because the radiation detectors do not look at a single value of 8 i t is necessary to integrate W(6) over the face of the detecting surface. This process involves only the Legendre polynomials, replacing each P^(cos6) by g^P^fcose), where g^ is a function of the geometry. Figure 2 illustrates the notation to be used for a circular detector. Let 0 be the origin, the source of radiation, C the centre of the detector face and A the point on the upper circumference of the detector face such that A is in the plane containing the z axis and the line C 0. has length R. The latter The angle between C 0 and the z axis is 0 ; the half angle subtended at 0 is a. Any point P on the detector face can have the co-ordinates § and 8, where <j> is the angle between P C and C A, and 3 the angle between P 0 and C 0. The Legendre polynomials must be integrated over a l l points P. To do this, one must know the angle 6' between the line P 0 and the z axis. The length of P 0 is R/cosg, and the length of P C is R tan 3. The z component of the point P is R cos8 + (R tan3)cos<j>sine. Dividing the 23 Figure 2 Solid angle correction 24 1 z component by the length of the line P 0 gives the cosine of 0 . cos0' = cosgcosS + sin3sin0cos(j>. 1.50 By the addition theorem of Legendre polynomials, when equation 1.50 is valid, one may write 1 Pv (cose ) = Pv (cosg)P (cose) +2 v 1 m v J r E y ' l ' , P"v (cosg)Pv(cose)cos m*. m = 1, fv+m)! m 1.51 Averaging the Legendre polynomial over the detector face yields a 7 - J •' d<j> ,/ d3 siri3 P (cose1) = P (cos0) =^ j P (cos3)sin3d3 , a 2TT O J O J 0 *> a where Q = 2ir / sin3d3 0 " = 2Tr{l-cosa)'. 1.52 The higher terms vanish because cos m<t> integrates to zero. Equation 1.8 is corrected for solid angle by rewriting i t as W(0)' = E BV UV FV gV PV (cos0) v V 1 a g = -j-, r- / P (cos3)sin3d3. v (1-cosa) ' v & v 1.53 J For second and fourth order Legendre polynomials, g2 = j cosa (1+cosa) 1 2 g^ = g- cosa (1+cosa)(7cos a-3). 1.54 25 1.5 Methods of Orienting Nuclei Having established that a collection of nuclei whose spins are not randomly oriented can produce decay products with a non-uniform angular distribution, we must now illustrate how the nuclear spins may be oriented. Before doing so, however, i t is necessary to explain some terminology. An ensemble of nuclei is said to be oriented i f the relative populations of the magnetic sublevels are not equal; that i s , i f a(M) ^ (21+1) ^. Furthermore, i f a distinction exists between up and down, that i s , i f a(M) ^ a(-M) for at least one value of M, the nuclei are said to be polarized. On the other hand, i f the nuclei are oriented, but with up and down symmetry (i.e. a(M) = a(-M) for a l l M), then the system is said to be aligned. In the case of alignment, the odd orders of the orientation parameters vanish, since they are composed of the expectation values of odd powers of M. It is often useful to note whether a method of orientation produces polarization or only alignment. In this work static nuclear orientation is employed; that i s , the nuclear orientation degeneracy of the Hamiltonian is broken and, following thermal equilibrium, the relative populations of the nuclear substates are given by the Boltzman distribution: a(M) « exp(-E(M)/KT). In this scheme, "reorientation" processes must exist to allow thermal equilibrium to arise. Nuclei caiupossess the.follgwingcelectromagneticemoments: electric charge, magnetic dipole moment (I >_ 1/2), and electric quadrupole moment (I >_ 1). The interaction between the charge and an electric field cannot orient the nucleus; the interaction between the magnetic moment and a magnetic field can polarize; and the interaction between the electric 26 quadrupole moment and an electric field gradient can align. Moments of higher order are "usually negligible. Let us examine the contribution from these moments to the Hamiltonian. Reviews of these methods are to be found in Blin-Stoyle and Grace (1955) and DeGroot, Tolhoek and Huiskamp (1965). i. Brute Force The most direct way to orient nuclei is to apply a large magnetic field. The additional term to the Hamiltonian is ->• ->E = -H-y = - ^ H-l = -Y HM 1.55 where H is the magnetic f i e l d , I the nuclear spin, and y the nuclear magnetic moment. The last line of equation 1.55 assumes that the quantization axis is the axis of the magnetic f i e l d . With this interaction the nuclear populations are a(M) = Z" 1 exp(BM) Z = E exp(BM) M = where 3 = yH/IKT. sinh 3(1+1/2) sinh;.-> 3/2 1.56 27 To achieve appreciable orientation B must be of the order of unity. -24 For magnetic moments of the order of the nuclear magneton (5.05 x 10 7 ergs/gauss) one requires H/T to be 3 x 10 gauss/kelvin. At room tempera9 ture, H would have to be about 9 x 10 gauss, a rather large value. Even -2 5 at low temperatures (T ~ 10 K) the field must be about 3 x 10 gauss. Dabbs et a l . (1955) used this method to orient ^^Tn. ii. SpincHaffii^tpnian - So far i t has been assumed that the atoms are free. This is not true in general, for in crystals the atoms may feel very strong forces. For a paramagnetic crystal with axial symmetry about the z axis Abragam and Pryce (1951) have derived the effective spin Hamiltonian X = [g„ H z Sz + g,(H S + Hy Sy )] + D[S -3 is(S+l)] + ASz z I + B(S I + yS Iy ) II i xx z xx 2 L 6 6 J + Q[I 2 L - j 1(1+1)] - ^-H-T . 1.57 The electron states are represented by an effective spin S, such that 2S+1 equals the number of spin levels. The f i r s t term is the interaction of the electrons with the external magnetic f i e l d , which has different strengths along and perpendicular to the z axis; the D term is the splitting of the electronic levels by crystalline fields; the A and B terms are the hyperfine coupling of the electronic and nuclear spins; the Q term is the nuclear level splitting caused by the crystalline electric field gradient; the last term is the external magnetic field acting on the nuclear magnetic moment. Different methods of orienting nuclei use various terms in this Hamiltonian. 28 iii. Gorter-Rose Method The method proposed by Gorter (1948) and Rose (1949) employs the hyperfine f i e l d in a paramagnetic material. A f i e l d of a few hundred gauss i s needed to saturate the electronic magnetic moments through the g terms of equation 1.57; then, the nuclei are made to follow the electrons by the A and B hyperfine terms. The advantage of this procedure i s that one can produce effective fields at the nucleus a thousand times stronger than the external f i e l d . Ambler et al_. (1953) have done the f i r s t experi- ments with this method. iv. Bleaney Method The electronic spins can be ordered by internal as well as by external means. In the method of Bleaney (1951) the D term in equation 1.57 splits the electron levels, and thus the nuclear levels through A and B. No external f i e l d is needed; only low temperatures are required. The f i r s t successes of this method were performed by Daniels et^ al_. (1951) and Gorter et_ al_. (1951). v. Pound Method Pound (1949) suggested using the nuclear quadrupole interaction. Like the Bleaney method, no external f i e l d i s needed, but one does require large electric field gradients (^ 10 e.s.u.). Uranium and neptunium have been so oriented by Roberts, Dabbs 131 and Parker (1958), and Shirley (1960). I has been aligned by Johnson, Schooley and 29 vi. Ferromagnets In ferromagnetic material, the strong exchange force between the electrons causes their spins to be completely polarized along some axis, say the z axis. One may then set and to zero. Marshall [1958) has shown that the hyperfine interaction behaves like an effective magnetic field at the nucleus given by H eff = -— AT p S+H. 1.58 As a result the nuclei within a domain are automatically polarized when the temperature is sufficiently low. method by orienting ^Ce Grace et a l . (1955) proved this in a single crystal of cobalt metal. In multidomain samples, an external f i e l d is s t i l l necessary to magnetize the whole specimen so that a single hyperfine field is present. Fields of a few kilogauss are required to achieve this in iron, whereas over ten kilogauss may be necessary to saturate cobalt. This method, which is the one employed in the present work, is not restricted to magnetic atoms. Samoillov et a l . (1960) oriented diamagnetic atoms in ferromagnets. The nucleus of a non-magnetic atom in a ferromagnetic host sees a hyperfine field caused by the polarization of the conduction electrons. The systematic behaviour of the field for different impurities in iron has been given by Shirley and Westenbarger (1965). Calculations of the field have been done by Campbell (1969). Using this effect, the number of nuclei which can be oriented is greatly increased. One is not restricted to using the hyperfine methods only for magnetic atoms; i t may be used for any atom which can be introduced into the ferromagnetic host. 30 vii. Antiferromagnets The case of antiferromagnetism is similar to that of ferromagneT tism, except that the exchange forces align rather than polarize the electron spins. Ear sa ssingle crystal "with;es^mmetry.tlow.er than cubic no : applied field is heeded; -thev electron" and nuclear spins align along a crystalline axis as the temperature is lowered. 54 Daniels, Giles, and LeBlanc (1961) observed orientation of Mn and Co in crystals of MnCl -4H 0, CoCl -6H 0, MnSiF -6H 0, and 60 2 2 2 2 6 2 Co(NH ) (S0 ) -6H 0. 4 2 4 2 2 viii. Dynamic Nuclear Orientation There exist methods for producing non-uniform populations of the nuclear magnetic states without low temperatures. In dynamic nuclear orientation, optical and microwave radiations excite atomic transitions; the subsequent emission of radiation by different channels can create nuclear orientation. For a complete discussion of these methods see Jeffries (1963). 1.6 Low Temperature Production A description of the low temperature apparatus is given in a later chapter; this section will only give an outline of the principle used for cooling. Preliminary cooling down to one kelvin is accomplished by liquified gases. Liquid nitrogen can cool to 77 K and liquid helium at one atmos- phere to 4 K. Pumping the liquid helium to a pressure of 1 mm of mercury 31 further reduces the temperature to 1 K. At this point, the system is ready for the final plunge. The last step is an adiabatic demagnetization of a paramagnetic salt. At this low temperature, the entropy remaining in the salt i s mainly that of the ionic spins. Since entropy is a measure of disorder, i t may be reduced by ordering those spins in a magnetic f i e l d . This process can be illustrated in figure 3. Initially, with no external f i e l d , the electronic spin states are only slightly split by the inter-atomic interactions; the splitting is much less than KT, and the states are very nearly equally populated, as indicated in figure 3 by the thickening of the lines. When the external field H i s applied, the levels are widely split; after thermal equilibrium i s re-established, the lowest levels are the most heavily populated. Thermal contact to the 1 K temperature bath is broken and the f i e l d i s removed adiabatically. In an adiabatic process, the entropy remains constant; that i s , the populations do not change as the level splitting is decreased. temperature is-produced : " f sp 7 A low '.J picduc^c The final temperature after demagnetization is determined by the inter-atomic interactions. It is possible to reduce these forces by separating the magnetic ions by non-magnetic ions and water of hydration. Some common salts used are ferric ammonium alum [FeNIfy(SO4)2"I2H2O], manganous ammonium alum [MnCNfy)2(SO4)2"ol^O], chrome potassium a'iumi'fCrK'CSO^) ^.lZH^O]and r cerium ^magnesium nitra'tecJge^g'j'CNOp • 2'4Hg'Of •:. I- C , , A 2«.:-o Further details regarding this and other methods of cooling can be found in the manual for low temperature physics by White (1968). 32 Figure 3 Level splittings and populations during adiabatic demagnetization 33 1.7 Other Experimental Methods In the standard method of orienting nuclei the anisotropy of the emitted radiation is observed as a function of temperature and/or angle. Equation 1.8 contains the temperature only in B terms, which are functions of 3 = uH^^/IkT. The hyperfine interaction yH ^^ may be so obtained. g In addition, once the B terms are known i t is possible to evaluate UF, which can give information about spin assignments and gamma ray multipolarities in the nuclear decay scheme. Other methods which can determine these quantities are: i. Angular Correlations An experiment with a basic idea similar to that of nuclear orientation is angular correlation. When two radiations are emitted in sequence the angle between them has a distribution W(0) = E F (1) F (2) P (cose) 1.59 V where the F terms are the same ones used in nuclear orientation. Spins and multipolarities can be found in this way. The advantage of this method is that i t is performed at room temperature, .i.^ i r. : A review of this method is given by Frauenfelder, and Steffen (1965) ."v.. . .- • ii. ' • . Perturbed Angular Correlations If in a sequence of nuclear transitions the intermediate state is long-lived, the precession of the nucleus in the hyperfine field can 34 be observed in the time dependence of equation 1.59. A periodic re-orien- tation combined with the exponential decay must be included in this equation to give the correlation as a function of the time between transitions. The precession frequency so determined gives the hyperfine inter- action. A f u l l treatment of this technique can be found in Karlsson ejt al. (1964). iii. Internal Conversion A nuclear gamma transition can occur not only by emitting electromagnetic radiation (the gamma ray) but also by emitting inner orbital electrons. This process is called internal electron conversion. The ratio of converted electrons to gamma rays does not depend on nuclear matrix elements but only on the transition energy and multipolarity. The theory of internal conversion is discussed by"Rose (1965) and experiments are reviewed by Ewan and Graham (1965) . iv. MBssbauer Effect The recoil-free absorption of gamma radiation (Mo'ssbauer , 1965) can be used to determine the hyperfine splitting of the nuclear levels which is seen as structure in the gamma peak. v. NMR Nuclear magnetic resonance observes the hyperfine splitting by the absorption of electromagnetic radiation at the resonant frequency •'di = yH ^/ftl. The f i r s t "bbse-r-vation'Of-.this'effect: dn ferromagnets was reported by- Portis and Gossard -(I960) who dejtee;teed an\ NMR. signal in rr cobalt.- The "signal".was- very strong duetto enhancement': in the domain walls".- - l •/rf,-.. ••-;>.::- - ... 35 Impurities in ferromagnets can be studied by NMR, but the concentration must be > 0.1% i n order to observe a measurable signal. vi. NMR/Nuclear Orientation Trace impurities in ferromagnets can be studied by a combination of NMR and nuclear orientation. A small coil i s placed around the sample so that i t produces an r.f. f i e l d perpendicular to the orientation axis. The thermal equilibrium distribution of nuclear populations i s destroyed by putting in radio-frequency power at the hyperfine resonant frequency. The observed effect will be a decrease in the anisotropy. This will be restored by spin-lattice relaxation upon removal of the r . f . field. This method has the advantage of the high accuracy found in resonance experiments, this accuracy being limited by the line width observed. Typically, hyperfine fields can be measured to MD.03%. The resonant destruction of nuclear orientation was f i r s t observed by Matthias and Holliday (1966). Some problems with the method are r.f. heating and the preparation of samples with the thickness of the skin depth. 36 CHAPTER TWO EXPERIMENTAL APPARATUS AND PROCEDURES 2.1 The Cryostat The low temperature apparatus i s described in detail by Gorling (1970). A brief description of the'cryostat which i s illustrated in figures 4 and 5, i s given here. Liquid nitrogen at 77 K and liquid helium at 4 K were contained in two glass dewar vessels. Immersed in the liquid helium were the superconducting magnets and the brass jackets of the inner cryostat. Two superconducting solenoids were necessary. The primary one provided the large field necessary for cooling by adiabatic demagnetization. It was made from niobium-titanium wire wound onto a solenoid of outer diameter 10 cm., inner diameter 5 cm., and length 18 cm. It provided a central field of 0.8 kG per amp and would take up to 62 amps, or 50 kG. The secondary magnet produced the field necessary to saturate the magnetization.of the ferromagnetic sample. It was constructed from superconducting niobium-titanium wire wound around a core of diameter 3.5 cm. and length 3.5 cm. The number of turns was 1700. Its maximum current was 25 amp (or 9.8kG). The inner cryostat consisted of three brass jackets the innermost of which contained the paramagnetic salt p i l l s and the radioactive specimen. The middle jacket held liquid helium which was kept at a temperature of 1 K by continuous pumping. This liquid was insulated from the 4 K helium by the vacuum (10 ^ torr.) in the outer jacket space. 37 PUMP LINES A | i LIQUID HELIU M ,\ fcT] RETURN LINE f ) OUTER DEWAR INNER DEWAR LIQUID OUTER J A C K E T I LIQUID NITROGEN I SCALE: CM I 1 L 0 5 10 15 Figure 4 Outer c r y o s t a t HELIUM 38 INNER JACKET PUMP LINE OUTER JACKET PUMP LINE BATH PUMP LINE N E E D L E VALVE OUTER J A C K E T SPACE 1K BATH INNER JACKET SPACE GUARD PILL SUPER-CONDUCTING COOLING M A G N E T M A I N PILL COPPER HEAT CONDUCTOR POLARIZING MAGNET RADIO-ACTIVE SPECIMEN SCALE: I 1 0 2 CM 1 4 l_l 6 8 Figure 5 Inner cryostat 39 High frequency vibrations could be a source of heating, either through mechanical flexing of the supports, or through currents induced by a conductor moving in a magnetic f i e l d . taken to reduce vibrations. The entire cryostat was supported on a dexion frame sunk into a cement base. the For this reason efforts were resonant frequency of the frame. The weight of the cement lowered To isolate the assembly from vibrations from the floor, the frame rested on springs or on a foam rubber pad,bothliof /whichfhad-equal effect. 2.2 Paramagnetic Salts and Thermal Contact i. Main P i l l In figure 5 i t may be seen that the innermost jacket of the inner cryostat contained two paramagnetic salt p i l l s . The primary p i l l was made of chrome potassium alum or cerous magnesium nitrate. The latter salt allowed lower final temperatures to be reached but the former had a higher heat capacity. With chrome potassium alum one expects to reach a temperature of ^12 mK, the N6el point for this salt, but with cerous magnesium nitrate one can obtain temperatures down to ^2 mK. The lowest temperature that was reached with the present system using the nitrate was 8 mK. ii. Guard P i l l and Shield The secondary paramagnetic salt p i l l consisted of manganous ammonium alum. The purpose of this p i l l was to reduce heat conduction down the supports to the primary p i l l . In later experiments a guard 40 shield (not shown in figure 5) was introduced. (0.1 mm.) This shield was a thin copper sheet wrapped around the p i l l and sample assembly such that i t was in contact with the brass wall of the guard p i l l but did not touch the main p i l l or the sample. The role of the guard shield was to adsorb the helium gas remaining in the inner jacket during demagnetization and to transfer the heat of adsorption to the guard p i l l . The two p i l l s were separated by a nylon tube. Short bits of nylon wire projecting from the guard shield kept the assembly from touching the wall of the jacket. The whole assembly was suspended from the top of the jacket by two 1.5 cm. lengths of stainless steel tubing, 0.56 mm. in diameter. iii. Thermal Contact Because a magnetic field of several kilogauss had to be applied to the radioactive sample to achieve nuclear orientation, and because the low temperature was produced by the removal of a magnetic field from the paramagnetic salts, the sample and the salts had to be spatially separated and yet remain thermally connected. This was accomplished by a folded copper sheet with one end embedded within the main p i l l and the other end extending 9?cm. below the p i l l . attached to the lower end of the copper. f o i l , i t was soldered to the copper. The radioactive sample was If the sample was a metal If i t was a crystal, i t was smothered in Apiezon N grease and wrapped within part of the copper f o i l ; needless to say, in the latter case a much poorer thermal contact was obtained. iv. P i l l Preparation The paramagnetic salt p i l l s were made by mixing finely ground salt with an equal mass of glycerine. The glycerine served to make thermal contact with the copper and to form a glass-like solidaat low temperatures. The folded copper fins f i l l e d the entire length of the main p i l l . The salt was packed between the folds to maximize the contact area, which 2 was about 900 cm . The total amount of salt and glycerine was about 120 grams. Outside the p i l l , the fins were pressed and glued together. The surfaces were previously covered with an epoxy paint to reduce heating from eddy 'Currents. 2.3 Low Temperature Procedure After the cryostat had been pre-cooled to 4 K by liquid helium the large superconducting solenoid would be charged to 50 amps (40 kilogauss) in five minutes. (the A superconducting wire across the magnet leads persistent switch) could be made normal (i.e. possessing resistance) by a heating wire around i t . While the magnet was being charged, the switch would be normal (an open circuit); to put the magnet into the persistent mode, the heater was turned off, the switch made superconducting, and the magnet short-circuited. After the power supply was turned off the current remained in the magnet but completed the circuit through the switch rather than through the power supply. The heat of magnetization in the paramagnetic salts was conducted away by helium gas. The pumping of the outer jacket space would be stopped, and the pumping of the 1 K bath would begin. By pumping the liquid helium in this space to 0.2 torr., the temperature was reduced to 1.05 K. Helium gas was added to the inner jacket space as required to maintain thermal contact with the bath. After the paramagnetic salt p i l l s had ceased to cool, as indicated by a carbon resistor embedded in one of them, the inner jacket space would be pumped down to 10 ^ torr. pumping was stopped. When this had been achieved, the The magnet was discharged by restoring the current from the power supply to 50 amps, making the persistent switch normal, and reducing the power supply current to zero in ten minutes. The 7 pressure in the inner jacket space normally f e l l to about 10 -6 to 10 torr. The polarizing solenoid was now charged and maintained i n the persistent mode. The counting of the gamma rays could then begin. 2.4 Gamma Ray Detection Once the radioactive nuclei had been oriented one could observe the anisotropic emission of gamma rays. To do this, devices sensitive to such radiation were positioned about the source. i. Types of Detectors Two types of gamma ray detectors were used. Initially two sodium iodide detectors were employed because of their high efficiency, but as the gamma ray spectra became more intricate the higher resolution of lithium-drifted germanium detectors became necessary. The sodium iodide detector consisted of a crystal of Nal(Tl) which was 2 inches in diameter and 2 inches long. It was in optical 43 contact with a photomultiplier (RCA 6342 A) biased at 1100 volts. Because this detector was extremely sensitive to magnetic fields, i t was shielded by mu metal and an iron tube along i t s entire length. r Such protection screened only the field of the polarizing solenoid; when the large magnet was charged these detectors were inoperable. The resolution photopeak of ^Co (full width at half maximum) of the 1.33 obtained in this system was 60 keV. of the two lines in ^Co MeV The separation is 159 keV; thus these lines are not quite fully separated. The Ge(Li) detectors consisted of a crystal of lithium-drifted germanium kept at the temperature of liquid nitrogen at a l l times. f i r s t of The these detectors had a trapezoidal crystal with an active area 2 of 11.2 cm and a drifted depth of 1.1 cm. 1000 volts, negative. detector was 3.9%. The bias on the crystal was 2 The relative photopeak efficiency for this Ge(Li) The peak-to-Compton ratio was 9 to 1, and the resolu- tion of the 1.332 MeV peak was 5.0 keV. The second Ge(Li) detector had a cylindrical crystal. 3.8 cm. in diameter and 4.1 cm. long. It was The bias was 3000 volts, positive. In this case, the relative photopeak efficiency was 6.7% with a peak-toCompton ratio of 25 to 1 and a FWHM of 2.4 keV. The relative photopeak efficiency of a Ge(Li) detector is defined in this work as the efficiency relative to that of a 3 in. x 3 in. Nal detector for the 1.332 MeV peak of 60co when the source is 25 cm. from the detector. 44 ii. Handling of Detector Output The pre-amplifiers of the detectors produced a voltage pulse _7 with a rise time of the order of 10 of 10 ^ sec. sec. and a decay time of the order The height of the pulse was proportional to the energy of the detected gamma ray. Sorting the pulses according to their heights therefore produced the energy spectrum of the radiation. Figure 6 shows the path of the detector signals. The pre- amplifier pulse was shaped and enhanced by a linear amplifier (Canberra Industries model 1410). The height of the output from this unit ranged from 0 to 10 volts, and was also proportional to the energy of the original gamma ray. Pulse shaping was achieved by a choice of appropriate differentiation and integration time constants to produce a bipolar output. For the Ge(Li) detectors, the best resolution was obtained using 2 u sec. for the f i r s t differentiation, 1 u sec. for the integration, and 7 u sec. for the second differentiation. For the Nal detectors, the bipolar pulse produced by double delay lines was used. The output from the linear amplifier was analyzed in two ways: by single channel analyzers (SCA) and scalers, and by a multichannel analyzer. The single channel analyzers produced a logic pulse of fixed height whenever the height of the input f e l l between two pre-set levels. The logic pulses were counted by a scaler (Canberra Industries model 1473). The linear amplifier signal could be sent to several SEAs, each set to select a different gamma ray. The resulting array of scalers was controlled by a timer/sealer (Canberra Industries model 1492) which stopped them after a pre-determined time. The scalers were used during an experiment to monitor the most interesting gamma rays. 45 DETECTOR PRE-AMP TIMER/ SCALER TIMING SCA LINEAR AMP SIGNAL • FROM SECOND LINEAR AMP DUAL INTEGRAL DISC. SCALER ROUTER JT. MIXER/ REJECTOR MULTI-CHANNEL ANALYZER Figure 6 Flow diagram of detector output 46 The more important path for the signals was the one leading to the Nuclear Data model 2200 multichannel analyzer. 512 channels and had a digitizing rate of 50 MHz. This unit contained With i t the entire spectrum could be recorded by pulse height analysis. Since two detectors were normally used during an experiment i t was desirable to be able to accumulate both spectra simultaneously in separate sections of the memory. Three additional units were necessary: a dual integral discriminator (Nuclear Data 531), a mixer/rejector (Nuclear Data 521) and a router (Nuclear Data 538). The analogue signals from the two linear amplifiers were presented to the dual integral discriminator. When either of the inputs exceeded a pre-set level, a corresponding logic pulse was sent to the router. The router generated signals for the analyzer to determine which half of the memory was to receive the analogue pulse. signals to the mixer/rejector. The dual discriminator also sent two One was the sum of the signals at the two analogue inputs; the other was the sum of the logic signals which went to the router. The mixer/rejector could therefore determine whether the analogue signal i t received was a coincidence between the two detectors. If so, that signal had to be suppressed. went to the analyzer input. The outputs of the mixer/rejector One Was the analogue pulse, which was de- livered to the analogue-to-digital converter, and the other was a logic pulse fed to the coincidence input of the analyzer. If this logic pulse was missing, the analogue signal was .r.ejiect'ed. The multichannel analyzer would accumulate in this way for a pre-set time. Afterwards the contents of the memory could be displayed on an oscilloscope and/or put onto paper tape in 40 sec. by a high speed Teletype punch. The memory would be erased and a new count begun. The 47 paper tape was then the only record of the last count. After the experi- ment, the contents of the paper tapes were transferred to magnetic tape and processed by the IBM 360 computer. iii. C Counting Procedure The two most convenient angles for observing the change in gamma ray intensity when the nuclei were oriented were along the axis of orientation (the axial direction) and perpendicular to i t (the equatorial direction) . Detectors placed in these positions measured W(0) and W(TT/2) respectively. Since the orientation axis was a vertical one, the axial detector was placed beneath the cryostat and the equatorial detector beside i t . Since the two Nal detectors were identical i t did not matter at which station either one was put. detectors. That was not true for the Ge(Li) The axial detector was 20 cm. away from the source while the equatorial was only 10 cm. away. It was therefore advantageous to put the more efficient detector at the axial position. The angular distribution If(8) is normalized to unity at high temperatures. One therefore measured W(0) and W(ir/2) by dividing the gamma ray intensity by i t s value when the specimen was warm. During the time preceding the demagnetization, as many counts as possible would be acquired so that after demagnetization, when the specimen was cold, the gamma ray intensities could be properly normalized. _2 After the specimen had been cooled to 10 K and the polarizing f i e l d had been applied, cold counts would be taken on the multichannel analyzer and on the single channel analyzers. The counts from the latter 48 were recorded in a log-book. A good indication of W(0) and W(TT/2) could be immediately provided by the single channel analyzers. Of special interest were the intensities of the ^Co peaks, since they acted as an in situ thermometer, as will be described. But because the scaler counts included not only the gamma ray, but also any background of the same energy, the final results were those given by the computer analysis at a later time. The procedure followed while the specimen was cold depended upon the type of experiment being performed. When only a small anisotropy was expected counts would be acquired for as long as the temperature was below a desired level. Count times ranged from 1000 to 4000 sec. After the specimen temperature had exceeded the limit (three to six hours after demagnetization) the system would be warmed to 1 K and another magnetization would be initiated. However, i f the experiment consisted of observing changes in the anisotropy as a function of some other parameter (applied magnetic field or temperature) several counts of 200 to 400 sec. would be taken at each setting. the specimen was too warm to be of use. Counting would proceed until Remagnetization and often a retransfer of liquid helium would be necessary. Warm counts would be acquired again during the few hours prior to the next demagnetization. They would also be accumulated at the very end of the experiment, provided that the liquid helium level was above the source. This was because the liquid absorbed a few percent of the gamma rays, and a noticeable increase in the equatorial intensity would be observed when there was no liquid between source and detector. 49 2.5 Paper Tape Analysis After the completion of a run the data consisted of a l i s t of scaler counts in the log-book and a great deal of paper tape. The infor- mation stored on the tape was f i r s t converted to magnetic tape for analysis by a computer. i. Spectrum Analysis The scaler counts taken during the .experiment gave the f i r s t approximation of the observed anisotropy. As has been mentioned the background included in the count could produce errors. A constant background tended to decrease the difference between the warm and the cold counts. But the background was not necessarily constant. Associated with the photopeak of a gamma ray is a broad Compton plateau of lower energy exhibiting the same anisotropy as the peak. Also, gamma rays emitted in one direction could be scattered, with reduced energy, into a different detector. This background would possess the anisotropy pertaining to the original direction of the gamma ray. The total background under a peak therefore depended in a very complicated way on the peaks of higher energy. The main component, however, was the Compton contribution. The anistropy of this unwanted segment could overwhelm that of the peak of interest, especially when only a small effect was present in that peak. In this event the scaler readings were not even a reasonable f i r s t approximation. The separation of the peak from the smooth background upon which i t stood was performed by our computer program SPECTAN reading the multichannel analyzer output. The program could also detect and correct 50 displacements of the peak positions. For nuclei of short lifetimes a correction for the decay could be made. By analyzing the warm counts f i r s t , the cold counts could be immediately normalized. The statistical error was calculated assuming that the number of counts in each channel had a variance equal to that count- (the. Poissonf distribution.),. Up to five peaks from one detector could be analyzed at one time. ii. Background Calculation The assumption used to find the background beneath a peak was that the background spectrum was linear over a distance of about twice the total width of the peak. This assumption was far more realistic for peaks taken with Ge(Li) than with Nal detectors because the peaks of the former were so much narrower. Since only relative intensities were required in this work, this procedure for analyzing spectra was found to be satisfactory in most cases. The method of calculation is illustrated in figure 7 . The f i r s t step was to sum the channels within the specified peak window which produced an uncorrected count equivalent to the scaler count. The second step was to skip several channels (the "gap" in figure .'7Q on either side of the peak window and to sum some more channels (the "width"). This second sum, divided by the number of channels contributing to i t (that i s , twice the "width") is the average background per channel between the two background windows. Because the peak window is symmetrical- ly located between the background windows, the background in the peak window is the average times the number of channels in the peak window. 51 • • • •« w e 1 1 D T H G A P 4-* WINDOW H J O 1 1 5 10 WINDOW GAP WIDTH L 15 s 7 = 2 = 4 TO 20 25 30 22 •I 22 UNCORRECTED COUNT YIN) Ns7 z 28 BACKGROUND = Y(N) 4/ 25 Figure 7 Background calculation Y(N) 22-7+1 2«4 52 The limits of the peak window, the number of channels in the "gap" and in the "width" were a l l variable. The corrected peak intensity i s the uncorrected count less the background. The variance of this quantity is the sum of the individual variances. These variances were easily estimated since the uncorrected count and the second sum both possess Poisson statistics. were taken as the sums themselves. The variances Since the background is proportional to the second sum, i t s variance equals this sum times the square of the proportionality factor. If a correction for the decay was to be made i t was applied to the background-free count. The time at which the spectrum had been taken was included on the magnetic tape, the h a l f - l i f e was given as input to the program, and the peak was increased to correspond to a time prior to the experiment. The variance was similarly adjusted. iii. Normalization The program SPECTAN would be given the spectra to be analyzed and was told the number of warm counts. It analyzed these f i r s t so that i t would have the average normalization count and corresponding variance before i t started the cold counts. The intensities of the peaks in the cold spectra were produced in the same way as in the warm ones, but a further step was added. intensity was divided by the normalization. The The statistical error in the normalized count would be found by using the rule that the relative variance of a quotient is approximately the sum of the individual relative variances. This assumes that the dividend and divisor are independent (which they certainly were). 53 iv. Shifting Peaks An additional feature of the program SPECTAN was the detection and correction of the peak position. Adequate results were obtained only i f the peak spanned tens of channels, for the window could be shifted by only an integral number of channels. This condition was met by spectra taken on a Nal detector, but usually not by those from Ge(Li) detectors. However, i t was more often the Nal detectors that tended to shift — .or::^ The program calculated the centroid of the peak window by giving to each channel a weight equal to i t s count. the f i r s t spectrum was used as the standard. The centroid of If a subsequent spectrum had the .centrbid--in-a;different channel, the peak window was shifted accordingly. The centroid was calculated for the new window and another shift was made i f necessary. The corrections continued until the cen- troid had the proper position within the window. There was of course the possibility that a divergent or oscillating situation would occur, thus a maximum number of attempts was specified beforehand. If centroid corrections were not desired, the maximum was put to zero. v. Program Output The computer print-out would show the results of these calculations. For the warm spectra one would be given the uncorrected and corrected counts, the variance for the corrected count, the peak window, and the number of centroid corrections. After the last warm spectrum the totals and averages would be produced. For the cold spectra one would be given the same data and, in addition, the normalization and 54 their standard deviation errors. If the output was to be subjected to further computation, some of these results would be written into the computer memory. For both warm and cold spectra the program would write the uncorrected count, the decay correction factor, the corrected count, and its variance. This data could be used by programs written for specific experiments. 2.6 Thermometry: ^Cobalt Experiments In nuclear orientation experiments the problem of measuring the specimen temperature is solved by using the anisotropy of some known radioactive material. Such an isotope i s ^Co dissolved in iron. The magnetic moment and hyperfine field are known, the decay scheme has been completely determined, and a large anisotropy exists at the desired temperatures. Additional advantages are that the isotope is reasonably long-lived (half-life of 5.25 years) and that the two gamma rays have identical anisotropics. This thermometer was f i r s t used by Stone and Turrell ((M)6'2). i. Decay Scheme The decay scheme for ^Co as given by Lederer, Hollander and Perlman (1968) i s shown i n figure 8. The ^Co nucleus decays by emission of beta particles to an excited state of ^ N i . Two gamma rays are promptly emitted as the ^ N i nucleus decays to i t s ground state. The energies of the gamma rays are 1.173 and 1.332 MeV, and both are pure E2 transitions. 55 Figure 8 Co decay scheme 56 Using equations 1.10 and 1.12 (or the tables of Yamazaki, 1967) we find that for both gamma rays U F 2 2 = -0.42056 U F 4 = -0.24281. 4 2.1 Anisotropy and Temperature n. The temperature dependence of the anisotropy depends on the hyperfine field of Co dissolved in iron and on the magnetic moment of the ^Go nucleus. The hyperfine field is given by Koster and Shirley (1971) as -287.7 kG and the magnetic moment is given by Shirley (1971) as 3.754 nuclear magnetons. The relative intensity of the gamma rays in the axial and equatorial directions can be calculated from equations 1.8, 1.9 and 1.56: 6 a(M) Z B v W(0) W(TT/2) uH/IkT Z exp(BM) sinh g (1+1/2) sinh B/2 (2v+l) 1/2 1 + BUF 2 2 I C(5,v,5;M,0)a(M) M 2 + BUF 4 4 4 2.2 57 In figure 9 are plotted W(0) and ture. At 12 mK we have W(0) = 0.60 and as functions of tempera- W(TT/2) W(ir/2) = 1.16. This anisotropy is large enough to serve as a good thermometer in this temperature region. iii. Sample Preparation The f i r s t experiments with this apparatus were carried out on ^Co in iron as a test of the system. A carrier-free solution of ^CoCl,, with an activity of 1 mCi was obtained from Atomic Energy of Canada, Limited. Solution containing 2 uCi was placed on an iron f o i l of dimensions 1.0 cm x 1.0 cm x 0.03 cm. After the solution had evaporated to dryness, the iron was heated at 700°C in a hydrogen atmosphere for three days. To test that the diffus- ion of cobalt into iron had been successful, and to remove any surface activity, the specimen was immersed in HCl for one minute. The activity after this operation was seen to be slightly lower than i t had been previously, but after subsequent acid washes i t was constant. iv. The ^Co Experiment The objects of the experiments with ^Co were to see what temperatures could be achieved, how long they would last, and how much current was needed in the polarizing solenoid to saturate the magnetization of the iron. These experiments were'done before the guard shield was used and before the Ge(Li) detectors were obtained. The salt used was chrome potassium alum. Two Nal detectors were used, one in the axial and one in the equatorial positions. Their distances from the source were 9.0 and 3.5 58 59 inches respectively. With a radius for the detector face of 1 inch, the solid angle factors of equation 1.54 axial: are a = 0. I l l g 2 = g 4 = 0.970 equatorial: 0.991 a = 0.286 g 2 g4 = 0.940 = 0.809 . 2.3 Counting times of 200 sec. were used for the cold counts, while the warm counts were usually accumulated for 1000 sec. An example of a 200 sec. warm count is given in figure 10. The current in the polarizing solenoid was increased from 0 to 5 amps at 1 amp intervals, and similarly decreased. each current setting. of the ^Co Counts were taken at In figure 11 is shown a plot of the anisotropy peaks as the current was increased and decreased. Saturation occurs around 3.5 amps. For this polarizing solenoid this current corresponds to a central field of 1.5 kilogauss. The quantity plotted in figure 14 is the anisotropy as is usually defined in the literature: _ W(TT/2) W(0) WO/2) ' Z ' 4 More will be said about this function in section 2'.,7i. Figure 12 shows the plots of the temperatures of the ^Co elapsed time since demagnetization. against The temperatures were calculated by 60 6 0 CO SPECTRA, 2 0 0 SECONDS 60 AXIAL 40 20J 2 Z < o CC ai a. w tz EQUATORIAL O O 300 200 100 h 40 80 120 CHANNEL Figure 10 160 NUMBER Sample ^°Co spectra 200 240 Figure 11 Co anisotropy versus polarizing solenoid current 62 211- RUN 1 17 IM 13 5 V) z > UJ UJ 25 RUN QC 2 I— < CC UJ I 21 UJ < 17 S 13 9- I "II. I I I 030 TIME Figure 12 100 AFTER 1:30 DEMAGNETIZATION 2:00 (HOURS) Co temperature versus time elapsed from demagnetization 63 taking a weighted average of the temperatures indicated by the axial and equatorial detectors, thus 2 ° = T 2 2T L °1 a 2 = 2 + + G L 7 2 °2 2 2 2, o 2 // R[ a 1 + a 2 ] 2 l 2.5 where a. is the variance of T. and a 2 I i s the variance of T. The error bars in figure 12 indicate plus or minus one standard deviation. Figure 12 shows that, at this early stage, temperatures less than 20 mK could be maintained for over an hour and a half and that lower temperatures around 10 mK could last a half hour. In more recent experiments using a guard shield around a main p i l l of cerous magnesium nitrate a temperature below 12 mK was maintained for 7 hours. 2.7 The Term "Anisotropy" The usual definition of "anisotropy" is given by equation 2.4. As an average of the effect seen by the two detectors i t i s undesirable because i t is not linear in both detectors. Even the difference WO/2) - W(0) places equal weight on both detectors, which may not be justified. For example, consider the case in which the fourth order term in equation 1.8 is zero. We then have W(0) = 1 + A W(TT/2) = 1 - 2 I A 2.6 64 where = ^2^2^2* Clearly i - this case a far better measure of the n "effect" is the term k^. If both W(0) and W(TT/2) have the same error, then A^ estimated by the axial detector has half the error as A^ estimated by the equatorial. The axial result should therefore be weighted four times the equatorial figure. For these reasons very l i t t l e use will be made in this work of equation 2 . 4 . The word "anisotropy" will instead be used to refer to a non-uniform angular distribution. The "effect" seen by a detector at an angle 6 to the axis of orientation is e(0) = W(6) - 1. 2.7 The "effect" is the increase of the gamma ray intensity relative to i t s warm value. A negative effect is exhibited by ^Co in the axial direction. 65 CHAPTER THREE ORIENTATION OF F e NUCLEI IN IRON 59 3.1 Introduction 59 The Fe nucleus has been oriented by Tschanz and Sapp (1970) in the hyperfine fields of paramagnetic Ce-Zn nitrate and Nd-Zn nitrate. An unexpectedly large anisotropy was observed, the temperature dependence of 59 which yielded a magnetic moment for the Fe nucleus of 1.1 +_ 0.2 u ^ . However, there are several uncertainties in this work, especially regarding the fraction of iron atoms at the double nitrate sites. A verification of their result by a different method would be useful. 59 decided to orient We Fe in iron since in this system there exist no problems about solubility. The59hyperfine f i e l d is -333 kG (Hanna et a l . 1960). 3.2 Decay Scheme of Fe 59 The decay scheme of Fe i s illustrated in figure 13. The parent 59 nucleus beta-decays with a 45 day h a l f - l i f e to two excited statesLof Co which decays to the ground state emitting gamma rays of energies 1.292 and 1.095 MeV. The life-times of the excited states should be sufficiently short (compared to the nuclear spin precession time) to neglect intermediate state reorientation. Tschanz and Sapp (1970) consider the possibility of some reorientation but reach no definite conclusion. In any event, com- parison between our results and theirs can be made under the assumption that no reorientation is present. The main features of the decay scheme in figure 13 are obtained from Lederer, Hollander and Perlman (1968). The half-lives of the inter- 66 1.292 M E V 1.095 M E V 0.0 Figure 13 Fe decay scheme MEV ' 67 mediate states, measured by Sidhu and Gupta (1967), Nordhagen, et^ al_. (1967) and Agarwal et_ al_. (1967) are 0.59 nsec. and 4-50 psec. for the 1.29 and 1.10 MeV levels respectively. Following Blair and Armstrong (1965) the 59 spins of the both excited states of Co are known to be 3/2. Internal conversion coefficients measured by Collin et^ al. (1964) indicate that both gamma transitions are predominantly E2. For the beta-decay the mixing of the Fermi (L = 0) and Gamow-Teller (L. = 1) contributions i s 0 P P unknown. 59 An estimate of the magnetic moment of the obtained from the shell model. neutrons. Fe nucleus may be The nucleus contains 26 protons and 33 The number of neutrons i s 5 in excess of the 28 necessary to close the li^j^ shell; these 5 are to be distributed among the 2 p ^ an< 2 ^ l f ^ 2 shells such that the total spin is 3/2. Clearly the configurations (2p ^ ) (lf ^ ) and (2p ^ )°(lf ^ ) have total spin 5/2. In the config1 4 2 3 3 2 J urations (2p y ) ( / ) > ^3/2 - 5/2- 3/2^ - 5/2" 4 3 2 5 1 5 2 3 2 5 lf 3 2 2 ( 5 lf ) a n d ( 2p ) ( lf ) w e a s s u m e 2 that the neutrons tend to pair within each shell such that there i s no resultant contribution from shells with an even number of neutrons. It remains to find the magnetic moment of the shells with an odd number of neutrons. For a single neutron with orbital angular momentum L and total angular momentum I, the magnetic moment i s p = -(1.91)u (L = I - 2~) N Y7f Ci-9i)y The single neutron states 2 p ^ an< 2 N ^ ^5/2 (L = I + j) t n e r e : L : o : r e . 3.1 have magnetic moments 68 -1.91 and +1.36 respectively. Chapter 12 of Rose (1967) discusses the situation of three fermions in equivalent orbits (that i s , states differing only in the Iz quantum number). In order to achieve a totally anti-symmetric wave function certain limitations must be placed on the allowed values of the total spin J. The resultant magnetic moment i s (J/I) times the single particle 3 moment. For both (2p^ ) a n 2 ^ ^^5/2^ 3 a t o t a l s Pi- °f 3/2 is allowed and n the magnetic moments are -1.91 and (3/5) (1.36) = 0.82 u N respectively. Although no definite value has been predicted this discussion does 59 indicate that the magnetic moment of the Fe.nucleus is likely to be small. 3.3 Expected Anisotropy- of 59 Fe 59 Because the spin of the Fe nucleus is 3/2, the angular distribu- tion function of equation 1.8 need be carried only to the second order. We have W(0) = 1 + A 2 = W(TT/2) = 1 - j A2 A =B U F 2 The quantity A 2 2 2 2 3.2 i s the most convenient measure of the anisotropy. For a pure E2 transition between states of spins 3/2 and 7/2, the coefficient F 2 is F 2 = -0.14286 3.3 69 The re-orientation following the beta-decay yields Fermi) and the = 0.2 (pure Gamow-Teller). = 1.00 (pure For the present we shall assume situation which produces the least anisotropy, that is, pure Gamow- Teller. The orientation parameter has an extreme value of 1.0 for com- plete polarization of the nuclear spins. The ultimate value for k^ i s therefore -0.0286. Of course, i f there i s any Fermi contribution to the beta-decay, the anisotropy will be increased. Assuming a working value for the magnetic moment as 1.0 and using the hyperfine f i e l d of -333 kG, one calculates the temperature dependence for k^ as in figure 14. It can be seen from this graph that the effect i s small even at temperatures of 10-12 mK which we expect to achieve. 3.4 The Experiment 59 The Fe source was prepared by Atomic Energy of Canada Limited. An iron f o i l 1.0 cm. x 0.8 cm. x 0.03 cm. was neutron irradiated until 59 the activity of Fe was 20 uCi. 60 A specimen of Co dissolved in iron was prepared in the manner described in section 2..'6. The activity of ^Co was 2 yCi. 59 Because the expected anisotropy of the Fe gamma rays is so small, i t was decided to use a salt p i l l of cerous magnesium nitrate. Lower temperatures should therefore be possible, but warm-up times would be faster since this salt has a lower heat capacity. The guard p i l l was made of ferric ammonium alum. The two gamma ray peaks of *^Fe are too close to those of ^Co to 70 Figure 14 A 2 versus inverse temperature for by Fe permit the use of Nal detectors. employed. Instead, the two Ge(Li) detectors were The axial detector was 25 cm. from the source and i t s detecting area had an effective radius of 2.3 cm. The equatorial detector was 8.5 cm. away, and i t s detecting area had an effective radius of 2.0 cm. The solid angle corrections were: a = 0.092 axial: g equatorial: 2 = 0.994 a = 0.236 g 2 = 0.959 3.4 Sample spectra are shown in figure 15. The count times were set at 200 sec. so that the temperature could be closely monitored. After the f i r s t demagnetization the samples warmed from 15 mK to 60 mK in one hour, but on later runs they remained below 30 mK for an hour. Afterwards counts were taken for only 20 min. before another magnetization was begun. In this way more cold counts could be obtained and less time was spent at the warmer temperatures. 3.5 Analysis After eleven demagnetizations were conducted there were 85 pairs 59 of spectra to be analyzed. It seemed that any effect present in the Fe gamma rays was indeed very small, and that i t would be necessary to reduce the statistical error as much as possible by combining spectra of the same temperature. It was hoped that a measurable temperature dependence could be found, for such a phenomenon would allow one to f i t the results to the 72 FE 1.10 10005 EQUATORIAL 800 FE 1.29 600H CO 400 1.17 CO 1.33 UJ Z Z200| < X o OC UJ OV) Z300{ O o AXIAL 240 180 120 60 40 Figure 15 Sample ~80 120 ••••• .«•••*• 1 •• • 200 160 C H A N N E L NUMBER Fe/ Co spectra t nrirnri'-T-ii 240 known temperature dependence of the B 2 term. If this could be done, one would have 3 = p H^^/IkT independent of the other terms in k^. An observable effect would also have made i t possible to carry out an NMR/ nuclear orientation experiment. Further, i t would even be possible, in principle, to find those other terms and determine the beta-decay mixing ratio. For these reasons the temperatures of the spectra were carefully evaluated. The f i r s t step was to analyze the peak intensities by means of the program SPECTAN, described in section 2.(5. Table 1 gives the output of this program for one pair of warm and one pair of cold spectra. The second step was to calculate the temperature of each spectrum by means of the ^Co intensities. First, a grouping of spectra was carried out. Consecutive counts with l i t t l e difference in ^Co intensities were put together; four or five spectra would be in each group, and there were 21 groups in a l l . Then W(0) and W(TT/2) were evaluated for the ^Co peaks and a weighted average was taken of the axial and equatorial measurements of the temperature. For the 17 groups with a temperature below 25 mK Table 2 shows the axial and equatorial intensities and temperatures, as well as the average temperatures. With the aid of Table 2 groups which f a l l within the same temperature ranges could be combined. Only then, when the temperature of each 59 spectrum had been so ascertained were the Fe peaks investigated. The standard deviation errors of the group temperatures were used as a guide to selecting the temperature ranges. few standard deviations wide. Each range should be a The ranges chosen were: 13 to 14 mK, 14 to 15 mK, 15 to 17 mK, 17 to 19 mK, 19 to 23 hiK, and 23 to 27 mK. One may 59 now find A„ for the Fe peaks. The average values within these ranges, 74 TABLE 1. Sample output of the spectrum analysis program for Uncorrected Corrected Peak Count Axial 1 1450 1182 (Warm) 2 870 872 3 525 421 4 420 398 Equator. 1 7239 5162 (Warim)) 2 3899 3569 3 2432 1606 4 1614 1467 Axial 1 1436 .990 .027 (Cold) 2 861 .•998.9 3 342 4 Equator. (Cold) Spectrum Norm 59„ , 60„ Fe and Co. Norm Error Window 1237 1.046 .034 ( 45, 56) .034 813 .983 .037 (186,197) .6650 .036 239 .568 .050 (100,110) 253 .602 .038 240 .603 .051 (217,228) 1 7453 1.029 .012 5171 1.002 .018 ( 38, 55) 2 3947 1.012 .016 3626 1.016 .019 (174,197) 3 2789 1.147 .022 1961 1.220 .036 ( 91,109) 4 1874 1.161 .027 1718 1.170 .032 (205,226) Peak identification: Error 1 2 Count 59 F e , 1 .10 MeV 59 F e , 1 .29 MeV 3 60 , Co, 1 .17 MeV 4 Co, 1 .33 MeV rn 75 TABLE 2. Cobalt temperatures from ^Co intensities during ^^Fe runs, Group W(-0) W(TT/2) T(0) mK T(TT/2) mK T (Average) mK 1 .696 +_ .017 1.122 + .013 15.2 +_ 0.6 16.6 +_ 1.5 15.4 + 0.6 2 .815 .018 1.105 .013 21.5 1.2 18.8 2.0 20.8 1.0 3 .880 .018 1.017 .012 28.5 2.7 24.9 2.8 26.8 1.9 4 .667 .015 1.168 .012 14.2 0.5 11.7 1.0 13.7 0.4 5 .666 .017 1.116 .013 14.2 0.5 17.4 1.6 14.5 0.4 6 .772 .018 1.094 .012 18.7 1.0 20.6 1.7 19.2 0.9 7 .758 .018 1.121 .012 18.0 0.8 16.7 1.5 17.7 0.7 8 .797 .018 1.092 .012 20.2 1.1 34.2 1.7 24.3 0.9 9 .726 .018 1.132 .012 16.4 0.8 16.0 1.2 16.3 0.7 10 .798 .018 1.101 .012 20.2 1.1 19.5 1.9 20.0 1.0 11 .728 .016 1.114 .011 16.5 0.8 17.6 1.5 16.7 0.7 12 .727 .015 1.141 .011 16.5 0.8 14.5 1.6 16.1 0.7 13 .635 .015 1.159 .011 13.1 0.5 12.7 1.0 13.0 0.4 14 .694 .016 1.158 .011 15.1 0.6 12.7 1.0 14.5 0.5 15 .8151 -016 1.076 .011 21.5 1.4 23.7 2.3 22.1 1.2 16 .649 ..017 1.156 .012 13.6 0.6 13.0 1.4 13.5 0.5 17 .668 .015 1.138 .011 14.2 0.5 14.8 1.0 14.3 0.4 76 as measured by the axial and equatorial detectors are given in Table 3. TABLE 3. Average A„ for Fe for each temperature range. 1.10 MeV Peak Temperature Range (mK) A 2 (Axial) 13-14 + .008 + .008 14-15 + .005 .008 15-17 + .005 17-19 A^ (Equatorial) -.0006 A 2 (Average) +_ .009 + .0037 -.018 .009 - .0067 .006 .008 -.020 .008 - .0085 .005 -.009 .012 -.032 .016 .020 .010 19-23 -.013 .006 +.0015 .008 - .0056 .005 23-27 + .011 .008 -.010 .011 + .0005 .00068 +_ .010 + .0043 + .0067 +_ .006 1.29 MeV Peak 13-14 +.0006 +_ .009 +0.008 14-15 +_.005 .009 + .018 .0096 + .0115 .0066 15-17 + 0006 .0075 -.011 .008 - .0025 .005 17-19 -.004 .013 + .010 .016 + .003 .010 19-23 -.0008 .00068 -.010 .008 + .0046 .005 23-27 -.003 .009 -.002 .011 .0025 .007 77 3.6 Results and Discussion For the size of the effects r e l a t i v e to the s t a t i s t i c a l errors no f i t t i n g t o the temperature dependence can be done. Instead, one can average A 2 over the lowest three temperature ranges (13 - 17 mK) with the result: A, (1.10 MeV) = -0.0038 + 0.0033 A2 (1.29 MeV) = +0.0044 +_ 0.0036 3.5 Because of the beta-decay mixing there i s no reason to expect the two gamma rays to possess exactly the same value of k^, however, they both must be negative. Also the minimum effect f o r pure Gamow-Teller decay i s the same for both peaks. over the two gamma rays. We are therefore j u s t i f i e d i n averaging The conclusion i s that to within an error of 0.003, A 2 i s zero. 59 An upper l i m i t may be placed upon the magnetic moment o f the nucleus. Fe By assuming that the beta-decay i s pure Gamow-Teller, we have at the mean temperature o f 15 mK: |B2 U 2 F 2 | '< 0.003 B 2 < 0.105 yH IkT y < 0.480 < 0.9 y. N 3.6 Although, t h i s r e s u l t i s just consistent with the r e s u l t of Tschanz and Sapp ( y = 1.1 +_ 0.2 y ^ ) i t does not explain the large effect which 78 they observed. At a temperature of 10 mK they find A 2 ~ -0.05. The hyper- fine field in their system is 594 kG, thus 10 mK corresponds to 6 mK in 59 our work. Furthermore, they estimate that only 30% of the Fe nuclei were at double nitrate sites and that the remaining nuclei were in pockets of frozen solution contributing nothing to the zero field anisotropy. 59 Their value of k^ z -0.05 is an average over a l l value of A 2 Fe nuclei, hence the for those nuclei in the sites is much larger s t i l l . The largest possible anisotropy occurs i f the beta-transition i s pure Fermi, in which case A 2 = -0.14 when nuclear spins are fully saturated. 59 In that event the upper limit imposed on the magnetic moment of Fe by our work is 0.38 u ^ , considerably less than that measured by Tschanz and Sapp. A 2 Since they determined the moment by the temperature dependence of their result is independent of the beta-decay mixing. Tschanz and Sapp measured the temperature from the susceptibility of the crystals. Since they used Ce-Zn nitrate for low temperatures (< 10 mK) and Nd-Zn nitrate for higher temperatures any f i t of their results to a theoretical curve requires accurate matching of the two temperature, scales. This is especially so because the overlap region contains the greatest variation in k^. Any mis-matching can greatly alter the value of the magnetic moment. Moreover the lack of agreement between the theoretical f i t s and their data implies that the temperature dependence is not contained solely in B 2> The most questionable assumption in the interpretation of their results concerns the effect from nuclei not in the double nitrate sites. If this effect had a separate temperature dependence, instead of the constant effect assumed, the observed anisotropy would not show the B depen2 dence . 79 These d i f f i c u l t i e s , the temperature uncertainty, the influence of the nuclei outside the sites and the non-FS^ behaviour of k^, make their experiment hard to interpret. Our measurement, however, i s easy to interpret but contains a null result. Lower temperatures might yield an observable effect, which would permit an NMR/nuclear orientation experiment. This would best resolve the problem of the magnetic moment and would also settle the other factors in A . 0 80 CHAPTER FOUR ORIENTATION OF RUTHENIUM NUCLEI IN IRON 103 4.1 Introduction 103 Several new examinations of the decay scheme of Ru (Avignone and Frey, 1971, Pettersson et_ al_. 1970) have settled one of the spin assignments, making i t possible to use the method of nuclear orientation to resolve the multipole mixing of the dominant gamma transition. In- ternal conversion studies have been unable to fix the mixing ratio because the two multipolarities involved have almost the same theoretical conversion coefficients. The nuclear orientation factor F, on the other hand, could be very sensitive to the mixing. According to Hansen (1958) Ru is soluble in iron to a few atomic percent, which is quite sufficient for the trace impurity needed to 103 orient Ru in iron. A small anisotropy [e = (1. +. 0.5)%] for Ru in iron at a temperature of 40 mK has been reported by Kul'kov et al_. (1965), and thus one would expect a much larger effect at 10 mK. The hyperfine field on Ru in iron is known from Kistner (1966). Thus anisotropy of a few percent would allow the NMR/nuclear orientation technique to be used to measure the magnetic moment with great precision. The decay scheme parameters could be found from the remaining factors in the effect. 81 4.2 Decay Scheme of 103 Ru Figure 16 shows part of the decay scheme for 103 Ru following Avignone and Frey (1971), Pettersson, et_ al_. (1970) and Raeside, et a l . (1969). The dominant decay mode i s by beta emission to the 537 keV 103 level of Rh followed by a gamma transition (of energy 497 keV) to the 40 keV level. The 497 keV gamma i s known from internal conversion data to be a mixture of E2 and Ml multipolarities. Avignone and Frey 2 report that Ml i s favoured, while Pettersson et_ al_. give 6 (E2/M1) = (1-17)%. Previously the spin of the 597 keV level had been known to be either 7/2 or 5/2, but both Avignone and Frey, and Pettersson e£ al_. have found an El transition between the 537 keV level and the 295 keV level of spin 3/2, which rules out the possibility of a spin 7/2 for the higher i i level. 3 The second most dominant gamma transition i s the 610 keV one between the 650 and 40 keV levels. The intensity of this gamma ray i s 6.2% of the intensity of the 497 keV gamma ray. Its multipolarity i s Ml or E2; no mixing ratios have been determined. 103 The Ru nucleus contains 44 protons and 59 neutrons. This number of neutrons i s 9 over the magic number 50. o r in a 2d^2 -^§7/2 sne H- The odd neutron is The Schmidt limits for the 2d^^ magnetic moment are -1.91 and +1.36 y^. For three neutrons in the l g 7 / 2 or bit, such that total spin i s 5/2, the magnetic moment i s (5/7) times the moment for one neutron in that shell, as discussed in section 3.2. The •z Pettersson et^ al_. (1970) i n i t i a l l y concluded that the 242 keV transition is M2 because of a plotting error. In an erratum they explain that the El assignment i s correct. 82 Figure 16 Ru decay scheme 83 3 resultant moment for the ( l g ^ 2 ) configuration is +1.06 u^. However, these numbers can only be considered as guides because the large number of nucleons makes a proper calculation very complicated. One can only conclude that the magnetic moment of the 103 Ru nucleus would be expected to be between -1.91 and +1.36 jl,.. N 4.3 Expected Anisotropy of the 497 keV Gamma Ray 103 The spins involved in the decay of Ru through the 497 keV transition are such that the second and fourth order terms of equation 1.8 will be present. For a parent spin of 5/2 the orientation parameters B2 and have maximum values of 1.3363 and 0.4629 respectively. The value of 3 which makes B2 equal to half i t s maximum is 1.0. 99 The hyperfine f i e l d on Ru in iron has been found by Kistner (1966) to be 500 kG. If one assumes a working value of one nuclear 103 magneton for the magnetic moment of Ru, then B = 1.0 corresponds to T = 7.3 mK. The other unknown factors in the anisotropy calculation are the Fermi/Gamow-Teller mixing ratio and the M1/E2 mixing ratio; table 4 gives the possible U and F coefficients. Since F^ = 0 for an Ml transition (which we expect may dominate) and since ^/B^ > 10 for 8 < 1.0, we may, at least temporarily, neglect the fourth order term of equation 1.8. The sign of the anisotropy effect is given by the sign of the F 2 term, which can be negative by virtue of the large interference term. For mixed multipoles the resultant F^ i s calculated by equation 1.13 and plotted as a function of the mixing ratios in figure 17. The range of F 2 for the 497 keV peak is +0.93 to -0.47 84 TABLE 4. U and F coefficients for possible transition modes in the 497 keV and 610 keV sequences in Transition Mode 103 Ru 497 keV Fermi 1.00 Gamow-Teller 0.65714 610 keV 1.00 -0.14286 0.87482 0.58029 Transition Mode Ml 0.13363 0 E2 0.32453 0.11756 M1/E2 Interference 0.69437 0 -0.43643 0.24939 -0.37796 0 -0.47755 0 Figure 17 F ? versus mixing ratio for 497 and 610 keV gamma rays 86 and i t s zeroes occur for 6(E2/M1) = -0.1 and 6(M1/E2) = -0.24. The 2 former value is just within the range 6 (E2/M1) = ( 1-17)% given by Pettersson et_ al_. (1970) . The anisotropy of the 497 keV gamma ray can be large or small, positive or negative. Of particular interest are the sign of the effect and i t s temperature dependence. There is a similar lack of information about the 610 keV gamma ray. It is known to be Ml or E2, or any mixture of these. The U and F coefficients for the transitions in this sequence are given in table 4. The anisotropy of this radiation may also be large or small, positive or negative. The coefficient F^ ranges from -0.604 to +0.416 with zeroes at 6(M1/E2) =0.28 and 6(E2/M1) = -0.50, as shown in figure 17. 4.4 Experiments i. Source Preparation The sources for the experiment were made by the same method used for the ^Co sources. To prepare a sample some ^^RuCl^ obtained from AmershamyiSeariee, Corp. of Arlington Heights, Illinois, was evaporated on an iron f o i l of dimensions 1.0 cm. x 0.8 cm. x 0.03 cm. Since Ru is more electronegative than Fe, the Ru atoms replace the Fe atoms i n the foil. The sample was furnaced for 72 hours at 900°C in a hydrogen environment to clean the surface and diffuse the Ru atoms into the iron. Afterwards, surface activity was removed by etching the specimen with 103 HCl. Because of the short h a l f - l i f e of Ru several samples had to be made during the course of the experiments. The activities of the specimens ranged from 5 to 10 uCi. 87 ii. Initial Runs The *^Ru source was soldered with a separate ^Co source to the end of the copper heat link. The main paramagnetic salt used was cerous magnesium nitrate with a guard p i l l of manganous ammonium alum. The guard shield was around the main p i l l and source. The counting in the f i r s t experiments was done by Nal and Ge(Li) detectors. These early runs showed that no large effect existed in the 497 keV peak at 10-15 mK, and that the Compton background from ^Co might mask any.small effect. The Nal detector was especially suspect of faulty background separation because the continuous spectrum on either side of the wide peak was very far from being linear. Fitting to poly- nomials was tried, but was finally abandoned in favour of a more direct way of determining the background, as discussed below. iii. Final Runs As i t was important to estimate accurately the amount of background originating from the ^Co spectrum, a separate run was undertaken using ^Co alone. This was necessary because the precise shape of the ^Co spectrum depended not only on the geometry and electronic settings but also on the temperature of the sample. The two Ge(Li) detectors were used, with the more efficient one in the more distant axial position. The run, with a 5 yCi source of ^Co in iron, covered the range of 8.5 mK to 15 mK with 18 pairs of spectra of count time 1000 seconds each. The next run was conducted under identical conditions, but with 103 the inclusion of an 8 uCi.source of Ru. The background in the 497 keV 88 window was 45% of the total count. A total of 73 pairs of spectra was obtained, at temperatures from 9 to 13 mK. The results of this run did show a very slight effect. However, in order to remove any doubts about whether the 60 Co might be causing this, a second run followed using the same geometry and apparatus, but 103 with a much larger Ru/Co ratio. The Ru source was 28 uCi and the ^Co source was 0.8 yCi. This ruthenium activity was considerably larger than would normally be considered desirable in view of heating effects i n the sample but proved to be tolerable for the purpose.of these tests. The background in the 497 keV window was 1% of the total. In this run 35 pairs of spectra were obtained in the temperature region of 10-12 mK. Figure 18 shows two axial spectra, one from the preliminary ^Co run and one from the second run with the very large Ru/Co ratio. The count time for each i s 1000 seconds. 4.5 Analysis The background run was analyzed f i r s t . The background in the 497 keV window was found to have a temperature dependence similar to, but less than, that of the ^Co peaks. In the ^ Ru/^Co runs, the number of counts in the 497 keV 3 window was f i r s t determined knowing the temperature from the analysis of the ^Co peaks, the corresponding background was removed and the decay correction applied. Statistical errors were also calculated. 89 8r 6 0 CO 1.173 BACKGROUND 1.332 co O _L z" z < 0 80 40 120 CHANNEL X 160 200 240 NUMBER u cc UJ 0. 1 0 3 RU AND 6 0 CO 0.497 MEV 1.332 t -:30h ' z' 1.173 o X 30 20\ 10| 0.610 40 Figure 18 80 Sample spectra from 120 160 Ru runs 200 240 90 i. Background Run The analysis of the ^Co background run determined the ratio of the continuum background within the window of the 497 keV ruthenium transition to the intensity of the ^Co peaks. This background consisted of true Compton counts from the Co peaks plus radiation scattered 60 from other directions. Ambient background would also be in this window. The ratio of this background within the 497 keV window to the ^Co peak intensity, denoted.as R, depends on the selected windows and i s quite different for the two detectors. . It i s the temperature dependence of R that i s desired, but since W(0) and W(TT/2) for the ^Co peaks are functions of temperature and since they are much simpler to measure, R was found as a function of the corresponding W. The f i r s t step in the analysis of the background spectra was to find the temperature dependence of the 497 keV background and to f i t i t as follows. Consider a gamma ray emitted in the direction 6, which i s scattered into the detector at a different direction, thus recording a spurious event, as i n figure 19. There i s a probability f(9) that Figure 19. Gamma Ray Scattered into Detector. 91 such a scattering occurs and that the scattered radiation f a l l s within the 497 keV window. Since the intensity of gamma radiation into the direction 9 is W(9), the number of background counts in the 497 keV window i s BG = A + J'f(9) W(9) d(cos9) 4.1 where A is a constant background. By virtue of equation 1.8, W(9i),egandbe writ.teniassaosumeof -.. Legendre polynomials of second and fourth order. The scattering func- tion f ( 9 ) may also be expanded into Legendre polynomials, but as a result of the integration in equation 4.1, only the zero, second and fourth orders survive. The temperature dependence of the background is contained i n B and B^, the ^Co orientation parameters. 2 BGfJD = BG [1 + aB (T) + bB (T>] Q The observed background/temperature 2 4 Thus 4.2 points can be fitted to equation 4.2 to find a and b for each detector. The constant BG is found from o the warm counts. Table 5 l i s t s the results of such a f i t for the back- ground run as well as the corresponding coefficients for the ^Co peaks. The ratio R could now be calculated. For every temperature, the coefficients in Table 5 could be used to find the 497 keV background and intensity of the ^Co peaks. In figure 20 are plotted the fitted R-W curves for each detector; also shown are the experimental points from which the f i t s were made. 92 93 TABLE 5. Coefficients of B„ and B. from Co background run Window a b Axial 497 keV Background -0..2657 -0,.1358 Axial Co* -0,.41702 -0..23606 0,.1600 -0,.0706 0,.20356 -0..08173 60 Equatorial 497 keV Background i3 * • i Equatorial 60 Co * * Includes solid angle correction. ii. Ru Runs 103 The h a l f - l i f e of Ru is 40 days, a time short enough that during the 48 hours of the experiment the source decays by 2.5%. Since effects of less than 1% were being sought, this decay had to be carefully corrected. The h a l f - l i f e given by Flynn, Glendenin and Steinberg (1965) is 39.6 +_ 0.1 days. The program SPECTAN summed the channels in the Co and 497 keV 60 60 windows, and calculated the decay factor. A second program used the peaks to determine the temperature, the corresponding background/^Co ratio, and thus the background in the 497 keV window. This was removed and the decay correction applied. in the 497 keV gamma ray peak. The result was the number of counts Co 94 iii. Statistical Error The statistical error in the background arises from the s t a t i s t i cal error in the ^Co peaks. Let nab'e the number of counts in these peaks, N the number at high temperature, R the ratio of background and ^Co intensities. The ratio R i s a function of W = n/N which i s approxi- mately linear over changes i n W of the order of /n/N. Therefore we may write R ~ a + b W. The background becomes BG = Rn 2 = an + bn /N = £ (n. + cn) 4.3 2 where c = aN/b. The variance in the background i s then VCBG) = (b/N) V(n + cn) . 2 4.4 2 However, the variance on the right hand side of equation 4.4 can be expanded thus: V(n 2 2 2 2 2 + cn) = <(ri + cn) > - <n + cn> 4 3 2 2 22 = <n > + 2c<n > + c <n > - <n > 2 2 2 - 2c<n ><n> - c <n> For a Poisson distribution of mean A . 4.5 95 A <n> 3 2 A + 3A + A A 4 + 3 2Z 6A° + 7A + A. 4.6 The mean A may be approximated by n; since this is a large number 3 (MO ) only the highest order need be retained. the order of N). (Recall that c is of Then V(BG) 2 3 2 2 ~ (b/N) [4n + 4cn + c n] = 4 b W Rn + na 4.7 The variance in the background is estimated by equation 4.7. 4.6 Results The corrected normalization counts should differ only by statistics. A test of the detector stability throughout the run could be provided by comparing the sample variance of the 497 keV warm counts to the variance expected from equation 4.7. These quantities agreed to within 1.5 standard deviations. l. The 497 keV Transition It was convenient to divide the cold counts from the f i r s t run into two temperature regions: 9 - 10 mK and 10 - 13 mK. run, which because of the stronger 103 For the second Ru activity did not achieve for long 96 a temperature below 10 mK, the counts in the one temperature region, 10 - 12 mK, were selected. The average effect of the 497 keV gamma ray in the detectors is shown in Table 6. TABLE 6. Effect in the 497 keV gamma ray Temperature (mK) 1 - W(0) 9-10 1st run 0.0101 +_ 0.0016 10 - 13 1st run 0.0059 +_ 0.0018 10-12 2nd run 0.0048 +_ 0.0007 WO/2) - 1 0.0025 + 0.0005 The results in Table 6 do not include measurements from the f i r s t run in the equatorial direction. The equatorial detector was found to have a small attenuation (>1%) when the large cooling magnetic field was on. No such effect was seen on the axial detector. Since in the f i r s t run a l l normalization counts were taken with this field on i t s results are suspect. As a further test for errors in the background correction a comparisom.'was made between the normalized count of the 497 keV peak corrected 97 by the R-W method and the same count corrected by the linear background method. For the f i r s t run there was agreement to within 0.008 with neither method producing consistently higher numbers. The errors for individual spectra were 0.011 for the R-W method and 0.015 for the linear method. For the second run, agreement was within 0.001 with individual errors of 0.004 for the R-W method and 0.005 for the linear method. The confirmation of the background correction by a different method of calculation indicates that both should be reliable. We there- fore conclude from the axial detector that the 497 keV gamma ray exhibit an effect of -0.004 to -0.005 in the axial direction at a temperature of 11 mK. The f i r s t run also indicated an even greater effect in the 9 - 1 mK region. However, because this effect was not confirmed in the second run and because the f i r s t and second runs did not agree entirely in the 10 - 13 mK region, the lower temperature result was not accepted as proven. Only the result of the second run with i t s large Ru/Co ratio was considered conclusive. ii. The 610 keV Gamma Ray The smaller intensity of the 610 keV gamma ray makes the statistical error in i t s analysis much larger. Also, the ratio of background to peak i s larger, creating more d i f f i c u l t i e s in the background removal. The results of analysis from the f i r s t run, where the background i s 90% of the total, are too uncertain to be reliable, but the second run, where the background i s only 18% of the total, does yield some useful data. 98 After correcting the counts in the 610 keV window using the R-W relations for this energy the 610 keV gamma ray shows the normalized intensities in the 10 - 12 mK region to be W(O) = 0.996 + 0.003 WO/2) = 0.999 + 0.002 . 4.8 With the possibility of systematic errors in equation 4.8 we conclude thaitethe : effectsiselesssthan 0.005;;" 1 iii. Magnetic Moment and Mixing Ratio With the negative effect seen in the 497 keV gamma ray in the axial direction one concludes that F i s negative, and hence that the 2 mixing ratio 6 must be negative, as in figure 17. 103 A lower limit on the magnetic moment of the Ru nucleus can be found by assuming the beta-decay and gamma-decay mixings are such that 2 ^2^2 ^ Sa m a x i m u I I K Using the extreme value 6 (E2/M1) = 0.17 given by Pettersson et a l . (1970), or 6 = -0.41, and the extreme case of pure Fermi beta-decay, |U F | < 0.3263 2 2 . 4.9 In order to observe W(0) = 0.996 at 11 mK, the value of B i s 2 2 0.004 0.3263 = 0.012. 4.10 99 For spin 5/2, has this value when y H ^ or eff > 0.098 y > 0.15 y 4.11 N where equation 4.11 uses the hyperfine f i e l d value of 500 kilogauss and a temperature of 11 mK. For the 610 keV gamma ray, the observation of an effect less than 0.005 when B = 0.012 indicates 2 B 2 U F | < 0.005 2 2 |F | < 0.416. 2 4.12 The E2/M1 mixing ratio for the 610 keV gamma ray is limited as follows: 6(E2/M1) > 1.15 or S(E2/M1) < -0.03. iv. 4.13 Summary 103 Although i t had looked possible that nuclear orientation of Ru in iron might resolve uncertainties in i t s decay and i t s magnetic moment, the small effect observed does not allow this to be done, other than to say that for the 497 keV gamma ray the E2/M1 mixing ratio i s negative and that the magnetic moment is greater than 0.15 nuclear magnetons. Some limitations may also be put onto the mixing of the 610 keV gamma ray. 100 CHAPTER FIVE SPIN-FLOP TRANSITION OF Mn C1 -4H20 OBSERVED 2 BY NUCLEAR ORIENTATION 5.1 Introduction 54 That the MnCl 2 Mn nucleus in a crystal of antiferromagnetic : 4H 0 can be oriented at low temperatures was demonstrated 2 by Daniels et al. (1961). The hyperfine field direction is defined by that of the electronic spins so that at low temperature the nuclear spins point along the axis of electronic magnetization. Turrell, Johnston and Stone (1972) have used the angular dependence 54 of the gamma radiation from Mn in CoCl • ^H^O to observe the 2 change in direction of the electronic spins as the crystal goes through the spin-flop transition. It should be possible to study the spin-flop transition in MnCl * 4^0 by the same method. 2 This experiment was carried out in collaboration with R.L.A. Gorling. 5.2 Molecular Field Model The dynamics of the electronic spins of an antiferromagnet can be described by the molecular field approximation of the Ne"el (1932, 1936), Bitter (1937) and van Vleck (1941). is given by Nagamija et a l . (1955). A thorough treatment of the model 101 The forces acting on the electronic spins are the exchange force, the anisotropy energy, and the external f i e l d . energy between spins ^ and The exchange is E. . = - h J. . Q + 4 IJ i j l v • £.) y 5.1 A simple model is made by letting J „ = J when spins i and j are nearest neighbours and J\ ^ = 0 otherwise. The spins may be considered to be members of two intersecting lattices A and B such that a l l the nearest neighbours of a spin in lattice A are in lattice B, and vice versa. For antiferromagnetism, J is negative. Equation 5.1 is minimized when the spins i n lattice A are antiparallel to those in lattice B. The anisotropy energy arises from the crystalline fields, and is often of the form E. = - D[S - \ S(S + 1)1 l IZ 3 2 J J L 5.2 For positive values of D the spins tend to align along the z axis, the easy axis of magnetization. The applied magnetic field adds the Zeeman energy E. = l * gB l & where B is the Bohr magneton and g the electronic g-factor. 5.3 102 Equations 5.1, 5.2 and 5.3 can be combined and summed over a l l the spins. Letting N be the number of spins in each lattice and q be the number of nearest neighbours, the total Hamiltonian, after dropping constant terms, i s X = ~£ J. -. s\ • = - 2 Nq J \ 3. -D IS - gg H • Z 2 • 3 B - ND[S 2 z + t. S^] - gBM • f§ + A . 5.4 In equation 5.4 a l l spins in lattices A are and a l l spins in lattice B are \ . Equation 5.4 can be re-written by introducing the molecular field H E = -2JqS/g6, the exchange f i e l d . The exchange f i e l d H~ = ^ E H A is the f i e l d produced by the spins in lattice A acting on the spins in lattice B . The anisotropy f i e l d , field. = 2DS/gB, is a fictitious Equation 5.4 i s written as X = ^ A * ^B -I C 4 + s2z ) --15 • c£ + A 5.5 Let H be along the z axis and consider two possible spin configurations. In the f i r s t case, the spins are antiparallel and the energy i s E(l) = - NgBS[H + H ] E A 5.6 E / S 103 In the second case, the spins are at an angle 6 to the z axis but have opposite projections in the xy plane. S • tn A B =S = S In this situation cos 2 0 2 (2 cos 0 - 1) 2 5.7a 2 S. + S = 2S cos 0. Az Bz 5.7b n The energy of this configuration i s E(2) = NgB[H (2 cos 0 - 1) - H 2 E cos © - 2H cos0]. 5.8 2 We choose the angle which minimizes E(2). This occurs when C 0 S 8 = 2f^H" 5 ' 9 and the minimum value of E(2) is E(2) = - NgBS[H + H /(2H - H^)]. 2 E E 5.10 The second case corresponds to the spin-flop configuration which occurs when the field applied along the easy axis exceeds a c r i t i c a l value H c = [H C2H - HA )] A h E When H > H , the energy E(2) < E(l). 5.11 104 The c r i t i c a l field in MnCl^ ' 4H2O has been measured with susceptibility and nuclear resonance studies by Gijsman et a l . (1959) to be 7.87 kG at 1.04 K. 7.50 kG. A linear extrapolation to zero temperature gives Metcalfe (1971) observes a spike in the thermal resistivity at 0.132 K when the field along the crystalline c-axis (the easy axis of magnetization) is about 7.0 kilogauss. Miedema et a l . (1965) using heat capacity measurements found their data consistent with the exchange field Hg = 12.8 kG and the anisotropy field H A = 2.2 kG. The advantage of the nuclear orientation experiment is that the spin-flop i s directly observed (assuming the nueleif-followttheeelectron spins). The interpretation in terms of H^ and H^ is made through equations- 5.9 and 5.11. 5.3 Decay Scheme and Anisotropy 54 The gamma radiation from oriented Mn nuclei is a suitable probe because i t s decay scheme is well established. The decay scheme:in figure 21 is from Lederer, Hollander and Perlman. The magnetic moment is 3.302y^ (Templeton and Shirley, 1967) and the hyperfine field on Mn in MnCl^ * 4^0 is deduced to be 645 ± 12 kG from the thermal capacity data of Miedema et a l . (1965). 54 The coefficients in equation 1.8 for U 2 = 0.82808 U F 2 = -0.59761 F 4 4 Mn are = 0.41785 = -1.06904 5.12 105 303 DAYS Figure 21 Mn decay scheme 106 The detectors used in the experiment subtended a half-angle at the source of 10°. The solid angle corrections are g = 0.974 and g^ = 2 54 0.916. The angular distribution of the gamma ray from Mn, with this solid angle correction, i s W(6) = 1 - (0.482) B P - (0.409) B P 2 2 4 . 4 5.13 Equation 5.13 is plotted for several temperatures as a function of cos0 in figure 22. As the electronic spins change direction, the nuclear spins follow by means of the hyperfine interaction. The direction of the spins can be found by locating the orientation axis relative to the axes of three gamma ray detectors. For each detector equation 5.13 finds the angle between the orientation axis and the line joining the detector to the radioactive source. It should be noted that i n the model in section 5.2 there are two spin lattices each with i t s own orientation axis. In the non-flop situation these axes coincide but are of opposite sense; the gamma radiation, depending on only the square of the cosine, cannot distinguish the two senses. In the flop phase, i t is not immediately obvious that the same symmetry exists. co-ordinates (0,<J)<) and The two orientation axes have spherical polar (0,ir +_<(>). Consider a detector along the axis with polar co-ordinates (6 , < | > ). The angle a between this axis and the axis '(9,40 i s d> ) cosa = cos0o cos6 + sin0 o sin0 cos(d>-o V T Y 5.14 108 For the detector along the z axis (the axial detector in the normal phase), 6 Q = 0 and a = 6 for both orientation axes. For a detector in the xy plane 8 Q = TT/2 and cosa = sine cos (<j>-<|> ) . The two orientation axes therefore differ only in the sign of cosa, and thus the gamma ray intensities from them are identical. For detectors in other positions this would not be so. 5.4 The Experiment i. Source The experiment was conducted upon a crystal of MnC^ • 4^0 which 54 had been grown from a solution containing some MnC^. The activity of the crystal was 15 yCi. It was covered in Apiezon N grease and fixed to the copper f o i l at the end of the heat link. Also soldered to the heat link was a ^Co-Fe f o i l of activity 2 yCi. The paramagnetic salt used was chrome potassium alum. ii. Detectors Three gamma ray detectors were employed. In the axial direction a 5 in. x 4 i n . Nal(Tl) detector (denoted C) was positioned 14 i n . from the source; the two Ge(Li) detectors were set up in the equatorial plane, both 4 i n . from the source. The crystals had been aligned so that the crystallographic c-axis (the .easy axis."of magnetization) was vertical; the axes of the Ge(Li) detectors were mutually perpendicular, with one (denoted detector B) 35° from the crystallographic b-axis, and the other (denoted detector A) at 125°. A l l angles quoted here have an error of 5°. 109 The analysis of the orientation axis is simplified by using three mutually perpendicular detection axes, but because each detector cannot distinguish the sign of cos0 an ambiguity exists in the location of the orientation axis. This ambiguity may be resolved by moving the equatorial detectors so that the angle between them is 45° (or 135°). In a subse- quent experiment the A detector was moved to make an angle of 125° with the B detector. The ambiguity was settled conclusively by the results of this run. Since the multichannel analyzer is able to route the signals from only two detectors at one time i t was necessary to alternate the second input between the A and B detectors. the f i r s t input. The C detector was always Sample B and C spectra are shown in figure 23. iii. Procedure Several normalization counts of 1000 sec. were taken before and after the experiment. After the demagnetization the current in the polarizing solenoid was increased to provide an external field along the c-axis of the MnC^^I-^O crystal?.; acquired at each current setting. Several counts of 400 sec. were The field in the centre of the solenoid had been calculated and measured to be 390 gauss per amp. The current was measured by the voltage across a resistor in the power supply for the superconducting magnets. 110 EQUATORIAL 0.835 MEV 80. 1.173 1.332 60] X10 40 "o 201 Z < I o cc 50 _L 150 100 200 250 C H A N N E L NUMBER UJ Q. v> 4 0 p Z AXIAL 3 0.835 MEV o o 1.173 X10 30h 1.332 20 10h 100 50 figure 23 Sample 5 4 Mn/^°Co spectra 150 200 250 Ill 5.5 Analysis i. Temperatures The program SPECTAN was used i n the normal way to find the r e l a t i v e i n t e n s i t i e s of the ^ C o 4 and ^ Mn peaks. copper heat l i n k , as given by the ^ C o The temperature of the thermometer, was .15 mK, but the MnC^^rL^O c r y s t a l was at a temperature of 65 to 55 mK, as measured by 54 the Mn anisotropy below the spin-flop t r a n s i t i o n . This higher temper- ature was due to the poor thermal contact and the high heat capacity of the crystal.?. The c r y s t a l continued to cool throughout the several hours of the experiment. ii. Spin-Flop Transition The spin-flop t r a n s i t i o n was analyzed by finding the orientation axis r e l a t i v e to each detector as a function of f i e l d along the c-axis. These angles were determined by locating the angle 6 corresponding to the observed r e l a t i v e i n t e n s i t y W(0) on the appropriate temperature curve i n figure 22. In the antiferromagnetic phase the c r y s t a l tempera- ture could be found from the observed anisotropy. Because the c r y s t a l continued to cool and because the temperature could not be measured d i r e c t l y i n the flop phase, the temperature was interpolated from those measured i n the antiferromagnetic phase before and a f t e r the c r y s t a l was i n the flop phase. The spin-flop effect i s shown i n figure 24 where cosa f o r each detector i s plotted against f i e l d applied along the c-axis. t r a n s i t i o n the cosine i s an average over a l l the spins. t r a n s i t i o n i s seen from figure 24 to be at H Through the The centre of the = 6.50 + 0.05 kG. 112 •8r DETECTOR A .4* O .21 < Z Q Z o CL CO UJ OC OC o CJ O JL .81 DETECTOR B • * .6 .4 UJ z 00 O O 1.0r DETECTOR C •8h .6 .4 7 FIELD ALONG C AXIS Figure 24 8 IN KG Cosa versus field along c-axis J 9 iii. Orientation Axis The orientation axis of the spins can be found from the angles i t makes with the three detector axes. The three angles are not indepen- 2 2 2 dent and can be fitted to the relation cos a. + cos a„ + cos a„ = 1 to A rs L yield the fitted values of these angles in Table 7. From these spherical polar angles for the orientation axis may be found as follows: 6 = a c tan<j> = cosa^/coscig 5.15 where <j> is the azimuthal angle measured from the B detection axis. TABLE 7. Measured and fitted angles of the spin orientation above the spin-flop transition Field Applied • along c-axis (kG) Measured Angles ct.. AA a_ „ B B (degrees) Fitted Angles a a. A B (degrees) c a_, a c Polar Angles* —-, r 6 <p • (degrees) 7.01 58 + 3 42 + 2 73+2 56.0 40.4 70.8 70.8 -36.3 8.97 60+2 42+1 67+2 59.5 41.1 65.6 65.6 -34.7 0=a ;< | > is relative to axis of the B detector. The b-axis in this system is at <j) = (-35+5)°. The error in fitted and polar angles i s +1° . 114 Without knowing the signs of the cosines in equation 5.15 one cannot determine the sign of <f>. In Table 7 the negative sign has been chosen because of the result of a separate experiment in which the A detector was at <j>^ = +125° relative to the B detector. With this geometry ± tan<j) = + cosa. - cosa,, costb. ~ A H A . cosag sintf)^ 5.16 A typical set of values from the scaler counts of this experiment, above the spin-flop, is cosa = 0.71 A cosa = 0.68 5.17 B Solving for. cf) by equation 5.16 yields = 62.5° or -29.5°. sign for <j> in Table 7 is therefore justified. The negative Thus the orientation axis is within the be plane, to within the accuracy of the detector settings. iv. Molecular Fields As the applied field along the c-axis is increased beyond the transition the spins tend to be pulled back to the c-axis according to equation 5.9. The molecular fields and can be found from this effect. From equation 5.9 one would expect that in the spin-flop phase cos8 is proportional to the applied f i e l d . By fitting a line through the two points above the transition and through the origin (figure 25) the value of 2Hg 21.4 +_ 1.0 kG. is calculated from the slope. We have 2H^ - From equation 5.11 and the transition field of = 115 FIELD Figure 25 ALONG C AXIS IN KG Cos0 versus f i e l d along c-axis above spin-flop 116 H = 6.50 +_ 0.05 kG, H (2H - H^) = H, 2 c A £ = 42.25 +_ 0.65 (kG) . 2 As a result, the fields H„ and H. are E A H A = 2.0+0.1 kG PL = 11.7 + 0.6 kG . E — 5.18 These values are in excellent agreement with the results of Miedema et a l . v. Summary Nuclear orientation provides a very direct means of observing the spin-flop transition in antiferromagnets. The transition in MnC^^r^O is seen to occur for an applied field of H at 50 mK. c = 6.5 +_ 0.05 kG There is good agreement with the simple molecular field description of the spin orientation above the transition. The molecular fields deduced from the experiments concur with previous measurements. The spins are seen to l i e in the be plane of the crystal. 117 CHAPTER SIX ORIENTATION OF 6.1 2 0 7 B i NUCLEI IN COBALT Introduction Nuclear orientation experiments have been carried out recently on bismuth dissolved in nickel and iron. 2 0 4 Bacon et al_. (1972) have oriented B i in iron and nickel and Kaplan et^ al_. (1972) have studied the same hosts. 2 0 6 B i in Bismuth i s soluble in nickel to a few atomic percent (Hansen, 1958) and only slightly soluble in iron (yQ.l' atomic percent). It has been reported that Bi in cobalt may be soluble to 1%, but this is unconfirmed, [(Vol, l'96j7). Since the solubility need only be a L small fraction of a percent to prepare a source for nuclear orientation experiments this may not be a problem. The purpose in studying Bi in cobalt i s to compare the hyperfine measurements for the ferromagnet hosts Ni, Fe and Co. The tendency of other non-magnetic atoms dissolved in these metals i s to exhibit a hyperfine field approximately proportional to the host, electronic magnetization. This i s consistent with the theory that the field arises from conduction electrons polarized by the host magnetization (Campbell, 1969) . Bacon e^t al_. (1972) measure the hyperfine field for Bi in Ni and Fe to be 325 +_ 35 and 1180 +_ 130 kG, respectively. Kaplan et al_. (1972) found the field in Ni to be 390 +_ 15 kG but could conclude for Bi in Fe only that the field was between 800 and 1000 kG. The field in cobalt is expected to be between those of Ni and Fe. 207 The magnetic moment of Bi has been measured using nuclear orientation by Kaplan et al_. (1973) and they estimate a value of 4.69 +_ 0.25p. 118 The hyperfine interaction then is quite huge. We might expect a very large anisotropy in the dominant gamma ray, which has M4 multipolarity. Problems with the experiment might.be the low solubility of Bi in Co and the reorientation of the nuclei in the intermediate isomeric state in the dominant decay sequence. 6.2 Decay Scheme 207 The Bi nucleus decays by electron capture with a h a l f - l i f e of 207 30 years to excited states of the Pb nucleus. The decay scheme i s shown in figure 26 following Lederer, Hollander and Perlman (1968). The primary decay is by f i r s t forbidden electron capture to the 1.633 MeV level of 207 Pb. This state has a h a l f - l i f e of 0.8 seconds, probably too long to ignore reorientation effects, but also too short to assume complete reorientation. It emits an M4 gamma ray of energy 10064 MeV. The secondary decay sequence i s by electron capture to the 2.330 207 MeV level of Pb followed by prompt gamma-decay to the 0.569 MeV level. The 1.761 MeV gamma ray is 94% Ml and 6% E2 (Alburger and Gunyar, 1955). 207 The 0.569 MeV level of Pb is populated by both these gamma 207 transitions and by a small amount of electron capture from Bi. The interpretation of the anisotropy of this gamma ray, which is dominated by the influence of the 1.064 MeV transition, i s rather complicated. 6.3 Expected Anisotropies The effect of the long half=life of the 1.633 MeV state on the anisotropy can be measured by observing the deviation in the measured effect from the one calculated assuming no reorientation. 119 Figure 26 Bi decay scheme 120 For the electron capture with angular momentum change L g c =2 between the states of spins 9/2 and 13/2 and for the M4 gamma transition between levels of spins 13/2 and 5/2 we have the coefficients U 2 = 0.91606 U F 2 = -0.85726 F = 0.39060 . 4 = 0.73396 4 6.1 Higher order terms are present but their effect is much smaller (^0.01). The maximum values of B and B are 1.6514 and 1.0644 respectively. The 2 4 term B is half i t s maximum for 8 = 0.7; using the magnetic moment of 2 207 4.69 y^ for the Bi nucleus and a working value for the hyperfine field in cobalt as 500 kG, then 3=0.7 corresponds to a temperature of 25 mK. At this point, W(0) = 0.37 and W(TT/2) = 1.34, easily observable effects. However, the relatively long h a l f - l i f e of the intermediate state could cause a change in the magnitude of the anisotropy calculated above. Two types of processes may operate (de Groot et a l . 1965). First, there are sudden violent disturbances due to nuclear recoil and electronic transitions in the atom after the electron capture. The effect of these disturbances i s hard to estimate theoretically but phenomenologically we can describe i t by introducing attenuation coefficients which may diminish (or enhance) the undisturbed anisotropy. Thus W(6) = E G V V B U V V P (cosO) V 6.2 where G^ allow for these changes. Secondly there will be relaxation due to spin-lattice interaction from an i n i t i a l population characteristic of the oriented Bi nucleus (but possibly attenuated by the preceding effects) to one characteristic of the 121 nuclei in the 1.633 MeV state oriented by the hyperfine interaction in the Pb atom. This type of process i s well known from nuclear magnetic resonance theory. It i s described by a relaxation time, T^, which i s typically the order of seconds or more at low temperatures. This process has been observed i n nuclear orientation, for example by Stone et^ a l. (1971) who measured the anisotropy of the gamma ray from the 40 sec. isomeric state 109 of Ag. An important feature of the spin-lattice relaxation i s that T^ is inversely proportional to temperature (except at the very lowest temperatures) . The 1.761 MeV gamma ray does not contain uncertainties due to reorientation i n the intermediate state, but i t i s a much less intense transition. It i s 10% as strong as the 1.064 MeV gamma ray and, being of higher energy, has a lower detection efficiency. The anisotropy parameters for this gamma ray are calculated using L g c = 1 for the electron capture reorientation, and 6(M1/E2) = -0.252 for the gamma transition. This mixing ratio i s chosen because i t 2 yields 6 = 0.064, consistent with Alburger and Gunyar (1955). The nega- tive sign i s selected because, as i s discussed in sections 6.5 and 6.6, we observed a positive effect from this gamma ray i n the axial direction, which result i s compatible only with the negative parameters are therefore U = 0.92496 2 gamma ray. 4 = 0.39406 The orientation terms The anisotropy U = 0.74948 2 F value. F = 0.59871 4 6.3 and B are the same as those for the 1.064 MeV 4 For & = 0.7 we have W(0) = 1.38 and W(TT/2) = 0.87. 122 The anisotropy of the 0.569 MeV gamma ray is d i f f i c u l t to estimate accurately because i t contains uncertainty not only from the reorientation of the isomer but also from the different feeding transitions. Nevertheless, because i t is such an intense peak i t s anisotropy can be easily measured to give information on sample saturation, as is discussed in section 6.4 6.4 Sample 207 The Bi activity was obtained from Amersham/Searle Corp. in 207 the form of carrier-free ^ ^ l ^ * solution.. Since BiCl^ melts at n 260°C and boils at 400°C i t would be impossible to prepare the sample by the regular method of evaporating the chloride on to a cobalt f o i l and furnacing at 700 - 900°C. The Bi ions, being more electronegative than Co, were allowed to plate spontaneously on to the surface of 0.2 g. of cobalt powder (99.99% Co) which had been introduced into the solution. The powder was used rather than a f o i l in order to increase the surface area. About 95% of the activity (^7 yCi) would be on the powder after 30 minutes. This method was suggested by R. Butters of the Metallurgy Department of the University of British Columbia. 207 The Co powder with the cm. in diameter and 4 cm. long. Bi was sealed in a quartz capsule 0.9 This length was necessary to allow the capsule to be sealed without heating the powder over 100°C. Hydrogen at 13 cm. of Hg was also in the ampoule. The Co and Bi were melted together by heating the powder to the melting point of cobalt at 1500°C. This temperature was achieved by encasing the quartz capsule within a block of graphite of length 5 cm. and 123 diameter 2 cm. and by heating the graphite to white heat by an induction furnace. The minimum power necessary to melt cobalt was determined before- hand with a non-active sample of cobalt powder. The entire length of the quartz tube had to be heated to prevent Bi from condensing on colder parts of the wall. Since quartz glass softens at 1500°C the internal pressure had to be properly adjusted. The original pressure of 13 cm. became 78 cm. at this temperature. The slight over-pressure pushed the glass walls away from the sample preventing the glass from encircling the cobalt. The graphite was maintained at 1500°C for several minutes, and was quenched to 100°C in 30 sec. by plunging i t into water. The bead of cobalt so produced had an activity of 3.7 yCi and was roughly cylindrical 207 in shape, 0.4 cm. long and 0.2 cm. in diameter. A Bi activity of 3.7 yCi in 0.2 g of cobalt corresponds to a Bi-Co alloy containing 10 ^ at.r% bismuth. No further annealing nor flattening was carried out because of the results of a study by Feldman et^ al_. (1971) which showed that for Bi in iron annealing at temperatures over 300°C tended to decrease the number of Bi atoms at lattice sites. The only further procedure applied to the sample was to remove surface activity by sandpapering. No noticeable decrease in activity was found. It was necessary that the specimen should be magnetically saturated in the f i e l d of the polarizing solenoid. Cobalt alloys generally require large fields for saturation and, in this case, there was a sizable demagnetization factor (^3.5) because of the shape of the specimen., The magnetization curve was determined by measuring the self-inductance of a small c o i l , wound around the specimen. The inductance was plotted against the field applied along the axis of the c o i l . Since the sample was approximately 124 cylindrical i t was oriented with the long dimension parallel to the coil axis and the applied f i e l d . so found. Figure 27 illustrates the saturation curve Also shown for comparison are W(ir/2) - W(0) for the 0.569 and 1.064 MeV gamma rays as a function of applied f i e l d , taken during the nuclear orientation experiment. saturation at 14 kG. A l l three curves are consistent with The superconducting polarizing solenoid in the de- magnetization cryostat could produce fields up to 16.5 kG and therefore would be adequate. 6.5 The Experiment A 1.3 uCi source of ^Co in iron was soldered to the copper link 207 next to the bead of potassium alum. Bi in cobalt. The paramagnetic salt was chrome Initially two Ge(Li) detectors were used and figure 28 shows a sample of the spectrum taken with a Ge(Li) detector. The i n i t i a l runs were performed to check the saturation of the sample and to measure the size of the anisotropy of the 1.064 MeV gamma ray. Also of interest was the sign of the effect on the 1.761 MeV gamma ray, since this determined the sign of the mixing ratio. The 1.064 MeV gamma ray showed a surprisingly small effect [W(0) = .82] at 11 mK. The 1.761 MeV gamma ray had a positive axial effect and although the statistical errors were quite large in this weak gamma ray, i t was obvious that i t too had a smaller effect than expected [W(0) = 1.1]. The reduction in the observed effects compared to the expected values could be explained by assuming that the sample was not a good solu207 tion of Bi in cobalt and that only a small fraction (M0%) were in the lattice sites. of the Bi To confirm this view the temperature dependence 125 Ar G A M M A RAY ANISOTROPY AT 11 MK 1.064 MEV o s £ I i * 0.569 MEV 1 2 , MAGNETIC SATURATION AT 300 K 1.0| i LU.8f 1.6| o UJ O . z ,< to. Q J 6 APPLIED Figure 27 Saturation of L_ 8 L 10 FIELD - KILOGAUSS 207 Bi/Co sample 12 14 8r 126 0.569 MEV GE(LI) DETECTOR X10 1.064 1.761 1.173 CO o 1.332 UJ • z z < I o cr, 569 UJ CL to • • MEV 40 Z3 NA I DETECTOR o o 30 X10 1.332 1.064 20 1.761 10 100 50 Figure 28 150 Sample ^ Bi/^°Co spectra 2 7 200. 250 127 of the anisotropics would be needed. Using the Ge(Li) detectors the 1.064 MeV gamma ray was observed within a temperature range of 11 to 70 mK. After demagnetization the temperature was 11 mK. Higher temperatures were produced in a reversible manner by restoring a small current (up to 2 amp) to the large magnet. At each temperature four 1000 sec. counts were acquired. Ten temperature points were taken, including two at the lowest temperature (zero field from the large magnet) taken before and after the others. The study of the temperature dependence of the 1.761 MeV gamma ray was conducted with a 5 in. x 5 in. Nal(Tl) detector in the axial direction at a distance of 25 cm. from the source. The increase in the count rate over that of the Ge(Li) detector was a hundredfold producing a tenfold decrease in the statistical error. Because of the poorer resolu- tion of this detector i t was necessary to reduce the activity of ^Co to 0.2 yCi so that a cleaner spectrum could be obtained for analysis of the 1.761 MeV peak. A sample spectrum from this detector is given in figure 28. The run to determine the temperature dependence of the 1.761 MeV gamma ray was carried out with the same procedure as was usedlfoflthe •!.061MeV study. Four 1000 sec. counts at each of ten temperature points were acquired. The temperature was increased by putting some current into the large magnet. A Ge(Li) detector was placed in the equatorial direction to observe the ^Co gamma rays for temperature measurement. These peaks could not be resolved by the Nal detector. 128 6.6 Analysis The gamma rays were analyzed by assuming the extremely simplified 207 model that a fraction f of Bi nuclei were in lattice sites of the cobalt host, feeling the f u l l hyperfine f i e l d , and that the remaining nuclei felt an average of zero f i e l d . [W(6) - l ] Thus = f[W(0) - 1] Q b s that i s , the observed effect is reduced by the factor f. 6.4 In addition the 1.064 MeV gamma ray may be altered by the reorientation as discussed in section 6.3. In non-magnetic metals there exists a relationship between the Knight shift and the spin-lattice relaxation time, T^. This i s the 'Korringa relation' (Korringa, 1950) o AH (TjT) YfjFjj) Here V = constant . 6.5 is the nuclear gyromagnetic ratio and the constant involves funda- mental atomic constants. The Knight shift depends on the polarization of conduction electrons by the applied magnetic f i e l d , and we might expect a similar relation to apply for hyperfine interactions in ferromagnetic metals. Thus we assume 2 2 (T-jT) Y H £ soeonstant . n ef 6.6 Data supporting this assumption has been given by Stone (1971) and theoretical 129 arguments supporting i t have been discussed by Campbellf(1969). Experi- 18 2 mental evidence indicates that the constant has a value ^2 x 10 G sec K. Most of the data is relevant to impurities dissolved in iron, butOCampbell predicts should be independent of the host material. If we apply the relation to the case of Pb in cobalt we obtain T T = 7 sec K, 6.7 1 where we have used the H £^ value for Pb in Co of 262 kG (Koster and Shirley, 1971) and assumed a magnetic moment of 2 y^ for the isomeric state. The hafftlife of this state is 0.8 sec. so that for the temperature range investigated spin-lattice effects may be ignored. In equation 6.2, i f the reorientation factors a n ^ G4 e a r gqual (to G, say) then [W(0) - l ] o b s = fG[W(6) - 1] . 6.8 With this very simple model the temperature dependence of the gamma rays could be fitted to the expected anisotropy to find the hyperfine field and the attenuation^fG. For the 1.761 MeV gamma ray the attenuation is the fraction of nuclei in lattice sites, and for the 1.064 MeV transition, i t is that fraction times the reorientation factor. The 1.064 MeV peaks were analyzed by the program SPECTAN in the normal way. Temperatures were measured by the ^Co thermometer. The 1.761 MeV gamma ray was studied using a Nal detector with such a wide FWHM that the linear approximation for the background was not considered suitable. 130 Instead a l l spectra of the same temperature were added and plots were made of the summed spectra. These were analyzed by fitting the continuous background on both sides of the peak with a flexible curve, drawing in the portion under the peak, and evaluating the remaining area under the peak by a planimeter. We estimate the total error from statistics, background f i t t i n g and planimeter reading to be vl.5% for each of the cold spectra. In figure 29 are shown the plots of the axial and equatorial effects from the 1.064 MeV peak and the axial effect from the 1.761 transition. MeV The equatorial effect for the latter i s not shown because i t was small (^ 1-2% at 11 mK), a result which is consistent with the value of f = 12% found from the axial effect. These data were fitted to the expected anisotropies (with solid angle corrections) with the attenuation and the hyperfine field as unknown parameters. The fitted values are given in Table 8 and the curves are shown in figure 29. TABLE 8. Fitted values of attenuation and hyperfine from the gamma rays . oft 207Bi. p Attenuation Hyperfine field kG ^ Deg. of Freedom 1.064 MeV Axial .176 +_ .005 560 + 30 6.1 8 1.064 MeV Equat. .186 +_ .007 880 +_ 65 13.5 8 1.761 MeV Axial .12 590 + 200 10.2 7 + .025 132 6.7 Discussion For the interpretation of the effect from the 1.761 MeV gamma ray, the model which assumes that there is no contribution from nuclei not in lattice sites is a very extreme one. However, our analysis indicates that this approximates to the actual situation. Of course, more likely there ise a large fraction of nuclei which feel weak fields (< 100 kG) and a smaller fraction ( M2%) which feel a much larger f i e l d of ^600 kG, of perhaps a distribution of fields of this order of magnitude. Thus we conclude from the f i t for the 1.761 MeV gamma ray that 12% of the nucleiaare in fields in the range 400 to 800 kG. For the 1.064 difficult. follows: MeV gamma ray the interpretation is more The significance of the fitted parameters in table 8 is:.as the attenuation factor determines the ratio of the observed and expected effects at saturation, whereas the hyperfine field determines the temperature range in which the saturation begins to occur. equality of the attenuation coefficients for the 1.064 MeV axial and equatorial effects indicates that the simple model in which taken as equal is reasonable. The near and G^ were We see that G = 1.6 ± 0.4, andenhancement. However, the discrepancy in the fitted values of the hyperfine field suggests that the previously ignored spin-lattice relaxation effects are important at the higher temperatures, altering the values of B^ and B^. The axial and equatorial effects would not be changed to the same degree because of the difference in the relative signs of the axial and equatorial directions. and in The observed discrepancy could be explained by a value of T^T of approximately 0.2 sec K, a much lower value than that estimated in section 6.6. 133 The nuclear orientation of 206 advantages over the present experiment. Bi in cobalt would have several Because the h a l f - l i f e i s 6 207 days instead of the 30 years for Bi one can achieve the same activity 207 with 1/1500 the number of atoms needed for. Bi. This reduction in concentration may allow for a larger fraction, perhaps 100%, of nuclei at lattice sites. Secondly, there are no long-lived intermediate states 206 ^0207 in the decay of Bi. A third point i s that as"'for h.> Bihitheedecays is by electron capture; thus there are no problems about heating of the source from absorbed beta particles. The major disadvantage is the short h a l f - l i f e of 6 days which would have made i t d i f f i c u l t for experimenting with source preparation. 207 However, now that the technique has been developed with Bi there 206 should be few problems in preparing sources with Bi. 134 REFERENCES Abragam, A. and Pryce, M.H.L. 1951. Proc. Roy. Soc. Ser. A, 205, 135 Agarwal, Y.K., Baba, C.V.K. and Bhattacherjee, S.K. 1967. A99, 457. Alburger, D.E. and Gunyar, A.W. 1955. Nucl. Phys. Phys. Rev. 99, 695. Ambler, E., Grace, M.A., Halban, H., Kurti, N., Durand, H., Johnson, C. and Lemmer, H.R. 1953. P h i l . Mag. 44_, 216. Avignone, F.T. and Frey, G.D. 1971. Phys. Rev. C, 4_, 912. Bacon, F., Haas, H., Kaindl, G. and Mahnke, H.E. 1972. Phys. Let. 38A, 401. Berstein, S., Roberts, L.D., Stanford, CP., Dabbs, J.W.T. and Stephenson, T.E. 1954. Phys. Rev. 94, 1243. Biedenharn, L.C. and Rose, M.E. 1953. Bitter, F. 1937. Phys. Rev. 54_, 79. Blair, A.G. and Armstrong, D.D. 1965. Bleaney, B. Rev. Mod. Phys. 25_, 729. Phys. Rev. 140, B1567. 1951. Proc. Phys. Soc. A, 64_, 315; P h i l . Mag. 42, 441. Blin-Stoyle, R.J. and Grace, M.A. 1957. Handbuch der Physik, Vol. 42 (Springer-Verlag, Berlin) p. 555. Campbell, I.A. 1969. J . Phys. C:Solid State Phys., 2_, 1338. C o l l i n , W., Daniel, H., Mehling, 0., Schmitt, H., Spannagel, G., and Subdhi, K. S. 1964. Z. Physik, 180, 143. Dabbs, J.W.T., Roberts, L.D. and Berstein, S. 1955. Daniels, J.M., Giles, J.C., Leblanc, M.A.R. 1961. Daniels, Phys. Rev. 98, 1512. Can. J . Phys. 39_, 53. J.M., Grace, M.A. and Robinson, F.N.H. 1951. Nature 168, 780. DeGroot, S.R., Tolhoek, H.A. and Huiskamp, W.J. 1965. In Alpha-, Beta-, Gamma-Ray Spectroscopy, Vol. 2, (K. Siegbahn, editor, NorthHolland Publishing, Amsterdam) Ewan, G.T. and Graham, R.L. 1965. In Alpha-, Beta-, Gamma-Ray Spectroscopy, Vol. 2, (K. Siegbahn, editor, North-Holland Publishing, Amsterdam). Feldman, L.C., Kaufmann, E.N., Mingay, D.W. and Augustyniak, W.M. Phys. Rev. Let. 27_, 1145. 1971. Ferentz, M. and Rosenzweig, N. 1955. Tables of F coefficients, Argonne National Laboratory Report #5324. 135 Flynn, K.F., Glendenin, L.E. and Steinberg, E.P. 1965. 22_, 416. Nucl. Sci. Eng. Frauenfelder, H. and Steffen, R.M. 1965. In Alpha-, Beta-, Gamma-Ray Spectroscopy, Vol. 2 (K. Siegbahn, editor, North-Holland Publishing, Amsterdam). Gijsman, H.M., Paulis, N.J. and Van der Handel, J. 1959. Physica, 25, 954. Gorling, R.L.A. 1970. Master's Thesis, University of British Columbia (Unpublished). Gorter, C.J. 1948. Physica, 14_, 504. Gorter, C.J., Poppema, O.J., Steenland, M.J., and Beun, J.A. 1951. 17_, 1050. Physica, Grace, M.A., Johnson, C.E., Kurti, N., Scurlock, R.G., and Taylor, R.T. 1955. Comm. Conf. Phys. Basses Temp. (Union Int. de Physique Pure et Appl., Paris) p. 263. Gray, T.P. and Satchler, G.R. 1955. Proc. Phys. Soc. A68, 349. Hanna, S.S., Herberle, J., Perlow, G.J., Preston, R.S. and Vincent, D.H. Phys. Rev. Let. 4_, 513. Hansen, M. 1958. Constitution of Binary Alloys. New York). 1960. (McGraw-Hill Book Co., Jeffries, CD. 1963. Dynamic Nuclear Orientation, (Interscience Publishers, John Wiley and Sons, New York). Johnson, C.E., Schooley, J.F. and Shirley, D.A. 1960. Phys. Rev. 120 , 1777. Kaplan, M., Johnston, P.D., A41, 315. Kittel, P. and Stone, N.J. 1972. Phys. Let. Kaplan, M., Johnston, P.D., communication). Kittel, P. and Stone, N.J. 1973. (Private Karlsson, E., Matthias, E. and Siegbahn, K. (editors). 1964. Angular Correlations (North-Holland, Amsterdam). Kistner, O.C. 1966, Perturbed Phys. Rev. 144, 1022. Kogan, V., Kulkov, V.D., Nitikin, L.P., Reinov, N.M., Sokolov, I.A., and Stelmah, M.F. 1961. Soviet Physics JETP, 12_, 34. Korringa, J. 1950. Physica 16_, 601. Koster, T.A. and Shirley, D.A. 1971. Hyperfine Interactions in Excited Nuclei, (G. Goldring and R. Kalish, editors, Gordon and Breach, New York) Vol. 4, p. 1239. Krane, K.S., Olsen, C E . Sites, J.R. and Steyert, W.A. 4, 1906. 1971. Phys. Rev. C, 136 Kul'kov, V.D., Kogan, A.V., Nikitin, L.P., Savin, E.P., and Stel'Makh, A.F. 1965. J. Expt. Theoret. Phys. 48, 122. (Eng. Trans, in Soviet Physics JETP, 48, 83.) Lederer, CM., Hollander, J.M. and Perlman, I. 1968. (John Wiley and Sons, New York). Marshall, W. 1958. Table of Isotopes, Phys. Rev. 110, 1280. Matthias, E. and Holliday, R.J. 1966. Metcalfe, M.J. 1971. Phys. Rev. Lets. 17, 897. Phys. Letters A, 36, 373. Miedema, A.R., Wielinga, R.F. and Huiskamp, W.J. 1965. Physica, 31_, 835. Mossbauer, R.L. 1965. In Alpha-, Beta-, Gamma-Ray Spectroscopy, Vol. 2 (K. Siegbahn, editor, North-Holland Publishing, Amsterdam). Nagamija, T., Yoshida, K. and Kubo, R. 1955. Neel, L. 1932. Ann. de Phys. 18_, 5. Neel, L. 1936. Ann. de Phys. 5_, 232. Advan. Phys. 4_, 1. Nordhagen, R., Elbek, B. and Herskind, B. 1967. Pettersson, H., Antman, S. and Grunditz, Y. 1970. Portis, A.M. and Gossard, A.C. 1960. Pound, R.V. 1949. Nucl. Phys. A104, 353. Z. Physik, 233, 260. J. App. Phys. 31_, 205 S. Phys. Rev. 76_, 1410. Raeside, D.E., Reidy, J.J. and Wiedenbeck, M.L. 1969. 347. Nucl. Phys. A134, Roberts, L.D., Dabbs, J.W.T., and Parker, G.W. 1958. Proc. 2 United Nations Int. Conf. on the Peaceful Uses of Atomic Energy, 1_5 (United Nations, Geneva) p. 322. n d Rose, M.E. 1949. Phys. Rev. 75, 213. Rose, M.E. 1965. In Alpha-, Beta-, Gamma-Ray Spectroscopy, Vol. 2 (K. Siegbahn, editor, North-Holland Publishing, Amsterdam). Rose, M.E. 1967. Elementary Theory of Angular Momentum (John Wiley and Sons, New York). Samoilov, B.N., Sklyarevskii, V.V., and Stepanov, E.P. 1960. JETP 11_, 261. Soviet Physics Shirley, V.S. 1971. Hyperfine Interactions in Excited Nuclei, (G. Goldring and R. Kalish, editors, Gordon and Breach, New York) Vol. 4, p. 1255. 137 Shirley, D.A. and Westenbarger, G.A. 1965. Sidu, N.P.S. and Gupta, U.C. 1967. Phys. Rev. 138, A170. Nucl. Phys. A91, 557. Stone, N.J. 1971. Hyperfine Interactions in Excited Nuclei, (G. Goldring and R. Kalish, editors, Gordon andiVBreach,'oNew York) Vol. 1, ;pvr>237. ' . . • f. • \ /Si Stone, N.J., Fox, R.A., Hartmann-Boutron, F. and Spanjaard, D. 1971. Hyperfine Interactions in Excited Nuclei (Goldring and Kalish, editors; Gordon and Breach Science Publishers, New York), Vol. 1, p. 351. Stone, N.J. and Turrell, B.G. 1962. Phys. Letters 1_, 39. Templeton, J.E. and Shirley, D.A. 1967. Tschanz, J. F. and Sapp, R.C. 1970. Phys. Rev. Lett. 1_8, 240. Phys. Rev. C, 2_, 2168. Turrell, B.G., Johnston, P.D. and Stone, N.J. 1972. State Phys. 5_, 197. Van Vleck, J.H. 1941. J. Phys. C: Solid J. Chem. Phys. 9_, 85. Vol, A.E. 1967. Handbook of Binary Metallic Systems, Vol. 2. Program of Scientific Translations, Jerusalem). White, K.W. (Israel 1968. Experimental Techniques in Low Temperature Physics, (Oxford University Press, London). Yamazaki, T. 1967. Nuclear Data, A3, 1.
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Nuclear orientation studies at low temperatures
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Nuclear orientation studies at low temperatures Daly, Patrick William 1973
pdf
Page Metadata
Item Metadata
Title | Nuclear orientation studies at low temperatures |
Creator |
Daly, Patrick William |
Publisher | University of British Columbia |
Date Issued | 1973 |
Description | The interaction between the nuclear magnetic moment and the hyperfine field in ferromagnetic and antiferromagnetic materials has been used to orient radioactive nuclei at low temperatures. The measurement of the resulting angular distribution of the gamma radiation has yielded a variety of information in both nuclear and solid state physics. The ⁵⁹Fe nucleus in an iron sample has been studied in an attempt to verify a report by Tschanz and Sapp (1970) who measure an unexpectedly large effect in the gamma ray distribution from that nucleus oriented in double nitrate salts. Since in their experiment assumptions have to be made regarding the fraction of iron atoms in lattice sites, the well-defined environment of the iron lattice was chosen for the present experiment in order to permit an unambiguous interpretation of the results. In contradiction to predictions based on the work of Tschanz and Sapp, a null result was obtained at a temperature of 15 mK, indicating an upper limit for the magnetic moment of ⁵⁹Fe of |µ| < 0.9 µ[sub N. The anisotropy of gamma radiation emitted in the decay of ¹⁰³Ru oriented in ferromagnetic iron was measured to determine the magnetic moment of the isotope and to resolve mixing ratios in the decay scheme. The small effect observed at 11 mK permitted a lower limit to be placed on the magnetic moment: |µ| > 0.15 µ[sub N]. The sign of the E2/M1 mixing ratio for the 497 keV gamma transition was determined for the first time to be found negative. The techniques of nuclear orientation have also been used to observe the spin-flop transition in antiferromagnetic MnCl₂H₂O. By observing the gamma ray anisotropy in the decay of radioactive ⁵⁴Mn we find that at 60 mK the transition occurs for a critical field of H[sub c] = 6.50 ± 0.05 kG applied along the crystallographic c-axis. In the spin-flop phase, the electronic spins point at an angle of 71° to this axis. These data allow values of the molecular exchange and anisotropy fields to be determined yielding respectively H[sub E] = 11.7 ± 0.6 kG and H[sub A] = 2.0 ± 0.1 kG in agreement with the values reported by Miedema et al. (1965). This experiment demonstrates that nuclear orientation is a viable technique for studying spin-flop transitions and moreover allows the measurement of the directions of electronic spin magnetization in the spin-flop phase. The nuclear orientation of ²⁰⁷Bi in cobalt was studied to determine the hyperfine field. An unexpectedly small gamma ray anisotropy was observed. However, saturation of the effect occurred at a relatively high temperature. These results indicate that only a fraction (~12%) of the Bi nuclei felt a large hyperfine field (~600 kG). Presumably this was the fraction of atoms in lattice sites. Comparison of the anisotropics of the 1.064 MeV and the 1.761 MeV gamma rays shows that the former effect is enhanced, possibly because of the long half-life of the 1.633 MeV intermediate state. |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2011-03-02 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0084949 |
URI | http://hdl.handle.net/2429/31939 |
Degree |
Doctor of Philosophy - PhD |
Program |
Physics |
Affiliation |
Science, Faculty of Physics and Astronomy, Department of |
Degree Grantor | University of British Columbia |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-UBC_1973_A1 D34.pdf [ 5.65MB ]
- Metadata
- JSON: 831-1.0084949.json
- JSON-LD: 831-1.0084949-ld.json
- RDF/XML (Pretty): 831-1.0084949-rdf.xml
- RDF/JSON: 831-1.0084949-rdf.json
- Turtle: 831-1.0084949-turtle.txt
- N-Triples: 831-1.0084949-rdf-ntriples.txt
- Original Record: 831-1.0084949-source.json
- Full Text
- 831-1.0084949-fulltext.txt
- Citation
- 831-1.0084949.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0084949/manifest