UBC Theses and Dissertations

UBC Theses Logo

UBC Theses and Dissertations

Muonium in some insulating oxides and diamond Spencer, David Philip 1985

Your browser doesn't seem to have a PDF viewer, please download the PDF to view this item.

Item Metadata


831-UBC_1985_A1 S63.pdf [ 6.61MB ]
JSON: 831-1.0060480.json
JSON-LD: 831-1.0060480-ld.json
RDF/XML (Pretty): 831-1.0060480-rdf.xml
RDF/JSON: 831-1.0060480-rdf.json
Turtle: 831-1.0060480-turtle.txt
N-Triples: 831-1.0060480-rdf-ntriples.txt
Original Record: 831-1.0060480-source.json
Full Text

Full Text

M U O N I U M IN SOME INSULATING OXIDES A N D D I A M O N D  by DAVID PHILIP SPENCER B.Sc, The University of British Columbia, 1977  A THESIS SUBMITTED IN P A R T I A L F U L F I L M E N T OF T H E REQUIREMENTS FOR T H E D E G R E E OF DOCTOR OF PHILOSOPHY  in T H E F A C U L T Y OF G R A D U A T E STUDIES Department of Chemistry  We accept this thesis as conforming to the required standard  UNIVERSITY OF BRITISH C O L U M B I A February, 1985 ®  David Philip Spencer, 1985  In piesenting this thesis in partial fulfilment of the requirements for an advanced degree at. the UNIVERSITY O F BRITISH C O L U M B I A , I agree that the Library shall make it freely available for reference and study.  I further agree that permission for extensive copying of this thesis for  scholarly purposes may be granted by the Head of my Department or by his or her representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission.  Department of Chemistry UNIVERSITY O F BRITISH C O L U M B I A 2075 Wesbrook Place Vancouver, Canada V6T 1W5  Date: February. 1985  Absnaci Muonium can be regarded as an isotopic analogue of the hydrogen atom, with a positive muon replacing the proton.  This thesis is concerned with searches for muonium in weak  transverse and zero magnetic Fields in a variety of non-magnetic oxide insulators with few nuclear moments, and in diamond, a semiconductor.  The technique of Muon Spin Rotation (uSR) was  used, which can be thought of in analogy with N M R and EPR.  It relies on measuring the  direction of the muon spin with time via the unique signature of the parity violating decay of the muon, in which its decay positron is emitted preferentially along the muon spin direction.  This  can provide a sensitive measure of the interaction of the muon spin with its environment. Muonium has been searched for successfully in a-quartz, fused quartz., hexagonal and fused germanium dioxide, magnesium oxide, and diamond; and unsuccessfully in tetragonal germanium dioxide, rutile, strontium titanate, zircon, and beryl. All the samples save the quartzes showed a large "missing fraction"; that is, not all the initial muon polarization could be accounted for.  This indicates that muonium was formed but rapidly depolarized.  The missing fractions  were found to be very strongly temperature dependent in strontium titanate but only weakly in rutile.  At room temperature, the observed muonium fractions were: a-quartz and fused quartz  (65+5)%, hexagonal germanium dioxide <2%, fused germanium dioxide (10±4)%, magnesium oxide (30± 10)%, and diamond (33± 4)%.  It is noteworthy that the fractions are very different in  quartz and germanium dioxide, despite their structural analogy.  The observation of muonium in  diamond represents the first time that muonium has been seen in a room temperature semiconductor. A spin-Hamiltonian new to uSR was found in low-temperature a-quartz and hexagonal germanium dioxide, in which the  hyperfine interaction is completely anisotropic.  In  low-temperature a-quartz in zero magnetic field this interaction gives rise to three signals whose frequencies are orientation independent, but whose amplitudes are strongly orientation dependent By studying the orientation dependence it was determined that muonium occupies a site which is identical with that seen for hydrogen in EPR studies. ii  This was the first time that muonium was  shown to behave just like an isotope of hydrogen in a solid.  In hexagonal Ge0 , eight frequencies 2  were observed at 6K, in marked contrast to the observation of only three frequencies in structurally analogous a-quartz. The interpretation of this is that there are three or four different sites, the exact nature of which has not been determined.  iii  lents  A.  B.  C.  D.  E.  Chapter I: Introduction  1  1.  The Muon and its Polarization  .'  1  2.  The Muonium Atom  5  3.  Muon and Muonium Fractions  8  4.  MSR In Solids  9  Chapter II: Experimental Apparatus  15  1.  Cyclotron, Beamlines and D.C. Separator  15  2.  Surface Muon Beams and Conventional Beams  17  3.  Data - Collection and Analysis  18  4.  Eagle: a Standard Surface Muon Apparatus  26  5.  Some Considerations in Apparatus Design  28  Chapter III: Muonium Hyperfine Interactions  38  1.  Spin-Hamiltonian I: Isotropic Hyperfine Interaction  38  2.  Spin-Hamiltonian II: Axially Symmetric Hyperfine Interaction  41  3.  Spin-Hamiltonian III: Completely Anisotropic Hyperfine Interaction  43  4.  Superhyperfine Spin-Hamiltonian  46  Chapter IV: A Brief Survey of Muonium Formation and Relaxation  48  1.  Gases  49  2.  Liquids  50  3.  Solids  51  Chapter V : Results and Discussion  56  1.  Quartz  56  2.  Germanium Dioxide  85  3.  Magnesium Oxide  92  4.  Strontium Titanate  96  5.  Rutile  103  iv  6.  Beryl and Zircon  106  7.  Diamond  106  F.  Chapter VI: Summary and Conclusions  112  G.  Appendix 1: A Muon Magic Beamline for T R I U M F  116  1.  Introduction  116  2.  State of the Art  116  3.  Alternatives  116  4.  Justification  118  5.  Design Considerations  119  6.  Conclusions  122  H.  I.  Appendix 2: Solving the Anisotropic Hyperfine Hamiltonian  123  1.  Solve H , the Hyperfine Part  124  2.  Solve the Zeeman Part  126  3.  Solve the Total Hj  127  A  Appendix 3: The Superhyperfine Interaction Hamiltonian  BIBLIOGRAPHY  146 147  vj.  List of Figures 1. A typical uSR time histogram  3  2. A typical uSR asymmetry plot  4  3. Muonium in fused quartz at 10 G ......  7  4. The muon signal in fused quartz at 85 G  10  5. M20 beam line  19  6. A typical uSR apparatus  .....  7. Comparison of W S F T and F F T  20 —  25  8. Interior of Eagle  27  9. Sketch of the variable collimator  29  10. Dependence of number and asymmetry of positron decays on energy  33  11. A N as a function of cutoff energy  34  12. B r e i t - R a b i diagram for the isotropic case  39  13. B r e i t - R a b i diagram for the cylindrically symmetric case  42  14. B r e i t - R a b i diagram for the completely anisotropic case  45  15. L o w - f i e l d beating of muonium in o -quartz  62  J  16. Anisotropy of splitting in quartz ....  .......................................................—.—............................— 63  17. Muonium oscillation at zero applied magnetic field in a -quartz  64  18. Fourier transforms of low temperature a - quartz data  66  19. Diagram of the hydrogen site in a - q u a r t z  70  20. High statistics muonium signal  72  21. Fourier transforms of quartz data, unapodized and weakly apodized  74  22. Fourier transforms of quartz data, medium and strong apodizations  75  23. W S F T of high statistics quartz data  ....................... 76  24. Relaxation rate versus temperature for a - quartz  —  79  25. Relaxation rate of the muonium signal in fused silica versus temperature . 26. Muonium signal in the perpendicular telescope in zero applied magnetic  vii  82 field  84  27. Relaxation rate of the muonium signal in fused germania versus temperature  87  28. F F T o f fused germanium dioxide in a 60 Gauss magnetic  88  29. F F T and W S F T of hexagonal G e 0  field  89  2  30. Muonium signal in room temperature M g O  93  31. Relaxation rate of the muonium signal in M g O versus temperature  .'.  94  32. Amplitude of the muon signal in SrTiOj versus temperature  97  33. M u o n relaxation rate in S r T i 0 versus temperature  99  3  34. Asymmetry plot of early time data in SrTiOj at zero and seven Gauss  101  35. Plot o f early time data in SrTiOj at 25 Gauss and in the perpendicular direction  102  36. Temperature dependence of the muon signal in T i 0  104  37. Plot o f the observed signal in T i 0  2  at 20 K and 295 K  105  38. Asymmetry in diamond at room temperature and 6 K  108  39. Relaxation rate of muonium i n diamond versus temperature  109  2  viii  List of Tables i. Some properties of the positive muon and of muonium  6  ii. Summary of results obtained  57  iii. Purities of materials used  58  iv. Comparison of muonium and hydrogen results in quartz  67  v. Frequencies observed in G e 0  90  2  viii  Acknowledgement'? I wish to thank my supervisors, D o n Fleming and Jess Brewer for their efforts, both mental and physical, over, the years it has taken me to complete first my research, and then my thesis, and I hope that they feel it was all as worthwhile as I do.  I must also single out John Coope for his  assistance in showing me how to do the calculations of Appendix 2, naturally any errors or omissions are my own.  Peter van den Bosch and Jon Nightingale were generous with their time  and patience in making the document itself look as attractive as I believe it to be. There are far more people I feel I should acknowledge than I could ever give individual attention to, so...  I wish to thank Jim Mcllroy and his "Beamline Boys" for their eternal  cheerfulness and willingness to jump to my aid at moments of crisis; D o n Arseneau, Masa Senba, Dave Garner, Dale Harshman, Carl Clawson, Eric Turner, Paul Percival, T o m o Uemura, R o b Kiefl, R o l f Keitel, G l e n n McRae, Randy Mikula, Bob Butters, T J . Ahrens, Joe Nagel, J. Weil, Y . Koga, and N . Drikitis. I will single out my parents, for whose assistance I can never possibly be grateful enough and my p a r e n t s - i n - l a w , who have been of great help in seeing me through this.  Most of all, to  my wife Sylvia, who has put up with so much over the years, who has helped me in so many ways, and to whom I dedicate this, my thesis.  1  A . C H A P T E R I: I N T R O D U C T I O N  1. T H E M U O N A N D ITS P O L A R I Z A T I O N The muon is a sub-atomic particle with about 200 times the rest mass of the electron, or about l / 9 t h the rest mass of the proton.  It has a charge of either +1  or - 1 , decaying into a  positron and two neutrinos, or an electron and two neutrinos respectively.  This work will be  concerned exclusively with positive muons.  u*  >  e  +  +  v  +7  (1.1)  u  e The muon has a mean life o f about 2.2 microseconds.  '  This may seem short, but its parent, the  pion, has a mean life of only 26 nanoseconds:  rr  +  > u* + v  (1.2)  u  The muon lifetime is long enough though so that detecting it and measuring the elapsed time until the observation of its decay positron is relatively easy.  Moreover, many interesting subjects can be  studied in the microsecond range available. The muon can be used in studying a wide variety of subjects, including charge exchange in gases [Fleming 84] \ relaxation dynamics in both gases [Senba 84] and liquids [Stadlbauer 84], magnetic fields and relaxation phenomena in solids [Heffner 84], diffusion on surfaces [Harshman 84), the Lamb shift of the muonium atom in vacuum [Oram 84], and is used in particle physics experiments [Carr 84].  The acronym invented to describe the field is u S R or M S R for M u o n Spin  Rotation (Resonance, Research...).  The technique and its applications have been reviewed by  Brewer et al. [75], Schenck [76, 81], and Heffner et al. [84]  and has been the subject of three  conferences, the proceedings of which have been published in Hyperfine Interactions [Gygax 79, Brewer 81b, Yamazaki 84].  ' T h e format which will be used for references is [first author's last name year of publication modulo 1900], unless the author's name has immediately preceded the reference in the texL in which case only the year modulo 1900 will appear.  2 What is it that makes the muon such a powerful tool?  A very simple property - it is  possible to determine the direction of the spin of a muon as a function of time.  This is due to the  non-conservation of parity in the weak interaction [Lee 56] which has two effects: the muon is produced 100% polarized in the direction of its momentum, and when it decays the positron produced tends to come off along the muon spin direction. In practice, due to quantum mechanics for one thing, it is not the spin of an individual muon which is followed, but rather the spin of an ensemble of muons. ensemble is generated by studying one muon at a time.  Generally, however, this  That is, in a typical experiment, a muon  stops in a target and begins precessing in a magnetic field applied transverse to the spin direction. 2  After some time it decays.  The direction in which the decay positron is emitted and the elapsed  time since decay are determined. Another muon stops in the target, and so on. After perhaps 100 million muons have been studied, the result is as shown in Figure 1, where the number of positrons detected in a particular direction is plotted against elapsed time until decay.  The obvious  exponential reflects the muon lifetime, while the oscillation is due to the precession of muons in the applied field . The amplitude of this oscillation as a fraction of the number of counts is termed 3  the asymmetry and is plotted in Figure 2  4  .  In magnetic materials such as iron, the muon may  also precess in a local magnetic field. This phenomenon was first observed by R.L. Garwin et al. [57], in testing parity conservation, an experiment suggested by Lee and Yang [56]. 5  Garwin et al. stopped muons in  two different targets: graphite and nuclear emulsion. The amplitude of the oscillations observed in the nuclear emulsion was markedly less than that in graphite.  Further experiments a year later by  Strictly speaking, the modern interpretation of the meaning of quantum mechanics prohibits speaking of such things as the precession of a single spin. However, its usage here reflects convenience as well as continued use within the field. Whenever the phrase "applied field" is used, it should be taken to mean "applied 2  3  magnetic  field".  This spectrum is a combination of spectra from two different directions, the left and the right, hence it is called a Left-Right spectrum. Backward-Forward and U p - D o w n spectra may also be formed. The directions used are usually from the muon's point of view. It was the second experiment to show parity violation, immediately following the classic C o /3 decay experiment of Wu [57] et a l ; indeed Garwin's paper immediately followed Wu's in Physical Review's Letters to the Editor. 4  5  60  •3  1600 r 8. t  1400 h  3 «T>  55' s  1200  h  o w  CO  1000 h  w  to I-* ON O  o c  Ok CO P O U  8  T I M E  (Microsec)  4  Figure 2 A typical uSR asymmetry ploL Aluminum at 60 Gauss.  5 Swanson [58] in a variety of materials revealed the same wide variations in the "residual" muon polarization, from about 10% in benzene to 100% in typical metals.  In non-magnetic materials,  such differences are now known to be due to the thermalization process of the positive muon in matter.  2. THF, MI IONIUM A T O M When a positive muon passes through, or stops in, matter, it may pick up an electron to formau e~ +  bound state called muonium (Mu)  6  '.  Since the reduced masses of muonium and  hydrogen are only very slightiy different, properties such as ionization potential and radius are virtually identical, as shown in Table I.  Consequently it behaves chemically like an isotope of  hydrogen but one with a far greater mass difference than is available from deuterium or tritium. When the muon captures an electron, their spins may be either parallel or anti-parallel.  If they  are anti-parallel, the muon's spin will oscillate at the hyperfine frequency of 4463 M H z (for muonium in a vacuum) and be unobservable by conventional M S R techniques. This results in only one-half the muonium ensemble being directly observable. Any muonium formed may be readily distinguished from unreacted muons and/or muons in diamagnetic environments, as the characteristic Larmor precession frequency of muonium in a weak transverse field is about 103 times larger than that of the muon.  This can be seen by  comparing Figure 2, a plot of the muon asymmetry in aluminum at 60 Gauss, to Figure 3, a plot of muonium asymmetry in fused quartz in 10 G .  Several phenomena can be observed:  1.  the probability of muonium formation;  2.  the rate at which the muonium signal vanishes;  3.  the mechanism for its disappearance; and  4.  the strength and directionality, if any, of the hyperfine coupling between the muon and its  'The spin of the electron may be either parallel to the muon spin, thus forming "triplet" muonium, or it may be anti-parallel, forming the so-called mixed muonium. 'The first written reference to this is in a paper by Friedman and Telegdi [57], which was published in the same volume of the same journal as Wu's and Garwin's papers. The dates indicate that the three papers were received at virtually the same time.  Table I Muon  Some properties of the positive muon and o f muonium  - Charge  =  +1  Spin  =  1/2  Rest Mass  = 105.6596 M e V = 206.7685 m e = 0.1126 m  P Magnetic Moment  = 4.4905x10 " e r g G " = 3.1833 <z  p  = 0.004836 u c  g-Factor = 2.0023318 = 1.000006 g Mean Lifetime Gyromagnetic Ratio  e = 2.197037 MS = 13.5544 k H z G "  Muonium  Rest Mass  1  = 3.18334 0.1131 m7 , , H = 207.8 m  Reduced Mass  = 0.9952 m . = 0.9956 M  Bohr Radius  = 0.5315 A = 1.0044 a o  Ionization Energy  H  H  = 13.54 eV = 0.99561.E  Hyperfine Frequency  H  = 2.8044 x 1 0 " rad s" = 3.1423 WOJJ = 4463 M H z  1  Figure 3 Muonium in fused quanz at room temperature in a field of 10 G .  8 election. These subjects, in solids, form the basis of this thesis and will be dealt with in more depth later. Even in the first ySR experiments, Garwin et al. [57], Friedmann and Telegdi [57], and Swanson [58]  observed a decrease in the initial amplitude of the muon signal which was  attributable to the formation of muonium. However, muonium was not observed directly until Hughes et al. [60] observed it in gaseous argon.  The first observation of muonium in a solid was  reported by Myasishcheva [68] in crystalline and fused quartz, in solid carbon dioxide and in ice. Muonium was not observed in a liquid until Percival et al. [76] observed it in water; the surprisingly late date was due to earlier workers having failed to remove all the dissolved oxygen.  1 MUON AND MUONTUM FRACTIONS After stopping in a given medium, the muons can be found in either a diamagnetic environment or a paramagnetic environment  For present purposes, all the muons in diamagnetic  environments, such as bare M or M u O H , will be treated as precessing at the same frequency, and +  the corresponding fraction of the stopped muons denoted as P ^ .  Those muons in paramagnetic  environments will be separated into two groups: muonium atoms in which the hyperfine interaction is comparable to that in vacuum, and is nearly isotropic; and muonium-like species, such as muonic radicals [Roduner 84] and anomalous muonium [Brewer 73, Patterson 84] in which the 8  hyperfine interaction is far reduced from that in vacuum and possibly highly anisotropic. corresponding fractions are denoted P  M u  and P ^ j * . u  The  These three fractions are calculated from  the measured initial asymmetries relative to a standard sample, typically aluminum, as in Figure 2, in which 100% of the muons are in a diamagnetic environment  Since none of the work of this  thesis observed any M u states, it will now be dropped from this discussion. From inspection, the asymmetry in Figure 2 is about 0.30, the actual value coming from a fitting procedure which will be described later.  The fractions P ~ and P , can be calculated from: x  D  8  Mu  Muonium itself is, of course, a radical.  9  P  D  =  D — A max A  2A, Mu A max  P. Mu ~  where A ^ , for example Figure 4, and A  M u >  (1.3)  for example Figure 3, are the measured initial  amplitudes of the corresponding uSR signals and &  mSiX  - for example Figure 2, is the measured  initial amplitude in the standard sample. The factor of two present in the definition of P ^  u  takes  into account the loss of one-half the signal due to the formation of anti-parallel muonium. The phrase "missing fraction" will occur repeatedly in this thesis, as it does in M S R literature. The missing fraction P  L for "lost", is that part of the initial asymmetry unaccounted  L >  for in the muon and muonium signals observed in a particular sample, 1- P - P Mu D ] 1  r  (1.4)  r  4. MSR IN SOLIDS The hydrogen atom, and therefore the muonium atom, represents the simplest atomic defect in a solid, as it has no electronic core.  This means that theoretical calculations of its  behaviour may be done more rigorously, possibly leading to the interpretation of more complicated defects.  Hydrogen is also a very common, and commercially important impurity in solids.  The  study of the muon can lead to a better understanding of the role of hydrogen as it can provide information otherwise unobtainable. Of particular importance is the role of hydrogen in the semiconductors silicon and germanium, where it is not known what effects, if any, the large amounts of hydrogen present has on the properties of these two materials.  While no hydrogen EPR signal has ever been  substantiated in them, muonium and muonium-like states are observed and have been studied extensively, as shown by many of the papers in the Proceedings of the Yamada Conference VII, also known as MSR83 [Yamazaki 84]. The study of muons in solids is pursued for two reasons: what they can reveal about the materials themselves, and what they can reveal about the role of hydrogen in these materials. One example of the former is the measurement of Knight shifts [Schenck 81], and of the latter, muonium in silicon [Patterson 84].  The measurement of the frequency and relaxation can be used  Corrected i  Spectrin  i tn  <=>  tn  11 to identify the static and dynamic environments of the muon in the lattice.  In the study of  diffusion, the muon is important as it has a mass between the ultimate light interstitial, the positron, and the proton.  For example, there is now evidence that quantum diffusion of the muon takes  place in low-temperature copper [Kadano 84].  Muonium can be considered in much the same  light as the muon, where the frequencies and relaxation rate provide information about the static and dynamic interaction of muonium with the lattice. Using M S R to study solids instead of the more conventional techniques of EPR and N M R has both advantages and disadvantages.  Some of the advantages are:  1.  Short lived species, with lifetimes down to 10 ns, can be observed;  2.  In principle, the fate of all muons can be determined;  3.  The large mass ratio between the proton and the muon allows isotopic effects to be measured more sensitively than with any other probe;  4.  The muon diffuses much more readily than the proton;  5.  Species are observed in M S R which are unseen by standard EPR;  6.  Relaxation can be observed in the absence of applied magnetic fields;  7.  Muons can be injected into and studied in any material, with no concerns about solubility;  8.  The muon is an infinitely dilute probe, as there is only one muon at a time, hence there are no probe-probe interactions.  Some of the disadvantages are: 1.  The short lifetime leads to lower resolution;  2.  The muon may not have reached its equilibrium state in the time available before decay;  3.  For chemical reactions, there is no opportunity to study the reaction products with techniques such as infrared spectroscopy or via classical chemistry;  4.  It requires a source of muons.  This last disadvantage is becoming less of a barrier with time as the field of M S R has grown much more quickly in the 1970's since the introduction of machines such as those at the "meson  12 factories" at T R I U M F ' , SIN , and L A M P F 10  11  each of which increased the number of muons  available by a factor of about 100 over that obtainable from the earlier machines.  There are also  other machines and more are becoming available all the time, some with even higher intensities. The technology of producing and transporting muons has also been improving, resulting in still further increases in availability. When this present work was begun in 1978, muons had been studied in a variety of solid targets, but most of the results were for muons in diamagnetic environments; there were still relatively few solid systems in which muonium had been directly observed and could be studied, consisting of fused and crystalline quartz, solid C 0 , ice, silicon and germanium [Brewer 75]. 2  It  was not clear why muonium formed in some targets but not in others, and many of the results were in samples of unknown purity.  There was a clear need for more systems in which to study the  behaviour of muonium, more understanding of why muonium forms, and more detailed studies of muonium in some selected systems.  Hence part of the motivation of the work of this thesis was to  do a survey of what were felt to be likely candidates for muonium formation.  Zero or low  transverse fields were used as the spectra are much simpler to interpret than those from the longitudinal fields used by Ivanter et al. [72], for example.  Higher transverse fields resulted in the  signals being lost due to the time resolution of the available apparatus.  As uSR generally uses  much lower applied fields than EPR and N M R , the Hamiltonians describing these properties had to be solved for their low field limits rather than their high field limits. The materials used in this present work were chosen largely because: 1.  they were not ferromagnetic nor were they expected to have significant paramagnetism;  2.  hydrogen atoms had been observed in similar materials [Pake 73];  3.  they were largely free of nuclear moments;  4.  they were available in reasonable purity;  5.  they were available as single crystals;  'rr/'-t/niversity A/eson .Facility originally, although by the time it was actually built there were four universities involved. °Schweizerisches Aistitut fur Nuklearforschung Los Alamos Meson .Physics .Facility l  11  13  6.  the samples had to be large enough in area and volume to stop a significant number of muons and  7.  the materials were non-metals. T h e absence of magnetism and nuclear moments was necessary as it was known from  previous work [Brewer 75] that muonium would be unobservable in zero or low transverse fields in such systems.  That is not to say that muonium might not be present but it would be depolarized  and hence unobservable, resulting in a missing fraction.  Impurities may trap and depolarize  muonium, possibly too quickly for the muonium signal to be observed.  It was expected that  muonium might be trapped at crystal boundaries and depolarized, and if anisotropics were present as first observed by Brewer et al. [73] in silicon, a polycrystalline material would broaden the signal considerably.  The area and volume requirements have declined significantly over the time this  work was in progress; the earliest experiments used a cube about four cm on a side, whereas later a 0.5 mm slice, 7 mm by 7 mm, was sufficient  Metals are not suitable as muonium cannot form due  to the "sea" of conduction electrons present screening the muon. One system, quartz, was studied in some detail, as it proved to exhibit properties previously unseen in uSR.  Another major motivation for the study of quartz was that it was well  characterized by the E P R of hydrogen in quartz [Isoya 83, Perlson 74].  This allowed a comparison  of the behaviour of muonium and hydrogen atoms, which had not previously been possible in any solid system.  Although muonium is clearly just like a very light isotope o f hydrogen (Table I), and  is known to behave this way in gases [Garner 78] and liquids [Stadlbauer 84], at the time this thesis was begun, there was no direct evidence for this in solids. structurally analogous to a - q u a r t z ,  Since hexagonal germanium dioxide is  it was studied as well, in the expectation that a site for  muonium could also be established, thus predicting the site for the hydrogen atom.  A s there is  little hydrogen in hexagonal germanium dioxide [Goodrum 72], the hydrogen E P R experiment might be too difficult, so /uSR could contribute data otherwise unobtainable. In addition to the insulators, muonium was also studied in diamond, a semiconductor. M S R studies had been carried out in its structural analogues silicon and germanium, but only one  14  experiment [Swanson 58] had been carried out in diamond. Since the start of this present work, the situation with respect to the state of knowledge of muonium in solids has changed dramatically due to a recent clever technical innovation at SIN [Kiefl 84].  By applying transverse fields as high as 12 kilogauss in an apparatus with high time  resolution, it is now possible to observe muonium in materials with appreciable nuclear magnetic moments, such as potassium chloride.  Some of their results will be discussed in Chapter V.  It  should be noted, however, that this new technique is not a substitute for, but rather it is complementary to, the zero and low transverse field work of this thesis.  15  B. C H A P T E R II: E X P E R I M E N T A L A P P A R A T U S  1. C Y C L O T R O N B E A M L I N E S A N D D C S E P A R A T O R The experiments described in this thesis were carried out at T R I U M F , a cyclotron on the U B C campus.  T R I U M F is a sector focussed machine, accelerating negative hydrogen ions to an  energy of 500 M e V .  Two electrons are then stripped off, leaving a positive hydrogen ion, or  proton, which travels down a (primary) beam line at about 75% of the speed of light  These  protons strike a "production target" of some material, generally beryllium or graphite, and undergo a variety of interactions, some of which produce positive pions.  Several nuclear reactions are  possible, for example:  ]p+?Be  These pions then decay to muons.  •?ir  +  +  1  8Be  (II.1)  The muons are transported down a (secondary) beam line to an  "experimental" target where they are stopped. For present intents and purposes a beam line is made of magnetic elements, dipole and quadrupole magnets, with the analogous optical elements being prisms and focussing lenses, respectively.  A s a prism changes the direction of impinging light, depending on wavelength, so  does a dipole change the direction of a charged particle beam, depending on the ratio of momentum to charge.  A quadrupole magnet behaves like a focussing lens except that a single  quadrupole only focusses e.g. horizontal.  in one  direction, e.g. vertical, while it defocusses in the  other,  Hence quadrupole doublets or triplets are used, one focussing in each direction,  resulting in a net focussing effect  Momentum and charge selection is performed by choosing only  those particles which come through a dipole in a particular direction. One of the difficulties with charged particles is that only the is selected.  ratio of  momentum to charge  This means that a particle with a charge of + 2 and a momentum of 2p will be bent in  the same direction as a particle with a charge of +1 and a momentum of p.  A much more serious  problem is that particles of different masses can have the same momentum/charge ratio.  Since  16  some of the other interactions which the proton may undergo in the production target generate copious quantities of positrons over a wide range of energies, this leads to the positive muon beam being contaminated with positrons:  In fact, a ratio of positrons to muons of 40:1  observed — rather more like a positron beam contaminated with muons.  has been  Since the experimenter is  generally interested in detecting positrons from muons decaying at rest this contamination can be quite a problem. One way of solving this problem is to use a Wien filter [Jackson 62], a device which uses crossed electric and magnetic fields.  The forces on a given particle are qy_xB for the magnetic  field, and q E for the electric field, where: q is the charge on the particle; v_ is the velocity of the particle; B is the magnetic field; and * E is the electric field. The fields are perpendicular, and the values of the fields are chosen such that only a particle o f a particular velocity is transmitted.  Since momentum has already been selected, only a particular  type of particle will be transmitted through this system. "D.C.  Separator"  contamination.  1 2  This device is frequently referred to as a  and in principle is capable of providing a muon beam 100%  free  of  W h e n operated at a sufficiently high voltage and magnetic field it will also rotate  the spin o f the muon relative to its momentum.  This separation can be done using crossed, static fields in a D . C . Separator, or with just a single oscillating field in an " R . F . Separator". This latter device relies on the fact that the cyclotron output has a time structure, putting out protons in a burst every 44 nanoseconds. A s a result, the pions, muons and positrons also have this feature. Since, for a given momentum, these three have different velocities, their time of flight to a particular point will be different If the R . F . Separator oscillates at the same frequency as the cyclotron and a suitable phase is chosen, the unwanted particles can be steered away. Since its location is fixed, it usually works for only a particular momentum or perhaps several particular momenta, whereas the D . C . Separator has a continuous range. 1 2  17  2. S U R F A C E M U O N BEAMS A N D C O N V E N T I O N A L BEAMS There are two different types of muon beams used in nSR, "conventional" and "surface". Conventional muons have historical precedence , but surface muons are somewhat easier to 13  describe, and were used for the bulk of the work described in this thesis. In the center-of-mass frame the pion decays to produce a 4.1 MeV muon and a muon neutrino as in eq. 1.2.  Since the neutrino always has its spin anu-parallel to its momentum  direction and since the spin of the pion is zero, the muon spin is also anti-parallel to its momentum direction, yielding 100% longitudinally polarized muons . 14  If this pion is at rest in the lab frame,  the result is a 4.1 MeV muon in the lab frame as well - a rather simple case of relativistic kinematics.  In the center-of-mass frame, this decay is isotropic.  Surface muons arise when a  low energy pion produced from proton bombardment in a reaction like (eq. II.l) stops at or near the sjufjace of the production target momentum of 29.8 MeV/c . 15  These muons come off with an energy of 4.1 MeV, yielding a  Some then travel down a beamline through the various magnetic  and electric fields present, and through a thin window to stop in the experimental target their low energy, they stop very-easily, their total range being of the order of 140 mg/cm  2 16  Due to .  This  corresponds to roughly one metre of air at a pressure of one atmosphere or to 0.5 mm of quartz. The spin of these muons is pointing anti-parallel to their momentum direction in the center-of-mass pion frame, which here corresponds to the lab frame. Conventional muons arise from pions which were of higher energy and escaped the production target completely to travel down a beamline where, due to their short lifetime, many of In fact, as an historical footnote, it could be said that Friedman and Telegdi's experiment in 1956 [57] used surface muons, and since Garwin used cloud muons, a case could be made for surface muons having historical precedence. The neutrino is said to have negative helicity, where helicity is defined as: 13  14  where o is the spin and p is the momentum direction, and is a good quantum number for neutrinos. Frequently other longitudinally polarized particles, such as muons, are spoken of as having helicity, but this is not strictly true. A typical nuclear science unit: 1 MeV/c = 5.3x10" kg-m/sec. Another typical nuclear science unit: the range in centimetres of a particle is naturally dependent on the density of the material, and in fact, the range in different materials scales fairly well with just the density, so that R(nuclear science) = R(length units)/p (density) is quite useful. 15  16  22  18 them decay in flight  In the lab frame, due to relativistic kinematics, those pions which decay  "forward" (in the centre-of-mass frame), that is, down the beamline, produce muons with very different energies from those resulting from pions decaying "backward" (in the centre-of-mass frame), or up the beamline.  However, if the initial pion momentum was high enough, those from  forward decays and backward decays will travel down the beamline in the lab frame, but two different types of conventional muon result forward and backward muons, with very different momenta. The spins of forward and backward muons point in opposite directions; the spin of a forward muon, like the surface muon, is anti-parallel to its momentum direction, and the spin of a backward muon is parallel to its momentum direction. There are, of course, muons with momenta between these extremes, but since their polarization is also intermediate, they are of lower utility. A typical momentum for a forward muon might be 170 MeV/c (kinetic energy = 94 MeV), for a backward muon 80 MeV/c (27 MeV). The corresponding ranges would be about 34 g/cm and 2  5 g/cm . 2  Clearly these ranges are far higher than those for surface muons, and so they are quite  inappropriate for targets such as low pressure gases or thin solids. In order to transmit surface muons, the beamline must have all its magnets set for the same momentum.  Conventional muons, however, require that the momentum selected be different in  different parts of the beamline in order to admit pions of one momentum and release muons of a different momentum.  With reference to Figure 5, the muons from pions which decay in the center  or decay section are "captured" by the quadrupoles and the second dipole, B2, is set so that the forward muons go down leg B while the backward muons go down leg A.  3. D A T A - C O L L E C T I O N A N D ANALYSIS Figure 6 shows a schematic diagram of a typical M S R apparatus.  The collimator is  typically about 2.5 cm thick lead or brass and is used to define the "spot size" of the beam. practice most of the muons coming out of the beampipe stop in the collimator.  In  The remaining  muons pass through and "fire" the M] and M counters, typically plastic scintillator, and most of 2  20  P.  M, M, collimator  Figure 6 Schematic diagram of a typical uSR  apparatus.  21 these muons stop in the target counter.  Some may miss the target or pass through it and strike the M  A "good" muon is thus a signal from M i and M and no signal from M . 2  3  3  It is then said  that Mj and M are in coincidence and M is in anti-coincidence, which may be written in Boolean 2  3  notation M ! - M - M . 2  eventually decay.  3  The muons which stop in the target precess in the magnetic field  and  17  Some of the resulting positrons pass through positron counters Pj and P . 2  In  the jargon of the trade, Pj and P form a positron telescope. A good event is thus M]-M -M -PrP 2  2  3  2>  with a delay allowed between M r M - M and P r P , but not otherwise. 2  3  2  The measurement of this delay is the key to doing the experiment, and the observed delay is recorded for each event . 18  Tne maximum delay allowed is controlled by a "data gate" which  opens when a muon is detected and closes at some pre-determined time later, typically eight microseconds or some 3.6 muon lifetimes. assumed it will never be detected.  If a positron has not been detected by then, it is  A l l of the work of this thesis used the time-differential M S R  technique [Brewer 75, Garner 79], which is the most commonly used, although time-integral techniques are extremely useful in some applications [Schenck 81, Heffner 84].  This discussion  applies principally to the time-differential technique. One of the things which may happen is .'that a second muon may enter during the time the data gate is open. This is known as "pile-up". If a positron is detected, it may not be possible to distinguish which muon it came from, so no events are accepted until the current data gate, which is re-started for every muon detected, has expired. Only after that will a new muon-positron pair be considered a good event  Consequently, if the incident muon rate is too high there will be no  good events at all. There is an optimum rate which is given by [Garner 79]: N=  JL  (II.2)  As drawn, the magnetic field is perpendicular to the page, so the muons precess in the plane of the page. Some time elapses between when the muon passes through the counter and when it stops in the target This is not usually important in condensed matter as it is very nearly the same for each muon. 17  18  where N is the number of incident muons per second and T is the length of the data gate in seconds.  A typical T is eight microseconds, which yields an optimum incident muon rate of 62,500  per second.  Further discussion of this can be found in Appendix 1.  The data gate is selected according to the type of experiment being done.  For example, if  the signal is long-lived, a longer gate might be used. The incident muon rate can be adjusted by changing the size of the collimator (see Figure 6).  The times of decay are recorded with a certain  coarseness, known as time per bin, which is also varied according to the type of experiment being done.  The data consists of a number of bins, perhaps 2000 per telescope, containing the number of  events for each time interval.  The data from each telescope then constitutes a time histogram such  as in Figure 1. After the data has been collected, it is fitted to an equation like:  N( t ) = N e  "t/r  0  where  N (t)  (1+^[A e i  i c o s ( u . t + 0j ) ] ) + B  is the number of events recorded at time t  (II.3)  from the muon stop, and the  parameters are defined as: N o is the normalization; r the muon lifetime; A ^ the initial amplitude of the ith signal; X ^ the relaxation rate of the ith signal; CJ^ the frequency of the ith signal; 0^ the phase of the ith signal; and Bg the time-independent background. In general, there may be more than one signal present  The fitting is via a chi-squared  minimization, where:  (II.4)  where:  x  2  is the chi - squared;  23  N .  is the Number of events observed in bin i;  N.  is the Number of events expected in bin /;  1 0  16  o ^ is the uncertainty in the number of events in bin /. Bin i corresponds to time t = i At, where A t is the time per bin. Since Poisson statistics are applicable, a ^ =v -N. _ )  (N.  X = £T 2  The ideal value for x  2  so  1 0 1 0  M  2  (IL5)  1 6  N,  1 0  when the equation describes the data perfectly and Poisson statistics apply,  is  *idea. = " -  < >  f  IL6  where n is the number of bins fit and f is the number of parameters used.  There is, however, a  distribution around this value due to statistical fluctuation. More importantly, if the model from which the fitting equation is derived does not describe the data correctly, the calculated minimum value of the x  2  may be much higher than the ideal value. In this case, the calculated value will  depend on the number of events collected, growing larger with higher statistics, while the ideal value does not change with the number of events collected. The minimizations done in this work used the general non-linear least squares fit program MINUIT [James 71] running mainly on a V A X 11/780 located at TRIUMF. provides best-fit values for the parameters of interest and their associated errors. M S R data, all the parameters are floating, but with limits.  This program When fitting  That is, they are allowed to vary freely  and independently between some upper and lower bounds. It is also possible to "fix" one or more parameters, if their values are known a priori. Some analysis was also done using Fast Fourier Transforms (FFT) [Brewer 82], which are particularly useful for data where several frequencies are expected.  One of the problems with  using an F F T is that in M S R data, the errors in the data are larger for longer times.  Thus, the  24  signal to noise ratio in the FFT deteriorates if longer times are transformed.  If an apodization  function is used which handles this correctly , the signal is unavoidably broadened. 20  problem with an FFT is that of "ringing".  Another  Since only a finite length of time is being transformed,  every signal will have a number of side-lobes [de Haseth 82]. small peak close to a large one.  19  This makes it difficult to find a  Apodization can also be used to suppress this ringing, and there  are apodization functions which minimize the broadening for a given amount of suppression [Norton 76].  Unfortunately, none of these functions corresponds to the form given above for  "correct" statistical weighting. Since a search for weak or closely adjacent signals formed a large part of the work of this thesis, an attempt to circumvent some of these problems was made by developing another form of analysis, one which uses a chi-squared niinimization, but finds the minimum as a function of frequency . 2:  That is, holding the frequency fixed at some value, the minimum in the chi-squared  with respect to the other parameters is found. By "stepping through" a range of frequencies, a table of these parameters and the chi-squared is generated as a function of frequency. some of the appearance of a Fourier transform, but strictly speaking it is not . 22  statistical weighting procedure is followed, and there are also no side-lobes.  This has  The correct  Another advantage  of WSFT is that if some parameters are known a priori, for example the phase, these can be fixed at the appropriate values, improving the signal-to-noise ratio.  Samples of WSFT and FFT are  shown in Figure 7, which is taken from the study of hexagonal G e 0  2  which will be discussed  below. The apparent signal-to-noise is obviously much higher for WSFT than FFT. eight clearly visible peaks in the WSFT spectrum, some of which are lost in the FFT.  There are WSFT can,  however, produce harmonics of the real frequencies, and its resolution tends to be poor.  The  Apodization is simply the multiplication of the signal by a function which varies in time. The approximate function is e [Brewer 82], which assumes that the background B is negligible, an assumption which is increasingly invalid at longer times. An earlier form of this technique was developed independently by Nagamiya et al. [75] to look for weak uSR. signals from negative muons. For a range from, say, zero to fifty megahertz in steps of 0.05 megahertz, some one thousand chi-squared minimizations must be done. This takes rather a long time, perhaps six hours of C P U time on a V A X 11/780, hence it has been dubbed World's Slowest Fourier Transform or WSFT. 19  / 2 t  20  21  22  25  250  V ki  1  i  .  -i  r  r  200  c p tr » 150  5 •3 K)0  10  in en o  50 r-  j  8.E-6  11  J  i_.  L  JLL  h  t> *  o 0, 4.E-6  2.E-6 t  Figure 7 Comparison o f W S F T (top) and F F T . Note the enhancement o f the signal-to-noise ratio of the W S F T spectrum over the F F T spectrum. F o r convenience, the chisquared was subtracted from the maximum observed value.  26 harmonics may be intrinsic to the method, while it may be possible to improve the resolution.  4. E A G L E : A S T A N D A R D S U R F A C E M U O N APPARATUS Most of the experimental work in this thesis was carried out in an apparatus known as "Eagle", the name being an allusion to both a fancied resemblance to the first lunar lander, and to its American origin at the University of California at Berkeley.  It is discussed in [Clawson 82].  Figure 8 shows the interior arrangement of this apparatus. The vacuum chamber is mounted between a pair of large coils in a modified Helmholtz configuration.  A cryostat capable of providing temperatures from 3 to 300K can be placed in the  center of this apparatus.  At present the apparatus is suitable only for surface muons as there is a  single thin muon counter". The eight positron counters are used to define four telescopes, which correspond to four histograms for the data. The collimator and muon counter are between the two counters of the Backward telescope so that a positron emitted by a muon stopping in the collimator or muon counter cannot trigger an event in that telescope — it is very unlikely that a positron could go through one counter and scatter backwards to go through the other.  Further, any positron from  the collimator or muon counter would have to pass through and fire counter B1 in order to trigger any of the other telescopes.  Hence, if the counter B1 is put in anti-coincidence with the other  telescopes, no events of this sort should be observed. One problem with the apparatus is that the vacuum chamber has nine holes for the light pipes to the counters, two holes for the muon degrader support", one hole for the collimator, one hole for a removable sample support rod, four large ports, and approximately eight kilograms of graphite, all of which lead to considerable vacuum problems. Originally the collimator was fixed in place and the entire vacuum chamber - including that of the cryostat - had to be vented in order to change i t  This operation could take as long as  "The very different ranges and energy deposition rates for surface and conventional muons necessitate very different thicknesses of counters, so generally a counter suitable for one type is not at all suitable for the other. Surface muons require a thin counter and all the muons which hit the target stop, so a veto counter is not used. Even with surface muons, it is possible to have a target too thin to stop all the muons, so thin plastic film may be inserted in front of the target 24  27  Figure 8 Interior of Eagle. The D (defining) counter is the muon counter and L . R. F . B (for Left, Right, Backward, Forward) are the positron telescopes.  12 hours until the cryostat was operational again, an unacceptable delay given the premium piaced on beam time (there is fierce, although rarely bloody, conflict for the limited amount of time available).  Needless to say, the collimator was not changed very often under these circumstances.  In fact once. The variation in the area of samples studied in this present work, from 0.5 cm to 5.0 cm , 2  made a more easily variable collimator was clearly desirable.  2  Accordingly, a brass bar with four  circular threaded apertures was mounted on a rod which was put through a vacuum seal in the bottom plate (see Figure 9).  The collimator to be used was chosen by simply varying the height of  the pipe. The apertures are threaded for two reasons: 1.  Smaller collimators could be installed ahead of time, perhaps even non-circular ones  2.  Threaded holes are believed to cause less scattering of the muons than unthreaded ones, or at  25  :  least those muons which would have been scattered are stopped in the threads instead.  5. S O M E C O N S I D E R A T I O N S I N A P P A R A T U S  DESIGN  In designing and choosing apparatuses for this thesis study, there were several criteria that had to be considered.  This discussion will be largely restricted to the use of surface muons.  These criteria are: 1.  the muon must be detected;  2.  the muon must stop in the target; and  3.  its decay positron must be detected. 1) The muon is detected by a counter, the active part of which is plastic scintillator,  typically either 0.013 cm or 0.026 cm thick . 26  allow easy detection of the muon.  The energy deposited must provide enough light to  Since the collimator is immediately ahead of, and somewhat  smaller than, the active area of the counter, no muons should be able to get into the target without This is of limited sharply defined spot dimensions. "Plastic scintillator is or 10 mil being the 25  utility as scattering from the counter, windows, etc., will blur any The tendency is towards a nearly Gaussian shape in two manufactured to round numbers in Imperial units, 5 mil (0.005") example here.  Figure 9 Sketch of the variable collimator. The pipe passes through a vacuum seal and the handle at the bottom allows the operator to adjust the height to choose the size o f the collimator. A friction collar is then tightened to hold it in place.  30 being detected.  The only way a muon might be missed is if two muons come within the recovery  time of the counter, which is about one hundred nanoseconds. 2) Since surface muons have a short range, it is quite easy to stop them in condensed matter targets.  Instead, the problem is to ensure that they do not stop before they get to the target  This necessitates minimizing the total thickness of the beam windows, counter, cryostat windows, heat shield, and magnetic shielding. For example in one run it was observed that 36% of the muons stopped in the 0.005 cm thick (40 mg/cm ) magnetic shielding before entering the target; 2  however, close to one hundred percent of the muons passed through the 0.0025 cm (20 mg/cm ) 2  magnetic shielding which was substituted . Strictly speaking, one should not add the thicknesses 27  of different materials in such a simple-minded way, but it was adequate for present purposes. The range can be determined more accurately from the Bethe-Bloch equation [Evans 55]. Another problem, particularly with surface muons, is multiple scattering by the beamline windows, etc., a pronounced effect due to the low energy of these muons. For example, changing the D.C. Separator windows from 0.0025 cm to 0.0008 cm thickness increases the luminosity (muons/sec/cm ) of the incident beam by a factor of about three. 2  The material in front of the  counter reduces the number of muons entering the apparatus, but the counter itself and the material between it and the target scatter the muons, so that some of them miss the target and stop in the target holder, the walls of the cryostat, or elsewhere. scatters muons through a mean half-angle of 6 ° .  For example, a 0.030 cm thick counter  Since the counter is 5 cm far away from the  target, the beam spot is at least 1.0 cm from scattering alone. By contrast, a 0.013 cm thick counter scatters them only 2.5°, for a minimum beam spot of 0.4 cm.  For small targets then, a thinner  counter is much preferred . 28  The final consideration is that the muons follow a curved path in the applied magnetic field.  This is particularly bad for surface muons which have an energy of typically 3.1 MeV (25.9  The dilemma was that if insufficient shielding were used, the magnetic field would not be low enough and the experiment would fail no matter how many muons stopped in the target The necessity for the use of zero field will be explained later. If the counter is too thin, insufficient light is produced and counter efficiency declines, resulting in a reduced signal and an increased background. 27  28  31  M e V / c ) after passing through the thin counter. calculated from H p =  The radius of curvature of such a muon can be  86 K G a u s s - c m [Trower 65].  For example, this means that in an applied  field of 480 Gauss, the radius of curvature (p) would be 176 cm, and if the distance from the thin counter to the target were 10 cm, the beam would be laterally displaced by about 5.7 mm. target is only 1 cm in diameter, this can be a serious problem. transverse to the muon momentum.  If the  This applies only when the field is  As mentioned above, a D . C . Separator may also be able to  rotate the muon spin relative to its momentum.  If the spin is rotated 90 degrees this results in  transversely polarized muons, which will precess in an applied magnetic field which is parallel to the muon momentum direction.  This allows the use of surface muons in higher fields than might  otherwise be feasible. 3) There are two considerations in detecting the decay positrons, the number detected and the asymmetry they show.  A s a consequence of the weak interaction, positrons emitted with  different energies also give different asymmetries.  This behaviour is given by the simplified  equations [Gamer 79, Sachs 75]:  d R U ^ l  *l[(3-2x)-P<1-2*>cos*]  =  = #-[  1 + aPcosfl]  (11.1)  where:  n = x  a -  2  ( 3 - 2 x )  l ^ ^ i  (118)  (H9)  and where: R is the probability of a positron being emitted at an energy x in a direction 8 from the muon spin direction; n is the number of positrons emitted at a certain energy; a is the asymmetry of those positrons emitted at a certain energy; P is the muon polarization at the time of decay; and  32  x is the energy/maximum energy of 52.8 MeV. The equations for n and a are plotted in Figure 10. The major point here is to note that the higher energy positrons also reflect a higher asymmetry. The quantity of interest is A N 2  29  , where A is the asymmetry, and N the number of  positrons detected of all energies. A =  ffi  d  d  x  N=/ndx  X  (11.10)  The A. defined earlier in eq. II.3 as the amplitude of the rth signal is closely related to this A, with only some small corrections for counter geometries and efficiencies being necessary. the standard sample ^  In practice,  is used as an empirical measure for this correction, simultaneously  m 3 i X  taking into account the fact that the muons may not be 100% polarized.  The quality of signal is  actually proportional to Aj/N, but usually the square is taken for convenience; also, the time it takes to do the experiment is proportional to A N , so it is a useful quantity. 2  This dependence can be  simply justified by noting that A N is the signal and since Poisson statistics apply, the error in the number of counts is just the square root of that number, therefore AN i— signal/noise = ^ = A j / N  (11.11)  Now consider the effect of rejecting all positrons below some certain energy E . From Figure 10 it can be seen that this would, of course, decrease the number of events detected, N , but would also increase A , the observed signal amplitude. possible to calculate A N as a function of E 2  2  A N A M -  By combining equations (eq. II.5) it is  c  (f-l)(3f  3  2  +f +f - H )  2  g ( j. 3 _ f Z_f _ 1 )  where f= the lower energy cutoff/52.8 MeV, the maximum decay energy.  (11.12)  This function is plotted  in Figure 11. The maximum is at f= 0.73 (38.8) MeV, and is 1.7 times that observed with no energy cut This does not necessarily apply when comparing different signals as a relaxing signal should not be treated in the same fashion. 29  Figure 10 Number n (top) and asymmetry a o f the emitted positrons as a function of energy, x = energy/52.8 M e V .  LO  r-  LO  LO CM  LO  r-  LO  in  CNJ  Figure 11 Theoretical plot of the square of the asymmetry times the number observed integrated from a lower limit to the maximum allowed energy, plotted as a function of the lower limit. The points are from a Monte Carlo simulation showing the effect of scattering. Note the decrease in A'N.  35  In practice, the simplest way to make such a cut is to place some "absorber" or "degrader" in the positron telescope as in Figure 8.  By choosing a suitable thickness it is possible, in principle,  to "range out" the lower energy, undesirable positrons.  Due to restrictions of geometry and cost,  the material of choice is carbon in the form of sheet graphite.  It can be calculated from range  tables [Berger 64, 66] that 18.4 grams/cm of carbon is required to eliminate positrons below 38.8 2  MeV.  This is all very well, but one is not interested in the total distance a positron travels, which  is what range tables give, but how far in a particular direction the positron travels.  Due to  scattering, this distance wili be less than the expected range so the absorber thickness must be reduced correspondingly, but more importantly the range is blurred by the random nature of the scattering process.  Thus, a particular lower-energy positron might escape an absorber which a  particular higher-energy one did not  Also, the higher-energy positrons are more likely to  undergo energy loss by Bremsstrahlung which can be a catastrophic process, whereas the range tables use only an average value for Bremsstrahlung effects. A Monte Carlo calculation was done some years ago by Leiss et al. [57] for exactly this interaction.  They calculated the probability of electrons escaping a thick carbon absorber in a  particular direction. Some of the resulting effects on A N are plotted in Figure 11. Comparing 2  the points to the line, it can be seen that the maximum corrected gain in A N of 1.25 is much less 2  than the maximum expected gain of 1.7. the 18 expected.  The maximum occurs at about 13 grams/cm rather than 2  A measurement was done with 7 grams/cm of carbon, all that would fit in the 2  available space in Eagle, which was found to increase the asymmetry by 30%, whereas the increase calculated from Leiss et al. is only about 25%.  This unexpectedly high value is probably due to the  sample thickness, as discussed below. For higher Z (atomic number) materials these "blurring" effects are even worse, and lower Z materials are not cheaply available in sufficiently dense forms, so graphite remains the optimum choice. As an historical note, perhaps the chief reason for the use of degrader in the positron telescope with surface muons was to prevent the scattered beam positrons from triggering the  36 positron telescope, thereby increasing the background. Since the beam momentum, 29.8 MeV/c, is less than the momentum of the desired decay positrons, they can be partially ranged out along with those, decay positrons reflecting lower asymmetry. With the advent of clean, separated beams, this practice should no longer be necessary. One of Tire side-effects of the dependence of the asymmetry on energy is the dependence, in some circumstances, of the signal amplitude on the sample thickness.  This makes it more  difficult to compare samples of different materials as almost invariably they are different sizes, densities and chemical composition . 30  Consequently, the most reliable telescope in Eagle for  asymmetry measurements is the Backward one (see Figure 8), as in this case the decay positron need only penetrate the small thickness that the muon itself travelled, thus reducing the unwanted sample dependence. Occasionally it is necessary to have zero applied magnetic field.  Since the ambient field at  T R I U M F is one to three Gauss, in virtually any direction, to obtain a reasonably good (less than 50 mG) zero field it is necessary to cancel this field.  There are two ways to do this: use magnetic  shielding material; or use a triaxial coil arrangement A disadvantage of the former is that the limited range of the surface muons means that little or no mass must be in the path of the beam.  In the restricted space available it was not  possible to use a cylinder with an open end toward the beam, so a thin magnetic shield had to be placed in the path of the muon. The commercially available shielding obtained had a thickness of .005 cm, or 40 mg/cm , which as mentioned in Section I.I.5, stopped some of the muon beam. 2  Thus it was necessary to make the shield thinner.  It was found that the material, a nickel alloy,  could be etched by a hot acid solution, reducing the thickness to 0.0025 cm, or 20 mg/cm , which 2  was thin enough. Another disadvantage of magnetic shielding is that while the shield is in place, no external magnetic field can be applied. It is sometimes desirable to turn on a small magnetic field to normalize asymmetries, for example.  Chemical composition may be important as the influence of Bremsstrahlung is strongly Z-dependent 30  The disadvantage of the coil method is that this ambient field changes with time.  It is  then necessary to monitor the ambient field, a task rendered difficult in Eagle as the sample is in a cryostat, an unsuitable environment for a temperature sensitive probe, and a large vacuum tank surrounds that, so any monitoring probe must be quite far from the sample.  In any future  apparatus, such as the Omni apparatus currently in the final stages of construction at TRIUMF, a provision should be made for a triaxial, temperature controlled, magnetic field probe.  38 C. CHAPTER III: M U O N I U M H Y P E R F I N E INTERACTIONS  1. S P I N - H A M I L T O N I A N I: ISOTROPIC H Y P E R F I N E INTERACTION The first studies of the hyperfine states of muonium were done by Hughes et al. [60] who were interested in measuring the muonium hyperfine frequency v in vacuum as a test of Quantum 0  Electrodynamics (QED).  They did this by using high-pressure argon gas which both stopped  their conventional muon beam and served as the medium for muonium formation.  By measuring  the frequency as a function of pressure, they were able to extrapolate to get a value of v  0  for  muonium in a vacuum. In such studies muonium is, of course, in an isotropic environment The isotropic hyperfine interaction is by far the most commonly observed in M S R , and will not be pursued much here, as it is well treated elsewhere [Brewer 75, Fleming 79].  The  Hamiltonian describing it is: H = g UnS-B • e B— — where  g  g  g n I-B + M ix— —  A I-S — — 0  (III.l) '  v  , g ^ are the electron and muon g-factors, respectively;  M g , M are the Bohr (or electron) and muon magnetons; y  S , I_ are the electron and muon spin operators; B is the applied field; and A = h^o is the contact hyperfine interaction. 0  In vacuum, p = 4463 M H z as given in Table I; the corresponding value for the hydrogen atom is 0  1420 M H z or 507 G , a unit more natural to EPR.  In this work all energies will be given in  frequency units. The solution of this Hamiltonian yields the familiar Breit-Rabi diagram shown in Figure 12.  When the electron spin is initially parallel to the muon spin, muonium is in a triplet state, but  when the electron spin is initially anti-parallel to the muon spin, the state is a mixture of singlet and triplet in which the muon spin oscillates at the hyperfine frequency. transitions, only four are allowed in a transverse  Of six possible  magnetic field by angular momentum  ©  o o o  K)  oo o  cs  oo o  o  o o o o o oo o o o r '  -  i  CM  l  ro  o o  l  Figure 12 Breit-Rabi diagram for the isotropic hyperfine interaction. The arrows indicate the allowed transitions. The numbers are the labels used in the text  40 considerations.  These four are indicated on Figure 12. If the field is in the z direction, the  observed muon precession can be written [Percival 82]: P  = P + iP xy  x y 1 , icJi t , 1 , i w t = ^-c e + ^-s e 2  2  , (III.2)  23 23  2  TTT  where the two precession frequencies CJ.J = 2itv „ are: u„='u.-Q  (III.3)  w  (III.4)  = u. + S2 , with  23  a>  1 ± cj ), where =-^(o)  +  (III.5)  fi  M  o; = ggMgB,  the Larmor precession frequency of an electron  C J ^ = g^u^B,  the Larmor precession frequency of a muon  g  = ^{[ul + 0)1] ~  (CJ,  +  2  CQ  4w 2 1 C =n 2  1 +  s = 1- c 2  ^  23  - G) } 0  )  2  (III.6)  0  (III.7)  0),  [  W +  2  +  U  2  0  ]  1  /  2  (III.8)  2  In the low-field limit, s =c =-|- and w is the hyperfine frequency, A =w /27r. 2  2  0  0  0  The well-known results are given in reports by Garner [79], and by Brewer [75], and will not be derived in detail here. 1.  For present purposes, the important features to note are:  two components oscillate at or near the hyperfine frequency p and are difficult to observe 0  directly, although this has been achieved with high time resolution [Holzschuh 81]; 2.  in zero applied field, the two low frequency transitions degenerate to one, with a frequency of zero, meaning the signal is constant in time, barring relaxation;  3.  for low fields (<20 Gauss)  31  this degeneracy still holds within the typical experimental  frequency resolution and a single characteristic low-frequency muonium signal can be observed, as exhibited in Figure 3, whose frequency is given by 1.39 MHz/Gauss times the  The field regimes assume the hyperfine constant is similar to that for muonium in vacuum.  31  41 applied field ; 32  4.  for medium fields (20<B<150 Gauss) the degeneracy breaks down and two signals are observed.  The average value of their frequencies, v. = 2ira>_, is again 1.39 M H z / G times the  applied field; and their splitting, 2Q, varies as the square of the applied field. This effect is usually referred to as the quadratic Zeeman effect  The measurement of fl is the standard  method of determining cj , since a  -4S2  ~  (o)  +  2  (CJ,  + CJ 4fi  12  3 )  2  + CJ  2  2  + 2w )  3  2  2  2. S P I N - H A M I L T O N I A N II: A X I A L L Y S Y M M E T R I C H Y P E R F I N E I N T E R A C T I O N In the case of an axially symmetric anisotropy, the Hamiltonian can be conveniently expressed as [Brewer 79b]:  H  =  9e B^*^" M  9  M M-'M  +  A  ° - ' -  +  A  A  I  z  S  (  z  I  I  U  0  )  where the symbols are defined as before, with AA the deviation from isotropy in the z (symmetry) direction, and the subindex z indicating the z component, The solutions of this Hamiltonian are given in detail by Patterson et al. [78], and the approximate solutions for small B are given in papers by Beder [80], Baryshevsky et al. [78] and Percival et al. [82].  The resulting energy level diagram is shown in Figure 13.  dependent on the angle between the applied magnetic field and the z-axis.  The details are  The observed muon  polarization is still given by equation III.2, but now with the precession frequencies CJ  12  = a). - S2 + ^ 5 C J  (111.11)  The occurrence of this low-frequency signal is the signature for isotropic muonium, and much of the work of this thesis consisted of searching for i t 32  -4000  1  0  1  100  1 200  1  300  1  1  400  1  L_ 500  600  Field (gauss)  700  L 800  900  1000  43  o> = u>- + S2 - 2 8CJ , where  (111.12)  23  = ^AA(3cos 0 - 1)  6CJ/2TT  (111.13)  2  with 0 and c defined as before, but with D replaced by u>^, where 0  = A + 4-AAsin26  W./2TT  (111.14)  0  This expression is only valid for u> _ > -~bu>. For zero magnetic field a>_ = 0, (HI-15) and the intensity I of this transition in the direction parallel to the initial spin of the muon is (111.16)  I - 75-sin 0 2  where <j> is the angle between the z axis and the initial spin of the muon.  The intensity in a  direction perpendicular to the initial spin of the muon is cos\//cos0 where $ is the angle between the direction of observation and the initial spin. The two main consequences of this Hamiltonian are: 1.  the observation of a low-frequency muonium oscillation in zero applied magnetic field, A A , whose amplitude is orientation dependent, but whose frequency is orientation independent; and  2.  a splitting of the two transitions in low applied fields 8CJ which is field independent in the low field regime, but strongly angle dependent, This splitting is in addition to the usual splitting seen in medium fields for isotropic muonium.  ^ 2 3  -  (*>u =  Thus  2 0 - 6a>  (111.17)  This allows the determination of the sign of A A , which otherwise would be difficult to determine.  3. SPIN - H A M I L T O N I A N III: C O M P L E T E L Y ANISOTROPIC H Y P E R F I N E INTERACTION Arguably, this section should be first, as the axially anisotropic and isotropic Hamiltonians are but special cases.  However, the order chosen reflects the increasing complexity as well as  historical sequence of observation.  The appropriate Hamiltonian is:  44  H  where A is now a tensor.  g  =  E  M  B - * - ~  M  9  M  M - ' -  - ' = ' -  +  (  I  I  U  8  )  Although more physically realistic, g is not taken as a tensor here, since g  none of the systems studied in this thesis had a sufficiently anisotropic g for this complication to g  be necessary at the low fields used. If the axes are chosen such that A is diagonal, H = q Mr,S-B-g  B— —  ^LL  M I'B + A I S + A I S + A I S ix— — xx x x yy y y zz z :  (111.19)  from which it is clear that equation III.l can be derived by setting A = A = A = A , and ^ xx yy zz equation III.2 by setting A = A , with A = A , and A A = A - A ^ ° xx yy xx zz xx 0  c  0  J  The eigenvalues and eigenvectors of this Hamiltonian can be found in Appendix 2, and the energy level diagram is plotted in Figure 14. The low field behaviour depends strongly on the orientation of the applied field relative to the principal axes of the tensor. The major new features are: 1.  With zero applied field, three frequencies are observable, given by i ( A - A ), i ( A  -  £ xx yy ^ yy A ) and i ( A - A ) as shown in Section A2.1, eq. A2.15. The signs are not observables in zv 2 zz xx V  y  M  6  zeio field. A useful feature is that these frequencies have a "sum rule"; that is, the sum of two of them equals the third.  Since the signs are not observed, the two lowest observed  frequencies sum to the highest observed.  This can be a helpful clue . 33  The three resulting  frequencies will be called a "triple". 2.  In zero field, the intensity pattern of these frequencies is strongly angle dependent, but the frequencies themselves are orientation independent  Orientation is with respect to the initial  muon polarization, as for the cylindrically symmetric case with zero applied field. This orientation dependence can, in principle, be used to determine the orientation of the principal axes of the tensor with respect to the crystal axes. When the first experiment on low-temperature a-quartz was done, the anisotropic Hamiltonian was completely unexpected, and this feature was one of the first things noticed. 33  re  3000  —  7  ? cr & P>  & re  N  8  I re  •§. n  a  a cr a>  re  3. 3 re i—•  a R  3 o  a.  8  100  200  300  400  500  600  Field (gauss)  700  800  900  1000  46 3.  In an applied field, the spectrum becomes much more complex, with as many as nine frequencies.  This case is not considered further here, since no experimental work was done  on it, but more theory may be found in Appendix 2. As shown in Section A2.3.C, taking the initial muon polarization as the z' direction in the lab frame, and if the site x axis is parallel to z', the only frequency observed will be v  = ^(A  -  A ), and the intensity of this frequency in the direction of the initial muon polarization will be 1/2. That is, the amplitude of the signal at this frequency will be 0.5, just that seen for muonium in a transverse field for an isotropic hyperfine interaction. and z' is 7 ,  , then the intensity of die v  On the other hand, if the angle between x  frequency will be icos 7 ,  . This is just the same  2  as it was in the case of the axially symmetric hyperfine interaction, the difference being that there are now three such frequencies, each varying as the square of the direction cosine between the site axes and the z' axis.  That is, l(v v  Jj-cos7 , 2  .  2  ) = Jcos7 , xy  2  y  'zz  , l(v v  ,  2  ) = icos 7 xv  2  , and l(v  'yz  v  ) =  yz  y  These intensities are all in the direction of the initial muon polarization. If there is  more than one equivalent site present, it can be much more complicated, as shown in Section A2.3.b.  It should be noted that the expression for intensities given in the paper by Brewer et al.  [81] is incorrect. 4. S U P E R H Y P E R F I N E S P I N - H A M I L T O N I A N  The superhyperfine interaction occurs when another nucleus interacts with the electron of the muonium atom. The Hamiltonian describing the superhyperfine interaction is: H = q M ^ S - B - q u I-B + I-A-S + J - A - S e B— — M M — — = — — = — N  3  -  (111.20) v  '  N where A and A  are both tensors and J is the spin operator of the nucleus.  For a spin one-half  nucleus in zero magnetic field, neglecting the states corresponding to singlet muonium, this Hamiltonian gives rise to a six-by-six matrix, as shown in Appendix 3.  This treatment is  completely independent of that carried out by Beck et al.[75] for the specific case of muonium in alkali halides.  47  The intensity behaviour is now tensorial in nature, this tensor being different for each of the three observable frequencies.  This is far more complex than for the case with just the  anisotropic hyperfine, where there was one tensor which described the intensities of all three frequencies simultaneously.  It should be noted that the sum rule for triples holds for the  superhyperfme interaction just as it did for the anisotropic hyperfine. A numerical example should illustrate, taking the observed values for the anisotropic hyperfine and the superhyperfme determined by Isoya et al. [83], for hydrogen in quartz. The reason for this choice will become apparent later.  At 80 K, the principal values of the hyperfine  tensor are, 51.9483, 51.7505 and 51.6981 mT, with the corresponding values of the " S i superhyperfme tensor being 0.068, 0.134 and 0.130.  The tensors will be assumed to be coaxial.  Substituting these into eq. A3.2, Appendix 3, gives the resulting transition energies: 0.0290, 0.1263 and 0.1553 mT. As noted, the intensity of each frequency is given by a tensor, defined by: Ijj=-u«I'U  (111.21)  where u is the direction of the initial muon spin, also assumed to be the direction of observation. If the polarization makes an angle 114° from the x axis, 24° from the y axis, and 90° from the z axis of the tensor, the corresponding intensities are: I = 0.0064,1 = 0.1456 and I =0.0535, where I2  an intensity of 0.25 corresponds to a fully allowed transition.  13  23  48 D. CHAPTER IV: A BRIEF SURVEY OF M U O N H J M F O R M A T I O N A N D R E L A X A T I O N Muonium's ionization energy of 13.6 eV is higher than that of almost all materials. has considerable kinetic energy available during its slowing down.  It also  Perhaps the question should  not be "Why does muonium form?", but rather, "Why doesn't it always form?" It must be noted that the term "relaxation" of the muonium (muon) signal is used to mean three very different processes: 1.  muonium (muon) may react chemically to form a diamagnetic (paramagnetic) compound in which the muon precesses at a vastly different frequency.  The muonium (muon) is thus  effectively removed from the ensemble in an irreversible manner. 2.  muonium (muon) may be dephased. that is, different members of the ensemble may precess at slightly different frequencies, or the frequencies may change as a function of time, but there is no chemical change.  The result is, though, a loss of signal amplitude. In principle this is a  reversible process. 3.  if the muonium electron undergoes spin exchange with some other electron, the signal will be lost, irreversibly.  The first mechanism can easily be distinguished from the last two.  In zero applied field, any  muonium that reacts to form a diamagnetic compound causes no loss of polarization, whereas spin exchange or local magnetic moments will cause a loss of polarization.  The application of a strong  longitudinal field, however, can be used to quench this loss of polarization [Brewer 75]. Although the work in this thesis is concerned almost exclusively with solids, most of the existing work on muonium formation has been done in gases and liquids, so this work will be briefly reviewed. Ultimately it is hoped that an understanding of muonium formation in gases, the simplest systems, will help lead to an understanding of muonium formation in liquids and hence in solids.  49 1. OASES Of the three phases of matter, gases are usually the easiest to understand. This is true for muonium studies as well.  In the gas phase the muon goes through a number of charge exchange  cycles where it captures and loses an electron as it slows from its initial energy of approximately 3 MeV  34  down to thermal energies.  The amount of muonium formed is thus the result of the energy  dependence of the relative cross-sections for these charge exchange processes [Fleming 82]. There are other considerations, including the possibility of hot atom reactions [Arseneau 84], but they are less important  In all gases thus far examined only two, helium and neon, do not show  large muonium fractions. The relaxation of muonium in gases, either pure gases or with added impurities, is due to the three processes mentioned above.  Chemical and spin exchange reactions of muonium are the  subject of considerable current interest [Senba 84, Fleming 82], with the lower limit on the rates of reaction being due to line broadening from inhomogeneous magnetic fields. This is much more troublesome for gas phase than for condensed phase work as the muons have a broad stopping distribution, since: 1.  they have a high initial energy and hence travel a long distance (approximately fifty centimetres) before stopping in the gas;  2.  the energy loss process has a statistical fluctuation; and  3.  the beam has an initial spread in energies due to the momentum acceptance of the channel.  The magnetic field produced by any finite set of coils varies spatially. This means that muons stopping in different regions will be in slightly different magnetic fields, hence precess at slightly different frequencies and the signal is dephased. This effect can be reduced by the use of large coils to produce a homogeneous field over a sufficiently large volume.  For example, the coils  presently used by the gas chemistry group at T R I U M F have a diameter of 150 cm.  The  diamagnetic muon signal in gases will also relax due to magnetic field inhomogeneities, although about one hundred times more slowly than muonium.  Other possible sources of muon relaxation  "Although surface muons have an energy of 4.1 MeV, the windows and counter reduce the energy to about 3 MeV.  50 include thermal muonium formation or processes involving n  molecular ions [Fleming 83].  2. L I Q U I D S  In liquids, the situation is much more complicated due to the collective processes.  It is  sufficiently complex to have been the subject of considerable debate in the literature [Percival 84, Walker 81] as well as in conferences, group meetings and the hallways of TRIUMF.  Two main  differences are observed for muonium and diamagnetic muon fractions between gases and liquids: a much lower muonium fraction and a much higher diamagnetic fraction are observed in liquids than in most gases; and there is frequently a missing fraction, of the order of 20% in water for example, whereas such a missing fraction is almost never observed in the gas phase [ Arseneau 84]. The earliest model attempting to explain muonium formation was the "hot atom" model which observed that "At a kinetic energy of 200 eV, practically all the muons are in the muonium state.", [Brewer 75] from which they slow down; but due to their high energy, they may participate in epithermal reactions which would be forbidden for thermal muonium. The observable fractions of muonium and diamagnetic muon would then depend on charge exchange processes, as in the gas phase, followed by possible chemical reactions before the muon or muonium had finished thermalizing. Such processes are relatively unimportant in the gas phase [Arseneau 84] but the large diamagnetic fractions seen in liquids relative to those in vapors suggests that they may dominate in condensed media. In 1978 a new model, the "spur" model, was developed [Percival 78] in analogy with the model proposed by Mogensen [74] for positronium formation. of a positron and an electron in a bound state.  Positronium is an atom conisting  This model is based on the premise that in its  slowing-down process the muon creates spurs consisting of electrons and radicals typical of the particular solvent  It would then be possible to form muonium by a thermalized muon picking up  51 one of these electrons, in contrast to the cyclic charge exchange process which takes place at higher energies. One of the big continuing mysteries in uSR is the origin of the missing fraction,  . This  is zero to a fair degree of precision in many materials, particularly in metals and gases at modest pressures [Arseneau 84], but in many cases it is not; indeed in some extreme cases neither muon nor muonium signals are observed at all! The reasons for this loss of polarization have been the subject of considerable speculation.  One advantage of the spur model was that it gives a plausible  explanation in terms of the radiolytic species formed in the spur, such as solvated electrons, which would be expected to depolarize muonium by spin exchange or chemical reaction [Percival 84]. No model which does not involve muonium or muonium-like states has yet been proposed as a mechanism for the creation of the missing fraction. For the moment the debate of hot-atom versus spur in muonium formation remains unsettled.  As an explanation for the missing fraction the spur model is very attractive, although  difficult to test conclusively. Thehot-atom model is even harder to test Relaxation processes in liquids are similar to those in gases except the stopping distribution is much narrower, so that field inhomogeneities are less important  The most important difference  between gases and liquids comes from the possibility that radiation-induced spurs contribute significantly to the muonium formation process in liquids.  3. S Q L L D S  a. Muptiium Formation The mechanism of formation of muonium in solids is even cloudier than in liquids. Again the muon has enough ionization energy to forcibly seize an electron from most materials.  Although muonium is not observed in metals due to the conduction electrons  present which screen the muon charge, it is observed in some insulators and some semiconductors; but often, as in liquids, there is a large missing fraction. works for liquids, then perhaps it holds for solids as well.  If the spur model  One might speculate that if the  52 missing fraction is due to muonium spin exchanging with a spur electron, then materials which allowed the electrons to diffuse rapidly away might, therefore have low missing fractions, a - quartz has large channels which might permit such diffusion, and it has no missing fraction. Ice, on the other hand, also has large channels, but it does have a missing fraction [Percival 85].  A more sophisticated approach along these lines might be profitable.  Typical results  prior to the present work may be found in a paper by Brewer et al. [75] and are mentioned in Chapter I.  In some cases, this missing fraction is easily explained as being due to muonium  being relaxed by nearby nuclear spins.  For example in CaF the superhyperfme couplings  between the muon's electron and the local  2  19  F spin 1/2 nuclei would produce such a large  number of different frequencies that the result would look like a fast relaxing muonium signal, or even, since the first few nanoseconds are not readily observed, no muonium would be seen at all, and a missing fraction would result  The superhyperfme interaction has long been the  bane of workers in this field as the majority of nuclei have magnetic moments, but recendy at SIN the use of strong magnetic fields to quench the superhyperfme interaction has begun to conquer this problem [Kiefl 85].  Missing fractions are still observed in their work, however,  so that mysteries yet remain. That the systems in which muonium is observed are not at thermodynamic equilibrium is clear from the fact that muon and muonium signals are present simultaneously. It is very unlikely that an equilibrium would give observable quantities of both muonium and diamagnetic muon signals in so many different solids.  Further, it is known [Isoya 83] that  hydrogen atoms are not observed in crystalline quartz, for example, when it is at equilibrium at room temperature, whereas muonium can be observed, so muonium cannot be at equilibrium. b. Muonium Relaxation Relaxation is a much more complex topic in solids than in gases or liquids.  In  addition to the sources of relaxation already mentioned, there are two more: nearby nuclei with magnetic moments; and the anisotropic hyperfine interaction, or AHI.  53 The former is the most common, these nuclei acting as local magnetic dipoles, producing what is termed random local magnetic fields, or R L M F . 35  This can be described in  a simple -minded way as the magnetic moments adding to the applied field at the muon or muonium, causing the muon spin to precess at a different frequency than another muon near a different magnetic moment  The superhyperfine interaction also occurs in the presence of  nuclear moments, but rather than acting like random dipoles, a very strong coupling takes place.  A great deal of work has been done on the importance of R L M F on the relaxation of  muons in metals [Gygax 84, Celio 84].  The relaxation, X in eq.II.3, if caused by R L M F is  related to the square root of the second moment [Abragam 61] <At> > 2  « M  2  adapted from a paper by Schenck [76], where M and f is the concentration of nuclei with spins.  N  n  ^jpf  (IV.l)  is the nuclear moment, J is the nuclear spin This assumes just one type of nucleus with  spins, but it may be easily generalized. The AH1, discussed in Section III.3, is more unusual; indeed, the work presented herein is the first muonium study to have correctly identified relaxation due to AHI.  That is,  earlier studies in crystalline quartz [Myasishcheva 68] and ice [Percival 78] had relaxations caused by AHI, but this was not recognized at the time. Relaxation of the muonium signal by AHI can be caused either by: 1.  Static muonium atoms in a large number of sites with different values for the hyperfine parameters; or  2.  Muonium which is not static, but moves among different sites which are not perfectly equivalent In the first case, the relaxation envelope will indicate something of the width and  strength of the anisotropics present in progress [Turner 85].  The detailed theory of this is not yet known, but work is  Qualitatively, if members of an ensemble are in sites with different  principal values of the hyperfine tensor, then they will precess at different frequencies, R L M F may be present, of course, in liquids and gases, but in these systems the rapid motions of the molecules average these fields to zero. 35  54 resulting in a dephasing.  This is true in both zero and transverse fields. The larger the  variation among the sites, the greater the width of the distribution, then the faster the relaxation will be. In the second case, if the muon spin evolves in time in one of the n sites for a time A t , then "hops" to one of the other n - 1 sites, the intensities of the three frequencies will be different  This, repeated for many muons, with random A t ' s and possibly more than one  hop, results in a dephasing of the signal.  Such a dephasing will occur even if all sites have the  same principal hyperfine values, so long as there is a difference in relative orientation, whereas such sites would not cause dephasing for sialic muonium. This is because it is essential to consider the amplitudes of the different frequencies not as they are observed, which is necessarily an ensemble average, but rather on a single site basis since the relative orientations of the sites must be considered. From looking at the transformation given in eq. A2.63, it is clear that, generally, the transition probabilities are different in two different sites.  Therefore,  if a muonium atom spends some time in one site, then "hops" to another, the result will be a dephasing of the muonium signal. Typically, the model of diffusion used herein will assume that the time spent between sites is much less than the time spent in each site. To see this more clearly, consider a simple example in which there are two sites which have the same principal values for their axially symmetric hyperfine tensors, but one of which has its z axis parallel to the initial muon polarization, and the other its z axis perpendicular to the initial muon polarization. In each of these sites two of the three possible frequencies collapse to one with frequency v, and the remaining frequency has a value of zero. In site 1, the zero frequency transition has all the intensity, as shown by eq. A2.62, and in site 2 all the intensity goes to the frequency v.  The spin of a muon of muonium in site 1 will not change  with time whereas in site 2 it will oscillate with a frequency v.  So long as there is no hopping,  the spectrum will consist of two sharp lines, one at 0, and one at v.  However, if the muonium  atom hops from one site to the other, that spin which was constant will now oscillate, and that spin which oscillated will become constant at whatever value it had when it hopped . A n  55 ensemble, each member undergoing hopping at random times, will undergo dephasing. This is true even in zero magnetic field, and for fields of a few gauss, the relaxation may be even faster.  Qualitatively, one can see that the shortest possible time for the signal amplitude to  fall to 1/e of its initial value is of the order of one-quarter of the period of the fastest oscillation.  Hence the relaxation rate  X = 4f  (IV.2)  where f is the highest frequency signal observed in the static case. If the rate of hopping increases sufficiently, the average of the tensors will be observed.  For example, if a three-fold axis is present, the result will be just the average of  the three tensors after they are rotated to the common crystal axis system. The final result is that an anisotropic hyperfine Hamiltonian becomes an axially symmetric one, with the unique axis parallel to the symmetry axis.  56 R CHAPTER V: REST JITS A N D DISCUSSION  The particular samples were chosen by the criteria given in Section 1.4.  Briefly, samples  had to be insulators with few or no nuclear moments and single crystals were preferred. This excludes most available materials, all the hydrides, halides, and alkali metal compounds, but still leaves a significant number from which to choose.  The original motivation for the study of quartz  was its natural helicity. The motivation changed to a study of the effects of anisotropic hyperfine interactions.  Germanium dioxide was studied due to its close analogy to quartz, and beryl and  zircon were chosen because they a T e silicates like quartz.  Rutile and magnesium oxide were  chosen because they are simple oxides, and hydrogen atoms have been observed in CaO [Henderson 72], which is very similar to MgO, and strontium titanate was chosen for its ferroelectric properties.  Diamond does not fit into this group, but rather it was chosen due to the  intense interest in M S R studies of semiconductors [Yamazaki 84].  The results on all these  materials are summarized in Table II,. with their sources and impurities given in Table III.  All the  samples used were synthetic, unless otherwise noted.  LOJIARIZ Quartz has a three-fold screw axis, the t or z axis, and there are two sets of three two-fold axes, a and B.  In this work the usual coordinate system [Nuttall 81] will be used, in  which the x axis is parallel to an a axis, and the y axis is not an axis of symmetry, but is orthogonal to x and z. a.  Asymmetry of Life  Quartz in its fused and crystalline forms forms the backbone of this thesis.  Quartz  has been studied with M S R for rather a long time; in fact it was one of the first solids in which muonium was directly observed [Myasishcheva 68]. because it occurs in right-  and left-  In the present work, it was initially used  handed crystalline forms.  This early study was  inspired by the two observations that life on earth uses essentially only L-amino acids, and that spatial parity is violated by the weak interaction.  Soon after the discovery of parity  57 Table D Summary of results in a variety of solids Sample  T(K)  muon fraction  (%)  muon  muonium  muonium  missing  relaxation  fraction (%)  relaxation  fraction (%)  (us )  (us )  1  a-quartz Fused S i 0  2  Tetragonal G e 0  2  Fused G e O j  4-295  35±5  6515  017  4-295  30±5  7015  017  8-295  70±5  <5  3015  295  76±3  .071.02  1014  2.811.8  1415  60  36±2  .021.01  1616  8.612.4  4816  914  4.411.0  6 Hexagonal G e 0  MgO  2  295  61±3  90  48±3  6  10013  32  52±3  4-295  60-70  BejAljSuO,,  295  ZrSiO« Diamond  Ti0  2  <2  3914  .161.02  1316  .351.06  3917  .091.01  30110  5.711.3  3 0 1 10  40±2  80  SrTiO,  1  <5  4816  16±5  <10  84115  295  40±5  <5  60110  295  48±2  .051.01  3314  1.8071.042  1914  6  21130  .121.03  2615  7.211.0  53130  Table m  Purities of materials used  58  Material  Manufacturer  Impurities  a -quartz  Sawyer Research Products Inc.  <0.01 ppm Fe, <0.3 ppm M n , 1-5 ppm Na, 5-20 ppm H , 0.2-0.5 ppm A l . 0.5-1.0 ppm L i  Adolf Meller Co.  <50 ppm Fe, Si, Ni, Mo, Cr, <30 ppm Co, B, <10 ppm Cu, Te. Mn, Mg, Be, Bi, 40- 64 ppm A l  SrTiO,  Adolf Meller Co.  <50 ppm Fe, Si, N i , Mo, Cr, <30 ppm Co, B, <10 ppm Cu, Te, M n , Mg, Be, Bi, 40- 64 ppm A l  MgO  Adolf Meller Co.  110-140 ppm Fe, 3-20 ppm M n  tetragonal Ge0  John W. Goodrum  40 ppm Fe, 30 ppm Cr, 80 ppm A l  2  hexagonal GeOj  John W. Goodrum  50 ppm W. 20 ppm Na  TiOj  59 violation, it was suggested [Vester 57] that there might be some connection between these two observations.  A large number of experiments have failed to produce conclusive evidence of  such a connection, but the present work was begun at a time when Garay [74] apparently had found one. The attempts to find a link have taken two routes.  One is to attempt through heavy  bombardment with particles of only one handedness, for example longitudinally polarized electrons, to differentially decompose a racemic mixture [Bonner 76].  The other is to use, say,  the positron, which is created in 0 -decay with a handedness, as a direct probe, looking directly, on an individual basis, at the interactions of the polarized particles with the medium. This discussion will be restricted to the latter approach.  An excellent collection of work on  the subject of the origins of optical activity in nature has been edited by Walker [79]. Garay did an experiment looking at positronium formation. The positron acts much like the muon in that it captures an electron whose spin is either parallel or anti-parallel to its own, forming the e e~ bound state called positronium (Ps). +  The two are called ortho (o-Ps)  and para (p-Ps) positronium, respectively. A positron undergoes mutual annihilation with an electron to produce gamma rays.  If the positronium is unperturbed the annihilation will be  between the bound positron and electron; p - P s having a shorter lifetime (T = 1.2xl0 s) than 10  o-Ps (r =1.4xl0 s) producing two gamma rays, whereas o-Ps produces three gamma rays. _7  In condensed matter, o-Ps usually transforms to p - P s , which then decays, but this still takes longer than the direct decay, and so the initial quantities of p - P s and o-Ps can in principle be measured from the observed lifetimes [Mogensen 74]. relative amounts of p -  Garay attempted to measure the  and o - Ps in single crystals of the L and D forms of some amino  acids, and claimed the ratio of p - to o - Ps was different in the two isomers. Garay suggested an explanation of this difference in terms of a "helical electron gas" model which said that electrons travelling in a handed molecule would tend to have their spins aligned in one way, say parallel to their motion.  In the same molecule, but of the opposite  handedness, the spins would be anti-parallel to their momentum. Thus the electrons would  60 show positive helicity in molecules of one hand and negative helicity in the other.  Since a  positron or a muon will tend to capture electrons whose velocity is parallel to their own , if 36  these electrons have helicity then the ratio of p - to o - Ps in the case of the positron, or triplet to mixed muonium in the case of the muon, normally expected to be one, would be different from one and would be different for the two molecular hands. The muon should be more sensitive than the positron to any such effect because: 1.  its polarization (helicity) can be reversed by changing from surface or forward muons to backward muons, thus providing a check on systematic errors;  2.  its greater mass means its direction at the time of capture will be more likely to be the original direction than for the lighter positron. That is, at lower energies they retain their helicity better than positrons. At the time it was not possible to observe muonium in amino acids or sugars, as they  have many nuclear moments, so quartz was chosen as it was known from previous work that it exhibited a large muonium fraction.  An additional feature of quartz is that its optical  handedness is reversed for the secondary optic axes relative to the principal optic axis. A high statistics experiment was performed to measure the muonium amplitude in a weak transverse magnetic field using forward and backward muons in lefthanded quartz, parallel and perpendicular to the principal optic axis.  and right-  The final result is a  difference in amplitudes of (1.08 ± 0.86)% where the error is, as usual in this thesis work, one standard deviation.  This number gives no particular confidence in saying that it is not  consistent with zero. A more detailed discussion can be found in [Spencer 79]. At about the same time other experimenters [Jean 79] repeated Garay's work with positrons and got a null result to better precision than his experiment  They stated that his  results could be explained by the presence of defects in his samples.  T h e energy difference will be smaller and the time in which they are close is longer. 3  61 b. Asymmetry of Muonium (in Quartz) i.  Room Temperature: Axial Hyperfine As is said to happen so often in science, serendipity took a hand, yielding results  unconnected with the original intent of the experiment.  When the first test runs were  performed, the muonium signal was observed, as in Figure 15, to have a distinct beat frequency of 0.4 MHz, whereas that expected from the quadratic Zeeman effect in the applied field of 6.3 G is only 0.018 MHz.  The initial explanation of this phenomenon  was that it was due to the intrinsic quadrupole moment [Baryshevsky 78, Beder 78]" of the triplet muonium atom interacting with an axially symmetric crystal field gradient [Brewer 79].  This can be described by the cylindrically symmetric anisotropic  Hamiltonian described in Section III.2. The first work of this thesis was to test the predictions of this model as given in Section III.2. One prediction was that the splitting which produces the beat should equal •^(3cos 0 - 1)AA (eq. 111.13), where 6 is the angle between the electric field gradient and 2  the magnetic field. This splitting should be independent of the strength of the applied field up to the region where the quadratic Zeeman term begins to contribute. true can be seen in Figure 16.  That this is  The model also predicted that in zero magnetic field there AA  would be an oscillation at a frequency dependent  -y  whose amplitude would be orientation  This can be seen in Figure 17, which shows two orientations, one in which  the t axis is parallel to the initial muon polarization, which shows no oscillation, and one in which it is perpendicular, which shows an oscillation. With the use of eq. in.16, the symmetry axis of the electric field gradient was found to be parallel to the c axis, the highest symmetry axis of the crystal. It is also possible to use the model to calculate the sign and estimate the amplitude of this splitting [Beder 80].  The sign may be determined by measuring the  "This intrinsic moment arises from a d-orbital admixture in the ground state of the triplet atom.  62  Figure 15 Data showing the distinct beating of the muonium signal in room-temperature quartz.  crystalline  63  3  COS  2  6  -1  Figure 16 Anisotropy of splitting in a - quartz. 6 is the angle between the c - a x i s of the crystal and the applied magnetic field. The round points are taken at 6.3 Gauss and the triangles at 3.3 Gauss. T h e error bars are smaller than the points. T w o triangles are plotted at 3 c o s 0 - l = O as the sign was not detennined. l  64  0.07  E E U3  <  E E  <  -0.07 0.0  1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  Time in lis (248 ns/bin)  Muonium oscillation in zero field room temperature a -quartz. W h e n the C - a x i s is parallel the initial polarization of the muon (top), the oscillation has zero amplitude. However, wher the e - a x i s is perpendicular to the initial spin, an oscillation occurs (bottom).  65 total splitting given by eq. 111.17, as a function of the applied field.  If they are of  opposite sign then as the field is increased the splitting will go to zero and then increase again, whereas if they are of the same sign then the splitting will increase monotonically. From the angular dependence it is clear that the sign of the splitting depends on the angle. For # = 0 °  the splitting, and therefore A A, was found to be negative.  This was  confirmed in a measurement by Brown et al. [Brown 80] as well as Barsov et al. [83] who measured the sign of A A for 8 = 9 0 ° .  However, the sign predicted by the model is for  A A to be positive [Beder 80]. It was then concluded that the oscillation was not due to the intrinsic quadrupole moment, but rather a more complex mixture of the intrinsic moment plus effects induced by the lattice [Weil 81b]. All these results are from room-temperature measurements.  Since there are  large channels along the t axis of a-quartz and given the cylindrical symmetry of the gradients about the c axis, it was thought that the muonium might be diffusing rapidly in that direction. It would then be interesting to see what happens at a lower temperature, where the muonium might diffuse more slowly. ii.  Low Temperature: Completely Anisotropic Hyperfine Accordingly, with zero magnetic field, a sample was cooled and a dramatic change  took place.  Below about 100 K, instead of the single, low-frequency (0.4 MHz) signal  seen at room temperature, three higher frequencies at 1.7, 6.2 and 7.9 M H z were seen. The frequencies were independent of sample orientation, but the amplitudes were strongly angle dependent as shown in Figure 18. The top spectrum shows the intensity pattern for the crystal y axis parallel to the initial muon spin direction, exhibiting the 1.7 and 6.2 M H z signals, and the bottom one shows the spectrum for the c axis parallel to the initial muon spin, yielding the 1.7 and 7.9 MHz signals.  This latter orientation is the one for  which no oscillation could be observed at room temperature. sum to the upper.  The lower two frequencies  More exact values, as given in Table IV, are: 1.692(11), 6.218(12),  7.899(10), and summing the lower two gives 7.910(16), in excellent agreement with the  66  4. 5 E - 4  -i  •  1  1  1  r  1  4. E - 4 3. 5 E - 4 3. E - 4 OJ  Y AXIS PARALLEL TO SPIN  2. 5 E - 4  3 O  2. E - 4  1.5E-4  1.  E-4  5. E - 5  0  /\  .0014  1  1  1  1  1  i  1  1  L A ^ V _  I  1  . 001 2 . 001  OJ  C AXIS PARALLEL TO SPIN  8. E - 4  3 D  c  6. E - 4  I-  _  4. E - 4  L  1  A , 2  1  3  4  5  Frequency  _  L  6  -  7  8  9  (MHz)  Figure 18 Fourier transforms of 6 K data in a - quartz showing the angular dependence or the signal intensities, with the frequencies constant  10  67  Table IV Predicted and Observed Hyperfine and Superhyperfme Transitions for M u o n i u m in a-quartz  2.26(7)  9.14(6)  11.4(1)  Observed M u o n i u m Frequencies ( M H z )  1.692(11)  6.218(12)  7.899(10)  Ratios  1.33(4)  1.47(1)  1.44(1)  Predicted Superhyperfme Frequencies ( M H z )  1.9  8.0  9.8  Predicted Intensities  0.026  0.0031  0.071  Predictions of Muonium Frequencies from Hydrogen E P R Data ( M H z ) [Isoya 83]  The observed hyperfine frequency for muonium is 4505.0 M H z at 30 K , [Holzschuh 81] whereas that predicted from hydrogen E P R data is 4561.8.  Angles of the principal axes for hydrogen and muonium, angles in degrees.  Hydrogen  113.75(19)  23.75(19)  90.00  Muonium s e l l  118(5)  28(5)  90(5)  M u o n i u m set 2  62(5)  152(5)  90(5)  68 observed value.  This looked like a quadrupole moment in a non-cylindrical field  gradient, and is described by the anisotropic Hamiltonian given in Section III.3 where each of the three sites related by the three fold axis has its own hyperfine tensor, with the same principal values.  This Hamiltonian has the same form as that observed for  hydrogen in a - quartz [Isoya 83]. These measurements represented the first time that a completely anisotropic hyperfine interaction was seen for muonium; and it is of some historical interest to point out that the first time an anisotropic spin Hamiltonian was seen for the hydrogen atom, it was also in quartz [Perlson 74]. If the site seen for muonium were exactly that seen for the hydrogen atom, then the frequencies arising from the anisotropies would be as given in Table IV, which also presents the observed frequencies, as well as the isotropic value of v measured by 0  Holzschuh et al. [81].  Note that the observed frequencies from the anisotropies are  about 1.4 times less than those predicted, and the isotropic term is 1.3% less than that expected based on the hydrogen results.  These differences are perhaps due to the  increased zero-point motion of muonium versus hydrogen, a trend which is consistent with Isoya et al.'s [83] observation that deuterium shows a larger shift than hydrogen from the vacuum value and also that the anisotropies are larger for deuterium than for hydrogen.  The value of v for both muonium and the hydrogen atom are higher than 0  the vacuum value; 0.9% for muonium, and 2.2% for hydrogen.  This suggests that the  lattice compresses the atom, as discussed in Isoya et al. [83]. As shown in Section A2.3.C and mentioned above in Section III.3, the observed M S R intensities can be used to deduce the orientation of the principal axes of the site with respect to the crystal.  However, as shown in Section A2.3.d, in zero magnetic field the  presence of the three-fold symmetry axis destroys all information about one of the angles, <f>, that which gives the angle around the t axis with respect to the y axis.  If a  magnetic field is applied, this information should be obtainable, at the cost of a very  69 messy spectrum with as many as nine muonium signals at different frequencies.  The  angular information can in principle be extracted from the splitting of these transitions as a function of the angle between the crystal axes and the field.  Unfortunately, this was  not realized at the time the experiments were done and so it remains as a worthwhile project for future studies.  Nevertheless, the angles of the site axes with respect to the c  axis were determined, and were found to be 118(5)°, 28(5)°, and 90(5)° at 6 K, which are within error of the hydrogen numbers which are 113.75(19)°, 23.75(19)° and 90.00° at 40 K [Isoya 83].  The hydrogen site is shown in Figure 19 and more details of the  hydrogen site can be found in the paper by Isoya [83].  There are two choices of sets of  angles due to the fact that it is the square of the cosine which is observed in zero field. These are given in Table IV.  One set does agree with the observed hydrogen numbers.  This was the first time that muonium was unequivocally shown to behave just like an isotope of hydrogen in the solid phase. iii Superhyperfine interaction For hydrogen in quartz, Isoya et al. [83] confirmed their choice of site by observing the coupling of a nearby " S i nucleus via the superhyperfine interaction.  The  expected angle the axis between the two nearest neighbour silicons makes with the c axis is 112° from the structure [Le Page 80], and the angle found is 116.5° which agrees well enough.  A n attempt was made to observe this interaction for muonium which would  have been the first direct observation with M S R of lines due to a superhyperfine interaction.  A broader discussion is available in a paper by Beck et al. [75].  Only 4.7%  of the Si nuclei are " S i , but there are two equivalent silicons in the site, so 9.4% of the sites should have a " S i and exhibit the superhyperfine coupling.  The calculation done in  Section III.4 for hydrogen can now be applied to muonium where the crystal z axis is parallel to the initial muon polarization, in zero applied field. It is necessary to scale the H atom values both by the ratio of the magnetic moments and also by the factor 1.4 observed for the anisotropic hyperfine frequencies.  The result is three peaks, at 1.9, 8.0,  h Figure 19  2  A  Diagram of the hydrogen site in a quartz, adapted from Isoya [83].  71 and 9.8 MHz, given in Table IV, whose intensities with respect to the major anisotropic hyperfine peak at 7.9 MHz, are 0.026, 0.0031, and 0.071 respectively, taking the natural abundance as well as the intensity tensors into account  These peaks are in addition to  those normally expected from the anisotropic hyperfine interaction in this orientation at 1.7 and 7.9 MHz.  As shown in Section II.5, the quality of a M S R signal is proportional to  A N , which means that about 250 times as many events would be required to see the 2  largest of these three additional signals as to see the one at 7.9 MHz. A series of runs at 6 K were made over a period of about 24 hours and the sum of these, corresponding to 100 million events, is plotted in Figure 20.  Unfortunately, when  a standard fit to the time histogram was attempted, it was found that the clearly visible relaxation envelope of the normal 1.7 and 7.9 M H z signals did not fit to any of the standard functions, Lorentzian, where the relaxation is e  or Gaussian, where the  -o t 2  relaxation is e  2  .  This can be seen by inspection.  Other functions were tried, at  considerable expense of real and C P U time, but none of them proved satisfactory.  In  spite of this, some fits were tried with an additional frequency to see if the x improved 2  noticeably. X  2  The best fit gave an additional frequency at about 8.5 MHz, and changed the  from 4246.6/323 degrees of freedom to 3269.8/320 degrees of freedom.  It was  thought that this improvement might be due to an artifact of fitting a shoulder of the poorly fit main peak, so as a check, the x for a frequency of 7.3 MHz, an equal distance 2  from the main peak, was found to be 4030.2/320.  A rule of thumb [Bevington 69] when  adding a parameter is that unless the change in x  2  exceeds the ratio of the x  number of degrees of freedom, it cannot be said to be significant ratio is 13.1.  2  to the  In the present case the  For three parameters, then, the additional parameters may be significant if  the change in x exceeds 39.3. 2  The observed change is 976.8. By comparison, the fit to  an assumed peak at 7.3 M H z changed the x  2  by only 216.4. This latter change is far less  than that observed for the 8.5 M H z fit but by the crude rule used might possibly be significant  This is a good indication that the 8.5 M H z peak is real, but the fit is still  i—..  to  •! • * — r  CM  O  eo pi CM  CM  .1 CM  o  c 2  oi  c  ~  a. E  oo 6 to  d  CM  *  o d  d  •. ' i m o d  _l o o d  J m o d I  o d  Figure 20 M u o n i u m signal at 6 K . Very high statistics. Note that the relaxation envelope is not Lorentzian nor Gaussian, as the signal relaxes at early times faster than later times.  73 appallingly bad. Fast Fourier Transform techniques were also applied to this high statistics data using various apodizations, and the results are shown in Figures 21 and 22.  The usual  ringing is so intense in both the unapodized and weakly apodized spectra in Figure 21, 38  that it would mask any additional frequencies out to 9 M H z . The strongly apodized spectrum (bottom Figure 22) is so broad that it also masks all frequencies out to 9 MHz. The spectrum with medium apodization (top Figure 22) indicates a peak at about 8.6 MHz.  That it is only a bump on a larger peak and shows up only under carefully chosen  conditions casts doubt upon its reality. A WSFT analysis was also done, as shown in Figure 23, and it, too, shows a peak at about 8.5 MHz.  This WSFT peak looks much more convincing than the FFT peak, but  it is still a small bump next to a huge bump, though, so this, too, is not conclusive. Moreover, the WSFT analysis is not truly independent of the normal fitting procedure, so it cannot be said that this 8.5 MHz signal has been seen by three independent methods of data analysis. That this peak is at lower frequency than the 9.8 MHz predicted by the hydrogen data (Table IV) is surprising, as the coupling is very sensitive (1/r distance from the muonium electron to the " S i nucleus.  3  [Isoya 83]) to the  Since the muonium has a larger  zero-point motion than hydrogen, < l / r > should increase, shifting the peak to higher 3  frequency.  No such peak is observed.  Isoya et al. [83] report that the coupling for  deuterium is about 10% smaller than for protium, as expected, counter to the result for muonium. The 8.5 M H z peak may be: 1.  an artifact of the fitting procedure due to the obscure nature of the relaxation;  2.  from a next-nearest or next-next nearest neighbour interaction, as seen by Isoya et al. [83] for hydrogen;  The apodization functions used are those recommended by Norton and Beer [76] as optimal in their suppression of sidebands relative to the broadening of the central maximum.  38  74  . 0 0 1 7 5 .0015 . 0 0 1 2 5 h .001  h  o  °- 7 . 5 E 4 \5.E-4  \-  2 . 5 E 4 \0 .001  8.E-4  r-  6.E-4 u it o a.  4.E-4  h  2.E-4  r-  Figure 21 Fourier Transforms of the high statistics 6 K a -quartz data. The unapodized spectrum (top) shows a great deal of ringing, as does the weakly apodized.  75  Figure 22 FourieT Transforms of the high statistics 6 K a -quartz data. The strongly apodized spectrum (bottom) is very broad, but the medium apodization gives a hint of a bump at 8.6 M H z .  76  1.6E5  0  2.5  7.5  10  12.5  15  17.5  20  22.5  25  9.5  10  Frequency (MHz) 500  400 -  O o 3  ~ o  300  I S  200  CO  cn ao  100  _J 5.5  6.5  I a7  7.5  Frequency  8  8.5  (MHz)  Figure 23 W S F T of the high statistics 6 K a - quartz data. The bottom plot is shown with a greatly expanded scale. For convenience, the chisquared was subtracted from the maximum observed value.  77 3.  from a different type of site from the major peak at 7.9 MHz.  The question would then remain: where is the peak expected from the superhyperfine interaction?  The relaxation seen is probably due to the R L M F from the " S i nuclei. It is  perhaps ironic that the difficulty in seeing the " S i superhyperfine interaction in quartz was due to the presence of the other " S i nuclei. It should be noted that the predicted frequencies and intensities are based solely on the data for hydrogen, and are quite sensitive to the values used. For instance, it may be that the 8.5 MHz peak is not that from the 7.9 MHz transition, as expected, but that from the 6.2 M H z transition, which had been expected to be at 8.0 M H z (Table IV). Although this peak was expected to have a much lower intensity, the larger zero-point motion might alter the interaction enough to allow i t expected increase in the frequency.  This would also then fit with the  Due to this aspect of the problem, any further search  for a peak from the superhyperfine interaction in quartz must also show the peak or peaks at more than one crystal orientation, allowing all three to be visible, thus making an unambiguous assignment possible.  A sample which was isotopically enriched in " S i  would be extremely helpful in any future search.  The maximum amplitude for the  occurrence of a single " S i nearest neighbour per site occurs when the abundance is 50%. Since the residual range of surface muons is about 100 mg/cm , for a 1 cm sample, only 2  2  100 mg of quartz would be required. Since silicon makes up about 50% of the mass, only 25 mg of " S i would be required. To make a sample in which changes in orientation were possible would require more material, but if three such thin slices could be grown in the appropriate orientation and arranged as three perpendicular faces of a cube, still only 75 mg of " S i would be required. Perhaps these could be grown as surfaces on the faces of a cube of unenriched quartz. If a spin rotator, such as that mentioned in Chapter II were available, then the changes of the crystal axis with respect to the muon spin direction could be accomplished by rotating the muon spin, allowing the use of a single thin slice.  78 iv Intrinsic Quadrupole. Moment Although the picture of an intrinsic quadrupole moment in an electric field gradient is physically appealing, a simple calculation by Beder [80] of where this moment is located shows that fully one-half is located quite far from the nucleus, beyond two Bohr radii.  At this radius the electron is certainly interacting strongly with the lattice and  the electric field gradient is changing rapidly.  Although the Hamiltonian is correct, the  physical interpretation is not as clear as for the case of a nuclear quadrupole in an electric field gradient, where the extent of the nuclear charge distribution is small compared to the distances to neighbouring atoms and more fundamentally, the quadrupole moment is not at all influenced by the electric field. A calculation by Weil [81] for hydrogen in quartz suggests that a 2p admixture might be sufficient to account for the observations. However, this is forbidden by the symmetry of the site for a static hydrogen atom.  If,  however, the atom is moving, even with just the expected zero-point motion, this admixture will be allowed [Isoya 83].  Goshen et al. [83] show that a d-orbital  admixture caused by the lattice is a reasonable possibility and would give this result even for a static atom.  More detailed calculations of the hyperfine structure  for  hydrogen-like atoms is given in a paper by Weil [81b], which does not deal with a particular lattice, but rather provides a framework for lattice calculations. v Relaxation in alpha-quartz As the temperature is raised the muonium begins to move, "hopping" from site to site.  There are three equivalent sites related by a 120° rotation around the three-fold  crystal t axis.  As previewed in Section IV.3.b, all three sites have the same three  characteristic frequencies in zero magnetic field, but they give different amplitudes due to the different orientations with respect to the muon polarization at the time the muonium goes to a particular site.  The result is a relaxation, and the observed rate as a function of  temperature is plotted in Figure 24.  As the temperature increases further the hop rate  begins to be comparable to the characteristic frequencies themselves, resulting in a peak in  79  O  o  CM  O  in  o o  o m  -I  o  CO  o  Figure 24 Relaxation rate versus temperature for a -quartz in zero field.  o  80 the relaxation at about 140 K.  From the rough formula, eq. IV.2, X = 4f gives 32 ns'\  and the maximum observed rate is 27+7 MS" , in very good agreement, especially since 1  there is no guarantee that the real maximum was observed. Above this temperature, the muonium hops so rapidly that it begins to average the sites resulting in a motional narrowing which reduces the relaxation with increasing temperature, as seen in Figure 24. The result is that the average of these three sites has cylindrical symmetry with the unique axis the three-fold axis.  This should sound familiar, as it is the observed  room-temperature situation of a slowly relaxing single frequency as in Figure 17. If the tensors corresponding to the three low-temperature sites are averaged to give the axially symmetric tensor, the result gives A A = -1.44 MHz, which has the observed sign, but a different magnitude, from the value observed at room temperature of -0.8 MHz.  Holzschuh et al. [81] calculate A A = -5.4 MHz, again with the observed  sign but yet another magnitude.  That the calculated numbers are different from the  observed one is not particularly surprising as the crystal lattice does change slightly with temperature, and as the observed frequencies are small differences between large numbers; extrapolating from 6 to 295K is almost certain to be incorrect signs agree is probably fortuitous.  That even the  Isoya et al. [83] have observed changes in the  hyperfine parameters for hydrogen, but the EPR spectrum was not reported above 100 K since the hydrogen signal vanishes, at too low a temperature to be useful for extrapolation. Holzschuh et al. [81] report changes in the hyperfine frequency u at or near room 0  temperature, too high a temperature to compare to the muonium low-temperature results.  The juSR results above 100 K and the site determination by EPR together  provide a nice example of the complementarity between M S R and EPR. vi Fused Quartz Some measurements were also done on fused quartz; even though it had been used many times before in M S R [Brewer 75] due to its large muonium fraction, this was the first study in which the relaxation rate as a function of temperature was studied in  81 detail.  The relaxation rate is shown in Figure 25, where it is seen to decrease  monotonically with increasing temperature, in marked contrast to the behaviour seen in crystalline quartz, shown in Figure 24. This fused quartz data was taken in an applied field of 10 G , and a typical plot of the room temperature data is in Figure 3. The model proposed to explain the relaxation is the same as that for crystalline quartz, except that in fused quartz a large number of anisotropic sites exist with different parameters, so each muonium atom precesses in its local environment, each with a different set of frequencies and amplitudes. The result is a relaxation as the phase coherence of the ensemble is lost. Below 50 K, the relaxation rate is no longer temperature dependent, so it is concluded that the muonium is static and so the relaxation rate represents the distribution of frequencies from the anisotropic hyperfine interaction.  At higher temperatures the muonium hops  and averages these sites which leads to motional narrowing, as was the case for the crystalline phase.  In principle it should be possible to calculate the diffusion rate from  the observed relaxation rates; however, the theory is not yet sufficiently advanced for this to be possible. vii Summary of MSR in  Quartz  The amplitudes of muonium in fused and crystalline quartz are within error of each other as shown in Table II. Neither exhibits a missing fraction, which is unusual in the solid state. The relaxation rate in fused and crystalline quartz may still be decreasing with increasing temperature at room temperature.  The levels of impurities are at about  the same levels in the fused and crystalline, so if impurities dominated the relaxation, then they should have the same relaxation rate. due to " S i nuclear moments.  This would also be true if the relaxation were  The observed rates are similar, with the fused quartz  having a slightly higher rate, 0.30(5) MS" compared to 0.19(2) us~ \ perhaps due to the 1  random anisotropics not being fully averaged.  Measuring the relaxation rates at elevated  temperatures might clarify this. A sample of natural quartz did show a faster relaxation, 0.43(3) MS" compared to 0.19(2) M S " , than the synthetic one used in these studies, 1  1  T  T  T  T  "4  ro (  srt)  CM  3 > B i iu o T v e x e i a n  *~  nfl  Figure 25 Relaxation of the muonium signal in fused quartz versus temperature.  83 presumably due to higher concentrations of impurities. a-quartz undergoes a phase transition at about 846 K [Frondel 62], and studying the muonium frequency and relaxation rate as it goes through this transition might be interesting.  It is unfortunate that even at room temperature the muonium is hopping so  fast that the direct hyperfine observations of Holzschuh et al. [81] only reveal the average. There are also a number of other forms of Si0 such as cristobalite, tridymite, keatite, 2  coesite and stishovite [Frondel 62] and from the point of view of studying muonium formation it would be interesting to see if these showed a different probability of muonium formation. They are not generally available as large single crystals, though, so that it might be difficult to see muonium if the completely anisotropic hyperfine interaction is present .  Fused silica can be prepared with various amounts of other  39  materials, such as N a 0 present 2  A series of experiments with various mixtures, as  suggested by A. M . Stoneham [84], might be revealing as to the mechanism of formation of muonium in solids, much as mixtures of liquids have done for the study of the formation of muonium in liquids and to belabour the point mixtures of gases have done for muonium in gases. As shown in Section A2.3.C, an anisotropy in the hyperfine Hamiltonian can lead to a signal amplitude in the direction perpendicular to the initial muon spin even though there is no applied magnetic field present tested i t  This result was surprising to me and so I  The data are plotted in Figure 26, and a signal is clearly visible. For certain  orientations of the crystal, in particular that with the principal axis of symmetry parallel to the initial muon spin, this signal will have zero intensity. Quartz could be considered a laboratory for M S R since it has been used in so many experiments.  It was the one of the first solids in which muonium was observed, the  first material used in a D E M U R (.Double Electron Muon .Resonance) experiment [Brown 80], the first material in which the completely anisotropic hyperfine interaction was 'For a powder in zero magnetic field this interaction would yield three frequencies with the same intensity.  3  84  V) TD  C  o o  CD W O  o (D  E  Figure 26 Muonium oscillation in zero applied magnetic field in a direction perpendicular to the initial muon spin.  85 observed [this thesis] and perhaps the first material in which the superhyperfine interaction was directly observed [this thesis].  This is also in some sense its role in EPR  where many different paramagnetic species have been seen, ranging from the atomic hydrogen mentioned frequently in this work, to germanium, lithium, [Weil 71] and aluminum [Weil 75].  An excellent review is given by Weil [84].  ?.. G E R M ANTI TM DTOXTDF, Germanium dioxide was studied due to its similarities to silicon dioxide. In particular, its hexagonal form is closely analogous to a-quartz [Young 62].  A tetragonal and a fused sample  were also studied. The purities of these samples are given in Table III. The results, however, are somewhat disappointing. As is so annoying! y common, there is a large missing fraction.  This cannot be due to nuclear moments, as only 7.8% of the Ge nuclei  ( Ge) have nuclear spins compared to the 4.7% " S i in quartz. One of the problems with the data 73  is that the hexagonal and fused samples were only about one cm in diameter, and the scattering from the thin counter and windows means the beam spot in the Eagle apparatus is necessarily larger than this.  This reduces the data rate as many muons miss the target, and increases the  background due to stops in other parts of the apparatus . 40  Moreover, the conditions were not  necessarily the same for all the runs, so detailed interpretation of the changes in fractions is not possible.  The exception to the foregoing is that changes from no muonium observed to some  muonium observed is significant No muonium was observed in the tetragonal sample, but a diamagnetic signal of 70 % was seen, leaving a missing fraction of 30 %. yellow-brown, although translucent  This sample was so impure as to be distinctly  These coloured impurities are believed to involve unpaired  electrons [Goodrum 72], so they would be highly depolarizing.  Taking the large missing fraction  Here the term background refers not to the background represented by the B term in equation II.3, counts which are independent of time, but rather to the more insidious one of spurious muon and/or muonium signals arising from muons stopping in materials other than the sample. It should be noted that muonium has never been observed in any of the materials used in the construction of this apparatus, so this background acts principally to reduce the sample muon and muonium signal intensities.  40  g  86 as an indication that muonium was formed and depolarized by these impurities implies that it might be fruitful to grow and study a cleaner crystal, if possible.  Thermal or other bleaching  methods should also be considered. Only a small amount of muonium was observed in the fused sample which otherwise behaves just like the fused Si0  2  sample, in that the relaxation rate decreases with increasing  temperature (Figure 27), although the error bars are much larger.  This relaxation can be  attributed, as in the case of fused silica, to random anisotropic hyperfine interactions. However, the relaxation rate in fused G e 0 is higher at all temperatures than it is in fused Si0 . 2  2  If the  relaxation is due to the random anisotropic hyperfine interaction, then this would indicate that the deviations from anisotropy were greater by a factor of 2-3 than in fused Si0 , which may be 2  related to the frequencies in hexagonal G e 0 , as mentioned below. The hyperfine frequency v „ 2  was measured in a field of 60 G , and found to be 4.0± 0.5 GHz, consistent with the vacuum value, and the value found in quartz. frequencies v  u  and u  2i  The Fourier transform is plotted in Figure 28, where the  are clearly visible, the value'for v coming from these two frequencies as c  substituted in eq. III.9 As was expected, the most interesting result is that for the hexagonal G e 0 , which was 2  polycrystalline, with impurities as in Table III.  Due to its structural analogy to a-quartz [Smith  64], it was expected to show similar behaviour, a large muonium signal in a site with a completely anisotropic hyperfine interaction.  What was observed was a small muonium signal and a more  complex zero-field frequency spectrum than anticipated. The ordinary Fast Fourier Transform (FFT) of the zero-field spectrum at 6 K showed several frequencies as in Figure 29 (top), but due to the low signal-to-noise ratio, the World's Slowest Fourier Transform (WSFT), as discussed in Section II.3 was invented and applied. The resulting spectrum is also shown in Figure 29 (bottom).  The frequencies, amplitudes, and phases  of the peaks are given in Table VI. The most dramatic observation is that among the eight largest peaks there are two "triples". For the completely anisotropic hyperfine interaction one expects, as is observed in low-temperature a-quartz, three frequencies where the two lowest sum to the  87  CO  CD LO CM  I-I—BB-  CD CNJ  • • " CD  LO  (0  S EH CD LO  , •  oo  |  -Q  ^  CM  CD  (_STY)  y  Figure 27 Relaxation of the muonium signal in fused G e 0  5  versus temperature. Compare with Figure 25  I  I  ID  CO  CO  1  1  1  UJ  UJ C M •  F  i  g  U  m  ID  I • 1  UJ COm  I rI 1  UJ • (O  I  L  i UJ •  o f fused G e O . at 60 G , room temperature. Notice the two frequencies, the lower one  being v „ and the upper is v „ . The splitung is TSl  89  1.E-S  8.E-6  h  6.E-6  h  tl o a.  4.E-6 h-  2.E-6  250  t  I  1  1  1  1  r  -i  B  •a  1  r  B  200  t> u m 3  _j  150  u  m o SO C  J 5  10  15  20  25  Frequency  30  35  40  45  50  (MHz)  Figure 29 FFT (top) and WSFT or 6 K high statistics hexagonal G e 0 data. The letters group the eight observed frequencies according to their assignment. See text. 2  90 Table V  Frequencies observed in G e 0  2  Frequency ( M H z )  6.690(13)  16.947(18)  23.705(24)  100* Amplitude  0.27(4)  0.21(4)  0.18(4)  Relaxation Rate ( M S )  0.28(6)  Phase (degrees)  31.(5)  1  B Frequency ( M H z )  9.148(8)  18.974(7)  28.177(20)  lOO'Amplitude  0.42(3)  0.43(3)  0.14(3)  Relaxation Rate  (us ) 1  0.17(3)  Phase (degrees)  21.(2)  Frequency ( M H z )  10.21(4)  100* Amplitude  0.28(5)  Relaxation Rate (*iS"')  0.98(26)  Phase (degrees)  13.(11)  Frequency ( M H z )  31.98(3)  lOO'Amplitude  0.67(6)  Relaxation Rate ( y s )  1.57(20)  Phase (degrees)  15.(5)  1  91 highest. This signature is observed in hexagonal G e 0 , only there are iffio. such triples. The sets 2  of frequencies are labelled A and B, as shown in Figure 29;there arc also two additional peaks labelled C and D. The sum of the two lower type A peaks is 23.64(2) MHz, which is in reasonably good agreement with the frequency of the highest type A peak, 23.71(2) MHz.  The sum of the  two lower type B peaks is 28.12(1) MHz, which is also in reasonably good agreement with the frequency of the highest type B peak, 28.18(2) MHz.  Both of these are slightly outside the error  bars, but that the near agreement could be a coincidence is beyond credence. then that there are two sites with different values for the interaction. intensity for each site, the relative intensity for B : A is 1.5(2): 1.  The inference is  If one takes the total  Whether this ratio reflects a  difference in abundance of the sites, or some sort of thermodynamic equilibrium value is not yet clear.  The frequencies seen are about four times higher than those seen in a-quartz, which is  consistent with the faster relaxation seen for muonium in fused G e 0 versus fused Si0 . 2  2  A least squares fit using MINUIT, as discussed in Section II.3, to all the 8 frequencies independently and simultaneously, gives essentially the same relaxation rate to the six members of the two triples. It is reasonable to suppose that the two triples would have similar relaxation rates, assuming this is due to the " G e nuclear moments, as in the case of " S i nuclear moments in low-temperature a - quartz. The other two large signals, however, show a much higher relaxation rate, but again these two relaxation rates are within error of each other. dependence measurements were done, as the signals were much too weak.  N o temperature From eq. I V . l , the  relaxation expected from the R L M F of " G e would be 1.3 times that seen for " S i .  In fact the  observed relaxations for the two triples at 6 K are rather a bit less than that seen in a - quartz at the same temperature.  Given that the two triples A and B correspond to two distinct sites, what can be  said about the other two frequencies C and D? They could also be members of a triple whose third member, at a frequency of 21.77 or 42.19 M H z , is too weak to have been observed. Alternatively, they could be from two different sites, both of which have cylindrical symmetry and so show only one frequency; or if the two sites were more anisotropic, then the other four frequencies are not observed due either to low intensity or too high a frequency.  Another  92 hypothesis is that the two frequencies are for the same cylindrically symmetric site and represent two of the three possible frequencies from a site where the isotropic part is only about 1% of the vacuum value - in analogy to anomalous muonium in Si, Ge, and diamond [Patterson 84], but with an even greater reduction in frequency.  In zero-field this last case would be indistinguishable  from the first It is possible that any of these three or four sites may be associated with an impurity, particularly in view of the very weak signals present  This is a well-known phenomenon in EPR,  for example the germanium- hydrogen- lithium site in quartz [Weil 71].  Such an impurity  associated site might show a much smaller hyperfine coupling than normal muonium, and if the impurity had a nuclear moment this might split the line, an unresolved splitting resulting in a broad line. No determination of the site orientation in G e 0 single crystal.  2  was possible, as the sample was not a  Rather it was polycrystalline with a small number of crystals, so it was also not in  the powder limit - in which case all members of a triple would necessarily show the same intensity. Experiments with a large single crystal should give much more information, but if it is not possible to grow such a crystal with sufficiently high purity, a coarse crystalline powder could at least test the hypothesis that triples are present  3. M A G N E S I U M O X I D E A single crystal of MgO was studied over a temperature range from 6 to 295 K. was observed at all temperatures.  A run at 295 K at 10 G is shown in Figure 30; P  M u  Muonium = 30(10)%,  P =40(2)%, hence P =30(10)%. D  L  The relaxation rate as a function of temperature, as shown in Figure 31, is not what one might have expected.  Any crystal has some impurities which act as depolarizing sites for  muonium, or indeed the muon, by being paramagnetic; by being chemically reactive; or by creating a locally anisotropic hyperfine environment  A host site would not be expected to have an  anisotropic hyperfine interaction, as the material has cubic symmetry.  It is true, however, that the  Figure 31 Relaxation of the muonium signal in MgO versus temperature.  95 local symmetry need not be as high as mat of the crystal, as is the case with a-quartz.  Any  nucleus with a magnetic moment, such as hydrogen, will also depolarize muonium by the superhyperfine interaction or dipole-dipole broadening.  The expected temperature behaviour of  such a system is for the low- tempera tare relaxation rate to be low, since the muonium diffuses only slowly and so is unlikely to find a depolarizing site during its lifetime; at higher temperatures the relaxation rate increases as the muonium diffuses faster and so is more likely to find these sites; and then at still higher temperatures the relaxation rate declines again as the muonium is moving past these sites so quickly, or is in them for such a short time, as to avoid being depolarized. In MgO, neither of the two elements has a large percentage of nuclear moments ( Mg 10.1%, 25  17  0  0.037%), and the crystal structure is highly symmetric, so that at low-temperatures, where the muonium is static, relaxation should be slow.  The simple scaling from eq. I V . l predicts a  relaxation rate only 1.5 times higher than for a - quartz, about 0.3 MS" . 1  No attempt was made to  account for the difference in structure, and hence distances from the muonium atom to the nuclei will not be properly accounted for. concentrations.  Also, relative abundances were used rather than  The relaxation seen is about 20 times faster than this prediction, a surprisingly  large disagreement  As the temperature increases, the relaxation rate first falls, then rises to a  value comparable to the low-temperature value, at about 200 K, then falls again.  That the  relaxation rate rises, then falls with increasing temperature is consistent with the expected behaviour due to impurities; what is puzzling is the high relaxation rate at the lowest temperature. This may indicate that even at 6 K muonium is diffusing rapidly and finding a trapping site different from that found at the higher temperature.  A n exploration of the region below 6 K  might be helpful. The observed hyperfine is considerably reduced from the vacuum value, measured in the same way as for fused G e 0 , being only 3.4(3)GHz, in agreement with a more recent measurement 2  from SIN, 3849 M H z [Kiefl], a 14% reduction from the vacuum value.  This means that the  electron wave-function has been expanded by the lattice, the opposite of what is observed in Si0 , 2  and to a far greater extent  Such an expansion might suggest that the muonium is in a site  96 surrounded by M g " ions, just as the site in a-quartz is surrounded by O " ions, the former being much more ionic than the latter due to some covalency in the S i - O bond, thus accounting for the difference in magnitude of the effect.  Since MgO has the same structure as KC1, it is interesting to  note that there is also a reduction in the muonium hyperfine constant in KC1 of about 5% [Kiefl 84].  If the site is the same for both, this is consistent since the Mg** would attract much more  strongly than the K*. However, no temperature variation of the muonium amplitude in MgO was observed in this study, such as has been seen in KC1 [Kiefl 85].  If this is due to a thermally  activated process in KC1, as suggested in the paper by Kiefl et al. [85], then perhaps a wider range of temperatures would show the effect The model is then that muonium hops rapidly between equivalent sites even at low temperatures.  It would not be difficult to study this material at lower temperatures where this  model predicts a decrease in the relaxation rate.  4. S T R O N T I U M T I T A N A T E  SrTiOj was chosen as it met the criterion for having a low concentration of nuclear moments ("Sr 7.0%, " T i 7.3%,  49  T i 5.5%), and it also shows ferroelectric behaviour [Sakuda 71].  While not a true ferroelectric, nor an antiferroelectric, still its dielectric constant e is a strong function of temperature.  It was thought that there might be some connection between muonium  formation and the dielectric constant, so this seemed an ideal test case, where e could be varied at will, merely by changing the temperature.  A single crystal was used, and the levels of impurities  are given in Table III. Strontium titanate exhibits tantalizing behaviour for the would-be observer of muonium. Generally, a muon signal amplitude is, of course, anti-correlated with a muonium signal amplitude.  Naturally, then, the 100% muon signal amplitude observed at room temperature  precludes a search for muonium at that temperature.  However, the muon signal was measured as  a function of temperature from room temperature down to 4 K. The observed amplitude changed dramatically below 60 K , as shown in Figure 32 with a minimum at about 32 K .  The muon  1  1  1  1  1  0  CD  —  CO  CD  —  CD  •  Li  cd •  CD  o  —  •  u 0)  6  •  Eo CD CM  o  • ' —  CD  1  co • CD  1  1  co  •  CD  &  •  1 CO  LO  CD  ^ A m p l i t u d e o f the muon signal in S i T i O versus temperature, normalized to the maximum s  observed value which is consistent with a 100% muon signal.  98 relaxation rate is shown in Figure 33, exhibiting a monotonic increase with decreasing temperature with a possible anomaly at 35 K, not at the temperature corresponding to the minimum amplitude. SrTi0 has been studied extensively with other techniques for two reasons: its phase 3  transition at about 105 K [Willemsen 76, Kinase 82], and its ferroelectric behaviour [Sakuda 71]. The muon does not show any remarkable behaviour in the vicinity of the phase transition, although a more accurate measurement of the relaxation might be of interest.  At the lower temperatures  where the ferroelectric properties of SrTiOj become noticeable, the muon fraction declines quite suddenly.  Two sets of data were taken under slightly different conditions and they agree in the  details of the remarkable structure of the muon amplitude curve. It is interesting to note that the minimum of the curve occurs around 32 K, which has been reported as the ferroelectric transition temperature [Cowley 64].  At any rate, the dielectric constant is a strong function of temperature,  increasing monotonically with a decrease in temperature throughout the region where the muon amplitude changes. As noted before, in a non-magnetic material the only known mechanism for a decrease in the initial muon amplitude is the formation of muonium or muonium-like species.  If the  muonium is not observed directly, then it has been rapidly depolarized. If muonium is formed, then reacts chemically to form a diamagnetic compound, this can be identified, since the observed muon fraction will depend both on the rate of reaction and the applied field, increasing with the former and decreasing with the latter [Brewer 75].  No such field dependence was observed. A  search was made for muonium in SrTi0 at several temperatures, particularly at the minimum in 3  the muon amplitude at 32K, but none was observed. If the loss of muon amplitude is evidence for muonium formation, why then is it a function of temperature?  The temperature dependence of the dielectric constant is the most obvious  source, but just how could it bring this about?  One might naively expect that since the attraction  between the muon and electron goes as 1/e (where e is the dielectric constant) muonium would be less likely to form where a high dielectric constant was present; that is, at low temperatures. However, it must be remembered that this c is the M k e, which does not necessarily apply on a  99  CD  CD 00  CD  QJ 3  (0 CD "3"  OJ  a 6 OJ  CNJ  CD  CD  CD  LO  CD  CD  CO  CD  CNJ  Figure 33 Relaxation of the muon signal in SrTiOj versus temperature.  —  100 microscopic scale.  Moreover, the same reasoning would apply to whatever the muon was seizing  the electron from. The ferroelectric transition is associated with a relative shift between the Sr  ++  and the TiCy~ ionic sublattices, so any muonium formation could be associated with this lattice change.  Perhaps "holes" are opened up which the muonium atom could occupy. Since the strontium titanate lattice is much more ionic, the electric field gradients should be  far larger than in quartz and germanium dioxide, so any muonium formed might be expected to have a hyperfine constant different from the vacuum value, and highly anisotropic as well.  If this  were so, then the conventional transverse field search for muonium used frequently in this thesis might well fail, if for no other reason than that the frequency might be split and shifted far from the expected 1.39 MHz/Gauss.  If the Forward-Backward spectrum taken at 32 K in a field of 7G  shown in Figure 34 is examined closely, a "dip" can be seen at early times (<300 nanoseconds), with an amplitude about half the missing amplitude. This would be consistent with the missing muon polarization being due to muonium formation followed by rapid relaxation of the muonium signal, with two exceptions: its presumed frequency is only half that expected in this field; and it is not seen in the Left-Right spectrum, as expected in a transverse field.  The dip is seen in other  fields as well, almost unchanged; for example at 0 G (Figure 34) and 25 G (Figure 35).)  It is not  seen in runs of other materials taken before and after these runs, hence it is unlikely to be due to a defect in the data collection system.  It is hard to say what this dip could be exactly, but almost  certainly it is a remnant of the disappearing muonium.  It has been suggested [Brewer 85] that an  extremely anisotropic hyperfine interaction might be responsible. It is worth noting that the recombination time of spurs formed in the slowing down of the muon can be expected to increase with the dielectric constant of the material [Mogensen 74].  If  the spur species are thermalized, then lower temperatures would also increase the recombination time. If there were sufficient time to permit an actual buildup of the radiolytic species, this might explain muonium formation and relaxation in SrTi0 in much the same way as the spur model does 3  in water, [Percival 84].  The result of a given experiment would then depend on variables like the  incident muon rate and the size of the sample.  It would be of considerable interest to study other  101  Fi&ure 34 Early time behaviour of the asymmetry in S r T i O , at about 32 K in zero applied magnetic field (top) and 7 Gauss. Note that the anomaly does not change its apparent frequency.  102  0  .1  .2  .3  4  T I M E  .5  .6  .7  .8  .9  1  (microseconds)  Early time behaviour of the asymmetry in S r T i O , at about 32 K in 25 Gauss applied magneuc field. The Backward-Forward spectrum (top) still shows the anomaly at the same frequency as jn 7 P T Ofieldand seven G . but the L e f t - R i g h t does not show it at all. N o r did it at 7 Gauss.  103 ferroelectrics to see if the muon shows similar behaviour to that discussed above.  For example, it  might be possible to test the hypothesized buildup of radiolytic species by varying the incident muon rate or simultaneously bombarding the sample with other particles, such as electrons.  5. RIJTILE In addition to its small concentration of nuclear moments, rutile (Ti0 ) was also studied 2  because hydrogen is a well-known impurity in it [Chen 84]. impurities as given in Table III, was obtained.  Accordingly, a single crystal, with  Rutile, like strontium titanate, showed interesting  behaviour in that the amplitude of the observed muon signal decreased with decreasing temperature as shown in (Figure 36.)  However, the effect is much less dramatic, and the small  dependence was not apparent when the sample was run, so a search for muonium was made only at the highest temperature, which now seems the least likely place.  N o muonium was observed.  That the dependence of amplitude on temperature is slight can be seen by looking at Figure 37, which presents asymmetry plots at 295 K and 20 K , in a field of 160 G . The data may indicate the formation of muonium which rapidly reacts to form a diamagnetic compound, like the O H " observed by Bates et al. [79], for example; or perhaps there is a time delay before the M u atom ionizes to u\ in the same way that Schumacher [80] suggests H atoms in T i 0 diffuse before 2  ionizing.  If the reaction rate is comparable to the triplet precession frequency and if the rate is  temperature dependent, the resultant muon amplitude will vary with temperature as well [Brewer 75].  Another signature of this effect is that the amplitude will decrease as the field increases  [Percival 78].  However, only one run was done at a different field, the room temperature search  for muonium at a lower field. Curiously, rather than showing an increased muon amplitude it shows a decreased amplitude. In principle this could be due to the 160 G field being sufficiently high that the muons are steered away from the target as mentioned in Section II.5.  The residual  energy of the muons after the thin counter is about 3 MeV which corresponds to a momentum of 25 MeV/c.  In the applied field of 160 Gauss, such muons would have a radius of curvature of 500  cm, and so over the 5 cm distance from the thin counter to the target, the muons would "drift" less  104  CD CD CO  CD LO CM  CD  ^  CD X CM  w 0J  CD LO  CD CD  s  CD LO  LO CM CM  CD  CM CM  LO  CD  LO  CM  CM  CM  CD  CD CD  LO CO  CM  ^ T t a f m u o n signal amplitude in T i O , plotted as a function of temperature. Note the expanded scale.  The asymmetry seen in T i O , at 295 K (top), and 20 K_ The solid plotted line is a fit to the amplitudes as plotted in Figure 3C  106 than 1 mm, much less than the size of the crystal, which was a 1.5 cm cube.  Thus, this effect  seems unlikely to be important The missing fraction of about 30% is presumably clue to depolarized muonium, although the source of the depolarization is not at all obvious. The only species with significant abundances of nuclear moments are " T i at 7.3% abundance, and T i , with an abundance of 5.5%, which, 49  although somewhat higher, are not so much different from quartz as to make the signal utterly vanish.  6. B E R Y L A N D ZIRCON As part of the survey of materials, B e A l S i 0 , or beryl, and ZrSi0 or zircon were also 3  2  6  18  4  studied briefly, beryl as a single crystal and zircon as an assemblage of single crystals.  Both were  of geological origin. Since both the zircon and the beryl are silicates, it was expected that they would behave somewhat like quartz in their ability to form muonium and perhaps in providing an anisotropic site for i t  The EPR of atomic hydrogen has been observed in a geological sample of  beryl [Andersson 74] at room temperature where the observed hyperfine interaction has cylindrical symmetry; the hyperfine constant is reduced from the vacuum value by 1%, and the anisotropy is about 0.06%. This is similar to the behaviour of muonium in room-temperature a - quartz. The beryl shows some sign of muonium, but the time available to study it was too brief to test more carefully, so no P j ^ is quoted. u  The zircon shows a large missing fraction, but no  muonium, and so the system may be worth further study with a larger, purer sample.  Further  pursuit of these two systems with uSR is desirable, if large samples of higher purity become available.  7. D I A M O N D  Diamond is a semiconductor like its structural analogues silicon and germanium, although less commonly used for that purpose.  Since muonium has been observed in silicon and  germanium [Patterson 84], it was hoped that a study of diamond would prove fruitful.  As  107  mentioned earlier, two different species of muonium are observed in silicon and germanium: an isotropic species with a hyperfine frequency reduced by about 50% from the vacuum value, and a species, called anomalous muonium, exhibiting cylindrical symmetry, with a hyperfine frequency greatly reduced from the vacuum value.  Anomalous muonium is more difficult to observe in  germanium than in silicon due to the G e nuclear moments [Estle 84]. 73  Diamond was very briefly  studied by Swanson [58], who saw a reduced initial muon amplitude. The first sample of diamond studied in this thesis was a large (about 1.5x1.5x3 cm) very impure intergrowth of crystals of the class known by the inelegant name of bort  It showed  absolutely no /iSR signal at all, the only non-magnetic material in which that had been observed. A second, smaller (1cm ) sample of higher purity was also examined but no muonium was seen, 3  and the muon signal was consistent with a muon background in the sample holder.  Some time  later, a small (7x7x0.5 mm), clean, Type IA diamond was obtained from De Beers Diamond Research Laboratory.  In the meantime, however, Holzschuh et al. [82] at SIN had studied a  somewhat larger diamond (9.2 g versus 80 mg), as well as diamond powder, and saw muonium at low temperatures, their signal vanishing well below room temperature. of interest to study the new IA diamond at low temperatures.  Naturally, therefore, it was  Again serendipity took a hand: the  cryostat was not yet cold when the time came to do the experiment and so the first run was done with the diamond at room temperature.  I then made the first observation of muonium in a room  temperature semiconductor (Figure 38). A temperature study down to 6K was also performed and the relaxation rate is plotted in Figure 39, along with those from SIN, which are remarkably different  Both the diamonds are of  the same type, with nitrogen as the principal impurity. Their explanation for their results was that muonium was converting to anomalous muonium, this rate increasing with temperature [Patterson 84b]. Patterson et al. saw anomalous muonium up to 1080 K [84b].  No anomalous muonium  was observed in the present sample, but in part due to the size and shape of the sample, it was put in in an arbitrary orientation for which there may be as many as eight lines sharing the  108  0.15  T i m e in /is ( 8  ns/bin)  re Asymmetry in diamond. Note that the room temperature (top) muonium signal relaxes  38  ,.  more slowly than the 6 K signal.  \  109  LO CO CD CD  CO CD  CM  W  CD  OJ  S  3  (0  CD  JH  LO  OJ  CD CD  0) E-  LO  00  CNJ  CD  (_srT)  y  "•TSata of the muonium signal in diamond versus temperature. The diamond shaped points •re the TRIUMF data, the squares are SIN data, a  110 polarization; so it. is not surprising that anomalous muonium was not observed, and this should not be taken as an observation of absence.  The isotropic hyperfine interval was measured to be  3693± 83 M H z at room temperature, which is consistent with the SIN result of 3711± 21 M H z [Holzschuh 82] measured at 6 K, so it would seem that both experiments are looking at the same type of muonium.  It is not at all clear however, how the SIN model, which was also thought to be  applicable to silicon and germanium, will mesh with the present thesis data.  Both diamonds were  of the same type (IA), but these classifications are rather broad and so they could have rather different levels of impurities.  It is usual that the sample with the highest relaxation rate has the  most impurities; however, since the relaxation rate in the SIN sample is lower at low temperatures and higher at high temperatures than the rate in the T R I U M F sample, as in Figure 39, it is rather difficult to assign this difference to impurities.  Looking at the T R I U M F data alone, what can one  speculate might be the source of the relaxation? moments seem unlikely to be a source.  There is only 1.1%  n  C present, so nuclear  Looking at Figure 39, a plot of the relaxation rate as a  function of temperature, one is most reminded of Figure 25, a plot of the relaxation in fused quartz, over the same temperature range, which was attributed to random hyperfine anisotropies.  If the  muonium in diamond were not strictly isotropic, but instead showed some anisotropy, perhaps with symmetry along the four <111>  axes, as is true for anomalous muonium [Patterson 84], then the  small applied field used would result in a large number of peaks for static muonium, perhaps eight, centered at the normal muonium frequency. relaxation.  A n unresolved splitting would give an apparent fast  If the muonium were diffusing it might average these multiple sites, resulting in  motional narrowing. This scenario does seem rather far-fetched, but a zero-field study at low temperature would quickly confirm or deny this model. The most obvious failure of the model is that, in contrast to the SIN model, it "explains" the T R I U M F results but not the SIN results. The EPR spectrum of hydrogen has been observed in diamond [Karaseva 83], where the species observed is isotropic and has a hyperfine constant only 0.4% less than the vacuum value. Neither of the two species of muonium correspond to this.  It has been suggested [Estle 85] that  the hydrogen atoms are not in the bulk of the diamond, but rather on the surface.  This would  Ill explain the hyperfine constant, but perhaps the spectrum should have resembled that seen for hydrogen on the surface of fused silica [Papp 75], showing a broadening due to anisotropics.  112 F. CHAPTER VI: S U M M A R Y A N D CONCLUSIONS The cylindrically symmetric Hamiltonian has now been observed for muonium not only for room-temperature a-quartz, but also for ice [Percival 82], and for anomalous muonium in the semiconductors diamond [Holzschuh 82], silicon [Patterson 78], germanium [Graf 79], gallium arsenide, and gallium phosphide [Kiefl 85].  In the semiconductors, the hyperfine parameter A is 0  sharply different from the vacuum value, and A A is comparable to A . The two insulators have 0  values of A comparable to the vacuum value, while A A is less than a 0.1% perturbation on A . 0  0  For both a-quartz and ice, the Hamiltonian z-axis has been shown to correspond to the principal symmetry axis of the crystal, although the signs of A A are opposite [Percival 82].  For the  semiconductors the axis of symmetry has been shown [Patterson 78] to be the four <111>  axes,  which complicates the spectrum even further. It is important to note that isotropic muonium is also observed in the semiconductors, with a hyperfine constant reduced from the vacuum value by about 50% [Patterson 78, Graf 79].  The  observation of two such widely different species at the same time has not, as yet, been satisfactorily explained.  The anomalous species is not observed for hydrogen, so it seems that it may be a  metastable state, decaying either to the normal atom or to a diamagnetic species. be due to impurities or imperfections in the lattice.  Its presence may  Sahoo et al. [85] now give fairly convincing  calculations to show that anomalous muonium is associated with a doubly positively charged vacancy. Their model works very well for silicon, but not as well for diamond. In this thesis work, a Hamiltonian new to jxSR was discovered, explained, and observed in two systems, a-quartz and hexagonal germanium dioxide; one with a completely anisotropic hyperfine interaction  This Hamiltonian has not been seen in any other /aSR study, and it is also  unusual in hydrogen atom EPR work [Weil 84].  This Hamiltonian allows, in principle, useful  geometric information about the site to be extracted from the spectra.  Since muonium atoms can  be observed in times shorter than those required for hydrogen atom studies, and since muonium may form in materials where hydrogen may not even be present muonium can be used to provide otherwise unobtainable information about the system and about what the hydrogen atom might do  113 in such a system, provided it really does act like an isotope of hydrogen. The study of muonium in low-temperature quartz herein provided the first evidence that this is indeed true in a solid. A summary of representative results obtained in mis thesis work is given in Table II.  It is  now clear that nearly every non-magnetic insulating or semiconducting material shows a reduced muon amplitude due to muonium formation. Chapter IV.  This agrees with the statement at the beginning of  That muonium is not directly observed with conventional techniques is easily  explainable for those materials containing nuclei with significant abundances of nuclear moments. That muonium should be unobservable in so many other materials is surprising.  It might be  suspected that impurities play a role, as they certainly do in the semiconductors silicon [Albert 84], germanium [Clawson 81] and gallium arsenide [Kiefl 85].  In a semiconductor, though, impurities  will, due to their effect on electronic properties, have an influence over a much longer range than they would in an insulator. Also, why would quartz, no purer than the other materials, show such a large, long-lived muonium signal?  Perhaps muonium in these others diffuses rapidly to  depolarizing impurities. Why then is a small, slowly relaxing, signal seen in germanium dioxide? Perhaps the observed muonium is trapped, not in some intrinsic site as it is in low-temperature a-quartz, but at some non-depolarizing impurity. Such a site may have been missed in other materials, such as zircon and rutile, as fewer events were taken. more vigorously due to the expected analogy with quartz.  Germanium dioxide was pursued  If this hypothetical diffusion were  strictly thermally activated, then it might be possible to observe muonium at lower temperatures. This might also be achieved by lowering the level of impurities. Both of these courses could be pursued, perhaps simultaneously. The nonobservance of muonium in these materials is reflected in the large missing fraction which is found even when some muonium is present quartz, where there is no missing fraction. fraction [Percival 85].  The only exception to this statement is  In ice, recent analysis shows that there is a missing  In this case, as for liquid water, it is suggested [Percival 85] that the missing  fraction is due to depolarization of muonium by spur electrons.  The present results for strontium  titanate might be construed as supporting the spur model, but such an interpretation is far from  114 unequivocal, and it must be stated that a general mechanism for muonium formation and the missing fraction in solids is not yet established.  It is likely that no single mechanism is sufficient  in contrast to the situation in gases [Arseneau 84].  The large diamagnetic signals seen in strontium  titanate, rutile and magnesium oxide may be related to the formation of MuO~ as suggested in the Fe 0 2  3  system by Boekema [84], or the H~ seen in MgO by Gonzalez et al. [81] in MgO, but this  provides no information on the missing fractions.  The very small muonium signals seen in fused  and hexagonal germanium dioxide compared to quartz is perhaps the biggest challenge to any theory attempting to explain the missing fraction. If it had only been one sample that showed this, it might be passed off as a poor sample; since two different samples, a fused sample and a crystalline sample, both have a large missing fraction, it must be explained in some more interesting way. Several suggestions have been made in the body of this thesis as to future directions for the study of muonium formation. Stoneham [83] has also made some suggestions as to how to pursue the problem. A few of these suggestions are: 1.  that the muon fraction may increase in alkaline earth fluorides close to their melting points due to anion Frenkel disorder;  2.  that  tne  - i Y [ / * T 3 t i ° should show interesting behaviour near the metal-insulator transition P  l  >  ra  in silicon - a study of the metal-nonmetal transition is currently planned for T R I U M F [Noakes 85]; 3.  and the muonium fraction should increase in materials in which the preferred u* site is blocked by hydrogen, as in the case of CaO with increasing amounts of water added to form Ca(OH) , and by sodium ions as in the case of (Na 0) (Si0 ) where he expects the muonium 2  2  2  x  y  fraction to increase with increasing amounts of Na. In the last case, since the muonium fraction is already so high, and since Na will depolarize muonium, as would hydrogen, it might not be feasible.  In the case of the CaO study, it is  encouraging to note that hydrogen atoms have been observed in this material [Henderson 72]. Stoneham does not address the missing fraction.  115 The latest weapon in the arsenal of muonium observation is the use of high transverse fields to decouple the superhyperfine interaction [Kiefl 84].  If the depolarizing impurities  hypothesized above depolarize via some magnetic coupling, then a sufficiently high transverse field should decouple the muonium and allow its direct observation.  If a chemical reaction is  responsible for the loss of signal, then the high fields would not help. High fields will also be of no service if the depolarizing mechanism is an anisotropic hyperfine interaction.  In fact, if  low-temperature a - quartz with an arbitrary orientation were placed in a high transverse field, six lines would be produced rather than the usual two [Kiefl 84] observed for other materials.  Given  the A N rule, this would mean nine times as many events would be required to achieve the same J  quality of signal.  The use of high transverse fields is by no means a panacea, but rather it is  complementary to work at low and zero fields. Even in the high fields used to date, 12 K G , there are still missing fractions, although it is more difficult for the fractions to be measured accurately. Muonium was found to be behaving as an isotope of hydrogen in a -quartz, giving further credence to the use of muonium as a proxy for hydrogen, particularly in systems or temperature regimes where hydrogen is difficult or impossible to study. The main difficulty in pursuing various of the experiments which have been suggested herein is that every proposed experiment is in competition for beam time not only with other M S R experiments, but with experiments in more esoteric branches of physics.  Hence the proposed  "muon magic beamline", as discussed in Appendix 1, which would allow up to five times as many experiments in a given amount of time, would permit the exploration of many new phenomena, including the questions raised in this work, than might otherwise be possible.  116 G . APPENDIX 1: A M U O N M A G I C B E A M L I N E FOR T R I U M F  1. INTRODUCTION A magic muon beamline is one which delivers a single muon on demand, whenever demanded, and delivers none at any other time.  This cannot, as yet, be done.  However, it is  possible to turn off a flow of muons either at a pre-determined time or on demand. When turned on, it is merely a matter of waiting until a muon comes along. The term "magic beamline" will be used to mean a beamline which can be switched on and off, with delays for both. It is envisioned as an "add-on" to a present beamline, probably either M15 or M20, associated with the D C separators presently there or possibly upstream.  Ideally, the muons rejected by the switch would  be sent to an alternate beamline.  2. STATE O F T H E A R T Conventional M S R allows the muons to enter the target at random intervals.  The rate  limiting factor is the necessity of distinguishing which positron came from which muon, meaning only one muon may be allowed in the target at one time, any "second muon" events must be discarded. The maximum muon rate is, therefore, N =  where T is the length of the data gate  used. The resulting maximum good muon rate N = N/e. That is, at the optimum incident rate, only 1/e (37%) of the muons are acceptable.  Therefore, N = g  If one had a magic beamline,  by turning it off every time a muon entered, and turning it back on at T, the maximum good muon rate would be N"  = 1/T.  This represents an improvement of a factor of 2e, so a magic beamline  could improve the real data rate by a factor of 5.4.  In a typical experiment, T=10 microseconds,  so the maximum rate would go from 18,000 to 100,000 per second.  3.  AT XERNATIVES At the present time, there are three other major ways of doing M S R with increased data  rates.  117  1.  Pulse M S R . Periodically, a pulse containing many muons enters the target and the positrons are detected. One of the problems is that the time of stopping cannot be well defined, another is that due to the large numbers of muons and other particles scattering through the positron counters there can be large distortions at the beginning of the spectra.  Probably the biggest single drawback  is that one needs a pulsed source of muons - one which has a period longer than the longest time range of interest  2.  Stroboscopic M S R . In this method, a transverse magnetic field is applied such that the muon precession frequency is an integer multiple of the frequency of the cyclotron. In principle, an infinite number of muons can be accepted.  One major disadvantage are that the magnetic field cannot be varied  smoothly, as most cyclotrons are unable to vary their frequency except by, for example, a 1:5 selector.  Another is the intrinsic linewidth of 1/7  for many experiments.  3.  , about 400 kHz, which is much too broad  Finally, M S R experiments in zero applied field are totally excluded . 41  Cern wire-chamber technique. Here, a wire chamber is used to monitor the decay positrons in order to determine their origin. In principle, it is also possible to monitor the incoming muon beam with chambers and determine where in the target they stop.  Then, since one can determine where the positrons  come from, it should be possible to logically "split" the target into many sub-targets only limited by the resolution of the chambers. This reduces the second muon problem noted above, as it less likely that two muons will enter the same small region of the target.  The limiting factors are the time required to  It might be possible to use a magic beamline to produce a time structure of whatever frequency is desired, thus removing some of the limitations mentioned. As the muons generated while the switch is off are thrown away, while with the regular stroboscopic technique all the muons generated are used, this might not be advantageous. 41  118 calculate the positions and the spatial resolution, which is limited by scattering.  Perhaps one  million muons per second could be tracked, but the spatial resolution might be only 5 millimetres, so a target 5 cm. on a side would be needed to get this high rate. of the interesting solid state targets are about 1 cm or less.  However, many  Again, though, at this rate, only  37% of the muons would be considered "good", so the good muon rate would be 370,000 per second.  Such a technique is not at all suitable for gas-phase work, as the stopping region is  quite ill-defined. Surface muons have only about 140 mg/cm range and are 100% polarized. For many 2  experiments their use is essential, as in thin solid-state or gaseous targets. means they are preferred for all M S R experiments.  Their high polarization  Using any of these 3 techniques with surface  muons may not be possible, as they scatter easily reducing the spatial resolution for (3), and their time structure is very broad, broadening the pulse in (1) and making (2) less sensitive.  4.  JUSTIFICATION  A magic beamline would work with surface muons for gas, liquid and solid targets in and out of cryostats just as conventional y S R now operates, but 5 times faster.  This makes it much  more flexible and much more powerful than any of the other methods for increasing the rate at which experiments can be done.  The availability of such a beam line would have greatly  facilitated the work done in this thesis.  Moreover, by increasing the real data rate by a factor of 5,  it will inevitably lead to the performance of experiments which are not presently feasible due to the high statistics required,  M SR at T R I U M F currently does about 2000 runs each year.  Some of the background present in M S R is due to muons which decay in the collimator or elsewhere.  With a magic beamline this source would be shut off and the remaining background  muons would then decay with the muon lifetime and not contribute to the background at longer times. The combination of higher data rates and reduced background would allow the time scale of M S R to be pushed further out, possibly rendering the measurement of 100 microsecond  119 relaxation times routine. If a beamline were split into two legs, and this switch were installed upstream of the split, it may be that rather than use a septum and have one magic leg and one ordinary leg, it might .be better to have two magic legs, where experiments take turns. For example, at this increase in data rate, experiments which presently take 1 hour for each run would take only 10 minutes. For some experiments, a change is made to the target for each run, for example, the temperature is changed, or the gas mixture might be slightly altered, meaning that at least 10 minutes downtime will occur. This leads quickly to the conclusion that for a magic beamline to be used to its full potential, it must be part of a two leg beamline.  Even if the switch were downstream of a split, during the  macroscopic downtime on the magic line, the other side would be an ordinary line. Of course, an even better scenario is one in which there are two legs which are magic simultaneously (a doubly magic beamline), and the rejected muons are just dumped. Naturally, the next step is to have two legs magic and one leg ordinary . . .  . If the T of interest is 10  microseconds, and the muon rate is 1,000,000 per second, it would be possible to service as many as 10 ports, all of them magic. A multiplexed muon beam!  5. DESIGN CONSIDERATIONS Clearly a magic beamline is desirable. How can one be made? muon beam through an estimated angle of 3 degrees.  It is necessary to bend the  This can be done by an effective electric  field of 400 kV (the potential in volts divided by the gap times the length) or magnetically with 5000 G - c m (the field times the effective length). Presently, at L A M P F , the TPC has an electrostatic separator for surface muons which can also be used as a switch. It does this by charging the plates of the separator to let the muons through then dumping the charge to steer the muons away. triggerable, switching in about 100 nanoseconds.  The on and off switches are  The problem in adapting it for a magic beamline  is that it can only be used with a frequency of 500 Hz. The application at T R I U M F calls for  120 100,000 Hz.  The limiting factor in the speed of switching and the frequency of switching is  dissipating the heat generated.  By accepting a slower switching time it would be possible to  increase the frequency, so, for example, if more of the charging devices were used, the same technology should work here.  This work, including the actual construction, should cost of the  order of $100,000. The present L A M P F device cost $20 - 30,000. Another method is to use a kicker magnet, a well-established technique used elsewhere to do fast switching on high energy proton beams.  The main difficulties are that although the  momentum of the surface muons is much lower, the required angle is much greater, and the beam envelope of a muon beam is much larger than for a proton beam. For example, the AGS fast kicker has a rise time of 150 nanoseconds, and generates 100 k G - c m , but the gap is only 0.5 inches vertically and 1.25 inches horizontally. A typical muon beam envelope will be at least 6 inches in diameter.  In addition, its repetition rate is only 45 Hz.  Fermilab estimates the lifetime of their 250 k G - c m kicker as greater than 10 pulses. 7  limit would be reached in 100 seconds of use on a magic beamline. hours.  This lower  10' pulses would take only 3  However, these fields are much larger than we would require, so perhaps both higher  repetition rates and longer lifetimes might be achievable. A report was commissioned by T R I U M F on the feasibility of using a 43.5 kHz kicker with a rise time of 20 nanoseconds as part of the injection to a kaon factory. The requirements for field and switching time were somewhat different from the application proposed here, but the repetition rate was nearly what is being envisioned. One of the results of the study was that it was predicted the tube would last about seven hours. At perhaps $40,000 per tube, this is not very attractive. A magic beamline with a switching time of 200 nanoseconds instead of the 20 nanoseconds for. the kaon factory might mean the lifetime would increase, but if only to 70 hours, it is still inadequate. A third method has been suggested by Syd Kreitzman of TRIUMF.  He suggests that if 2  sets of saddle coils producing a field normal to the beam were positioned successively along the beam line, each with a circularly polarized R F magnetic field, by suitable choices of their phases it should be possible to transmit the muons (phase=180°) or bend the muons away (phase=0°).  121 That is, a muon passing through the first set would experience a magnetic field bending it upwards, say, then the second set would either bend it downwards if the switch were on. or upwards if the switch were off.  Switching is then merely a matter of changing the relative phases of these 2 sets.  The switching time should then be very quick, as no energy changes are involved. A disadvantage of this is that it does not appear that any such device has as yet been tried anywhere. Saddle coils with R F fields are used commonly in N M R , but typically to produce fields over a much smaller volume. In such a large volume as, say 15 cm diameter by 20 cm long, a field of 125 G might be difficult to establish at the 1-50 MHz required. The actual frequency is not so important but the lower limit is determined by the switching time required, and the upper limit is set by the requirement that the field not change too much while the muon passes through. Surface muons only travel 8 cm per nanosecond.  For example, a  coil 20 cm long, with a frequency of 50 M H z would mean the field would change direction by about 46 degrees, reducing the effective field by 15%. The capital cost is estimated to be less than $150,000. An advantage of the first two methods is that they bend the beam in only one direction, so it may be possible to take the rejected muons and send them out an alternate final leg on the same beamline.  This would probably require a septum magnet to "amplify" the bend.  This is not  possible with the third method. The triggerability and rise/fall times of the switch are important considerations.  If it is  not triggerable, that is, if the on and off intervals must be preset and only variable macroscopically, this would decrease the efficiency somewhat perhaps a factor of two less than the triggered version. The rise/fall times must be less than the mean time between muons which should be at most two microseconds (M15 should have a total flux of 1 million muons/second), and preferably they should be kept to of the order of 200 nanoseconds or less.  If the mean time between muons is  increased, then the rise/fall times could be allowed to be longer, but all this is at the cost of a loss in efficiency.  One of the limits as to how fast the muons will stop or start after the switch is  flipped off or on is the time-of-flight of the muons. As noted above, their velocity is only 8 cm  122  per nanosecond, so a switch 8 metres upstream of the experimental target would mean the experiment would not notice the change until 100 nanoseconds later. time required to tell the switch to flip, and the switching time itself.  This is in addition to the This gives some motivation  for keeping the switch close to the target  Development of the magic beamline requires work in several areas.  All three of the ideas  need to be explored, as the D C - t y p e switches are presently used at frequencies  100-1000  times  lower than required and the R F - t y p e is untested and may not be feasible, although a bench test of the basic idea should not be too difficult Work should be done on the optics of these systems as applied to an existing beamline, either M 1 5 or M 2 0 .  M 1 5 may be easier due to the better optics design.  Some work should be done to see the effects of the switch on the M S R spectra, in particular the problem of an incomplete switch, i.e. if the muons are still transmitted at a level of, for example, 1 0 % , what are the consequences to the maximum good muon rate and the background?  A magic muon beamline would increase the real data rate of a typical conventional M S R experiment by a factor of 5 at a cost far less than that of an ordinary beamline.  123 H.  APPENDIX 2: SOLVING THE ANISOTROPIC HYPERFINE  HAMILTONIAN The Hamiltonian i s given by: H  T  "VM *" 1  =  where g  +  9  M  ?  ?  e B '  +  (A2.1)  and g are the muon and e l e c t r o n q - f a c t o r s e '  M  respectively, and  M are the muon and Bohr ( e l e c t r o n ) magnetons, b  I and S are the muon and e l e c t r o n spin and £ i s the h y p e r f i n e It can be s p l i t H  = H  T  + H  A  operators  tensor.  i n t o two p i e c e s , (A2.2)  B  where #  = I'A«S and  A  "B  =  " /i M ' 9  M  I  I  +  9  (A2.3) M  §  e B '  ?  (A2.4)  To solve i t , f i r s t choose as a b a s i s s e t , where e, u i s the order  of l a b e l l i n g ,  |X> = j7^( | aa>  +  | 00>)  |Y> =  \aa> - |/5/3>)  |Z> = ^(  |a/3> t | j3a>)  |0> =  |-a/3> - |j8a>)  (A2.5)  Note: 21  |X> =  |z>  21  y  |Y> = |0>  21 |Z> =  |X>  |0>  X  =  |Y>  y  1  y  1  1|0>  21 y  i | 0>  21  1  |X> = |Y> z  -= i | Z>  21  |Y> = |X>. z  I|Z> -=  21  X  -=  1  1|Y>  21  X  21  1|x>  21  X  ==  -i|Y>  2  I  z  -i|X>  2  I  z  |Z> = - | 0 > |0>  =  -|z>  1 24  2 S | |X>  2 S l |x> -= - i | 0 >  2 S  |X> =  |Y>  |Y> = - I 0 >  2 S | |Y> -= i | Z>  2 S | |Y> =  |X>  |Z> =  2S  |Z> -= - i | Y >  2S  |Z> =  |0>  |0> =- i |X>  2 S  |0>  |Z>  x  2S 2S  2 S  x  -  |Z>  y  X x  |0>  2  y  1  |X>  = -|Y>  2S  y  y  zl  z z  =  (A2.6) I«B  = I  B x x  + I B y y  + I B z z  I-B|X>  = i ( B |Z> + i B |0> 2 x y  I-B|Y>  = i(B  |0>  I•BIZ>  = i(B  |X> -  iB  I-B|0>  = i(B  |Y> -  iB  1  1  S-B  1  1  2. x '  = S B x x  + S B y y  S-B|X>  = i(B  S-B|Y>  = i ( - B |0> 2 x  S-B|Z>  = -j(B |X>  S-B|0>  = 4"(-B |Y> 2 x  1  1  1. In  + iB  the axis  y  iB  -  1  x  |Z> + B  |X>)  1  |Y> -  B  |X> -  B  |Z>)  + B  lY>)  |0>  + iB  1  SOLVE / / ,  1  z  1  |0>) (A2.7)  + S B z z  |Z> -  x  + B |Y>) z  y  IZ> + B  1  iB |Y> y  +  iB  y  +  THE H Y P E R F I N E  system  i n which  |X>)  1  B |0>) z  |X> + B 1  z  z  (A2.8)  |Z>) 1  PART A is diagonal  //»=A I S + A I S + A I S A xx x x yy y y zz z z Now t h e m a t r i x e l e m e n t s  can e a s i l y  be w o r k e d o u t :  1 25  X  V  >  (  " 4 ^ =  472  A  >  x x l ^  ( A  A  -  A  xx  i  ^  >  z z  )  y y  " y y  +  A  +  A  | Z  (  a  a  >  +  A  x x l  Z  +  l ^ >  t  l  t  a  >  a  a  >  -  A  y y l  a  a  >  +  A  z z l ^  >  )  )  = -J-(A - A + A ) |X> 4 xx yy zz 1  (A2.10)  = E |X> x  Y  V  >  =  4^  [ ( A  472  =  ( A  4-(-A  =  A  4  +  xx- yyH^> A A  xx  A  " y y " + A  xx  yy  z z  + A  |aa>-(A  Z 2  )  (  l ^ >  x x  -A  y y  )|aa>-A  ) |y> zz 1  (A2.11)  y  Z  >  " 475  ( A  xxl^  +  =  472  4  ( A  ( A  A  +  a >  +  A  A  A  A  A  A  a  >  y y l  y y ~ zz  y y "  p  y y l  +  x x l ° ^ +  xx  xx  )  (  A  " a  ^  a  >  zz' ^ A  "  z z l ^  a  >  )  +  l**>  )  zz |Z> (A2.12)  = E |Z> z  "  A  l  0  >  "  4 7 2  (  A  A  =  x x i ^ x x l  472 " xx (  4 " (  A  A  xx -  a  +  A  /  A  a  >  3  >  +  A  yy^  t t  >  "  A  " y y l ^ > A  y y ~ z z  y y "  A  )  (  +  l ^ >  Z I *> A  A  Z  A  z z ^ -  a  >  " )  \fia>)  zz |0> )  = E |0>  (A2.13)  Q  which  |^>]  - |aa>)  = E |Y>  ^AI  z z  c a n be w r i t t e n  i nmatrix  notation a s :  126  A  xx  -A +A yy zz -A  0  4  xx  +A  yy  0  0  0  0  +A  zz A  xx  +A  0  C  0  0  yy  -A  0  zz -A  xx  -A  yy  -A  zz (A2.14)  A lucky  choice  of  Check — the  trace  Check -  A  let jh  E  2.  E  3.  E  X  ,  Q  -  E  = A  -  E  "  E  -  E  "  E  "  E  0  energies  three  as  for  the  then  given  (  I  A  (  A  xx  . (-2A xx + 2Az z )  2- " xx  0  |(2A + 2A ) 4 yy zz  2-  0  •i(2A + 2A ) 4 xx yy  I  E  6.  E  2.  S O L V E T H E ZEEMAN  z  The m a g n e t i c  field  (  (  A  ( A  A  yy  xx  +  +  4  5.  z  get  i " yy  XX  E  y  are  yy  o  ->  .  levels  •  isotropic  case.  by:  yy  XX  z  0  - j A , just  4.  y  = A  ^^  3  one a t  y  X X  0.  = A  The t r a n s i t i o n 1.  =  xx  1  at  basis.  +  A  z > 2  A  zz>  A  z z  +  A  zz>  +  A  yy>  )  (A2.15)  PART part  is  given  by: (A2.4)  This  gives  127  (7  2  2 z ^B  2  e  -7  z 'e  M  )  u  (w 7 - 7  is  )  x e  V  -4B (7 +7 ) 2 y 'e 'n  (7-7  )  -2 i B x (7' e+7  )  M (A2.16)  where 7 = g M„ and 7 'e e B 'M  = q M  " M M  y  SOLVE THE TOTAL H  3.  n  S u b s t i t u t i n g eqs.  A2.14 and A2.16 i n t o A2.2 and  m u l t i p l y i n g both s i d e s by 4 f o r convenience, 4tf = T  A  xx  2B 2B  -A  yy  +A  z (7 e- 7 (7  zz  M  -7  -A  ) )  (7  +7  y 'e  M  xx  2iB  x 'e ' u  -2iB  2B (7 - 7 ) z 'e 'u  )  +A w  yy  (7  y e  +A  2B (7 - 7 ) x 'e 'u  zz  -7  'n  -2B ( 7 + 7 ) x 'e 'ju  )  2iB (7 +7 ) y 'e 'n  - 2 i B (7 - 7 ) -2B (7 +7 ) y 'e 'n x 'e 'u A +A -A xx yy zz 2B (7 +7 ) z 'e ' M  2B (7 +7 ) z 'e ' M -A  xx  -A  -A yy zz (A2.17)  which i s r a t h e r unpleasant, however a n a l y t i c a l expressions f o r the energy can be found e a s i l y f o r c e r t a i n s p e c i a l cases.  1 28  a.  Special ~ i.  A  xx  2B  z  -A  case:  Finding  yy  +A  the  2B  zz  (7 - 7 ) ' e 'u  -A  0  z  eigenvalues  (7 - 7 'e +A  xx  B = B = X ^  0,  4*  )  LL  yy  +A  zz  0  0  0  0  of  H„  0  0  0  0  A +A -A xx yy zz  2B  2B  -A  (7 z  )  +7  'e  z  \u  (7 +7 ) ' e 'u -A  xx  -A yy  zz  (A2.18) To  solve,  set  |4#  -  T  X^| = 0  (A2.19)  One s e t o f s o l u t i o n s i s : (A - A + A - X)(-A + A + A xx yy zz xx yy zz = X  +  2  (A =  X  2  X = ;{2A =  A  xx  -  zz  *  The  other  (A  + A yy  =  X  2  xx  +  zz  ±  A yy  A  zz  xx"  2  ~4A  z  A  yy  the  A + xx  zz zz  of  (A  2  ) 2  z  xx  A yy  4 B 2 ( 7 -7 )2 z ' e 'u  A ) zz  2  zz  -  x x  -  A  -  2  « ^ ^  +  -  4 B 2 ( 7 -7 )2 z ' e 'u  4B ( 2  T  -7  z 'e 'u  2  »]  xx  )  )2  1 6B| ( 7 ^ )  +  2  A  zz  -  + A ) yy  2  ]1/  2  }  (A2.20)  1 / 2  J '  2  0  +  /  2  i  is:  A yy (A  )  y y  =  levels:  solutions  X + A  A ) yy  " £ ( 7 ^ )  +  * „ > '  -XM-A xx  -  4(A  +  energy  A  set  -  2  [< xx-  +  A  X + A zz  [4A  [ ( A  2A  -  -  A + A )(-A + A + A ) yy zz xx yy zz  gives  ' " T«*„  xx  -  2A  z z  which E  X(A  X)  X)  - 4 B 2 ( 7 -7 ) z ' e 'u  -  4B ( +7 )2 z ' e 'u 2  7  2  =  0  129  X = ^(-2A ±[4A z z  =  A  " zz  1  [ ( A  - 4A  z z  V  +  xx  2  z z  +  +  4  4(A - A x x  i(-A  + [(A  x x +  A  y y  )  2  zz  E, = { { - A  - [(A  x x +  A  y y  )  2  2Z  E  3  =  A2.21 E  12  =  1 3  -  E! [ ( A  E  A  E  [ ( A  ~  - E I  ) 2  +  4  - E  - A ) xx yy  [ ( A  xx  +  A  yy  ) 2  +  4  B  2  [ ( A  ~  E  xx  + A  yy  ) 2 + 4 B  e  z  /  2  ~  E  2  ( A 2  2  %  ) ]  l  2  e +  7 ) l  /  2  2 2 )  }  l / 2  M  '  >  (A2.23)  }  energies:  )  2  ]  1  /  2  (  A  2  '  2  4  )  2  + 4B ( -7 ) ] z 'e ' u 2  2  1  /  2  7  z  (  7  e  +  M  7  )  2  ]  1  /  2  (A2.25)  }  [ ( A  2  + 4B ( -7 ) ] z 'e 'u 2  2  1  /  2  7  z  { 7  e  + 7  M  (A2.26)  ) 2 ] l / 2 }  3  = i{2A - [ ( A -A ) 4 zz xx yy  -  B|(7  1  2  E (,  +  E2 4  ]  4B (7  +  ^ ^ M  B  = -J-{2A + [ ( A -A ) 4 zz xx yy  '23  4  2  /  3  = |{2A + 4 zz  E 14  +  )  1  »]  2  xx" yy  E1 -  =  (  1 6B| <  +  and A2.23 l e a d t o the t r a n s i t i o n  " I E  Z V^  B  2  )  y y  xx  + A  yy  ) 2  +  4 B  2  + 4B ( -7 ) ] z 'e 'u 2  2  1  /  2  7  z(7  + e  2  l / 2  7 ) ] M  (A2.27)  }  EI  = 1{2A - [ ( A -A ) 4 zz xx yy + [(A +A ) xx yy  2  2  + 4B ( -7 ) ] z e 'u 2  2  1  /  2  17  + 4 B ( +7 ) ] z 'e 'u 2  2  7  1 / 2  }  (A2.28)  1 30  1  34  E,  = E= i[(A  2  +A  xx  yy  )  2  + 4B2(7  -7  z 'e ' M  )  2  ]  1  /  (A2.29)  2  Check: f o r B =0, these agree with the z e r o - f i e l d z  expressions. Now, i f the a n i s o t r o p y i s s m a l l , A , A and A are ' ^ ' xx yy zz J  approximately  equal, and i f B ( 7 - 7 ) i s much smaller than 2  A , x x  which of these t r a n s i t i o n s  without,  0  fJL  i s "normally"  observable,  f o r the moment, c o n s i d e r i n g t r a n s i t i o n  intensities?  That i s , what i s the order of magnitude of E ^ j ? A - A or l i k e A + A ? xx yy xx yy E  1 2  E „ 3  Is i t l i k e  i s c l e a r l y s m a l l , and t h e r e f o r e should be observable. i s c l e a r l y l a r g e , and t h e r e f o r e l i k e l y Observable E  1 2  E  1 3  E  3  2  unobservable.  Unobservable E,4  E « 2  E„ 3  which i s q u i t e a n a t u r a l grouping, as i n z e r o - f i e l d the s t a t e 4 corresponds the  to the s i n g l e t and the s t a t e s 1,2,3 to  triplet. i i.  F i n d the T r a n s i t i o n  Intensities  To f i n d the i n t e n s i t y of the E  1 2  transition,  first  find  the a p p r o p r i a t e e i g e n v e c t o r s . To f i n d  s u b s t i t u t e E, i n t o the top l i n e of the s e c u l a r  equations corresponding to A2.18.  131  (A  xx  -A  yy  +A  -A - [ ( A -A ) zz zz xx yy  + 2B_( -7 )C 7 e  M  2  2  2  + 4 B ( 7 ~7 ) ] z 'e M  1 / 2  )C  = 0  2  1  (A2.30)  2B  C, =  c  A  xx  ^A  - [ ( A  yy  xx  ^A  yy  )  2  + 4B  2  1  ]  /  2  2  (A2.31) where B = B .  C  c,  = ±  A l s o , C?+C| = 1  4B {A  1-1/2  :  ^ A -[(A ^ A )2+4B ] xx yy xx yy 2  (A2.32)  + 1 1 / 2  }  2  2B  = +  A  xx  -A  - [ ( A -A ) + 4 B ]T 7 2 yy xx yy 2  4B  2  -1/2  2  (A2.33)  + 1  { A - A [(A-A)2+4B ] xx yy . xx yy 2  1  /  2  }  2  Look at l i m i t s : as  B approaches 0, C  2  approaches 0  C, approaches ±1 as  B approaches », C  2  approaches ±—^  d  approaches  To f i n d \[/ , s u b s t i t u t e E 2  2BC,+(-A 1  A r  C l  _  2  +^  from A2.21 i n t o A2.18.  +A +A -A + [ ( A - A ) + 4 B ] ^ C xx yy zz zz xx yy 2  1  2  2 2  2  2  -A - [ ( A -A ) +4B ] xx yy xx yy _ _ -  From the o v e r a l l  2  = 0  l / 2  normalization,  „ c  (A2.34) 2  1 32 {A C, = ±  xx  -A  2  yy  4B  A  2  - [ ( A -A ) +4B ] xx yy  - A  l / 2  }  2  -1/2 (A2.35)  + 1  :  - [(A  - A  )  + B ] 1/2,  2  2  C, = ± 2B A X  xx  - A  + 4B l]/ 2  2  2  - [ (A - A ) xx yy  yy  -1/2  n  + 1  (A2.36)  2B Now,  f o r \j/ as 2  B approaches 0, C  approaches ±1  2  C, approaches 0 B approaches °°, C  approaches  2  approaches  +^  So for B=0, a |1> = |2> o \2>  (A2.37)  = |1>  z  where a i s the P a u l i spin matrix, For B=°°, a_  1> =  7  l ( a  z  | X >  1( | Y >  -  -  72  a |Y>) z  |X>)  1> ff  zl  2>  =  7? zl U  = |2>  X  >  +  a  zl  y  >  )  (A2.38)  Since there i s no matrix element f o r a t o connect l e v e l s 1 z and 2 f o r B=°°, and s i n c e  1 33 a  and a  x  of E , B=0.  2  cannot connect  y  1 and 2 f o r any B , the i n t e n s i t y z  i s equal to zero.  However, t h e r e i s an i n t e n s i t y f o r B  To g e n e r a l i z e , set  4K  A  xx  - A  •=K. yy  "1/2  :  2K X> + 4K ]1 "7~ 1-[1+4K  + 1 2  {!-[1+4K ]  1 / 2  }  1  2 2  2  Y>  2  (A2.39)  -1/2 2  {1~[1+4K ]  1 / 2  }  r  2  2  1 -[1+4K ] 2K  W  7  1 / 2  |X> + |Y>  (A2.40)  > =  Check: does < ^ | ^ i  0 ? Yes, so a l r i g h t  2  so f a r ,  In g e n e r a l ,  a_|1> = N,  w  2 R  2  1-[1+4K ]  2  l  Y  " l  >  x  >  1 / 2  "1 -[ 1 +4K 1/2 ] Y> + 2K 2  a |2> = N z 1  2 2  1-[1+4K ] 1/22  2K  <2|a |l> = N T N , z  2  1-t1+4K ]  N  1N  T72 " —  2K  4K -1-1-4K +2[1+4K ] 1/22  =  (A2.41)  X>  2  2  2  2  T72  {!-[ 1+4K ] "}2K 1-[1+4K ] 1/2 2  =  -N,N  2 2  {1-[1+4K ]  1/2  }K  1 34 (A2.42)  K  Now, :  N,N  +  4K  2 +  =  2  2  1+4K ]  1 / 2  }  2  -1/2  N,N  2  {1-[1+4K ] 2  1-[1+4K ]  1 7 2  2  2  1/2  2  } 1  1  /  2  2  }  -1/2  1/2 x ' x+1  2  1/2  }  1+4K -[1+4K ] 172  ;  2  2  ., 1-y aside: = — — = y -y  1-y 1 —r-r* —r = — -y(1-y) y 4  a  N,N  1 / 2  2+8K -2[1+4K ]  K{!-[1+4K ]  N,N  2  r(x+i ) i x  2  2K  =  1 / 2  4K'  2 + x + x  aside: =  2  + (1-[1 4K ]  K  = -  2  -  2  [1+4R ]  (A2.43)  T72  S u b s t i t u t i n g A2.43 i n t o A2.42, we get |1> = - T-pr ' ' (1 4K )  (A2.44)  <2\a  z  2  1 / 2  +  To check, <1|o  I2>  =  NTN  2K  2  2  2  1-[1+4K ] 1 1  1 / 2  1-[1+4K ] 2K  1/2-  131  To [A  find  xx  +A  substitute  yy  -A  zz  -{-A  zz  E4  +[(A  into  3 of  2  +A  xx  line  2  1  2  2  From  B  = B  2  *  z  (7  e  +7  ju  1  1  +A  yy  +[(A  (A2.45)  2  7 1 / 2  xx  -1/2  ;  +  (A2.46)  2  }  4B2  * xx  -1/2  + 1  +A  2  yy  ) +4B|]  1 / 2  1  }  2B-  X A  xx  +A  yy  +[(A  xx  +A  2  ) +4Bf]  (A2.47)  1/2  z  yy  i l L {A  xx  +A  yy  +[(A  xx  A  xx  +A  yy  4K  +[(A  xx  + 2  +A  2B  X  0  2  i u  {A  =  2  2  =  21  C +C =1  U„„+A„„+[(A„„+A„„) +4Bl] xx yy xx yy  C,  C  1  2  ).  the normalization,  C,=±  }]C +2B C ' ' 1  z  • [A +A +[(A +A ) +4Bl] / C +2B xx yy xx yy z  Where  l / 2  ) +4B|]  yy  A2.18.  yy  ) +4B|] *  1 / 2  ;  2  +A  yy  ) +4Bi]  T72  }  ^  n-1/2  2  Z> +  |0>  z  -1/2  2  "-T7272 {1+[1+4K13 "}  +  1  , /  X  2K  2  T72 1+[1+4K1]  |z>  +  | 0 >  (A2.48)  1 36  where K-  A  B  2  + A xx yy  Similarly,  1^3 =  {1+[1+4KJ] 4K  1/2  }  2  1"  +  1 / 2  1+[i+4KJ]  7  2K  1 / 2  Z> - |0>  :  (A2.49)  Since i b , and \p are composed only of |X> and | Y>,  and s i n c e  2  the only a which connects  |X> and |Y> i s a , z  i t f o l l o w s that  for any magnetic f i e l d p a r a l l e l t o the z a x i s there never be an E  s i g n a l observed  1 2  will  i n the x or y d i r e c t i o n s .  Since \£ and \p are only made up of | Z> and |0>, and s i n c e u  3  the only a which connects frequency  E  3 U  them i s  , even the high  s i g n a l w i l l have no i n t e n s i t y  i n the x or y  directions.  b.  Solve So  f o r three  sites  f a r , only the simple case of a s i n g l e s i t e has been  considered.  Now c o n s i d e r what happens i f there are three  s i t e s r e l a t e d by a one hundred twenty degree r o t a t i o n around a symmetry a x i s ?  The eigenvalues  and t r a n s i t i o n  intensities  must be c a l c u l a t e d f o r each s i t e s e p a r a t e l y , then the s p e c t r a summed to give the e x p e r i m e n t a l l y observed To do t h i s ,  i t i s u s e f u l to take the magnetic f i e l d  spectrum. i n terms  of the e x t e r n a l c r y s t a l axes, x', y' and z', and transform it  t o the i n t e r n a l axes, x, y and z.  There i s a  transformation three  from x',  s i t e s , D(0),  y',  z' to x, y , z f o r each of  D(+120), and  w i l l be used to represent  D(~120).  a l l three  The  the  notation  D(n)  cases.  D(n) = cos ( <j>+n )cos\p -  sin(tf>+n)cos^'  cosflsin (<j>+n) sin*//  sint9sin<//"  + cos#cos (4>+n) sin\//  - s i n (4>+n) sinxp  -cos (</>+n) sin<// - cosflsin ( 0+n)cos\//  sinflcos^  + cost9cos(0+n)cos^  sint9sin (0+n)  cos0  -sintfcos (<p+n)  (A2.50) 6,<t>, and  i// are  d e f i n e d c o n v e n t i o n a l l y as i n C l a s s i c a l  Mechanics by G o l d s t e i n  [50].  What i s B in terms of B' B = 'B  X  ?  DB' [cos (0+n )cos\//-cost9sin ((j>+n)sin^]B' +  [ s i n (<p+n )cos\//+cost9cos (<p+n) s i n ^ ] B  y  + sin0sin\//B' 2  B  [-cos  (<j)+n) s i n ^ - s i n t 9 s i n (4>+n) s i n ^ ] B '  y + [ - s i n (<t>+n) sinv//+cos0cos (0+n) c o s ^ ] B  y  + sint9cos>//B'  sint9sin(0+n)B  B  x  - sinflcos(0+n)B  y  +  sin0B  z  z (A2.51)  138  One c a n t h e n s u b s t i t u t e  t h e s e v a l u e s f o r t h e components of B  i n t o e q u a t i o n A 2 . 1 7 , and g e t H_ axes.  i n terms of the  external  139  O • u a * o — u 0 c u • 1 * — c — — ^  (A2.52)  K  • C  c » c — • c • • m — o  n o i u —  —  —  c  N C *•» C D  • •  c  o -*« i> to  ~ a { • l  B  o a u " c • • •> c o • o • «— c — — a ei - ~ i  v  o  •  B  C •  « ~  U —  to •—• o •» •* O C 01 — — o  c  -» «  o c o + — • • a • o — u c »  •*  — c •  « N B  • — c i m «» f- o  •  •  —  o  — . N  c c n • • •» • • to •—• •—• o e n u o «> •o o c  a u r « .  • c «l  I  j- to  n O •K  c  ^4 • K m to  c -» 01 • • o V  »o  c * _u1 • -  c au > 1 e so 11 O f- u  o a u  m o o  «M  i  U  B  •»  «o o c u •> c »— » N c » n V) n « o •» O c tt> u u o> n n c i o *> u •  +  B  c •* C • — ~ «0 - » K n o c o  -» c «0 «~ c  ~  •  •  B  4  ~  _  —»  _ c — •» c • m c a — • o ~* f 10 ~- U •> • <t> c «• » t> n ~* B c o •» © — — u — o •>  M  •  > N  a  • M» B  10  —  O  C  U  — c  •  -r< a tn •—  «i  c  m  —  w- to -~ o a u %  K  •  •o o _ _ I — - c ~ — a to  CO —  o U ~ c • • —  c —> • • N • a — -» w c o -~ o »  » —  o o  a "  B  •*  « ao  — •0 o u  u -~ c — c •  • • —  o u «i i  i  •  c  1 40  There are three of these corresponding of n. and  to the three  Each must be solved s e p a r a t e l y , and  intensities calculated.  then be the sum  The  of a l l t h r e e .  i n t e r p r e t i n g experimental  the  experimental  To use  frequencies spectrum  task  for the e x i s t e n c e of computers.' The Hamiltonians  spectrum can be p r e d i c t e d .  6,  will  i t f o r the purpose of  data would be a daunting  n u m e r i c a l l y s o l v e d f o r a given  values  <j>, 4/, and  the  can  but  be  experimental  T h i s can be done as an  iterative  process, d r i v e n by the disagreement between the p r e d i c t i o n and  the o b s e r v a t i o n , e v e n t u a l l y r e s u l t i n g  \p which should s u c c e s s f u l l y p r e d i c t the  i n a set of 6, <j>,  experimental  spectrum.  c.  Time dependence of the p o l a r i z a t i o n The  for a single  site  time dependence of the p o l a r i z a t i o n can be w r i t t e n  as: <a  >(t) = u « t  M • uX  f  U£  (A2.53)  where the elements of M are given  by:  M (f) = M (0) 2 + rr L I <m I a, I n> I I <n I o. I m> I {cos (CJ N n<m f I mn 'mn f i  f i  1  1  1  11  n  = 0  6* mn  = arg <mlo- |n> r  The  d i r e c t i o n of i n i t i a l  'mn  f  mn 3  e  -  1  1  t + rj  1  ) - cosrj 'mn (A2.54)  i  mn  1  1  p o l a r i z a t i o n u^ and  the  direction  of o b s e r v a t i o n u^ are known i n the e x t e r n a l frame. consequence of t h i s e x p r e s s i o n  f o r <o  >(t) u  u_£ i s p e r p e n d i c u l a r  to  f  One  i s that even i f  }  0  141 u-, <a  >(t) * 0 provided that u. i s not along a p r i n c i p a l  axis. When H_ i s solved i n the iX>, |Y>, jZ>, |0> b a s i s , the four e i g e n f u n c t i o n s which r e s u l t are of the form \p  = c, |x> 1n  < U V  x a  +  1  *n  1 1  C  C  1m 1n  +  C  Y >  =  z'V  C  3nl  C  +  2m 2n  C  Z >  +  C  C  4nl  3m 3n  0 >  C  +  +  C |Y> 4n  C  (A2. 55)  O  C  i f m equals n. |X>  S  '  Im4n  = 1  = - i C 4 n |X>  ff  +  i f m does not equal n  3n >  n  2nl  = 0  n  K  y  =  C  (A2. 56)  C | Z > + C |0>  +  1n  - iC |Y> 3n  2n  iC |Z>  +  2 n  iC |0>  +  1 n  |X> + C |Y> - C | Z > - C |0>, so 1n  2n  4n  3n  (A2. 57)  * * * * 1 \b > = C C + C C + C C + C C n 1nr3n ^ r r i 4n. ^3m^1n 4iri2n "fc "fc "ft "fc <\b l a 1 \b > = C C . _ - i C _ C~ + i C , C , + i C . C . m' y *n 1m 4n 2m 3n 3m 2n 4m 1n <\b 1 a Y  m  x  1  1 Y  c  L  _ i  1  Y  1 m  c  <\i/ l a | Yii > = m z n Y  *  c  + c  *  ic  1c  c - c c - c c ^2nT 1n ^3nT 4n ^4nT3n  ^1nT2n  1  i  0r  (A2. 58)  In the case where B' = 0 the only <i//_ | oj | \p^> remaining a r e <^ \o |tf,> 3  =  x  <^2  1  0 „  1 =  |o„|^ > 3  |iK> = 1 = <i|/„ I a |\// > x  2  <tf«|a |^,> = i  = -<^, | a | K >  i  = -<i// | a | <// >  < ^  3  1 »  | ^ 2 >  Y  =  y  y  2  3  < ^ 2 k | ^ i > = 1 = < i h | a |\// > z  z  <^4IO l^3> z  All  T?  mn  =  2  -1 = «// | o _ | t f , >  a r e 0, and  3  (A2. 59)  1 42  M(t)  =  M  (0)  xx  +  -^cos^^t  M  yy  ( 0 )  +  f  ^  l  ^  M  zz  (0) + +  2  C O S w  gcosco^t  34  t  (A2.60)  Simplifications: in  the absence the ^  2  4  ,  the M f o r each  u>^ a n d t o ^ a r e c o m p a r a b l e  hyperfine they since  frequency  and hence  site  is  t h e most  interesting  oscillates,  result  part  t o t h e vacuum  not r e a d i l y  c a n be d r o p p e d f r o m f u r t h e r  which  M(t)  an a p p l i e d f i e l d ,  same;  since  The  of  observable,  consideration;  i s g o i n g t o be t h e  the M^(0) can a l s o  part  be d r o p p e d .  is:  =  1 2  C O S C J  13  0 0  0  t  5  c o s w  23  0  t  COSU  t  2 }2  (A2.61)  1 43  In the s p e c i a l case where the i n i t i a l p o l a r i z a t i o n and d i r e c t i o n of o b s e r v a t i o n c o i n c i d e , that  i s , u* =  and  where u, M are i n i n t e r n a l c o o r d i n a t e s , <o >(t) = u«M«u f o r one s i t e .  (A2.62)  u  d.  Spectrum f o r Three The observed spectrum  Sites f o r the case of three s i t e s  r e l a t e d by a t h r e e - f o l d symmetry a x i s can be c a l c u l a t e d f o r t h i s simple case.  Using the D's as d e f i n e d e a r l i e r , and  with u' the p o l a r i z a t i o n d i r e c t i o n  in external coordinates,  simply average the s p e c t r a from the three s i t e s .  That i s ,  <a >(t) = ^[D -U'-M-Do-u' + D -u'•M-D -u' 0  +  +  + D..u'-M-D.«u'] where D  0  = 5(0), D  (A2.63)  = D(+120)  +  and  p.  =  D(-120)  <a >(t) = u  {[cos(</>+n)cos<//-cos0sin(0+n)sin\//]x' n=0,*120,-120 + [ s i n (</>+n)cos\//+cos0cos(0+n) s i n ^ ] y ' 2  + sin0sin\//z ' } K  1  3  l  +  {[-cos(</)+n)sin\/'-cosesin(</>+n)cos\/']x  + [ - s i n (</»+n) sin^+cos0cos (<f>+n) c o s i l y  1  2  + sin0cos\//z ' } K 3 2  + {sin0sin(c/>+n)x' - sintfcos (</>+n )y' 2  + cosflz*} K where  1 2  (A2.64)  - -^-cosa>^ 31 K  23  =  K  12  =  Clearly,  22"  C O S W  COSW  23  12  t  t  f o l l o w i n g some a l g e b r a i c m a n i p u l a t i o n ,  (A2.65)  1 44 2  2  2  2  2  2  <a >(t) = l [ { | s i n 0 s i n ^ / ( 2 z ' - x ' - y ' ) + | ( x ' + y ' 2 ) } K  1 3  + {|sin t9cos i//(2z' - x ' - y ' ) + | ( x ' + y ' ) } K 2  2  2  2  2  2  2  2  2  2  2  2  2  + {|cos t9(2z' -x' -y' )+|(x' +y' )}K ] 12  2  2  But x' +y' +z'  2  2 3  (A2.66)  =1  2  2  2  2  <a >(t) = • -{[sin 0sin \!/(3z' -1 ) + 1 - z ' ] c o s c j t 4  1 3  2  2  2  2  + [sin 0cos \//(3z' -1 ) + 1-z' ] c o s c j t 23  2  2  2  + [ c o s 0 ( 3 z ' -1 ) + 1-z' ] c o s w t }  (A2.67)  12  Note: sinflsini// = COS7  , A. A  sinflcost// = c o s 7 , y y  cos0  = COS7  ,  (A2.68)  zz the  d i r e c t i o n c o s i n e s , the angle between the c o r r e s p o n d i n g  e x t e r n a l and i n t e r n a l axes. 2  <a >(t) = l { [ c o s 7 u  So  2  x x  2  , (3z' -1 ) + 1-z' ]coso> t  + [COST  13  2  2  yy  + [cos7  2  2  , (3z' -1 ) + 1-z' ] C O S C J 2  2 2  9  zo  t  2  ,(3z' -1)+1-z' ]cosw t} 1 2  (A2.69)  Note the l o s s of a l l <f> i n f o r m a t i o n and the l o s s of s i g n i n f o r m a t i o n f o r the d i r e c t i o n c o s i n e s .  T h i s means that the  s i t e cannot be determined from zero magnetic  f i e l d data  alone, and a l l the i n f o r m a t i o n which can be o b t a i n e d i s obtained from a s i n g l e o r i e n t a t i o n of the c r y s t a l . In  the presence of a magnetic  field,  a l l the angles  w i l l have an i n f l u e n c e and so they can, i n p r i n c i p l e , be determined, but at a cost of having a spectrum of from 3 to 9 lines.  I f the a p p l i e d magnetic  symmetry a x i s , t h i s s i m p l i f i e s back t o the 3 l i n e s  field  i s along the  the spectrum c o n s i d e r a b l y ,  found i n zero f i e l d . .  Unfortunately, for  the same symmetry reasons as f o r zero f i e l d , dependence i s again the  lost.  the <j>  I t seems that i n t h i s  instance,  information can only come from a complex spectrum.  146  I. A P P E N D I X 3; T H E S U P E R H Y P E R F I N E I N T E R A C T I O N H A M I L T O N I A N For a spin 1/2 nucleus in zero magnetic field, with an anisotropic hyperfine interaction, where the axes of the superhyperfine tensor coincide with the axes of the anisotropic hyperfine tensor, the Hamiltonian yields the matrix: -2A  xx  0  -iA  -2A •A  XX  0  N  iA zz 0  -iA  0  -A  N zz  N  zz  o  0  o  iA  - 2 A. yy  0  o  -2A  o  N iA "xx  yy  -A  zz  iA  yy  N  N "yy -iA  yy  -iA  N xx  -2A  zz  N xx  N xx  0 0 -2A  zz  J  (A3.1) where the basis is: | X a > , |X/3 •>, | Y a > , | Y/3 >, | Z a > , |Zj3 > and X , Y , and Z defined as for the anisotropic hyperfine interaction in Appendix 2, and a and 0 are the spin states of the nucleus. By rearrangement, the matrix can be made block-diagonal as follows: -2A XX •A  N  iA A 0  zz  ..N -iA zz -2A •A  yy  iA 0  N  N  N  -lA -2A  XX zz  XX  0  0  0  0  0  0  0  0  0  yy •A  yy  N  A  0  -2A  •A  XX 0 0  0 0  0 0  •A  N  iA A  zz N  yy  N  -iA -2A  zz yy  •A  N  iA  A  N  yy ..N -iA XX -2A zz  XX (A3.2) The basis is now | X a >, | Y a >, | Z/3 >, - | X/3 >, | Yj3 >, | Za >. so the spectrum will have three pairs of degenerate levels.  This has two identical blocks and  Abragam 61 A. Abragam, The Principles of Nuclear Magnetism, Oxford University Press, London 1961. Albert 84 E. Albert, A. Moslang, E. Recknagel, A. Weidinger, Hyp. Int. 17-19 (1984) 611. Andersson 74 L.O. Andersson, 18th Ampere Congress, Nottingham, (1974) 129. Arseneau 84 Donald J. Arseneau, David M . Garner, Masayoshi Senba, Donald G . Fleming, J. Phys. Chem. 88 (1984) 3688. Barsov 83 S.G. Barsov, V.G. Baryshevskii, A.L. Getalov, V.A. Gordeev, S.P. Kruglov, L A . Kuz'min, S.A. Kuten, S.M. Mikirtych yants, G.V. Shcherbakov, Leningrad Nuclear Physics Institute preprint #879, Leningrad, (1983) [Russian]. Baryshevsky 78 V.G. Baryshevsky, S.A. Kuten, Physics Letters 67A (1978) 355. Bates 79 J.B. Bates, J.C. Wang, P.A. Perkins, Phys. Rev. B 19 (1979) 4130. Beck 75 R. Beck, P.F. Meier, A. Schenck, Z. Physik B 22 (1975) 109. Beder 78 Doug Beder, Nuclear Physics A305 (1978) 411. Beder 80 Douglas S. Beder, Phys. Rev. B 21 (1980) 3861. Berger 64 M.J. Berger, S.M. Seltzer, NASA SP 3012 1964. Berger 66 M . l . Berger, S.M. Seltzer, N A S A SP 3036 1966. Bevington 69 P.R. Bevington, Data Reduction and Error Analysis for the Physical Sciences, McGraw-Hill, New York, 1969. Boekema 84 C. Boekema, Hyp. Int 17-19 (1984) 305. 147 Bonner 76 W.A. Bonner, M.A. Van Dort, M.R. Yearian, H.D. Zeman, G.C. Li, Israel J. Chem. 15 (1976/77) 89.  148  Brewer 73 J.H. Brewer, K . M . Crowe, F . N . Gygax, R.F. Johnson, B.D. Patterson, D . G . Fleming, A. Schenck, Phys. Rev. Letters 31 (1973) 143. Brewer 75 J.H. Brewer, K . M . Crowe, F . N . Gygax, A. Schenck, in Muon Physics vol III edited by Vernon W. Hughes and C.S. Wu, Academic Press 1975, 4. Brewer 79 Jess H . Brewer, Douglas S. Beder, David P. Spencer, Phys. Rev. Letters 42 (1979) 808. Brewer 79b J.H. Brewer, D.P. Spencer, Hyp. Int. 6 (1979) 181. Brewer 81 J.H. Brewer, D.P. Spencer, D . G . Fleming, J.A.R. Coope, Hyp. Int 8 (1981) 405. Brewer 81b J.H. Brewer, P.W. Percival, eds. Muon Spin Rotation Proceedings of the Second International Topical Meeting on Muon Spin Rotation, Hyp. Int. 8 (1981) 307-853. Brewer 82 J.H. Brewer, D . G . Fleming, P.W. Percival, from Fourier, Hadamard and Hilbert Transforms in Chemistry, ed. A . G . Marshall, Plenum Publishing Corporation, 1982, 345. Brewer 85 J.H. Brewer, Priv. Comm., 1985. Brown 80 J.A. Brown, S.A. Dodds, T.L. Estle, R . H . Heffner, M . Leon, D.A. Vanderwater, Solid State Communications, 33 (1980) 613. Carr 84 J. Carr, G . Gidal, B. Gobbi, A . Jodidio, C.J. Oram, K..A. Shinsky, H . M . Steiner, D. Stoker, M . Strovink, R.D. Tripp, Hyp. Int. 17-19 (1984) 887. Celio 84 M . Celio, P.F. Meier, Hyp. Int. 17-19 (1984) 435. Chen 84 Y . Chen, R. Gonzalez, K.L. Tsang, Phys. Rev. Letters, 53 (1984) 1077. Clawson 81 C.W. Clawson, E E Haller, K . M . Crowe, S.S. Rosenblum, J.H. Brewer, Hyp. Int. 8 (1981) 417. Clawson 82 Carl William Clawson, Ph.D. thesis, Lawrence Berkeley Laboratory, University of California, 1982. Cowley 64 R.A. Cowley, Phys. Rev. 134 (1964) A981.  149  de Haseth 82 J.A. de Haseth, from Fourier, Hadamard and Hilbert Transforms in Chemistry, ed. A . G . Marshall, Plenum Publishing Corporation, 1982, 387. Estle 84 T.L. Estle, S.L. Rudaz, E. Holzschuh, R.F. Kiefl, B.D. Patterson, W. Ktindig, K.W. Blazey, Hyp. Int 17-19 (1984) 623. Estle 85 T.L. Estle, Priv. Comm. 1985. Evans 55 R.D. Evans, The Atomic Nucleus (1964) 586. Fleming 79 D.G. Fleming, D . M . Garner, L.C. Vaz, D.C. Walker, J.H. Brewer, K . M . Crowe, p. 279 in Positronium and Muonium Chemistry, ed. H.J. Ache, American Chemical Society, Washington, D.C. 1979. Fleming 82 Donald G . Fleming, Randall J. Mikula, David M . Garner, Phys. Rev. A26 (1982) 2527. Fleming 83 Donald G . Fleming, Randall J. Mikula, Masayoshi Senba, David M . Garner, Donald J. Arseneau, Chem. Phys. 82 (1983) 75. Fleming 84 D.G. Fleming, D.J. Arseneau, D . M . Garner, M . Senba, R.J. Mikula, Hyp. Int 17-19 (1984) 655. Friedman 57 J.I. Friedman, V.L. Telegdi, Phys. Rev. 105 (1957) 1681. Frondel 62 C. Frondel, The System of Mineralogy, Volume III, Silica Minerals 7th edition John Wiley & Sons, Inc. 1962. Garay 74 A.S. Garay, L. Keszthelyi, I. Demeter, P. Hrasko, Nature 250 (1974) 332. Garner 78 D. M . Garner, D . G . Fleming, J.H. Brewer, Chem. Phys. Letters 55 (1978) 163. Garner 79 D . M . Garner, Ph.D. Thesis, University of British Columbia 1979. Garwin 5750 Goldstein R.L. Garwin, LClassical . M . Lederman, M . Weinrich, Phys. Rev.1950,107. 105 (1957) 1415. H . Goldstein, Mechanics, Addison-Wesley,  1  Gonzalez 81 R. Gonzalez, Y . Chen, M . Mostoller, Phys. Rev. B 24 (1981) 6862. Goodrum 72 J.W. Goodrum, Journal of Crystal Growth 13/14 (1972) 604. Goshen 83 S. Goshen, M . Friedman, R. Thieberger, J.A. Weil, J. Chem. Phys. 79 (1983) 4363. Graf 79 H . Graf, E. Holzschuh, E. Recknagel, A. Weidinger, T. Wichert, Hyperfine Interactions 6 (1979) 177. Gygax 79 F.N. Gygax, W. Kundig, P.F. Meier, eds. Muon Spin Rotation Proceedings of the First International Topical Meeting on Muon Spin Rotation, Hyp. Int. 6 (1979) 1-450. Gygax 84 F.N. Gygax, A. Hintermann, W. Rtluegg, A. Schenck, W. Studer, A.J. van der Wal, Hyp. Int. 17-19 (1984) 177. Harshman 84 D. R. Harshman, R. Keitel, M . Senba, E.J. Ansaldo, J.H. Brewer, Hyp. Int. 17-19 (1984) 557. Heffner 84 Robert H . Heffner, Donald G . Fleming, Physics Today 37 (1984) 38. Henderson 72 B. Henderson, H.T. Tohver, Phys. Status Solidi 51 (1972) 761. Holzschuh 81 E. Holzschuh, W. Kundig, B.D. Patterson, Helvetica Physica Acta 54 (1981) 552. Holzschuh 82 E. Holzschuh, W. Kundig, P.F. Meier, B.D. Patterson, J.P.F. Sellschop, M.C. Stemmet, H . Appel, Phys. Rev. A25 (1982) 1272. Hughes 60 V.W. Hughes, D.W. McColm, K. Ziock, R. Prepost, Phys. Rev. Let. 5 (1960) 63. Isoya 83 J. Isoya, J.A. Weil, P.H. Davis, Journal of Physics and Chemistry of Solids 44 (1983) 335. Ivanter 72 I. G . Ivanter, V.P. Smilga, Sov. Phys. JETP 34 (1972) 1167. Jackson 62 J.D. Jackson, Classical Electrodynamics, John Wiley and Sons, 1962,412. James 71 F. James, M . Roos, MINUIT, C E R N Computer 7600 Interim Programme Library (1971).  151  Jean 79 Yan-Ching Jean, Hans Ache, in D.C. Walker, Origins of Optical Activity in Nature, Elsevier Scientific Publishing Company, Amsterdam, Netherlands, (1979) 67. Kadano 84 R. Kadano, J. Imazato, K. Nishiyama, K. Nagamine, T. Yamazaki, D . Richter, J.M. Welter Hyp. Int. 17-19 (1984) 109. Karaseva 83 L.G. Karaseva, T.A. Karpukhina, B.V. Spitsyn, Russian Journal of Physical Chemistry 57(2) (1983) 302. Kiefl R.F. Kiefl, Private Communication. Kiefl 84 R.F. Kiefl, E Holzschuh, H . Keller, W. Kundig, P.F. Meier, B.D. Patterson, J.W. Schneider, K.W. Blazey, S.L. Rudaz, A.B. Denison, Phys. Rev. Lett. 53 (1984) 90. Kiefl 85 R.F. Kiefl, J.W. Schneider, H . Keller, W. Kundig, W. Odermatt, B.D. Patterson, K.W. Blazey, T . L Estle, S.L Rudaz, submitted to Rapid Communications, Physical Review B (1985). Kinase 82 W. Kinase, N . Ohnishi, Phys. Letters 90A (1982) 101. Lee 56 T.D. Lee and C.N. Yang, Phys. Rev. 104 (1956) 254. Leiss 57 J.F. Leiss, S. Penner, C.S. Robinson, Physical Review 107 (1957) 1544. Le Page 80 Y. Le Page, L D . Calvert, E J . Gabe, J. Phys. Chem. Solids 41 (1980) 721. Mogensen 74 O . E Mogensen, J. Chem. Phys. 60 (1974) 998. Myasishcheva 68 G . G . Myasishcheva, Yu.V. Obukhov, V.S. Roganov, V.G. Firsov, Sov. Phys.-JETP 26 (1968) 463. Nagamiya 81 S. Nagamiya, K. Nagamine, O. Hashimoto, T. Yamazaki, Phys. Rev. Lett. 35 308. Norton 76 Norton, Beer, J. Opt Soc. Am. 66 (1976) 259. Noakes 85 D. Noakes, Priv. Comm. 1985. Nuttall 81 R.H. Nuttall, J.A. Weil, Can. J. Phys. 59 (1981) 1696.  152 Oram 84 C.J. Oram, J.M. Bailey, P.W. Schmor, C.A. Fry, R.F. Kiefl, J.B. Warren, G . M . Marshall, A. Olin, Phys. Rev. Lett. 52 (1984) 910. Pake 73 G.E. Pake, T.L. Estle, The Physical Principles of Electron Paramagnetic Resonance, W.A. Benjamin Inc., Reading, Massachusetts 1973. Papp 75 N . Papp, K.P. Lee, J. Mag. Res. 19 (1975) 245. Patterson 78 B.D. Patterson, A. Hintermann, W. Kundig, P.F. Meier, F. Waldner, H . Graf, E Recknagel, A. Weidinger, T. Wichert, Phys. Rev. Letters 40 (1978) 1347. Patterson 84 B.D. Patterson, Hyp. Int. 1 7 - 1 9 (1984) 517. Patterson 84b B.D. Patterson, E Holzschuh, W. Kundig, P.F. Meier, W. Odermatt, J.P.F. Sellschop, M . C . Stemmet, Hyperfine Interactions 1 7 - 1 9 (1984) 605. Percival 76 P.W. Percival, H . Fischer, M . Camani, W. Ruegg, A. Schenck, H . Schilling, H . Graf, Chem. Phys. Lett. 3 9 (1976) 333. Percival 78 P.W. Percival, E Roduner, H . Fischer, Chem. Phys. 3 2 (1978) 353. Percival 79 Radiochimica Acta, 26 (1979) 1. Percival 82 Paul W. Percival, Jean-Claude Brodovitch, Kenneth E Newman, David P. Spencer, Chem. Phys. Letters 9 3 (1982) 366. Percival 84 P.W. Percival, J.C. Brodovitch, K . E Newman, Hyp. Int. 1 7 - 1 9 (1984) 721. Percival 85 P.W. Percival, K . M . Adamson-Sharpe, J.C. Brodovitch, S.K. Leung, K . E Newman, Chemical Physics 9 5 (1985) 321. Perlson 74 B.D. Perlson, J.A. Weil, J. Mag. Res. 15 (1974) 594. Roduner 84 E Roduner, Hyp. Int. 1 7 - 1 9 (1984) 785. Sachs 75 A . M . Sachs, A . Sirlin in Muon Physics vol II edited by Vernon W. Hughes and C.S. Wu, ' Academic Press, 1975.  153  Sahoo 85 N . Sahoo, K.C. Mishra, T.P. Das, Private Communication, 1985. Sakudo 71 T. Sakudo, H . Unoki, Phys. Rev. Letters 26 (1971) 851. Schenck 76 A. Schenck p. 159 in Nuclear and Particle Physics at Low Energies, Ed. J.B. Warren, Plenum Publishing Company, New York, 1976. Schenck 81 A. Schenck, Helvetica Physica Acta 54 (1981) 471. Schumacher 80 Rolf Schumacher, Ber. Bunsenges. Phys. Chem. 84 (1980) 125. Senba 84 M . Senba, D . M . Garner, D.J. Arseneau, D.G. Fleming, Hyp. I n l 17-19 (1984) 703. Smith 64 G.S. Smith, P.B. Isaacs, Acta Cryst. 17 (1964) 842. Spencer 79 David P. Spencer, Donald G . Fleming, Jesse H . Brewer, Randy J. Mikula, in D.C. Walker, Origins of Optical Activity in Nature, Elsevier Scientific Publishing Company, Amsterdam, Netherlands, 1979. Stadlbauer 84 J.M. Stadlbauer, B.W. Ng, Y . C . Jean, Y. Ito, D.C. Walker, Hyp. Int. 17-19 (1984) 715. Stoneham 83 A . M . Stoneham, Priv. Comm. 1983 Swanson 58 R.A. Swanson, Phys. Rev. 112 (1958) 581. Trower 65 W. Peter Trower, University of California Radiation Laboratory report #2426 (1965) 15. Turner 85 R.E. Turner, Priv. Comm. Vester 57 F. Vester, Seminar at Yale University, Feb. 7,1957. Walker 79 D.C. Walker, Origins of Optical Activity in Nature, Elsevier Scientific Publishing Company, Amsterdam, Netherlands, 1979. Walker 81 D.C. Walker, Hyp. Int 8(1981) 329.  154  Willemsen 76 Herman W. Willemsen, Robin L. Armstrong, Peter P.M. Meincke, Phvs. Rev. B14 (1976) 3144. Weil 71 J.A. Weil, J. Chem. Phys. 55 (1971) 4685. Weil 75 J.A. Weil, Rad. Eff. 25 (1975) 261. Weil 81 J.A. Weil, Hyperfine Interactions 8 (1981) 371. Weil 81b J.A. Weil, Can. J. Phys. (1981) 59 (1981) 841. Weil 84 J.A. Weil, Phys. Chem. Minerals (1984) 10 (1984) 149. Wu 57 C.S. Wu, E. Ambler, R.W. Hayward, D . D . Hoppes, R.P. Hudson, Phys. Rev. 105 (1957) 1413. Yamazaki 84 T. Yamazaki, K. Nagamine, eds. Muon Spin Rotation and Associated Problems, Hyp. Int. 17-19 (1984) 1-1034. Yasunaga 70 Hitosi Yasunaga, Journal of the Physical Society of Japan 28 (1970) 454. Young 62 R.A. Young, B. Post, Acta CrysL 15 (1962) 337.  


Citation Scheme:


Citations by CSL (citeproc-js)

Usage Statistics



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"
                            async >
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:


Related Items