ON THE THEORY OF RADIATIVE BETA PROCESSES FOREST GENE HESS A thesis submitted in partial fulfilment of the requirements for the„ degree of Doctor of Philosophy in Physics We accept this thesis as conforming to the standard required from candidates for the degree of Doctor of Philosophy Members of the Department of Physics The University of British Columbia June, 1955 ABSTRACT This thesis is devoted to some aspects of the theory of the weak continuous gamma radiation (often called 'Internal Bremsstrahlung') which accompanies beta processes, i.e. negative and positive electron emission and orbital electron capture. Whenever a beta process is accompanied by this gamma radiation, i t will be called a "radiative" beta process, otherwise "radiationless". The results presented in this thesis go beyond those obtained by other authors in two respects. In the first place, the radiative beta emission probability is calculated for an allowed transition taking into account an arbitrary mixture of all the five beta interactionsc Previously, only the theory for the case of pure interactions has been carried out. In the calculations, as in previous ones, Coulomb effects have been neglected. In the second place, the radiative K capture probability is calculated for an all owed transition taking into account again an arbitrary mixture of the five beta interactions, and, in addition, Coulomb effects. Previously, only the case of pure interactions with the neglection of Coulomb effects has been considered. In order to take Coulomb effects into account, a "semi-relativistic" approximation for the solutions to the Dirac equation with a Coulomb potential has been developed. It turns out that taking Coulomb effects into account reduces the probability of radiative K capture by an order of magnitude. THE UNIVERSITY OF BRITISH COLUMBIA Faculty of Graduate Studies PROGRAMME OF THE FINAL ORAL EXAMINATION FOR THE DEGREE OF DOCTOR OF PHILOSOPHY . FOREST GENE HESS B».A„. (British Columbia) 1949 M.A, (British Columbia) 1951 MONDAY, June 27th, 1955 at 3:00 p„nu i n ROOM 303, PHYSICS BUILDING COMMITTEE IN CHARGE: HoF.- Angus, Chairman External Examiner - L„EvHv Trainor of W<, Opechowski G.M.- Shrum J,B. Warren F„A. Kaempffer B,Co Binning D. Derry T„E. Hull A.Do Moore Queen's University LIST OF PUBLICATIONS Some Directional Correlation Functions for Successive Nuclear Radiations, Canadian Journal of Physics 30, 130, 1952'. Alternative Solution to the Ehrenfest Problem, American Mathematical Monthly, 6JL, 323, 1954« THESIS ON THE THEORY OF RADIATIVE BETA PROCESSES This thesis i s devoted to some aspects of the. theory of the weak continuous gamma radiation (often called "Internal Bremsstrahlung") which . accompanies beta processes:,, i.e. negative and positive electron emission and orbital electron capture. Whenever a beta process i s accompanied by this gamma radiation, i t w i l l be called a "radiative" beta process, otherwise "radiationless"'.-The results presented i n this thesis go beyond those obtained by other authors i n two respects. In the f i r s t place, the radiative beta emission probability i s calculated for an allowed transition taking into account an arbitrary mixture of a l l the five beta interactions. Previously, only the theory for the case of pure interactions has been carried out.- In the calcula-tions, as i n previous ones1, Coulomb effects have been neglected. In the second place, the radiative K capture probability i s calculated for an allowed transition taking into account again an arbitrary mixture of the five beta interactions, and,. in addition, Coulomb effects. Previously, only the case of pure interactions with the neglection of Coulomb effects, has been considered.. In order to take Coulomb effects into account, a "semi-r e l a t i v i s t i c " approximation for the solutions to the Dirac equation with a Coulomb potential has been developed. It turns out that HTaking Coulomb effects into account reduces the probability of radiative K capture by an order of magnitude. GRADUATE STUDIES Field of Study: Physics Electromagnetic Theory Theory of Measurement Nuclear Physics Quantum Mechanics Spectroscopy Special Relativity Theory of Solids Group Theory Low Temperature Physics Quantum Theory of Radiation Theoretical Nuclear Physics Theoretical Physics Seminar Other Studies: W. Opechowski A.M. Crooker K.C. 'Mann G. M„. Volkoff A.M. Crooker W. Opechowski R.W.- Maurer H, Koppe J.M. Daniels F.A. Kaempffer W. Opechowski W. Opechowski Network Theory Integral Equations Complex Variable Theory A.D. Moore T.E.- Hull D'„. Derry ACKNCWLEDGEMENTS I am grateful to Professor VL* Opechowski for suggesting this problem and for his continued interest and.valuable advice throughout the performance of the research. I am also indebted to the National Research Council of Canada for financial assistance in the form of a Studentship and Fellowshipo Table of Content8 Chapters Paffe Introduction and Sumoaiy .............. ..»* 1 I Outline of Theory for allowed Transitions 8 A. Non-Radiative Processes 8 1. Electron Iiknission 8 2. Positron Emission and Orbital Electron Capture 13 B. Radiative Processes 13 1. Radiative Electron Emission 14 2. Radiative Positron Emission 17 3. Radiative K Capture 18 II Radiative Electron and Positron Sales ion 23 III Radiative K Capture.*. 31 A. Harrison and Sohiff results - 31 B. Radiative and Radiationless K Capture 35 C. Expression for L 44 D. Calculation of L in a "Send-Relativistio" Approximation 49 E. The ease of Z*119, A- \ 68 Tables 111.1 D+ of Eq. (in.D«15a) and Q+. of Eqs. (III.D,12b) and (III.D.Mb) . 7 2 131 111.2 Integration results for Cs 74 Table of Contents (Continued) Graphs Facing Page 32 ' • II.l Qlt Qg of Eq. (11.16) for P 29 HI.l Integrand of Eq. (III.D.15b) for Cs 1 3 1, Ey=102 kev 63 131 111.2 Integrand of Eq. (III.D.15b) for Cs , E y = 154 kev 63, 131 111.3 Integrand of Eq. (III.D.15b) for 0s , E ? - 205 kev • .63 131 111.4 Integrand of Eq. (III.D. 15b) for Cs , E f _ - 308 kev 63 i n .5 Comparison of the results of Saraf , Morrison, and 67 Schiff, and Eq. (III.D.27) III.6 Comparison of nnon-relativistic M, nsend.-relativistic" and "exact" integrands far Z= 119, A- £ 68 Appendices Page A Evaluation of Traces 75 B Table of Functions 77 Evaluation of B ^ i n (III.C.16) 80 Rietoann Surface and Hyper geometric Fli. 81 C Normalized Solutions to the Dirac Eq. vdth a Coulomb potential 83 ffave Functions Used in the Thesis 87 Normalization Procedures 88 Comparison with Other Solutions 92 The Semi-Relativistic Approximation 94 D D e t a i l s oP Ck. LI 1°* References * 108 I INTRODUCTION AMD SUMMARY Every beta process, i.e. electron or positron emission or orbital electron capture, has a certain probability of being accompanied by a weak continuous gamma radiation* This radiation is not due to de-excitation processes in the nucleus before or after the beta process since such pro-cesses give rise to monochromatic radiation. Rather, i t is attributed to that caused by the changing dipole moment of the atom when the electronic charge, involved in the beta process, is suddenly shifted; i t occurs "during" the beta processo Such radiation has been termed "internal bremsstrahlung". Whenever a beta process is accompanied by internal bremsstrahlung, i t will here be referred to as.a "radiative" beta process, otherwise - "radiationless". The measurement of internal bremsstrahlung is complicated by the fact that i f any deexcitation processes occur, the resulting monochromatic radiation will probably "smother" i t . In this connection, internal bremsstrahlung was first observed by Aston.(27) in his measurements on the electron emission of fta ET . The bremsstrahlung from electron capture was first observed by Bradt, et al (46) in FeSS. Let us first consider radiative electron and positron emission. The theory for these processes has been given by Knipp and Uhlenbeck (36) and Block(36) for an allowed transition and for "pure" interactions only. In what follows, "pure" interaction means any one of the five well-known beta interactions. Chang and Falkoff (49) and Bolgiano, Madansky and Rasetti (53) extended the calculations to first and second forbidden beta transitions. In all these calculations Coulomb effects have been neglected. Goodrich and Payne (54) measured the internal bremsstrahlung r>3<2-aecompanying the electron emission of r « They found that their results agreed quite well with the theoretical results for an allowed transition mentioned above. The calculations given in Chapters I and II of this thesis take into account an arbitrary mixture of the five beta interactions and the results show that a non-zero contribution to the probability arising from cross-terms between scalar and vector, tensor and axial interactions, is present. The calculations are carried out for an allowed transition with the neglectibn of Coulomb effects for electron and positron emission. The results' show that the pure interaction terms of the probability are the same for electron and positron emission; on the other hand, the cross-term contribution changes sign from one case to the other. The explicit expression for the photon energy spectrum is given by Eq. (II.16) and the corresponding curve's plotted in Graph II.l Because of the agreement of Goodrich and Payne's results with the theory including only pure interactions, i t was felt that an upper limit on the magniture of the interaction constants apoearing in the cross-terms could be obtained from the condition that the contribution of these terms does not destroy the agreement. However, after the calculations were carried out, more accurate experimental results appeared in the paper by Liden and Starfelt (55), who criticized Goodrich and Payne's result and concluded that there is a considerable discrepancy between theory and experiment. Consequently, the program connected with the upper limit on the cross-terms had to be abandoned because a proper comparison could not be made until Coulomb effects were taken into account. In any case, the results do extend previous calculations and should be of some 3 use for comparison with calculations including Coulomb effects for small values of Z. The calculations of the radiative K capture probability were carried out for the following reasons. Prior to the experimental work of Saraf (54)» experimental results agreed with the theoretical results of Morrison and Schiff (40) for an allowed transition and pure interactions neglecting, among other things, Coulomb effects* These authors calculated the ratio v/^ i k j ^ as a fn. of the photon energy, where vv^A is the probability per second that a photon of energy t i c k in the range dOicPi) is emitted during the K capture process and we is the probability that a radiationless K capture occurs per second. Before Saraf«s work this ratio was measured on an arbitrary scale as a fn. of -4 the photon energy,, Since' the ratio is of the order 10 , one can say that Wc is practically the same as the total probability for K capture. Saraf, making use of this fact, was able to measure the actual value of this ratio as a fn. 131 of the photon energy. He found that for radiative K capture in g -^Cs agreement existed between the Morrison and Schiff theory and experiment for high photdn energies but that definite disagreement outside of experimental error existed for low photon energies, the disagreement increasing the lower the photon energy. This disagreement may be due to the neglection of Coulomb effects in Morrison and Schiff's theory or possibly due to the neglection of the radiative \> capture from the L shell which may be of considerable magnitude for the second order radiative process* The theoretical program which was started to explain this discrepancy was to calculate first the radiative K capture for an allowed transition taking into account Coulomb effects and an arbitrary mixture of the five beta interactionsj then, to do the same for radiative p capture from an L shell (the results for the radiative S capture from the L shell are expected to be much smaller than those for radiative s capture from the K shell). t Because of the complexities involved in taking Coulomb effects into account, this program has been carried out only for radiative K capture. It turns out that taking Coulomb effects into account reduces the probability by an order of magnitude and hence increases the disagreement mentioned above. It follows that the radiative /> capture probability should be investigated. In this connection, i t should be mentioned that Glauber and Martin (54) have published a short note in which they state that they have calculated radiative S and |» capture probabilities for an allowed transition taking Coulomb effects into account. They obtain explicit results for he which was studied experimentally by Madansky and Rasetti (54) who also found a disagreement at low photon energies with Morrison and Schiff's theory. Glauber and Martin find that taking radiative capture into account explains this low energy discrepancy. • However, the information in their note is not sufficient to compare their method (which is entirely different from ours) and their results with those presented in this thesis. An outline of the contents of this thesis on radiative K capture is given below. In Chapter I, the probability for radiative and radiationless K capture is formulated by means of the Dirac hole theory along with a discussion of the intermediate states involved. The contents of Chapter II have already been mentioned above in connection with radiative electron and positron emission0 In Chapter III, Section A the results of Morrison and Schiff are obtained in detail. These authors only outlined their calculations verbally without showing explicit details; for this reason, Section A should be useful i f one desires a detailed understanding of the assumption they make. Essentially, they have neglected Coulomb effects entirely, representing the electrons and neutrinos with Dirac plane waves and, in particular, the K electron by a Dirac plane wave with momentum zero and energy fnc* i.e. a particle at rest. In Chapter III, Section B, the radiationless K capture probability Wc is calculated for an allowed transition taking into account Coulomb effects and an arbitrary mixture of the five beta interactions. It has been obtained already by other authors but is given i n this thesis as a preparation for the more complicated calculation of the radiative K capture probability vv^ J/i . In obtaining. Wc the well-known assumption for the electron and neutrino wave fns. at the nucleus: where is the electron or B B u t r i n o momentum, R the nuclear radius, is made. Also, only.the components of the wave fns. having an angular dependence given by the zero order spherical harmonic are retained. Both these approximations are common to the allowed transition in the radiationless case. Also in Chapter III, Section B, the radiative K capture probability W^Jk for an allowed transition with Coulomb effects, is reduced to a 'sum' over all the energies possible for an electron in the intermediate state. In obtaining the wave fns, for the electron and neutrino at the nucleus the same two assumptions as made in the radiationless case above are used. Although these assumptions are valid in the first order radiationless case i t does not immediately f o l l o w that they are also valid for the second order radiative case. In particular, because of the "sum" over intermediate states for energies up to infinity, the assumption iB not always true. Also, contributions from the neglected components of the wave fns. are possible and i t has not been shown that such contributions are negligible. This would require a separate investigation in itself. However, i f we define an allowed transition to be the one for which only those components of the Dirac wave fns. are retained which have an angular dependence given by-zero order spherical harmonics, then this last objection disappears. Also in Chapter III, Section C, the "Sum" over intermediate states mentioned above is considered. It reduces to a sum over all discrete energies and an integral over all continuum energies. Arguments are presented to show .that the sum over discrete energies should be negligible and it is then omitted. The remaining integral over continuum energies is too complicated as i t stands to be evaluated^ within a reasonable length of time0 In Chapter III, Section D, the integral is carried cut for gg Cs'3' by using what we shall call a "semi-relativisitic" approximation for the solutions to the Dirac Equation with a Coulomb potential. This approximation preserves the normalization of the wave fns. and simplifies the integrand to the extent . that i t can be plotted and evaluated numerically. The integrand is plotted for Cs'3' for four different photon energies in Graphs III.l to III.4. The integration results are tabulated in Table III.2 In Graph III.5 the ratio Wj^k/tsf, is plotted in comparison to' Morrison and Schiff results. For values of Z other than that of C& graphical evaluation of the integral would again be necessary; however, the most complicated parts of the integrand have been tabulated in Table III.l for any value of Z. In Chapter III, Section E, the validity of the "semi-relativisitic" approximation is investigated. It happenls that for It* HI the integral over continuum energies can be easily plotted without making a "semi-relativisitic" approximation. In Graph III.6,,the integrand is plotted for Z - 119 and compared with the curves of the integrand for Z = 119 in which "semi-relativistic" and "non-relativistic" approximations have geen.made. In the "non-relativistic" approximations Schrodinger wave fns. are used. The integrand and the "semi-relativistic" integrand have the same g B n e r a l shape and i t is anticipated that for ss^'*' these two curves are raich closer to-gether; however, the "non-relativistic" integrand does not have even the same shape as the other two - especially at low photon energies where the in-tegrand is largest. "In fact, the value of the integral for the "non-relativistic" approximation turns out to be much larger than for the "semi-relativistic" approximation, and, consequently, the value of the ratio W^V^c i s» i n t h e "non-relativistic" approximation, much closer to the experimental one. This is a ratledf curious result because, from the theo-retical point of view, the "semi-relativislic" approximation should give better results than the "non-relativistic" one. 8 CHAPTER I Outline of Theory for Allowed Transitions In this Chapter a brief review of the theory of electron and positron emission and orbital electron capture for an allowed transi-tion is presented, only those essentials necessary for an understanding of this thesis being discussed. For more details, the reader is referred to the longer articles of Rose (55)* Konopinski (55), Konopinski and Langer (53), Blatt and VJeisskopf (52), De Groot and Tolhoek (50). The radiative processes, i.e. those for which internal bremsstrahlung occurs, are also considered and expressions for their probabilities are formulated along with a^dlccussion of the intermediate states. A, Non-Radiative Processes It is the purpose of this Section to show how the matrix elements for the non-radiative electron and positron emission and orbital electron capture processes are obtained. These matrix elements are then used in Section B for the radiative processes. In the treatment of the non-radiative processes given below, we shall adopt the convention that an anti-neutrino is emitted in electron emission and a neutrino is emitted in positron emission and orbital electron capture. The results for the opposite convention in which a neutrino is emitted in electron emission and an anti-neutrino in positron emission and orbital capture, are the same. We shall also use the word "lepton" when referring to an electron, positron, neutrino or anti-neutrino. 1. Electron Emission The basis for the theory of. beta-decay, as initiated by Fermi (34), is the assumption that the emission of leptons can be treated in a manner 9 analogous to the emission of electromagnetic radiation, .as explained in Fermi's original paper,, As a result, the interaction matrix element is assumed to be: <FiH.io)i t Jarnti, ify'AM it o-» where: if/ is the wave -f h. for an electron ib & positive energy state, <p is the wave -fn. for a neutrino in a negative energy state, Tp\ > are initial and final nuclear states respectively, * is the symbol for the adjoint (i 0e. transpose and complex conjugate) of a matrix, J l ,SL/n are operators operating respectively on lepton and nucleon wave fns., ( ) n means the value of the fn. in the brackets, taken at the th position.of the n nucleon, jdtm is taken over the volume containing the n**1 nucleon i 0 e 0 the nuclear, volume. . ( I i ) In order to write the theory of electron.emission in a symmetric form, we shall say that the initial state of the system 0 consists of the initial nucleus and a neutrino in a negative energy state and the. final state of the system F consists of the final.nucleus and the emitted electron By restricting the interaction matrix element to be relativistically invariant, it can be shown that on quite general assumptions there are only five possible choices for XL , In general, some of the components of -fl-are "large" or non-relativistic whereas the.remaining components are "small" or relativistic of order v/c where, v is the nucleon velocity, the "large" and "small" components giving rise to different selection rules* For an allowed transition, only the "large" components are used, the resulting expression for the interaction matrix element being no longer reiativistrically invariant,, The "large" components are: 0 called the Scalar interaction (I.3a) 1 11 i. Vector " (1.3b) /30s " " Tensor M (l»3c) 5 s " " Axial » (I.3d) lfi*s : " 11 Pseudbspalar » (l.3e) where the matrices /3 , <? , V5- are defined in (k.(e>S > C A - 0 « Note, in particular,. that all the interactions are Hermitian. The Pseudoscalar interaction (l«3e) is actually, a relativistic one of order v/cj however, it is the largest component of its operator Si- because the others are zero and is retained for this reason*. There are nine different matrices given in ( lo3). In order to write a linear combination of the five interactions, i t is convenient to define the following quanitities0 r - 1, 2, 3, 4, 5, 6, 7, 8, 9 A r = fir , / ^ . ^ J . ^ ><*l, ^ 3 . Lpr6 Or z c s , c y , c r, c T, c T, e A, c A, c A, c p, a . * ) With these quantities, the interaction matrix element can be expressed as follows: ( Fi H.|o) = & z i c r £ / ( r A r n , ( * • ; A ; ^ ) ^ A I T T S ) It can be shown that the C's are realo In (1.5) we have imposed the res-triction on the C's, letting G be the coupling constant of the beta interaction. An allowed transition is characterized not only by the neglection of the "small" terms in obtaining the interactions (1.3) but also by an approximation concerning the matrix element of the emitted leptons at the n nucleon in (1.5). In this approximation i t is assumed that the matrix : element involving the lepton wave fns. is slowly varying over the nuclear a.nd volume and hence can be taken outside the nuclear volume integral^sum over the nucleons in (1.5)* The lepton wave fns. which are largest at the nucleus and the largest terms in the expansion of these wave fns. about the origin are then used in the extricated matrix elemento Thus, for an allowed transition, (FlH-10)- frCC^'A"^ {K (r.7a) where: R is the nuclear radius and /A- = t . (r.7b) The nuclear matrix element (1.7b) is assumed to be zero unless a neutron changes into a proton. If we neglect nuclear charge affects and describe the leptons by plane wave solutions to the Dirac Eq. then, for an allowed transition, we approximate expu LjC-r* by 1 at the nucleus and so (lo7a) becomes (FlMo)- CrVi C r ( V A r i l ) /Ar (r.8) r / where ^,t\ aro the matrix amplitudes of the electron, neutrino wave fns„ respectively. Approximating the exponential by unity at the nucleus i s equivalent to expanding i t i n terms of spherical waves and noting that the largest term of this expansion at the nucleus i s the L— 0 term or — ^ — • ^r| } i,e„ that part of the plane wave fn. having zero orbital angular momentum. If, on the other hand, we neglect r e l a t i v i s t i c effects and yet tcUce into account nuclear charge effects by describing the leptons by pherical wave solutions to the Schrodinger 3q<, with a Coulomb potential, then D these solutions (for energies i n the continuum) vary as r near the origin, the largest one being that for J= 0 as i n the plane wave caso 0 We see from the above considerations that an alioitrod transition i s one for which the r e l a t i v i s t i c terms of the interaction hamiltonian are neglected (excepting pseudoscalar interaction) and one for which' the emitted leptons carry off zero orbital angular momentum. Definite selection rules exist for allowed transitions; ho\*ever "forbidden" transitions obeying different selection rules do occur, their theory being obtained from a study of the terms neglected i n the allowed transition,, For further details, see the references l i s t e d at the beginning of thi3 Chapter,, The nucleons and their wave fns« w i l l be treated non-relativistically c Consequently, for Scalar and Tensor interactions we shall use the relations f& =/* JP* (1.1) for the nuclear matrix elements, essentially replacing fi by unity. Such a replacement i s permitted because i n the representation used i n this thesis, the diagonal terns of /3 which are t l connect the "large" non-relativistic components of the nucleon wave fh-s„ whereas the' diagonal terms which are -1 connect the !'small" r e l a t i v i s t i c components which we are neglecting,. 13 2. Positron Emission and Orbital Electron Capture In the theory of positron emission or orbital electron capture, an electron in a negative continuum or positive discrete state is captured by the nucleus with the emission of a neutrino,, From the theory of second quantization, creation and annihilation operators are the adjoint of each other and so the interaction harailtonian for the processes under consideration is the adjoint of the hamiltonian for electron emission. As a result, the interaction matrix element for these processes is obtained by taking the complex conjugate of (1.7)> giving (F!H + 10) = G - £ I C r ( 4 > * A r 4 0 R / A r <*.!©*) where is zero unless a proton changes into, a neutron. In (1.10a) <p is the wave fn. for a neutrino in a positive energy state and ljf is that for an electron in a negative continuum or positive discrete state for positron emission or orbital electron capture respectively., B. Radiative Processes. In this section, the probabilities for radiative electron and positron emission are formulated for use later on in Chapter II. In the calculations of Chapter II the lepton wave fns0 are approximated by plane wave solutions to the Dirac Eq„ neglecting nuclear charge effects (Z = 0) so that the final formulae and the discussion of intermediate states given in this section will apply only for this approximation. Also, the formula for the radiative K. capture probability is presented,, This thesis, from Chapter III on, is devoted to the calculation of this probability, taking into account both relativistic and Coulomb effects i.e. by using solutions to the Dirac Eq. with a Coulomb potential for the lepton wave fns; consequently, the formula developed in this Section for this process is considerably more general than those for the radiative electron and positron emission in which plane wave fns. are used. The intermediate states for this process are also discussed, relations between their matrix elements being used to simplify the radiative K capture formula. I . Radiatiive Electron Tftrri-ssion This process requires an intermediate state. Let us designate the states cf the system as follows: 0, Initial state consisting of nucleus of charge %mt and neutrino with energy - E j j (C-^O), £^ is the energy of the anti-neutrino, lt intermediate state consisting of final nucleus of charge Z I and emitted electron with energy £ s F, final state consisting of final nucleus of charge Zt-J , an electron with energy E, and a photon with energy E f s c k . Conservation of energy between 0 and F states tells us that the available nuclear energy, W, is given by and defining W e = £ + £ y ( I l l b ) we also find that the energy difference between the I cJid 0 states is given by The probability per sec, Pi for the emission of an electron with energy E, an anti-neutrino with energy £7- , and a photon with energy C r is obtained from the formula a well-known result of time dependent second order perturbation theory. Here P is -obtained by integrating over the directions of emission of the photon, anti-neutrino and electron, by summing over the polarizations of the photon, and by summing over all possible states of the anti-neutrino and electron having the energies a.n.4 £~. & is the Dirac delta fn. I* H E r + E j + E ^ - w ) , (r./jb) enforcing the conservation of energy (i.lla); it is used here because i t prevents any confusion which might arise otherwise in obtaining the density of^inal states. The matrix elements in (l„13a) are given by (II H_|O) = G E C r (K*s A" a> /Ar (r./f a) where: -e is the electron charge (C-^o)> t£ , & iV a r e t n e m o m 9 n ' t a °f t n e photon and electron in its intermediate and final state respectively, is the unit polarization vector for the photon, ar-. the Dirac plane wave amplitudes for an electron in a positive energy state E and a neutrino in a negative energy state ~E— , ^ satisfying the Eq. ( c Z - f + /3mc*)£, - E t , (X.is) a»U t[ satisfying the sane Eq. with m a 0 and £ , f> replaced by (1.14a) i s obtained from ( l 0 S ) , va l id for an allowed transi t ion. (I. 14b) •is taken from Heitler (54) Po 144, Eqs. (25) and (27). The momentum conservation i n (l»14b) permits only two energies for the electron i n the intermediate state; these are £TS - i [ <cfc>*+ (me*)1]1*; . <r- / 6 ) being determed by (1.14b). In ( l „13a) the sum over intermediate states I, using (1.12), becomes a sum over the energies (1.16) and the spin orientations of the elctron wave amplitude £ s satisfying (1.15) with E replaced by E*s i nd by We want to calculate'-the probabil ity per sec. , S ( k ) d k that a photon of energy %ck ^ n the range d(Hck)±s emitted during the electron emission. To obtain i t , we must multiply (1.13a) by the density of states of the emitted anti-neutrino and electron and integrate over a l l Ep. an<4 ET and also just multiply the result by the photon density of states. Thus, carrying out the integral over E p f i r s t and converting the remaining integral over E to one over Wft defined i n ( I . l lb ) , we f ind that Ury(Kc) ii^tr where ^ The formulae (1.13) and (1.17) are used in Chapter II to calculate 17 the probability for radiative electron emission for an allowed transition, neglecting nuclear charge effects. 2 . Radiative Positron Emission In this subsection i t is shown that the probability for radiative positron emission, to the same approximations as in subsection 1 above, can be obtained from the probability for radiative electron emission given in subsection 1 naerely by replacing the electron mass m by -m In the results,, This fact is used in Ch. II. For radiative positron emission an electron in a negative energy state jumps to an intermediate energy state with the emission of a photon and is then captured by the nucleuw with the emission of a neutrino0 The formulae (1.13) are valid for this case i f we replace v by >%e by-eaxi& use the matrix elements ( F | n t | r ) - O E C , . (/l"V4's) / A ' ( r . i f i ) < r i M , i o > - - e [ i = i * ] * ( t ' , * -T- *K-K where \\' is the wave amplitude for a neutrino in a positive energy state E^ lv\A 4' -*-s the" wave amplitude for an electron with energy —t momentum ~f> , iT *\nd fy being the po^ Ltron energy and momentum, r ^ i t i ^ CT.IIT3.) is obtained from (r. I0±) (cfff.*))uid £<]. (r.lgb) from Heitler (54), p. 144, Eqs. (25) and ( 2 7 ) . Note that the positron and photon momenta are conserved in (l018b) exactly as in (I.14b). The wave amplitude satisfies (coT- (-£) -h /3n^ic' l)i ' = -£ £' (r.n) E} f» being positron energy and momentum. Eq. (1.19) is the same as Eq. (1.15) with m replaced by - m. If' satisfies (1.19) with m a 0, the result being the same as Eq. (1.15) with m s 0. IS It is easy to see now that, by taking the complex conjugate of the matrix elements (1.18) and replacing /m by -m , e by -e , one obtains the matrix elements (I.1A). • A glance at (1.13a) shows that i t is the square of the absolute value of the sum of the matrix elements over the intermediate states which determines the desired probability and so tha fact that the matrix elements (1.18) are the complex conjugate of the matrix elements (I.14) is irrelevant to the probability; also, since £ is squared, its sign is irrelevant to the probability; hence, only the o f m by -m e replacement fep m is reliant to the probability. This is the result which we set out to prove in this subsection. 2. Radiative K Capture In accordance with the hole theory, two types of process es are possible in the transition from the initial state of the system consisting of the nucleus Nz and Z atomic electrons to the final state of the system consisting of a nucleus N^ -y , one unoccupied K state and 1-1 atomic electrons, a photon and a neutrino. These processes are: I. The K electron makes a transition to either (a) an unoccupied discrete state or (b) an unoccupied positive continuum state with the emission of a photon. The electron is then captured by the nucleus with the emission of a neutrino. II. An electron in either (a) an occupied discrete state or (b) an occupied negative continuum state is captured by the nucleus with the emission of a neutrino. The K electron then jumps into the remaining hole with the emission of a photon. Let us label the initial state of the system by 0, the final state by F and the intermediate states by 1(a), (l)b, 11(a), 11(b) corresponding, to the processes ;;iven above. The conservation of energy between initial and final states tells us that the available energy W is given by: vrhere: is the energy of the nucleus , E is the total energy of the K electron, £ is the neutrino energy t i c k is the photon energy. The probability per sec; **^ A^ , for the emission of a photon with energy Ack in the range d( c A ) is given by the formula * ft (zw)3 * } ' v ' I \ r (FlH » | D ( r H 4 r 1 0 ) r-r (F/"y I ttXl | H j O ) I a result of second order time dependent perturbation theory. The symbol I represents all the intermediate states for the processes 1(a) and 1(b) above; similarly II represents the intermediate states for 11(a) and 11(b). The probability is obtained by summing over the polarizations of. the photon and integrating over the directions of emission of the photon and neutrino. S v repressnts the sum over all possible neutrino states having the energy Ev and 59 is the sum over the two K electron states, giving the average results for the capture of either K electron. The neutrino density of states does not appear because the neutrino wave fn. will be normalized per unit energy interval, automatically accounting for the density of states. Tho matrix elements for the transitions are: ' where: J ^ j - .- Dirac Coulomb wave fn. for the electron i n intermediate states 3(a) or 1(b), ^ = Dirac Coulomb \*ave fn. for K electron, <P - Dirac wave fn„ for neutrino with positive energy } Dirac Coulomb wave fn. meaning a solution to the Dirac Eq. with a Coulomb potential. (F/Kr/1T), (T / H + |o) are given by (1.22a), (1.22b) respectively with I replaced by II on the right hand side of the Equations. The expression (l.22a) i s taken from Heitler (54), Eq. (21b) on p. 143, that for (I.22b) from (I.10a). The formula (1.21) can be simplified by means of the relation £ ( F j ^ 4 J D j ( £ i H ^ = y ( F/Hj r*yra . lMC) + y (p I rir UTa^falH^lQ) ' which we shall now prove. The matrix elements on the l e f t are given by (1.22) with .1 replaced by D and the sum over D means the sum over a l l electron wave fns belongin; to a l l discrete states. Obviously, the sum over D includes the sum over the intermediate states l a i.e . over a l l the unoccupied discrete states, so that i n order to prove (1.23), we must prove that the sum over D for discrete states which are occupied by electrons i s the same as the sum over the intermediate states I l a . Note that the process represented by D i s the same as l a except that we are now concerned with occupied states 0 Thus, using for the initial state energy, Wj, bein;; the nuclear energy and the total K electron energy, then for an occupied state we have E 0 - viz + r s tick, EB being the total energy of the occupied discrete state. The correspon-ding intermediate state in the process 11(a) has energy Now, making use of (1.2$), we find that Furthermore, the matrix elements of D and Ila for the occupied state under consideration are the same, just interchanged; however, the processes represented by D and Ila are initiated by different electrons. In ofcder to , make the process Ila the same as D, the electron initiating Ila must be interchanged with the elctron initiating D, such an interchange producing a minus sign because the wave fn 0 for the system of electrons is anti-symmetric. This minus sign cancels out the difference in the energy denominators given above and hence (1.23) is proved. An argument similar to this one is &Lven by Heitler (54), bottom of p. 213, in connection with the Comptoh effect. The argument given hero could be made much more straightforward by using the . formalism of the second quantization,. In the same way, one can prove that r* (PiH4..cXciH,i<rt _ £ (FfH+UbXrHHylo) ^£-» (rfHj|n.b)(gb,H^,o) c E C - £ T 0 ~ X b . E x b - C 0 J i b £jrb -£0 (T'24) where the sum over C on the left is the surat over all continuum energy states for the electron in the intermediate state. Thus, using (1.23) and (1.24), we may write vrhere we now understand the sum over I to mean the sum over all discrete states an4 all c o v x t i n a u m stales for the electron in the intermediate state I, no matter whether these states are occupied or not. The matrix elements in (1.25) are s t i l l given by (1.22) and since we no longer need to distinguish between I and II, it is convenient to replace by and Jet" where E is the energy of the electron in the intermediate state« M^dk is calculated approximately in Ch.III taking into account Coulomb effects. The radiationless K capture probability , yi is also calculated there to the same approximation and is used for the ratio vv. , such a ratio being a more reliable result than the approximate result for alone. Let us now write down the formula for u/c which is used in Cho III. It is obtained from first order time dependent perturbation theory, i.e. „ a M io J&vK l(P» H + | W | x ' | ( r z 7 ) C t A • the matrix element being given by (1.22b). 5„ a.nA S0 have been defined already after formula (le21); also, no neutrino density of states appears for the same reason as in (lo25)o US Chapter II Radiative Electron and Positron Emission In this Chapter, the calculations for the probabilities of radiative electron and positron emission are shown i n detail for an allowed transition. The lepton wave fns. are approximated by Dirac plane waves and hence nuclear charge effects are neglec'ted,. i.e. Z=0 „ Vie shall take into account a mixture of the five beta interactions and i n this sense extend the results of Bloch (36), Knipp and Uhlenbeck (36) and Chang and Falkoff (49), who have calculated these probabilities to the same approximation as mentioned above, but for pure interactions only. The results show that non-zero contributions to the probabilities appear also i n the cross terms between scalar and vector interactions and i n those between tensor and axial interactions. The calculations for radiative electron emission are carried out i n detail, the results for radiative positron emission being obtained from them by replacing m by -m, as shown i n Ch. I. •Vte shall f i r s t calculate P (E, £Vt t~\ from the formula (1.13a) using the expressions 1(1.14) for the matrix elements and (1.12) for the energy denominator0 The sum over the intermediate states I i s the same as the sum over the energies ±(£ s( and the spins of the electron i n the intermediate state. It i s performed by f i r s t rationalizing the energy denominator so that £ 5 * " ' appears, then by replacing Eb in the matrix element by ^s\& where using (I.14b), H3 being independent of s } and finally by carrying out the sum over the'electron spin, i0e<> the sum over the four intermediate states <£s belonging to ± | E S J , using the relation where 1 is the unit matrix.' This last relation is the closure relation for the complete set of column matrices . The result is: 5 ^ - U / e £ * - U / / with c* + (mc^ . (IT. 3b) Let us now write down some relations of use in connection with operations and averages of anti-neutrino states. Using the projectipn operator which is 1 for positive energy states and 0 for negative energy, states and the closure relation for* the \\'s just as in (II.2), we find the relation = ATT ($* x i i ) (p. 5-*} Z5 for the integral over all directions of emission of the anti-neutrino and the sum over the spins of the anti-neutrino in a positive energy state £-. In obtaining (II.5), we have used the results jdJl- j'-Pp = 0 j 3 a C«n*t»rvt vector In performing the sum over the two spins of the emitted electron in a positive energy state, we shall also make use of the well-known result (see Heitler (54), p. 108, Eq. (l4b)« Using the relations (II.3), (11.5) and (IL.6). in (1.13) and (le14) we find: ^/^V|E C F I M; | , ) t r. ' l ' - ' r tr = (an-Kc)* & where H5 = c*. [?+ tit) tfi ^ H . (lT.7b> A considerable amount of alegebra is required to simplify the product of the matrices appearing in the square brackets of (II.7a), especially for arbitrary values of r ;nd f' . The algebra is given A pp. D. The result is: *• 3L £ ^ . We have yet to perform the operations to find P of (I.13a). Collecting all the quantities affected by these operations in (II.7a) and (II.8), we may write i-Bii^-THF-rtf i-B£(r3.?k)(fe-k) f ^ ^ - ^ ^ ^ . f t Where the B's and a's are unaffected by the operations (II.9). Slimming (H.10)over the polarizations of the photon and integrating over the directions of emission of the electron, we find that (II.10) reduces to k k as shown i n App. D. In (II.11),9 i s the angle between p and A . If we now integrate over a l l directions of emission of the photon i n (II.11) then a l l terms containing ( A-• R ) disappear and we are l e f t with the result: jdSLYj<l/LEILR LLX.U) = frr3- ^LiAi^UB^lB^Kco^^B^^s^H^ ^ The integrals required to evaluate (II.12) are given below. By-substitution one can show that hence, by expressing the fn . in the curly brackets of (LTol2) i n povrers of _ c/>«o*0 -£" and using the integrals we can evaluate (11,12). In connection with the results (11.13), we must remember that when m i s replaced by -in, the second integral , for which the logarithm results, i s unaffectedo Collect ing (II.7a), (II.S), (11.10), (ll<>12) and (11.13) together with ( l 0 13a), we obtain it Pi£, £>, Ev) . if £ Cr i, (A' (/»")' -It is convenient at this point to notice that the correspond-5 ing result for poitron emission, obtained from (II.14) by replacing m by -m everywhere but in the logarithms, is the same as (ll„14) but for a sign change in .0 4 A A /S , Th$ term later gives rise to the cross terms between the beta interactions. Another m will appear in the expression for S(k) from the lower.limit on the integral in (I.17a). but this m is obviously the same for positron and electron emission. In conclusion, the only difference between the electron and positron emission, will be a change in sign for the cross terms for P ( £ - , ^-g-, £ y ) and S ( k ) d h , Let us proceed to calculate S(k)4k for electron emission. From the formula (l„17a), i t is seen that an integral over We has yet to be carried out. The required integrals are given in App., D. In App6 A it is proved that T ArAr' - 4- Kr> (t.tn) and J/i A l\ - 0 except for the cases: (lL./4) Graph II.l <}lf Qg of Eqs. (11.16) for P' 32 (3.84)0^ ("pure" term) (3.84)0^ (cross term) The factor 3.84 was used originally to compare the curves with Goodrich and Payne1 s (54) results. It was determined so that the ordinate for the "pure" term was 146.7 correspon-ding to a photon energy of 485 kev. 2.00 . » 4o© 660 £>00 looo (II.15b) and (11.15c) show that cross-terms arise between the scalar and vector, tensor and axial interactions respectively as in the well-known case of radiationless electron emission (cf. De Groot and Tolhoek (50)). Collecting our results together with the approximation (1.9) for the nuclear matrix elements, one obtains SLK)dh = where: x = £-K <^=—„ K » fWlCb) one1*- m c * - ' / The result (ll 0l6c) agrees with that of Chang and Falkoff ( 4 9 ) , who have carried out the integration over We to obtain Q/ only, i.e. they considered pure interactions only0 As before, the results for positron emission are obtained from (II.16) by changing; the sign of tb p ".ross-terms exactly as in the radiationless case treated by De Groot and Tolhoe.k( 50). 32 A graph of the curves Qt and Q2 for P for which [A/ = 1-71 -mev / w c ' 32 is given in Graph (II.l). The reason for choosing P is that the radiative electron emission of this element has been the object of several experimental investigations (see the Introduction)0 Chapter III Radiative. K Capture • • •», In this Chapter the radiative K .capture probability i s calculated for an allowed transition taking into account an arbitrary mixture of the five beta interactions and, i n addition, Coulomb effects. In order to take Coulomb offects into account, a "semi-relativistic" approximation for the solutions to the Dirac equation with a Coulomb potential i s developed,, The calculations are complicated to the extent that, i n order to obtain the radiative K capture probability for a given photon energy and nuclear charge Z, i t i s necessary to evaluate an integral graphically. The most complicated parts of the integrand involved are tabulated, by means of suitably chosen dimensionless quantities, for various photon energies and for any value of £ so that the required integral, can be evaluated f a i r l y easily for any value of 7j - This integral i s then evaluated expl i c i t l y to find the radiative K capture probability for S5- Ci>3t (neglecting "screening"). Before carrying out the main calculations, we shall f i r s t derive Morrison and Schiff's (40) result for Z-0 , The reason for doing this i s that these authors have just outlined their calculations i n words without showing any details e x p l i c i t l y . From the derivation, we can see exactly what the assumptions and procedure are. A. Morrison and Schiff Results. The main assumption i s that the nuclear charge can be neglected; consequently, the leptons can be represented by pL?ne waves. The K electron i s represented by a plane wave inLth aero momentum ( E K -wic*). 3 i The required formulae for the radiative and radiationless K capture probabilities are: These expressions are taken from (1.25) and (1.27). Since plane waves are used for the leptons, the neutrino density of states has been included. The available energy W i s given by W = £ y + *ck (TIC.A.3) from (1.20). From (1.22), we see that the necessary matrix elements are given by (*r-A-<M>) ( F I H + I I ) . Q J*r l n * A r s ) t where: /\r i s one of the five beta interactions £ f> a.re the wave aaiptitude and momentum for the electron i n the intermediate state 0 • £ satisfies (1.15), the same Sq. with/* = 0 and ^ the same Eq. with/»n=0. From (lII.A.4a), the intermediate electron energies possible are (jtt.A.O Let us evaluate Wc f i r s t using the. matrix element (III.A.4b) with | replaced by £ o . The sum* over the..two neutrino states belonging to the energy , the suia ovov tho two Y. electron states and the integral over the neutrino directions arc carried out by standard procedures (see Heitler (54)). The required relations are: = v r - i XY [±r^ * "^J Using this result with ^ - A " - = j (jr. A.7) i n (HIc-A.2), ( i l l . A . 4 ) , one obtains * " ? ^ I K r • ' toj^A i s obtained as follow;:,. The sura over intermediate states i n (III.A.1) i s carried out f i r s t by siunoin^ over the two spins of the electron having the two energies (lIIoA*5). In order to do this the energy denominator, given by (1.26) with EK = / m c * 2 i s rationalized no that E*~ appears i.e. by multiplying numerator and denominator by E- fak ~<**c*)- then, wherever E*\ appear i n the numerator, i t i s replaced with \\ £ where using the momentum, conservation i n (HI.AJ+a); f i i u l l y , the sura i s carried out usin* Spins , , , r , trie uni'o matr?jx. The re stilt i s : y (F IH+iiXriHyio) r ztrtc 1^ Cre jf\r Squaring this result and performing the K electron and neutrino operations with the help of (ill.A.6), one obtains a result independent of photon diredtions and polarisations so that the photon operations merely introduce a factor 2( TT ). The result i s , using (III.A.7), and so, from (III.A.S) and (III.A.9), wo find Morrison and Schiff's result: This ratio can also be obtained from the results i n Section III.B, which take Coulomb effects into account, by netting Z ^ 0 . However, the value of the ratio so obtained i s smaller by a fiaator two than that given by the Ilorrison and Sch; "." formula (lII.A O10) O It would thus appear that the Morrison and Schiff formula i s not quite correct even for the case 2, - 0 , owing to the very crude assumptions on which i t i s based. 35 B. Radiative and Radiationless K Ca-pture. In this Section, explicit expressions for the matrix elements appearing in w^ fc a.n<l */c of (I„21) and (1.27) are presented, and the sums over the spins and polarization along with the integrals over the solid angles are performed,, As a result, a simpler expression for v/^ A. . and a final expression for \J are obtained. Calculation of VJP The calculation of W is outlined below. Although the procedure followed is essentially the same as that, of De Groot and Tolhoek. (50), it is presented here mainly as an introduction to the more complicated calculation of v j ^ k . From (l027), the probability for the radiationless electron capture per sec. is given by f. where tyK1<P being the wove fns. for the K electron and neutrino, evaluated at the nuclear radius R. The choice of the electron and neutrino wave fns0 is' important, the same type of choice being used in the calculation of . For the K electron, the only wave fns. possible are those from App. C for which and a - ti< , corresponding to the two spin orientations. The K electron wave fns. w i l l be approximated by setting the small radial component G equal to zerio. Thus \*here: i s a four component column vector with unity i n i t s m^n row and zeros elsewhere, Y : -J4-JT F^(R)IR supplies the dependence on the nuclear radius. The neutrino wave fn. belonging to the continuum i s obtained by the following argument. We look for those wave fns. which are largest at the nucleus. The large and small radial.components of the neutrino wave fn., denoted by f and g respectively, vary at most, as (^j:)^ (/JpP«^ ) where A = j + for the neutrino ( 2 = 0 ) 0 At the nuclear radius then, the well-known approximation t e l l s us that the largest wave fns. are those for which j i s smallest i . e . . J = a n d hence K = £/ . It can be shown (cf. Marshak (4?-)) that at the nucleus | varies as ] ^ whereas varies as so that JJ^ may be neglected with respect-to On the other hand, §~ varies as ^ , and |~ as so that we can also neglect ^ with respect to <J" . Hence, the neutrino wave fns. used hero, talcing into account a l l neglections mentioned, are; <PmW - Y: N (A^.C.) (nr.0.5) where yrn's 1,2,3,1} corresponding to j * K = + I, p - ± ^ (4=0) * \ ' " -37 A/ = JL R I » R I • 7Tti 3c 3 ' Note that the sum 5 . appearing in ( i I I o B . l ) denotes the sum over these four states. The matrix element appearing in (III.B.2) now has the form, R using (III.B.3) and (III.B.5). Substituting this expression into (III.B.2) and ( I I I o B . l . ) and carrying • ut the sums 50 i*</ S with the help of the relations? A. one obtains This result is independent of the neutrino direction of emission, so that the integral over the solid angle merely contributes a factor 47T . Using the relations • A r < , - A r V, x/r'/tf - A R Nx (nr.B.s) A/A/ * = _ f i - (5«c ( B - B . S " > one can reduce the trace operation in (ill.3.7) to the evaluation of Now, from App. A, 3d ind TA. f\r'hr/6 - 0 except for cr03s-terias between the scalar and vector, tensor and axial interactions, i'or which (.sr. Q. io b) respectively. Substituting the above results into (III.B.7) and making use of the relations (1.4) and (1.9), one obtains This result nay be obtained froiu De Groot and Tolhoek*s (50) Eq. (61) by multiplying their result by t T to include the integral over neutrino directions and by to include the averaging factor for the K electrons. 3 9 Calculation of w^dh The radiative K capture probability per second is obtained from the formula * * (i»rJ J * y * ; * ' r ^ r - ^ given by Eq. (1.25). The summation over I denotes the sum over all the discrete and positive and negative continuum states corresponding to all the solutions of the Dirac Eq. for the elctron in the intermediate state» In (III.B.12) the energy denominator is given by (1.26) and the matrix elements by (1.22). In the following, we shall label the wave fn , appearing in (1.22), by tye . The K electron and neutrino wave fns. are given by (III.B.3) and (III.B.5) respectively. in the intermediate state The wave fns. possible for the electron^are obtained as follows. It is assumed that at the nucleus, ±E- <^ I (ur.b.13) n as before so that, using a similar argument as in (III.B.4) for the neutrino, we find from A^.O that the largest wave fn. is that for which ji » } K = */. We shall also assume that the largest components of the electron wave fns belonging to these quantum numbers are F G and shall neglect Cr*, F. As in radiationless transitions, the components of the lepton wave fns. which are retained here are these whose angular dependence is given by the spherical harmonic of order zero and those neglected are theones whose angular dependence is given by the spherical harmonic of order one. A word about the assumptions is necessary. For radiationless transitions, the above assumptions are valid and give good agreement with experimental results. It must be realized, however, that radiative processes are second order processes, including a summation oyer the intermediate states of an electron which can have energies up to infinity, so that the assumption (III.B.13) is not always true. In the radiationless case, i t can be shown that the components of the wave fns. neglected supply a negligible contribution .to the probability. It does not follow that this is s t i l l the case for the radiative transitions and a detailed investigation of this point should be made. Of course, the assumptions made may be regarded as defining, by analogy with the radiationless case, an allowed radiative beta transition. We shall use these four wave fns. all belonging to the same energy E in the discrete or continuum case (for the K-state & -0 )t i Expressions for the matrix elements appearing in (III. B.12) may now be given. From Eqs. (III.B.3), (III.B.14), (1.22) we find that / where C K \ n B£K/ 0 use being made of the integral" J • • kr ' 4/r in obtaining (in.B.15). Similarly, from Eqs^ (III.B.5) and (III.B.14) -where * ' " [ 0 ««> & ; < J If we now substitute the results (III.B.16) and (III.B.15) into (III.B.12) and carry out the summation over the four intermediate electron states belonging to the energy E, using the unit matrix, we get y (FiH^irxriH,)o) _ e r i ^ i i j V i y * G- 21 c lA>r ' 1 f r - £0 e L * J r where 71 Or* Q-2Z (3BT.B- f7fc> which, from (III.B.16£) and (HI.B.15b), is a diagonal matrix. Squaring the sum of the matrix elements in (III.B.17) and summing over the two K electron states and the four neutrino states with the help of (III.B.6) one obtains / 41 Let us now obtain an expression for the trace in (III.B.18). Using the result that antioommutes with ^ and I > one can show that fk (i±£)xek = L^. (mr.B.H) Also, since ( ° ° ) and Q is diagonal, Tshe re from (III.B.17b), (III.B.16b) and (III„B.15b). With the help of the relations given in (III.B.8), one can new show that the trace operation in (III.B.18) is reduced to the evaluation of which is obtained from (III.B..10). 13 The above results are independent of the directions of emission of the photon and neutrino and also Independent of the polarization of the photon; hence the remaining operations to obtain ik merely multiply our results by a factor 3. (tf.rr)3'. Collecting all the results together into (ill.B.12) and using the substitutions (l 04) and (1.9), one obtains where, frc '37 (r.io) , (m.Rizb) ' 4 4 C. EXPRESSION FOR L The expression for the radiative K capture probability obtained in (III.B.22) contains in one factor the same photon energy dependence as in 10 Morrison and Schiff 1s expression (III.A.^ and in the remaining factor an additional, dependence due to the presence of IM , where L, as defined in (ill.B.21), is given by B£K = / *n 1 ssair FK*crxLr (nr.cn) The summation dummy E runs over all possible energy values for the intermediate electron, i.e. all negative continuum, discrete, and positive continuum values. It is the object of this Section to evaluate the radial integral B£K and to write out explicitly the summand of (ill.C.la). In performing the "sum" over E, we shall neglect the contribution of discrete states and ' approximate Gg- (fO, Several arguments are presented to justify these approximations in the hope that everything up to the end of this Section can be considered accept-able for an allowed transition* A more questionable assumption is then made in Section III.E where L is calculated in what we shall call a semi-relativistic approximation. We shall neglect the summation over discrete levels for the V following reasons. First, the Morrison and Schiff result'* of Section III.A make use of intermediate electron wave functions belonging to the continuum only. Second, i t is shown in App. C that even when the nuclear charge is taken into account, Q£ is zero for a discrete state in the non-relativistic l im i t . Thus, we may conclude that the greatest contribution to L should come from the "summation" oviir the continuum levels» 4-5" Our neglection of the discrete sum means that we are disregarding those i n t e r m 3 d i a t e states arising when the K electron jumps to an empty discrete level with the emission of a photon and is then captured with the emission of a neutrino and those when a bound electron, other than * the K electron is captured with the emission of a neutrino, the remaining hole being filled by the K electron with the emission of a photon. These transitions are labelled l(a) and 11(a) in Ch.I. Let us now evaluate B£K exactly. The integrals required for the electron in intermediate discrete state are evaluated' in App. B and will not be discussed further0 The radial wave fns. for the electron in an intermediate continuum, state and a K electron, as obtained in Appendix Co, may be written in the form: F / f r ) - Kti.Ej e" *" Car.a*) where ( cf J V C P ' * ' > ct> * " Substitution of the expressions (III.C.2) and (HI.C.3) in (HI.CLla) provides one with integrals of the type given in App. B, Eq. (B.*0) and evaluated exactly there. Using the result (B^ along with the definition one finds that ft = - J ^ e l - a W 1 A//V0 e l " J - / ( ^ i r V f ) [ A ) + UV*)] + •h^K-i tej2£)i uik) + UC-rl] j (jK.C.Hb) In order to obtain an expression for G£ at the nuclear radius, we shall assume that as mentioned previously,, With this assumption we need retain only the lowest power of />R/*r in the expansion of &e(B) in (lII.C,3s). Thus (R) = ajI^Tw *!(*,£) e ^ " { A"* - l l ^ l ^ l ) ^ ' 1 1 (TJTX.5) R where N is given by (lII.C.3b) Assembling the results (III.C.4) and (III.Co5) into (Ill.C.la) with the help of the definition one obtains i.e. the integrand is the same for each integral with the exception that 8 is -1 for the first and for the second. Because of the change in the sign of S , a different fn. of E is integrated over the negative range compared' to that over the positive range. This fact makes it difficult to apply contour methods for the evaluation of L. Note also that no density of intermediate states is required in the integrand since our continuum wave fns. are normalized per unit energy .interval. 4 7 Using the relation and the energy-momentum relation, i t is easy to show'that. Q(E) is real. Let us now write the integral over the negative continuum in (III.C.7) as one over the positive continuum. Using l£i=JE, one readily obtains This*is the expression for L which xre shall use in the next Section where L is calculated in a semi-relavistic approximation. The integrand of (III.C.&) is finite everywhere except possibly at those points, where the energy denominators vanish. No such zeros will occur over the energy range ^c* to oo provided |E"K -tick | * /mci . This relation restricts the photon energy to the range besides the range 0 fick <~ W (UL-C.n) which already exists physically because the photon energy cannot exceed the available energy W for the transition. Since for Cesium, from Saraf (54), the restriction (lII.C.9a) is automatically satisfied; hence, no difficulties in the calculation of L for Cesium arise from the energy denominators, If, on the other hand, i t happendd that 48 then poles would be present and the integration over such poles would have to be carried out using the Cauchy principal value for the integral. D. CALCULATION OF L IN A SEMI-RELATIVISTIC APPROXIMATION The understanding of this Section and the purpose of the next may-be enhanced by the following outline. The solutions to the Dirac Eq. with a Coulomb potential are given in Appendix C. These solutions consist of four components column matrices where each component has an angular dependent factor given by a spherical harmonic and a radial dependent factor given by either F(r)/r or G(r)/r. In the following we shall refer to just F and G as the radial solutions to the Dirac Eq. For positive continuum energy states, the F's are the "large"radial solutions and the G's the "small" in the sense that the G's go to zero as the energy E approaches /wo*1 whereas the F's do not. On. the other hand, for negative continuum energy states, the G's are the "large" and the F's the "small" radial solutions in the sense that the F's go to zero as E approaches -/wic* whereas the G's do not. We shall also speak of the radial solutions to the Schrodinger Eq. with a Coulomb • potential; as above, such solutions, when divided by r, given the complete radial dppendence of the Schrodinger wave fns. . . Rules are formulated in Appendix C which, when applied to the large normalized radial solutions of the Dirac Eq. give the corresponding normalized radial solutions to the Schrodinger Eq., and when applied to the small radial solutions give zero. The rules are formulated for both discrete and continuum "of energy values. After a discussion below/the salienttfeatures of these rules, it is shown that in a non-relativistic approximation i.e. one in which these rules are applied to the radial solution, only the negative energy continuum levels (positrons) contibute to L in Eq. (III.C.8). so L, as obtained in the last Section, is too complicated to be evaluated by any simple means. It is not even easy to determine roughly how the integrand varies as a fn. of E except that i t is finite everywhere. In order to study this question and also to evaluate the integral, i t v/as decided to plot the integrand (for Z = 55, i.e. Cesium neglecting screening) up to some value of E large enough so that the remaining•area to infinity could be evaluated by a rapidly convergent series, the main area being obtained graphically. However, even to plot the integrand within a reasonable length of time, it is necessary to make some approximations. It appears that the most consistent and simplest way of doing this, at. the same time including the positive energy contribution to the integrand, is to make a "semi-relativistic" approximation based on the above mentioned-rules and to adopt as a criterion for a good approximation the requirement that the normalization of the wave fns. be not unduly affected. A plausible argument is given in Appendix C to show that this requirement is satisfied for the semi-relativistic approximation used. A further check on the validity of the semi-relativistic approxi-mation is provided in the next section. It arises from the fact (discussed below) that when A is g - for which value Z, = =H9, a fairly simple exact expression for the integrand of L in (III.C.8) can be obtained. In the next Section, the exact integrand is plotted for this extremely high Z value and compared with graphs of the integrand using: (a) The semi-relativistic approximation for the same Z value and (b) A non-relativistic approximation for Z = 119, obtained by applying the rules mentioned above to the radial solutions of the Dirac Eq. which appear in L. Definite agreement exists between the general shape of the "exact" and "semi-relativistic" curves, but for high energies there is quite a large discrepancy between the ordihates0 It would be surprising i f this discrepancy s t i l l existed to the same extent in the case of Cesium whose nuclear charge is less than half the value considered; however, a proper investigation of this point appears to be quite complicated and has not been the carried out. Oh/other hand, the "non-relativistic" integrand does not have the same general shape as the other two curves and consequently i t appears that a considerable error would be present in the results i f the ."non-relativistic" approximation had been used. Difficulties in Evaluating T. . Let us consider some of the difficulties in evaluation L without making any approximations outside of those made in the previous Section. The main difficulty is finding some convenient expression for the hypergeometric fns. of the form . F ^ + i j r i a A . i A + i ^) (cf (nr.c.4a)) (nr.o.u) where The usual expansion of (HI.D.la) in powers of does not converge forUN(see(B.«)), a situation which occurs for large /> ; thus, several expansions of (HI.D.la) over the range of j> from 0 to « are required i f we try to calculate L (or just plot the integrand) by means of series. Whether'the resulting series converge rapidly or not, has not been investigated. , Besides the series for the hypergeometric fn. a series is also necessary for the complex gamma fn. appearing in (lII.C.3b). k procedure for evaluating L for small photon energies is to first expand s/i Af / h r appearing in the radial integral Bfc-K (ill.Colb) into a series of powers of fc before carrying out the radial integral. Thus, substituting the series S i n h r = ( M r ) ™ into the required type of radial integral • hr * (see (B.dO) ) gives using the same method for evaluating the integral (B.Ao). The advantage of this procedure is that, after using the relations (B-^ ) and (B.J5"), one finds the second parameter in the hypergeometric fn. to be a negative integer, -2n, or zero and the fn. itself is then represented by a finite series (from (Bji)), containing 2n terms for the coefficient of k*"} The integration over E must then be carried out (probably by means of other expansions or else graphically) to obtain L as a power series in k, the coefficients also being fns. of k since k apoears in the energy denominators of L (iil.Cola). 53 This procedure has not been used in this thesis, but it would have to be used if values of L for photon energies smaller than 100 kev in the case of Cesium were desired; however, such energies are below those occurring in the measurements ofSaraf (54). The reason why the method that is used to evaluate L breaks down for these small energies is that the integrand of L is plotted (after being approximated) and L, which is then obtained graphically, decreases for small photon energies to a value comparable to the errors in plotting. It may well happen that special approximations will be required to overcome this same sort of difficulty in the procedure just outlined. Before proceeding to the semi-relativistic approximation, let us first note that i f the quantity <2/\ appearing in the parameters of the hypergeometric fn. (HI.D.la) is an integer, then an explicit finite expression for these fns. can be obtained from the integral representation (B-/3 ) by integration by parts. Now, in our case so that S.t\ , considered as a fn. of Z , is an integer only when Z =<?> 111,137 giving 22\ SjL^ljO respectively. The case Z-0 is of no interest to us and neither is Z = 137 for which theintegrand becomes infinite, but the case ?' - Hi, A = is because the integrand of L as given by (III.C.8) can be obtained exactly and plotted in comparison with a graph of the integrand for the same value of 2 using the semi-relativistic approximation. Such a comparison provides a check on the validity of the semi-relativistic approximation and is carried out in Section E. 5 H-The Semi-Relativlstic Approximation The rules by which the non-relativistic normalized radial solutions of the Schrodinger Eq. (including the nuclear coulomb potential) are obtained from the normalized radial solutions of the Dirac Eq. are listed below for continuous and discrete states. What we call the "semi-relativistic" approximation is then presented and used to evaluate L. Rule for Positive Continuum States (E ^ /wc*): Consider our "large" radial solutions F for these states in Appendix C and regard them as fns. of the momentum f> using £ - ¥ Jjcrjl + (*«.t*)x . Then, in order to obtain the radial solutions to the Schrodinger Eq: (a) Replace E by the first term in its expansion in powers of p i.e. t ^ c / (b) Set A = IK | everywhere. (BCD. 3b) Wherever E appears, it is divided by cb so that, using the relation — C ^ -+• y /mc we ean see what is neglected in this ratio when rule (lII.D.3a) is applied. In obtaining the radial solutions to the Schrodinger Eq., we also write K - i f l or K = -I for K«+C|+-l) or Ks-Q-M[) respectively, both cases producing the same radial solution to theSchrodinger Eq. possessing an orbital angular momentum JL . If the rules (HI.D.3) are applied to the "small" radial solution G for the positive continuum, one obtains zero. Rule for Negative Continuum States (E — - /™ c*): In this case, we again regard the radial solutions as fns. of ^ using Then to obtain the corresponding radial solutions to the Schrodinger Eq: (a) Replace E by the first term in its expansion in powers of f> i.e. t - -/wc* (B) Set *-IK| everywhere. (M-OAb) This rule is obtained as follows. It is shown in Appendix C that for the negative continuum, the G's are the "large" solutions and the F's the "small" in contrast to the positive continuum case. In particular, the negative continuum solutions are related to the positive continuum solutions by the relations F*He | ,* ,2) - e l V G Y - 2 ) (zr.ort) considering the wave fns. to be fns. of E, K, Z, . fi4^ is some constant phase factor. The R.H.S. of (III.D.5) shows that the negative continuum states on the L.H.S. represent the motion of a positron in a positive energy state (with K replaced by -K ) In malting a non-relativistic approximation,the R.H.S. of (5a) must give the radial solution to the Schrodinger Eo.. for a positfcon i.e.e from rule (III.D.3) we just put in the appropriate places, either sign of K giving the same non-relativistic result. The relations (III.D<>5) then show that the rule (III.D«4) must be used radial r i d i a t to obtain the non-relativistic^wave fns. from our negative continuum^ wave fns. JT6 Applying the rule (III.D.4) to our "small" solns F for the negative continuum, the result is zero. This statement also follows from (lH.D.5bj and the fact that G for positive continuum states goes to zero in the non-relativistic limito Rule for discrete states The normalized non-relativistic radial soins for the discrete case are obtained from our radial solutions in Appendix C by setting everywhere and £T = /wc* (BE. C>. l>b) in just the square root factors /"X, which appear in front of the radial solutions. Thus /A" and hence G go to zero. Note however that the quantity t/AB appearing in the wave fns. does not go to zero but J A B ~> ™ C where n is the non-relativistic energy quantum number, The result (HIcD.6c) is obtained by first writing exactly for a general discrete level; then using (III.D.6a), one obtains (HI.D06c). Let'us now discuss these rules and, in particular, how a non-relativistic approximation would affect L in Eq. (ill.C.la). First of all, the rules were determined in each case so that the "large" radial solutions gave the non-relativistic radial solutiorfs for an electron or positron depending on' the energy region. When the same rules are applied to the "small" radial solutions we get zero0 Now L depends only on the "small" radial solutions G~ in the positive continuum and the "large" radial solutions G"in the negative continuum (neglecting the contibution of the discrete states), so that in a complete 57 non-relativistic approximation only the negative continuum states would contribute to L. It is the purpose of the semi-relativistic approximation to include the positive continuum contribution to L. Semi Relativistic Approximation L depends on the wave fns. of the K electron and the electron in the intermediate state, for both of which |K|=I. If we set everywhere in these wave fns. the calculation of L is simplified considerably because, as already mentioned, explicit finite expressions for the hypergeometric fns. and gamma fns. can then be obtained. Setting A = I in the K electron wave fn., one obtains the substitution (III.D.6c) with n = I ; substituting these results in the electron wave fn. (HI.C.2a), we find that the wave fn. contains in one factor the non-relativistic radial dependence and in the other factor the constant N0 which can be taken outside the integrand of L and which then cancels out in the ratio wKik/'wic provided the Kelectron wave fn. is approximated the same way in Wc. Since we are only interested- in this ratio, no further approximations are made for the K electron. Note however, that we shall s t i l l let relativistically, in the energy denominators of L in (III.C.8) since no simplification is required there. The semi-relativistic approximation that is applied to the radial solutions for the electron in the intermediate state is carried out as follows. The square root factors which appear in front of F, G (App. C) arte left unchanged; in the remaining fns, we set }=/ (If. P. 8b) and, considering tnese remaining rns. as fns. of p using replace E by tiie first term in its expansion in powers of p, i.e. as in rules (III.D.3) and (III.D.4). If we also approximated E in the square root factors, we would be carrying out a non-relativistic approxi-mation with the resulting neglection or the positive continuum contribution to L. Levus call the radial fns. resulting from the application of this semi-relativistic approximation, semi-relativistic fns. It is shovaa in Appendix C, that, if we look for general radial solutions to the Dirac radiai Eos. having A = |K( , vie find that such solu-tions are possible i f , and only i f , the Dirac radial Eqs. themselves are approximated, the semi-relativistic fns. being solutions of these approximated f Eqs. which we shall refer to as semi-relativistic radial Eqs. These semi-relativistic Eqs. are given in Appendix C and are used to show, by a plausible argument, that the normalisation of tne semi-relativistic fns. is not unduly affected. A special advantage of the semi-relativistic fns. is that, by the choice of suitable units, the most complicated fns. appearing in the integrand of L need to be tabulated only once as far as Z is concerned, since the resulting table can be used for any Z. It should be mentioned that only the approximation a=| is needed to simplify the integrand of L to the extent where it can be plotted within a reasonable length of time; however, it is not easy to make any statements about the normalization of the resulting wave fns. and also the most compli-cated parts of the resulting integrand of L must be recalculated for each different value of Z. It would be possible though to check the resulting integrand far Zs119, 1-3 with the exact integrand for this ease just as we shall do in the next Section for the semi-relativistic approximation, but, because of time limitations, such a comparison has not been carried out. The Semi-Relativistic Integrand of L We shall now write out explicit^ the integrand of L as given by (III.C.S) using the semi-relativistic approximation. In the next Subsection thiB integrand is plotted as a fn. of the energy E for several values of the photon energy. 7?cA -Carrying out the substitutions (HI.Dofl) one finds I x Le-£K+*tk ' E + eK-Hck / where: ^ p s ' m h >ry and the K electron energy, E K , is obtained from L^. (EL. 0.7). The most difficult quantitities to obtain explicitly in the intewgrand are QCt/^c1) of (UI.D.9d). Qt-«*c*) cannot be obtained from by replacing m by -m or y by -y in the resulting expressions. This is apparent when one notices that the factor ( / + iy ) and not (1+ iiy ) appears in the second term of the right factor of Q in ( )/ Let us now obtain exploit expressions for the Qis. 60-• i Using (B.|<?), one finds that U(K) in (III.D.9e), for be only, cajr^written in the form: Uik) 4- x _ ( — ^ - ' ^ ' J ' . _ _ / ) ( *.„.,„> Before proceeding, let us define the dimensionless quantities: X= , £ hom Cm-0.9 Ci) , ^=Z-«r^.c (nc.o.ll) Then, forming <*£h) + u(-fc) and substituting the result into (lII.D.Qd) xvith T = + / , one obtains where c.c. meaning complex conjugate. Q so defined turns out to be real and positive. A simpler procedure is now used to find Q (-/wic1) from (HI.D.9d). Using the relation (B. lb ), one can prove that (i+Cy)ulk) - ti-iy) ~ ^ j £/* +i(r-M)]'fL'Ll* -i(p+Kr\)T11 j • tm.D.13) from which one can obtain where 1 •' t , n ( 4 x ) * ± / - - cc I (ar.D.j*b) Q - so defined turning out to be real and positive 6/ As seen from the Eqs. (III.D.9), (III.D.12) and (III.D.14), the integrand of L in the semi rel. approximation is quite complicated. In the next subsection we shall plot it up to a certain value of E and then evaluate the "tail" integral by a rapidly convergent series. Before doing this, let us first write the integrand of L in a more conven* ient form. Defining n, = - T e^ T T / q. = TT\coth(rr/a) ± i] fa.D.ir*) sinh rr/q 1 - Z«X far. m i C a. (nr. P. wc) we find from (III.D.9), (HIoD.12), (III.D.14) that I - - N» JL (nr.D.iffJ) In the formula (15b) the D's, and Q's, and energy denominators OS T are all p\§itive so that I is positive or negative depending on whether the negative or positive continuum states respectively contribute the most to the integral. The calculations performed in the next two subsections show that the negative continuum states contribute•the most (as we should expect from a non-relativistic approximation) giving I positive. Tabulation of the Integrand of I I of Eq.. (III.D.15b) must further be written as where J ° ^ = / 7^ ^ £ Lmz.PMi>) being the total integrand in curly brackets in <III.D.15b). I-. is obtained 3t . . x graphically in this Subsection and 1^ , the "tail" integral in the next. It turns but that for Z = 55, tie value ' ea = 3.5- (m-D.n) is sufficiently large for I A to be evaluated analytically, using series expansion of the integrand. The quantities D and Q appearing in the integrand of I ; depend on q'and x only and are tabulated in terms of these variables in Table III.l. "The values of D and Q, given in Table III.l, may be used for any value of Z provided one uses the dimensionless units q and xv With these values, the integrand of I, for different values of Z can be plotted fairly quickly and I ( evaluated graphically. The series expression for the !!;tailn integral I A given in this Section is also valid for any value of Z. Consequently, the results given in this Section can be used to calculate the radiative K capture probability for any value of Zf even though just the results for Cesium have been obtained explicitly here. For Cesium, about 25 values of q were chosen between 0 and 10 for each value of x equal to 0.5, .75, 1, 1.5 in order to obtain enough points to plot the integrand accurately. Far large values of the integrand, the accuracy of the calculations was to. three significant figures; how-ever, for values of £ near £0 the integrand is small in some oases and con-sequently the accuracy decreases to two significant figures. Checks on the calculations were made vrherever possible. In par-e ticular, the integrand of I for f = 1 or q - 0 must be obtained separatly ' A 4 3 from the calculations by a limiting process (the results of which are pre-sented later on in this Subsection), the plotted integrand approaching this point as though it were extrapolated; also, the integrand of I ( for energies near £0 can be checked by the series expression for the integrand of I 4 to an accuracy within 5$. In Graphs III.l, III.2, III.3, III.4 the integrand of I, is pre-sented for x =.5, .75, 1, 1.5. From these Graphs, we see that for small electron energies the positive continuum contribution to the integrand is the significant one even though the largest overall contribution to I comes from the negative energy continuum. As we shall see, the same situa-tion occurs for Z=119 in both the "semi-relativistic" and "exact" curves of the integrand but not for the "non-relativistic" integrand. The values of I, , obtained graphically, are presented in Table III.2. Certain precautions are necessary in calculating the Q»a because we must stay on the correct sheet of the Riemann surface associated with these complex fns. For this reason, the formulae used to calculate the Q*s are given explicitly below. In order to express the Q's in terms of real quantities, we define 1+ L x) = C e z - ~ i L6„ r+ e •0* tzn0+r <j + x 0^ 0 L k T£ b' = e +-G_ + \U r+/r_ ) which, when substituted into (III.D.12b) and (III.D.14b) give Q_ -eT^sinb' (nr. o.jq) (nr. P.a.o) 4f The ranges for the angles 9 + , Q_ given in (III.D.18) are determined as follows. Consider the first term in each of the Q»s in (III.D.12b), (III.D.14b) (the o.c. terms are treated in exactly the same way); the multivaluedness arises from those factors in these terms which are raised to pure imaginary powers. As shown in Appendix B, the multivaluedness appears as soon as we make use of the explicit expression (B.19) for the hypergeometric fns. According to the discussion given in Appendix E, we must choose the ranges of 8+ and 0_ in such a my that 0 A S a requirement satisfied by the ranges in (III.D.18). When q is zero (E= mc*), one can prove from the definitions (III.D.18) and (III.D.19) (respecting the ranges on the e*s) that: X 4- Z * By substituting (III.D.21) into (III.D.20), one finds the values of the Q»s for q= 0. The "Wl" Integral 3jg (mo zi) <The "tail" integral is obtained here as a series in powers of 1/q. The results (III.D.24) and (III.D.25) are valid far any value of Z. Let us first expand the expressions (HI.D.12b)t (III.D.14b) for tbs Q*s. In order that the expansion of the complex binomials in the for-mulae converge absolutely and uniformly over the range of 1^ (3.5 & l ^ x> we must have For Cesium, the value of q corresponding to the energy £0 is 1 GS and since we see that the above relation is satisfied for all q> qc as required,, Expanding these binomials to 4 terms, one finds the Q»s, correct'to 2 terms, to be given by r L ' J Combining these expansions with the D's of (III.D. 15a) and making use of the expansion one obtains an absolutely and uniformly convergent series for the products DQ appearing in the integrand of (IU.D.15b) „ Since also the other parts of the integrand can be expanded in absolutely and uniformly convergent series of powers of 1/q. Consequently, all these series can be multiplied and added together to give a convergent result for 1^ . Writing 1^ in the form the series (III.D.24b) being correct to the terms shown,. Integration of (III.D.S4a) gives l - S J i j l i - i A - A i ] (w.D.zir) one finds where the setels (HI.D.25) converging more rapidly then the series (III.D.24b). Provided \ or q e is suitably chosen, this result is quite general. For Cesium, 1^ is tabulated in Table III.2 using (III.D.22). From the results for I, and 1^ , I = 1, + 1^ is obtained and tabulated in Table UL.2 along with I/x. Using the relation with ~£ given by the series (III.D.24b), one oan check the tabulated values a<f of the integrand of I, for values of q near q^ -8.36. The ratio wj^ dk/wg In the radiationless probability wc, the quantity • ro m (nr. c . i i ) for either a serai-relativistiv or a non-relativistic approximation, (a different N0 for each case), and also Ey= W, the available energy, since no photon is emitted. Collecting the formulae (III.B.ll), (III.B.2Z), (III.D.lbd) and the above results together, one obtains U r r * V I w / (no1)* Onf. O. 173) where Y being unity i f just a "pure" beta interaction is considered. The result (III.D.27a) differs from that of Morrison and Schiff (40), who considered "pure"', interactions only, in the quantity appearing in the l.sr-g —— Photons per disintegration (xlO ) per unit kev energy interval 8 x . Photons per disintegration (xlO ) per unit kev energy interval 131 1 • 0.5-Morrison and Schiff l o o Graph III.5 A comparison of the experimental results of Saraf (54) with the theoretical results of Morrison and Schiff (40), Eq. (III.A.10), and the results of this thesis, Eq. (III.D.27), The results of this thesis are represented by x*s. The ordjUaates of these four points must be divided by 100 for a comparison with the continuous curves on the same scale. 2.00 square brackets. In Graph III.5 our results are compared with the Morrison and Schiff formula. Saraffs (54) experimental curve is also shown. A glance at these curves shows that the theory of radiative K capture alone ,! will not explain the experimental results. Also, we see that taking Coulomb effects into account lowers the ratio w^ dk/xvc by an order of magnitude in the low photon energy region. E. The Case of Z » 119. A = t? As mentioned previously, the integrand of L as given by (III.C.8) can be obtained exactly for this case. We must remember that the integrand of (III.S.8) has been obtained by neglecting contributions of discrete l states and by approximating the electron wave fn. at the nucleus; hence the statement that this integrand can be obtained exactly for Z= 119 means that no further approximations are required. A brief outline showing hew; the exact integrand is obtained is presented here. The calculation of the integrand in a semi-relativistic approximation presents no difficulties since it is completely similar to that given in the last Section for Cesium except that Z must be changed from-55 to/119. The D»s and Q»s of the last Section can be used without recalculation. The exact and semi-relativistic integrand are plotted in Graph III.6 for x =.75, i.e. a photon'energy of 332 kev. Along with these two curves the non-relativistic integrand is drawn for comparison. It is. corresponding obtained by replacing the radial solutions to the Dirac Eq. by the/radial solutions to the Schrodinger Eq.; in this non-relativistic approximation, only the negative continuum states contribute to'L. ' It is perhaps too much to ejfrpect very good agreement between any of the approximate cufeves with thB exact one for such a high value of Z. We can see though that the general shape of the exact and semi-relativistic curves is the same, neither one agreeing in shape with the non-relativistio curve; however, there is a discrepancy between the ordinates of the exact and semi-relativistic curves which cannot be overlooked for this value of 2 . It is assumed that for Cesium, whose nuclear charge is less than half that for this case, this discrepancy decreases to a much smaller value. In the semi-relativistic approximation, we have approximated A = Jl-(l*)x 6y . I in certain places. As mentioned in the last paragraph of Section III.D before the Subsection on the semi-relativistic integrand of L is presented, it is not necessary to mate the above energy approximation although certain advantages in calculating the integrand of L do result from i t . If the resulting integrand for just the approximation 2=1 were also plotted in Graph IIi;6 for Z = 119, then we could definitely state that the discrepancy hetween such a curve and the exact curve must decrease for Z decreasing because the approximation A s 1 is so much better for small Z. It does not follow that this approximation would be better than the semi-relativistic one, because the approximation on A and E might have to be made together. In fact, as shown in Appendix C, if one looks for solutions to the Dirac radial Eqs. for which /1= 1, then solutions are possible if , and only if , the approximation on E is also made. It appears that the normalization of the»Jsejni-relativistic radial solutions is assured (Appendix C); how-ever, nothing can be said about the normalization for theAs/only approxima-tion since we do not have the Eqs. which the resulting fns. satisfy and the fns. themselves are too complicated to easily evaluate the required nor-malization integral. In the remainder of this Section the details involved in calcu-lating the exact integrand for S = 119 are outlined. Using (B.18) in (III. C.4a) with 2=-^ , one obtains q, x, y being given by (III.D.ll) and (III.D.9c). Substituting (III.E.2) into (III.C.6), one obtains where (ar.£-where /3 = 6++L 4 - ft JU >;/r_ 7 0 CUT- ET. 3 c) ( A * #n (M.D.IS))} (nr.tr.s Using (B.3), one finds far' N in (III.D,3o) 7T e Substituting (III.E.3) and (III.E.4) into (III.C.8) v/ith the definition Q ' ( E ) - Q(fiT) (nr.c.4) one obtains 2. TT t< (BL.£. <a.N where I, , x ( * f J (LLL) JL (nr.£. &M It is the integrand of I^which Is compared in Graph III.6 with the integrand of I (III.D.15b) for Z=119. The nuclear radial dependence present in (6a) disappears in the ratio wKdk/wc for Z-119,^= \ and so is not taken into aocount in the comparison (Cf. (III.D.15d)). TABLE I I I o l A table of the Quantities ^ (q) (0^-q^lO) defined by Eq.(lII.D.15a) and tyq,x) (0±"q>10) (x • .5, .75, 1, 1.5, 2) defined by Eqs.(III.D. 12b) and (III.D.14b). D+(4) D_(q) Q+(q,.5) QJ<1,.5) E/mc^ 0 6*2832 .2640 3„915 1.0000 . .2 6.2832 .2641 3.845 1.0032 •4 6.2832 .2617 3.601 1.0128 .6 6.2834 .00019 •2546 3.114 1.0287 .8. 6.2856 .00242 •2419 2.448 1.0505 1.0 6.2948 .01162 .2249 1.770 1.0779 1.2 6.3167 .03355 .2060 1.210 1.1105 1.4 6.3555 .07235 1*1477 1.6 6.4095 .12629 .1697 .4619 1.1893 1.8 6.4807 .19754 1.2347 2.0 6.5S70 .28378 .1396 .2423 1.2835 2.2 6.6665 .38331 1.3355 2.4 6.7776 .49442 .1159 .1186 1.3901 2.6 6.8987 .61553 1.4472 2.8 7.0284 • .74525 .09738 .06284 1.5064 3.0 7.1656 .88241 .08965 .04699 . . 1.5674 3,5 7.5346 1.25145 .07383 .02434 1.7271 4.0 7.9321 ' 1.64893 .06178 .01365 1.8947 4.5 8.3501 2.06695 ..05240 .00820 2.0683 5.0 8.7830 2.49983 .04494 .00519 2.2465 5.5 9.2271 2.94392 .03895 .00343 2.4283 6.0 9.6801 3.39694 .03406 .00236 2.6130 7.0 10.6054 4e32217 .02672 .00121 2.9887 8.0 11.5486 5.26547 .02152 .000688 3.3705 9.0 12.5033 6.22016 .01769 .000420 3.7566 10,0 ^13.4672 7.18400 .01479 .000268 4.1458 7<2 TABLE III.l (continued) 0. Q (q,.75J <Mq,.75) ft (q.,D +• 0 .3663 1.793 .4204 .7343 .2 .3702 1.843 .4306 .7682 .4 .3805 1.953 .4577 .8723 .6 .3899 2.039 .4914 1.0255 .8 .3828 1.986 .5162 1.1848 1.0 .3655 1.752 .5203 1.2751 1.2 .3383 lo402 .4999 1.2554 1.4 .3082 1.042 1.6 .2779 .7410 , .4180 .8560 1.8 .2498 .5181 2.0 .2249 ..3618 .3327 .4711 2.2 .2029 .2543 2.4 .1839 .1814 .2664 .2470 2.6 U673 .1314 2.8 .1527 .09675 .2174 .1336 3.0 .1399 .07231 .1977 .1001 3.5 .1142 .03740 .1591 .0517 4.0 .09490 .02096 .1311 .0288 4.5 .08014 .01254 .1098 .0172 5.0 .06851 .007920 .0934 .0108 5.5 - .0805 .00712 6.0 .05176 .003576 «0701 .00494 7.0 .04054 .001839 .0547 .00256 8.0 .03256 .001033 .0437 .00142 9.0 .02671 .0006232 .0358 ..000856 10.0 .02225 .0003981 .0299 .000542 TABLE III.l (continued) <1 ^(q.1.5) Qjq,1.5) <Mq,2) £_(q,2) 0 .3947 .1425 .2803 .03794 ' .2 .4106 .1508 .3128 .03990 .4 .4565 .1764 .3526 .04644 .6 .5286 .2244 .4202 .05828 .8 .6211 .2946 .5167 .07722 1.0 .7165 .3929 .6413 .10587 1.2 .7941 ..5100 .7892 .14758 1.4 .8292 .6181 .9449 .20559 1.6 .8121 .6773 1.0817 .28007 1.8 .7535 .6620 1.1697 .35801 2.0 .6737 .5846 1.1772 .42099 2.2 .5921 .4792 1.1072 .44492 2.4 .5180 .3753 .9923 .42269 2.6 .4544 .2870 .8629 .36897 2.8 .4012 .2178 .7436 .30397 3.0 .3567 .1652 .6412 .24201 3.5 .2741 .'08573 .4580 .13058 4.0 .2186 .04733 .3456 0 7186 4.5 .1792 .02790 .2728. .04170 5.0 .1500 .01735 .2223 .02554 5.5 .1856 .01640 6.0 .1103 .007678 .1578 .01100 7.0 .08486 .003860 .1189 .005452 8.0 .06738 .002154 .0933 .093000 9.0 .05489 .001289 .0754 .001780 10.0 .04566 .0008182 .0624 .001123 TABLE'III.2 A tabulation of I (III.D.15b), I, (IU*.D.16b), and I4(III.D.25). x .1, 1^ I,, I/x .5 .0040 .0154 .0194 .0388 .75 .0728 .0244 .0972 .130 1 .163 .034 .197 .197 1.5 .383 .057 .440 .294 75 APPENDIX A - EVALUATION OF TRACES Each of the (4x4) Dirac matrices and 07 matrices can be expressed (see A&) as the direct product of any two of the (2 x 2) matrices: s . - ( o ? ) s,-c::) s^tri) s 3 = (rj, uo p the last three being the \auli matrices satisfying the relations , 5;s,= - V * = ' f S* C ' j ' k ; C y c l ' C s,A = J e 5 L = 0 j , K . 1 , 1 , 3 Using this property along with the properties of the direct products, i„e. (*,*b,)(ixKbi) = (a.l^b.bj) (/A.3i^ T \ ( a . x b ) = TVia TVt b (A. 3 b ) where a's and b's are (2 x 2) matrices, and with the relation TA. A 5> - T*. B A (A. «•) one can evaluate the required traces in a simple and direct manner,, An illustrative example is presented at the end of this Section. An alternative procedure would be to express the <J matrices as products of <X matrices e.g. <T = -L<>Cxoi.3 f and then make use of the general rules given in Heitler (54) for the traces of products of Dirac matricesj however, I prefer the method presented here. The interaction matrices A are: r = ' , a, J , * , R> S 7, I , 1 A r * - A r (fi.Sj.) A r» A , ' - A<r, , /S^x, A ^ , <r, ; 1 ^ } <r3 > i / 3 r j- = / ^-^b) Using (A2) and (A3), one can show that yj- - + i ' * l * z « 3 = -(5,xi 0) (A.f») One can also show that /3 anticommutes with o(. and ^ and commutes with <r By writing doxra. all the products required and using (A.2), (A.3), (A.4), (A.6), one can show that: T A ArAr' = 4 o " r r , and T i A r A r ' / 3 = 0 except for the cases T A A ' A a / 3 - T/iA*A ' j8 = 4 (Scalar, Vector interactions) ^ ^ J ) T r v A ^ A " 1 * 3 ^ = ^ A M + 3 ^ ' * t / 5 - f (Tensor, Axial interactions) /Y\ ~ 3, 4 , 5". Tftrample: - - T/t ( s, s0 s3 s,) x ( h Sj *.*) (using (A.3a)) = - U (?,s3 s,) TA. s;5^ (U5il1£ (A»3b)) r , r (using (A.2)) = 0 since 7^.S3=0 from (A.2). APPENDIX B Table of Functions In the table below, z is a complex number, R(z) is the real part of z, and any other variables are complex numbers unless otherwise specified* The arguments of the hypergeometric fns. will always be shown so that no confusion should arise between them and other fns. 1 * denoted by the'same letter F. Gamma Fn. o ru) rti-2) = _ E — , (B.3) Sfh 7TZ + 1 Jo (BJL) - n! ^ n integral or zero (B.2) Beta Fn* ' r(X-h/3) ' o Confluent HyperReoaetric Fn* ^Series expansion: (B.4) 5) convergent for finite z. Integral expression: • B(«,*-«)F(«,/3n) = /0V*t"~' (i>-t)*-*-''tit (B,6) R(/&)>RL«}>0 (MacRobert (54), p.346, Cf. also Magnus (49), p. 88.) Useful relatione: (B.7) d. F ( < c , / s ; h ) = Fl<* + t,P + i >D (B.8) r<L JTH = <*.[F(*+i ,/3 ; vr) - F j jrr>] (B.9) ( (B.7), (B.8) from Magnus (49), p.87; (B.9) proved by substitution of (B.5).) (ar-«)F^ 4 - * R*+»,r+i;*> - yF6x,r ;2> (Magnus (49), p.87) (B.lOa) / 3 . £ r . ( * , / 3 j &) - R * + / , / 3 ; £ ) ] = -H F ( o t ( f r o m (B.8), (B.9)) (B.lOb) Asymptotic expressions: • F(*tfl ; I ) - X. (G; t & - J Wn«re, f Q r l a r g Q ^ JL =. ZK£i £i.+ 0(^)J ( B. l l a ) i re*-*) = ^r"*e*L'"-°<£)] ( B - l l b ) „;>• (Sommerfeld p.795-6) Hypergeometric Fn. Hypergeometric Fn»i Series expansion: If «" or /d is a negative integer or zero, the series is finite. The infinite series converges for . (Magnus (49), p.7) Integral expression: 3^,r-(J)FK(J ,» iO»J * («-*> (/-**) * f O (B. 13) R / r^M/»)>0 (Magnus (49), p.8) The integral on the right is multivalued; in order to obtain the sereis (B.12), we must stay on the Riemann sheet for which (See the discussion of this point below.) Trouble might arise in the integral for z real and greater than 1; however, this situation does not occur in our calculations. • Useful Relations: Fix.fi, r,l) = F(&,«•>*; ?) (B.14) F(«,fl,r\l)'(i--l)BF(p.r-*,it-,i:() (Magnus (49),p.8) (B.15) ix-tiy) F(x-iy,*\, ; l(k)) (*-»;Q F U + ».y XW, Z*i-k)) _ ( B . 1 6 ) Here A , y, p, fck, and/4 are real. (B.16) is obtained from the relation and (3.17), the last relation being obtained by subtracting the 1st and 2nd relations from, the bottom of p.9 in Magnus(49) from each other . Special Cases: (From (B.13)) (B.17) F = ! L (I-2)"°t- I J 0£-RL*)C\ (B.18) Evaluation of B^ in (Ill.C.lbl To evaluate B ^ , it is necessary to evaluate integrals of the type /•to |(fc) = / r * A " ' e-CA* ^ ^ r - M ^ J i - F ( x t _ i / 3 z A f l . i i g r ^ j r ( B < 2 0 ) where yU , A are real and positive and f, k, /3 real. In order to integrate (B.20), first substitute the integral expression (B.6) for the confluent hypergeometric fn., then interchange the integrals - possible because of the presence of the decreasing exponential, and carry out the integration over r using (B.l) and finally, compare the result with (B.13) and (B.14) to obtain f(k)= SUX) — r (« »+*p,a*,*V-l i^TTTT)) . (B.21) This result has been obtained previously by Hulme et al (35). The radial integrals involving the electron in an intermediate discrete state are obtained as follows. Using the relation (B.7), the required integrals can al l be expressed in the form: />*-< e-(/<<-*-<A)r F ( ^ | 4 ^ , ; a | r ) ir (B.28) where A , , Y are. real and positive and cC is aero or a negative in-teger (the special case 'in which 1 is not required because of exter-nal zero factors). Since et is zero or a negative integer the series for the confluent hyporgeomtric fn. is finite; hence, a term by term integration is possible, yielding after comparison with (B.12): / U A . T 1 « ^ K ^ * * + , i F ~ - f t ) (B.23) Riemann Surface and Hypergeometric Fn. Explicit expressions for the hypergeometric fns. used in this thesis are given in (B.17), (B.18), (B.19). Because of the presence of non-integral powers in these expressions, they are multivalued. It is the object of this Subsection to find out which value we must use. Z in these expressions is of the form I = J '' % / + i(^-x) (B*24) from which l-i -/- < (? + x\ 1 + iCl-x) (B.25) U3ing the definitions (III.D.18) without the restrictions on the ranges of the G's; we find for (B.25) ' . (B.26) The quantity (B.26) is the only tiling in our hypergeometric expressions which is raised to a non-integral power and hence it causes the multi-valuedness. In order to obtain tns series expansion (B.12) for our nypergeometric fns., we can expand the binomials involving (B.26) provided we take Lim (,-i)r- lv =| 2 -5>0 for V a complex number. This relation defines the sheet of the Riemann surface on which we must remain; thus, using (B.26), we find (,-*)* e-"<*++V and since the value of this quantity has to be unity as Z (or q) goes to zero,we .must have 9++e„-K) f o r < J - > ° - (B.27) This condition defines a sheet of the Riemann surface and it is satisfied by the range chosen for the 0»s in (III.D.18). The same condition (B.27) is obtained again for the complex conjugate or (B.26) which also appears in our Equations in Section III.D. 8 * APPENDIX C. THE NORMALIZED SOLUTIONS TO THE DIRAC EQUATION VvITH A COULOMB POTENTIAL. The. normalized solutions to the Dirac Eq. with a Coulomb potential, (.- i i L 4 + p +fi y = E y/ i c . i i ) where: = _ i V? ^ (C. I b) are given in the first Subsection below along with the special wave fns. used in the 1 thesis* These solutions have been carried out in detail, although only the details necessary for an understanding of the thesis are presented here. In the next two Subsections, normalization procedures are discussed along with a comparison of the solutions with those given by others. In the last Subsection the 11 semi-relativistic" approximation is presented. ' Summary of Normalized Solutions to the Dirac Equation with a Coulomb Potential Following the procedure given by Kramers (38), one finds that the components of the wave fn. ^ can be separated into angular and radial dependent factors. Their are two independent solutions characterized by the quantum numbers j, m where j repres-ents the total angular momentum and m its z, component. These solutions are: These "two solutions are distinguished by the quantity K = ±( j + ^ ), the -f-sign denoting the solution on the left and the - sign the one on the right. In (C.2), the Y«s are spherical harmonics depending on the polar angles ©, <P and are defined by the Equations•• The fns. F,G depend on the radius vector r only; their superscripts denote to which sign of K the wave fn. composed of them belongs. These radial fns. satisfy the Eqs. the sign of K determining the superscript of F, G. By setting E = EV+ mcz, one can show from the Eqs. (C.4) that for positive energy states and for either sign of Z, the G»s are negligible with respect to the F's in a non-relativistic approximation. Thus, only the first two components of the wave fns. ( G. 2 ) ^ necessary in the non-gs-• (cs) relativistic limit so that, i f these wave fns., in this limit, describe an electron with orbital angular momentum X, then we must have j = for K =+(J + £) wave fns. t-tr\ forK=r-(J + -|) wave fns. in the non-relativistic limit. In solving (C.4) for normalizable solutions, one obtains the following results: Discrete Solutions (O^E^mc^): EC*.'/**)- r ~ T ? *i'»o,/,*,_•_-• & ( r ) « M i l A " U y j 5 r ) A { { + - {-] where: (c.b) /f c /i c F, G are normalized so that (C.7) as discussed in the Normalization Subsection below„ In the non-relativistic limit, it can be shown that besides (C.5), the non-relativistic total quantum number. Continuum Solutions (lEl^ mc*1-): F M r j E u ^ f t > ^ L ^ 1 . ( ^ f e T a # f , , r [ &0>=-<5 where: t —p , X rea-l Here F, G- are normalized per unit energy interval, i.e. if F, G and F', G* belong to the energies E and K* respectively, then or where 4£ is any small energy interval containing E», 8T7 Wave Fns. Used in. the Thesis (c.io) In obtaining F*"f G^one uses V JHt*?;, 7 = 1 , *-&=0 and (B.5 ) In (C.fe). No solutions for K= -1 estist since F~, G" are zero. Continuum Solutions for K,*£ j,. j= Far K =tl, j= 'jf , m = £ £ , these solutions for an electron, have exactly the same form as in (C.10) with F +, G r given by (C.8). For K =-1, 3 = -g, we have (cu) where F", G are obtained from (C.8). The wave fns. of the neutrino have exactly the same form for K s ± 1 as those for the electron waive fns. of the continuum just presented except that the nuclear charge Z and the mass m are set equal to zero in t tte radial solutions F, G, n Normalization Procedures In thi3 Subsection, the normalization procedures are outlined and carried out in some detail for continuum states. Essentially, the procedure followed is the one given by Sommerfeld although the method originated with Kramers(38). For the normalization of discrete and continuum wave fns. the following integrals must be evaluated: (y*A>lr = lull* / " ( i F l 1 * Ifrl^Ur H (discrete) Cc.\z) (w'*ll>J.r¥ = lMr/ a 8(- r /*F + <Lr =S(E-E)( continuum) (Cl3) where IT is the normalization factor. In obtaining the radial integrals (C.12), (C.13), ti» orthogonality of the spherical harmonics (C.3) has been used. In (C.13) the primed quantities belong to the energy E'. Instead of using (C.13) to find N, the following integral J*E fy'*lpd? = \U\X IIEI*?'* F + G'*G) IT* / (c . / 4 ) is used where AE is any small energy interval containing E». The basic formula used to evaluate these normalization integrals is obtained from (C.4) and its complex conjugate Eqs. for an energy E*• If we call the two Eqs. in (C.4) (a) and (b) and the corresponding complex conjugate Eq»s. for an energy E», (c) and (d), and form the quantity (a)F«*r (b)G'*> (c)F+ (d)G then we can obtain the relation S"7 . i d .. /£Z«')(F ' V+&'«6.) (CIS) from which using the fact that F(0) = G(0)= 0. The quantity in square brackets on the right is evaluated at the radius R. From the Eq.(C.16), we Bee that the m normalization integral on the left can be evaluated from the asymptotic values of F, G for large R on the right. (C.16) has been obtained without any restrictions on £ and E*, so that even for the discrete states we can Consider the solutions F, G as fns. of a continuous variable E and make use of a Taylor expansion in powers of (E»-E) to find I'tf + Wlr - itc. Lfr [ r & + I* FJ R (en) In obtaining (C.17). we have used the result Urn [ F*& + & * F ] R = 0 By using the asymptotic solutions for F and G and carrying out the limit on the R.H.S. of (C.17), one can evaluate the discrete normalization integral. For details written in English, see Hill and Landshoff (38). 10 Let us carry out the normalization of Continuum vTave fns. in detail. For convenience, we shall use (in the next few paragraphs only) the quantities: A*|C " and the unnorma.11 zed solutions: (>(/») = -*J17T U n W * {{ -T l (CU) K+ to evaluate =/ (c./9) where F, G are given by (C.18) and A In (C.19), Tf* is contained in the small interval & W. From (B.ll) we find, for largef> f run-is*) (C.iO) Cc.2/i) and defining = re \-;6& ^ a e (cub) we obtain, from (C.21) and (C.18), where (zrt^l** are Thus, from (C.18) and (C.22), we nave, asymptotically, i ~ i n -&JL&-r- a. Jw + o" a l e a c<?5 ^ - $v & Cr ^ ZiSfw^S *T e * SU \ From (C.23), we find (.C.13) using - ( n - H - O y O for large f> and — i — for W near Wf. Using • Heitler (54), Eq. 1, p. 66^6 obtain finally . fr.a*) = <rT*.*r*e7SirAfl . Lei?) Thus, from (C.25), (C.21b), (C.19), and using (C.26), with (C.20) and (C.18), we fing the normalized solutions (C.8) for F, Go Comparison With normalized Solutions Of Other Authors 12 Same of the other authors solve the Dirac Eq.: where / t u t all quantities appearing in (C.27) being defined in (C.l). Let us in-vestigate how to obtain our solutions from the solutions Ip"'. Consider the operator (C. 272.) Making use of the relations one can show that 0<*t- = -<*t0, Ofi-y&O and hence " 0H»=HD , where H is our Hamilton!an, given in Eq. (C.l). Multiplying (C.27a) by 0 from the left and using the relation (C.29), one can show that Of 7 is an eigenfunction of H to the eigenvalue 3. Thus, assuming if are normalized, we must have: (1*8) (C.30) where e is same plase factor which is the same for all four components. The result (0.30) has been used to compare the normalized solutions to the Hamiltonian (C.27b) given by other authors with the solutions given here. The solutions given in (C.2) and (C.6) for discrete states have been cheoked with those of Kramers (38), who uses the same Hamilton!an H, and with those of Hill and Landshoff (38), who use H» (0.27b), and agreement has been found except far a common phase factor e ^ (C.30)» Rose (57) has tabulated the continuum solutions without stating which form of the Hamiltonian he uses; however, if we assume that he used the Hamiltonian H* (C.27b), then agreement exists between the continuum solutions (C.2) and (C.8), for both positive and negative energies, and his solutions except for a common phase factor e t u S in (C.30). THE SEMI-RBLATIVISTIC APPROXIMATION We shall first study how to obtain the normalized radial solutions of the Schrodinger Eq. with a Coulomb potential from our "large" radial solutions to the is given by the "radial solution" divided by r. We shall find that when the rules obtained for the "large" radial solutions, are applied to the "small" radial sol-utions, the result is zero. A semi-relativistic approximation based on these rules is then formulated and the normalization of the semi-relativistic wave fns. studied. Discrete States: Kramers (38), p0 311, states what substitutions are required for discrete states. A restatement in terms of our solutions and notation is as follows: The normalized discrete radial solutions of the Schrodinger Eq. are given by Dirac Eq» For convenience the quantities F^ G are called the radial solutions to of the u/a.i/e the Dirac Eq. although the complete radial dependenceAis given by these quantities of the w±ve £*>. divided by r. Similarly, for the Schrodinger Eq., the complete radial dependence^ (C.3»b) the complete radial dependence being FmjL/r. This solution is obtained from either (i.e. for K=± (j+i)) of our "large" discrete radial solutions of the form. by first setting E = mc^ in just the second factor J~& then everywhere else setting A = iKl (instead of A * ^K^-ft*)*" )* and n = n' + |K| , where K is either K = £-f- 1 or K = -£ (corresponding to those solutions for which K = 4- ( j. jg) or K = -(j +• )^ respectively), either choice producing (c.l) with orbital angular momentum it. One may check these statements with the help of the relations (B.I 0 ) for the confluent hypergeometric fns. and the exact result which, when substituted into JAB with ^ = IK| , /n * lK| (Cf gives Possihle minus sign factors may appear in obtaining (C^l) but these are irrelevant and hence omitted. Applying the above rules to our "small" radial solutions G in exactly the same way, one obtains zero because the factor J/\ in front of G goes to zero. • Positive Continuum States For the case of positive continuum states, the radial solution to the Schrodinger Eq. is given by X (from Sommerfeld, pp. 115»-127) which is normalized per unit energy intervale Again, the complete radial dependence is E^/r. This expression is obtained from our "large" continuum solutions in a manner somewhat different from that for the discrete states. The expression (C.32) is a positive energy solution, so we must put £=.+ 1 in our expression for F. How, consider the large radial solutions F as fns. of p using then, in order to obtain the corresponding radial solution to the Schrodinger Eq: (a) Replace E by the first term in its expansion in powers of p, i.e. E - > A n c A (C.33i) (b) Set ^-|K| everywhere. (C-33b) Wherever E appears, it is divided by cp so that, using the relation we can see what is neglected in this ratio when rale (C.33a) is applied. Both values of K = l\ 1 or K=-J? produce the same result (C.32) having or-bital angular momentum Jt. Again, tne relations (B.10) are useful for cheeking this procedure. Also, possible minus sign factors have been omitted. If these rules are applied to our small solutions G in exactly the same way, one obtains zero because the factor (J E| -mcA) * in front of G goes to zero. Negative Continuum States We shall first show that for negative continuum states the G's are the "large" solutions and the F*s the "small". . Let us write our radial solutions as functions of E, K, 2. Then, for negative energy states, E=-|E| and S- -1, and we have for an electron, F(-|E| , K, Z), G(-|E|, K, Z) It is not difficult to prove by direct substitution in J, G (or from the Eqs. (C.4)) that these radial solutions are related to the positive energy radial sol£uions according to the relations: F(-|E|, K,. Z)=et-4> G* (|Bl, -K, -Z) (C.34a) G(-|Elf K, Z)s F* (|E/f -K, -Z) (C.34b where qp is real and constant. The parameters on the S.H.S. of (C.34) show that the radial solutions there represent a positron in a positive energy state. Since far a positive energy state, F is large and G is small (independently of the sign of Z), then the relation (C.34) tells us that for an electron in a negative energy state and represented by the radial fns. on the L.H.S. of (C.34), the G»s are "large" and the F»s are "small". In order to obtain the non-relativistic radial solution from the large solutions G for the electron in a negative energy state, ire see from (C.34) that this is completely equivalent to obtaining the non-relativistic IS radial solution for the F»s for a positron in a positive energy state. If tire replace z by -Z in the radial Schrodinger solution (C.32), then we , obtain the non-relativistic radial solution for a positron in a positive energy state. This solution is then obtained from (C.34b) by setting, just as for the electron in a positive continuum state, both signs of K giving the same result. The corresponding rule far the "large" radial solution G for an electron in the negative continuum then, from (C.34b), must be to set E^lEl-^-mc*- and A=|K| (C3f) in exactly the same places as for the positive energy ease. As before, the relations (B.10) can be used to check this procedure and possible minus sign factors have been omitted. If the rule (C.35) is applied to the "small"1 radial solution HtjUn exactly the same way, one obtains zero because the factor (|E|-mci )^ in front of F goes to zero. The Semi-Relativistic Approximation for Continuum States. Let us study the Dirac radial Eqs. to find what effect setting A = |K| and E-mc*1 in some of the places mentioned above has on them. The exact radial Eqs. are: from (C.4). We shall not approximate the square root energy factors appearing in front of F, G because , for the "small'1 solutions, these factors go to zero in a non-relativistic approximation and we do not want such a result in the approximations which we shall make. F ^ V I E M ^ G * r, * f\ Let Substituting (C.37) into (C.36), we find that F, , G satisfy ••'[A; ^ J * - . ^ ] & . these Eqs. s t i l l being exact. Comparing (C.37) with our solutions in (C.8), we see what the expressions for F( , G, are and , in particular, from which parts of Eq. (C.38) the quantities t , K, Z0CE/op, Ztfmc/p seem to arise. Le€ us now force A to be |KJ, i.e. the smallest total power of r ( for series solutions for F, G in (C.3/) is |K| . This requires that f, = z : « c c . 3 7 ; ' (at) 00 Cc.ii) with-non-trivial solutions far a„, b Q. Substituting (C.39) infco(C38) and equating the coefficients of like powers of r, , we find far r , the Eqs: cb 100 and since we require non-trivial solutions for a , b c in (C.40), we must have ( A + K K A - K ) * - + (cm) Wow, since \ « |K(, the L.E.S. of (C.41) is zero. In order that (C.41) .be satisfied we must have either Z=0 (irrelevant) or E=«£mc'L • Conversely, if Essemc*-, we must have A="|K| (since 0 from (C.37)). These substi-tutions constitute the semi-relativistic approximation. If we set in (C.38) /\ = |E| and E = 6me4 we find the approximate Eqs. . |j£ = -t zjms. C£±i2 and, from the definitions (C.3?), we find our semi-relativistic solutions $, h given by and satisfying the Eqs: thenunnormalized" solutions^ ,|f to these Eqs, being the semi-relativistic fns. obtained from our exact normalized solutions F , G by setting >= | K ! and E = rac everywhere but in the square root energy factors in front of the solutions. ( c . < u ) 1 , > (cm) (cit) 101 From these considerations, we see that in order to find the radial solutions to the Dirac Eqs. for X = |K| we must approximate the Eqs. by-setting E =6 mc* in appropriate places and conversely, the solutions to the approximated Eqs. being the semi-relativistic radial solutions. Normalization of the Semi-Relativlstic Radial Solutions For E near Sme*" , our semi-relativistic "large" radial solutions become almost equal to the normalized Schrodinger radial solution , the "small" solutions going to zero, and consequently, the semi-relativistic fns. , 5/ are normalized in this limit. This statement is obvious when one realizes that the only difference between the"large"semi-relativistio radial solution and the normalized Schrodinger radial solution lies in the square root energy factor in front of these solutions. If we use Eqs.(C.44) in exactly the same way as in the Normali-zation Subsection of this Appendix, we find that the exact normalization of 5-, if can be obtained from the relation -LA. is** + JrVj = C£-£'> + h"is\ the primed quantities belonging to the energy E». Considering only the positive energy case, the negative energy case being similar, we find for the normalization integaal: /on ,2. 7A£T « taking E , + raci outside the last integral (permissible because A E is small). The first integral on the R.H.S. is unity as one may check by following through the normalization procedures in this Appendix for the exact Dirac radial solutions in the continuum setting A S | K | and E — mc1" everywhere but in the square root energy factors. The second integral on the R.H.S. is too complicated to obtain even an order of magnitude for it ; however, the following argument seems to show that it is negligible. Remember that when E is near me , becomes the normalized Schrodinger radial solution and & goes to zero. This means the second integral on the R.H.8. of (C.46) is negligible when E is near mcf. Let us now compare the integrals involving the & 's only, and show that the one on the R.H.S. of (C.46) is negligible with respect to the one on the L.H.S. The factor in front of the R.H.S. integral is of the order Z l« X ~ (C~'° for E near mc*", this factor decreasing as E increases. The R.H.S. integrand 103 is finite everywhere as a fn. of r since the radial dependence jf-/r varies with r as ths non-relativistic Schrodinger radial solution divided by r does. If we carried out the energy integral first in both the integrals under consideration, the resulting integrands would have exactly the same energy dependence but would differ in their radial dependence by the factor r. lor this reason, the variation of the integrals with respect to the energy would be the same and since ths second integral on ths R.H.S. is negligible for E near mc^ , it must s t i l l be negligible for higher energies especially since the factor (C.47) multiplying it decreases as 1/E3* • 10H-Appendix D Derivation of Eg. (II.8) The object i s to simplify the expression <H, + weK.?A(<4f.)?.?fc<„(Wt) " (0>|) For convenience, the following definitions and units are used i n this appendix Also, the following relations w i l l be used: H a» £* • - " (D.3e) R s r H + « ( D . j f ) r • <* = ° (p. 3 M being any vectors. A l l these relations but the last follow from the commutation relations for the Dirac matrices. Expanding (D.l) i n powers of We, one obtains (H s+W^ x ( » s + W e) = H S X H S + l/v/ e ( M f xHs) + W^X ( p . f i ) where 105 IB £ using (Do3a), (D.3c). Simplifying the coefficient of We by the use of (D,3e), (D„3f), (D»3g), D.3d), (Di3h), one finds M + X H i . - £«it££L'+ ztfg? >. i-t . (P.s) The f i r s t term i n (D„4a) can be written i n the form E using (Do3f), (D.4b). The terms of (De6) are simplified as follows: the f i r s t term reduces to using (D.2), (D.3C), (Do3e) ^\3t\\ ; the second tern reduces to . using (D,3a), (D B3b), (D„3c) and (D„3h)j and the last term reduces to Co.3a),(o.3e) usinc aid (D.3c). Collecting the formulae (D.4), (D.5), (D«6) and. (D.7) together and changing back to the original units (0,2), one obtains (II.8). Integrals involved In obtaining Eq. (i l . l l ) In order to carry out the operations , we shall let be fixed in direction along one axis of an orthogonal » _ » • —* set of unit vectors u, , ux ° The momentum f> will be in a direction subtending the angle & with the U3 axis and an azimuthal angle <f to the axis. £^ can be either TT or . Subscripts 1, 2 , and 3, will refer to the components of a vector along the corresponding axes. The sum over the photon polarizations is carried out first, then the integral over just the angle <P for the directions of emission of the electron, The integrand turns out to be either independent of <^ , in which case a non-zero result occurs, or linear in sin<f or cos<P , in which case the integration over <p is zero. Two examples are sufficient to understand the procedure. The expression, when integrated over <f gives zero. Also, using in the expression one finds ^ —•> — » since ft is along the t/j axis. l o l Integrals connected with Eq. (II.16) The integrals oyer WQ which are required can be written in either of the forms j ( y - x ) A y r^Tl ( 0- 8b) by means of the substitutions 4 The change of variable y = cah u U - L .^ + vy is enough to express the integrands of (D08) in.a recognizable form, i.e.i as a product of U and a polynomial in cosh u and i i ' n h * An integration by parts ^ ets rid of the U term, leaving an integration of a polynomial in cosku and i i n f t U which can be evaluated by standard procedures The results are: (p.i) (D. 9 b) 108 References Aston, G. H. (27). Proc. Cam. Phil. Soc. 23, 935 (1927). Blatt, J. M. and Weisskopf, V. F. (52). "Theoretical Nuclear Physics", John Wiley and Sons, 1952. Bloch, F. (36). Phys. Rev. 50, 272 (1936). Bolgiano, P., Madansky, L., Rasetti, F. (53). Phys. Rev. 89, 679 (1953). Bradt, H. et al (46). Helv. Phys. Acta 19, 222 (1946). Chang, C. S. W. and Falkoff, D. L. (49). Phys. Rev. 76, 365 (1949). De Groot, S . R . and Tolhoek, H. A. (50). Physica 16, 456 (1950). Fermi, S. (34). Z. Phys. 88, 161 (1934). Glauber, R.J. and Martin, P. C. (55). Phys. Ret* 95, 572 (1954). Goodrich, M. and Payne, W. B. (54). Phys. Rev. 94 , 405 (1954). Heitler, W. (54). "The Quantum Theory of Radiation", 3rd Ed., Oxford, 1954. Hill, E. L. and Landshoff, L. (38). Rev. Mod. Phys. 10, 87 (1938). Hulme, H. R., Mc Dougall, J., Buckingham, R. A., Fowler, R. H. (35). Proc. Roy. Soc. A, 149, 131 (1935). Knlpp, J. K. and Uhlenbeek, G. E. (36). Physica 3, 425 (1936). Konopinski, E. J. and Langer, L. M. (53). Ann. Rev. of Nucl. So. 2, 261 (1953). Konopinski, E. J. (55). Ch.X of K. Siegbahn»s "Beta and Gamma Ray Spec-troscopy", North Holland Publishing Co., 1955. Kramers, H. A. (38). "Quantentheorie des Electrons und der Strahlung", Leipzig (1938). Liden, K. and Starfelt, N. (55). Phys. Rev. 97, 419 (1955). Mac Robert, T. M. (54). "Functions of a Complex Variable", 4th Ed., lilac Millan, 1954. Madan&ky, L. and Rasetti, F. (54). Phys. Rev. 94, 407 (1954). Magnus, W. and Oberhettinger, F. (49). "Formulas and Theorems for the Special Functions of Mathematical Physics", Chelsea, 1949. M&rskak. R.EL. (4a). Pfr^s. Rev. 61,411 (1942.). Morrison, P. and Sehiff, L.I. (40). Phys. Rev. 58, 24 (1940). Rose, M. E. (37). Phys. Rev. 51, 484 (1937). Rose, M. E. (55). Ch.DC of K. Siegbann's "Beta and Gemma Ray Spectroscopy", North Holland Publishing Co., 1955. Saraf, B. (54). Phys. Rev. 94, 642 (1954). Sommerfeld, A. "Wellenmeehanik", Ungar Publishing Co,, New York.
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- On the theory of radiative Beta processes
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
On the theory of radiative Beta processes Hess, Forest Gene 1955
pdf
Page Metadata
Item Metadata
Title | On the theory of radiative Beta processes |
Creator |
Hess, Forest Gene |
Publisher | University of British Columbia |
Date Issued | 1955 |
Description | This thesis is devoted to some aspects of the theory of the weak continuous gamma radiation (often called 'Internal Bremsstrahlung') which accompanies beta processes, i.e. negative and positive electron emission and orbital electron capture. Whenever a beta process is accompanied by this gamma radiation, it will be called a "radiative" beta process, otherwise "radiationless". The results presented in this thesis go beyond those obtained by other authors in two respects. In the first place, the radiative beta emission probability is calculated for an allowed transition taking into account an arbitrary mixture of all the five beta interactions. Previously, only the theory for the case of pure interactions has been carried out. In the calculations, as in previous ones, Coulomb effects have been neglected. In the second place, the radiative K capture probability is calculated for an all owed transition taking into account again an arbitrary mixture of the five beta interactions, and, in addition, Coulomb effects. Previously, only the case of pure interactions with the neglection of Coulomb effects has been considered. In order to take Coulomb effects into account, a "semi-relativistic" approximation for the solutions to the Dirac equation with a Coulomb potential has been developed. It turns out that taking Coulomb effects into account reduces the probability of radiative K capture by an order of magnitude. |
Subject |
Beta rays |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2012-02-07 |
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.0085928 |
URI | http://hdl.handle.net/2429/40533 |
Degree |
Doctor of Philosophy - PhD |
Program |
Physics |
Affiliation |
Science, Faculty of Physics and Astronomy, Department of |
Degree Grantor | University of British Columbia |
Campus |
UBCV |
Scholarly Level | Graduate |
AggregatedSourceRepository | DSpace |
Download
- Media
- 831-UBC_1955_A1 H35 O6.pdf [ 4.99MB ]
- Metadata
- JSON: 831-1.0085928.json
- JSON-LD: 831-1.0085928-ld.json
- RDF/XML (Pretty): 831-1.0085928-rdf.xml
- RDF/JSON: 831-1.0085928-rdf.json
- Turtle: 831-1.0085928-turtle.txt
- N-Triples: 831-1.0085928-rdf-ntriples.txt
- Original Record: 831-1.0085928-source.json
- Full Text
- 831-1.0085928-fulltext.txt
- Citation
- 831-1.0085928.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:
https://iiif.library.ubc.ca/presentation/dsp.831.1-0085928/manifest