ES flu 051 any SUFERCQNBE he Eh {. we ’ (‘ E i Exes \- ”dam-lg we h‘ 5 till . u 32%;.” ‘ Y “‘1 :a 3: I‘: i‘. an}. 1... L LIBRARY Michigan State University __.- _ This is to certify that the thesis entitled Theory of Transition Element Superconductors presented by Paul B. Williams has been accepted towards fulfillment of the requirements for ? Ph.D. Physics degree in M/% Mdor professor DateWO 0-169 ABSTRACT THEORY OF TRANSITION ELEMENT SUPERCONDUCTORS by Paul B. Williams A theoretical expression for electron-phonon coupling is derived and then used in an attempt to explain differ- ences in superconducting transition temperatures for groups VB and VIB transition metals. The bare ionic potential, based upon the Herman- Skillman atomic structure calculations, is_used to determine unscreened coupling. Screening is then calculated using the Bardeen self-consistent field method modified for exchange. The net result is an expression for the screened electron- phonon coupling, gK K" exhibiting a positive dependence upon the electronic density of states at the Fermi level. Such a dependence is due entirely to exchange and is thought to stem from the preference of conduction electrons to group between nearest neighbors in some resonant structure remnant of atomic "d"—electrons. A suitable computer program is developed to determine the product of coupling times phonon density of states as a function of phonon frequency, a2(cn)F(¢n). The method first uses interatomic force constants to determine the phonon Paul B. Williams density of states, F(m0, via the Born-von Karman method. The expression for screened electron-phonon coupling is then incorporated to include umklapp as well as normal scattering processes. Finally, within the currently accepted theory of superconductivity, a2(m)F(cn) is used to determine an expression for the superconducting transition temperature of tantalum and tungsten, THEORY OF TRANSITION ELEMENT SUPERC ONDUCTORS /U9 Paul B? Williams A THESIS Submitted to Michigan State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Department of Physics 1970 ACKNOWLEDGMENTS I wish to thank Dr. Michael Harrison for his excellent guidance and encouragement throughout the course of this research. Special thanks to Dr. ‘William.Hartmann for providing a copy of the Gilat and Raubenheimer program, G(NU), and for some very helpful advice in the area of accoustical properties of solids. Also, I wish to thank Jean Strachan for a truly outstanding Job of typing. 11 TABLE OF CONTENTS ACKNOWLEDGMENTS . . . . . . . . . . . . . . . . . . 11 LIST OF TABLES . . . . . . . . . . . . . . . . . . IV LIST OF FIGURES . . . . . . . . . . . . . . . . . . V Chapter I INTRODUCTION........_......l II THE'I‘HEORY.’...............14 III THEGAPEQUATIONS............ 29 Iv THE ELECTRON-PHONON COUPLING CONSTANT . . 46 v EVALUATION OF MATRIX ELEMENTS . . . . . . 66 VI THEORETICAL INVESTIGATION OF SUPERCONDUCTING TAN'I‘ALUM AND TUNGSTEN . . 7% .VII CONCLUSION...............91 Appendices A EVALUATIONOFInlm...........96 B EVALUATIONOFINTEGRALS.........100 c chPtJTERPRoanN.............105 LIST OF WRENCESO O O O O O O O 0 O O O O O O O O 139 111 Table 1. 2. LIST OF TABLES Page IsotOpe effect (experimental for some simple and transition metals . . . . . . . 5 Summary of physical parameters influencing superconductivity . . . . . . lO Character table of the cubic group O . . . 20 List of potential coefficients and force constants for tantalum.. . . . . . . 84 List of potential coefficients and force constants for tungsten . . . . . . . 85 A vs Kf for tantalum.and tungsten 92 Figure 1. 2. 5. 6. 7. 8. 10. 11. LIST OF FIGURES Qualitative behavior of superconducting transition temperatures through the ‘ periodic system . . . . . . . . . . . . Areas of maximum electron density within the B~C.C. primitive cell for electrons having group symmetry. . . . . . . . . . . . . Important diagrams contributing to self enem O O O O O O O O O O O O O O O O O Contours in com lex plane showing some .polesatoddffl .......... Theoretical density of states for tungsten from band structure data . . . Density of phonon states for tantalum . Density of phonon states for tungsten . Coupling times density of states for tarltalumoe_eoeee'oeeoeoeee Coupling times density of states for" tlmsa ten 0 C O O O O C C C O O O O I O O A vs Kf for tantalum. . . . . . . . . . A vs Kf for tungsten. . . . . . . . . . Page 12 23 35 38 45 87 88 89 2‘88 CHAPTER I INTRODUCTION The first successful microscopic theory of supercon- ductivity was presented by Bardeen, Cooper, and Schrieffer (3.0.8.) in 19571 and then later refined into a more ele- 2 and Bogoliubov’. In gant and accepted form by valatin simplest terms, the theory assumes attractive Cooperll elec- trons forming pairs of opposite momentum and spin which re- side in a condensed ground state '15; :> . A reduced hamiltonian (fired ) which considers only the elec- tron pairs is used in conjunction with a parameterized trial ground state wavefunction subject to the constraint that the expectation value of the total number of particles r is e b . Let d som fixed num er No Ck? an Ckl be creation and annihilation operators for electrons in state k(spin 1/2) and state -k(spin -l/2) respectively. Then Hr“? El: 6 le-(T jI-4C—k J'CH ”1+2 was C—kl Cm Ir, IE kZSCkTC -|= N with ’fl’l + N. CT |< |W>:fl k V1 + XL?— The parameter X“ in the ground state wavefunction is then determined by normal variational methods for a minimum ground state free energy (see Schrieffer5). The key results are that the pairing phenomenonis effective within a small range of energies of approximately 2"“ mo ( COD =- Debye frequency) about the Fermi energy and that quasiparticle excitations from the ground state always have an energy greater than or equal to an energy gap A. The B.C.S. result for A is, in the weak coupling limit (small . -1 0) V) AEZ’fiwDe /N( \/ 2 Assskfig where N(o) is the density of electronic states at the Fermi energy and'v is an effective pairing potential acting only within the range of energies about the Fermi level where pairing is approximately assumed to exist. Although the B.C.S. theory was highly successful in predicting qualita- tive features of superconductors, it soon became apparent that a more powerful theory was due. One had to consider the retarded nature of the electron-phonon interaction and the screened coulomb repulsion between electrons before serious attempts could be made to calculate quantitative properties of superconductors from first principles. A particularly useful feature of the B.C.S. theory is the isotope effect. If one considers the isotopes of any particular element, each arranged in its own lattice, the normal frequencies of vibration will be (.0200 k/|Vl in the harmonic approximation where M is the mass of the isotope in the lattice in question and k is an effective spring constant common to all the isotopes of that partic- ular element. Since the transition temperature To satis- fies Eq. 4, then the ratio of the transition temperatures between any two isotopes of a given element is %=\/m— . This important result is known as the isotope effect. Al- though it does not include the effects of mutual coulomb repulsion between pairs, it is supported quite well by the simple metals including such strong coupled superconductors (V large) as lead and mercury. Hewever, transition.metal superconductors are a different story altogether. They seem to exhibit little or no isotope effect and agreement exists in too few cases to produce any firm conclusions. The results of measurements of the superconducting tran- sition temperatures of various isotopes is shown in table 1 where TOO M-F 6. It would seem that for the transitidfi metals, the B.C.S. theory fails because either the electron pairing potential is not phonon mediated but possibly due to some overscreening by nonsuperconducting electrons in different bands as suggested by several auth- ors7 or that some electronic mechanism, intermediate to the superconducting electrons and the ions, serves to alter the bare electron-phonon coupling constant in ways more complex than previously described. Not long after the advent of the B.C.S. theory, theo- reticians attempted to apply Green function methods to the superconducting many body problem. Their goal was to dee velop a theory that would predict the superconducting pro- perties of metallic superconductors from the most basic properties of the metals themselves. Ordinary Green Table 1. Isotope effect (experimental) for some simple and transition metals.* Simple Transition Metal 8 Metal 8 Pb 0.48 Os 0.21 H8 0.50 Ru 0 Sn 0.51 Mo 0.3} T1 0.49 Zn 0.50 Cd O.#O * The quantity 8 is defined with TC~M‘3. functions soon met with failure due to infinities resulting from summing over certain Feynman. diagrams of a given order. The culprits are those diagrams where two elec- trons of opposite momenta and spin repeatedly scatter against each other. These difficulties were soon sur- mounted, however, with the introduction of "anomalous prop- agators" (see Gorkov8) which destroy and create Cooper4 pairs in the ground state. The entire Green function prob- lem was reformulated into a matrix form (Nambug) and was successfully applied by Eliashberg10 to the problem of strong coupled superconductors. Strong coupling was particularly troublesome because a large electron-phonon coupling constant would tend to produce such energetic quasiparticle excitations that the resulting lifetimes would be considerably shortened (be; cause of increased available phase space for their decay) to the point where the width (due to uncertainty considere ations) of the quasiparticle levels is comparable to the energy of the excitations themselves. The problem was acute until an important result by Migdal11 10 was applied by Eliashberg yielding a set of coupled integral equa- tions which form the basis of the theory as it stands today. Important coulomb effects were soon incorporated12 giving a set of equations (known as the Eliashberg gap equations) which should enable one to calculate the super- conducting properties of metals and semiconductors to a relatively high degree of accuracy. These equations form a vital part of this research and, therefore, will be de- veloped in detail later. Use of the Eliashberg equations to find the supercon- ducting transition temperature Tc was successfully applied to aluminum, sodium, and potassium by Carbotte and Dynes13. Essentially, they began by using the Heine-Abarenkov form 1“ from which they of the electron-ion pseudopotential could obtain the electron-phonon coupling constant. The necessary phonon characteristics were represented through the use of a Born-yon Karman force constant fit to phonon dispersion data in the high crystalline symmetry directions where the number of nearest neighbors considered depended upon the complexity of the phonon dispersion curves forthe superconductor in question. The Eliashberg equations were then solved self consistently on a computer yielding re- sults that were quite successful. The transition temper- ature for both nonsuperconducting sodium and potassium was determined to be less thanAESOK. For aluminum, the pre- dicted transition temperature differed by less than 10% from the experimental value. In addition, the first strong coupled superconductor to be investigated using the Eliashberg equations was lead which has the pecularities as described earlier. Experimental results for lead strong- ly support the Eliashberg form as opposed to previous mod- els. Finally, little or no progress has been made in deter? 'mining the superconducting properties of transition metals from first principles. Their electronic band structure is of such complexity that it is difficult, if not impossible, to make use of pseudopotential methods to determine values of electron-phonon coupling. If the Eliashberg equations are considered a truly accurate set of solutions to the problem or superconduc- tivity and if superconducting properties can be accurately measured, then it would be logical to ask if these equa- tions could be applied in reverse to determine some of the more elusive properties of solids from their superconducting properties. Such a problem was undertaken by McMillan15 who applied rit to a variety of materials including transi- tion metal superconductors and their alloys. The primary objectives werean empirical electron-phonon coupling con; stant >\ and a coulomb coupling constant ,5”: as 16 defined by Mbrel and Anderson . From these empirical con- 2 stants the quantity (the average squared elec- 17 was found for several tronic transition matrix element transition metal superconductors. All of this was accom- plished as follows. The average electron-phonon coupling 2 0C (w) as a function of frequency and the phonon density of states F(w) are determined as the product OC 2(a)) Fan) directly through the insertion of experimentally determined superconducting properties into the Eliashberg equations. Then the quantity (:(2; is: <6»: N(o) where )( ::MZU/‘-£L£¥2,ififouz)F:2CL» «02>: wa oc2(w)l——(oe) co fiwflz(w)F-(w) The results are summarized in table 2. In this scheme McMillan has demonstrated that, for group'V and VI B.C.C. transition metals (except magnetic chromium), the factor 2 N(O) remains essentially constant 2 and that To depends upon M <60 > (the stiff- ness of the lattice). These results are uncommon and, as he pointed out, no satisfactory theoretical explanation has been given for such phenomena. The present theory of superconductivity seems to be satisfactory with the possible exception of the transition 10 :m.m w.5 m5 0n.0 00.8 mod mH.5 pm $5.H :.m OHH HN.0 mm.H NHH 0:.n CH 0.N 5.m 0mm mON.0 mn.H mm: ©H.H H4 n.© m.m: 0mm ®:H.0 0m.0 0mm NH0.0 3 w.® 0.:N OHM m5m.0 mw.H 00: mm.0 0: H.© 0.5 05H 55.0 0.0 mmm w#.# we N.5 0.5 0mm Hm.0 w.5 55m mm.m DZ 0.: m.n 0mm HM.H m.m man on.m > -Amm\>ov ANM\N>TV. Aeoum>o\moumumVAmxoofios\bsv Axmv Ax v 53302 58 «>53. A32 > e we deem: mfimufi>fiuosccooammsm wcfiocosawcfi maouoemamm Hmofimhnn mo assessm .m oHme 11 metals. These metals can be among the best superconductors and, when properly alloyed, the transition temperature and critical field can be greatly increased (i.e. Nb3 8n and V381) which seems to indicate that physical features which give rise to superconductivity are very sensitive to the electronic structure in the case of Nb and V. If one plots: the transition temperature Tc versus the number of valence electrons per atom for transition metals and compares it with a similar plot for simple metals, as was done by, flatthiasle, (Figure 1), there appears to be little agree- ment as to how'l‘c varies with the number of valence elec- trons per atom. This certainly suggests that, as in many other phenomena, transition metal superconductors must re- quire a modified or perhaps totally different treatment' _than was used for the case of simple metals. In fact, if one compares the transition temperatures of the group V transition metals with those of group VI transition metals, it becomes apparent that, with the addition of one more non- core electron per atom, the excellent superconducting prep- erties of group v are lost in group VI.’ Assuming that the theory of superconductivity is correct, it is then logical to try to find a means to determine the degree of electron- phonon coupling and the best value for the density of states (including only those electrons that can take part in pair- s-dl IMPLE ETAL S M _+ ANSITION METALS PR S'IVLBW NON j‘l Q \Z‘ .3 4 6 8 N 0 ‘0 ¢ ATOM inq tr sition VALENCE ELE CTRONS PER 13 ing) for these two groups. It certainly would not be unrea- sonable to suppose that not all of the remnants of the s and d atomic states are itinerant. Perhaps a sizable fraction are not available for pairing but do act as an intermediate device to enhance coupling between the paired electrons and the ions. This would most likely reduce the density of. available states but, at the same time, could substantially increase the effective electron-phonon interaction. Such a concept is worthy of consideration and, therefore, will be explored in as much detail as present knowledge of the group'V and VI transition.metals will allow. CHMPTER II THE THEORY This problem is primarily concerned with the elec- tronic structure of the group V and group VI transition metals. As completely free atoms, transition elements are characterized by two partially filled atomic subshells (Bass, #dSs, or 5d6s) as opposed to the usual ggg_(or com- pletely filled shells as in the case of the inert gases). The number of electrons residing in both the s and d sub- shells determines the chemical and metallurgical charact- eristics of the element in question. For groups V and VI, the atomic Epnfigurations are: ELI V-3cl “rsz Cr—3d S L+ES' ' 5' Nt-LH s [Via-LN 55 I 3 2 W L) G's-Ed 6S -5d 652 group V' group V14 The objective is to determine how these electrons in par- tially filled subshells (valence electrons) interact with valence electrons of adjacent atoms in a metallic solid. It is often stated that the unfilled s and d subshclls transform into overlapping s and d bands in the solid state. moreover, a frequent assumption is that the relative propor- 14 15 tion of electrons in each of these bands is the same as it was for the partially filled atomic subshells. A more rea- listic viewpoint, however, is one proposed by several au- thoritieslg’20 in this area which is that the valence electrons collect into two fairly distinct groups. One such group is quite localized and may possibly be repre- sented by some linear combination of the former d orbital wavefunctions. The second group would be an itinerant con- ductionelike band which could probably be represented by plane waves properly orthogonalized to both core states and the localized states of the first group. It seems plausible that these two bands overlap with the localized band being the narrower of the two. Moreover, the number of electrons per atom residing in each of these bands would, in all like; lihood, not be in the same proportion as the relative num- bers of electrons in the s and d atomic subshells respec- tively. Such are the foundations to be later utilized in setting up an electronic model of the group V and VI tran- sition'metals. If one is to base a model on the concept of two over- lapping bands where one is localized and the other itinerant and conduction-like, then it is fruitful at first to inves- tigate how such a configuration might affect the pairing interaction in superconductors. The conditions for pair- ing are most favorable with a strong electron-phonon l6 coupling constant and a high available electronic density of states at the Fermi surface. The question then arises as to whether all_of the valence electrons are available for pairing. Perhaps either the more localized electrons or the conduction electrons do not take part in the pairing process at all and hence one must use a density of states that accurately reflects this reduction of available elec- tronic states. Certainly, then, if a given proportion of the valence electrons do not contribute to the pairing pro? cess, transition metal superconductivity may be enhanced through a mechanism.other than a high electronic density of states. Another possibility would be to examine how these non superconducting valence electrons might tend to modify ion (phonon) coupling to the superconducting electrons." If this modified coupling were quite sensitive to the number of electrons in localized bands, then the difference of a valence electron per atom between transition groups V and VI could account for the gross contrast between their ability to superconduct provided the extra electron in group VI were in the localized d band. The model to be used will be one in which nd valence electrons per atom reside in localized d-like band states and n0 valence electrons per atom form an itinerant conduc— tion band. The question naturally arises as to which band 17 contains the superconducting electron pairs. In the past, it has been quite customary to consider the itinerant elec- trons as those available for pairing 21’22. The localized states are then usually treated as either intermediate states for electron pair scattering by the lattice ions or as states acting directly between the pairs by providing an additional attractive coupling mechanism. McMillan15 seems to avoid this question altogether when he uses a density of states based upon the electronic specific heat 8 . Al- though the argument for d-like electrons over-screening the conduction electrons would enhance superconductivity,‘ it makes little sense to consider electrons with a suppos- edly high effective mass as capable of screening the effec- tively lighter conduction electrons. However, it would not suffice to say that electrons with high effective mass can- not be screening agents in the case of the ions due to the fact that they are also strongly coupled (source of high effective mass). Whether or not screening of one band by a second could enhance the electron-ion interaction of the first should be analyzed. These arguments give rise to two possibilities: (l) The conduction band electrons form superconduct- ing pairs. The d-band electrons screen the motion of the ions in such a way so as to enhance electron-phonon coupling in group V and diminish it in group VI. The difference in r1- l8 superconducting properties between groups V and‘VI is attributed to the difference of one electron per atom in the d-band. (2) The d-band electrons form superconducting pairs. The conduction band electrons serve as a screening mechanism to alter electron-phonon coupling and reduce the effective coulomb repulsion between d band electrons. The high den- sity of states of the group V d band at the Fermi level greatly assists superconductivity whereas a small density of states in group VI all but destroys it. Again, there is a difference of one electron per atom in the d bands of group V and group VI. One of the most crucial portions of the model will be the structure of the localized electrons. One must know their density of states distribution, energies, and, if possible, their wavefunctions. unlike core electrons, the localized electrons cannot be associated solely with a par- ticular icn within the crystal. The characteristic that distinguishes them from conduction electrons is that their distribution varies grossly within the structure of the‘ unit cell. Localized electrons will tend to bunch up mainly between nearest neighbors and next nearest neighbors_respece tively. Such electrons should'besensitive to thetnotion of the ions even though they cannot be associated solely with any one particular lattice site. Since both group V and 19 group VI transition metals have BCC crystal structure, it is best to examine splitting of d-orbital wavefunctions in a cubic field before attempting to use them in constructing localized band states. Basis functions from which the localized states can be formed are as follows: xy=\/% {fies-Kerb? YZ :V‘WQEXVGDCPH )2/1(e,d>)} z xm 3X‘1/ (e are the orthonormal spherical har- 2. monics forming a basis for the angular part of atomic d wavefunctions. These new basis functions are all mutually orthogonal and therefore must span the same space as X . The purpose of the particular basis becomes quite clear when one examines the full cubic group OhfiTible 3) . Z 2 Z The functions X —->/ and 3 Z —R2 form a basis for the two diminsional representation E 8 ant?< y) y Z J andZ X form a basis for the three diminsional represen- 20 MN macaw sonao>GH on» uH H cams: Hadinno asoam oHoso HHsm moss H- H H- o n me said: H H- H- o n Ha Nada H- H- H H H we I I q I 0 0 m H N m mm «NM N» «x H H H H H H 4 mm soHuwpsomoamom mHmom zoo om mom mom m mmoHo *.o moonw oano on» no oHnmu nopomaono .m oHooB 21 tation 1—2 8 . This means that under spatial opera- tions which preserve the symmetry of the full cubic group, these basis functions do not mix but transform only into linear combinations of the members belonging to their own particular basis. For small cubic crystal fields, the five-fold degenerate atomic d-levels would be split into a doubly degenerate E- level and triply degenerate T 'level 8 28 respectively. Stronger cubic fields (as in a crystal) will further split these levels in every atom of the crystal thus forming two distinct bands. These bands are in accord with Fig. 2 showing electronic structures within the BCC primia tive cellen. These structures as shown in.Fig. 2 can be inferred from examination of the angular dependence of X2, Y2, and zx for the T28 case or 12-Y2 and 323-32 for the E8 case pro- vided that the wavefunctions are formed from proper linear combinations of the basis functions sc.as to preserve cubic symmetry. Such wavefunctions will be referred to as hy- bridized T or hybridized E respectivelyzn. Since, in 8 both groups V and VI,there are insufficient valence elec- 28 trons to fill both T28 and BS states, it is natural to de- termine which of these bands would lie at the lowest energy and hence begin filling first. Because it is generally true that nearest neighbor correlations are stronger than next 22 .hhumfifihm mzoum.~mmv 0N9 mnw>oa moonuoon How HHoo T>HuHEHHm .U.U .m on» :HnuH3 huHmoop conuome Enfiflxmfi mo mmoH< am am». .m Tasman Itux (”0' ’fi 5‘? pl. 1}- A 25 nearest neighbors and so forth for more distant ions, then the states representing correlations between nearest neigh- bers might be expected to receive electrons at a lower en- ergy than the next nearest neighbor correlated Eg states. Although band structure data is seriously deficient for groups V and VI transition metals, one can investigate the band structure of one particular metal (tungsten) in hope of getting a handle on the essential features of band structures of chemically similar materials. Mattheiss25 preposes that the band structures of the BCC transition metals (group V'and VI) should be quite similar to one another with the Fermi energy being determined from the total number of valence electrons. Upon examination of the theoretical density of states curve as plotted by Mattheiss for tungsten (Fig. 5), there appears to be four fairly dis- tinct peaks with the lower three appearing quite similar. Furthermore, it seems that for up to six valence electrons per atom, the lower portion would provide a sufficient num- ber of states to accomodate all of the valence electrons. Although this problem will be reconsidered later inthe discussion of the filling of the bands, it will suffice to say now that the T states will be used to represent the 23 localized electrons and the E8 states can be disregarded for the use of six or less valence electrons per atom. .... 2” Proper tight binding wavefunctions representing T28 states are 8 .. J —‘ 15 £520 A1553— ‘LW 26 tar-mm) 1 where the sum.over 5 is over the entire lattice and 49-400 >yz and AHA-Zx 2 The convergence factor }1(h) satisfies the normalization ‘9 2 . condition that j f0) Ir‘J- :: l and is best guessed at from the atomic d-electron radial wave- 26 functions or the band structure itself. If .F(V9 is normalized as it is here, then these wavefunctions represent a continuum of states in k-space and, hence, must be used with an electronic density of states. Later it may become necessary to modify the spin density of //<:(I) to reflect the occurance of antiferromagnetiesgiructure in this band. Such a modification can be included in the spinor factor :><§I) ‘which may be set to vary spa- tially from spin up to spin down with a period one half of that of the lattice. Representation of the conduction states should be much simpler. For one thing, the density of these electrons should not vary grossly within the unit cell since they are 25 free electron like and hence are most easily represented by some variety of plane wave. Ordinary plane waves would suffice but should require that a large number of thembe used in order to cope with strongly varying potentials throughout the crystal. On the other hand, augmented plane waves (APW) are particularly advantageous in the sense that they are ‘most frequently used in band structure calculations making it easier to apply the results of such calculations to this particular investigation. Orthogonal plane waves (OPW) of- fer a third alternative which should be most compatible with using tight binding functions to represent the localized states. E's-i: 7434:) — \__i}____ _ V—Nzg|k> Ck—L-SCIS -"C 3' WCk +' where ‘ 6k are the Bloch energies relative to the Fermi level and V is the screened coulomb interaction as a func- tion of the distance between two electrons. using Nambu's notation, define new field operators in terms of column and row matricies . w_ E“ SDI: Ct C .k: C‘kil' 3 l; lgl jkl where the; arrows T and ‘L denote spin up and down re- spectively. Next using . 5-2 :(3 _i , express the Hamiltonian in terms of these field operators noting that the spin is handled automatically. 2 32‘ :;5:Cn:n=+C:C.kfi1)—§:ek fl“; gj+ZeK , ,Zik‘( 1”“? kka5 kkscc :gfl,( “W” ”)(CETCW‘CELCM) :kZig ...(a-°y+% ) (fig % kk’. Z|k>CH3Cmgl@1st 2!; lg’s.s'l -1Z<‘<-U<<+l|V|k>l 1&sz -L.f Clan CanI + CH: («1+qu ’JCwCELCIJCk 11%;}qu =;z<_k—L|\k>(SI£lOE@fififi), k,k'.L 33 or, neglecting the term Zék (unneeded in Dyson self energy), H =25 %+Ggqé+;€s@w%~) fig; 1% 2% 6% Z«+%gE§E-U :1; where the brackets <--- > denote the grand canonical aver- age, T is the flick time jjfering operator)6 : T ‘i ' + CP%‘(T):C H cp‘fum T: C ‘H’T(a.1¢+a_$e) , and 74,: 211'th . In the Na‘mbu scheme the finite temperature electron propagator becomes ;) h Gash): —; {ENC-$05 > } “r 3S .>Hm>wuommmmu.mwmmwooum coconm can nfioasoo wumoflcwfl mmcaa a>m3 can pounce .wmuoao «How on mcwusnfluusoo mfimummwo ucmuuomfiH O I l O O O O O I Q Q ‘ ‘\ .0 00-.....- .m ousmflm . as _ :éfCMT fl-éHTCPIJomiHqugawchT) 7>o :. £1,3qu eiH‘T 475* (0)6"HT4>:I°°‘(08>4 T, Té 0 gamma "”e' H" o >‘e wane» Z < I e :H I > . ’r‘>° dz; Lnem‘rZ—QK—FHecHTfiAO) €-'iHT¢;,(o)\I> " ' Z-.<¢I (MI) “‘9 ifevnTZCIT :2, Z&; I 6 PH I > . 7no 2! ’ev‘ZfiIEF r18HTé°>l>AT * ‘ Z *9 f‘c} 14786 JE: +( E EJIZ‘ d)?‘ 0)“) <5‘ 493413“) cl TH”, Z (ilee ’H|‘> ‘ 2g. >2ngKHEVEJKII-(flw-(E‘EEI) z Z;€"E‘ x 4|¢gé°4><| flaw» i> (PE .3 z _-L “'J—l’z- Nc ov P8 37 Define the spectral density A (k _Iév)=(1;(_§:)ZJKJI 43:”; >LE86V‘E +3) ‘9 O where (1 =k_‘k/ . k-IZ,£VA)=[°A¢(JS'EIV) 7—73.7— —‘ 3; $77-14? 10 Now substitute Eq. 10 into Eq. 7 and transform the sum over h’ into an integral representation using the J I residues of poles within “an. infinite contour C to pro- Ject out terms of the infinite sum, i.e. integrating over the contour will give infinite series of residues of the ‘poles contained within thecontour C (see Fig. ll). Functions tanh —E-t2—— and (1+Ct‘921 , one has the residues for values of F2‘ - o dd ¢ 1T pz _Lt 2—ch _+c;_‘_“_'.__ 1 ResI1+ 6* ”mm: 73— 7;- 62/: 42/2( RestahIILzz_ —Lm (Z- L 11:; 2"? (2-75“ 626 —6_‘ e Now, making use of these residues, let LWhI—, Z to get (tam-2%: gtfdtzIQIZIépzfaRt-m Av+\fl<_k’)ttnkgjgz n 1 {(1 w..-z -‘VX1+€”}+( I Mr 2+ lWI‘HCin 1}} 38 m_x< 44mm .m\: 660 um modem mfiom mcfi3onm wasam Xmamfiou ca musoucou .v mhsmflm 39 ,/l I :When the contour is deformed to C , two additional poles must be included whose residues are RQSI Gaze) _ GIEIIwhtv) W " mam-£17 :‘GIk: wni'V)_Cr(k.iwni’1/) 1+6 “M " 1- e“ m, usinsGIk’,w'+iO—G(k’w'-i8)=z «'lmGQéw') and the fact that the contribution on the circle at in- ’fAAk E. I finity vanishes, the self energy becomes 2 (kLLM-‘1472;loilmc(kw , l 1 , _ {(fwh‘ w’ “ W)+(éwn-w’+ 10(1 +3 WBVJ'VIV-zgk ”4%?de Q 1'; ‘V I". _ z z _’) . 12 OEGQ, m )1+F(ek- UM v’Bo-Igggtl Nk IsVM‘V Letting i a)“ —-> a) +1 5 , the self energy can be. expressed in terms of the real analytic variable 00 (see Schrieffers, P. 197). Its most general form is :05» “0):{1’ZIkuw8w1+¢ (kw) ' This expression must then be substituted into the free energy yielding a set of four coupled integral equations. Several steps can be taken to simplify these equa- tions. First, set the quantityé’flgflU) =€k+X(k)w) and let AQSJU); CMQQ’VZ (kW) . Then the following criteria can be used to reduce the algebra considerably: (l) X(k,w), being due to nomal coulomb pro- cesses, is slowly varying for changes in (J) of the Debye frequency and,hence, independent of CL) . (2) ZLkJ») and (Pagan) are important ‘mainly at the Fermi level. Since they vary quite slowly with k , one can set |k‘:‘<}_ . (3) Since transition metal superconductors are us- ually impure, it will suffice to determine the spherie’ cally averaged self energy because impurity scattering reduces the effects of isotropy. (4) The chief contribution to the sum overE is when 6.2/3 )6 (DD . Following essentially the technique of (SSW) the sum over E can be replaced by continuous integration with the following definitions It- 0 41 2:... 121%,ggg—TTI1‘3 [ale-:34 “k 15 where m dek= If dk’ . Since the chief contri- ’ I bution to the sum over k is near 6k =5 ’R 0),; , the limits on the integral over Ek’ can be approx- imated at + 00 respectively and, hence, by residues fdéfl' C(19)) oériTrL—EUIZM’) 7)1_¢(w’:-] 16 k z szI‘Z‘Iw'F CM with the one pole being above the real axis. By avera- ging Eq. 12 over the spherical part of k (and remem- bering that ((1).. —9 (Di-i5) , Eq. 15 and Eq. 16 can be applied to yield for the phonon part of the self energy at k 2k} III r—f Ive—em 4521; [Ad (lg-Ky) 'ngch Zia-wLV+zJX1+e‘PV'T+¢U-I”'*W‘XP e "55] I I zIr—Ifi—I II I" I—t—EEI 42 where (Us: w+(—1):'V+ i 5 . Equation 17 can be simplified with the transformation of the negative part of the integral over to the positive portion. Then :35“ w) :EEWkIWWFKWWQ {Irv—If III III-III w where the upper signs belong to factors multiplying the 0'2 component only, lower signs belong to factors multiplying the I compOnent only, and W9=M7flz W712i: )NlS-Is, twiwcs gm; 19 where the Fermi-Dirac distribution functionf-(w) The coulomb self energy is )=jg?\<)a‘_[g$ImEvg G(ls'w'kgjwkk'flanhgg'. ' 20 An approximation is developed for the 0'; component of E c. in (SSW) and, hence, only the necessary re- sults will be stated here. #3 These are the death-N(oflfifie £333 Ly Ia.» where the "potential" Lax; is approximately U = \é__ I c’ 1. + N(O)\/C (.h (6%0 00,) and is the average of the coulomb interaction over the Fermi surface. This finally completes the develop- ment of the finite temperature gap equations. It is apparentmthat, even in reduced form, they are intractable and, thus, some further approximation is needed to bring forth a reasonably accurate set of solutions. In practice, the gap equations can be selved by com- puter for the zero temperature case. Finite temperature properties can often then be accurately found from the zero temperature solutions by rules obtained with exper- imental evidence. From Eq. 18 and Eq. 20, the zero tem- perature gap 1equations become. AIwFawo Jaw’ RQ{V_TA_(_Z‘”_')_ }K( on ('w) _N__(_o) AI (10') I 2“» U R {Maw 21 22 23 MI and [1" Zhui] (mi-[211‘ Rafi/“7%) \(_(w‘,w) 24 where Am”: (hwy/Z (10) and N(o) is the electronic density of states at the Fermi level. 45 . .mump wusuosuum Icahn Eouw soummgsu How moumum mo huwmcoc Hmowuouomna .m enamam .mecmmo>m. >ommzu - h. . n._ n._ 3. 3 ho no O 3 _l a l. .v 8 mm 3 m d , o E N. W. m w. 0 H . 0. ON NOLV 83d SUBBOAH 83:! 83111.8 NldS CHAPTER IV THE ELECTRON-PHONGN COUPLING CONSTANT Derivation of the finite temperature Eliashberg equa- tions begins with the irreducible self energy {III-é- ErfiGIk’MIo;{ZnggIzQII-I’I III +VI—IEI 30 Green where G and D are the electron and phonon Nwmbu functions respectively,‘v is the screened coulomb repul- sion, and 3kg)!" is the screened electrOn-phonon coupling constant. The most important and also most elusive quan- tity to determine theoretically is 61%,; . This cone-"1 stant will be derived as rigorously as possible utilizing Bardeen's self consistent calculation including exchange32 and other necessary modifications needed to tackle the previously discussed model of groups V and VI transition elements. Consider the displacement operator of an ion at point due to a lattice wave of wave vector db and po- . .,, i" m 40.9. Sfimzzefi'fi—(dewa “Ian e B!" where gang» is the polarization unit vector for wave- 47 vector % and polarization index 06:1,2 D 3 The normal coordinates are given as QT: :VZMwT‘ a“? (Iaj Tc where “1‘ and 0.1:: are the creation and annihilation operators for phonons of wave vector x and polarization index DC . Thus the lattice displacement operator at lattice point becomes €13," SR m =2- 4‘er écrce (av: +£_T¢)+hc Now assuming that \/Cr) is the external potential at distance I: from the center of a single ion, the poten- tial at point B due to displacement of the ion by the amount 83 will be approximated as V“ “’83) . Then the change in lattice potential at point I for a wave or wave VQG‘COI‘ i and polarization index ‘3: 18 Z [Sfim' {-WII III] M where M is summed over the entire lattice. Substitut- ing Eq. 3 into Eq. 4 gives 48 , ITEM \/;=‘§1’ IANME 60?: e (QTId-Tt)§3awr'fim) +hQ This is the so-called. rigid ion model of Nordheim”. Let - (Lin—BM) \/(r-B_mI- ZM Q . .L Then, eel! e‘IQj) Em Iéf szI “I” “'“I‘ IEVI‘L-IIE E) + Inf Since 2; Gui-L). Bm_ =ZN 6b$+§h where GUN is the Ihth reciprocal lattice vector. Then gym-1h ‘4' " 2mm I‘I“+“I‘Z‘4+II INT“ ( The total electron-phonon interaction Hamiltonian for scattering of an electron from state k)$ to state k S is Just the matrix element of VT, be- tween states 195 and E,S multiplied by CEIS Ck; which destroys an electron in state .3 and recreates it in IS 3 . So one can then write \0 (A 49 He-Ph —gflE¢ CESCKS (a? ‘+a'"f‘ “I letting I<'-IS/=°_1_)+.Gh glqgf‘i‘I/‘zwr—WIIs-EIQQILE (ISI Q. - )r Subscripts i may be replaced by kk since q is unique for a given kt; and the sum over reciprocal lattice vectors is inferred in the sum over k _|_< be- cause Gh is also unique for a given 13 E in .umklapp processes“ . The best way to determine 5kg“ is to find [VII k EL: <“v O h '9‘ -I j which gives MOSIS’Iec: §I§J§°‘V “T“ 3 . To begin the Bardeen calculation it is necessary to return to 10 H 1 12 13 5O ‘\€;é19VAEEEEQ&Eé;éste+UCfifIzéggfgliag;1f5C7IEEISIIIC 55 01‘ IgffI/z—N—IM—Zyze I”; TVI/(r-BwIJIIIC ‘ whose phonon number eigenvalue has been determined in phonon occupation number space. Let ice“ upon a conduction electron state ”5) V19- :{v‘m’ILeTQw +\[I"dlk-oo I» where \/‘“’— . Doing the same to a VJ Is>= Z{\/_ d- electro on state I_I<_> (hd) (Md) e IISIIIQ I’Vvl Ii operate 49>} u H 16a 16b 7 , H 51 where \/‘I-—a;:d):<.|<+ j+gnI\/¢1;)¢II£> 18a \[_.f““)=<|_<-I-Gn Wang) , This represents the action upon a core electron by a 8b H lattice wave of wave vector i and polarization at . To be self-consistent, one must also determine how the electron states are perturbed by such a lattice wave. For the screening effect, expand the perturbed electronic states Ik» and |K>> in terms of I I? no @1454anWQWIroIk-fig} H 9a H IK>>IK>+XIA (K Is + I 334 Ing— Q} 9" h “12" n ‘32 (.3 I IS “I to first order where a, b, A and B are small pertur- bation coefficients. These latter two equations, when summed over all of the states of the band in question, are the net response of that particular band to first order due to the motion of the ions plus the screening effect of the other band. Next, call the Hartree-Fock a potential energy of the d plus conduction-electrons \4‘. 0 Operating with \é upon a perturbed state using 52 WII>>=EXIEIIIIS>>->IE>I +€28K’){|I>>->lK>J V; |p>=e2XIKII>I IQ sflkImsms>->II>I 20b with the specification that the eigenvectors span spin space as weall seHilb rst space and that the sums include spin. Call the Harte -Fock potential energy in the per- tu Mb dl th ”\4; .. n I; Il_<>>= e ZEII<>|D¥<>IE>>I KZZIKKKIIfiIKIIIIsI-IKIVRIQIKII III©=€ZXKI<I§>IKI§IK>|K>§ MEEIII'IIII’IIKI-IEPIIK>IE>>} 53 is the effect upon a particular bandstate perturbed by the bare ion and screened by the other non-core electrons. Now substitute Eq. 19 into Eq. 21 and subtract Eq. 20 to \; |I>>=(\I*-\;°I\k>> i.ezwgag‘ghmwQIR'IQmTIMIIR‘IIWQ 4:;th (IE-1'94 R"Ik’>+I>oD-II’IIII€IR'IE‘I'EII‘II {a}. New R“ III» +ayIIII+Ar‘KI a -{A"ZD.(KII|K>>+ 814mm bl rap—GI] “2 n5 5n and an identical expression with K and k, a and A, and b and B interchanged. Equation 22 represents the change in potential energy of an electron within a given band interacting with a bare lattice wave of wavevector a and polarization as when properly screened by all of the other electrons in both bands. It is now necessary to solve the Schrodinger equation for the perturbation expansion coefficients. For the un- perturbed lattice, Schrodinger's equation is I—IOIDIVD IIgI =EOIISIIK> 23s for the conduction electrons and HoIK>+IAIKI=EOIKIIIs> for the d-band electrons where \VL includes all of the potential fields acting upon an electron within a particu-' lar band in an undistorted crystal. Now, for a lattice perturbed by a wave of wavevector Cb ‘ and polarization 0C , the Schrodinger equation ficomes H. |k>+(\;+vg +\/,IIIg>> =EIIsIIIs» 2“ for the conduction electrons and HOIK>>+(\;+\;.+\;IK>>=E> as 55 for d-band electrons. Substitution of the perturbation expansions Eq. 22 into equations Eq. 24 and then sub- tracting Eq. 23 and including the assumption that E(k) - Eo(k) and E(K) = EO(K).gives ‘ Ham—Lime»tetanus—T9} +\/,, fia1.(kn)lk+fb+Q>>+bT(k"‘lbib-Cw} + :{Wm’lkjg+g,>+\[‘l"lk-i-Gn>} +_ af‘lm \[mc [kRTZGQJerWfl n: M ‘1“ +bgqu'fwdlk)*\4:d'k23-ZGS}] + ea;(k)lé}ll£n)+bj,(EMF-1&4 R"|k> af(gn)]|D -e‘ 231; )ka «E + lfilfi‘l|1§>+b*f«“)+LTAEA+B§.(En) who KIR‘HLGrG; + 81m (KI R‘1 “.6 are] k> .eag ”WWW R-1\k>+B’;.(K;)\K+1+Gny(m\h wows ug— W} .5 56 for conduction electrons and an identical expression with K and-k, A and a, and B and b interchanged for the s -- d-i—ban‘d ‘ electrons. Now operate with bravectors lts> + meg’kaW‘cMQhKE‘R1|k+ib+§rv>|k> + LllkakfipGnlflélRam-Gab] .92 ;( k’)[a’fr(lg’m)l_|<’> + aidk'ka'TQrKk’l R‘Wkflk’tycé 57 5%. (km)|E- $199] +6 ;(K')[A}(Kn) I19 +81} (Km\lk>] +8 gm [A34 Km1IK> +B}i|@m\\K+$+Q»> +ch(_K;n)} WM K K— swam- 19.)} Now substitute these expressions for the coefficients into Eq. 26 and Eq. 27 and solve. The result is 8 '0 MLK.l<_+ —(K+§+Q,.KK+TQM "31'9““ MKK K ”i GmIa 1 WORK?" G JKKW‘QRK KleIfi KKK) -IK>I M + Lug qug +(3°Q(IK+i+G»I) mK Kira.» J, . ’_ .. El“ ) HK m QIK- $39)] +€Z’) WK r9 |llg> —(K+$G»Kk+i+G~JR IIQID) 6O +WKETQKK’MM lfiflp‘m ‘II£>) KEW‘KKI’K GIIK> ~0é+l+§nI ’3 62 Next, multiply by Q‘a‘BM and sum over all lattice sites a . (G‘U‘B :3 ‘b Ewe <1lvv=N 3* Equation 34 can now be used in Eq. 16 and Eq. l8 andthese equations in turn are substituted into Eq. 32 and Eq. 33. using Eq. 13, the expression for the electronéphonon coup- ling constant is EKK 3a.,«=‘\/2§_—'EKTY€,¢'|K’>B *{gfiwmwgfllRWK+1+G>1K> “ll§+°K-+GK>)} \ 63 ‘{( gkfiékai QMQ+TQMKJS 32 GMIFflIOID —\K’>)} *{fiWWanKKIR“|1K’-fl)-CK>IK>- —(k+1+QKKKIR”1|k)|K'- K993] + 62 SK)" mm KG KKKWKKK) ~(K+}*GKI|K>)} {(ng'i-KK-mKKK *GKKK IR"| K KKKK ~IK'+‘])+GKE 64 K; K. W) (ISquQKKK-wAR'flDID K —(K+K+QK|\K )3 K‘KKKKK—a.) '<‘-<*°fi*G)|Q§|R"|K>\K’—) ,5 for conduction electrons coupled with a similar expression for the-Ku‘d-Kbandm electrons. Only the first term in Eq. 35 represents the unscreened electron-phonon_coupling con- stant whereas the remaining terms under the sum are en- tirely due to screening. The problem now is to determine the bare ion poten- tial \/fikfl . This potential may be broken down into two parts, one of which is the nuclear coulomb potential and the other being due to the screening effect (including ex- change) of the core electrons. Since the core electrons all lie within closed shells, their potential is spheri- cally symmetric about the nucleus and hence the vector notation (r) may be dropped. Since five or six atomic electrons enter band states in the transition metals in 65 question, then the remaining core will be quite small and should resemble their atomic counterparts quite closely. In fact the use of atomic potentials to represent core po- tentials may very well be a better approximation for tran- sition metals than for normal metals. The calculation of these potentials is largely unnecessary due to the efforts of Herman and Skillmanas. Essentially, they have deter- mined a set of self-consistent Hartree-Fock potentials for neutral atomic elements. It is from these numerical poten- tials and radial densities that the necessary core poten- tials must be determined. Up until now the possibility of phonon induced inter- band transitions has not been considered. ‘The reason so is that such transition probabilities would be proportional to matrix elements of the type and . Since the potential V6.) and hence its gradient is spherically symmetric, these ‘matrix elements can be non zero only if the bra and ket contain at least one identical spherical harmonic. However, in specifying the wavefunctions, such a possibility was eliminated. Thus, for this particular model, interband transitions are ruled out. CHAPTER V EVALUATIW OF MATRIX EIIEMENTS The purpose of this section is to develop specific expressions for matrix elements of the type < 302’) l| 143:!) |.g(r))|fi(:9> and (EMIVVYHI 31m) where WV) . is the gradient of the bare ion. potential. The form of the wavefunctions will be linear combinations of atomic orbitals and plane waves from which any ortho- gonalized plane wave can be formed. Part A - Fourier Coefficients of Wavefunctions Consider a linear combination of atomic orbitals MM 1 Z ik‘U ~ A v- ="'" e l (II-1’ ) 1 kL) W 3 4)" m I where N is the total number of lattice points, r: is the Jth lattice point, and 43 is a given atomic wavefunction having quantum numbers n, l, and m respectively. The Fourier transform is m . tf’I Mm CLG)-‘=éi}nsf€ A (NA: 2 66 67 Performing the integration, this coefficient then becomes Mm '(Jr+f)r M ‘ Q, (E) firm—XE ’Iuiia) 3- Mm Mm - Ck (P):"_‘(q (ZHP JHIUD) 9r) 4 where M if- I: 13:36,) =f E (PMMQ‘H’: « 5 Orthogonalized plane waves may now be treated using COK‘E’izifmfeiP-t 4%“! ‘ 6, aw “7. ‘fr Mk 0C8): 7211?); e- e _techrt-r n mgr'nfik “Mr-133“ The integrations are then performed to give .1 . W V‘, %n a! (m G(Ekg—Tjs he; [1 X1“?- “090 CEU”. 8 hca-i(-ohl 68 Likewise, proper linear combinations of d-‘states having T6 symmetry possess the coefficients #122 “2- -Z i Cx y:"(E) Ck _ C4—— ) 9a '5 V2" n21 h2-1 CZEF-gfi r—CL- 3 9" and C1121 “2-1 + C }_ 9° C7649 V‘é‘ . To evaluate IhUM , let (phlyir :RMU') ( and €{E-r: 5‘77 Z(Z(“)+1)V:, 3‘.(PV‘)X(Q’) 11 ('a o 10 Transform from the primed to unprimed coordinate system letting 6? be the angle between 2' and Z. Then (see EdmondsBB), ( , ~ ' Y.0(6’)=XX.M(9'¢) (fitm'fi! PM (COSGP) 12' m':-( 69 Thus ( ( (U- 12€5A=V91T Z} (2‘ .)1/+1 Lars”)! I'D.“ ((056?) 13 '34 «7.1,, :]3/[R“ (HjtJPfly (e,¢)>’(a,¢)r Séheclegié 4'. Performing the integrations over angles, the final result ‘ _ ml 1.?fim,)'-"47((:h1)foh)/P\ C c as o J R (r) 14M '_2 4 b- 11; If the radial wavefunctions are assumed to have hydrogen atom-like form, then KW 2mm ‘ 15 (‘0 h ffe i f h e t on i were/(istheectvemassoteelcrz, EM s the state energy, and the coefficients, [3“ are set to give the best fit to the currently accepted radial wave- 70 functions for the atomic state in question. Equation 15 can then be substituted into Eq. 14 to yield the required 1 MM (see appendix A for the first few co? efficients). I Part B - Evaluation of Screening Matrix Elements A typical screening matrix element is of the form M. E: {R (h’)'uk> ’IE> +Lb ~ag+cy$3 +{E 1 (aft—1,9,3} {\D ‘. +1? QmKK- VG”! R’llkflk’) +|k"l'Q~>>}] 76 Where 0 is the volume of the unit cell and the sums are over all occupied states. Substituting in the correct ex- pressions for the matrix elements from chapter'v and using C,<>r>--(N1’TW5M 10’, d (SEMVGM alga” wit él‘ H<_+:b*§.> +L‘1TQZ;( — k’) Eliza/mtg“ X Ban E10; wit" {2: We )1 m we mks s5? -I(k+1,+G,\I(E)I(k}I(E+T§»D 5e.,£rmlk_k+é -GB law-firm} 77 www- ,- am)1 S — 4 131+ Q, '{r\/ 1- COS 2(Sén‘1Ke72k Where [0(E}) is the density of states at the Fermi energy and Kc represents the correlation momentum. Such momentum represents a.minimum scattering momentum associated with a correlation hole about each electron. Electrons at relatively large distances from each other are effectively screened so small momentum transfers do not take place between them. Substituting these integrals into Eq. 6 yields 3m ‘3' [1+ ;§B&in d {'21:- + ' '1 (8k:$§%|2)(“%fil%%)} IT:+Q,II<,(FV(E 9‘“ (Nil for the electron-phonon coupling. =00: 61)::ng \O 81 An expression for the superconducting transition tem- perature is derived using a technique similar to that used by Carbotte and Dynesl3. The desired quantity is the coup- ling times density of states function which is ccz(w)F(w)= $52622 Biflg It, '%*§n|| 2%) "m [141“sz where the sum over modes of vibration is implied. A com- puter program developed by Gilat and Raubenheimer39 pro- vides RWz‘i'j—EL Zen? SW14” w) n in histogram form for the phonon densityof states F(m)Am . Comparing Eq. 11 with Eq. 10, one finds that, for each value of the phonon wavevector q, the proper_value of c2(mo must be computed for all possible umklapp processes and then multiplied by the number of states placed in each histogram channel for that particular q. Thus, the histogram will then reflect a2(co) F(m) Am for channel width Am. The coup- ling factor for phonon wavevector q and the nth umklapp 82 process is then EPH’I‘IIQIEC—LE} l we“ 12 A computer subroutine, EPH, which calculates Eq. 12 is added to the Gilat and Raubenheimer program (see Appendix C). Another subroutine, SWEEP, which constructs the phonon density of states histogram, is modified to find the coup- ling factor for all umklapp processes involving phonon wave- vector q. In the same subroutine, the number of states placed into each histogram channel, F(m)AA(n, is then multi- plied by EPH for all umklapp processes to construct a new histogram containing c2(m) _F(co) Mm) in each channel. The main program.then calculates A: 2[é‘ma(wlf(w) AI» 13 for comparison with McMillan15 together with the transition 83 temperature using €‘[,\3'_}%L(tig.ié+2>.\))\) } T: 365 for strong coupled superconductors. The value of )3 repre- sents effects of electron-electron repulsion whose values are taken from ncMillan. Plots of both F(m0 Am and c.2(co) F(ai) 4(a) are constructed and shown in Tables 6 through 9. Two crucial quantities are Kf and Ké. Upon examina- tion of calculations for Kc for sodiumuo, one can reason that,for the greatly increased density of electrons in a transition.metal, Kc = Kf should be a fairly good approx- imation. On the other hand Ff is freely varied about its free electron value to see what results can occur. Since, for the case of a real metal, Ff. changes from point to point on the Fermi surface. Other parameters utilized by the program.are the Debye temperature on, the inter- atomic force constants, and the coefficients for the func- tions representing the ionic potential. Such data is 84 Table 4. List of potential coefficients and force constants for tantalum. .— Ionic Potential Coefficients Force Constants and Density of States (dynes/cm) ‘mr 1J725xlo9 ,1 16983 AA a1 -9.774x10‘2 31 11201 10 a2 3.894x10 ,2 1182 a3 -4.77ox102° 32 1u83 a, 3 .oo8xlo30 "3 3546 a5 -1.069x10”'O a -5#27 he 3 a6 2.227x10 Y} 1943 a7 —2.698x1058 ,4 3577 an -718 n(er) 1.52 x1034 ,4 983 ”'5 '1’93 35 812 ”6 -3705 86 134 "7 558 B7 -237 Y7 106 -683 85 Table 5. List of potential coefficients and force constants for tungsten. Ionic Potential Coefficients Force Constants and Density of States (dynes/cm) E 1.780x109 ,1 23000 a1 1.061x10”1 31 19200 a2 3.926x1010 ,2 #7300 a3 -4.795x102° 32 800 a, 3.023x1030 "3 3200 a5 -l.076x10“O B3 1100 a6 2.248x10u9 y} 4900 a7 -2.734x1058 o(E) 2.56 x1033 86 summarized in tables # and 5. The final results include Figs. 8 and 9 of which Fig. 9 agrees quite well with an experimental investigation by_ L. Y. L. Shenhl. The dimensionless coupling parameter, A, and the superconducting transition temperature stemming from it show reasonable values which upon freely varying Kf, can be made to agree with experimental values. 87 .ESHMNMM—m “Mam—ownmwvsosonm mo muflmsoo o.m on 06 on o.~ o._ - u q d , q .m ousmwm 0.. O.N QM O.¢ I.-zu 3,0: mm 88 OK .soummcsu Hon mmumum sosonm mo husmsma 2.. so}... .5 magmas 0d 0.0 03V QM O.N O.— . sq I1 q . hflJ-I -o._ H. W H O h. .6 n1 -o.~ 2. .Edasussu Hem moumum no wuwmcop mesa» madamsoo .m munmflm .2223... 89 0.0 0.? 0.» 0.N 0.. - Id I d (093 Die Zl-O' “(and (m) as 90 OK 0.0 .soummssu How madman mo muflmsmv moawu mswamsou .2 N.95... 0d 0.¢ 0.m 0.N u - 0.. .m ousoqm ..0 «.0 mac . v.0 m.o ed 5.0 (cos bio 2|_OIX)("’).-! (")z" CHAPTER‘VII CONCLUSION The objective of this research was to examine, theo- retically, transition elements in order to determine whether or not their superconductivity is due to a phonon‘ mediated interaction between electrons. The task was to find an expression for electron-phonon coupling from an appropriate model and, then, with the theoretical apparatus used by Carbotte and Dynes13 for simple metals, find the superconducting properties of the transition metal in question. It may be recalled that all of the conduction electrons were assumed to lie in a single band similar to the free electron case. Although such simplifications cannot be expected to yield exact expressions, the results can, nevertheless, tell a great deal about the system at hand. First to be considered are the coupling times phonon density of states curves for tantalum and tungsten (Figures 8 and 9). Comparing these two figures to the corresponding phonon densities of states shows about twice the effective coupling at lower (transverse) frequencies compared to that at the highest (longitudinal) frequencies for both metals. in the case of tantalum, the curves compare favorably with 91 92 Table 6. 1 vs K for tantalum and tungsten. r Tantalum Tungsten Kf (x108) 1 Kf (x108) 1 1.601 32.8 0.75 0.405 2.0 124.34 0.8 0.342 2.1 5.414 0.9 0.266 3.3 2.358 1.0 0.213 5.5 2.207 1.4 0.183 4.0 1.933 1.779 0.164 5.0 l .041 2.0 0.154 6.0 0.671 Tc=4.48°K 0.69 Tc=0.012°K 0.29 93 experimental results obtained by L. Y. L. Shennl. Thus, it appears that the qualitative features of the theoretical screening are basically correct for tantalum. Moreover, treating all conduction electrons as members of a single band with no distinction for s-like or d-like states has not produced any serious difficulties. The curves (Figures 10 and 11) giving the dimension- less coupling parameter, 1, versus the Fermi momentum Kf’ require deeper interpretation. For tantalum it is apparent from the results that in the theory one must use a value of K1. about twice as large as would be expected for an elec- tron gas of similar density (1.6 x 108 qm-cm/sec) to achieve the correct values of 1 (0.69) and, hence, Tc (4.48OK). The situation with tungsten is somewhat reversed. Here the theory requires a value of K? of about one half that expected in an electron gas of similar density (1.8 x 108 qm-cm/sec) to get the correct 1 (0.29) and To (0.012OK). Thus, if the theory is essentially correct for group V transition metals like tantalum, then what pr0per- ties of group VI transition metals such as tungsten prove to be incongruent with these calculations? It must be remembered that all electrons were assumed to be capable of pairing and no distinction was made between s-like and d-like carriers. One such interpretation is that the conduction electrons in tungsten do separate into both an 100.0 10.0 1.0 0.1 0.01 1.0 0.3 0.1 0.01 94 F -l _ -'- ----- 1.0 2.0 3.0 ”.0 5.0 .0 Kr Figure 10. A vs Kf for tantalum. \ I 0.0 0.5 1.0 1.5 2.0 Kr _ Figure 11. A vs Kf for tungsten. 95 sélike and a d-like group. If only the more weakly coupled srlike electrons could form Cooper pairs then both the abnormally low value of Kf required in the theory and the low superconducting transition temperature can be Justi- fied. The assumption that the theory is correct for tanta- lum is supported, in addition to agreement with the Shen results, by reasoning that all of the conduction electrons can form Cooper pairs and that strong coupling results from the tendency of these electrons to form into a group of d-like electrons strongly correlated to the positions of the ions together with an itinerent s-like group. In conclusion, a theoretical investigation of the superconducting properties of transition metals using the same phononrmediated theory successful for simple metals hasfyielded credible results. APPENDICES APPENDIX A EVALUATION or I "1‘“ (1) 1a Ocase I?F%°6p= ZV—[Rnfl JOHN-fich- v- - :ZV—[mk W:D: y x (s h Pyjrzcjr T’Y‘ co C K‘- ’ :Z-Pvi; Dhofirgi gag?! X c (“H-ah 1% P )J S “(( Vz/Em 96 97 00 \+2 (PBPF:V—; {{Rno5;;1).'(¢+z Jx \__J:Z Fir—1. (2) 1= lease I 1('r'='¢»,.>==zv3—w\/;1_zg;'M-"‘ Wm» [RhéflJ'fiPHV-zdr 0 —ZV3~—WV1:—mPM(' cos agiDhix (*0 1° -r‘fl' 2.;1 co. - Eif-Me ‘5' SWAN-gr”? r V—T—rE—mcarfllfl him "" 1‘59»- :2\//37T(1‘)’ MN P(c.osep):DM" u-1)!{i+2./ “Cr:- ”a 1:ng 5! 221%;LZ)’ ]( 2/151) {W3 [(‘+Z)|(H3)I L J+(‘L)3( L212 P5’1 JZUPJ'JW‘T' ztu3- 2E“): (__P2. P2. ans] (3) 1=2caae hZM Meg—zviz—eaxl a / Rum ja( Fr) r201 " __ 51r(2-m! 'W“ — ‘ Z (awn): (COSGPZ—an ['3[:£11€“‘V:Zfi$‘ can,» 4* __ ‘- ‘ZA Efil _1 range ‘52 SLhFY‘J" P0 +‘F1f:€*zeo¢-ZCOSFPJ:‘| 99 (2%)) .. m. 30.1.1)“ 2) Hi LZéK-hz)! J“; 3&6 S‘jh’kflz); {E2 45259} «H! Elfiwafi It +617, _P_____“’1 [ii-3M3! [ng'fi—wos (2,”; m.) (~P?‘ —_2__/‘E )m ‘4“: 33a+z)l(.+3)!)g +(_() }62/4 5.27 P 3-2 ' “(PS-r3), J I 7! Z‘fi‘fi" {(P 2 _. Z/f‘nEth)t-r3 Ifggr Z. .5'____77(2‘”7)' [3%le 992R: i APPENDIX B EVALUATION OF INTEGRALS Part- I Exchange Screening [00:11) : 1+ km with EL” a_ function of V E— scalar k'only d \5 . Thus, 1181113 ‘ _ '3 dEk' "_ EU???) ENC) (“a S “93:5,Cb Cb and . I S Lkevr‘a‘b = Q3: Where due to the facet H( k" is small only near‘the Fermi surface ,. A k’ is the effective A316 _1__ . ’ E(l$'+cb)- E(K) ‘k‘ E ‘2 becomes 100 101 [Eds-2} T) HE) rEZAkScierJ¢ 1 443319232”) W lé erk'cose) :(ZflHeflEi/fj dkalecicishe mvm k2»: REZkk’c :6- over the range of integration lc' - k, (211) (9(5) _szkeda) :lfT—l—Vk 60(1— case) 7T _(_21)3I0(Ek)(h( (-1-C059)| —‘iabk V 99 3 2 243372 C(Ekh (1- C0560 \ 102 Finally, L L’ _L __ E(L ‘1) HL’) IL -L'\ '— ( WE $2;ka P) (h (1- cos ZPLSCHWc/ZKY) ) Where the correlation hole, Kb is given by 3° 9"- .S (n-1(Kc\Z Kr) 103 Part- II NOrmal Screening using cylindrical mrdinates take 3 axis 'It 2 /E(_L’- ‘jfi‘ ELL’) 00 and E(L)= 39m; . .::l+.(_ h) I‘ELFEJG9 5L2: CD ’62 (2 +ODJZ‘rFf-rf’z+z‘9 2,2)1/2 _ VMZ‘IT) £500? 27]”:ij zléflmei _i)_d_Z(l 87> (“(z'L— an) +415?“ -”z°ok,+°p) imp?) ) %% 16% =£2W{§"” 22W i} 4% (gm—L? (.1 2-H LKJ 105 mm A ¢Jc~uq. u iwuvmuag upsc»qrq_.xc~.c1_vpasacu «CH Ac.cmu¢vp¢§acu nwfim mm Ac.cmuq.pqzacu pcfim AUH .xcfiv chacu ~m- ~m Am.ouq\m.ou¢.x¢.~.mu.r_«.m_.x~Vpcyacu so” A\\\ m**A:z.\A:2Vc ac ><1a<1mm .xo_.fiI_vp<5acu ccfi o~ Am.ouc.n >Im.xm.m.our.u c:r.xm.d.¢u.u :Ir.xm.rpr.qp.xdvp_::u u_ :2c 2? > . .n_.> c2< Am_.m#.>u p q~ .Accc>.<»u uaq whzmzupchm zcszcu 2b m><33< Acuchmzuspcv >42C u:» c mfi u «A— X\23§§Cu —~ ¥._.~.XX.EX\(x\2CE§CU c~ 5x.pzbaaz.cm<.u2:4c>.u;c.¢Hu.<4\mx\2cs§CU o >.>u\dx\ZC5§CU a CC.J<§<.Tum?.QM.:.x3.hau.~Q\wx\ZC$ECU w 22<\mx\2c::cu «c .ccwyv_»u.A~Vzupca zcwvzufipc c Am~.m~.>u zchmzuzpc m Amfivmqsq zcbmzuzbc q .oqvxo.,m_.rcc.Am_vmcc.am~vMac.Am~v>..mvc.Amv;< chmzuyhc ¢¢~Vu4khb.accm_vupc.Accw_vho UCFUqu q u m23492~ Ch A:zvc zqaccoo Quipuzzua2 cuub>ua c o NHszmm¢ 106 pzucp ._m_. hzbza quc_ .cn_# Ccua flaphch: m . quq::2 chhcchuppzucp Ccua sq Cd qucz u csbcz ccmfincfpcz q¢ the :2< 4p; uc mchmzuzhc nib uQ< gzpcz C2< czpcz r¢ ~¢ Fmahu uung:cg czbbkcga uk<>phg< cq cm to >m cmpquchp .v»2u§:cc ac :— pun. :uua u><1 mu» ;m cuu¢4Qu2 Zuuq u>IO .xc~. bcfiocu vccm w AI. .¢.$~ur . v u Amv>zc .xc_v »($ch «com , Alp .¢.m_um . v u A~V>IO .xcpv pcyocu _ccn A uzckcu>zucpu :u. .xcfiv pcygcu cccm AIIH\ w4u22<19 >gzu:cuau *u>- .mpmv p<§ocu qu Anr~chc_.xc_.fi1_vp<:acu «mac ~m m\\ vquZZu chqvk<5acu m_~_ quc~. n Ax<2~ .2_§_ .xx .272 .qu7. Irrvpc§ocu -~_ qm A mp2;;<:u: .A22VC<1c czpp:qzcu 3cm c cpguc cc Am\c n~**c~ uc wk»2:. cugu xcs 2x Amxc mpa¢cp 2». zqzcc»m_r uc IFCHE 4u22<1g u <>c cc .zx .<>c .o__ szoa c: .zx .<>C .o__ Ccua czcg cpqc car Ccua . 2c zulka:u 2» :«ua a: CH m»2<»u2cu muzcu u: .c2 1 cuz A 44mg ppz:\w§c»q me .3: wugpb rn V chmzuspc 222 . com acu «u c2< ecu qcu mu .ua:p ugaspm acu _ n qua: cm: .FZZ .auuz .cm- bzpoa cu: .222 .auwz .cm- Ccua czqu «kc: CPA C0 «(0 CO mo #0 no No #0 c0 om mm XXXXXXXx xxxxxxxx xxxxxxxx xxxxxxxx (hp 465 c.~ n no xxam n xcyhfi cmr CF cc c.~ n on yx¢¢ H 5.12:, car C» cc w.c u an Xv. H XCTP‘I aumz.Aq_r.r~r.m~rv CH cc c.». H AEVFYPU cs~cz.~ n 5 ¢c¢ CC c.c n Ahvahu c.cuh~vhko €5»CZ.~ H ~ CF CC pmzcc ;.co qmturqc¢¢._¢muc2ZCU ufcv Czpu oquZZ uczu:cum FZFCQ Ian: uI» ud C» v» c u wz mcuau *cm2* uc mmau 2» «Pg uC pm .uc. 2F A:Zvc CZWQCHM hochv c..w + <>c\zx u ¢4 uxu.~_m thaa au.fik.uufi.,hvzu»cav.-_m pzpaa A¢.Pn~.fihvzuhcav.~>~m thoa Acuz ._n~ .A~.xa. .~c_n pzhoa Aacflm .~¢V ukpaz .c_.m<:< .A>vvq.w¢mv Ch cc mm. 5.~ u x cmr cc «(a ¢r~ x<:fi.fiufi cwr cc h¢m .Hulw.c+xxu.mq ~n~ c.m\u::u n u4<1< _H U.~H§x cc CC o_~ L m- Ar.— u zyv Upaca uuolp CHEH a: zuycqa 3c: uzck 2~:C44paa Q ~.— u ~c~ pzhaa xqsh .2—2» .yy .222 .unZ ._p.~ bzhaa ~c¢ ¢~_ r n zpzw yum mfl_ ~c¢.~mm.~mmAzthvup cmmm «#P xqufi u xczp cum" -_ cmrm.cmr~.cmrfiax<2b.u~ ydi—7X\&Qnuf:o chFTchfpfix cmm 110 cc- cc— NCH cc~ mcm «£— me— «(a cc~ om~ mm~ cm~ mmg cm— mm~ fimfi cm~ o¢~ «q. b¢~ (a. m¢~ ¢¢_ mq_ m¢~ fic— <;.~HH cm CC mc.>¢.ch.lzbzhvu~ ~¢.~¢m.n¢mfiuaq¢ $9 qu U \l L cwm Chm 9 UV OK on «Tc ma .Amcchcou Ic_I3 chpzugupqpm uzp uc 51cm zup .qcpcc;¢ ulh c A 222 .~ u 2 .Ar.z.>uv .wccm czhaa A 222 .~ n 5 .Am.zc>uc .mccm thaa . 222 .~ u s .,_.2.>u. .ficcm szaa cccm szaa .222.~u2..:.>V.Ar.#ns.Asccc.:.x.1.b.w2 .mcfi pzpaa on.¢r.¢rfifiicycup Mrcc.mcc._cc.cccuau 4;¢m (cm *U aw .AmvBCJqu IcpIB vabzufupdhm ulb uC Ezcu ¢u— 4~ ~b~ cfifi wh— m>~ («by m>~ ~h— ehfi UH cw .AmV3CJ4cu chIs .Amvchgcu IcpIz cmcpzuzupqpm c1» uc :ocu ¢u_ .om_m Hypaa rcc~ bzpaa Acf o~nfioasvdhcv u:2~»zcc d u szh cc ck cc ~ n FPEF zy.aoc.c>.oc.oc.-c c» cc m¢1~u2~5~c x<2H c.m\uq:c n uc:u XX$N H xx ficp szaa .ZHEH .mz .xx .5! .m~fi~ hzbaa aQ\A1c Auaqhzcuh cq.c¢m.cqm quczcuh Ac.»c.ua cqfi cc (cm (cc c¢C$U Amvhzufiuhqhm uI» uC Ezcu zuh Aquhccgt uIH ccc.c¢>.c¢(1~u2_:»cu_ u:zphzcc u:2cp7cc >c.>c~.bcfiAuQc.c .cm~ Ic23a ughph .rm~ IUZZQ some .flccc .oomc Alcozc up u:z~p2cc .n.¢:.~.q>c.<2pcz.hhccha3 44cc mcc_ pzpaa caux.».Hpcaacnuc.ru.pcphu , <>c+couucau <2.#n» cncm cc c.m\<>CIncau A~.<7.~.<>C.<§~CZ.QFUV HQB 44cc rmfim pzpoa d>c .mc~> pzpaa uthp .rmpo szaa AccH thaa momr.comr.ccnc5<»c7cup c.c u Aficho c.cuAfic»»c <4 .<< c» mc cm mFqumzco ucacu mupga»»4:z FUZCC uZPHICQUZM czu op.cm.o szaa cubicaxc¢mm¢.~\c+ccuucau ..c<»c+ucccuucuc , <2.~u_ qcm cc c.cuucuc c.cn E<4u c.N\<>Clneau A—Iooc..cv hcga Azcc >c chum uzpp:cac:w uhzh ccm comm —~om UUUUU . 11.3”.» *9 114 mhm cbm mwm NNN Fhm Chm 0cm «cm 5cm rem Ucm ccm mcm «cm ch com 0mm cwm mmm qwm mmm me 0mm ocm «cm >¢m {cm mcm vammc cmcccc cth zux<3a< .uch qcnwc uIF 2F FFccv<2cF>ch Fcha¢c.chF.xa c~.c~nF FFcF cc Fcha¢c.anchc «n.01HF chF cc Fchaac.anvaa cc.wan och cc Fcha¢c.¢nFcha cc.m¢uF mch cc Fchatc.ancha Fc.cuuF wch cc Fchotc.ancha Um.mmuF cch cc Fcha¢c.chcha Fm.m¢nF mch cc vax0$comnFvaQ F¢.cruF cch cc Acha¢c.¢nFcha nm.¢ruF mch cc Fcha¢c.ancha cw.crnF mch CC ,chaac.¢uAcha Fch cc WU LcI FFcF chF 0coF qch bcoF (ch mch ccoF rcoF Ach Fch L L New Fem ccr 00m 116 mom can can can Nam quuc.Fnch dch cF.qch.cFFc.FI:ch c.cncqu co 23:»u1 FFVFcanFcrcc Fm~ FFchchchcc .r\CrC\FFF.11.FFV‘:FVHFFV Fcc r.Fu_ Fm» cc chc>nFrcc> ch\FFm.FchlFF.F.c chxFAm.Fc¢q*c.m+1F.chqchgcunAFcc> Fq¢._cczucu Jcchc> chxFFm.chuFF.Fc FUM\FFA.Fc¢¢c.m+FF.FV F<.cczucu cccc cc cw cc Fchc.uz.FFccch co cF cc FFrcc.u2.FFccch Z\hx\ZC§§CU >.>u\¢x\?csscc cc.w<5<.qum2.cv.:.xo.bac.Fa\rx\zc:fcc 22<\mx\2C2:CL FmF .ch mirc chmZuSFc Fch v<§< chmzuch FmF.4>.Fanmcc.Fvamcc.FnFcFcc chuzcch Frvqc.Fcccxa.Fvaz.FmF.>.FmF.chu chmzuch FmF.mF.<.Fchu.chc chmzuch Frccomcc.Fcc.Cvcuxu uZthcac:u czm ZQEFLQ c0 117 mmm cmr rmm mmr cmm 0mm car bmm can mar can mam me cam oFm cFr hFr cFm dFm FFr cFm 0cm dcm New (cm mcm Uficm rcr chFcchccncc :.Fuw mm cc cm.mm.¢mFchc:chcch chc.¢cccccc gcdc FF.cc 2.FnF (F cc Fccc.mF.F.r.u.>cc zucFu 44cc F1.Fcc 2.Fnc mF cc _4.FnF «F cc F<.cczu4ucgcc c.Fuc.>u\¢x\2C§ECL cc.u42<.Icm:.cw.:.x1.Fcc.Faxcxxzcsscc Fch> .FmF.ch>c chvzuch FmFVVnFFc> z.an chm cc Frcc.vxa FMF¥.rFL + FT.FVLV$FF.UVCIc.ncaFoccx0 FFFfioccczFx.mvc + Fx.¢cL$Fc.m.Lc¢FF.rcClc.mcaFcc.x0 FFx.rcc$F5.mcL¢FF.dFCIcoFctchch FF¥.¢cc$Fc.quIc.Fc¢chcxQ FFF¥.¢VU + Ffi.d.cc¢FF.¢VCIc.mc¢chcxa “cc FFFX.UFU$FT.NFU + FxomcL*F7.JFLV$FF.FFUIC.RV¢FGUVXO FFFx.mcc¢Ffi.FFL + Fx.FcL*Fc.vav$FF.mcclc.mc¢F>mch FFFxomcU*Ffi.FcL + Fx.FCL$FT.KVLC¢FF.UVCIc.mc¢chcxa , FF¥.UFU¢F7.FVCIFX.FVC*FT.wcLIc.AV¢F¢UVx3 FFF¥.UVL + Ffi.mch*FF.FcLIc.nc¢FcUcxa WAFx._.c+1w.cccc$1c.uccuc.mc¢1mgcxa + FFF¥.¢FU¢F5.WVU + Fx.mcu$Fc.chv¢FF.FVUIc.mc$Fc¢cxa FFFx.¢VU*Fc.FcL + Fy.Fvc¢FT.¢VLV$FF.mVLIc.mc¢Fw¢cxa FFFx.va¢FT.F.o + Fx.FVL$F7.r.Lc¢FF.¢cCIc.mc$chcx0 FFFx.NFO#:,o+:.L + FxcdetFF,.mFU.$FFonVL I ComctFF+le + FFxova$FficmVL*FFoqu I con$Fm¢va HmCC< CC<+ c<+ FF.FVcxo 11:.mchAPIVCm I Fs.rcwI1FImcchAn~cxo ..c.¢.mIAFI¢ca I AsIVCUIFF.qcmcaFXIrccacxccxa AvIrcuIAF.«chcx.qchA~ccx0 AwI¢cvacFI¢cu¢F¢ccxc urc FAwImcuthIACV I .:.ncu¢,F.mcuc*1x.Fch.Fccxo .17ImchAFIch I 1:.FCVICcImcucIAxImccmccccxa F.5Imcu¢.F.ch I Ac.FVVIAFIKCUVI1x.m.cIAomcxa FFT.UVU#FF.FFV + F7.FLU*FFIWFMC$FUEVxQ F.1I¢cm¢,cIm.u I AwImcm¢.F.qcmcI.y.FchcFucxa AAFIF.m*.c.¢cu I Aw.FcuIAFIJVUCIAXImchAcmcxa FFfi.NVU$FF.FFU + Ffi.FvvtFFochcanochzFocvx0 AAsImcvacFI¢cu I F:.¢cm¢AF.mcmcahx.mchAchxa I Fa.nca¢.FIncVIAXI¢chAqqcya umcc< - cc< I FFIfivanFIccq Ffi.mcv * FFIFFU * F1.NVL ¢ FFcch + U FFfi.m.v*FF.mcv + FcIVCm¢FFIKVUV t FXINVL t chcxa + q FF1.¢VW*FFIAVU + FcIAthFFIGVUV t Ffirvxa + r a F I—'(\,("\1 < + +-++ +<++ CI++ + + HRFQU'ifiq FFfionwsFFoqvv + Ffi.¢.M$FF.F.UV * FXIFVU ¢ Frrcxa + FfoFcu * FF.F.m¢ Fx.qco t Fmrvxa + Fc.«.m ¢ FFIrcm¢ Foacxcuccq FFT.chc ¢FF.N.m+FF.mvuzFF.Fvwv$FxIrcL$mech +FFT.FVU> tFF.FFU+F5.FFV#FF.FFUF$F¥.KFL¢FFRFXQ +FF7.FFU( #MF.NFM+FfiomvmaFFIMFvc$FXIF.L*decxa U +FXIRVU$Ffi.mVU¢FFINVU¢FGFVxa q +FF3.VVV#FF.FFV+F1Iva$FForcvc*Fchxo +Ffi.mcmr #FF.mvv*FmFFxQ +FXIRVLaFfi.chtFFIFvchIan & +FF1.mcm*FF.FcU+F5.ch¢FFomcvc¢FXIFVL¢F Fovxa +FfioFFUtFF.va¢FFVxQ HFFoivq c mccc I accq I FFIFccxo I 1F:Imcc¢.x.m.c I .XImchA:.«chc.mc¢AF~cxa I \Cf‘dO‘ 123 cxc Ohd chd Nfic {Nd mhc cbq mwc mbd th ehc occ mcq bcq ccc mcq #cq m¢¢ mcq ch ccc omc cmq hmc (mq rmc mmq ch cmd o¢¢ FAAFy.q.cI,_.mchAx.mchcfi.ccuc¢Ammcchc AFx.m.cI1:.mcc.¢.F.qcc¢Acmcxccac.nI1AXI¢CcIAqucccq I,AF.«.cIAFacxaIc.mI.chxacIAFIchIchcxacIc.nI~ A.FXIFchFTIFCcIAXImchAs.FcchAFImchANFanc +Fxoch¢Ffiorcu¢FFoFVL$FIFVxQV¢CIQIU AAA¥.mch,w.mcoc¢A¢cchAFIKCcIAmcxac¢c.mIFAx.mccq IFsImccaficcxaIFAyIACcIFw.m.cc¢.FImcc¢Amcxchc.¢Ir AFAXIFCthfi.FchAYIFchFcIFCchAFIFchFocxam IAXINVcIFF.mchcF.acc¢.nchaIAx.FchA:.Fch IAAFIFCcI.q.onFFIFch.Fcx0..ICImIFaanFIFcc ch FIFux mF cF c» cc rnx cF mFIqFImFFFIFch mF FIFufi «F cF cF cc Fun FF RF.FF.mFFrIFch FF cIFnF cm cc m$*FFoRFU*C.mIcoFHFFo¢.U CF FchcL*FF.m.U¢c.&HFF.¢vm FchFv$FFoFCMIFFIFthFFIAVLHFF.rFL FFIAVC¢FFIva+FFIFcL¢FFIAVUHFFIFVM N*$FFonw$c.mIcoFuFF.ncc FF.FVC¢FF.FFv¢c.AHFF.nFm FFFVC$FQVUCLHFF.FFU FFFVC*FQVZFUHFF.FVU NIF H F cF CC m LLa U U 23:»ua «FF CZZFFPCU cmF FF.ficx\2cyscc NXIKcINFchISxxcx\zcsycc EXIFZFQQZIcwq.u534c>IuccII>u\cx\7c:scc cc.m<5<.ccmz.cv.:.xQIFacIchmx\zcy§cu chchchFchCVQuu3v c7FF3cocrm czc ZQChua FFc>\22q¢> n FFcc FF.yc>u¢Fx.Fch>u> zIFnychc qucIFcc Ffi.xc<¢FF.xc>c+xux ZIFumecc c.cnx ZIanFFcc cIcn> ZIFchFcc Z\>x\2cy§cc >I>u\cx\ZCS5Cc ZZ<\AX\ZCSECL Fch>IAmFInFVQIFchcIFmFIAch.AmFIch>c chuzu§Fc FcI cucc x Icmchzccu FcF .mImFum IcIFch .xm.F<§ccc chchcc.F«.2:.FcchcFFc chmzuEFc FmFIRFc>c chmzuch Fnchqzc chwchFc .rchACccxc chmzuch Fchmcc..chmcc.,chFcc..ch>..mc4<.1ccm~.cn24cc Fm$$FJFFCC+m$$chmCC+m**F4vFCCFFQCUHCQC . cm cF cc Fcc> .mch Fchc oFIcFIFanxIF4c>ch 5mm cmm mmm mmm me Fm OF L C 127 omc cmc ch cwc ccc ccc mwc mwc ch ccc ccc ccc hcc ccc mcc ccc ccc mcc Fcc ccc Omc hmc cmc mmc cmc mmc mmc ch cmc cmc cmc hmc cmc mmc mm H 1 cc n 4 CNN C» Cc RBI H x3 cIc n c mcc u c mcq u 4 CNN CF cc coc\F<>c¢«2< I x<§cc u 3c . a: u zez ow + «22 n «2 F I m2: H n2< CIF I <>C\x<§c u «22 ccc¢23 + Fcc> u x.rmhommc.thc CH CC NIF H 551 ccc CC c.o u 5cm FFV;<\CU<¢:$c.cHFcc FCUE I cmgcc*2243 H cc VCI n U0 F3*2CU H mm Acmx I rm4<¢c.mc¢zcc n mq Fch4<+Fchqc¢uch<¢mcc¢CIN n «cc . Arcc<¢cc<1;cflv I FcvaLIZu I c¢AZ< I x<2uc M 2m 21 I 34cc u 3c omIAzz n A2 F I A22 I A24 c.F I <>cxx<2c I A22 A2 u F2 czc¢x3 I FCC; n xczc CUI UAc ccm CqN 130 GUN UmN cUN MUN AUN FUN GUN ocN ch NcN ch UcN ch ch FcN ch OmN qu NfiN UmN cMN mMN ANN FMN CNN OAN cAN «NAN NAN cm222vIFLwaaccuuuFC FA: .Frux .ch>qaccc .cxu .AcF FFFQQ m cF cc Fm: .anx IFxc>cocqccq 2ch2u5Fc FFEuchZI4z.chCISFCZI> Fcucqcc CIF VF can: chZFQQ uzrgc> uIF JLCZIUUFCIUZDJC>IF3.3cchmrcAc.FAcqcoFFVCQIC0c .F.C> ICCCF FZFQQ 22m I c2:cc> uch F I 2F22 I A: Jucz :\2:m n 52¢ FZFQQZIFcrICNC CF cc x? H 23 2.cc I 2:u n 22m FFZCQFC I 24cc n Fcmch FcmFFcI2ccCIccuuFF2cFFu FFN*¢x3Ix3*23IA*¢23v*U I Fchzsccq I IFoFcc IFoFcccca:m.FcAccaFcIFCAc:4c>.FFIU c2uch c2q uu:c<>2uch zucFu CUCC cc Fc.I 23C2auFFc:_<> FFIchuFFCc2Fc 2IFnF ccA cc Ucz:cc c24<>2mch cNchFFF2F c2< anczucu rocc F3CCc2 u:2FF2cc AFuFFIFIXCI FFIFIXCII2FUIF_.xchUCLHFc.xcc FTIXCT¢ZFUIFFIFIXFd$UCUHAF ZIan cUF cc cNF .ch .CFF FEV UF AFqu.FIFc< FXIFIFF<*ZFUIFXIFVd$JCCHFxIFyd FXITF<$ZFUIFXIFIFF<$UCLHAN ZIFux UAF CC AFuFFIFchq FFIFIXFchZFUI: .IxccafldCUHFchcF‘ Fficyv<¢ZFUIFFIFIXC<*UCCHAN 2IFux UcF cc CCO CF CC FCIFccuAF FFIFIFCzuch 7C UCTCCa acFIcFF IUF4cIFccc\FF.:._<>VCUTF: uF rmc .Umc .rAc FFFIFC:44>C uF AAc IcrN .crN FNIchUIFFc:Cc>IFFIFc:_<>c uF FAc .UAc IUAc FFIFC cF ZIFnF cUc cc c.anc.FC< 2.Fuc ch cc 2.FuF ch cc AuFFch ij Icoc IFOU FEC uF chF2ucFu FcCC<2 ccA CF cc c.A\FAF.:;<>IFFc4;c>cn:\FFCcc<>cucccch ccc .ccu Iccm AFFc:cc>ch cUU IcNU .cNU FcIcc.UIFFccc<>IFFc24<>ch u22FFZCC :anch:c<> crU IcrU .UAU Fch:;<>I: ccc .crU .ch FFxc44<>Iquch ccc CNc ccc CUc Ccc Ufc ccc «Ac UAc MAc AAc FAc cAc UFc CFc CCc Uom Com CoU qu Cam UNU CNm ccm CUm Ccm crm UAm CAU Coc ccc CNc 135 NNFCchF ANFCCFuF FNFCCFUF CNFCCFUF ochoFuF choCFuF NcFCCFUF choCFUF mchcFuF chcoFuF mcFoCFUF mchchF FcFCCFcF chCCFUF OUFCCFuF «UFoonF NUFooFmF cUFCCFmF UUFCCFLF cUFCCFcF mchchF NUFCCFUF FUFocFmF cUFcoFuF ochchF chCCFuF NcFocFUF chcoFmF UcFCCFuF chcoFmF mcFCCFmF AcFoCFmF FcFCCFUF ochoFuF ochoFUF cmFCCFcF ZCFFqFCQ UC CCCZ< u: szmCC FA$$AFIA¢$FFCUFQCUHN ch .cco .cFo FAFC uF c2< czFU chcccch 7&ZNUQ FFIFchFXVFIIFXIFFCHFFIIFvc 2Ian Uqc Cc CICHFcIch 2.FnF Uca cc FCchanyc: 2IFux Ucc cc 2IFnc Ucc cc chqu acFuc> Fccu<fi IFF3 chF<2 2CFF\> FIC .UFxc x m22FFZCC FFvunxczc AU. FUIFU .ZFCCIZICC I: CSIFHF AU CC c.cnxc2C Cf\Co—HC CZ¢CSHCZ FIc2cIch F<20cc FCZVC cc FCJQIFFIXCF.\\\\C F<23CA FCFFCZVC ZCFUzuch FZILZFCFILC Fcau uZFFrcqc2U czu FquF .FchIchc cF cc FFch .FcUc.ccc cF cc FxAFn2Fm F\FFHUCC AU FU cc Fc CF rCF ACF FCF ch CUT ccc CFO UAO CAO 137 FNCFCcuIFFCZcho 2CFU2u2Fc FclacqInc.zchc.uchCIo.A.c Fccc 42 xamFFC;A caFux CEIFuF c cc xchFccc 1F..c..cc Fccc 41cc xchcha FFIIc.>. Fcca chc xchFCJc FFIFIcIFC FCCc 4cc Fcca 422C U>I>u> x<2IFuF N cc xcmFFccc .cu> xawFFC4c FFIICIICC Fcca 44 .cchcccuchcFIIFchumx :FchIc.cchncxm FIcFuI$c.cccAua>U GIchIcuIFch c.a\c.UnFIcFuI FIFx<2ch2anq2 FUIUUIUU FCIcFIxFFF2 :c .zchpjxpccmpc cum—445aCZcch. b\> _.c u xoqz _ .uaxq x 2c capacia >czuccuoaczoqv pcaacu mcfl xamhhcqa A¢.m~u .n be ac u:;<>clcfi. H>~ccwu~ paxup .Aofio.cc~c c» cc coo cwficchmh haxu» .Acwc.ccv c» cc cro m>_cc»mh mepn7aw umo qhfioo_uw p\~»nvcu cmo 137 xaubbc4a xamFPC4a xamFPC4d xawhh04a xmmth4a xam»»C40 xaaPPC40 xawhp040 xamPFC4a xamPHC4a xamkbC4d xamPPC40 xam44040 xamFPC4a xawppo40 xauhPC4Q cupp04o xauFPC4a xamPHC40 Ancpcqu.1~vzukc0 zcamzuzpc 14104<.du.zuhca.uxu.c:0.N.ecmgc.wm<:c.cv Ian uzaercoc:m czu zereua Aaxm.0>m.~.x.>c pc4a 44cc 1.xcsceqc4ueicyceqc.u*c.cac\aavc¢c.au> Cit—"X cz.Hn_ m cc AH..1c..c4 22240 .e_. Fc4a 44qc a~.~.c.>v PC41 44v 4:40 44+>u> x<5.1n» s cc .cn> 11..c..cv Fc40 44cc 1_.x..c. ec40 44cc A_.x._.c. FC40 44 .Cthw o.uxc.mureca3 c.a\c.unelcau1 _+Axhpkzd:c mIH uc hacm muEFP CZH4QZCU uaqm czpu u IQXaluapa¢~roc.cnlcxu czwzuuaum uCZu1>ceccu w c.cmpu1@.ec