Zero Field E S R of Magnetic Impurities in Superconductors: a Novel Technique for Measuring Magnetic Penetration Depth and Application to High T Y B a C u 3 0 _ 5 c by Tamar Pereg-Barnea B . S c , The Hebrew University, 1999 A THESIS S U B M I T T E D IN P A R T I A L F U L F I L M E N T O F THE REQUIREMENTS FOR T H E DEGREE OF M A S T E R OF SCIENCE in T H E F A C U L T Y O F G R A D U A T E STUDIES (Department of Physics and Astronomy) We accept this thesis as conforming to the required standard T H E U N I V E R S I T Y O F BRITISH C O L U M B I A October 12, 2001 © Tamar Pereg-Barnea, 2001 2 7 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. Department of Physics and Astronomy The University Of British Columbia Vancouver, Canada Date dci-aber V2~, tool Abstract ii Abstract The Meissner effect in superconductors is characterized by an exponential decay of externally applied magnetic fields into the sample with a length scale A, that is called the magnetic penetration depth. A is directly related to the superfluid density of the superconductor and is therefore one of the most important probes of the intrinsic properties of superconductors. This work describes the development of a new technique for measuring the penetration depth at microwave frequencies and low temperatures. The surprising combination of two seemingly unrelated experiments has led to a direct method of measuring the magnetic penetration depth in the London limit of high T superconductors. Samples of high T Y B C O with uniform doping of magnetic ions (Gd Yi- Ba2Cu 0 -6 where x w 1.28%) were prepared. The magnetic ions embedded in the superconductor serve as magnetic field probes through a zero field Electron Spin Resonance (ESR) experiment. Measurement of the imaginary part of the magnetic susceptibility, x"(o;), through the sample's microwave surface resistance, reveals relatively sharp E S R transitions with integrated intensity that is proportional to the number of ions that are exposed to the field. W i t h a uniform ion distribution, the number of ions that are exposed to the field is proportional to the penetration depth A. The spectrum of is obtained from surface resistance measurement in a bolometric apparatus. Our new method was successfully applied to overdoped and optimally doped Y B C O samples and absolute values of the a-b plane penetration depth at T = 1.2 are given. c c x x 3 7 Table of Contents iii Table of Contents Abstract ii Table of Contents • iii L i s t of Tables v List of Figures vi Acknowledgments C h a p t e r 1: Introduction and Motivation C h a p t e r 2: C o n d u c t i v i t y a n d Surface Impedance 2.1 The Two Fluid Model 2.2 Conductivity in the Two Fluid Model 2.3 Surface Impedance vii 1 4 4 4 5 C h a p t e r 3: Penetration depth 3.1 Meissner effect and the London theory 3.2 London Penetration Depth 3.3 Magnetic Penetration Depth 9 9 11 12 C h a p t e r 4: E l e c t r o n S p i n Resonance 4.1 Magnetic Resonance 4.2 Magnetic Ions in a Crystal: Fine Structure 4.3 Effective Spin Hamiltonian 4.3.1 Semi-Classical Approach for Deriving the Effective Spin Hamiltonian 4.4 Zero field E S R and the E S R Spectrum 14 14 16 17 18 20 C h a p t e r 5: The Bolometric Experiment 5.1 The Measurement Principle 23 24 C h a p t e r 6: A E x t r a c t i o n Technique 6.1 The Samples 6.1.1 Sample Preparation 6.1.2 Gadolinium Concentration Determination 6.2 R and x" 6.3 Background Subtraction 6.4 Spectrum Fitting and Extraction of A 27 27 28 30 31 32 32 s Table of Contents iv Chapter 7: Results and Analysis 7.1 Twinned Optimally Doped Crystal 7.1.1 Oxygen Ordering 7.2 Detwinned Overdoped Crystal 7.3 Corrections for c Axis Currents 7.4 Error analysis 34 34 37 40 42 45 Chapter 8: 47 Conclusions References 49 Appendix A : Cavity Perturbation Technique 53 Appendix B : Surface Current in a Uniform magnetic Field 55 List of Tables v List of Tables 7.1 7.2 7.3 7.4 7.5 Crystal field parameters for the two E S R bands in an optimally doped Y B C O , given in G H z E S R Transitions in a twinned optimally doped Y B C O sample Crystal field parameters for the E S R spectrum in an overdoped Y B C O crystal, given in G H z E S R Transitions in a detwinned overdoped Y B C O sample The ESR-extracted penetration depth A of an optimally doped twinned sample and an overdoped detwinned sample, at T — IK. Both samples have 1.28% G d concentration 38 40 42 42 44 List of Figures vi List of Figures 2.1 A schematic picture of the complex conductivity of a superconductor as described by the two fluid model in different temperature regimes. 6 4.1 E S R transitions in a spin system 15 4.2 E S R transitions at a constant frequency and varying D C magnetic field for a spin | system 21 4.3 Theoretical tonian f° system with spin | and fine structure Hamil- r a U = B° 0° = B ° ( 3 S 0 2 2 2 2 Z - S ) 22 2 5.1 T h e electromagnetic fields configuration in the bolometric apparatus. 5.2 A typical surface resistance spectrum of a Y B C O sample at different 25 temperatures as measured in the bolometric apparatus 26 6.1 T h e Y B C O unit cell 29 7.1 T h e surface resistance of a twinned optimally doped sample with 1.28% Gadolinium . 35 7.2 7.3 T h e magnetic susceptibility x" of Gdo,oiYoQgBa2Cu^OQ^2 T w o different crystal field environments 7.4 T h e fitted E S R spectrum of an optimally doped Y B C O sample with 1.28% G d 39 7.5 Surface resistance of an over doped Y B C O sample with 1.28% G d . . . 41 7.6 T h e fitted E S R spectrum of an over doped Y B C O sample with 1.28% Gd B .„1 36 37 43 A typical high T sample in a magnetic c field 56 Acknowledgments vii Acknowledgments I would like to thank D r . Walter Hardy for the great idea behind this work and for his guidance and support. It was a pleasure to work with and learn from Walter at all times. Special thanks to Patrick Turner, who is my partner in this project, for his devotion to the work and for teaching me a great deal about bolometry and microwave techniques. T h i s work has also benefitted from the help of Richard Harris who was always willing to discuss physics, consult and even wrote the first spectrum analysis Mathematica program. T h e past year in the U B C superconductivity lab was challenging, interesting and very pleasant, thanks to everyone in the lab. I would like to thank D r . Doug B o n n for co-supervision, D r . Ruixing Liang for the high quality samples and chemical support and the other members of the group: D r . D a v i d Broun, Pinder Dosanjh, Jennifer Babiak, A h m a d Hosseini, Saeid K a m a l and Darren Peets who have contributed to this work directly and indirectly through discussions, editorial assistance or simply by making it fun to be a part of the group. We have also received help from D r . Andras Janossy whose previous study of G d doped Y B C O was the basis of this work and the correspondence with him, in person, was enlightening. Eventually I would like to thank my family for their support and my dear husband, A m i r , for being a source of comfort and inspiration. Chapter 1. Introduction and Motivation 1 1: Introduction and Motivation Since 1986, when Bednorz and Miiller first discovered high temperature superconductivity in Lcii.gsBao.isCuO^, the topic of high T superconductivity has remained one c the most interesting open questions of condensed matter physics. High T superconc ductors are not only superconducting at temperatures much higher than the boiling temperature of nitrogen, which opens up the possibility of many new applications, but they are also very different from the well known "conventional" superconductors. While previously existing superconductors were well understood by the B C S theory which is based on a phonon mediated s-wave pairing mechanism [1], the high T c superconductors show many deviations from this behavior. T h e London penetration depth, XL, provides also colloquial useful information about the nature of the superconductor. It will be shown in chapter 3 that 1/A^ oc n /m* s where n is the density of the supercurrent carriers (the supercarriers) and m* is their s effective mass. T h i s result motivates the study of the penetration depth, since XL and its temperature dependence provides access to two of the superconductor's intrinsic properties, the supercarrier density and the energy spectrum of the quasiparticle excitations. It is well established that except in special cases, superconductors have an energy gap, A , in their density of states spectrum. T h i s can be seen, for example, in S T M experiments. In the language of the B C S theory, superconductivity is the result of pairing between electrons (Cooper pairs). T h e gap energy is a measure of the energy required to break a pair. A t low temperatures, the system is condensed in its paired ground state (also called the condensate) and excitation cannot be created by collisions. T h i s leads to frictionless motion of the condensate. of the gap in k-space, A(fc), However, the shape may vary from one material to another. main differences between the high T c One of the and the conventional superconductors is the symmetry of the superconduction gap. A conventional superconductor has an s-wave gap symmetry, meaning that A is the same in every direction: A ( & ) — A n . More generally, an s-wave gap is often taken to be any dependence A (A;) that has the full symmetry of the lattice. T h e high T superconductors were found to have a d-wave gap c symmetry: A oc k% — ky. This fact is reflected in many of their temperature dependent properties. Some of the first evidence for the d-wave pairing symmetry was observed Chapter 1. Introduction and Motivation 2 by Hardy et al. at U B C [?], where measurements of A ( T ) at microwave frequencies showed a linear temperature dependence. T h e linear temperature dependence of A indicated a node in the gap function, a direction in which the gap is zero. Although nodes can exist in a gap with anisotropic s-wave symmetry, it is rather unlikely. In addition, measurements sensitive to the phase of the gap confirm that indeed the gap has predominantly d-wave symmetry. Another interesting feature of the high T c superconductors is their quasi two di- mensional structure where most of the conduction takes place in the copper-oxygen (CuC>2) planes. T h e carriers in the planes are believed to be holes rather than electrons and the number of holes per unit cell, x, is controlled by the structure of the material outside the C u - 0 planes. For example, in Y B C O , where the O3 sites form chains in the b direction, the overall oxygen content determines the average chain length which in turn affects the number of supercarriers in the planes. In 1989, Uemura et al. demonstrated a universal relation between the superfluid density, n , and s T. c It was found, in several high T superconductors and for a range of doping levels, c that T varies linearly with n c s up to a certain optimal density at which T c saturates and then is suppressed [3]. These results were obtained through muon spin rotation (/xSR) experiments, where the muon spin relaxation rate, a, which is proportional to 1 / A , was measured. Its limit at zero temperature, a(T 2 —> 0), is proportional to the carrier density at zero temperature and a sharp drop in <r(T) indicates T . c universal doping dependence of T This reflects certain intrinsic properties of the C u - 0 c planes. In the last decade a considerable amount of effort has gone towards measuring microwave properties of the H i T c superconductors. A wide variety of microwave data is available on both the real part of the conductivity and the London penetration depth at different frequencies and temperatures (see for example [4, 5, 6]). However, the limit as T —> 0 is still not well tested. In particular, an absolute measure of Xi(T) is challenging to measure. A n y volume exclusion technique that aims to measure Ax, to within 10% would require an accurate measurement of the sample thickness to better than 1/1000 and very precise control of demagnetization effects (see [7] for more details). In our lab, A Ax, C O of high purity YBaiCu^O-j-i, ( Y B C O ) is measured through a cavity perturbation method that is described in appendix A . T h e difficulty of measuring the sample thickness restricts the penetration depth to be measured only up to an additive constant, i.e the temperature dependence AAx,(T) = Ax, C O — Ax,(T& ) ase is resolved, but an absolute value of Ax(7], ) cannot be obtained. Therefore some ase Chapter 1. Introduction and Motivation 3 other method is required. Moreover, one cannot extract the complex conductivity a(u, T) from the surface impedance without knowing A ^ T ^ e ) . In order to complete the data, a value of A ( T = 0) is inferred either from fxSR experiments [8, 9] or from infrared reflection [10]. Values obtained from far infrared spectroscopy, for example, are A a = 1600A and A = 1000A for b YBa^Cu^OQ.^. However, the infrared spec- troscopy method is extremely demanding and cannot practically be used for a variety of samples. It should be mentioned that pSR experiments determine A i n the vortex state as the length scale over which the magnetic field decays from the normal core to the superconducting bulk. It is now known that the effective A is field dependent, which is interesting i n itself, but this means that the to the Meissner state A. In addition pSR fiSR values do not correspond directly measurements are limited to a geometrical average between A and A&. O n the other hand, extracting A from infrared reflection a involves performing a type of Kramers Kronig integral of the data from zero to infinity frequency. In practice, the measurements exist only over a finite range of frequencies and some extrapolation is required. This extrapolation reduces the accuracy and sometimes underestimates the importance of low frequency spectral weight. The purpose of the new technique described here is to provide an absolute measurement of XhiTbase) at a l ° w base temperature. T h i s method involves preparing a superconducting sample ( Y B C O ) that is doped with magnetic ions (Gadolinium) and therefore cannot be applied to existing samples. However, we assume that the doped samples have the same penetration depth as undoped Y B C O crystals with the same oxygen concentration. T h i s assumption is based on the low concentration of G d ions in the sample and on the detailed chemistry of G d substitution in Y B C O . In particular the Gadolinium substitutes only for Y t t r i u m , has the same charge and is located at a symmetric position between the CuO% planes and not on the planes. Finally, because the lattice parameters and very different from Y B C O T for GdBa Cus0 s c 2 7 are known and not {a = 3.92A, b = 3.85A, c = 11.75A and T = 94K [?]), it c is clear that G d is a rather weak perturbation on Y B C O . T h i s thesis contains some introductory material on superconductivity, microwave spectroscopy and E S R experiments (chapters 2 - 4 ) , description of the new method and its application to Y B C O (chapters 5 and 6), and finally the results are given and discussed in chapters 7 and 8. Chapter 2. Conductivity and Surface Impedance 4 2: Conductivity and Surface Impedance The theory of high T superconductors is yet to be established. Nevertheless, some c simple phenomenological descriptions of superconductivity are available and can be applied to high T materials. c 2.1 The Two Fluid Model The two fluid model of F. and H . London provides a basic picture of superconductivity without assuming a specific mechanism that leads to superconductivity. Any theory of the mechanism (such as the B C S theory) can be combined with this model to complete the picture. The two fluid model divides the charge carriers (originally thought of as electrons) into two groups. One group contains the normal carriers which do not take part in the superconductivity. The other component, called the superfluid, contains the supercarriers, gapped excitations (cooper pairs in the B C S language) that condense below T and exhibit frictionless motion. The supercarriers, c also referred to as the condensate, are responsible for the superconductivity while the normal carriers exhibit scattering processes. The two groups may exchange particles as long as the total spectral weight ne /m* (where n is the carrier density and m* is 2 the effective mass) is conserved. We shall follow the description of this model by Tinkham [11]. Let us denote the fraction of carriers that are normal by x and the fractional number of supercarriers by (1 — x), so that the total number of carriers is conserved. The normal fraction varies with temperature with x —> 1 above T and x —t 0 at T — > 0. c 2.2 Conductivity in the Two Fluid Model Within the framework of the two fluid model, let us use the simplest transport picture, the classical Drude model. In the Drude model under the influence of an applied electric field E, the charge carriers travel uninterrupted for an average time r. The free motion ends when a collision of an unspecified nature occurs and the carrier is scatted. The equation of motion of a charge carrier with effective mass m* and charge Chapter 2. Conductivity and Surface Impedance 5 —e i s : 771* r m*r + = -eE. (2.1) Assuming an A C field with angular frequency u and time dependence of e * leads to iw a complex conductivity a that is defined as J = aE and is given by: , . N <T(UJ) = a x - ne r/m* ne r/m* . 2 %a 2 = — . 2 = 1 1 + IUJT 1 - 2 1 + IUJT) „ 2.2 U> T 2 2 where n is the total number of carriers. In a normal metal at microwave frequencies and below, OJT <C 1 and o is effectively real: a « cri = A simple picture of a superconductor is obtained by taking the limit r —>• oo. In this case, the conductivity away from u — 0 is primarily imaginary and the real part <TI shrinks to a delta function at zero frequency. a>) = ^ W ) - i ^ (2.3) TO* TO*CJ where n becomes n , the number of supercarriers. In the presence of both super and s normal fluids the conductivity is the following sum: o = o~ -\- o~ = o~\ ia s n a.iu) = ^5{u) m* + 2 s TO*O; n TO* — T TO* , . ne n e aao; = - £ - + " (2.4) 2 1 + 2 UJ T 2 o;r 1+ w r (2.5) 2 2 2 2 2.6 Using (for example) the phenomenological Gorter/Casimir T dependence of x: x = ( ^ ) , leads to negligible contribution of normal carriers to the conductivity at T « 4 0.2T . Therefore at any finite low frequency, o w — ia . Figure 2.1 shows o\ and cr , C 2 2 schematically, at different temperatures. 2.3 Surface Impedance Microwave properties of superconductors such as the surface impedance can shed useful light on the intrinsic properties of a superconducting sample. Given that the shielding currents are confined to a very thin surface region, it is convenient to introduce the surface impedance, Z . It is defined by the ratio between tangential s magnetic and electric fields at the sample's surface and can be measured in microwave cavities (see for example [4, 5, 6, 12]). When the magnetic field is applied in the x Chapter 2. Conductivity and Surface Impedance T Figure 2.1: >T 6 C A schematic picture of the complex conductivity of a superconductor as described by the two fluid model in different temperature regimes. T h e upper figure 2 shows the normal limit where the lifetime r is finite and a = " , e / T / 2 m = ? e / y^il — iur). T h e middle figure shows the T —¥ 0 limit where the real part shrinks to a ^-function around u = 0 and the imaginary part is proportional to K T h e small peak in o\ at high frequency corresponds to excitations above the gap. T h e bottom figure shows an intermediate temperature, in the superconducting state. Chapter 2. Conductivity and Surface Impedance 7 direction, the electric field is induced in the y direction. T h e surface impedance is defined as: Z = ^ = R + iX S S where R is called the surface resistance and X s (2.7) s is the surface reactance. T h e relation s between Z and the conductivity o can be found through the following classical skin s effect calculation ([13]). Consider a semi-infinite conductor (or superconductor) in the Z > 0 space with uniform complex conductivity cr. A n oscillatory external magnetic field that is uniform in the x-y plane and varies in time is applied in the x direction: H(f, t) = H (z)e~ luJt x where taking the real part to obtain the physical field is implicit. T h e electromagnetic fields will obey Maxwell's equations and Ohm's law: BR VxH = J, V-B = 0, V x £ + - J = aE. = 0, (2.8) at Since H is applied in the x direction its solution for the equations above varies only in the z direction. Expressing the above equations in terms of the magnetic field, H, and taking V x V x H = - V ^ dz 2 and dt = -iuH (2.9) x v ' yields: d {-^ + iHo<rw)H {z) = 0 2 (2.10) x with general solution H (z) = H Qe lkz x X k = ipoau. In a normal metal, at 2 where microwave frequencies, a is effectively real so that k = (1 + and H (z) = H^e'l-i. (2.11) x where 8 is the classical skin depth: S = A / — — . 2 However, this is not the case in a superconductor where a is primarily imaginary for low frequencies. If we assume a = —io~2 then k = fiQO2C0 and the field decays as H oexp(—|) with the penetration 2 x depth: A = J—J—. (2.12) T h e electric field E obeys the same equation as H in the y direction and therefore has the same spatial decay. Eliminating the exponents, the equation for V x E is Chapter 2. Conductivity and Surface Impedance reduced to: dE ~d~z y dH . ~^ ~dT ^ 8 T = Q i k E y ~ ^° = i (- ) ujHx 2 13 and the surface impedance follows: Z = — s ] - - = - I — = /-—. ] In general the conductivity is complex: a — a\ — ia 2 * = V 2RT^I) X s i = s and we obtain 2KT^I) In the normal state o\ >• o" and therefore i ? = X 2 (2.14) ] = \J^-- s ( 2 - 1 5 ) In the superconducting state at low frequency <r becomes much larger than o\ as soon as T drops below T 2 c and the surface resistance and reactance reduce to: * = ^ . (,.7) Combining this result with our definition of the penetration depth cr = 2 l/\ u/i 2 0 gives the relation between the surface impedance and A and o\: R = 1 s lp uja 0 l _ 1_ 2 2 x 3 _ = 5 ^ A V _ v X = fJJ^J s = „ uX. 0 (2.18) 0-2 Microwave techniques that measure surface reactance give direct access to the penetration depth and therefore to the superfluid density, as we shall see later on. Measurements of the surface resistance (combined with A) can give the real part of the conductivity Oi. Chapter 3. Penetration depth 9 3: Penetration depth The magnetic penetration depth has been introduced in the previous chapter as the skin depth of a superconductor in a time varying magnetic field. In this chapter we show that the exponential decay of magnetic fields (constant or time dependent) is closely related to superconductivity. In particular, the decay constant, A is closely related to the superfluid density. 3.1 Meissner effect and the London theory The Meissner effect i n conventional superconductors was first observed experimentally by Meissner and Ochsenfeld, in 1933. In 1935 F . & H . London showed that the phenomena of superconductivity is closely related to near-perfect diamagnetism and explained the Meissner effect theoretically. In the Meissner effect, the superconductor expels any applied magnetic field from its volume through a shielding supercurrent traversing its surface. Type I superconductors (such as many elemental superconductors) are materials that have only two phases separated by a critical field H . The c normal state, above H , in which the magnetic field fully penetrates the sample and c the Meissner state where effectively no field penetrates the sample. Type II superconductors (for example many conventional alloy superconductors and the high T c materials) are materials in which a so called vortex state exists. In these materials, below the critical field H i the sample is in the Meissner state and above a higher c critical field H 2 the material is normal. In the region between i f C c l and H 2 there C is a vortex state i n which the field penetrates the sample in the form of vortices. In this state the bulk superconductivity isn't destroyed and the vortices can form a two dimensional lattice. The London theory describes type I superconductors in a magnetic field below H or a type II superconductors in magnetic fields smaller than c Hci- The theory begins with minimizing the free energy of a superconducting sample in an applied magnetic field. We assume that the thermodynamic state that minimizes the free energy includes a steady surface current. The free energy is: •F ~ j F dr + E)~iri ~t~ s ^magnetic (3-1) Chapter 3. Penetration depth 10 where F is the free energy of the charge carriers in the superconducting condensate s [14], E kin is their kinetic energy associated with the constant current and E gnetic ma depends on the magnetic field configuration. Let us define v(f) to be the velocity of the charge carriers and j the resulting supercurrent as: j (r) = n ev(r) where n is s s s s the density of the carriers in the condensate and v(r) is their velocity. The constant flow of currents results in a kinetic energy of: E = / \mJn,ir tm (3.2) and the vector r spans the sample volume. The magnetic energy is given by: ^magnetic =kl - - &if (3 3) The relation between the current j and the magnetic field H is given by Maxwells equation: s V x # = — VxB = J. (3.4) s Mo The total free energy is then: T = / F df + f[B + A | | V x B\ }df J 2fi J 2 2 s (3.5) 0 where the constant A^ is defined by AL = [ - ^ - ] 1 / 2 . (3.6) Minimizing the free energy with respect to the magnetic field B will determine the field configuration inside the sample. To minimize T we vary B —> B + SB. The first order variation in T is: 6F=— Mo f df\B-SB + X y xB-V J 2 L x8B}. (3.7) Using integration by parts we can replace the term V x B- V x SB by — V x V xB-SB — V B • SB to obtain: 2 ST = — fdf\B + V B] SB. Mo J L J 2 (3.8) Chapter 3. Penetration depth 11 Taking 5T = 0 gives the L o n d o n equation: B + \]y B 2 = 0. (3.9) In the above calculation we have shown that the existence of steady currents without driving electric fields leads to the London equation. T h e L o n d o n equation describes perfect diamagnetism as will be shown in a simple example in the next section. 3.2 London Penetration Depth The London equation shows that the magnetic field in the sample decays exponentially into the sample with decay constant XL- T h e detailed form of B depends on the geometry of the sample and the applied field. Let us consider a simple example: an infinite superconductor is placed i n the z > 0 half space, so that its surface is the x-y plane at z = 0. T h e magnetic field is applied in the x direction. T h e magnetic field is uniform on the superconductor's surface and varies only i n the z direction. In this case the London equation reduces to: ^ ) + A i ^ ^ = 0. (3.10) B (z) = B (0)e~ ^ . z T h e solution is an exponential decay: x L Note that the above x calculation assumes the bulk superconductivity is not destroyed by the magnetic field and therefore the theory is valid for weak fields. T h i s simple example can easily be generalized to other sample shapes since A is often small compared to the macroscopic sample dimensions. Except for the Meissner effect itself, another important prediction is the value of XL. We have previously shown that X L = f^*^] ^ ; 1 2 in other words XL gives the ratio of number of carriers per unit volume to their effective mass: n e /m*, also 2 s called the oscillator strength since it is the pre factor of the delta function of o\ in zero frequency. In the limit of low temperatures we expect all charge carriers to be part of the superconducting condensate. In the context of the B S C theory [1] the carriers are pairs of electrons (known as Cooper pairs) and therefore we would expect n s = \n , e the pair's effective mass should be close to 2m e and their charge is 2e. In metals like aluminum one would find A ~ 500A. In the high T cuprates the carriers c are believed to be holes and their concentration depends on oxygen content. In an optimally doped single crystal Y B C O A was measured to be between 800A and 2000A by pSR [8, 9] and far infrared [10]. Chapter 3. Penetration depth 3.3 Magnetic Penetration 12 Depth It is important to mention that the above calculation is only valid in the local limit, also referred to as the London limit. A n assumption has been made here that both the magnetic field B{r) and the current J (f) vary slowly in space. In order to define s what is meant by "slowly" another parameter, £, the coherence length, has to be introduced. So far we have considered a supercurrent that is carried by a condensate without introducing any mechanism for the superconductivity. However, the mechanism is important for determining the relevant length scales in the problem. In 1957 J . Bardeen, L . N . Cooper and J.R. Schriffer laid the fundamentals of the theory of superconductivity known as the B C S theory [1]. The details of the theory will not be given here and can be found in reference [14]. The theory describes superconductivity as the result of pairing between electrons with equal and opposite momentum, the Cooper pairs. The ground state of the system (also called the condensate) is separated from the excited states (quasiparticles) by the energy gap, A . The Hamiltonian of the system contains the pairing potential and the energy of the quasiparticles (above the condensate) is given by: E = y 2 ( £ 0 + |A(£;)P where c(k) = (3.11) e k —ep is the energy of a free electron with respect to the Fermi energy and A(fc) is the energy gap. A t low temperatures the system contains a condensate of paired electrons which exhibit frictionless motion. Excitation of quasiparticles involves pair breaking and therefore requires energy that is larger than the energy gap. This reduces substantially the probability of collisions since at low temperatures such energetic particles are not available. The gap A may vary in space and the coherence length, £, is defined as the length scale in which the gap can vary. A n estimation of the coherence length can be made in the following way. The relevant kinetic energy range of our system is of the order 2A around the Fermi surface: E - A < F 0 <E +A F 0 (3.12) where Ao is the maximal gap in case of non uniform gap function. In this range of energies we can approximate k « hp + 5k and E the relevant momentum range is 8k = Ep + Vph5k and therefore This momentum span corresponds to a Chapter 3. Penetration depth 13 spatial span of Sx ~ ^ and the coherence length is approximately: hvF T h e London limit holds when A > > £ . In this limit the correlation between electrons has a short length scale and their response to electromagnetic fields is local. In high T superconductors such as Y B C O , £ is of the order of several lattice constants and A c is a few thousand A so that the London limit is valid. In other cases (elemental metals for examples) £ may be considerably larger. In this case the free energy calculation should be modified to introduce a non-linear relation between J and A. Chapter 4. Electron Spin Resonance 14 4: Electron Spin Resonance Our measurements of A are obtained by applying an rf field to magnetically doped Y B C O and observing Electron Spin Resonance ( E S R ) transitions in a dilute array of magnetic moments. T h i s chapter will describe briefly the basic theory of magnetic resonance in general and electron spin resonance in particular. 4.1 Magnetic Resonance A n y physical system with quantized magnetic moments can exhibit magnetic resonance. O u r interest is in the spin magnetic moment of the electrons in the ion. However, the following description of the magnetic resonance effect is general, where the magnetic moment gu-sS can be replaced by g/^sJ where J is the total angular momentum in the system, J — L + S +1 where L is the angular momentum, S is the electron spin and I is the nuclei spin. Consider a system of free spins with magnetic moment PBQSmagnetic field all energy levels are degenerate. A p p l y i n g a uniform magnetic field in the z direction lifts the degeneracy completely. T h e system's new eigenstates correspond to the Zeeman levels with energy —fj,BgmH projection of the spin on the z axis. |m) are the eigenstates of S , z Without any external where m = —S...S is the T h e new eigenstates of the system, labeled where S \m) z = m\m). A n y magnetic interaction, such as the crystal field, may lift the degeneracy in some way that would change the eigenstates of the system (for example, the crystal field lift the degeneracy of different spin orientations but does not split the Kramers doublets). In general, the new eigenstates can be described as linear combinations of the \m) states. s' (4.1) :=-S 771= For simplicity we will assume that the energy levels \xb ) are very close to the \m) n states, and will use the \m) states instead. T h i s assumption is valid for the G d ions in Y B C O but is not at all necessary for the following discussion and the generalization is straightforward. In any case, the system has quantized energy levels with lower degeneracy than without any external field. A t zero temperature only the lowest en- Chapter 4. Electron Spin Resonance 15 Figure 4.1: E S R transitions in a spin system. The degeneracy between the different spin orientation was lifted due to an applied magnetic field ergy states are populated and at finite temperature each state is populated according to the Boltzmann distribution: N . e 3 = Y^<r EilkBT _ -E /k T l = B z (4.2) i where Nj is the probability of finding an ion in the j th state with energy Ej and z is the canonical partition function. Magnetic Resonance occurs when the system absorbs or emits an energy quantum (a photon) that is equal to the spacing between two energy levels, provided the transition is allowed by quantum mechanics. A photon can be fed into the system by an alternating magnetic field of the form Hf r = Hocos(ujt). This field is described by an additional perturbative term in the Hamiltonian: r i = r i o + iiB9S-H (4.3) rf where T-LQ describes the Hamiltonian of the system without applying external A C fields and HsgS is the magnetic moment). When the field is perpendicular to the principal spin axis, this term can be written as a linear combination of spin ladder operators. For example when H is applied in the x direction and the principal spin axis is z we can write: S-H = H ^Re(S e iut rf x0 + + S-e^) (4.4) where S ± = S ± iS . The new eigenstates of the system can be found from time x y dependent perturbation theory and can be written as a linear combination of the Chapter 4. Electron Spin Resonance 16 eigenstates of UQ with time dependent coefficients: m)) = ^2o (t)e- H \m) iE (4.5) kt m m where E m are the eigen energies of the unperturbed Hamiltonian, 'Holm) = E \m). m Perturbation theory shows that each state \<f>) would oscillate between two or more \m) states. T h e transition rate (number of transitions per second) is given by Fermi's golden rule: irH 2 where f{u) is a shape function normalized to unity and concentrated around u ~(Ei — Ej). = T h e shape function is a result of spin relaxation processes (spin-spin interactions a n d / o r spin-lattice interactions) and is usually a Lorentzian. fiij is defined as pij = gfJ,B(mj\S y\mi) (x or y depend on the direction of the applied rf field) and Xt provides the selection rules for allowed transitions. 4.2 Magnetic Ions in a Crystal: Fine Structure A magnetic ion i n a crystal experiences magnetic interactions via two mechanisms: the motion of charged particles around it (mainly electrons) and the magnetic moments (associated with the spin degree of freedom) of electrons and nuclei in the crystal. The first interaction is the fine structure of the crystal and the second one is the hyperfine interaction. O u r subject of interest is a magnetic ion substitution in an otherwise non-magnetic crystal, therefore we will only consider the fine structure. A general Hamiltonian for describing this magnetic system is: u = E t fc=l ^ - T ~ k + i { r k ) l k ' § k + v { n ) } + ( k<j 4 7 ) j k where n is the number of electrons associated with the ion, Pk is the fcth electron's momentum, describes the strength of the spin-orbit coupling of the fcth electron and V(rfc) is the potential energy of the fcth electron in the electrostatic field of the environment. T h e second sum is the Coulomb interaction between electrons. Notice that in this picture the electrons are divided into two groups. T h e first contains the electrons that belong to the ion, which will be treated quantum mechanically, and the second contains environment electrons, that are expressed in V(rk) and are treated classically. Further justification and details of this picture can be found in Griffith Chapter 4. Electron Spin Resonance 17 1961 [15]. 4.3 Effective Spin Hamiltonian The above description is a good approximation for the magnetic ion in the crystal but the form of the potential V(rk) is unknown and may have a very complicated spatial dependence. The symmetry of this Hamiltonian is very useful for further analysis. V(rk) and hence Ho commutes with all the elements of the point group of the ion's site in the crystal. The effective spin Hamiltonian expresses the energy levels of the system in terms of spin operators. It is equivalent to the full Hamiltonian of the system within the subspace of a given electronic shell (constant radial wave function R (r) and total n angular momentum). It has the compact form: % ( - ) = 4 8 PA where O is a Stevens operator of order q. The Stevens operators are functions of v q S ,S ,S+ 2 Z and 5 _ and are the operator equivalents of the spherical harmonic func- tions. The B are numeric coefficients that are called crystal field (CF) parameters. v q This sum is an expansion of the crystal field in the ion's environment in terms of a complete set of functions. We shall see that the number of elements in the sum is finite and depends on the symmetry of the crystal and the magnetic moment of the ion. The first step in building an effective Hamiltonian is to prove that the system can be described by a finite number of parameters. One can then choose an appropri- ate method for calculating them theoretically, but typically the C F parameters are measured through E S R experiments. The theoretical effective spin Hamiltonian with measured C F parameters can predict the magnetic field and frequency dependence of the E S R transitions as well as their intensities and show very good agreement to the experiments. The full proof for the C F theory is given in [15] and will not be repeated here. The basic scheme is to describe the different ionic states as Slater determinants of oneelectron states and impose the site symmetry on the Hamiltonian matrix elements. The C F Hamiltonian can then be expanded in terms of spherical harmonics. Its matrix elements between ionic states are then written in terms of matrix elements of spherical harmonics between the one-electron functions, that are also described Chapter 4. Electron Spin Resonance 18 by spherical harmonic functions. Group theory arguments show that the number of non-zero matrix elements is limited by j where j is summed over all irreducible representations of the site point group and rij is the number of times this representation occurs in the relevant representation (according to the magnetic moment). The number n is also the number of free parameters in our problem and therefore the number of C F parameters in the effective spin Hamiltonian. 4-3.1 Semi-Classical Approach for Deriving the Effective Spin Hamiltonian Let us assume that the crystal environment of the magnetic crystal is represented by a static charge distribution p. The potential energy V(r) is given by: ( -io) 4 n o = - / nr^<*«. J I it — f] where r is the position of the electron and R is a general point in the environment. We define the environment as the space outside the magnetic ion and therefore if both r and R are measured from the nuclei then |f| < \R\. We can expand the potential in Legendre polynomials: J \R- k 0 0 ^ = E^+r *( p C0S («)) (4.ii) where Q is the angle between R and f. We proceed by expanding each P further in k terms of spherical harmonics. It is convenient however to use a different basis set, Z , instead of the regular spherical harmonics Y : im lm Zio — Yio, Z[ = y2(Y _ m l m + Y ), lm (4.12) ^Im ~ f \ / 2 ( V i _ — Yim) m The above functions are all real, with Z and Z being functions of cos(ra0) and c sin(ra</>) s respectively. The vector space spanned is the same as for the usual spherical Chapter 4. Electron Spin Resonance 19 harmonics. The expansion of the Legendre polynomials in terms of Zi is: m P (cos(Q)) = Z (l)Z (2) k ka (4.13) ka a where (1) refers to the coordinates (9,4>) of the general point R and (2) refers to (9, 4>) of the electron position r. The sum over a is over all the Z\ m I including Zio and Z\ s m functions of the same for m ^ 0. Inserting this expansion into our expression for V(r) (by performing the integration over R) we obtain the potential V in terms of the electron coordinates, (r, #,</>): oo k=0 a where ^ k o i are characteristic numbers of the environment and are defined by: J k a ~ 2k + lJ R^ d K ^- ) 4 15 Since the ionic wave functions are Slater determinants of one-electron wave functions (the spherical harmonics \l,m) times a radial function) all matrix elements of V{r) are of the form: M = (Z, m\Z \l, m) times the average radius (f). The average ka radius is only a numerical constant since the effective Hamiltonian is defined within the subspace of the same radial dependence. Different I functions are not included since we consider only electrons outside of the core, within one unfilled shell. The 4> dependence of the spherical harmonics shows that the matrix element M is nonzero only when k < 21. The parity operator is a member any site symmetry group and represents time reversal symmetry. This imposes an additional condition, that k must be even. This reduction of the number of C F parameters can be applied to every crystal. More symmetry in the system reduces the number of C F parameters even further. Any Z ka function that violates the one or more symmetries of the system do not appear in the expansion. The above matrix elements depend only on the angular momentum degrees of freedom and can be described as equivalent operators (the Stevens operators) whose derivation can be found in [16]. A full list of the Stevens operators up to the sixth order can be found in [17]. It is important to mention that even if the semi-classical picture that is given here is not accurate, the expansion of V(r) in Stevens operators is still justified since it is an expansion in an orthonormal basis that spans the chosen subspace. The cancelation Chapter 4. Electron Spin Resonance 20 of many terms due to symmetry and orthogonality remains valid since the symmetry of the site can be imposed on the elements in the sum. However, the coefficient of each operator in the expansion may very well be different from what is calculated i n this simple picture. The C F parameters can be corrected in several theoretical ways [18] and through experimental measurements (such as the work described here). 4.4 Zero Held ESR and the ESR Spectrum This section will discus the principle of E S R experiments i n general with emphasis on zero field E S R . Most E S R experiments are performed with H f at a constant r frequency so that E S R transitions can only occur at the relevant energy hu. In other words, the perturbation to the Hamiltonian is proportional to Re(S e luJt + + S-e~ ), llJt where ui is constant. A t this frequency transitions will occur only if AE = E -Ei = hw. (4.16) j To observe the various transitions between different pairs of states the energy levels must be shifted so that the difference AE matches the applied frequency. This is done by applying a D C magnetic field, usually parallel to the principal axis of the ion's spin. In this way each energy level E , that corresponds to a state \m) (an m eigenstate of the S operator or a mixture of states that is dominated by \m)) will z gain an additional energy of — PBQ^H. This field would also lift the degeneracy associated with time reversal symmetry (the Kramers doublets of| ±TO)).The full Hamiltonian is then: H = H 0 + g/jLBH^ReiS+e™ + S_e~^) + gp H' S . B z z (4.17) The magnitude of the field H' can be changed so that each transition will occur at Z the frequency u but at different H' field as shown in figure 4.2 for a spin | system. Z Another way to observe the transitions is to take the zero field E S R approach. In this situation the experiment is performed without any applied parallelfieldH ' . Z In order to see all the possible transitions, it is necessary to scan a wide range of frequencies. The observable quantity is the absorbed power in the E S R transitions or the imaginary part of the magnetic susceptibility, x" • These quantities can be inferred from the energy change rate (from Fermi's golden rule). ^ = WijiNi - Ni)hw = ^ M V M A M l - e~^)hw (4.18) Chapter 4. Electron Spin Resonance 21 H ' Figure 4.2: E S R transitions at a constant frequency and varying D C magnetic field for a spin | system. The splitting in H' = 0 is the fine structure (where | ± ) states are degenerate. When H' is varied, the levels evolve according to their magnetic moment g(j, m. z z B 22 Chapter 4. Electron Spin Resonance 2 4 6 Frequency ( B ° ) Figure 4.3: Theoretical x"i ) f ° system with spin | and fine structure Hamiltonian Ho = B§0g = B%(3S - S ). The allowed transitions, ± ± ± § , ± | -> ± | and ± | —>• ± | , occur at 2B , 4B and 6 5 ° respectively. A Lorentzian shape function of width 0.06B has been assumed u 2 ra 2 Z 2 2 2 where Wij was defined in equation 4.6. From the relation ^p- = we get the \H OJX" 2 Q susceptibility: = ^ k l V ( ^ ( l - e"^) (4.19) where the average number of spins in the i t h level, iVj, is given by the MaxwellBoltzmann distribution: Ni = N-ez (4.20) Ei/kBT where N is the total number of spins per unit volume and z = partition function. The total x" 1S a s u m e~ / is the Ej k£lT of shape functions, centered around the frequencies of allowed E S R transitions with an integrated intensities of | ( / / j j | A ( l — 2 7 i e B ). Figure 4.3 shows a typical susceptibility spectrum for a Hamiltonian "H = k T 0 B0 2 = 3S — S and spin | . In a \ spin system the energy levels are E 2 2 2 m = B (m — | • | ) and the E S R transitions appear at hh) = E +i — E = B (2m + 1) 2 2 i.e at 2B%, 4B% and 6B°. m m 2 Chapter 5. The Bolometric Experiment 23 5: The Bolometric Experiment Performing zero field E S R on magnetic ions embedded in a superconducting sample involves measuring )d'(cu) at a low temperature over a wide range of frequencies. In the next chapter it will be shown that the susceptibility and the surface resistance are related through: X » where Rf SR = T ^ T - (5-1) is the contribution of the magnetic impurities to the sample's surface resistance and AL is the London penetration depth. This relation was used to extract the susceptibility from surface resistance measurements. The experimental part of this project was performed by Patrick J. Turner in an existing bolometric apparatus [19] which has been developed in the U B C superconductivity lab as a sensitive probe of the microwave surface resistance of high T samples. Measurements can be carried out c over a range of temperatures (1.3-10K)and over a broad band of microwave frequencies (0.1-21GHz). The choice of performing zero field E S R through a broadband measurement of surface resistance was made for several reasons. The first, practical, reason was the existing surface impedance measurement facilities in our lab. A n early idea was to measure the real part of the susceptibility, x!{T), through surface reactance measurement at very low frequency (10kHz) where the imaginary part, x", is negligible. This would be essentially a measurement of the D C susceptibility of the spins that are exposed to the magnetic field and therefore the magnitude should be proportional to the penetration depth. The effect of magnetic impurities on R (T) is an upturn S at low temperatures since X oc y/\ + x' • This upturn is a result of the 1/T paras magnetic temperature dependence of x'(T). However, the first order contribution to the D C response comes only from the —\—*+\ transition between the spin levels of the ion, while the higher spin transitions only contribute in the second order to the susceptibility. D C measurement of x' would show the expected upturn, but the G d contribution to this effect is actually quite small and might not be easily distinguished from other magnetic contributions (such as BaCuO flux on the sample's surface). A n alternative approach was therefore taken. It was a fortunate coincidence that an apparatus existed in the lab which is capable Chapter 5. The Bolometric Experiment 24 of measuring the frequency dependence of R over the range where the zero field E S R s peaks exist. The A C response of the sample, x"{ )i involves first order spin transitions UJ at certain frequencies (the ESR transitions). The frequency of the transitions together with their intensity provides the information needed to determine the number of ions exposed to the field (see chapter 6). The bolometric experiment was found to have both sufficient sensitivity and appropriate range of frequencies to reveal the ESR lines for this particular system. Our knowledge of the relevant frequency range was based on the crystal field parameters as measured in conventional E S R [20, 21, 22, 23, 24]. Furthermore, zero field E S R (rather than conventional ESR) is done in the same field configuration as other microwave measurements at U B C . The samples are usually very thin in the c direction compared with the dimensions of the a-b plane. The field is applied along the a (b) direction with the currents running mainly in the a-b plane. In this configuration the magnetic field is hardly distorted by the sample. On the other hand, conventional E S R must be done with external magnetic field along the c direction. In this configuration the field is distorted by the relatively large area of the a-b plane and the resulting demagnetization effects are much larger. A n additional advantage of the chosen method is the ability to isolate the E S R lines from other features in our data. The susceptibility x"( ) is extracted from u the measured surface resistance as explained in section 6.3. The first step in this extraction is to separate the G d contribution to the surface resistance, Rf (u>), SR from the superconductor surface resistance, i?f (u>). For Y B C O , Rg (u>) is known c C to be a smooth monotonic function of frequency, that can be fitted by a low order polynomial function. A n example for this functional form is shown in figure 5.2. On the other hand, Rf SR is a sum of sharp Lorentzians with widths of about 250MHz. The total surface resistance, as measured by the bolometric experiment is the sum of the two, but the two components can be easily separated thanks to their different shape. 5.1 The Measurement Principle Microwave measurement of the surface resistance of Y B C O samples have been performed in the U B C superconductivity lab since 1989. The typical measurement technique is cavity perturbation (this method is described in A). A resonator is designed to have a resonance frequency with high quality factor (Q). The sample's surface impedance, Z , is measured at this particular frequency over a wide range of tempers atures. Different resonators in different frequencies were designed, and R of various s Chapter 5. The Bolometric Experiment side view 25 end view Figure 5.1: The electromagnetic fields configuration in the bolometric apparatus. The sample and reference are the black rectangles Y B C O samples were obtained [7, 4]. The bolometric technique was developed to complete the resonator data over a broad band of microwave frequencies in order to show the detailed frequency dependence of the surface resistance. The basic idea behind the bolometric technique is to measure the sample's surface resistance through its power absorption in a microwave magnetic field. This is done using a reference sample (in our case a high purity Ag:Au alloy) of known surface resistance. The sample and the reference are inserted into a shorted transmission line in two positions that have nominally the same magnetic field configuration. The thermal response of the reference and the sample due to their rf absorption are measured as a function of microwave frequency. As shown in appendix C, since the two samples experience the same tangential magnetic field on their surface, it follows that the ratio of power absorption is also the ratio between their surface resistances: _ psample psample (5.2) s = preference R f Te ' ' erence Using the power ratio, the sample's R is inferred from the known R of the reference. s s Implementing this measurement requires a setup with a microwave magnetic field that is uniform over the dimensions of both samples, is equal in the positions of the sample and the reference, and whose frequency can be varied continuously . The experiment is designed as a rectangular termination of a coaxial transmission line. The field configuration near the end of the line is shown in figure 5.1. A t the termination, the center conductor is a broad flat rectangle. This provides a symmetric Chapter 5. The Bolometric Experiment Figure 5.2: 26 A typical surface resistance spectrum of a Y B C O sample at different temperatures as measured in the bolometric apparatus magnetic field that is reasonably uniform in the regions of the sample and reference. T h i s design supports a non-resonant T E M mode with a magnetic antinode and an electric node near the sample and the reference. T h e sample holder is designed as a cold finger, a technique that was first described by Sridhar and Kennedy [25]. This technique provides an isothermal connection (sapphire plate) between the sample and the thermometer and a weak link to a pumped helium bath. T h i s permits temperature oscillations of the sample, while ensuring that the time constant for the stage to cool, following power input, is not too long [19]. Chapter 6. X Extraction Technique 27 6: A Extraction Technique Having introduced the bolometric measurement technique and the concept of ESR we are now ready to combine these two types of experiments into a new direct method of measuring the London penetration depth of superconductors. The main idea here is to dope the superconductor with magnetic ions that serve as magnetic indicators. Ions that reside close to the sample's surface, within the penetration depth, are exposed to the magnetic field. Only these ions participate in the E S R process and therefore contribute to the line intensity of the various transitions. Ions deep inside the superconductor experience no magnetic field and therefore will not participate in the ESR process. The resulting frequency-integrated intensity of each ESR transition is proportional to the effective number of spins that are exposed to the field which can then be related to the penetration depth A. The following sections describe how this principle was applied to Y B C O with G d impurities. 6.1 The Samples The material of interest is YBa2CuzOrs- The magnetic dopant chosen for this study was gadolinium (Gd) in the yttrium site, so that the doped crystal is Gd Yi^ Ba2Cu^0 x x 7 This substitution is ideal for several reasons. First, the yttrium site lies halfway between the copper-oxygen bilayer (see the unit cell of Y B C O in figure 6.1). As indicated by the relatively low conductivity along the c direction [26], transport is primarily in the copper-oxygen planes. Choosing a dopant site outside of the planes reduces the effect of the substitution on the scattering mechanism in the d and b directions. We expect the extracted A to be similar to that of an undoped superconductor. Gadolinium is a rare-earth ion with an oxidation number of +3 (the same as yttrium) and substitutes only on the yttrium site. Another advantage of the G d ion is its E S R spectrum. The Gd +3 ion has an exactly half filled f-shell: the / orbital has angular momentum L = 3 and therefore 2 x (2 x 3 + 1) = 14 states. The shell is half filled by the seven remaining electrons (after 3 electrons have left the atom). Such a half-filled shell has total orbital angular momentum of zero. As a result, the ion's magnetic moment comes solely from the spin of the seven electrons and the system is a nearly free spin | system. Thus, the total angular momentum is S — \ and there is only one Chapter 6. A Extraction Technique 28 angular momentum multiplet, instead of the usual J = L + S.. .\L — S\ multiplets. This simplifies the E S R spectrum which has only three main lines, corresponding to the transitions: ± | —> ± | , ± | -» ± | and ± | —V ± | . Moreover, the splitting of the levels (resulting from the fine structure in the Y / G d site) are of order 1-15 GHz so that all of the transitions fall within the frequency range of our bolometric apparatus. Our present study has benefitted considerably from previous work on G d doped Y B C O by A . Janossy et al. [21, 22, 23], S.Pekker et al. [27] and A . Rockenbauerei al. [24]. Their work was conventional E S R in applied fields of aligned Gd Y^ Ba2Cuz,Oi-s x x powder samples. These experiments measured the C F parameters of the G d site for several oxygen contents in a variety of underdoped samples. In this work, we use the conventional E S R measured C F parameters as first guesses for our spin Hamiltonian. 6.1.1 Sample Preparation The gadolinium doped samples were prepared by Dr. Ruixing Liang of the U B C Superconductivity Lab. The procedure is very similar to the one used for making high purity Y B C O crystals [28, 29]. Growth of Y B C O crystals typically begins with powders of I2O3, BaCOs, and CuO. In the case of gadolinium substitution, a gadolinium-oxide powder is added (Gd^Oz). The ratio of gadolinium-oxide to yttrium-oxide sets the nominal Gd concentration. However, the precise concentration in the crystal has to be determined by other means. The rest of the process is identical to Y B C O single crystal preparation. A ceramic crucible (BaZrOz) filled with a mixture of the above powders is placed in a furnace and then heated up to 1000°C for 12-24 hours until the powders melt completely. The heating is done near the edge of the furnace and next to a quartz rod that serves as a cold finger to create a horizontal temperature gradient towards the center of the furnace. The temperature gradient, of 2 — 5°C/cm, improves the growth stability and enhances the growth rate. Once a uniform melt is obtained, the temperature is lowered slowly and at about 985°C crystals start to form near the coldest spot on the crucible. Post annealing is done slightly below growth temperature in pure oxygen flow to release residual strain and defects. A lower temperature anneal is used to set the oxygen content, which is a function of the annealing temperature and oxygen partial pressure [30]. The crystals grow naturally twinned so that the a and b directions are mixed (the Chapter 6. A Extraction Technique 29 Figure 6.1: The Y B C O unit cell. Gadolinium substitutes into the yttrium site. The left picture shows one unit cell and the right picture focuses on the C u - 0 chains along the b axis. Chapter 6. X Extraction Technique 30 so-called twins). The twins are removed by heating the sample to 200 — 400°C and carefully applying pressure along one direction in the a-b plane. This procedure is called detwinning. Another important issue is the crucible material. The BaO-CuO melt is highly corrosive and reacts with most materials typically used for crucibles. Extensive study has led to the currently used BaZrOs crucibles. These crucibles replaced the widely used Y S Z (Y2O3 stabilized Zr0 ) to create the high purity new generation crystals 2 of 99.995% atomic purity [?, ?]. 6.1.2 Gadolinium Concentration Determination The G d concentration in the crystal is not necessarily equal to the concentration i n the starting oxide powders and has been found that the actual G d concentration can be substantially higher than in the starting powders. This difference in concentration is a result of growing the crystal from a liquid melt. The G d ions migrate preferably into the crystal rather than the melt, and as a result the G d concentrations in the crystal is higher than initially expected. One way to determine the G d concentration of the samples is by the Curie-Weiss susceptibility in the normal state using a SQUID magnetometer. The samples were measured in a field of 2T applied parallel to the a (b) direction and the magnetization was measured for temperatures between 100 and 300K, above the critical temperature of the superconductor, T . In this field direction and temperature range, the c susceptibility of Y B C O is constant. The resulting magnetization is mainly due to the G d paramagnetic contribution to the D C susceptibility and therefore obeys the Curie-Weiss law [31]: •x{T)= +e 0 J- (6.1) Weiss where C is the Curie constant and is given by: C = NS(S i y + )g B 3k B and 9 is the constant Y B C O susceptibility. The data is fitted to equation 6.1 and the constant C is determined. A ceramic powder sample of Gd Y\- Ba Cuz07-s x x 2 with known x = 0.01 was used as a concentration standard. The concentration of Gd in the powder is known since its growth procedure does not include separating the crystal from the melt in any stage (thus all of the input materials end up in the sample and nothing in wasted). It was found that Tweiss = —0.5 ± 2K and g = 2.03 ± 0.01. Chapter 6. A Extraction Technique 6.2 31 R and x" s In a superconducting sample with magnetic doping the observed surface resistance has two components: one comes from the lossy part of the diamagnetic shielding current that is assumed to be independent of the magnetic doping (for small doping) and is a monotonic smooth function of frequency (as shown for a typical Y B C O sample in figure 5.2 in chapter 5). The other contribution comes from the transitions between the G d spin levels (the so-called E S R transitions). The magnetic susceptibility x" and the E S R component of the surface resistance Rf are related in the following SR way: Dtotal " = th x TDSC \HQUJA = r>ESR fh (6.3) \PQUIX where Rl° is the measured surface resistance and R resistance. This can be derived from the expression for surface impedance Z = tal is the superconductor's surface s c s s R + iX , as is now shown. A t low temperature (in the superconducting state) s s R <§C X and we can make the following approximation: s s Z =R 2 X + 2iX R 2 S 2 S s « -X + 2iX R . 2 s s (6.4) s On the other hand we can modify equation 2.14 to include the magnetic response of the G d ions: Z, = { . %iUj) )* with p. = /xo(l + X' ~ ix") (6.5) where the permeability p replaces po to account for the magnetic moments in the sample. In the superconducting state o% ^> o~\ at any low frequency, and therefore Z reduces to 2 Z] = ^ [ - ( 1 + X>2 + i ( ( l + X>i + x V ) ] . (6.6) 2 2 a Comparing equations 6.4 and 6.5 we obtain: Aj^a+X-) and R, = f ^ ± . [ i l + X > 1 + X „ a 2 l (6.7) For a dilute array of paramagnetic moments we can assume x ' < l - In our system, a rough estimation of x' was done in the following way. The real and imaginary parts of the susceptibility are related through Kramers-Kronig relation. Since the x " ( ) w depend on frequency only in the shape function, x'(u;) for each transition would be the integrated intensity of x" times the Kramers-Kronig transform of the shape function. For a Lorentzian shape function with width n, the maximum of the transformed Chapter 6. A Extraction Technique function is l/2r). Therefore, the maximal value of x' 32 1 S : xV)maz = (6.8) where A is the integrated intensity of the largest E S R transition. In our system this gives x' ~ 10~ and therefore we can approximate 1 + x' ~ 1 in the above expressions. 7 Finally, with a — 2 M O ^2 w e obtain: R, = + (6.9) and result 6.3 follows. 6.3 Background. Subtraction The extraction of x" from our measured R„ involves subtracting the superconductor's surface resistance (R ) from the total measured surface resistance. This is done by s c s fitting a polynomial function to the surface resistance data between the E S R peaks (at places that are far enough from any Lorentzian tail) and subtracting this fit from the total surface resistance. This fitting procedure was found to be the best solution for separating the two contributions. The functional form of R ° is unknown but it has s s been seen in various samples to be a smooth function of frequency. Another approach would have been to subtract R data from a different sample ( Y B C O without G d s doping) but so far this approach has failed since for the samples tried the background level is sample dependent, not surprising since it is very sensitive to surface quality. The next chapter discusses the results and the data is presented in figures 7.2 and 7.5 as the measured surface resistance with a polynomial curve fitted to the background between the E S R peaks. 6.4 Spectrum Fitting and Extraction of A Once the background is subtracted, the E S R spectrum is ready to be analyzed. We have obtained Rf SR a factor of irp A. 0 = i?* ' — R^ . ota c Dividing this data by frequency gives x" P to U In order to extract A one must first calculate x" from the Hamil- tonian of the system. Our strategy was to use C F parameters that were measured in conventional E S R experiments [24, 27, 21, 22, 23] as a first guess and make some adjustments to fit our own data. To do this a Mathematica notebook was written by Chapter 6. X Extraction Technique 33 Richard Harris (see appendix B). The algorithm receives the crystal field parameters and returns the susceptibility function. The procedure is the following: • Build a matrix presentation of the Hamiltonian given parameters, Tio = ^ in the subspace of spin | . p q BgO^, • Diagonalize %Q to find the eigenvalues which are the system's energy levels. • Calculate all matrix elements of S or S , according to the polarization of the rf x y field. Every matrix element that is non zero is an allowed transition and would be observed as a peak in the spectrum at fko = AE • For every allowed transition, the transition energy is calculated by Eij = Ei — Ej and the intensity of the transition is calculated from equation 4.19 using the eigen values of T-L . 0 • Build the x"^) spectrum as a sum of shape functions (Lorentzians) centered around the calculated transition energies. The width of the peaks is taken from the experiment, but the total intensity was calculated before. • After the spectrum is calculated it can be compared with the experimental spectrum. If the calculated resonances are not exactly at the same frequencies as the experimental ones, the C F parameters are adjusted to fit the data and the procedure is repeated with the new parameters. • The unknown parameter A is chosen to give the best agreement between the calculated and experimental spectrum. This is done by fitting the experimental data with a sum of Lorentzians and obtaining the area of the largest transition, A , where the background subtraction error is lowest. A is extracted by comparing this area, A, with the theoretical one: A = |/xoA||/%|(iVi — Nj) where the relevant matrix element, partition function and energies were calculated by the program discussed above. Chapter 7. Results and Analysis 34 7: Results and Analysis Before performing any experiment, a prediction of the ESR spectrum was generated. The crystal structure has a nearly tetragonal symmetry (with about 1% difference between the d and b dimensions) and very high anisotropy in the c direction as compared to the a-b plane. Considering all the Stevens operators from 0 ° and higher, the following reductions should apply. The system has spin | with no orbital angular momentum and therefore the highest order Stevens operator that should be considered is 6 (since k < 2S where k is the order of the operator and S is the total angular momentum, k is also even). In addition, since the a-b plane has approximately | rotational symmetry, all of the operators that violate this symmetry (the O .) should have negligible C F coefficients. The operators that survive this selection 2 are: O®, 0 ° , 0\, 0$, 0\ and 0%. Given the dimensions of the unit cell, we expect that the diagonal terms, 0%, be larger than the off-diagonal terms. Combining this with the fact that this sum is an expansion, the dominant term in the Hamiltonian is 0 ° . The expected spectrum would contain three main transitions, since the eigenstates of the system are the Kramers doublets, | ± m >. The spacing between the E S R lines of a system with Ho oc 0% is constant since 0% = 3S — S (as was shown in figure 2 2 4.3). However, in our system the spacing between lines may not be constant due to additional C F parameters (mainly 0 ° ) . 7.1 Twinned Optimally Doped Crystal The first measurement was performed on a twinned 1.28% G d doped YBa,2CuzO§m crystal. The surface resistance results are shown in figure 7.1 along with a polynomial fit to the background. A spectrum proportional to x" w a s obtained by subtracting the background fit and dividing by frequency as shown in figure 7.2. Note that the spectrum contains six sharp E S R peaks. This is in contradiction to the predicted spectrum. Any reasonable combination of the above operators could not yield such a spectrum and the solution to the problem was eventually found elsewhere. Chapter 7. Results and Analysis 35 Figure 7.1: The surface resistance of a twinned optimally doped sample with 1.28% Gadolinium - Gdo.01yb.99-Sa2Ou3O6.g2- The graph shows the full surface resistance (points) with a polynomial fit to the background between the E S R peaks (line). Chapter 7. Results and Analysis 0 2 4 6 8 10 12 14 36 16 18 20 22 Frequency(GHz) Figure 7.2: The magnetic susceptibility x" of Gdn.oi^o.99Ba Cu 0 . up to a factor of npoX obtained from the surface resistance measurement by subtracting the background (-Rf ) and dividing by the frequency. 2 c 3 G 92 Chapter 7. Results and Analysis F Band 37 G Band Figure 7.3: T w o different C F environments. T h e G d ion may have four (G-band) or three (F-band) nearest neighbour chains. T h e stars represent a chain, open circles are holes in the C u - 0 plane and the black circle represents a missing hole (electron) 7.1.1 Oxygen Ordering Research by Rockenbauer et al. [24, 27] has shown that such a spectrum should be fitted by two sets of C F parameters. A t an oxygen doping of 6.92 the C u - 0 chains along the b axis are not completely full. A s a consequence, some G d sites have four nearest neighbour chains and some have less. A t this oxygen content, 8% of the chain oxygens are missing and approximately 32% of the G d sites have only three nearest neighbour chains (neglecting the sites with two and less chains). T h e other 68% of G d sites have four nearest neighbour chains. This results in two different C F environments for the two different G d sites as shown in figure 7.3. We will denote the spectrum of the ions with four nearest neighbour chains as the G-band and the spectrum of the ions with three nearest chains as the F-band. T h e spectrum of the optimally doped sample (shown in figure 7.2) was fitted with two sets of C F parameters. T h e first guesses were taken from the work of Rockenbauer et al. on yo.999Gdo.001 Ba2Cu^Oe,7' [24], and some adjustments were made to fit our data (see 7 the fitting procedure in chapter 6). T h e two sets are presented in table 7.1 with the C F parameters of ref. [24]. T h e transitions were identified as summarized in table 7.2 and the fitted data is shown in figure 7.4. In this case, the function x"{ ) u 1 8 a Chapter 7. Results and Analysis 38 sum of the two bands: (7.1) where r is the weight of the F band, relative to the G band. In the spectrum, r can be extracted from the relative intensities of two peaks that represent the same transition in the different bands. The Hamiltonian of spins with four nearest neighbour chains is described by the G-band C F parameters and therefore, these spins will contribute to the intensity of the Gi transitions. In the same way, the intensities of the Fi transitions are proportional to the number of spins in the F band: A F Nf - Nf r(e Nf - Nf (7.2) where Afj is the integrated intensity of the i —> j transition in the F-band and Ef is the ith energy level in the F band and (3 = 1/ksT. This ratio is directly related to the oxygen content of the crystal: 4 + 3r (7.3) 4(1+ r ) ' Using the integrated intensities of the two largest transitions in the optimally doped sample's spectrum (which correspond to the ± | —>• ± | in the G and F bands) gives an oxygen content of 6.94, in reasonably good agreement with the expected content of 6.92, from the knowledge of the annealing temperature and pressure. G-band F-band -590 (-598) -533 (-541) B'i 0(0) 0(0) s° 4 Bf Bi -3.1 (-3.1) -3.1 (-3.1) 12.7 (13.2) 12.5 (17.2) 0.167 (0.167) 0.167 (0.158) r, F : G Ratio 0.34 Table 7.1: Crystal field parameters for the two E S R bands in an optimally doped Y B C O , given in GHz. In brackets - the C F parameters as measured by Rockenbauer et al. in conventional E S R on a powder sample of 10.9990^0.001-6020^306.77 [24]. The extracted penetration depth from this fitting is 1270A and since the crystal is twinned, this is a weighted average between A and a Another contribution to this extracted penetration depth comes from A . The correction is discussed in section c 7.3. Chapter 7. Results and Analysis 5xl0 9 lxlO 1 0 1.5x 1 0 2 x l 0 1 0 1 0 5xl0 Frequency (Hz) 9 39 1 x IO 1.5 x l O 10 1 0 2x IO Frequency (Hz) 5xl0 9 lxlO 1 0 1.5 x l O 1 0 2xl0 1 0 Frequency (Hz) Figure 7.4: The fitted E S R spectrum of an optimally doped Y B C O sample with 1.28% G d . The spectrum is a superposition of two bands. Figure (a) shows The G-band fitting, figure (b) shows the F-band fitting and figure (c) shows the sum. The dots are experimental data points and the line is the theoretical spectrum generated by the Mathematical program. 10 Chapter 7. Results and Analysis Label Gi G G> G 2 2 3 Fi F 2 Frequency (GHz) 14.04 5.86 6.90 1.05 13.02 5.20 40 Transition ± 5 / 2 -»• ± 7 / 2 ± 3 / 2 -> ± 5 / 2 ± 1 / 2 ->• ± 5 / 2 ±1/2 ±3/2 ± 5 / 2 -> ± 7 / 2 ± 3 / 2 ->• ± 5 / 2 Table 7.2: ESR Transitions in a twinned optimally doped Y B C O sample. Transitions of the G-band are labeled G\ and F-band transitions are labeled Fi 7.2 Detwinned Overdoped Crystal In order to verify the two band theory another measurement was performed on a Gdo,oi sYo,g Ba Cus06^Q sample. To further simplify matters, this sample was de2 72 2 twinned. This oxygen content was chosen in order to achieve a spectrum containing predominantly one band, the G-band. A t this doping, the chains along the b axis are almost completely full, we expect a ratio of 1:0.04 between the two peaks observed earlier. Figure 7.5 shows that our expectations were verified and the spectrum contains only three main peaks with only a small contribution of the Fi transition, while the other F-band transitions are not visible. Another interesting feature of the data is the difference in spectra between the two directions. This crystal was detwinned and therefore we can use this data to separate the measurements of A and At. The spectra show a difference in the intensity as a a result of the different penetration depths in the a and b directions. In addition, the shape of the spectrum in the two directions is different. The a direction spectrum has three main peaks and in the b direction we see four main lines. The additional line (marked as G' ) is the result of an off diagonal transition, ± | —>• ± | . Such 2 a direction sensitive feature can only result from an operator that violates the | rotational symmetry in the x-y (a-b) plane. Operators of this type are the result of the crystal's orthorhombicity and the relevant operators are: 0\,0\ and 0\. The dominant contribution is of course the first orthorhombic term in the expansion: 0\ = \(S%_ + S ) = \{Sl — Sy). The coefficient B\ of this operator was adjusted 2 (using a Mathematica simulation) to simultaneously fit the data in both directions. This transition was also visible in the optimally doped crystal but since the sample was twinned, the orthorhombic nature was not resolved. The determination of the C F parameter B\ can therefore only be done on detwinned samples. The crystal field parameters are given in table 7.3 and the identified transitions are given in table 7.4. Chapter 7. Results and Analysis 41 1000 10 Frequency(GHz) 15 1000 10 15 Frequency(GHz) Figure 7.5: Surface resistance of an over doped Y B C O sample with 1.28% Gd. Figure (a) shows the spectrum of the a direction and figure (b) shows the b direction. In both figures the measured surface resistance is shown (dots) along with a polynomial fit to the background (line) Chapter 7. Results and Analysis 42 The data with the theoretical fit for the d and b directions are shown in figure 7.6. G-band F-band B\ 26.7 (0) 0(0) -583 (-598) -530 (-541) r, F : G Ratio B\ B\ -3.1 (-3.1) 12.7 (13.2) -3.1 (-3.1). 12.8 (17.2) 0.04 0.167 (0.167) 0.159 (0.158) Table 7.3: Crystal field parameters for the E S R spectrum in an overdoped Y B C O crystal, given in GHz. In brackets - the C F parameters as measured by Rockenbauer et al. in conventional E S R on a powder sample of yo.999Gdo.ooijBa2Cu3O6.77 [24]. Label G G G' G Fi F 1 2 2 z 2 Frequency (GHz) 14.04 5.86 6.90 1.05 13.02 5.20 Transition ± 5 / 2 -» ± 7 / 2 ± 3 / 2 -> ± 5 / 2 ± 1 / 2 -»• ± 5 / 2 ± 1 / 2 -> ± 3 / 2 ±5/2 ±7/2 ± 3 / 2 -> ± 5 / 2 Table 7.4: ESR Transitions in a detwinned overdoped Y B C O sample. Both the a and b spectra show the same transitions except for the G transition that is only visible in the b direction 2 7.3 Corrections for c Axis Currents Corrections to the above figures should be made because of c-axis contributions to the effective penetration depth. In our experiment the magnetic field is applied along the b (d) direction and the currents flow mainly in the a-b plane. However, a small portion of the current flows in the b-c plane, along the c direction. As a result the measured penetration depth, A / / , is a mixture of the in plane penetration depth, A e a (A/,), and A . A method used for separating the two contributions was employed at the c U B C superconductivity lab [26, 32] in overdoped Y B C O crystals (x = 0.95 - 0.99). This technique involves a cavity perturbation measurement of R . The sample is then s cleaved into 2 or 3 long platelets along the b (a) direction and R is measured again. s The increase in the observed surface resistance can be attributed to the additional b-c planes in the platelets and therefore to the c axis transport. In the case of G d doped crystals, this measurement has not been done yet and the corrections given here are based on previous measurements of the c penetration depth. Chapter 7. Results and Analysis 10 43 15 Frequency (GHz) Figure 7.6: The fitted E S R spectrum of an over doped Y B C O sample with 1.28% Gd. The spectrum has essentially only one E S R band, the G-band. Figure (a) shows the spectrum of the a direction, where the currents are running in the a direction and the magnetic field is applied along the b direction. Figure (b) shows the b direction, where the middle main line is split into two peaks by the 0\ crystal field operator. The dots are experimental data points and the line is the theoretical spectrum generated by our Mathematica simulation. Chapter 7. Results and Analysis 44 We used the value of A = 11,000±400A that was measured in YBazCuzOs.w at 1.2K c [32]. The first, twinned sample had dimensions of x = 820//m and c = 19/im where x is the length of the sample along the measured direction. The second, detwinned sample has dimensions of a = 770/xra, b = 910/xm and c = 23pm. A simple, straightforward correction assumes that the G d ions that are exposed to the field near the b-c plane are in a volume of 2 x a x b x A : c +cA aX a A // = c (7.4) e A more accurate calculation([33]) of a field profile in a sample of rectangular cross section uses the anisotropic screening equation: ,8 H A; dz dH A; = H. dx 2 2 2 a (7.5) 2 The solution to this equation is: A // e ac 4 a + c it 2 T 1 J tanh(a ) - odd n>0 with n \ n Ctn a n P n tanh{(3 ) 2A 1 + a 1 + ~ 2X r Table 7.5 summarize the results for X i p ane (7.6) Pn V (7.7) a / rvnXaf V 2 (7.8) c of the twinned crystal and A and A^ for a the detwinned crystal and their corrections in the two methods. Xpiane, twinned sample A , detwinned sample Xb, detwinned sample a extracted A(A) 1270 ± 130 1684 ± 175 1347 ± 140 simple c correction 1019 ± 140 1355 ± 187 1069 ± 150 full calculation 1025 ± 140 1410 ± 187 1100 ± 150 Table 7.5: The ESR-extracted penetration depth A of an optimally doped twinned sample and an overdoped detwinned sample, at T = IK. Both samples have 1.28% Gd concentration. Chapter 7. Results and Analysis 7.4 45 Error analysis Errors in the extracted A come from both experimental errors and the function fitting procedure. The experimental errors in the bolometric apparatus are estimated as 510% in the raw R data [19]. Another source of error is the Gadolinium concentration s in the crystal that is estimated to be 1.28 ± 0.02%. In the fitting procedure there are two sources of error, the first is the background subtraction in which we fit a polynomial function to the data between E S R peaks. This fitting is quite accurate away from the peaks and only affects the tail of the Lorentzians. After the background Rg is subtracted, the data is divided by the frequency to give a quantity proportional C to x" • F ° the small peak the background fitting inaccuracy is enhanced by a large r factor of order ^ since the frequency is low and the actual line intensity may be changed by up to 50%. However, for the largest peak the background subtraction error is small compared to the intensity of the line and is divided by a large factor since the frequency is high. The best fit parameter A is chosen to fit the largest transition and therefore we neglected the background subtraction error. The second fitting error is the actual curve fitting of the theoretical data (a sum of Lorentzians) to the experimental data. A n estimation of this error is given by Origin (the data analysis program) as the errors of the fitted Lorentzian area. This error is estimated to be less then 2%. The errors mentioned above are uncorrelated and therefore add in the following way: where Ned is the gadolinium concentration, R is the surface resistance and A is the s intensity of the largest ESR transition that is used to fit A. This leads to the following relative error in A L (before c axis correction): (^^j = v (10%) + (2%) + (2%) = 10.4% / 2 2 (7.10) 2 Another error comes from the A value that is used for the c axis correction. We C can calculate the error of the corrected A A (or A^) using the simple formula for the c axis correction: A„ = A E / / - -A (7.11) C where the error in A / / was calculated in the previous paragraph and A = 11, 000 ± E C Chapter 7. Results and Analysis 46 400A. The total error in A is then given by: S\ = S\ eff + -5\ . a c The errors for each value of A are given in table 7.5. (7.12) Chapter 8. Conclusions 47 8: Conclusions The m a i n goal of this work was to develop a reliable method for measuring the absolute value of the low temperature limit of the magnetic penetration depth, \L(T —>• 0). T h i s goal was convincingly achieved through a surprising and unconven- tional combination of two experimental methods, E S R and broad band microwave surface resistance measurements. T h e new technique was applied to Y B C O crystals and the results that were obtained are in reasonable agreement with uSR [8, 9] and infrared [10] measurements. However, there are still some open issues regarding the reliability and the errors of the figures given. Most concerning is the difference between our measured A , X of b a the overdoped detwinned sample and X i of the optimally doped twinned crystal. p ane We expected A ; p a n e to be intermediate to the values of A and A& (when measured in a the same sample) since it is a weighted average of the two directions. In practice, Xpiane was found to be smaller then both A and A&. T h i s may indicate a change in a penetration depth due to the different superfluid density but may also be the result of a difference in G d concentration between the two samples. T h e later possibility is currently being checked out. Several directions of study are possible for further work. Performing a series of measurements on Y B C O with different oxygen contents will provide information on the dependence of the superfluid density and the crystal field parameters on doping. The results can then be combined with T C measurements of the same samples and presented as a Uemura plot [3] (when we plot T vs. 1 / A ) to verify the dependence of 2 C T C on superfluid density. Another product of this series of measurements would be the C F parameters for each doping. This would complement the work of Rockenbauer et al. [24] in single crystals (rather than aligned powder) and extend it to higher doping levels. Another useful measurement that can be done with the existing samples is the caxis penetration depth A . T h i s can be done by cleaving the sample as described c in reference [26]. T h i s measurement will provide, in addition to the value of A , the c appropriate c-axis correction to our X /ba However, by far the most significant result of the success of this method is that it allows us to measure A absolutely, on crystals that are essentially identical to the Chapter 8. Conclusions 48 crystals used in the lab's surface resistance measurements. In turns this allows us to extract the microwave conductivity o~\ from the surface resistance values. In the past all of our conductivity values have had relatively large uncertainties because of uncertainties in A . It may also be possible to use the same measurement principle in other superconductors. The feasibility of the technique depends on finding the right atomic substitution and knowing the fine structure of the ion at the relevant crystal site, at least to first order. Knowing the C F parameters, one can then calculate the frequencies of the various transitions and design an experiment to cover the relevant range. As described in chapter 5, the experiments were done in a bolometric apparatus [19]. Further development of the bolometric experiment is under way. The frequency range will be extended to higher frequencies (up to 50GHz) and to lower temperatures (down to lOOmK). With these experimental improvements we will be able to measure A L closer to the limit of zero temperature, and the broader band of frequencies may become useful in measurements of other materials. References 49 References [1] J . Bardeen, L . N . Cooper, and J . R. Schriffer. Theory of superconductivity. Physical Review, 108:1175-1204, December 1957. [2] W . N . Hardy, D. A . Bonn, D.C. Morgan, Ruixing Liang, and Kuan Zhang. Precision measurement of the temperature dependence of A in YBa2Cus06.95'. Strong evidence for nodes in the gap function. Physical Review Letters, 70(25) :39994002, June 1993. [3] Y . J . Uemura, G . M . Luke, B . J . Sterlieb, J . H . Brewer, J . F. Carolan, W . N . Hardy, R. Kadono, J. R. Kempton, R. F. Kiefl, S. R. Kreitzman, P. Mulhern, T. M . Riseman, D. L i . Wlilams, B . X . Yang, S. Uchida, H . Takagi, J . Gopalakrishnan, A . W . Sleight, M . A . Subramanian, C. L . Chien, M . Z. Ciepalk, Gang Xiao, V . Y . Lee, B . W . Statt, C. E . Stronach, W . J. Kossler, and X . H . Y u . Universal correlation between t and n /m* (carrier density over effective mass) c in high-7j c s cuprate superconductors. Physical Review Letters, 62(19):2317-2320, May 1989. [4] A . Hosseini, R. Harris, Saeid Kamal, P. Dosanjh, J . Preston, Ruixing Liang, W . N . Hardy, and D . A . Bonn. Microwave spectroscopy of thermally excited quasiparticles in YBa Cu306.99. Physical Review B, 60(2): 1349-1359, July 1999. 2 [5] S. Kamal, Ruixing Liang, A . Hosseini, D.A. Bonn, and W . N.Hardy. Magnetic penetration depth and surface resistance in ultrahigh-purity YBa2Cu30y_5 crystals. Physical Review B, 58(14) :R8933-R8936, October 1998. [6] S. Kamal, D . A . Bonn, Nigel Goldenfeld, P.J. Hirschfeld, Ruixing Liang, and W . N.Hardy. Penetration depth measurements of 3d xy critical behavior in YBa2Cu306.95 crystals. Physical Review Letters, 73(13):1845-1848, September 1994. [7] D. A . Bonn and W . N . Hardy. Microwave surface impedance of high temperature superconductors. In Donald M . Ginsberg, editor, Physical Properties of High References 50 Temperature Superconductors, volume V , chapter 2, pages 7-98. World Scientific, 1995. [8] Jeff E . Sonier, Jess H . Brewer, and Robert F. Kiefl. psi studies of the vortex state in type-ii superconductors. Reviews of modern physics, 72(3):769-807, July 2000. [9] J . L . Tallon, C. Bernhard, U . Binninger, A . Hofer, G . V . M . Williams, E . J. Ansaldo, J . I. Budnick, and Ch. Niedermayer. In-plane anisoropy of penetration depth due to superconductivity of cu-o chains in Y B a 2 C u 0 _ 5 , Y Ba Cu 0i5_,5, 3 7 2 4 7 and Y B a C u 0 . Physical review letters, 74(6): 1008-1011, February 1995. 2 4 8 [10] D . N . Basov, R.Liang, D.A. Bonn, W . N . Hardy, B . Dabrowski, M . Quijada, D . B Tanner, J.P. Rice, D . M . Ginsberg, and T. Timusk. In-plane anisoropy of the penetration depth in Y B a C u 3 0 _ and Y B a C u 0 8 superconductors. Physical 2 7 x 2 4 review letters, 74(4)-.598-601, January 1995. [11] M . Tinkham. Superconductivity. Gordon and Breach, science publishers, 1965. [12] D . A . Bonn and W . N . Hardy. Microwave properties of high temperature superconductors. Physics in Canada, 56(5):243-246, September/October 2000. [13] John David Jackson. Classical Electrodynamics. John Wiely & Sons, Inc., 3rd edition, 1998. [14] P. G . De Gennes. Superconductivity of Metals and Alloys, advanced book program. Persus books, 1999. [15] J . S. Griffith. The Theory of Transition Metal Ions. The Camridge University Press, 1961. [16] K . W . H . Stevens. Magnetic Ions in Crystals. Princeton university press, 1997. [17] A . Abargam and B . Bleaney. Electron Paramagnetic Resonance of Transition Ions. Clarendon press - Oxford, 1970. [18] M . Divis and V . Nekvasil. Semi-empirical and ab-initio calculations of the crystal field interactions in rare earth cuprates. Journal of Alloys and Compaunds, 1, 2001. References 51 [19] Patrick J . Turner. A bolometric technique for broadband measurement of suface resistance and application to single crystal Y B a C u 2 3 0 .99. 6 Master's thesis, The University of British Columbia, October 1999. [20] D . Shaltiel, S.E. Barnes, H . Bill, M . Francois, H . Hagemann, J . Jegondaz, D. Lovy, P. Monod, M . Peter, A . Revcolevschi, W . Sadowaski, and E . Walker. Unusual behavior of gd esr in single crystals of G d Y i _ B a C u 3 0 6 + y y 2 x with x = 0.1 — 0.8 and y = 0.03 — 0.06: Evidence for a magnetic interaction in the superconductors. Physica C, 161:13-20, August 1989. [21] A Janossy, A . Rockenbauer, and S. Pekker. G d 3 + electron paramagnetic reso- nance in the quazi-two-domentional antiferromagnet YBa Cu306+ : G d . Phys2 y ica C, 167:301-306, 1990. [22] A Janossy, A . Rockenbauer, S. Pekker, G . Oszlanyi, G . Faigel, and L . Korecz. Epr in gd-doped Y B a C u 0 6 + : local method to determine oxygen ordering. 2 3 y Physica C, 171:457-464, 1990. [23] A Janossy, L . - C . Brunei, and J. R. Cooper. G d esr determination of the spin 3 + susceptibility in G d : YBa Cu306+ high-temperature superconductors. Physical 2 y review B, 54(14):10 186 - 10 191, 1996. [24] A . Rockenbauer, A . Janossy, L . Krecz, and S. Pekker. Electron paramagnetic resonance spectroscopy of gd doped Y B a C u 0 6 + in a magnetic field. Journal 2 3 x of Magnetic Resonance, 97:540-552, August 1992. [25] S. Sridhar and W . L . Kennedy. Novel technique to measure the microwave response of high t between 4.2 and 200k. c Review of scientific instruments, 59(4):531-536, April 1988. [26] A . Hosseini, Saeid Kamal, D . A . Bonn, Ruixing Liang, and W . N . Hardy, c-axis electrodynamics of YBa Cu30,5_7. Physical review letters, 81:1298-1301, March 2 1998. [27] S. Pekker, A . Janossy, and A . Rockenbauer. Ordering of oxygen into chains and distribution of hole density in YBa Cu306+ - Physica C, 181:11-17, 1991. 2 x [28] Ruixing Liang, D . A Bonn, and W . N . Hardy. Perfecting the growth of ybco crystals. Physics in Canada, 56(5):243-246, September/October 2000. References 52 [29] Ruixing Liang, D.A. Bonn, and Walter N . Hardy. Preparation and x-ray characterization of highly ordered ortho-ii phase Y B a C u 3 0 . 5 single crystals. Physica 2 6 C, 336:57-62, 2000. [30] P. Schleger, W . N . Hardy, and B . X . Yang. Thermodynamics of oxygen in YxBaaCusCv Physica C, 176:261-273, 1991. [31] Charles Kittel. Introduction to Solid State Physics. John Wiely & Sons, Inc., 6th edition, 1986. [32] C . C . Homes, S. Kamal, D . A . Bonn, R. Liang, W . N . Hardy, and B.P.Clayman. Determination of the condensate from optical techniques in unconventional superconductors. Physica C, 296:230-240, 1998. [33] David Broun. The Microwave Electrodynamics of Unconventional Superconductors. P h D thesis, Pembroke College, Cambridge, January 2000. Appendix A. Cavity Perturbation Technique 53 A: Cavity Perturbation Technique Cavity perturbation is a particularly useful tool for measuring the microwave electrodynamics of small single crystals of high T superconductors. This appendix explains c the basic idea of cavity perturbation technique. In general a resonance cavity has a complex eigen frequency that is composed of the frequency u (the real part) and the inverse quality factor Q (the imaginary part). * = +55) (A.l) In an empty cavity the resonance frequency is a function of the cavity's geometry and the quality factor is a function of the losses in the cavity walls. When a sample (metal or superconductor) is inserted into the cavity the resonance frequency changes. The change is due to the new boundary conditions for the electromagnetic fields and resulting new field configuration as well as the loss on the sample's surface (the sample's surface resistance). In order to measure the temperature dependence of the sample's surface impedance the following steps are taken. The losses in the cavity walls are kept constant throughout the experiment, so that any change in Q is caused by a change in the sample's surface resistance. This is done by keeping the cavity at a constant base temperature and heating only the sample above this temperature. The field configuration in the cavity should be constant so that the changes in the resonance frequency would only be due to the change in the sample's penetration depth. This is done by keeping the sample at one place and working at a resonance frequency that is well separated from other resonances. When the above steps are taken, the changes in u are only a function of the changes in the sample's surface impedance. The extraction of the sample's surface impedance is based on the relation: Auoci j J AZ J . 2 S (A.2) sample The integrated surface current is assumed to be unchanged and therefore the change Appendix A. Cavity Perturbation Technique 54 in surface impedance can be resolved: Aw = Au + iA^- = %KAZ = %KAR - KAX . a S (A.3) S The proportionality constant K is a function of the cavity and sample geometry and can be calculated for simple geometries or measured using a reference sample with known properties and similar dimensions. Achieving the absolute surface impedance of the sample in this way is not trivial. The relevant origin of the complex frequency, (Do, is of the cavity with a sample of no surface resistance and zero penetration depth. This limit is of course not accessible by experiment. The measured Q of the cavity without the sample can be taken as a good zero point for the surface resistance R since then the sample's contribution s to the loss is zero. For the surface reactance, this does not work since the origin of X is a sample with zero penetration depth ( acavity with no sample is the opposite s situation). In this approach (taken by our group) one can only measure AX . S Another approach is to use the fact that at the normal state R = X s s [?]. The frequency origin is set in the normal state and the ^ origin is set well in the superconducting state where the loss is orders of magnitude smaller. This method yields an absolute measurement of Z but has large uncertainties. s Appendix B. Surface Current in a Uniform magnetic Field 55 B: Surface Current in a Uniform magnetic Field This appendix explains equation 5.2 of chapter 5: psample preference _ psample s R ference /g L \ r e j\ ' ' where the reference and sample experience the same magnetic field H. As discussed in chapter 3, the classical skin effect leads to an exponential decay of A C magnetic fields into the sample. In a normal metal this decay is proportional to i{i+i)r/s e to e~ l r x w h e r e j the skin depth and in a superconductor the decay is proportional s where A is the magnetic penetration depth. In both cases the magnetic field deep inside the sample is essentially zero and the currents flow on the surface within the small penetration depth or skin depth. This results in a surface current, J , (the s current density integrated over the thickness A or 8) that is equal to the applied magnetic field. To obtain this consider the geometry described in figure B . l where the magnetic field, H, is applied parallel to the crystal's a axis and the currents are flowing in the a-b and a-c planes perpendicular to the d direction. Now let us integrate Ampere's law, V x H = J, over the rectangular surface S that is parallel the a-c plane and has an edge I (see figure B . l ) . The upper edge of S lays on the upper a-b plane, where the supercurrent is flowing, and the lower edge is well inside the sample where both the magnetic field and the current are zero. The integration gives: j"(V x H)-dS = J J-dS'= J -l s (B.2) where the last equality is consistent with our definition of the surface current per unit width (in this case the current is flowing mainly in the b direction and the width is along the d direction). Using Stokes theorem we can simplify the left hand side of the equation: f(VxH)-dS=(f J H-dl (B.3) JdS where dl runs over the contour of the surface S. Since H is applied in the d direction and is zero in the middle of the sample, the only contribution to this integral is from Appendix B. Surface Current in a Uniform magnetic Field Figure B . l : A typical high T 56 sample in a magnetic field. T h e magnetic rf field i is c applied along the a direction the current is flowing mainly in the a - b plane along the b direction. T h e field and current deep inside the sample are zero. the part of the contour along the sample's surface (in the d direction) and therefore: H a •I — J s • I =$> H a = J (B.4) s W i t h this result we can now calculate the ratio of absorbed power in the sample and reference. T h e power absorption of a sample per unit area in an rf magnetic field is given by •Psample — ^R J s where R s g — — RH 2 S (B-5) is the surface resistance and J is the current density. Since the two sampl s experience the same magnetic field, result B . l follows. es
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Zero field ESR of magnetic impurities in superconductors...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Zero field ESR of magnetic impurities in superconductors : a novel technique for measuring magnetic penetration… Pereg-Barnea, Tamar 2001
pdf
Page Metadata
Item Metadata
Title | Zero field ESR of magnetic impurities in superconductors : a novel technique for measuring magnetic penetration depth and application to high Tc YBa2Cu3O7-[delta] |
Creator |
Pereg-Barnea, Tamar |
Date Issued | 2001 |
Description | The Meissner effect in superconductors is characterized by an exponential decay of externally applied magnetic fields into the sample with a length scale A, that is called the magnetic penetration depth. λ is directly related to the superfluid density of the superconductor and is therefore one of the most important probes of the intrinsic properties of superconductors. This work describes the development of a new technique for measuring the penetration depth at microwave frequencies and low temperatures. The surprising combination of two seemingly unrelated experiments has led to a direct method of measuring the magnetic penetration depth in the London limit of high T[sub c] superconductors. Samples of high T[sub c] YBCO with uniform doping of magnetic ions (Gd[sub x]Y[sub 1-x]Ba₂Cu₃0 ₇-[sub δ] where x ~ 1.28%) were prepared. The magnetic ions embedded in the superconductor serve as magnetic field probes through a zero field Electron Spin Resonance (ESR) experiment. Measurement of the imaginary part of the magnetic susceptibility, x"(ω), through the sample's microwave surface resistance, reveals relatively sharp ESR transitions with integrated intensity that is proportional to the number of ions that are exposed to the field. With a uniform ion distribution, the number of ions that are exposed to the field is proportional to the penetration depth λ. The spectrum of x”(ω) is obtained from surface resistance measurement in a bolometric apparatus. Our new method was successfully applied to overdoped and optimally doped YBCO samples and absolute values of the a-b plane penetration depth at T = 1.2 are given. |
Extent | 3972010 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | eng |
Date Available | 2009-08-06 |
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.0085714 |
URI | http://hdl.handle.net/2429/11792 |
Degree |
Master of Science - MSc |
Program |
Astronomy |
Affiliation |
Science, Faculty of Physics and Astronomy, Department of |
Degree Grantor | University of British Columbia |
Graduation Date | 2001-11 |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-ubc_2001-0492.pdf [ 3.79MB ]
- Metadata
- JSON: 831-1.0085714.json
- JSON-LD: 831-1.0085714-ld.json
- RDF/XML (Pretty): 831-1.0085714-rdf.xml
- RDF/JSON: 831-1.0085714-rdf.json
- Turtle: 831-1.0085714-turtle.txt
- N-Triples: 831-1.0085714-rdf-ntriples.txt
- Original Record: 831-1.0085714-source.json
- Full Text
- 831-1.0085714-fulltext.txt
- Citation
- 831-1.0085714.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-0085714/manifest