T H E O R E T I C A L C A L C U L A T I O N O F M U O N SITE IN YBa Cu Ojr 2 By Qiang Li B. Sc. (Nuclear Physics)Peking University A THESIS SUBMITTED T H E INPARTIAL REQUIREMENTS MASTER FULFILLMENT O F F O RT H E D E G R E E O F O F SCIENCE in T H E F A C U L T Y O FG R A D U A T E STUDIES PHYSICS We accept this thesis as conforming to the required standard T H E UNIVERSITY O FB R I T I S H June 1990 © Qiang Li, 1990 COLUMBIA 3 In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of British Columbia, I agree that the Library shall make 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. Physics The University of British Columbia 6224 Agricultural Road Vancouver, Canada V6T 1W5 Date: / / Abstract The muon is a useful probe of magnetic fields in superconductors, but knowing the field seen by the muon is often of limited value until we know where the muon is in the crystal lattice. In this thesis I employ two independent theoretical methods to search for candidate muon sites: the potential energyfieldmethod, which seeks the minimum of the electrostatic potential of the and the magnetic dipolarfieldmethod, which compares the calculated magnetic field (due to the host electronic, atomic or nuclear dipolar fields) with the observed local fields at the muon. ii Table of Contents Abstract ii List of Tables v List of Figures vi acknowledgement 1 2 3 4 viii Introduction 1 1.1 History of Muon Site Determination 1 1.2 Space Group 3 1.3 Muon Spin Resonance/Rotation/Relaxation (uSTZ) 4 Methods of M u o n Site Determination 7 2.1 Electrostatic Potentials 7 2.2 Dipolar Magnetic Fields 13 2.3 Conclusion 15 Explicit Calculations of M u o n Sites In Y B a C u 0 2 3 6 16 3.1 Potential Maps of A F M 16 3.2 Dipolar Magnetic Field Calculations 36 3.3 Comparison With Experimental Data 39 Calculations on Y B a C u 0 2 4.1 3 43 7 Muon Site Decision Experiment 43 iii 4.2 Potential Map of Superconductor Y B a C u 0 4.3 Dipolar Magnetic Field 46 4.4 Conclusion 47 2 3 7 45 Bibliography 50 A 53 Programs iv List of Tables 1.1 Transformations pmmm of Y B a C u 0 7 1.2 Transformations pl23 of Y B a C u 0 3.3 Properties and coordinates of ions 16 3.4 The magnetic field at z=0.03c 38 3.5 The magnetic field at z = 0.19c 38 3.6 Muon candidate site 39 4.7 The positions and properties of Y B a C u 0 4.8 The nuclear magnetic moments of ions in Y B a C u 0 a ; 2 2 4 3 3 • 6 2 3 45 7 2 v 5 3 47 List of Figures 2.1 Potential Vi(r) vs R 10 3.2 Two dimensional potential map at Z=0.02c 19 3.3 Two dimensional potential map at Z=0.03c 20 3.4 Two dimensional potential map at Z=0.05c 21 3.5 Two dimensional potential map at Z=0.07c 22 3.6 Two dimensional potential map at Z=0.10c 23 3.7 Two dimensional potential map at Z=0.15c 24 3.8 Two dimensional potential map at Z=0.17c 25 3.9 Two dimensional potential map at Z=0.19c 26 3.10 Two dimensional potential map at Z=0.20c 27 3.11 Two dimensional potential map at Z=0.26c 28 3.12 Two dimensional potential map at Z=0.30c 29 3.13 Two dimensional potential map at Z=0.35c . 30 3.14 Two dimensional potential map at Z=0.38c 31 3.15 Two dimensional potential map at Z=0.40c 32 3.16 Two dimensional potential map at Z=0.42c 33 3.17 Two dimensional potential map at Z=0.45c 34 3.18 The contour potential map at Z=0.03c 35 3.19 The contour potential map of Y B a C u 0 at Z=0.19c 37 3.20 Contour plot of magnetic field at 0.03c 40 3.21 Contour plot of magnetic field at 0.19c 41 2 vi 3 6 3.22 The crystal structure of antiferromagnetic YBa Cu 06 42 4.23 The contour plot of the potential map at Z=0.19c 48 2 3 4.24 The crystal structure of superconducting Y B a C u 0 7 2 vii 3 49 acknowledgement It has been a pleasure to work in this group from which I have benefited and enjoyed so much. I would like to express my thanks to Rob Kiefl who has been not only helpful to my thesis work but also a good teacher for a broad area in /iSR. I would like to say that I have very much enjoyed working with Tanya Riseman, a helpful colleague and an interesting person to talk with. All I have achieved from my work should be credited to my supervisor Jess Brewer. I am grateful for what he has done, although I have a little trouble with his American accent. The data employed in this thesis were provided by several people; I very much appreciate those involved. viii Chapter 1 Introduction 1.1 History of Muon Site Determination The first /x site determination was performed in a single crystal of gypsum (CaS0 • + 4 2H 0), in which the magnetic field due to the adjacent protons at fixed locations is 2 well known. If the fi were mobile in the crystal, the dipolarfieldsof the protons would + simply cause relaxation of the muon spin. However, if the muon replaces protons at fixed lattice sites, the Muon Spin Resonance (pSR) results should show exactly the same behavior as the proton NMR. The experiment showed the split-frequency muon precession expected for muons occupying proton sites, and thus gave the first positive identification of the /x site in a crystal. [1] + During the following years the fiSH technique has been applied to metals, oxides, magnetic materials and superconductors to explore the internal magnetic and electronic structure. But without knowing the exact locations of the fx in the samples the + information often has limited value. The local magnetic field at a stopped fi + site was measured in a single crystal of ferromagnetic cobalt as a function of temperature between 4 and HOOK by Graf. [2] The data are consistent with a relatively smooth temperature dependence of the hyperfine field only if it is assumed that the fi + is at an octahedral interstitial site; the calculated local field at the tetrahedral site would imply a discontinuous hyperfine field. These facts lead to the conclusion that octahedral interstitial sites are preferred 1 Chapter 1. Introduction 2 by the p . + The longlived metastable states of the u in the magnetic oxide a—Fe203 have also + been studied. [3] Three separate frequencies were seen. This indicates that the muons are localized at three different local energy minima in the unit cell. Other muon site determinations have been performed on Cu, [5] Fe, [6] [7] and Si. [8] Lattice distortion plays a key role in determining the muon site in Alkali fluorides and semiconductors. One of the most interesting magnetic phenomena in a superconductor is the well known Meissner effect. When a type I superconductor is cooled in a magnetic field below its transition temperature, the magnetic flux is abruptly expelled, except for a thin surface layer where the field decays roughly exponentially on the scale of the London penetration depth Aj,. In a type II superconductor the magnetic flux penetrates into the interior of the sample in an array of "vortices" each carrying the flux quantum —• —* ip when the magnetic field is between H \ and H . As a result, the internal magnetic 0 c c2 field in the superconductor is inhomogeneous. The pSTZ technique can measure the mean magnetic field at the muon site (through the p dispersion (through the p + + precession frequency) and its relaxation rate). Experiments on superconducting GdBa Cu307 were performed in the temperature 2 range from 30mK to 130K by A. Golnik. [9] In zero external field two different frequencies were distinguished, which correspond to internal magnetic fields of 33 and 52 mT, respectively. This shows that muons find two magnetically inequivalent stopping sites in the sample. Two types of muon stopping sites are in principle possible: interstitial sites and oxygen vacancy sites. So far there is no solid evidence to show the p TRIUMF [46] group predicts that p + + site in YBa Cu 07. 2 3 Our own stops primarily in a site close to the C u 0 plane. 2 The PSI [48] and Tokyo [47] groups claim that the p stopping site is close to an oxygen + Chapter 1. Introduction 3 on the CuO chain. The only point of agreement is that the muon site is l A away from some oxygen ion. So far no rigorous theoretical prediction has been made, due to the difficulties involved in evaluating electrostatic potentials and magnetic fields inside the superconductors. The current theoretical method is to apply the potential equation which best generates the fi sites in YBa Cu306 to the case of YBa Cu 07. During the round-table + 2 2 3 discussion at the Vth International Conference on pSTl in Oxford (1990), it became clear that one reason for the differences of opinion about the muon site in Y B a C u 0 2 3 6 is the uncertainty regarding the orientation of Cu magnetic moments which should lie in the Cu-0 plane oriented along an a or b axis. The initial approach taken to locating the muon site focussed on the antiferromagnetic (AFM) material YBa Cu306, which has a magnetic structure consisting of strong 2 nearest-neighbor correlated spins in adjacent copper planes which are aligned antiferromagnetically. This special case offers some hope of calculating the internal magnetic field at potential muon sites theoretically. The simulated results can be verified by comparing them with the experimental data. The first study of local fields at muon sites in antiferromagnetic YBa Cu 06+r 2 3 was done by Nishida. [14] The experimental data show that there are at least two independent fj, sites. + 1.2 Space Group A perfect crystal is constructed by an infinite regular repetition in space of identical structural units or building blocks. A space group, which consists of a Bravis lattice and a basis, displays all the symmetries of a crystal. Through the symmetry of the space group, one point generates several equivalent points at symmetric positions in a Chapter 1. Introduction 4 X -X -X X X -X -X X Y -Y -Y Y -Y Y Y -Y Z -Z Z -Z -Z Z -Z Z. Table 1.1: Transformations pmmm of YBa Cu 0*7 2 3 unit cell. Therefore, a space group fully represents a crystal structure. [10] [36] [13] Neutron-diffraction experiments [12] have shown that YBa Cu 07 has a pmmm 2 3 structure and YBa Cu 06 has similar structure but with 4 transformations instead of 2 8 for YBa Cu 07. 2 3 3 i.e., every input ion site will generate 8 equivalent ion sites for YBa Cu C>7 and 4 equivalent ion sites for YBa Cu 0*6, due to its antiferromagnetic 2 3 2 3 structure. In summary, the transformations of YBa Cu 07 are pmmm and orthorhombic. Its 2 3 startpoint, the origin of the cluster, is (0.0,0.0,0.0) which can be specified by users. The transformations of Y B a C u 0 are pl23 which is termed for convenience and tetrago2 3 6 nal. Its startpoint is (0.0,0.0,0.0). Through such transformations, (1,1,1) represents a set of 4 sites at (1,1,1), (-1, -1,1), (1, - 1 , -1) and (-1,1, -1) for Y B a C u 0 and a set 2 3 6 of8sitesat(l,l,l),(-l,-l,-l),(-l,-l,l),(l,l,-l),(l,-l,-l),(-l,l,l),(-l,l,-l) and (1, -1,1) for Y B a C u 0 . 2 1.3 3 7 M u o n Spin Resonance/Rotation/Relaxation (pSTt) The muon possesses a magnetic moment and a spin of ~; hence it exhibits Larmor precession in a tranverse magnetic field. After stopping inside the target, the implanted Chapter 1. Introduction 5 X -X X -X Y -Y -Y Y Z Z -Z -Z Table 1.2: Transformations pl23 of Y B a C u 0 2 3 6 highly spin polarized positive muons, precess in the local field and decay (fi —* e -++ e + + with a mean lifetime of 2.198/is. The decay positrons are emitted in directions v correlated with the direction of the muon spin. Inhomogenenities of the local field will cause loss of muon spin ensemble polarization, termed depolarization, by dephasing. In this case information about the internal or external field distribution can be derived from the experimental distribution of precession frequencies which are proportional to the local field. In general, the precession frequency u> and depolarization rate a are given by f^H w = (1.1) where M = ((H — (H)) }, u = {H) and 7 is the gyromagnetic ration of the muon. 2 2 M Specific equations have to be derived for specific cases. In this thesis two experimental techniques are referred to zero applied field (ZF) fiSH and Level Crossing Resonance (LCR). In ZF-fiS7l the fi + undergoes Larmor —* precession around the local magnetic field. For a local magnetic field given by H = (H ,H ,H ) x y z the muon spin initialially polarized along the z direction evolves as P,(i) = ^ + = cos 9 + sin^Bcos^nHt) 2 ^ cos( Ht) Hx2 Hy2 yfl (1.2) 6 Chapter 1. Introduction where 6 is the angle between H and the z direction. The internalfieldcomponents are assumed to have a gaussian distribution where A = 7 (if,- ). The average of the polarization in the z direction is 2 2 2 P (t) = J J J piH^p^piH^cos^ z + ain^cos^H^H Evaluation of this equation leads to the well-known Kubo-Toyabe formula for the evolution of a static spin in zero field: P*(t) = | + f (1 " 7 „ A V ) « p ( - ^ A * ) 2 2 2 2 The application of LCR was first suggested by Abragam [15]; it quickly became one of the most valuable pSTZ techniques. The level-crossing occurs when the field strength is adjusted to match the energy splitting of the muon's Zeeman levels to those of the magnetic levels of the neighboring host nuclei (in this study, 0 ) , which are determined 17 mainly by their quadrupolar interaction energies. Details of the LCR technique can be found in references [16] [17] [18] [19]. Chapter 2 Methods of Muon Site Determination In this chapter two theoretical approaches are introduced to search for candidate muon sites: electrostatic potential mapping and dipolar magneticfieldcalculation. 2.1 Electrostatic Potentials The principal interactions of the u in ionic lattices are the static long-range Coulomb + force and the overlap force, both of which are envisioned mainly as simple two-body forces; electronic polarization effects ( in which an ion may be polarized by the resultant electric field due to other ions in the crystal) are considered to be less important. The u interaction potential is the sum of all two-body interactions between the muon and + other ions. The concept of an additive two-body interaction is only an approximation to reality; the wave mechanical calculations made by Lowdin [22] showed that an appreciable fraction of the lattice cohesive energy cannot be represented in terms of two-body interactions. Moreover, the effect of electronic polarization cannot be represented by a two-body interaction. However, the polarization effects are largely suppressed by the high symmetry of the structure so that the approximation is adequate for most purposes. The Coulomb potential can be employed to describe the potentials of ions at distances much larger than their ionic radii; in insulating YBa2Cu30 (as in most ionic x oxides), this is a fair description of all the positive ions. In metallic YBa Cu 07, atten2 3 tion should be paid to the screening effect of conduction electrons on the field due to 7 Chapter 2. Methods of Muon Site Determination 8 ions fax away from the muon sites. For nearby ions, in this instance O , the calculation 2 - is more difficult. I resort to a familiar approximation, namely the Morse-like potential, which is frequently used for quantum mechanical calculations in solid state physics to calculate the potentials between / i and 0 : + V(r) = i where A = Ze and r is the /x —O + + 2- 2 - e x p ( - ^ ) - L # , B (2.3) distance. The first term arises from the Coulomb interaction between fi and 0 ~. The second term is the Morse-like potential, which + 2 represents the major part of the 0 potential. The semi-empirical constants B 2 - and p are determined by experiment, and r is the equilibrium fi —0 " bond lengh. The + 2 0 third term arises from the interaction of /x with the induced electric dipole moments + by the muon on the O 2 - ions, p. The induced moment was determined from the following procedure. A unit dipole po was placed on a test O 2 - ion giving a dipole moment of pi = poi on this ion, where i is a unit vector along x. The electric field at the test ion produced by these other dipole moments is given by E j = (Ei , E , E ). x iy This procedure was repeated with lz unit dipole moments along the y and z axes of the test ion, p = poj and p = p k; the 2 3 0 corresponding electric fields calculated at the test ion are given by E = ( £ 2 1 , E , 2 and E = (E ,Ez ,Ez ). 3 3x y z The net dipole moment at the O 2y E) 2z ion can be written 2 - This moment is induced by the electric field due to the point charges of all the other ions Eo plus the electric field due to induced dipole moments on all the other O ions. Thus, p = a [ E + (l/poXPxE! + p E + p,E )], 0 0 v 2 3 where c*o is the polarizability of 0 " in Y B a C u C v The last two equations can be 2 2 3 Chapter 2. Methods of Muon Site Determination 9 written as a set of three simultaneous equations: p x p v P x = CiolE + (l/p )(p E = a [E = ao[E -r(l/p )(p E 0x 0 x v x lx + (l/po)(p E x 0x If E and (l/po)(p Ei + p E 0 0x 0 0 x + p E -rp E )], (2.4) + p E +p E )], (2.5) y lv v 2x 2y + pE l2 v z z 3x 3y + p E )}. 2i z (2.6) 3i + p E ) could be derived from a lattice-sum calculation, 2 z 3 then these equations would be solved for p , p„, and p . x z In order to fix the parameters in Eq.(2.2), I simplified the model. First, no consideration was given to any screening effects, which would introduce a screening factor to the first term. However, one reference is given here for those interested in solving screening effects. Given the electron density, the charge distribution in the p —0 ~ + 2 could be estimated as for OH", [20] and the screening factor could be fixed. Second, the dipolar part is suppressed by the high symmetry of the structure, so it was neglected. Recent ZF-pSTZ and wTF experiments on 0-doped YBa Cu 07 [21] place the 17 2 3 p approximately 1 A away from some oxygen ion; LCK-pS1Z measurements suggest + specifically the 0(2, 3) oxygen in the C u 0 planes. Previous ZF-pSTZ experiments 2 [23] [24] on A F M Y B a C u 0 provide some constraints on the local field. Based on 2 3 6 this foreknowledge, and assuming that the muon's site preference will not be much affected by oxygen deficiency in YBa Cu 06, I narrowed my search down to a small 2 3 ensemble of candidate sites within one unit cell. This allowed me to fix certain necessary parameters in the semi-empirical equations with which I calculated the crystal potential in superconducting YBa Cu 07 to better locate the muon sites, 2 3 V;.( ) = 2 + !2l r (19.4-25r[A]) r r where Q is the valence of the ion and r is the ion's hard core radius. The first term 0 0 Chapter 2. Methods of Muon Site Determination 10 Figure 2.1: Potential Vj-(r) vs R is the Coulomb potential and the second term is the repulsive overlap potential. This choice of parameters ensures that the muon-oxygen bond length is 1 A. The total potential at the /x is in the form + v; = £ v ; ( r ) ot t=i The summation has been carried out in two regions defined by a spherical boundary withfiniteradius. Upon the asssumption of zero net charge in the unit cell, the coulomb contribution of ions outside this sphere goes to zero as long as the distance from that unit cell to the fi + site is long enough compared with its dimensions, because this unit cell can be considered as a point with zero net charge. The contributions from included ions are simply summed. The exponential term is convergent. Therefore the Chapter 2. Methods of Muon Site Determination 11 convergence of this sum is ensured by assuming that the net charge of the unit cell is zero. Another general method for calculating the interionic potentials, based on the density functional approximation to the energy of an electron gas, is introduced here but not used in this thesis because of the difficulties of finding the electro densities in YBajCuaCv Gordon and Kim [26] assumed that the interatomic interactions between all charges must be evaluated from the additive atomic densities. For a pair of ions, AB, the total density PAB{ ) R p (r) and PB^)', A 1S assumed to be the sum of the separated ion densities the total energy is therefore given by EAB = E[p ] = E[p + p ] AB A B so that the interaction energy can be written as E INT = E[p + ] A PB - E[p ] A - E[p ) B This is the essence of the approximation to the interionic potential. The energy of an isolated, closed-shell ion is given by E[p) = C„ J[p(r)f dr + C, J M r ) ] * ' * - Z f p(r)/r<ir 3 \j 3 j y r ^ * ' +1*[*•)]*•)* in which Z is the nuclear charge, p is the ion density and C = (3/10)(3TT ) / 2 k 2 3 and C = e -(d/tyS/v) ' . 1 3 The first four terms represent the kinetic, exchange and coulomb energies respectively, while the final term is an estimate of the electron pair correlation energy. The correlation energy density e [p(r)] is simply an interpolation between the high and low density c limits. Chapter 2. Methods of Muon Site Determination 12 The Coulombic interaction of this pair is given by V. = Z Z /R A I - Z B B flHh -Z A TIB J I^ d r J riA + / / '" 2 J ( r i ) r J '" 'Wr ( r ; 1 2 (2.7) The electron-gas contribution to the interatomic interaction is V = g fUpA(r) + p (r)]Eo\p (r) B - (r)E [pA(r)] PA The final form for Ei EUR) = Z Z /R A nt - Z B - G p {r)] B {r)E [p {r))}dr PB G B is thus, I Hdnljr, B J _ Z r J lB PB -p (r)E [pA(r)] G I SSiDldr, A + J {[PA{T) + (r)]E [p (r) A + A G + A - r 1A + / / " ^ " J J r B ^ d r , d r 2 12 p (r)] B PB(r)E [p (r))}dr G B in which EGW)] and r , r \ , r \ A B and r i2 = C [p(r)} 2/3 k + C [p(r)} e 1/3 + *Mr)} are functions of the internuclear separation, R. Here only nearest-neighbour interactions axe included in the lattice summations. The next-nearest-neighbour interactions lead to an overestimate of the cohensive energy. Rae [27] pointed out that the exchange energy in the equation includes a selfenergy contribution which, though negligible for an infinite electron gas, is significant for a small, finite number of electrons and leads to an overestimate of the exchange energy in the interaction of two light atoms or ions. A correction has been made by replacing the exchange term by a modified contribution, V (GKR) e = V (GK)[l e - 8/36 + 2S + 1/3S } 2 4 13 Chapter 2. Methods of Muon Site Determination in which 6 is a solution of (4iV) = 6 (1 - 9/86 + 1/46 ) -1 3 3 and N is the number of electrons. It has been concluded that for solids such as the alkali and alkaline-earth halides and the alkaline-earth oxides, which are largely ionic, the modified electron-gas approximation is a reliable non-empirical method for the calculation of interionic potentials. [28] [29] [30] One different potential assumption is also used by W.K. Dawson et al. [32] They assume that the muon-oxygen bonding potential has the form v = (t + where r is the / x — O + 2- 1 4 4 e y distance in A and a, b, and c are semiempirical constants (a = +0.359,6 = 0.0531, c = +8.32). The calculations are carried out within a finite sphere of radius on the order of 50A. Convergence of the lattice sum is ensured by assuming that the net charge of the unit cell is zero. 2.2 Dipolar Magnetic Fields In the zero field (ZF) case, [31] the local magnetic field at the / i is usually split up + into two contributions: H = Hkf + Hdip where the dipolarfieldcontribution is from the O 2 - interaction. The hyperfine field H/i/ results from the contact interaction between the fi and any polarized electrons + density at the muon site. When covalent bond effects occur as in muon-oxygen bond formation, Hhj is mostly due to these effects. Generally HH/ = 3f[n (r*/i) — n~(r )], + M 14 Chapter 2. Methods of Muon Site Determination where n (r ) is the density of spin-up electrons at the / i and n~(r ) the corresponding + + M M spin-down density. In the case of a—Fe20*3 the hyperfine effect is due to electron spin transfer into unoccupied metal 3d-orbitals, which causes spin polarization of the oxygen 2p-orbitals. These 2p-orbitals overlap with the muon s-orbitals, resulting in a nonzero spin density at the muon site. In order to estimate H/,/ for a muon participating in a "muoxyl bridge", Sawatzky [33] pictures a simple arrangment of one oxygen ion surrounded by metal ions; only those metal ions which form a direct link with the muoxyl bridge were taken into account. For this simple structure H ^ / is assumed to be the sum of the contributions of the linkages: H^/ = C • - A )cos d + A „]S,2 2 1t 2 i where S, is the unit vector of the magnetic moment of the metal ion, C is a constant, 9 is the angle between the metal-oxygen and muon-oxygen directions and A „ and A „ 2 2 are the magnitudes of the spin polarization of oxygen orbitals with er and 7r symmetry. The spin density at the muon is proportional to (A 2 a — A )cos 6i + A . 2 2 2 ir T Here I neglect Hhf in the case of antiferromagnetic YBa Cu 06, because of too 2 3 small electron spin density at muon site, and assume that the field at the muon arises mainly from the dipolar magnetic fields of surrounding atoms. The dipolar field is H , = (^j where H ; is the field generated by the i [3(/i, • r\)f, - /«,•] t h magnetic moment, r,- is the vector distance from that magnetic moment to the muon and fa is a unit vector in the direction of the magnetic moment. The resultant magnetic field at any point is the sum of such contributions over the lattice. It is calculated by summing explicitly over the contributions from lattice points within a certain radius R (the "Lorentz sphere") and replacing the summation over points beyond R by an integral. The accuracy of this Chapter 2. Methods of Muon Site Determination 15 method obviously increases with the radius R; an estimate of the error may be obtained by recording the variation of the calculated result with the value used for R. M . Bonn [34] has proved that a simple procedure for evaluating a sum over a simple lattice is to obtain by direct summation the contributions by the lattice points X(l) within a certain radius R and to replace the summation over the points beyond this radius by an integral. 2.3 Conclusion From the potential maps I was able to narrow down the areas of candidate fi sites in + the sample and then apply the dipolarfieldcalculations on these areas. The candidate muon sites are then those points which fall into the potential minima and have the correct magnetic field. Chapter 3 Explicit Calculations of Muon Sites In Y B a C u 0 2 3.1 3 6 Potential Maps of A F M Fig. 3.22 depicts a unit cell of the YBa Cu 06 crystal. The positions and properties 2 3 of the ions of YBa Cu 06 are shown in Table 3.3. The ion positions are taken from 2 3 T. Siegrist. [35] Hard core radii are from Handbook of Chemistry and Physics. [36] The valences of barium and yttrium can only be 2+ and 3-f, but Cu can be 1+ or 2-h Copper ions on the Cu-0 chain are thought to have no magnetic moments; they are therefore 1+ ions, and those on the Cu-0 plane are 2+ ions. [37] [38] [39] Cu2 has a 2 slightly smaller radius than Cul. ion Cul Cu2 Y Ba 01 02 03 X-axis (a) 0.5 0.5 0.0 0.0 0.5 0.0 0.5 Y-axis (b) 0.5 0.5 0.0 0.0 0.5 0.5 0.0 Z-axis (c) 0.5 0.143 0.0 0.314 0.347 0.12 0.12 valence haxd core radius 0.74A 0.7A 0.893A 1.34A 1.2A 1.2A 1.2A 1+ 2+ 3+ 2+ 222- Table 3.3: Properties and coordinates of ions in Y B a C u 0 . Here a = b = 3.856A amd c = 11.666A. Origin is at Y(0,0,0). 2 3 6 The calculations are performed over one unit cell at the center of a finite sphere of more than 400 ions generated by a software package, ORTEP. Twenty planes intersecting the z axis at even intervals are chosen in the unit cell, and the potential is calculated on each plane. The deepest minima of the potential maps are the expected 16 Chapter 3. Explicit Calculations of Muon Sites In YBa2Cu C*6 3 17 muon sites. The exact positions can be decided by finely adjusting z values. The following are the potential maps at constant distance along the z direction. The two dimensional potential maps on each plane are represented graphically by single characters which stand for the value of the potential at that point, using the table below to decode the potential map. Generally * represents the highest potential, capital letters are the next highest and numbers are the lowest. The length scales are in units of the lattice constants a = b = 3.856A and c = 11.666A. Coppers are at the four corners of the unit cell, oxygens are at the centers of edges and Yitrium is at (0,0,0). The graphs that follow indicate that there is a very large increase in potential energy close to the ions and that the potential has (naturally) the symmetry of the crystal. Surprisingly, a large region close to the Cu-0 chain in YBa Cu 06, which might make 2 3 an attractive home for a muon, doesn't have a potential minimum. Instead, the minima occur in two planes at 0.03c and 0.19c. Chapter 3. Explicit Calculations of Muon Sites In YBa2CusOs Value of Potential (F) F>36 35 < F < 36 Character Printed 10<F<11 9 < F < 10 A 9 1< F < 2 0< F<1 -1<F<0 1 0 a -26 < F < -25 F< -26 z j * Z 18 Chapter 3. Explicit Calculations of Muon Sites In YBa Cu 06 2 Z = 0.20000E-01 ?5 YrrQ^S Y=0 3 Y=0.5 X=-.5 .6 66 55555 "15555555 j555566666666666666666666( 66666555 >555666666666677777777766( "6666655 .J5555555 555555555 555566666666777889AAA98877] 666666" 5.555555"" "5566666666778ACGJMNMJGCA81 6666666. .. 5555555 556666666789DKV*******VKD9876666< 66666678AGU***********UGA876f" 66666?8AH***************HA8/t 66666679FZ***************ZF9766 6666678CQ*****************QC8?c 6666679H*******************H97( A "66677BN*******************NB7J 66678CU*******************V 6666/8DZ*******************l__ 666678DZ*******************ZD876 666678CU*******************UC87c 666677BN*******************NB?7t 6666679H*******************H976( 666667oCQ*****************qC87f 6666679FZ***************ZF97( 6666678AH***************HA87( 66666678AGU***********UGA876( 6666666789pKy*******VKD98766( 66666666778ACGJMNMJGCA877666t, __,566666666777889AAA98877766666l 5555566666666667777777776666666666555555555£ _5555556666666666666666666666666665555555555£ 51 555555556666666666666666666666666§5555§5555555 55J "5555555556666666666666666666665? X=0.0 55555555555555555555555555555555555555555555555555555555555 X=0.5 Figure 3.2: Two dimensional potential map at Z=0.02c Chapter 3. Explicit Calculations of Muon Sites In YBaiCuzOe 20 30000E-01 t X«-.5 555566666667 .555566666677 6666Z79CJY*********YJC977 5555: ::::~~" 881*************188? ~ 55566666. 9FX**+**********TF9 BK***************KB "6 CP***************PC 6 66678DS***************SD8 "6678DS***************§D8 -6678CP*************** 66678BK*************** 666??9FT*************T 6666678BI************* [BJ A X=0.0 666??§BHP*******Pfi ..6666778ACGKN0NKGCA877 666666677789ABBBA987776.. "666666667777888777766666 666666666667777766666666 566666666666666666666666 556666666666666666666666 5555566666666666666666665 666666 X=0.5 Figure 3.3: Two dimensional potential map at Z=0.03c Chapter 3. Explicit Calculations of Muon Sites In YBaaCuzO^ 5 5 4 4 4 4 5 5 5 5 5 5 5 5 6 6 6 6 6 7 7 9 E Q * * * * * * * * * * * q E 9 7 7 6 6 6 6 6 5 5 5 5 5 5 5 5 4 4 4 4 5 5 5 X=0.0 5544445555555566666779DN***********ND97766666555555554444555 Figure 3.4: Two dimensional potential map at Z=0.05c 22 Chapter 3. Explicit Calculations of Muon Sites In YBavCusOe X=-.5 S*** **** ***** ***** ***** x=o.o X=0.5 Figure 3.5: Two dimensional potential map at Z=0.07c Chapter 3. Exphcit Calculations of Muon Sites In YBa^CusOs • 0.10000 Y-0.5 Y-0 . -0.5 ******************* ******* X=-.5 ******YF97 ******************* ******]tr ******* ****** { *****************U ******* *****************! *****Q|r 9EQ****** ****QFA7665555555556AQ***************QA6555555555667AFQ**** ¥ KHE*: X=0.0 AFQ***** ****Q *****QE9766555555557I*****************I7555555556679EQ****** ******KB76655555556AU*****************UA65555555667BK******* ******pn86655555556E*******************E65 ******YF9765555555 ^^^citc^^^G^Y655555557I****^***** **^* **** * I?o555555679G******** X=0.5 : ! , C Figure 3.6: Two dimensional potential map at Z=0.10c Chapter 3. Explicit Calculations of Muon Sites In YBa2Cu30e Z - 0.15000 Y=-0.5 Y=0.5 *********** X -.5 Y=0 **********H8 **********£7 *********VB7 *********Kgg •*******YC7 ******** *******H *****XF9 ***S^^7 v J** ************* ** ?***************** w***************y 51* ************* * j i ~ n************* iN***********N s *********** v********** 69K********** 7CV********* ********* H******** 9FX****** §§66515. S*** ***** 58G******* 7EY******* 90******** ********* H********* 8K********* 8K********* X=0.0 7H********* 456D********* 44590******** 457EY******* 458G******* 4568EU***** IS*** 8ACEI ••SIB7,,. ****UFg6 *******G85 *****YE *******0 ******** ******** ******** ******** ******** ********n *******09 ******YE^ ******Gr ****U **SIB ECA86554 ?65555"" " 155 - f u t 55 . JDGJK ***SIB8766 *****XF976 *******g ******** •*******VC7 *********Kg *********VB76 •*****»***E76 **********H86 **********X86 5556A0*************0A6 "5581***************I8 BV***************VB6 7F*****************?? j*****************j K*****************R BIS**** X****** H******** ********* 7CV***+***** gK********** V********** *********** *********** ! *********** X=0.5 a S Figure 3.7: Two dimensional potential map at Z=0.15c Chapter 3. Exphcit Calculations of Muon Sites In YBaaCu^O^ Z « 0.17000 : Y=-0.5 Y=0 Y=0.5 #********G8665555554557CN* X=-.5 *********E8665555^ ********QB76b55555544568DM*******MD^ ********H9665555555445569ELTZ*ZTLE9655445555555669H********* *******p7665555555444^^ ******nfi§6655555555544455 *****MC8766555555555444445556665554^ m jfrp***** J8 JGDA876 ***PGB8 *****MC87 ******Qf)8 *******NB76 ********H96 ********QB765 *********E866 *********G86 *********H$6 j7AJY*******YJA] J7CN*********NC1 >8DP*********PD{ x=o.o 678ADGJK 8BGP**** 8CM****** 8DQ******* 7{3N******** 9H********* BQ********* ********** ********** * * * * * * * * * * X=0.5 Figure 3.8: Two dimensional potential map at Z=0.17c Chapter 3. Explicit Calculations of Muon Sites In YBa^Cu^Oe Figure 3.9: Two dimensional potential map at Z=0.19c 26 Chapter 3. Explicit Calculations of Muon Sites In YBaaCu^Oe Figure 3.10: Two dimensional potential map at Z=0.20c 27 Chapter 3. ExpHcit Calculations of Muon Sites In YBaiCu^O^ 555555555555566666666789CIQZ***ZqiC9876666666655555555555555 Figure 3.11: Two dimensional potential map at Z=0.26c 28 Chapter 3. Exphcit Calculations of Muon Sites In YBaaCuzOs L 30000 AL****** X« - . 5 9HZ***** 7CN***** **** ^^f ^||f If ppp 6 56l66677777776f66666 578BGP*******PGB r9DN***********N )ET************* s***************sp E?******************* * * * * * * * * * * * * * * * * * L j******************* j******************* >666666667BR*******************RB766666666 " 366666667BR*******************RB26666666' )6666o667AQ*******************QA76666666 >666666679J*******************J976oo666o >666666678E*******************E876666666 7AL*****************LA76 566666 78DS***************SD87c 6l9ET*************TE97f 6779DN***********ND9?7c 566666678BGP*******PGB8766e 555555! mm .CRJD96! ***0E96! ****NC7! ****ZH9( *****}• *****] 566666666667777777666( 566666666666666666666( 566666666666666f 555555555555-, 1555555555555J )555555555555? 1555555555555! Ml X=0.0 59EO**** rcR***** )HZ***** IT.****** jff****** x* 0 . 5 Figure 3.12: Two dimensional potential map at Z=0.30c Chapter 3. Exphcit Calculations of Muon Sites In YB&iCuzOe Z - 0.35000 Y=-0.5 **********$ *********TQ *********K? Y=0.5 * * * * * * * * * * * X=-.5 T********** K********** g£********** *********C6 K********* ********K *******QA ^q******** ******PB6544555555555555555555555555555555555554456BP******* X=0.0 55 ...3765 **SIB7 ***•*! ******p ******* ********K *********06 *********£? *********T9 **********g **********£ BIS*** ****** p******* ******** K********* ********** K********** 7********** *********** *********** Figure 3.13: Two dimensional potential map at Z=0.35c X=0.5 0=X * * * * * * * * * ********* ********n ********T *******rt(] ******** *****rt ***Aia89 .frfr- "9 ******** ******** n******* 7******* rjrt****** o****** "DM**** " "A** sssgfl "If 0=X ***A1Q89 *****rt06 *******£ *******rt< ******** ********flfl *********Q S'0=A vS9§aiA** rt**** ****** cm****** 1******* flfl******* Q******** 9-0-=• 0008S*0 Chapter 3. Explicit Calculations of Muon Sites In YBa^Cu^Oe 0.40000 0 " **ZMDB655455555555555555555555555555555555555555545568D8Z*** X=-.5 **RIB?E vqJD96S IFC975E X=0.0 >55§55 55! 55J **RIB75555^ Figure 3.15: Two dimensional potential map at Z=0.40c Chapter 3. Explicit Calculations of Muon Sites In YBavCu 0e 3 Figure 3.16: Two dimensional potential map at Z=0.42c 33 Chapter 3. Explicit Calculations of Muon Sites In YBaaCu^O^ = 0.45000 : - -0.5 Y=Q ****MC8665555555555555555555 ***UIB866555555555555555555555555555555555555555 34 Y=0.5 X=-.5 X=0.0 56 555555555555555555555! 15f 55 ]66§§555§55§5§^ VQJEA7665555555555555^ Figure 3.17: Two dimensional potential map at Z=0.45c Chapter 3. Explicit Calculations of Muon Sites In YBa2Cu 0 3 6 35 -0.6 H 1 1 1 1 1 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 Figure 3.18: The Contour Potential Map of Y B a C u 0 at Z=0.03c. Coppers are at the four corners. Larger number lines represent lower potentials. 2 3 6 Chapter 3. Explicit Calculations of Muon Sites In 36 YB&iCu 0s 3 The contour maps in Fig. 3.18 and Fig. 3.19 show the inverted muon potential energy — U(x, y) in the a — b plane (Cu ions at the four corners). Larger-numbered contour lines represent lower potential energy and indicate likely fi sites. The scales + are in units of the lattice constants (a, 6, c). Yttrium is at the center. Fig. 3.18 describes the plane 0.03A away from Y plane and between the CuC>2 planes - and Fig.3.19 the plane 0.8lA below the Cu02 planes —i.e., toward the CuO chains. These contour maps confirm the strong fi+—O 2- attraction that binds all muons closely to oxygen ions. For convenience, I define the z = 0.12c copper oxide plane as a reference plane; here yttrium is at (0,0,0). The first muon site (0.45a,0.06,0.03c) is 1.07A "down" away from 02 in the plane. The other site is (0.3a,0.06, 0.19c) which is 1.125A "above" away from 02 in the same plane. There are another six sites at symmetric positions on these planes. 3.2 Dipolar Magnetic Field Calculations For antiferromagnetic YBa Cu 06 the large Cu ion paramagnetic moments dominate 2 3 the dipolar magnetic field, which is thus ~ 2000 times larger than that from the nuclear moments. The magnetic structure consists of strong nearest-neighbor correlated spins in adjacent copper planes which are aligned antiferromagnetically. The moments are constrained to lie in the tetragonal a — b plane; I assume that they point along either the a axis or the 6 axis. From neutron scattering measurements, the ordered moment on the Cu ions in the Cu02 plane is estimated to be 0.66 /XB[40]. J.W. Lynn[41] reports that at lower temperatures the Cu moments in the oxygen-deficient CuO chain layers also order antiferromagnetically, with a moment that can be quite substantial (0.5 / / B ) ; however, in my calculations I only consider the contributions from Cu moments in the Cu02 planes. 37 Chapter 3. Explicit Calculations of Muon Sites In YBa^CuzOe Figure 3.19: The contour potential map of Y B a C u 0 at Z=0.19c 2 3 6 Chapter 3. Explicit Calculations of Muon Sites In YBa^CusOe 38 The calculated results are in Table 3.4 and Table 3.5. Here the magnetic field B = xjB 2 x + B 2 v + B \ z X-axis 0.05 0.10 0.45 0.05 0.50 Y-axis 0.10 0.00 0.00 0.00 0.05 Z-axis 0.03 0.03 0.03 0.03 0.03 Field(Gauss) 308.88 318.57 307.86 287.17 302.79 Table 3.4: The magnetic field at z=0.03c. X-axis 0.25 0.30 Y-axis 0.00 0.00 Z-axis 0.19 0.19 Field(Gauss) 268.45 308.83 Table 3.5: The magnetic field at z = 0.19c. Only the points at which the local magnetic field has a magnitude of around 300G are listed here. Contour maps of the magnetic field strength on these two planes are shown in Fig. 3.20 and Fig. 3.21. Chapter 3. Explicit Calculations of Muon Sites In YBa^Cu^Oe 3.3 39 Comparison W i t h Experimental Data Experiments [14] on A F M YBa Cu 06 below 100K obtain a dominant ZF-fiS1Z signal 2 3 at about 4 MHz, which indicates an internal field of about 300 Gauss. Meanwhile, the LCR experiment shows that the muon sites in YBa Cu307 are 1 A away from the 0 2 in the C u 0 plane. If we assume that the muon occupies the same site regardless of 2 oxygen deficiency, the overlap between the regions with potential minima and those with the magnetic field within uncertainty will be the candidate muon sites. Table 3.6 lists the muon sites in YBa Cu3C*6. 2 X-axis a 0.45 0.30 Y-axis 6 0.00 0.00 Z-axis c 0.03 0.19 Distance A 1.070 1.125 Field(Gauss) 307.86 308.83 Table 3.6: Muon candidate site in Y B a C u 0 . Here a = b — 3.856A and c = 11.666 2 A. 3 6 Chapter 3. Explicit Calculations of Muon Sites In YBa Cu 0 2 eft? 3 # 40 6 $ 1- 2. 00E+02 too. 0 2= 2. 50E+02 99. 6 3- 3. 00E+02 96. 7 4- 3. S0E*02 92. 1 5« 4. 00E*02 87. 3 6- 4. 50E+02 8 1 . 6 7- 5. 00E+02 74. 9 8* 5. SOE + 02 70. 3 9- 6. 00E*02 66. 3 H 0 - 6. 50E+02 62. 8 11* 7. 00E + 02 59. 4 12» 7. 50E+02 56. 2 13- 8. 00E+02 52. 7 14- 8. 50E*02 49. 5 15- 9. DOE+02 46. 7 16* 9. 50E+02 44. 3 17. 1. 00E+03 42. 0 - 1 8 - 1. 05E+03 39. 9 37. 9 19- 1. 10E+03 20- 1. I5E+03 36. 0 4 00. 0 99. 9 99. 1 97. 7 96. 0 93. 6 90. 6 88. 3 86. 1 84. 0 81. e 79. 6 77. 0 74. 5 72. 1 70. 0 67. 9 65. 8 63. 7 61. 7 .5 Figure 3.20: The contour plot of magnetic field at Z = 0.03c. It shows only a quarter of the plane. Coordinates are in' the unit cell length. Chapter 3. Explicit Calculations of Muon Sites In YBa Cu 0 2 3 6 Figure 3.21: Contour plot of magnetic field at 0.19c Chapter 3. Explicit Calculations of Muon Sites In 42 YBa Cu 0 2 3 6 Cu1 Figure 3.22: The crystal structure of antiferromagnetic Y B a C u 0 2 3 6 Chapter 4 Calculations on Y B a C u 0 2 4.1 3 7 M u o n Site Decision Experiment ZF— and WTF—fxSTZ measurements were made on normal and 0 enriched samples 1 7 of YBa Cu 07 at a temperature of 103K. In the latter samples, 38 % of the naturally 2 3 abundant and spinless 0 was replaced by 0 . At this temperature the correlation 1 6 1 7 time for muon hopping is much longer than the muon lifetime so that the nuclear dipolar fields seen by the muon are essentially static. The ^SR time spectra for the 1 7 0 enriched sample are compared with corresponding spectra in the pure 0 control 1 6 sample. The depolarization rate in the 1 7 0 substituted sample is enhanced in both ZF and wTF(24G), which suggests that the muon is close to an oxygen ion. The rms internal field was estimated by fitting the wTF relaxation function to a sum of two Gaussians. The fitted results yield an rms internal field (i? ) ^ along the applied field 2 1 2 direction equal to 1.31G in the unsubstituted YBa Cu C"7 and 2.65G in the 0 enriched 1 7 2 3 sample. If the muon is much closer to a single oxygen (i.e. if it forms an hydroxyl-like bond) such that the dipolar fields from other oxygens can be neglected, one obtains a powder averaged (B ) ^ 2 Z X 2 due to a single 0 equal to L\ =3.64G. The 0 nuclei then 1 7 1 7 0 exert an effective static dipolarfieldon the muon given by Bdi = (^JN/ 3)S [3(f R p g -q)f—q] where S is the component of nuclear spin along the electric field gradient direction q q and r is the unit vector between the muon and 0 nucleus separated by distance r. 1 7 43 Chapter 4. Calculations on YBa Cu 07 2 44 3 Averaging over all angles between Bdip and z yields A = \{n~t /r*)[S(S 0 N + l){3(f • qf + l}] ' , 1 2 from which r = 1.01 A can be estimated provided that the muon is bonded to a single oxygen and the bond axis is perpendicular to q. LCR was introduced in order to decide to which oxygen the muon is bonded. Principally, LCR is based upon the idea that when two systems are prepared so that an energy splitting in one system matches a splitting found in the second system, an exchange of energy or polarization may occur between levels of the two systems. In this study, the system consists of a spin 1/2 muon in an external field interacting with a spin 5/2 nucleus 1 7 0 in an electric field gradient plus external field. Whenever the Zeeman energy of the muon matches the quadrupolar level splitting of 0 , there may 1 7 be a resonance. The spin Hamiltonian for a p in the crystal is: + H = Hdip + HQ + H H = dip H q = H ^^[s .I-3(s .r)(I.r)] 47(2J-l) (4.9) M M - " [ 3 ( 7 2 / 2 ) + ^ 2 " ^ 2 ) 3 ( 4 = (u„ + uo) • H , z where Hn (4.8) z ' 1 0 ) (4.11) e x is the dipolar interaction between muon and 0 , HQ is the quadrupolar 1 7 p interaction of 0 , Hz is the combined Zeeman energy of the muon and the 1 7 external field, r is the separation between nucleus and p , s + M magnetic moments of /x and 0 and 6 = (V + 1 7 xx — Vyy)IV . xz 1 7 0 in the and I are the nuclear Here VJj is the electric field gradient (EFG) tensor. Due to the two splittings between |l/2) and |3/2) and between |3/2) and |5/2), two resonances are expected at approximately B = ^^QQ/iy^) r 2B , assuming 6 = 0. r ^d Chapter 4. Calculations on YBa Cu 07 2 45 3 From the data the best estimate of the average quadrupolar parameter is e qQ/h =6.6MHz, 2 which is very close to that measured by NMR for the 0(2,3) [42] planar sites in the absence of the muon[42]. These results suggest that the muon is bonded to an oxygen in the C u 0 plane at a distance of l A . 2 4.2 Potential M a p of Superconductor Y B a C u 0 7 2 ion Y Ba Cul Cu2 01 02 03 04 X(a) 0.0 0.0 0.5 0.5 0.5 0.0 0.5 0.5 Y(6) 0.0 0.0 0.5 0.5 0.5 0.5 0.0 0.0 Z(c) 0.0 0.347 0.5 0.143 0.312 0.12 0.12 0.5 3 Valence 3+ 2+ 2+ 2.3+ 2222- Hard core radius 0.893A 1.34A 0.74A 0.7A 1.2 A 1.2A 1.2A 1.2A Table 4.7: The positions and properties of YBa Cu 07. Here a = 3.856A, b = 3.87A and c = 11.666A. 2 3 Eq. 2.1 is used to generate the potential maps of YBa Cu 07, which have a structure 2 3 similar to those of YBa Cu 06 except that there is no potential minimum between two 2 3 Cu—0 planes. Furthermore a minimum-potential area appears around the bridging 2 oxygen 04 due to the two extra negative ions on the Cu chain, which decrease the potential around that region. Further studies of region will be conerned with the dipolar magnetic field. The plane at z = 0.19c has a potential map similar to those of YBa Cu 06. The 2 3 contour plot of this plane shows that there is also a very strong interaction between the muon and the ions. Muons are evidently bonded to the oxygen ions in the copperoxygen plane. 46 Chapter 4. Calculations on YBa^CuzOj 4.3 Dipolar Magnetic Field For YBa2Cu 07 the dipolar fields are due to weak nuclear moments, primarily from 3 ^ C u and C u nuclei but with minor contributions also from Ba and Y nuclei. For 6 5 17 0-doped samples the substantial contributions from 0 nuclear moments must also 1 7 be included. The spin-3/2 copper nuclear moments will precess about the electric field gradient, which is along the c-axis for Cu ions in the C u 0 plane and along the chain 2 for Cu ions in the CuO chains. This quadrupolar precession is assumed to average out all transverse field components; those remaining are all treated as static (i.e., any muon-nuclear "flip-flop" dynamics are neglected). Within these conditions, the nuclear moments are assumed to be oriented along the electric field gradient and the resultant local fields are characterized by a gaussian width A as usual. 47 Chapter 4. Calculations on YBa CU3O7 2 Ion ^Cu Cu 0 Ba Ba 6 5 1 7 135 137 89y Nuclear Magnetic Moment pi, Abundance % 2.226 69.01 2.385 30.91 38 -1.8937 0.8365 6.59 11.32 0.9357 -0.13682 100 Table 4.8: The nuclear magnetic moments of ions in YBa2Cu30 . x The calculations indicate that the muon relaxation rate of 0 doped sample should 1 7 be larger than that in the unsubstituted sample, which is consistent with the experiments. The minimum potential region at z = 0.19c yields a theoretical value of A consistent with those measured experimentally. Other regions with potential minima don't have dipolar magnetic fields of the right magnitude. Note that no consideration has been given to any lattice relaxation or hopping of muons, which would significantly affect the relaxation. 4.4 Conclusion The combined results of potential maps and dipolar magnetic field calculations indicate that the most likely muon site is on the plane at z = 0.19c, which is between the Cu—0 2 plane and the Cu—O chain and 0.84A "down" from the Cu—0 plane, which is indicated 2 on Fig.4.4. All the other regions either have minimum potentials but wrong dipolar magnetic fields or have the correct magnetic field but don't fall into the minimum potential areas. In conclusion, the most probable muon site is at (0.3a, 0.06,0.19c) and corresponding symmetric sites. Chapter 4. Calculations on 48 YBa^Cu^ 2. 00E-01 4. 00E-0I 6. 00E-0I 8.00E-01 1.00E»00 1. 20E*00 l.40E«00 1.60E»00 0.6 0.4 H i.80E»oo 2. OOE*00 2. 20E»00 2. 40E«00 2. 60E+D0 2. 8DE«00 3. OOE'OO 3. 2DE*00 3. 40E'00 3.60E*00 3. BOE'OD 4.O0E*00 0.2 0.0 -0.2 94.5 100. 0 34. 1 100.0 — 33. 9 99. 9 93. 6 99.9 93.4 93. 8 93. 2 99. 7 92.9 99. 6 92. 7 99.5 92. 4 99. 4 92.0 99. 2 91.6 98. 9 98. 4 69. 9 97.7 89.2 96. 3 85. S 93.9 81. 2 89.9 68.5 77. t 41.1 48. I 24.3 29.4 13.4 16.6 H -0.4 -0.6 -| 1 1 1 -0.6 -0.4 -0.2 0.0 1 0.2 r 0.4 0.6 Figure 4.23: The contour plot of the potential map at Z=0.19c 49 Chapter 4. Calculations on Y B a C u 0 7 2 3 Cu2 ^ s ite Cu1 Figure 4.24: The crystal structure of superconducting YBa Cu 07. 2 3 Bibliography [1] A. Schenck and K . M . Crowe, Physics Review, 26:57, 1971. [2] H. Graf et al., Phys. Rev. Lett, 37:1644, 1976. [3] E . Holzshuh, A.B. Dension, W. Kundig, P.F. Meirer and B.D. Patterson Phys. Rev., B27.-5294, 1983. [4] C. Boekma and P.L. Lichti, Phys. Rev. B., 31 1233, 1985. [5] S.R. Kreitzman et al., Phys. Rev. Lett. 63:1865,1986. [6] B. Foy, N. Heiman and W.J. Koss\er,Phys. Rev. Lett, 30:1064, 1973. [7] Y.J. Uemura, W.J. Kissler, X.H. Yu, J.R. Kempton, H.E. Schnoe, D. Opie, C. E . Stronach, D.C. Johnston, M.S. Alvarez and D.P. Goshun, Phys. Rev. Lett., 59, 1045, 1987. [8] R. Kiefl et al., Phys. Rev. Lett, 60, 224, 1988. [9] G. Golnik et al., Phys. Lett. A, 125:71, 1987. [10] C. Kittel, Introduction to Solid State Physics, 1976. [11] J. Pearson, Handbook of Lattice Spacings and Structure of Metals, 1963. [12] Y. Le Page et al., Phys. Rev. Lett, 36 3617, 1987. [13] J . M . Tranquada, D.E. Cox, W. Kunnmann, G. Shirane, M . Suenage, P.Zolliker, D. Vaknin, S.K. Sinha, M.S, Alvarez, A. J . Jacobson and D.C. Johnston,Phys. Rev. Lett, 60 156, 1988. [14] N. Nishida et al., J. Appl. Phys. Pt.2, 26:L1856, 1987. [15] A. Abragam, CR. Acad. Sci. Ser. 2, 299 95, 1984. [16] S.R. Kreitzman, Hyperfine Interaction, 31:13, 1986. [17] R.F. Kiefl, Hyperfine Interaction, 49 233, 1988. [18] S. Kreitzman, J.H. Brewer, P.R. Harshman, R. Keitel, D.LI. Williams, K . M . Crowe and E.J. Ansaldo, Phys. Rev. Lett, 56 181, 1988. 50 Bibliography 51 [19] R.F. Kiefl, Eyperfine Interaction, 32 707, 1986. [20] J.B. Bates, J.C. Wang and K.A. Perkins, Physical Review B, 19:4130, 1979. [21] Proceedings of the 5th International Conference on Muon Spin Rotation, 1990. [22] P. Lowdin, A Theoretical Investigation into Some Properties of Ionic Crystals, 1954. [23] J.H. Brewer et al., Phys. Rev. Lett, 60:1073, 1988. [24] R.F. Kiefl et al., Physica C, 162-167:161, 1989. [25] B.G. Dick and A.W. Overhauser, Phys. Rev., 112 90, 1958. [26] R.G. Gordon and Y.S. Kim, / . Chem. Phys., 56:3122, 1972. [27] R.I.M. Rae, Chem. Phys. Ze«,18:574-7, 1972-73. [28] W.C. Mackrodt and R. F. Stewart, / . Phys., C12431, 1979. [29] C.R.A. Catlow, K . M . Diller and M.J. Norgett, J. Phys. C, 10 1395, 1977. [30] W.C. Mackrodt Superconductivity, 343, 1988. [31] P.F. Meier, Hyperfine Interaction, 8 591, 1981. [32] W.K. Dawson et al., J. Appl. Phys. 64:5809, 1988. [33] G.A. Sawatzky and Van der Woude, J. Phys., C647, 1974. [34] M . Bonn and K. Huang, Dynamics Theory of Crystal Lattice, 1956. [35] T. Siegrist, S. Sunshine, D.W. Murphy, P.J. Cava and S.M. Zahurak, Phys. Rev. B 35:7137, 1987. [36] The Chemical Rubber Co., Handbook of Chemistry and Physics, 1974. [37] Y. LePage, W.R. Mckinnon, J.M. Tarascon, L.H. Greene, G.W. Hull and D.M. Hwang, Phys. Review. B, 35:7245, 1987. [38] J.M. Tranquada et al., J. Appl. Phys., 64:6074, 1988. [39] F . Beech, S. Miraglia, A. Santoro and R.S. Ruth, Phys. Rev. B, 35 8778, 1987. [40] J.M. Tranquada et al., J. Appl. Phys., 64:6074, 1988. [41] J. W. Lynn et al., J. Appl. Phys., 64:6065,1988. Bibliography 52 [42] E . Oldfield, C. Loretsopouls, S. Yang, L. Reven and H.C. Lee, Phys. Rev. B, 40, 6832, 1989. [43] W.W. Waren et al., J. Appl. Phys., 64 6081, 1988. [44] C H . Pennington D.J. durand, C P . Slichter, J.P. Rice, E.D. Bukowski and D.M. Ginsberg, Phys. Rev.B, 39 2902. [45] F . van der Woude and G.A. Sawatzky, Phys. Rev. B, 4 3159, 1971. [46] J . Brewer, MSR group, TRIUMF, Vancouver, Canada. [47] N. Nishida et al., Dept. of Physics, Tokyo Institute of Technology, Ohokayma, Meguro-ku, Tokyo, 152. [48] G. Solt et al., Paul Scherrer Institute, CH-5232 Villigen PSI. Appendix A Programs To determine the candidate muon sites a calculation consisting of three steps was performed. The first step is to generate a cluster of lattice points by using X T L P L O T ; the second step is to produce potential maps on the planes at equal separation along the z—axis in a crystal unit cell to decide possible muon stopping regions; and the last step is to calculate the magnetic field on the different planes to verify that those possible muon stopping regions have the correct magnetic field. This initial introduction to the programs will emphasize the relationships between the programs. Details about the procedures followed in each program can be found in the corresponding manuals. X T L P L O T provides information to another program ORTEP, which actually produces the model cluster. One of 230 different space groups has to be specified by users in oder to generate a crystal structure. 10 space groups are available in the files USR2:[ALANA]*.SPG on the TRIUMF VAX cluster. A special space group can be generated by editing a *.SPG file in which the specific transformations are indicated. Different radii have been employed to distinguish different ions on the plot because of the difficulty of labeling the ions. In this study two kinds of units are used: unit cell length and actual length. Positions are in the unit cell length, radius are in actual length. P O T E N T M A P program generates electric potential maps of a crystal structure over the x — y plane at the fixed z values specified by the users, minima of which are likely muon sites. The potential maps are characterized by letters; each letter 53 Appendix A. Programs 54 represents potential levels, which can be decoded by the table in the documentation. POTENTMAP's preselection of candidate muon sites is very useful, because each calculation of the dipolar magnetic field consumes at least half an hour of CPU time on the VAX. From the potential maps it will be easy to decide where the magnetic field needs to be calculated. The Morse potential is used to calculate the interaction between muon and ions. No consideration was given to any lattice relaxation or charge screening effects. The calculation involves more than 400 ions in the cluster. Some modifications of POTENMAP can be made to generate output in a format suitable for input to PLOTDATA which can produce beautiful contour graphs or density maps. Both methods were employed in this thesis. The lettered maps show pseudo-threedimensional structure and the contour map gives more detail of the two dimensional potential at every z value. CLUSTSUM calculates the dipolar magnetic field caused by other magnetic dipole moments at a site. Because of the limit of the size of the cluster generated by X T L P L O T , the contribution from those outside the cluster will be replaced by an integral. The sum will vary with the radius of the cluster. POTENMAP and CLUSTSUM read in cluster files generated by X T L P L O T . Because the cluster usually involves many more ions than on one plot, no plot should be drawn after one sees the prompt XTLPLOT). Also two kinds of clusters are optional: "unit-cell" shaped, to draw the entire structure of a unit cell; and "sphere", to look at the environment of a particular position. In the case of this thesis, the sphere is usually used. For POTENTMAP, a sphere centred on (0,0,0) is produced. For CLUSTSUM, a sphere centred on the candidate muon position is required. The cluster file created will be in the proper format to be read by the two programs POTENMAP and CLUSTSUM. The only thing that must be done before they are run is that the appopriate ion property be included in lieu of the "**•***" y O U w j U faA Appendix A. Programs 55 in the newly created file. For the program POTENMAP, this property would be the ion's hard core radius in your chosen length unit; for CLUSTSUM, it would be the ion's dipole moment, in Bohr magnetons or nuclear magnetons, depending upon which dominates the field. In the case of Y B a C u 0 , Bohr magnetons are used; nuclear 2 3 6 magnetons are used for YBa Cu 07. 2 3 POSITION, PLOTPARAMETERS and XTLMAIN, subroutines of X T L P L O T , play the roles of reading in ion positions, plotting parameters and lattice constants respectively. It is possible to revise the programs for a batch job. One example is usr2: [qiang.thesis.program]li.for, which calculates the magnetic field at 100 evenly distributed points on the plane specified by the user in YBa Cu 06 with the 2 cluster generated automatically. 3
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Theoretical calculation of muon site in YBa₂Cu₃O
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Theoretical calculation of muon site in YBa₂Cu₃O Li, Qiang 1990
pdf
Page Metadata
Item Metadata
Title | Theoretical calculation of muon site in YBa₂Cu₃O |
Creator |
Li, Qiang |
Publisher | University of British Columbia |
Date Issued | 1990 |
Description | The muon is a useful probe of magnetic fields in superconductors, but knowing the field seen by the muon is often of limited value until we know where the muon is in the crystal lattice. In this thesis I employ two independent theoretical methods to search for candidate muon sites: the potential energy field method, which seeks the minimum of the electrostatic potential of the μ⁺, and the magnetic dipolar field method, which compares the calculated magnetic field (due to the host electronic, atomic or nuclear dipolar fields) with the observed local fields at the muon. |
Subject |
Muon spin rotation Magnetic dipoles |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2010-10-21 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0302373 |
URI | http://hdl.handle.net/2429/29415 |
Degree |
Master of Science - MSc |
Program |
Physics |
Affiliation |
Science, Faculty of Physics and Astronomy, Department of |
Degree Grantor | University of British Columbia |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-UBC_1990_A6_7 L53.pdf [ 4.07MB ]
- Metadata
- JSON: 831-1.0302373.json
- JSON-LD: 831-1.0302373-ld.json
- RDF/XML (Pretty): 831-1.0302373-rdf.xml
- RDF/JSON: 831-1.0302373-rdf.json
- Turtle: 831-1.0302373-turtle.txt
- N-Triples: 831-1.0302373-rdf-ntriples.txt
- Original Record: 831-1.0302373-source.json
- Full Text
- 831-1.0302373-fulltext.txt
- Citation
- 831-1.0302373.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0302373/manifest