{"http:\/\/dx.doi.org\/10.14288\/1.0085185":{"http:\/\/vivoweb.org\/ontology\/core#departmentOrSchool":[{"value":"Science, Faculty of","type":"literal","lang":"en"},{"value":"Physics and Astronomy, Department of","type":"literal","lang":"en"}],"http:\/\/www.europeana.eu\/schemas\/edm\/dataProvider":[{"value":"DSpace","type":"literal","lang":"en"}],"https:\/\/open.library.ubc.ca\/terms#degreeCampus":[{"value":"UBCV","type":"literal","lang":"en"}],"http:\/\/purl.org\/dc\/terms\/creator":[{"value":"Ho, Ernest Chun Yue","type":"literal","lang":"en"}],"http:\/\/purl.org\/dc\/terms\/issued":[{"value":"2009-07-30T23:10:24Z","type":"literal","lang":"en"},{"value":"2001","type":"literal","lang":"en"}],"http:\/\/vivoweb.org\/ontology\/core#relatedDegree":[{"value":"Master of Science - MSc","type":"literal","lang":"en"}],"https:\/\/open.library.ubc.ca\/terms#degreeGrantor":[{"value":"University of British Columbia","type":"literal","lang":"en"}],"http:\/\/purl.org\/dc\/terms\/description":[{"value":"The rate of the nuclear reaction \u00b3He + \u03bc- \u2014> \u00b3H + v\u00b5 + \u03b3 has been calculated\r\nby two different methods: the elementary particle model (EPM) approach\r\nand the impulse approximation (IA) approach. The exclusive statistical radiative\r\nmuon capture (RMC) rate for photon energy greater than 60 MeV,\r\n\u0393rmcstat, using the elementary particle model approach is found to be 0.2113 s- 1\r\nand the ordinary muon capture (OMC) rate is 1503 s-1. Several trinucleon\r\nwavefunctions from different types of realistic nucleon potentials are used for\r\nthe impulse approximation calculation and it is found that the capture rates\r\ncalculated via the IA exhibit slight model dependences, possibly arising from\r\ndifferences in binding energy predictions, nature of potentials used or partial\r\nwave properties. The impulse approximation version of \u0393rmcstat ranges from\r\n0.1313 to 0.1387 s-1 and the corresponding OMC rate ranges from 1260 to\r\n1360 s-1. The difference in reaction rates between IA and EPM is larger in\r\nRMC due to some of the extra Adler and Dothan terms. Therefore, a meson\r\nexchange current (MEC) calculation of RMC seems necessary to account for\r\nthis discrepancy. [Scientific formulae used in this abstract could not be reproduced.]","type":"literal","lang":"en"}],"http:\/\/www.europeana.eu\/schemas\/edm\/aggregatedCHO":[{"value":"https:\/\/circle.library.ubc.ca\/rest\/handle\/2429\/11546?expand=metadata","type":"literal","lang":"en"}],"http:\/\/purl.org\/dc\/terms\/extent":[{"value":"3107819 bytes","type":"literal","lang":"en"}],"http:\/\/purl.org\/dc\/elements\/1.1\/format":[{"value":"application\/pdf","type":"literal","lang":"en"}],"http:\/\/www.w3.org\/2009\/08\/skos-reference\/skos.html#note":[{"value":"R A D I A T I V E M U O N C A P T U R E B Y 3 H e by Ernest Chun Yue Ho B . Sc., The University of Br i t i sh Columbia, 1998. A T H E S I S S U B M I T T E D I N P A R T I A L F U L F I L M E N T O F T H E R E Q U I R E M E N T S F O R T H E D E G R E E O F M A S T E R O F S C I E N C E in T H E F A C U L T Y O F G R A D U A T E S T U D I E S (Department of Physics and Astronomy) We accept this thesis as conforming to the required standard T H E U N I V E R S I T Y O F B R I T I S H C O L U M B I A A p r i l 27, 2001 \u00a9 Ernest Chun Yue Ho, 2001 In presenting this thesis in partial fulfilment of the requirements for an ad-vanced degree at the University of Br i t i sh Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for extensive copying of this thesis for scholarly purposes may be granted by the head of my department or by his or her representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. Department of Physics and Astronomy The University of Br i t i sh Columbia' Vancouver, Canada Abstract ii Abstract The rate of the nuclear reaction 3 H e + yT \u2014>\u2022 3 H + f M + 7 has been calculated by two different methods: the elementary particle model ( E P M ) approach and the impulse approximation (IA) approach. The exclusive statistical ra-diative muon capture ( R M C ) rate for photon energy greater than 60 M e V , T j J ^ , using the elementary particle model approach is found to be 0.2113 s _ 1 and the ordinary muon capture ( O M C ) rate is 1503 s _ 1 . Several trinucleon wavefunctions from different types of realistic nucleon potentials are used for the impulse approximation calculation and it is found that the capture rates calculated v ia the IA exhibit slight model dependences, possibly arising from differences in binding energy predictions, nature of potentials used or partial wave properties. The impulse approximation version of ranges from 0.1313 to 0.1387 s _ 1 and the corresponding O M C rate ranges from 1260 to 1360 s _ 1 . The difference in reaction rates between IA and E P M is larger in R M C due to some of the extra Adler and Dothan terms. Therefore, a meson exchange current ( M E C ) calculation of R M C seems necessary to account for this discrepancy. Contents i i i Contents Abstract i i List of Tables v List of Figures vi Acknowledgements viii 1. Introduction 1 1.1 What and why of radiative muon capture by 3 He? 1 1.1.1 What? 1 1.1.2 Why? 2 1.2 From ordinary muon capture to radiative muon capture . . . . 3 1.2.1 Ordinary Muon Capture 3 1.2.2 Radiative Muon Capture 6 1.3 Summary 17 2. The elementary particle model 18 2.1 Introduction 18 2.2 Kinematics and nuclear form factors 18 2.2.1 Nuclear form factors 19 2.3 Results 22 2.3.1 Ordinary Muon Capture 24 2.3.2 Discussion 27 2.4 Summary 28 3. The impulse approximation 29 3.1 Introduction 29 3.2 \"Internal\" and \"external\" degrees of freedom 30 3.3 Wavefunction 31 Contents iv 3.3.1 Antisymmetr izat ion and normalization of wavefunctions 32 3.4 Hamil tonian 33 3.5 Nucleon form factors 33 .3.6 LA o E P M translation . . 34 3.6.1 Operators 38 3.7 Results . 42 3.7.1 Ordinary M u o n Capture 42 3.7.2 Discussion 46 3.8 Summary and conclusion 50 Appendices A. Handling of radial wavefunction integrations in the IA . . . 53 A . l Introduction 53 A . 2 Integration of type: C(p, r) = ftp(p, q)ji(qr)q2+adq 53 A . 3 Integration of type: f((p,r)C(p,r)ji(sr)p2dpr2dr 54 A . 4 \"Proofs\" of the accuracy of spline interpolation 54 A . 5 Integration cutoffs 55 B. Second order terms for IA E P M translation 56 C. Wavefunction characteristics 60 C . l Channel specifications 60 C.2 Mode l dependent quantities of wavefunctions 60 D. Notations and conventions 62 Bibliography 64 List of Tables v List of Tables 2.1 A set of ( r 2 ) 2 values. The ( r 2 ) notation indicates that the above quantities are root mean square values of nuclear radii . Experimental data from Ottermann et al . [14] ( 3 He) and Beck et al . [13] ( 3 H) 20 3.1 Various O M C and R M C statistical rates when gP or FP is at the P C A C value. The numbers in the second column are the values obtained by using the \" | \" prescription up to O(-jjfr-) terms. The numbers i n the th i rd column are obtained v ia the correct approach up to 0(-^-) terms. The fourth column has values of the O M C rates using \" | \" prescription up to 0((^-)2) terms. Numbers in the fifth column are the values obtained by the correct approach up to 0((-^)2) terms; the sixth and rightmost columns contain r ^ ( \u00ab > 5 M e V ) and r \u00bb S t (K > 6 0 M e V ) respectively 46 C l Specifications for 22-channel wavefunctions 60 C.2 Some important quantities of trinucleon wavefunctions. B i n d -ing energy Eb in M e V . P{S) denotes the probability of S-wave and so on 61 List of Figures vi List of Figures 1.1 The Feynman diagram for ordinary muon capture 3 1.2 The Par t ia l Conservation of A x i a l Current Hypothesis 5 1.3 The vir tual pion exchange diagram, a consequence of the P C A C 7 1.4 The external radiating diagrams 8 1.5 Covariant derivative of the external weak hadronic currents illustrated. The Aa part, in one photon l imit , only affects the non-radiating weak hadronic current (i.e. the muon radiating diagram) 9 2.1 The E P M description R M C photon spectra using the full Adler Dothan amplitude for various values of Fp. The \"Kl ieb and Rood\" values are taken from the relativistic calculation of ref. [5] which are not shown in ref. [4] 25 2.2 Sensitivity of T ^ ( K > 60 M e V ) using E P M wi th respect to FP. 26 2.3 Photon polarization for various values of FP 27 3.1 Impulse approximation for O M C by 3 H e illustrated. The blobs represent the strong force that binds the nucleons. Please note that one has to sum all the contributions of the above process from every nucleon in order to get the complete rate 30 3.2 R M C photon spectra from two E P M calculations (one wi th the full Adler and Dothan amplitude and the other wi th only gauge invariant terms) and from I A calculations using various model potentials. A l l wavefunctions used have 22 Faddeev components and the permutation is projected on the same set of states. The infrared divergent part is not shown 43 3.3 Photon spectra from I A calculations for various values of gP. Wavefunction derived from Bonn-A potential is used for cal-culation 44 List of Figures vii 3.4 Same as figure (2.2) but for I A calculation. Using Bonn-A potential, r \u2122 ? ( \u00ab > 60 M e V ) = 0.1387 s\" 1 45 3.5 Plot of [1]\u00b0 vs. s for different nuclear potentials 51 3.6 P lo t of [a] 0 ' 1 vs. s 52 3.7 P lo t of [a] 2 ' 1 vs. 5 52 Acknowledgements viii Acknowledgements I would like to thank my supervisor Dr . Harold Fearing for his encouragement, guidance and support during my two years at T R I U M F . M y appreciation also goes to the kindness of Dr . David Measday of U B C \/ T R I U M F and Dr . Helmy Sherif of the University of Alber ta for being my second readers and giving me lots of constructive comments. This work would not have been possible without the tri-nucleon wavefunctions and instructions on how to use them given by Dr . Wolfgang Schadow. He, and Dr . Tae Sun Park, also provide me a lot of expert advice on setting up my own home L i n u x computer on which some of the present work is performed. Finally, I wish to thank Dr . Javed Iqbal for his introduction to me my present supervisor and both Michael Forbes and Alexander Busch for wri t ing two very nice WF\u00a3fi.2e thesis templates. I have chosen to use the one writ ten by the former because I received it a few days earlier than the latter one. To study in spring is treason; And summer is sleep's best reason; If winter hurries the fall, Then stop till next spring season. A Chinese poem for bibliophobes. Poet unknown. English translation by Yutang Lin [1]. ix 1. Introduction 1 1 Introduction 1.1 What and why of radiative muon capture by 3 He? 1,1.1 What? Radiative muon capture ( R M C ) by 3 H e is the process 1 and the corresponding non-radiative process is called the ordinary muon capture ( O M C ) by 3 H e . Here \/ i ~ stands for the muon that is captured; and 7 are respectively the muon neutrino and photon emitted. The capture takes place from a muonic atomic orbital , similar to an electronic state except that it is about 200 times smaller. Whi le it is in this orbit the muon can decay in the normal way with rate of 0.46 x 10 6 s _ 1 . Thus, as one wi l l see later, that muon capture is a fairly rare occurrence. Both capture processes (1.1) and (1.2) can be studied in at least two levels: the level of the nucleus and the level of the nucleons2. A t the nucleus level, one takes 3 H e and 3 H as whole entities and calculates the capture rate based on phenomenological nuclear form factors. This method is called the elementary particle model ( E P M ) and wi l l be discussed in chapter 2. One way to visualize processes (1.1) and (1.2) at the nucleon level is to regard the capture taking place on the constituent nucleons. B y recognizing 1 The \"breakup\" reactions: 3 H e + ^ ~ -\u00bb d + n + v^ + 7 and 3 H e + \/i~ -> p+n + n + 1^ + 7 are not considered here. 2 not to mention the level of quarks! 3 H e + pT -> 3 H + Vf, + 7 (1.1) 3 H e + p~ - r > 3 H + i\/p (1.2) 1. Introduction 2 that 3 H e ( 3 H) is some antisymmetric combination of | ppn) (| pnn)) and its permutations, one sums the following reaction that happens inside the nucleus to obtain a first order approximation to the capture rate: P + H~ ->\u2022 n + + 7 (1.3) for R M C and P + H~ -> n + Vp (1.4) for O M C This method is called the impulse approximation (IA) and w i l l be dis-cussed in chapter 3. Notice that since one is summing the contributions of the rate from constituent nucleons, the interactions between the nucleons are ignored by this method 3 . 1.1.2 Why? There are several motivations to calculate R M C (and O M C ) rates. Below are some major ones: 1. Experiment E592 at T R I U M F [2, 3] is an experiment to measure the photon spectrum of the radiative muon capture by 3 H e . A n accurate theoretical prediction seems necessary to interpret the experimental results. 2. W i t h the advent of computing hardware and software and methods of trinucleon wavefunction calculation, a more precise impulse approxi-mation seems overdue. Indeed, the most recent (and the only other) impulse approximation calculation was done about twenty years ago by K l i eb and Rood [4, 5] yet the Hamil tonian they were using d id not seem to completely satisfy some low energy theorems and they also handled some proton momentum terms in a crude way. The present impulse approximation calculation w i l l hopefully straighten out some of the errors and improve the approximation they made. 3 Please see section 3.1 for justification on ignoring the interactions between the nucle-ons. 1. Introduction 3 1.2 From ordinary muon capture to radiative muon capture Since there are many similarities between radiative muon capture and ordi-nary muon capture 4 , it is worthwhile to take a step back and have a glimpse at the simpler problem of ordinary muon capture. The relationship between different form factors of O M C derived from various hypotheses wi l l eventually be applied to the problem of R M C wi th some modifications. 1.2.1 O rd ina ry M u o n Capture In ordinary muon capture , the leading order Feynman diagram is shown in figure (1.1) F i g . 1.1: The Feynman diagram for ordinary muon capture. In figure (1.1), Pi stands for the four momentum of in i t ia l nucleus or nucleon (i.e. 3 H e for the E P M or proton for the IA) and Pf is the same quantity for the final nucleus or nucleon (i.e. 3 H for the E P M or neutron for the IA) ; \\i is the four momentum of the muon and v is that of the neutrino produced. The above diagram follows from the Fermi's original current-current coupling of the weak interaction. M O m C = ^ ^ ^ ^ - ^ d r o o i c (1-5) 4 See Mukhopadhyay [6] and Measday [7] for reviews on the topic of ordinary muon capture 1. Introduction 4 where the leptonic current J k p t o m c [ s jleptonic = u{u){la{l - ^)}U^) (1.6) and the hadronic current ^hadronic * s hadronic = *(Pf)Wu{PJ ( 1 7 ) The Wa is the weak hadronic vertex, which can be parameterized by four form factors. Wa = G ^ + G M ^ ^ + G ^ V + G P T 5 \u2014 (1.8) 2,Mn m where 75 = ^ V i V (1.9) A l l the GiS (i = V, M, A, P) are functions of Q2, the square of the momentum transfer at the hadronic vertex. 5 Note that Mn denotes the mass of the nucleus in E P M but the mass of nucleon in I A 6 . Relationship between various form factors The Isotriplet Vector Current Hypothesis ( I V C ) 7 stipulates that Gy(Q2) should be of the form, GV{Q2) = efGfv(Q2) - e^Q2) (1.10) and GM(Q2) GM(Q2) = HfGfM(Q2) - HiG^iQ2) (1.11) 5 Indeed, the form factors also carry the Pf (Pj) dependence if the initial (final) hadron is off-shell, but this dependence wil l , as usual, be neglected. 6 Please see appendix D for definition of m and for all other expressions that are not denned. 7 I V C is a statement saying that the vector current Jia = 4'W^'0'(Q)vector7\"**, * = +,0, - , $ = (| 3 H e or p), | 3 H or n))T, forms an isovector. Since r ' already transforms like a vector during isospin rotation, all that is required is W , Q(Q)Vector forms an isoscalar. The magnitude of it can be found for the i = 0 case. 1. Introduction 5 where Gy (Gy) is the electric form factor of the in i t ia l (final) particle and G'M (GM) is the magnetic form factor of the in i t ia l (final) particle. Notations of ei etc stand for charge and anomalous magnetic moment respectively. These form factors are supposed to be found experimentally from electron scattering experiments on the respective particles. The Conserved Vector Current Hypothesis ( C V C ) states that the vector current is conserved. This statement can be writ ten in equation: Q \u2022 ^hadronic = 0 (1.12) where Jhadronk ^S the vector part of the hadronic current. Notice that equation (1.12) is automatically satisfied (up to a small isospin breaking) given the form of Wa in equation (1.8). The Partial Conservation of Axial Current Hypothesis ( P C A C ) relates the divergence of the axial vector current to the pion field. A pictorial rep-resentation is shown in figure (1.2). Fig. 1.2: The Par t ia l Conservation of A x i a l Current Hypothesis Algebraically, it i s 8 : Q \u2022 Jhadronic = \u00abir\"4<\/ I = a \u00bb m \u00bb T O 2 l n 2 < \/ l * r ( Q ) l \u00bb > ( L 1 3 ) mTt W 8 It is a slight abuse of notation not to differentiate between Ja (a current operator) and ( \/ | Ja | i) (a current amplitude) but the context in which the above notations are used will usually clear the confusion. 1. Introduction 6 where 4>{Q) is the pion field which operates between the in i t ia l and final states and an is the pion decay constant (av = 0.9436). In the E P M , i is 3 H e and \/ is 3 H while i n the I A , i is the proton and \/ is the neutron. Note that (f\\JAQ)\\i) = -\u2014Gw{Qr')u{P,)>j*u{Pi) (1.14) vn \u201e\u2022 where Gn(Q2) is the pion-i-f coupling constant. Since, Q \u2022 ^hadronic = u(Pf){2Mn GA + GP^}^u{Pi) (1-15) equating (1.13) and (1.15) using (1.14) one wi l l obtain a relation between GA and Gp, which is: Gp = J ^ O * { l + e(Q 2 )} (1-16) Note that a\u201e disappears by virtue of the relation GA(0) = - o \u00bb G w ( 0 ) (1.18) This is just the P C A C at Q = 0. The e part is small and is approximately a constant over the range of Q2 concerned [4, 5]. Therefore the dominant part of Gp looks as if it were caused by a vir tual pion exchange (Figure (1.3)). Please do not confuse e wi th e, the latter wi l l be defined as the photon polarization four vector later in the chapter. 1.2.2 Radiative Muon Capture Let us shift the focus to radiative muon capture. There are several differences between O M C and R M C . These include: 1. Unlike the case of O M C where Q2 is a constant, Q2 for R M C varies approximately between \u2014 ra2 and m2. Moreover, there are two kinds of momentum transfer at the hadronic vertex: QL, signifying momentum transfer when the muon is radiating, and QH, the momentum transfer 1 . Introduction 7 Fig. 1.3: The virtual pion exchange diagram, a consequence of the P C A C when one of the hadrons is radiating. Specifically, QL = Pf-Pi = fl \u2014 K \u2014 V QH = (p} + n)-Pi . = p - v (1.19) where n, u and K are the four momenta of the muon, neutrino and photon respectively. 2. Gauge invariance (GI) . Let M = e \u2022 M be the R M C amplitude (e is the polarization vector of the emitted photon), gauge invariance means that K \u2022 M = 0 where \u00ab is the momentum of the photon. Gauge invariance essentially says that the longitudinal photon polarization is unphysical and would not affect physical observables. O f course, there is no such thing as gauge invariance in the case of O M C since the photon is absent i n the first place. 3. Using minimal substitution Qa -> QQ \u2014 eAa, C V C and P C A C now read (compare equations (1.12) and (1.13)): (Q - eA) \u2022 J h v a d r o n i c , R M C = 0 (1.20) (Q \u2014 eA) \u2022 ^ h a d r o n i c . R M C = (1-21) Note that AaJa l i terally means \"sticking a photon on the vertex of current J Q \" . 1. Introduction 8 Our aim is to construct an R M C amplitude that is gauge invariant and .satisfies both the C V C and P C A C . This method of constructing the R M C amplitude, to be shown below, is due to Adler and Dothan [8]. The method essentially checks the P C A C , C V C and G I of the currents wi th a photon attached to either the muon, proton ( 3 He), neutron ( 3 H) or the pion and adds in counterterms so that the sum of al l the currents (plus counterterms) satisfies G I , C V C and P C A C . Firs t consider the four external radiating diagrams in figure (1.4) and let the sum of their hadronic currents be ^hadronic.RMC-F i g . 1.4: The external radiating diagrams Let Mi be the value of the diagram on the upper left hand side of the figure. Enumerating the diagrams clockwise, one has, after suppressing the constant ^ \u00a7 V u d for convenience: M x = \u00ab H 7 Q ( l - 7 5 ) \u00ab ( ^ ( ^ \/ ) ^ Q ( Q H ) 5 F ( \/ P i -M 2 = u(u)la(l-15)u(fi)u(Pf)QfSF(\/Pf+ fi)Wa(QHXPi) : ^ h a ( l - 7 > ( ^ ( P \/ ) { m , _ ( Q V _ K ) 2 1. Introduction 9 {-ie){2<^ - *) \u2022 eGft^uiPi) M4 = u(v)7a(l-75)SF(li- jk)(-ie 4)u{n) u{Pf)Wa{QL)u{P^ (1.22) where SF is the Feynman propagator for spin | particles and Qi(j) = ie-i(f) ~p-aXpKpe\\, Hiy) being the anomalous magnetic moment of the in i t ia l (final) hadron. U p to one photon emission, (QQ \u2014 eAa)^hadronic.RMC c a n be calculated schematically as shown in figure (1.5). F i g . 1.5: Covariant derivative of the external weak hadronic currents i l -lustrated. The Aa part, in one photon l imit , only affects the non-radiating weak hadronic current (i.e. the muon radiating d i -agram). 1. Introduction 10 The result of the calculation 9 is: (Q \u2014 eA) \u2022 <\/hadronic,RMC = fi(P,){G? JQH + G\" JQ^5 + G ^ L F S F ( ^ ~ A ) Q i u ( P 0 + *{P\/)QtSF{ff+ \u00a3){G\u00bb Q\" + G\" JQHJ5 + G \u00bb ^ V M P 0 -m 2 , - (C} w - K ) 2 ml-Qjj m u(Pf){Gv jQL + 2MnGLA75 + G ^ 7 5 M P ) -G ^ 5 i l ! 7 R } \u00ab ( P 0 + fi(P\/){G^^^i^K\/,eQ}u(PO (1.23) where e has been defined in equation (1.17). Equat ion (1.23) can be reduced to the following fo rm 1 0 : (Q \u2014 eA) \u2022 ^ hadronic.RMC = s (P\/K (G? - Gv) > + (G? - G f r \/ 7 5 } \u00ab ( P i ) -u(Pf){GLMia^^ea + G L P ^ ^ M P i ) + u{Pf){2MnGHA+GHP^-}^SF{Pi- \u00a3)Qiu(Pi) + u(Pf)QfSF(\/Pf+ \/ ; ) { 2 M n G 2 + G | ^ } 7 5 u ( P , ) -m m 2 - (Q\u00ab - K ) 2 m \\ - Q \\ m u(pf){2MnGLA + G ^ y y p , - ) + m fi(i'\/){G?f^\u00ab^\/e\/,ett}tt(Pi) (1.24) 9 O n e calculates (Q - eA) \u2022 J h a d r o n i c . R M C i n o r d e r t o find t h e counterterms needed to satisfy equations (1.20) and (1.21). 1 0 T h e photon normal izat ion factor ^ k V w i U he dropped for convenience as it wi l l always appear in the phase space factor when one calculates the capture rate. 1. Introduction 11 The two tricks involved in the above calculation are the observations that {ft- fi)SF{P>- A) = (A- A)-B , M ft- fi - Mn = i + MnSF{Pi- fi) SF{Pf+ fi)(Pf+ fi) = i + MnSF{Pf+ fi) (1.25) and the Dirac equation for the in i t ia l and final hadrons. Notice that the second last term of equation (1.24) has already been worked out for the case of ordinary muon capture. In other words, equations (1.20) and (1.21) imply that, u(P,){(G\u00a3 -G$)t + {GHA - GLA) tf}u{Pi) - u{Pf){GLMia^-^-eQ +GpQ\u2014i^}u{Pi) + u{Pf){2MNG% + GP'C^h5SF{Pi- fi)Qiu{Pi) + m m s(P,)Q,SF(f\/+ mM\u201eG\u00abA + G g ^ h 5 \u00ab ( P , ) - ~ \/ l K e ) 2 (Q - eA) \u2022 A J h a d r o n i c , R M c = o 7 r m 3 T ( \/7 | {Q) \\ i) (1.26) provided the relationship between GA and G p is enforced by equation (1.16). In the above equation A J h a d r o n i c R M C is the extra piece of R M C hadronic current whose terms are to be derived from the Adler and Dothan procedure. {fl I 4>{Q) I i) is related to the pion photoproduction amplitude. The Born terms of {fy | {Q) | i) are given by: <\/7 I .1 \u00ab\u2022> = u ( P \/ ) { - ^ G f 7 5 } ^ ( ^ - fi)Qiu{Pi) + u{Pf)QfSF{Pf+ \/ , ) { - ^ L G f 7 5 H P ) + rn v \u2122 n \u201e L -i {-ie{2QL - \u00ab) ) \u2022 e 7 5 } ^ ( P i ) -u{Pf){-2-^G>nL{2Q\u00bb \u2022 e ) 7 5 } \" ( ^ ) (1-27) 1. Introduction 1 2 In the above equation, G'\u00a3 = dG^qT^\\Q^=Q2L- The first three terms correspond to the three external radiation diagrams and the last term is the gauge term to ensure the amplitude is gauge invariant when the pion is on shell. A look at equation (1.26) will reveal that the third, fourth and fifth terms almost fit the first, second and third term of the above expression respec-tively by virtue of equation (1.16) yet there are still two problematic spots: missing a -u(Ps)^{2Q\" .\u00a3)7\u00bb\u00ab(P<) ~ - u { P I ) ^ ^ ^ ^ { 2 Q \u00ab \u2022 O r y * ) = - u ( P \/ ) ^ ^ P ^ ( 2 Q f f \u2022 O T M P O = + 2MnG'iS\"'K)){2QH \u2022 e)75u(Pi) + 0(e2) term in the fifth term of equation (1.26) and the mismatching of momentum transfer between G\u201e and Gp there. The first problem can be solved by some rearrangement of terms in equation (1.26) and for the second one, one can expand G% as G% = G% + G'?{Q\\ - Q2H) + 0(Q4). After all the work, equation (1.26) can be written as11, u(Pf){(G$ - Gy) 4 + (GHA - GLA) muiPt) - u(Pf){GLMio\"e^eQ _ G \u00a3 ^ L \u00a3 7 5 W p . ) + u{Pf){2MnGHA + G^WSrift- A)QMP) + m m ml - (QH - K)2 ml-Q2H rn u(Pf)2MnGt{2%'K\\2Q\u00bb \u2022 e)^u(Pt) - u(Pf){\u00a3-^(2Q\" \u2022 e)75}u(P) + \u00ab ( P \/ ) { G ^ ^ ^ \u00ab a ^ ^ e a } \u00ab ( P 0 + {Q- eA) \u2022 AJhadT0IiiCtRMc = a\u201eml(f1\\4>{Q)\\i) (1.28) and (\/7 | jn(Q) | i) as <\/7 I J\u00bb(Q) I 0 = fi(P\/){-\u2014G^Sr'Pi- \u00a3)QAPi) + 771 TT u(Pf)QfSF(\/Pf+ m ~ G ^ } u { P i ) . + TTl-jr 1 1 Notice the gauge condition of the photon is chosen to be such that K-e = 0; therefore, Q L c = QH \u2022 e. 1. Introduction 13 ( - i e ( 2 Q L - K ) ) - e 7 5 } u ( P l ) -u{P{){-2-MllG'^2QH \u2022 6) 7 5}tt(P.-) + 7 8Mfi) (1.29) Assuming the linearity of Gn(Q2) over the range of Q2 concerned, the last two terms of equation (1.29) can be combined to form That is to say the C V C and P C A C conditions for R M C are: u(P,){(G\"v-0Lv) it + (G\u00ab-GLA) ^ W P J - i l P j K G i i ^ l - e , m mv - LJL u(Pf){\u00a3-^t(2QH \u2022 6 ) 7 5 } \u00ab ( P i ) + ^ Z ) ^ ! ^ ^ ' 4 ^ } ^ ) + (Q \u2014 e^4) \u2022 A Jhadronic.RMC provided equation (1.16) holds true. Now turn to gauge invariance, let Mi = e \u2022 Mj and if one then calculates K \u2022 M{ the results of calculation would be: K \u2022 Mi = e i u i ^ i l - ^ u ^ u i P ^ W ^ H P i ) K-M2 = -efu(u)7a(l-75)u(^u(Pf)Wa(QH)u(Pi) K-M3 = u(v)ya(l - y5)u(n)u(Pf){-\u2014% l m l - ( Q H - K y (-ie)(2QH \u2022 K)GHP}^u{Pi) m K-M4 = -u{v)la(l - 7 5 ) ^ ) u ( P \/ ) W Q ( < 2 I > ( P i ) (1.31) The tricks used in the above calculation involve the following: 1. Introduction 14 1. Using the antisymmetry of al3x wi th respect to switching indices to conclude that a\u00aeXKpK\\ = 0. 2. Using equations like Pi fi + Mn fi MP) -2Pi \u2022 K 2Pi \u2022 K- fi Pt + Mn fi u{Pi) -2Pi \u2022 K = u(Pi) (1.32) to simplify expressions. Since ei \u2014 ej = e, i t is easy to show that, \"t^Mi = u(u)ja(l - y5)u(n)u(Pf){(Gy - Gy)ya + i= l (GZQg - GLMQ%) 2mMnG'A(l + s) 1 2Mn a + - \u00b0 A h 7 + ml-Ql 2 9 ^ Q H ^ + ^Vw*) a-33) To satisfy gauge invariance up to O ( K 0 ) , one needs A M such that 4 K-iJ^Mi + AM) = K \u2022 M i = l = 0 + O(K2) (1.34) Postulate12 A4 Q adronic,RMC = -2u(Pf){G'vl\u00bb + G\\1rf}{g\u00bbaQH-e + Kae>>-K\u00bbe\u00ab}u(Pi) - u(Pf){GLMiaa^ + G LXj5}u(Pt) 2,Mn m 1 2 K l i e b and Rood did not have - u ( P f ) { 2 m \" f ^ } + e ) 2 ^ e Q a H i 5 + G'v(uf -Hi)^-i) (1-36) Putting it back into the left hand side of equation (1.30), one would get for the left hand side: -u(P,){2MnG'A(2QH \u2022 < r ) 7 5 } \u00ab \u00ab ) + u(P,){GLv^^-i^< ^^(P,) + 2 ^ ( 2 Q \" . . - ( 1 + \u00a3 ) Q ^ ) . _ \" rnl-Ql \u00a3-^(2QH \u2022 e ) 7 5 W ) (1.37) m _ 1 3 Equation (1.30) determines terms up to 0(Q\u00b0KN) and equation (1.34) determines terms up to O(K0Qm). Here K\u00b0 means \"photon four momentum up to order zero\" and n, m are arbitrary positive integers. 1. Introduction 16 which equals --^Ln(Pf)2MnG'A(2Q\u00bb \u2022 \u20ac) 7 5u(P t) + ^PfHGy^^ ^ K ^ M P O - u(PM\u00a32\u2122Mn^QH'eQW + \u2022 eh5}u(Pi) m~ \u2014 Q L m m (1.38) Using to eliminate G'A in favour of G'w and e in expression (1.38) one would end up having, ^ W } \u00ab ( P i ) + C(c: 2) (1.40) The first term is just the right hand side of equation (1.30). The magnitude of the second term can be approximated as follows: ^ < \u00b0 > < ^ ' + <&fc> ~ ^ ( 0 ) ( ^ 5 ) 2 + ^ ) ) < 3 2 ~ 2KG\"Q4 (1.41) keeping i n mind that a* ~ 1 (0.9436, to be exact); in the above equation, G\", \\G\"\\ \u00ab 1 is the small second derivative of either GA or G\u201e and K is some constant with magnitude \\K\\ ~ 1. The third term is annoying but the Adler and Dothan procedure is not sufficient to determine uniquely counterterms to get r id of it. To see this point, notice that a term in A Jh adronic, R M C W ^ the format -u(P^GLvYa^=^za^KCe,}u(Pi) (1.42) 1. Introduction 17 with YA = -^-^ or (and probably among many others) can both do the trick but both have some undesirable pole behaviours. Notice that the latter is formally a term of O(KQ) but it is exactly these terms that the Adler and Dothan procedure cannot determine uniquely. A more detailed study of these structure terms should probably involve a model that incorporates the exchange currents of mesons. Please see section 3.7.2 for possible implications of these undeterminable terms on the resulting capture rates. O(KQ) terms aside, the Adler and Dothan procedure does indeed generate terms of A J h a d r o n i c as in equation (1.35) that satisfy the gauge invariant requirement up to O{K0) and P C A C up to terms of 0{eG\") and O{Q0). 1.3 Summary After a brief introduction of the process of ordinary muon capture by 3 H e and its various hypotheses involving the O M C hadronic current, the correspond-ing radiative process is compared against the ordinary one. For the R M C amplitude to be gauge invariant and its hadronic currents to satisfy the C V C and the P C A C , an extra current A J \u00a3 a d r o n i c R M C must be added in addition to a naive insertion of photons to external particles. The Adler and Dothan procedure is employed to generate those extra terms as far as possible. This procedure guarantees the uniqueness of those extra terms up to O(K0) by gauge invariant requirement and O(Q0) by the C V C and the P C A C except the pion pole terms which are treated exactly. Terms of O(KQ), however, are not determined uniquely by the procedure and these terms can only be revealed by a detailed examination of the structure of the weak hadronic vertex. 2. The elementary par t ic le model 18 2. The elementary particle model 2.1 Introduction The elementary particle model is probably the simplest method which can be used i n the calculation of the R M C rate. It was first used by K i m and Primakoff [9, 10] in calculating the beta decay of complex nuclei. A s its name suggests, the central idea of the method is to completely ignore the internal structure of both the 3 H e and 3 H and regard them as \"elementary particles\". Since both 3 H e and 3 H have spin | , and isospin \\ up to a small isospin breaking, the calculation is relatively simple compared wi th other methods.. The drawback is one sti l l has to rely on experiments done on these nuclei for the phenomenological form factors parameterizing the reaction. For example, one can perform electron scattering experiments for determining the vector current form factor and beta decay experiments for the axial current form factor. Previous E P M calculations of R M C by 3 H e include Fearing [11], K l i e b and R o o d [4, 5] and Hwang and Primakoff [12] but only the calculation of Kl ieb and R o o d included some Adler and Dothan terms. Moreover, some physical quantities were obtained by Kl ieb and Rood v ia a non-relativistically reduced amplitude. The present calculation differs from those previous ones in that the full Adler and Dothan terms are included and the resultant amplitude is not non-relativistically approximated in the hope of getting a more accurate reaction rate. 2.2 Kinematics and nuclear form factors For the elementary particle model description of R M C by 3 H e , the degrees of freedom are the four momenta of 3 H e ( P \u00bb H E ) , 3 H (P>H), photon (K), neu-trino (v) and muon (\/J,), together wi th their respective spinors. The four-2. The elementary particle model 19 momentum conservation is written as: PSH + V + K = P3He + P (2.1) The quantity that is of interest is the differential capture rate, which is given b y 1 dT = E y d ( P W 5 p a c e ) T r ( p M t M ) ( 2 2 ) photon polarization for R M C , and the total capture rate Tomc = j d(phase space)^ Tr ( p i k T ^ M 0 \u2122 ) (2.3) for O M C . Here M (Momc) is the R M C ( O M C ) amplitude that was discussed i n chapter 1. One just needs to recognize the following name changes u(Pi) -> u(Psne) u{Pf) -> u{Pm) Mn \u2014> Mt on the hadronic vertex d -\u00bb\u2022 Fu i = V, M, A, P (2.4) for M (denote it as M(Psiil P ^ j e , Mt)). The name change from G \u2014> F is merely a change in accordance with the usual convention. p is the density matrix which describes the ini t ia l spin configurations of the muonic atom. The detailed implementation of p w i l l be discussed in the next section. d(phase space) (d(phase space) o m c ) is the differential phase space factor of the emitted particles for R M C ( O M C ) . The detailed calculation of these wi l l appear in the next section. Also note the integral sign of equation (2.2) does not apply to K since what one needs is the differential rate 2 . 2.2.1 Nuclear form factors One can conveniently parameterize the nuclear electromagnetic form factors linearly in terms of Q2 when Q2 is small, as in the case of R M C by 3 H e . FyHe(Q2) = 1 + ^r 2)3He,electricC} 2 + 0 ( Q 4 ) (2.5) 0 P M H 6 ( Q 2 ) = 1 + ^ 2 ) 3 H e , m a g n e t i c Q 2 + 0{Q4) (2.6) 1 Despite its notation, d { P h a s ^ P a c e ) j s actually a differential of angular variables. 2 One integrates K in equation (2.2) to get the total capture rate of R M C . 2. T h e elementary par t ic le mode l 20 and analogous equations also hold true for the case of 3 H . The above equa-tions relate the electromagnetic form factors to (r 2 ) where ( r 2 )^ is some form of root mean square charge or magnetic radius of the nucleus which are sup-posed to be found from experiments (see, for example, refs. [13, 14, 15, 16] for related experiments). {r2)h ( f m ) He, electric 3 H e , magnetic 3 H , electric 3 H , magnetic 1.976 \u00b1 0 . 0 1 5 1.99 \u00b1 0 . 0 6 1.68 \u00b1 0 . 0 3 1.72 \u00b1 0 . 0 6 Tab. 2.1: A set of ( r 2 ) 2 values. The ( r 2 ) notation indicates that the above quantities are root mean square values of nuclear radii . Exper-imental data from Ottermann et al . [14] ( 3 He) and Beck et al [13] ( 3 H ) . B y the I V C discussed in chapter 1, one has Fv = 2FvHe-Fvn FM \u2014 fi3neF^e - \/ ^ H - F J 1 (2.7) where \/ U 3 H e is the anomalous magnetic moment of 3 He wi th a similar definition for 3 H . Their values are A^He \u2014 \u20148.3689 n.m \/ i 3 H = 7.9173 n.m (2.8) Equations (2.7) imply that Fv(0) = 1 FM(0) = \/ i 3 H e - \/ \/ 3 H (2.9) The functional form of Fv and FM are chosen to be Fy = l + UyQ2 0 FM = 0u3He - M\u00bb H )(1 + \\RMQ2) (210) 0 2. T h e e lementary par t ic le m o d e l 21 with Ry = 1.94 fm and RM = 1-72 fm. These values are taken from a grand average of various experiments and (with RA) they produce \"good\" numerical values for O M C . The function FA(Q2) is not so clear cut since there does not seem to have any direct experimental determination of that form factor. One has to resort to the impulse approximation to estimate the dependence of FA on Q2 over the range of Q2 concerned. Delorme [17], using impulse approximation, expressed the nuclear form factors in terms of nucleon form factors. This method is also employed by Kl i eb and Rood [4, 5] and Congleton and Fearing [18, 19] in their E P M calculations. Fy(Q2) = gv(Q2)[l}\u00b0 (2.11) FM(Q2) = 3{9v(Q2) + 9M(Q2)M- - gv(Q2)[i}\u00b0 -^ ( Q 2 ) ^ ] 1 ' 1 (2.12) FA(Q2) = gA(Q2)[a]~ (2.13) M2 FP(Q2) = VgP(Q2) + \u00a7V2^gA(Q2)[o]2l (2.14) where [ l]\u00b0 , [a]+, [a]~, [ff]2'1 and [ i P ] 1 , 1 are the impulse approximation reduced matr ix elements resulting from integration of internal momenta. They wi l l be defined in section 3.6.1. The Q2 dependence of FA is given by FA(Q2) = 9A(Q2M-(Q2) , . FA(0) 9A(0M-(0) { \u2022 } Notice that Fy(Q2) + FM(Q2) -> 3{gv(Q2) + gM(Q2)}[o}- (2.16) when [zP] 1 ' 1 \u2014>\u2022 0. Since [ i P ] 1 ' 1 is indeed small, one can use the above l imit ing equation as an approximation and put it back i n equation (2.15). The result is: FA{Q2) = 1-Q2\/M2 2 FV(Q2) + FM(Q2)} ( 2 U ) FA(0) 1 - Q2\/M\\ Fv(0) + FM(0) * { - ' where MV = 0.710 G e V 2 and M\\ = 1 . 0 8 \u00b1 0 . 0 4 G e V 2 [18]. Th i s follows from the expansion g{ = ^ ( 0 ) \/ ( l - Q2\/M2)2 i = V, A. Upon expansion to 0(Q2), 2. T h e elementary par t ic le model 22 FA can be written as FA = F A (0 ) (1 +jU*Q2) (2.18) FA(0) = 1.212 \u00b1 0 . 0 0 4 [18] (2.19) w h e r e R\\ = Fv%ttZ{(o)Rlf + KBI - -k) w h i c h m a k e s ^ = i - 7 0 3 f m -The P C A C hypothesis discussed in chapter 1 relates FP to FA- The relationship is shown as equation (1.16) with GA replaced by FA and Gp replaced by Fp. e 2.3 Results When al l the relevant Feynman diagrams are evaluated and the gauge free-dom is exploited so that the photon is made transverse (i.e. K \u2022 e \u2014 0 and 0 = 0), the relativistic amplitude M (M = \u00a3 4 = i + J l e p t o n i c - A J h a d r o n i c , R M C in equation (1.35)) can be written in the following non-relativistic like form, which operates on the product space (see ref. [20], chapters 3 and 4 for in-formation on direct products of matrices) of the hadronic spin (via a) and leptonic spin (via di): M(P^P^Mt) = W^Vud(l-c?L-0){f1oL-ex + f2a-ex + ifz^x x &i \u2022 \u00b0 + -exa- s + -^-v \u2022 ex + 2m 2m . fl -. - _ , \/8 _ \/9 -. , i\u2014s x eA \u2022 cr + \u2014OL \u2022 so \u2022 eA + \u2014 aL \u2022 au \u2022 ex + 2m 2m 2m \/io -. \u2014 _ . . \/ n - s - - . -o \u2022 sv \u2022 ex + i^z^s x v \u2022 ooL \u2022 ex + 4 m 2 4 m 2 ' \/12 _ \u2014 _ , \/ l 3 _ -a \u2022 vv \u2022 ex + \u2014d \u2022 VoL \u2022 ex + Am2\" \" A ' 2m ' i-p^sx ex-oaL \u2022 s + ip^sx v \u2022 ov-ex} (2.20) 4mz 8 m J where S = U + K and N ' is the normalization factor for the particle spinors, 1 N ' = P^ + Mt 2\\ 2Pb (2.21) 2. T h e elementary par t ic le mode l 23 The coefficient of | in N' comes from the normalizations of the photon and the neutrino. The photon polarization unit vector A = 1,-1 is defined as, k x e A = -iXnex ( 2 . 2 2 ) That means the Sphoton polarization m equation ( 2 . 2 ) should be replaced by S A = - I , + I - \/t) z = 1 \u2014 15 are some very complicated scalar functions independent of a, CL and t\\ and will not be denoted. , Define d(phase space) to be the differential phase space factor which has the form \u201e , \\ i \/ \/r\\\\ 1 9 d3p3H d3u 1 d3k d(phase space) = C \\ 4>^{0) (2TT ) 3 (2TT) 3 \u00ab ( 2 ? r ) 3 (2TT)4~~K2dKdk ( 2 . 2 4 ) where cos(f?) = 0 \u2022 k and - Mm + Mt) - m2 - 2mMt [ K ' \u00b0> \" 2\u00ab(1 - cos(0)) - 2(m + M t ) l ^ 5 j expressing the four-momentum conservation of the 5 function. Define p^ - as the elements of the density matrix which takes care of the initial spin of the muonic atom. Since p is a density matrix, it has the property Tr(p) = 1. p is chosen to be diagonal in the | \/ , fz) (coupled lepton (s),) spin and nuclear (sn u ci e a r) spin state; f = s^ + s n u ci e a r) basis rather than the | s ,^ 5 n u c iea r ; 5^, s\u201e u c l e a r) because of the hyperfine splitting of the muon in the muonic IS orbit. In the special case when the initial spin mixtures 2. T h e elementary par t ic le mode l 24 are equal, p{j = \u00b16^ and is basis independent. In that case, equation (2.2) reduces to dTstat _ s-^ f fd(phase space) dK ~ x=k+J dK iTr(M(PsH, P^ e, M^M(Pm, Pane, Mt))} (2.26) The quantity d I j ^ g t is commonly called the differential statistical rate and the corresponding photon spectrum is plotted in figure (2.1). The total R M C rate is obtained by \/ ^dn. The E P M description of R M C statistical rate for photon momentum higher than 5 MeV 3 is T^(K > 5 MeV) = 0.8189 s _ 1 (2.27) 60 MeV) = 0.2113 s _ 1 (2.28) Figure (2.3) shows different photon polarization 5 P 7(K) predictions vs. FP 2.3.1 Ord ina ry M u o n Captu re The calculation of ordinary muon capture (OMC) is also done for comparison. Since R M C is the major topic of this work, the presentation below will be brief. Let M0\"10 be the amplitude for the ordinary muon capture as shown in section 1.2.1. The momentum transfer squared at the hadronic vertex 3 The K integration cannot go to zero because of infrared divergence. 4 The two step process fj,~ \u2014> e~ +ue+vli fo l lowedbye~+ 3 He \u2014> e - + 3 H e + 7 produces photons up to 53 M e V . 5 To obtain \u00a3P 7 (\u00ab), simply replace Ephoton polarization i n equation (2.2) by EA=- 1 I + 1(-I) ( A + 1 ) \/ 2-2. T h e elementary par t ic le mode l 25 40 60 photon momentum [MeV] Fig. 2.1: The E P M description RMC photon spectra using the full Adler Dothan amplitude for various values of FP. The \"Klieb and Rood\" values are taken from the relativistic calculation of ref. [5] which are not shown in ref. [4] Q2 = -(!\u00a7S^)2 = -0-2763 fm\" 2 . Upon expanding M^P^, Pme, Mt) and rewriting in a non-relativistic like format, one has: M \u00b0 - c ( P 3 H , P 3 H e , M t ) = K'^Vud(l-oL-v){gl+g2a-v + g3aL-o} where gx = Fv(l + v M t + pOj mv M 2 M t ( M t + P \u00b0 I ) 92 \u2014 (Fp \u2014 FA \u2014 Fy \u2014 FM) M + P4 #3 FA - {Fv + FM) v Mt + Pb (2 .29) (2 .30) (2 .31) (2 .32) 2. The elementary particle model 26 Fig. 2.2: Sensitivity of r\u2122(K > 60 M e V ) using E P M wi th respect to FP N ' = 1 M t + \\ 2 2 P & The differential phase space factor for O M C d(phase space)^ is: (2.33) d(phase space) c Ci^(\u00b0)|2f7r7 ^ ^ ^ ^ + \/0 = C|^(0)| 1 P\u00b0 V 2 (2TT) 2 m + M t where i \/ is evaluated at, due to the 5 function m 2 + 2 m M t v -(2.34) (2.35) 2(m + M t ) The value of P ^ for E P M description of O M C is therefore: P\\ = s]M? + v2 (2.36) using the value of v from equation (2.35). The statistical rate for O M C is obtained via equation (2.3) with pij = \\8ij. The result, wi th FP set at the P C A C value, is: TZt = 1503 s- 1 (2.37) 2. T h e e lementary par t ic le m o d e l 27 i * photon momentum [MeV] F i g . 2 .3 : Photon polarization for various values of Fp. with F y ( - 0 . 2 7 6 3 fm- 2 ) = 0.8267 F M ( - 0 . 2 7 6 3 f m - 2 ) = -14.067 F A ( - 0 . 2 7 6 3 fm~ 2) = 1.050 (2.38) which should be compared wi th the most recent calculation by Congleton and Fearing [19] who obtained a value of 1497 s _ 1 . The closeness of the above result to theirs is not so surprising given that the values of the formfactors (i.e. FV, FM and FA ) at the Q2 for O M C are approximately the same as the ones they were using (which also gives a reason why a linear parameterization of formfactors in section 2.2.1 is sufficient). Also interesting to mention is the most recent experimental result of statistical O M C rate of 1496 \u00b1 4 s _ 1 [21]. 2.3.2 Discussion Below are some points to note wi th regards to the E P M description of R M C . Please refer to section 3.7.2 discussions on the comparison between E P M and I A . 2. The elementary par t ic le mode l 28 1. The photon spectrum (figure (2.1)) and total R M C rate agree wi th the results of K l i eb and Rood [4, 5] who obtained a total rate of 0.814 s _ 1 via a non-relativistically approximated amplitude (also see figure (2.2) for variation of R M C statistical rate vs. FP) , showing that the extra terms not included by Kl i eb and R o o d do not contribute significantly to the photon spectrum. 2. The photon polarization P 7 ( K ) obtained in this work does differ from that of Kl i eb and Rood quite significantly (see figure 2.3). The reason behind that is while a non-relativistically reduced amplitude used by Kl i eb and Rood produces the correct spectrum within a few percent, it cannot produce P 7 ( K ) accurately. Fearing [22] noted that the first order contribution of P 7 ( K ) actually comes from 0(-^) terms in the squared Hamil tonian but K l i eb and Rood compromised P 7 ( \u00ab ) ' s accuracy by truncating many 0(-^) terms when they squared the already non-relativistically reduced amplitude (see section 3.7.2). 3. The Adler and Dothan terms have a great influence on the resulting photon spectrum. They contribute about 10% of the spectrum at the high photon energy range (K > 60 M e V ) . The large size of G'v and G'A because of the large root mean square nuclear radii of 3 H e and 3 H compared to proton or neutron is pr imari ly responsible for the contri-bution. 2.4 Summary In the framework of the E P M description of R M C by 3 H e , the \" P C A C \" value of the exclusive statistical capture rate, TRS^(K > 5 M e V ) , v ia a fully Adler and Dothan amplitude, is 0.8189 s _ 1 and the O M C statistical rate is 1503 s _ 1 . Whi l e these results agree wi th K l i e b and R o o d [4, 5] and Congleton and Fearing [18, 19] respectively, K l i e b and Rood's photon polarization P 7 ( K ) disagrees wi th the result,of this work, possibly due to the fact that they threw away some terms in the squared amplitude which contribute significantly to P 7 ( \u00ab ) . 3. T h e impulse approx ima t ion 2 9 3. The impulse approximation 3.1 Introduction In this chapter, the method of the impulse approximation w i l l be used to investigate the R M C by 3 H e . Even though the elementary particle model is a convenient approach to the reaction, it does not give the detailed mecha-nism of \"what is happening inside\". Indeed, each of the 3 H e and 3 H nuclei is a three-nucleon bound state. Perhaps it would be a good idea to look at the nuclear reaction using a nucleon perspective. The idea is, using a t r in-ucleon wavefunction, to integrate out al l the degrees of freedom relevant to the internal coordinates of the trinucleon system and to extract quantities that are dependent on the E P M degrees of freedom. In other words, one wants an expression similar to M(P3ne, P S H , Mt) in the E P M but one with al l the information regarding the internal degrees of freedom hidden in some formfactors resulting from internal integration 1 . The first step is to regard the three bound state nucleons as free nucleons and the Hamil tonian of the nuclear reaction as the sum of the Hamil tonian of its nucleon counterparts. This procedure is called the impulse approxima-tion. Naively speaking, the approximation of bound state nucleons by free nucleons in this case is a good one since the binding energy (~ 8 M e V ) is about 0.3% of the mass of the nucleus. However, the impulse approxima-tion deliberately neglects the contribution to the R M C from meson exchange currents between the nucleons, which may be large. Figure (3.1) shows the impulse approximation for O M C by 3 H e . 1 The method of finding this expression will be discussed in section 3.6 3. T h e impulse app rox ima t ion 30 F i g . 3 .1 : Impulse approximation for O M C by 3 H e illustrated. The blobs represent the strong force that binds the nucleons. Please note that one has to sum al l the contributions of the above process from every nucleon in order to get the complete rate. 3.2 \"Internal\" and \"external\" degrees of freedom In describing any non relativistic three-particle system, there exists a coor-dinate transformation which separates the centre of mass coordinates and the internal coordinates of the three particles. This enables one to treat the three-body problem as an effective two-body problem. Let ki be the (three) momentum of the of the ith nucleon and define: P \u2014 ka + hp -t- fc7, Qa = \u2014 3^\/9 \u2014 3^7, (3-1) Pa = \\kp ~ where (a, B, 7) are any set of numbers obtained from the cyclic permutation of (1 ,2 ,3) . 2 2 From now on, when no confusion arises, g^and <\u00a3,will sometimes be denoted as q and 3. T h e impulse approx ima t ion 31 It is now obvious that P is the centre of mass momentum vector and qa is the momentum of particle a (spectator) with respect to the centre of mass momentum of the other two particles (subsystem), while pa is the momentum of particle \/? wi th respect to particle 7 (or vice versa). Also note that the Jacobian of the transformation I d^'kp'kS1} I is 1. 3.3 Wavefunction The tri-nucleon momentum space wavefunction, provided by Schadow [23, 24, 25], is obtained by solving the Faddeev equation (see, for example, ref. [26]) with various model potentials. It can be written as a sum of channel wavefunctions where each channel represents a specific angular momentum and spin eigenstate. \\^) = J2MP,) = | V'Fad) + P | ^Fad) (3 .4) where | ipvaA) is the \"Faddeev component\" of the wavefunction obtained by solving the wave equation wi th one particular choice of nucleon labelling scheme, and P is the necessary permutation operator to make the full wave-function antisymmetric. The normalization is defined as: (ib\\tb) = 3 { # a d | rb) = 1 (3 .5 ) The fact that the wavefunction is antisymmetric with respect to nucleon labellings makes the second expression of the above equation legitimate. It is also numerically possible since it consists of a finite number of terms. The first expression is, however, numerically impossible in general since permuting the nucleon label ( P | ^Fad)) w i l l project the Faddeev component onto a set of eigenstates which is usually much larger, i f not infinitely larger, than the original (finite) set of eigenstates comprising the Faddeev components. It is therefore necessary in practice to restrict the set of projected states to one wi th a finite number of elements. Suppose one defines Pl^Fad) to be the permutation projected on a restricted set of states, then the nu-merical normalization of the wavefunction (V'lV')num = ( (V 'Fad l+^Fad lP^d^+PIV 'Fad) ) (3-6) wi l l not be 1 because of this restriction. For example, a 8-8 channel wave-function (i.e. | i\/>Fad) has 8 channels and the permutation is projected on the same 8 channel) has a normalization of 0 .93 and a 2 2 - 2 2 channel wave-function typically has a numerical normalization around 0 . 99 (see appendix C ) . 3. T h e impulse app rox ima t ion 33 3.4 Hamiltonian The effective Hamil tonian is defined by: 3 (3.7) where Hi and T7 are the weak interaction Hamil tonian and the isospin low-ering operator acting on ith nucleon of the trinucleon system. The relevant Feynman diagrams are exact duplicates of the E P M ones except that the 3 H e is replaced by a proton and the 3 H is replaced by a neutron. Since the wavefunction is fully antisymmetric with respect to interchange of particle labels, one can write equation (3.7) as where j is 1 or 2 or 3, provided H is acting on the above wavefunction. The relation between different nucleon form factors is an exact duplicate of the nuclear ones. In other words, al l the C V C , P C A C , I V C hold for the nucleon case once one replaces 3 H e with a proton and 3 H with a neutron. One has, 3 ri = Y,TrHi = ZrjHj (3.8) i=i 3.5 Nucleon form factors 9i(Q2) = 9i(0)(l + lr2Q2) (3.9) where i = V,M,A and 9v{0) = 1 9M{0) = Hp ~ gA{0) = -1 .267 \u00b1 0 . 0 0 3 5 [27] (3.10) with r2v = 0.576 fm 2 , r2M = 0.771 fm 2 and r\\ = 0.433 f m 2 [18]. 3. The impulse approx ima t ion 34 3.6 IA E P M translation The a im of this section is to find the correspondence between the E P M amplitude and the I A amplitude. . Let Mia be the R M C transition amplitude obtained by the impulse ap-proximation. This means Mia is the I A equivalent of M(Pme, Pm, M t ) in the E P M . More explicitly, one would use the following formula in order to have the quantity ^ (compare wi th equation (2.2)): fK = t J ^ J ^ e space) ( 3 . n ) and = j d{phase 5pace) o m c Tr(pMr c t M^ c ) (3.12) for O M C . Define M.epm as 4: Mepm = (27T)45 (4)(P3H + (y + K - fj) - P3He)M(P3He, P^ M t ) (3.13) It then follows that the I A equivalent of M.epm, M.ia, would be: Mia = 3\/(2ir)W>(k'a + {y + K - \/i) - fcQ)(27r)4^(^ - ^)(2TT) 4^ 4>(A:; - fc7) ^tu{k'Q, le-, fc7)^3He(fcQ, kp, ky)u)(ka, kp, ky, fcjj, k\u00b0)M(k'a, ka, Mp) d4k' d4k-( 3 1 4 ) In words, the above equation merely says the capture takes place on (spec-tator, in the sense of equations (3.1)) nucleon a 5 and so one integrates out the degrees of freedom of the two subsystem nucleons. The 3 comes from the antisymmetrization of the wavefunction. u>(ka,k0,kj,kp,k\u00b0), loosely speaking, specifies the \"dispersion relation\" of the subsystem particles in the phase space. There is a degree of ar-bitrariness in choosing the functional form of ui since al l the nucleons are 4 The S function merely enforces the four momentum conservation of the reaction. 5 That is why M(k'a, ka,Mp) appears. 3. T h e impulse app rox ima t ion 35 off-shell and therefore the non-relativistic wavefunction cannot enforce the energy-momentum relation. Suppose one were to integrate 6 Mepm over the phase-space of and P^, that is, one does the following: \/ d4Pm d4Pme e P m ( 2 7 T ) 4 (27T)4 ( 3 - 1 5 ) The corresponding thing in I A is: (3.16) can be writ ten as: (3.16) = 3 J(2n)35{3){k'a- ka + v + it - fl)(2Tr)5(k'\u00b0 - k\u00b0a + u + K - ^) (2ir)3*~~~~$, - k0)(2Tr)6(k'\u00b0 - *J)(2,r)3*~~~~(*; - k,)(2v)5(k!\u00b0 - k\u00bb) ^n(k'a, k'p, k'JipiHeiha, h, ky)u(kay kp, ky, k\u00b0p, k\u00b0)M(k'a, ka, Mp) . 1 1 ^ (27T) 3 (27T) 3 ( 2 7 T ) ( 2 7 T ) ' Using the inverse of equation set (3.1) to transform variables from { \u00a3(0,0,7)} to {P,pa, qa} and from to + hp + one has: (3.16) = 3f(2n)W\\hp-P) + fa-?a + P+ii-fl) (2n)5(k'\u00b0-k\u00b0a + u + K-tJ,0) (27r)35(3)(i(P' - P ) _ I ( \u00a3 _ + (pfa _ pQ))(27r)a(A;,0 - fc\u00b0) ^ ^ ( J * - P) _ I ( \u00a3 _ &) _ _ P-Q))(27r)5(fc;0 - fcjj) ^ H ( ^ , % , * > H . ( * a , * \/ \u00bb , f e r ) \u00ab A f | || ^ ' ^ ' t y | d{P,pa,qQ) d{P\\pa,^) f ^ f ^ f ^ d r ^ d ^ f p ^ dkl_ (2;r)3 (2TT)3 (2TT)3 (2TT)3 (2TT)3 (2TT)3 Ja^ (2TT) 6 One does this integration to get a relativistically invariant quantity so that a direct comparison of the same physical quantity obtained by the IA and the E P M is possible. 3. T h e impulse approx ima t ion 36 Again , using the multivariable version of one of the basic properties of 5 function: 5(n)(x) = | det(T) | 8{n\\Tx) (3.19) where n is any positive integer and T is any non-singular linear transfor-mation and | det(T) | is the absolute value of the determinant of T . We have, (3.16) = 3 j\"(2TT)3<5(3)(P' - P + v + K - \/2)(2TT)<5(\u00a3' - E + v + K - \/j\u00b0) ^ T T ) 3 ^ 3 ^ - pa)(27r)3<5(3>(& - qa + 2-(u + K - p)) (2n)5(k'\u00bb - k\u00b0p)(2n)5(k'f - k^t^PV\u00ab,O^He(Ap\u00ab,9a) - - - d3P' d3 The assumption of impulse approximation is that the two subsystem nucleons are both on-shell; therefore, u>(ka, kp,ky,kp, k\u00ae) can be writ ten as: u(ka, fy, \u00a3 7 I Arj, fc?) = (2kp)(2k\u00b0)(2Tr)6(kj - M2)(2ir)5(k2 - M p 2 ) (3.21) Using the above functional form of ui and integrating equation (3.20), we have: (3.16) = 3l{2TT)3S{3\\P'-P + u + K-fi)(2Tr)S(E'-E + u + K-fi\u00b0) (2ir)36Wtfa - pa)(2Tr)35Wtfa - qQ + \\{V + K - fi) ^ H . ( ^ . P a . fa)1\u00bbH*(P,Pa, Qa)M(k'a, ka, Mp) ^ l ^ f ^ ^ l ^ o f ^ d E ^ d E _ (2TT)3 (2TT)3 (2TT) 3 (2TT)3 (2TT) 3 (2TT) 3 (2TT) (2TT) 1 ' It is tempting to associate the centre of mass momentum of the three nucleons P(P') of I A with P W B H ) of E P M and E = k\u00b0a + ifcjj + fc\u00b0 wi th P ^ e o f E P M and E' = k'\u00a3 + k'p+k'\u00ae wi th P ^ o f E P M ; however, there are s t i l l two problems which need to be resolved: 1. What is the dispersion relation of the struck nucleon? 3. T h e impulse approx ima t ion 37 2. Do (E,P) and (E',P') satisfy the \"on-shell\" condition? It turns out that if one puts the struck nucleon on-shell, the answer for ques-tion (2) is negative; one could have adjusted the energy-momentum relation of both the struck incident and final nucleon in a way such that (E,P) and (E',P') are on-shell but the latter approach has the disadvantage that the off-shell nucleon wi l l no longer satisfy the Dirac equation, the equation that has been used extensively to obtain the Adler and Dothan terms in the I A amplitude. Let us see the extent the on-shell condition is not satisfied by (E,P) and (E',P') if we put the struck nucleon on-shell. One wants to know whether P 2 = = Ml and P ' 2 = P4 = M 2 . Now, P = \"^\"^ ki \u2022 kj = 3 M 2 + \u00a3 ki-k = 3 M 2 + \u00a3 (M2 + \\(ki - fc,)2) + 0(fc 4 ) (3.23). Except for the difference in binding energy, the error in assuming that P 2 = M 2 is of the order of the square of the relative momenta of the in i t ia l nucleons. Even for the case of 3 H , the error is at most 0((j%-)2) which is about 1%. The approach that is taken is to regard al l the nucleons as free nucle-ons and assume (to a considerable accuracy) that (E, P) and (E',P') is a four-vector satisfying the on-shell particle criterion. Then, noting that ipme(P,fa,q'a) = V*He(&>&)> o n e c a n associate \u2022M(P3H,P3He ,M t ) o 3 | ( 2 7 r ) 3 ^ 3 > ( \u00a3 ^ V'\u00bbH(&\u00bbOV'\u00bbHe(Pa,\u00aba) The right hand side of the above expression is just Mia defined previously. We have successfully separated the centre of mass coordinates from the inter-nal coordinates of the I A amplitude. In order to actually do the calculation, one also needs to expand those 6 functions in (3.24) in terms of angular mo-mentum eigenstates. To this end, one has (after setting \/7 = 0 and denoting 3. T h e impulse approx imat ion 38 v + k = s): (2*)*8Mtfa-pa) = (2n)M ~ P a ) }2(-iyS^Yi0i0(p'a,Pa) (3.25) Pa 2 (47r)5(2\/1 + l)(2i2 + l) \/ h (2\/3+ 1) \\ 0 ( 2 7 r ) \u00bb * < \u00bb > ( & - & + |s) = \u00a3 ^ \u00bb=l -3 \u2022 iix-i a+i\u00bb y,i.-\u2122. (^ > g a ) y - \u2122 . ( S ) ' I^(^^(^ .^(l^dr (3.26) since {pa,qa | -0) = \u00a3 i c ^i c(p, \u2022 1 f + (3.29) (3.28) (3.27) (3.30) 3. T h e impulse approx ima t ion 40 * * 9 - f f <* { - ^ [ ( ( 1 , 1 ) 2 \u00ae l ) i ] +^[((1,1)1 + ^ [ ( ( l , 1)0 \u00ae l)x]}v xs-o + *^ | [ ( (1 ,1 )1 \u00ae l)o]v \u2022 5 (3.31) q.va-u o 1 ) 2 \u00ae l ) i ] - ^ | [ ( ( 1 , 1 ) 0 0 -\u00a3 ^ [ ( ( 1 . 3 ) 2 \u00ae l ) i ] } a { - ^ [ ( ( 1 . 1 ) 2 \u00ae l ) i ] +^[((1,1)1 \u00ae l ) i ] -^ [ ( ( 1 . 3 ) 2 \u00ae l ) 1 ] } a . i m - s + { - ^ [ ( ( 1 , 1 ) 2 \u00ae l ) i ] - ^ [ ( ( 1 , 1 ) 1 \u00ae l ) i ] -^ y f [ ( ( i > 3 ) 2 \u00ae i ) 1 ] } \u00ab ? - H ^ + ^ | [ ( ( 1 , 3)2 \u00ae 1)!](7 \u2022 six \u2022 sv \u2022 \u00a7 + i ^ | [ ( ( l , 1)1 \u00ae l)o]\u00ab \u2022 v x 5 (3.32) Note that while a on the left hand side acts on the spin of the spectator nucleon, a on the right acts on the entire nucleus; u and v are mutually commuting vectors that are not concerned wi th the internal momentum (i.e. not p nor q) and commute with o. In this notation, [l]\u00b0 = [((0,0)0 \u00ae 0) 0]; [of- 1 = [ ( ( 0 , 0 ) 0 \u00ae IK], [of- 1 = - [ ( ( 0 , 2 ) 2 \u00ae IH], [ \/ ( 2 L 1 + l ) (2 \/ 2 + l ) (2 \/ 1 + l ) 0 0 J ^ ' ^ ' ^ ^ V ^ r j ^ * (^ (^ )ll{(n\u00b01,L1(75'>75) \u00ae y,f q)) \u00ae ( 1 0 \u00ae T d ) } e | | \u00bb c ( J ) ) \u00ab ( X ) | | ( l 0 \u00ae f ) | | i c ( J ) ) ( 3 . 3 3 ) Notice a specifies the dimension of the matr ix element. When a is not shown explicitly on a reduced matrix element, a = a; that is [((a,b)c d)e] = [((a, b)c[a] ~~~~ d)e]. F(h, l2; a, b, c; i'c, ic) is defined as: F(li,l2;a,b,c;i'c,ic) c a b - l \/ Yff*(v, x Y Y \u00b0 r ( v , x j l f t^t ) , X)dvdx J_C_1 \\h+h+b j h h b 1 ATTK ' \\ a c I J Z 0 h a \\ \/ (2\/ 2 + l ) (2a+Tj 0 0 J V 26 + 1 6hi< ( 3 . 3 4 ) being some dummy angular variables. The notations used here are the same as Br ink and Satchler [28] (also see ref. [29]) except that the Clebsch-Gordon coefficients are denoted by J M J \\ J 2 M i M 2 as opposed to ( J M | J x J 2 M i M 2 ) . Please note that the definition of reduced matrix ele-ment is: (-1) 2J\" J M J ' J \" M ' M ' \u201e (J||7>||J') = ( J M I T j n M \u201e | J ' M ' ) ( 3 . 3 5 ) <*c(^)IK(>i 1 ,\u00a3 1 (p ' .P) \u00ae g)) <8> (lo \u00ab\u00bb T d ) } e | | i c ( ^ ) > is the spin and angular momentum part of the reduced matr ix element between the helion and tr i ton channels. Its calculation is tedious but standard. (i(J)\\\\{(YLL1(P'>P) \u00ae Y?A?,Q)) \u00ae ( lo \u00ae Td)}e\\\\ic(J)) = 3. T h e impulse approx ima t ion 42 4 7 r ( - l ) \u00a3 ' + 2<+c + L ' + S ' + i + S < + d { ( 2 \u00a3 + 1 ) ( 2 \u00a3 , + 1 ) ( 2 < S + 1 ) ( 2 < s , + 1 } < | l | r \u00ab , | | | > < J L ' 011^,^(^,^)110 Xr> (3.36) where Td = 1 0 for d = 0 and ,p)\\\\0 L) (3.37) (^ cC^ ilio ~~~~ T^ J|\u00abc(X)> is the isospin contribution of the reduced matrix ele-ment, which equals B y doing all the procedures mentioned in this section, one can match the I A amplitude piece by piece wi th its E P M counterparts and thus a direct comparison between each piece would become possible. Upon doing all the appropriate couplings and integrating the internal coordi-nates, one would obtain the right hand side of equation (3.24) which is M{a. One then obtains ^ v ia equation (3.11). Figure (3.2) shows the statistical I A R M C photon spectra from wavefunc-tions of various model potentials. Figure (3.3) shows the photon spectra at different values of gp and figure (3.4) plots the statistical R M C rate vs. gp. Since the IA treatment of O M C is less complicated than R M C , all the nucleon momentum (q) terms in the I A O M C amplitude are expanded to 0(\u00a3j). Note the four-momentum transfer at the hadronic vertex Q2 for O M C is (3.38) 3.7 Results 3.7.1 Ord ina ry M u o n Captu re Q2 = - ( l ^ 1 ) 2 * -0 .964m 2 : M\u2122(k'Q,kQ,Mp) = 9*v\u00ab(L {gi +g2v -i> + g3aL \u2022 a + \u2014(g4o \u2022 q + 3. T h e impulse approx ima t ion 43 > s \u2022 a 0.014 0.012 0.01 0.008 0.006 0.004 0.002 X # + EPM calculation for comparison + * + EPM calculation with only GI terms x + IA calculation (Bonn-A) x + IA calculation (Bonn-B) x IA calculation (CD-Bonn) x * IA calculation (Nijmegen) x + IA calculation (Paris) \u2022\u20227\u00a7|i:i;, x IA calculation (AV14) \"\u2022(Js^ .. x Klieb and Rood values \u2022 20 X C X 40 60 photon momentum [MeV] 80 100 Fig. 3.2: R M C photon spectra from two E P M calculations (one with the full Adler and Dothan amplitude and the other with only gauge invariant terms) and from IA calculations using various model potentials. A l l wavefunctions used have 22 Faddeev components and the permutation is projected on the same set of states. The infrared divergent part is not shown. <75 \u2022 q + giOL \u2022 qo \u2022\u00a3> + g8aL \u2022 aq- \u00a3> + g9a -i>q-i> + iglQD x a \u2022 q) + 7^ (<7n 60 M e V ) = 0.1387 s \" 1 . = 9P 97 = 9M( 9s 99 \u00a710 Su \u00a712 9l3 914 rrr mu ) ~ 9 A mu 2M] .rn mu mu ~M{MJ ~ mf + 9AMI Zmu2 4Mp 3 mu m 9v2~MI + 9MMI = 9A 2 m 2 M 2 -9A -gp -gp-2 m 2 M l m2u M l m2u (3.40) 3. T h e impulse approx ima t ion 46 To perform the internal integration, one uses equation (3.28) to (3.32) and the formulae of appendix B to couple the internal operators with the ap-propriate spherical harmonics from the 8 functions. The final product after internal integrations wi l l then be a amplitude that, similar to the R M C one, depends on the external E P M degrees of freedom. The only traces left by internal integration are the \"formfactors\" that are functions of s. In the case of O M C , s = V and therefore s \u2014 Use equation (3.12) to integrate externally with the phase space and get the I A O M C rate. The statistical rates for both the ordinary and radiative muon capture are shown on table (3.1): Potentials O M C rate (statistical) ( s - 1 ) Bonn-A 1368.6 1368.1 1367.8 1357.9 0.6114 0.1387 Bonn-B 1341.3 1340.8 1340.4 1330.7 0.6033 0.1366 C D - B o n n 1336.1 1335.7 1335.3 1326.1 0.6022 0.1364 Nijmegen 1298.0 1297.5 1297.1 1288.2 0.5894 0.1334 Paris 1270.9 1270.4 1270.1 1260.3 0.5800 0.1313 A V 1 4 1271.0 1270.6 1270.2 1260.0 0.5796 0.1313 E P M 1503 0.8189 0.2113 T a b . 3 .1: Various O M C and R M C statistical rates when gP or FP is at the P C A C value. The numbers in the second column are the values obtained by using the \" | \" prescription up to 0{^-) terms. The numbers i n the th i rd column are obtained v ia the correct approach up to 0{-jtj-) terms. The fourth column has values of the O M C rates using \" | \" prescription up to 0((-\u00a3-)2) terms. Numbers in the fifth column are the values obtained by the correct approach up to 0((jjr-)2) terms; the sixth and rightmost columns contain T ^ ( \u00ab > 5 M e V ) and T^(K > 60 M e V ) respectively. 3.7.2 There are several points to note: Discussion 3. T h e impulse approx ima t ion 47 1. Using the correct approach to handle the momentum terms instead of the traditional \" | \" approach 7 decreases the O M C statistical rate by a mere 0.5 s _ 1 in al l I A calculations up to first order in nucleon momen-tum terms. In the prescription that is used here, there are 24 reduced matrix elements instead of 3 in the \" | \" prescription. The decrease in capture rate is due to the non-zero value of those extra \"internal formfactors\" which can be attributed to the non-vanishing amplitude of the non-S wave components of the tri-nucleon wavefunction since the \" | \" prescription is exact for S waves. The smallness of the effect is pr imari ly due to the minute contribution of the P-state wavefunction to the trinucleon wavefunction. Please note that the effect is a genuine effect and not a effect caused by numerical calculation. To prove this, a two channel Yamaguchi wavefunction consisting solely of S-waves is used to gauge the numerical uncertainty in wavefunction integration. There is an increase of the rate (due only to numerical integration) by 0.1 s _ 1 for the correct approach. This difference is much smaller than the difference (stemming from errors in both the numerical integration and the | approach) of the calculations of other regular 22-channel wavefunctions. Even though the \" | \" approach is exact for S waves up to first order, it is no longer valid at second order and thus a correct prescription for second order nucleon momentum terms is essential, as can be seen by a 0.6% decrease in capture rate when the correct pre-scription is used for second order nucleon momentum terms. 2. The determinable Adler and Dothan terms have a significant influence on the E P M calculation of R M C rates but their effect on the I A cal-culation is much smaller. They make the ratio of the I A R M C total capture rate to the E P M R M C total capture rate about 10% smaller than the O M C counterparts 8 . For O M C , the ratio of the I A calcula-tion to the E P M calculation ranges from 85 to 90% but in R M C the same ratio is just around 75%. If only the gauge terms 9 were added 7 The \" j \" prescription, first used by Peterson [30] in his O M C calculation, replaces all nucleon momentum terms q by | . 8 Notice that the influence of those extra Adler and Dothan terms other than those GI terms on the IA amplitude is minimal. One can barely distinguish the two IA photon spectra (one with the full Adler and Dothan amplitude and the other with only GI terms) visually when they are plotted on figure (3.2) 9 The gauge terms are the terms in A J f a d r o n i c R M C (equation (1.35)) that are explicitly 3. T h e impulse app rox ima t ion 4 8 but not the full Adler and Dothan terms, the ratio would be about 83% (see figure (3.2)). One might blame the discrepancy on the failure of the IA in describing R M C by 3 H e but that seems unlikely given its good predictions for O M C . Indeed, the Adler and Dothan procedure seems suspicious when it is applied to both the E P M and I A descrip-tion of the 3 H e \u2014>\u2022 3 H transition, in particular the E P M one when the derivatives of formfactors are large 1 0 . Since the magnitude of the Adler and Dothan terms are large in E P M due primarily to the large size of the derivatives of the nuclear form factors, it is important that al l the extra terms (to a sufficient high order) should be determined in order to have the full effect of those terms on the resulting E P M capture rate. However, this does not happen: the terms of O(KQ) are nor determinable (or at least cannot be determined uniquely) by the Adler and Dothan procedure. Some of those undetermined terms of O(KQ) carry coefficients of magnitude of ( \/ X 3 H E \u2014 n*ii)RM and these terms are actually first order terms disguised as second order and might have a considerable influence on the resulting spectrum. A s a result, a full meson exchange current ( M E C ) calculation seems necessary to tell whether the inadequacy, if there is any, lies wi th the I A or the E P M or both. Congleton and Truhl ik [31] calculated the M E C contribution to the simpler problem of O M C by 3 H e and found that I A + M E C predic-t ion of the exclusive O M C statistical rate agrees wi th both the E P M calculation and experiment [21]. 3. The E P M calculation of both the O M C and R M C rate and the I A calculation of O M C rate using wavefunctions of non-Bonn potentials agree wi th that of K l i eb and Rood . However, the I A R M C spectrum is about 4% lower 1 1 than theirs. A t first it looks a bit contradictory that a more or less the same I A O M C rate but a slightly lower R M C spectrum is obtained in this work but a closer look reveals that K l i eb and R o o d used a lot of approximations evaluating the reduced matrix elements for the R M C spectrum which they did not use for O M C . independent of K (this independence does not include the pion pole terms). They are the terms determined by gauge invariance alone. 1 0 In their original paper, Adler and Dothan applied only the procedure to the case of R M C by a free proton. 1 1 This already takes into account the fact that they took the value of C = 0.965 while C = 0.9788 is used in this work. 3. T h e impulse approx ima t ion 49 (a) Instead of evaluating the reduced matrix elements directly in terms of s, they expressed s in terms of an infinite sum of spherical har-monics of v and k and an artificial cutoff was imposed on them. (b) They did not fully square the resulting matrix element: only prod-ucts of any two most dominant terms and products of one domi-nant and one small term were considered. (c) In expanding the plane wave exp(is \u2022 r), they only included the term having jo(vr)jo(Kr) and they used this approximation as the premise in deriving several relationships between various reduced matrix elements for R M C . None of these approximations is used here. Some terms in A J ^ \" a d l . o n i c R M C that are present here but not included by Kl i eb and R o o d also tend to decrease the resulting photon spectrum. 4. The Bonn type potentials seem to give a higher (and perhaps better) R M C and O M C results than the other potentials. To analyze this properly, let us take a look at the three dominant reduced matr ix ele-ments: [ l ]\u00b0 , [a] 0 ' 1 and [a] 2 ' 1 . A l l the curves produced by the non-Bonn potentials seem to be a bit separated from the curves of the Bonn type potentials, especially i n the region when s is large. For [1]\u00b0 the prob-lem may be partly associated with the numerical normalization since [1]\u00b0(0) should be one in principle. However, even though one takes this into account (say, for example, scale the Paris potential curve so that it agrees wi th al l others at s = 0) the value of [l]\u00b0 is st i l l smaller than that of the Bonn type potentials, as can be easily seen from the fact that the fractional deviation is larger at large value of s. This seems to have to do wi th the binding energy differences by the Bonn and non-Bonn potentials: the non-Bonn potentials generally underbind the trinucleon by about 1 M e V . Congleton and Truhl ik [31] explained since [1]\u00b0 ~ l \u2014 const(r2)s2 (To see this, expand j0(^sr) in polynomials, note const > 0) and ( r 2 ) scales like the inverse of binding energy so one can expect a lower value for [1]\u00b0 when a wavefunction with a lower binding energy is used. For [ E P M translation 5 6 B. Second order terms for IA <-> E P M translation The second order terms of the I A amplitude after coupling wi th the S func-tions are shown below. Note rf1^ is defined as V2a + 1 I 0 11) (B1) I A after coupling and reexpressing in E P M format qxa-vq-w <-> { - ^ ^ [ ( ( 2 , 2)2 l ) x ] - ^ ^ [ ( f r 2)1 \u00ae l ) x ] -^r 7 ( V ) [((2, 2)0 \u00ae 1):] + ^ ^ [ ( ( O , 2)2[2] \u00ae l)x] + 1 5 ' 7 ( 1 2 1 ) [ ( ( 2 , 0 ) 2 \u00ae 1 ) 1 ] -9 ^ W W , 0)0[2] \u00ae l ) ^ x w \u2022 a + ^ \/ l 7 7 ^ ^ 2 ' 2 ) 2 \u00ae ~ 7 6 ^ 1 ) [ ( ( 2 ' 2 ) 1 \u00ae + yf^V^CCO, 2)2(2] \u00ae l)i]}w; x 5 \u2022 \u00ab\u2122 \u2022 5 + {\\\/ |(^ ( 12 1 )[((2 ) 2)1 \u00ae l)x] + 7?(V )[((2, 2)0 \u00ae 1):] -^ [ ( ( O , 2)2[2] \u00ae l)x])}v xs-aw-s-l-V(\\\\(2,2)l\u00ael)0}w-v + ffi\\(2,2)l\u00ael)0]v-sw-s (B.2) b. Second order terms for I A <-> E P M translation 57 q-vq-wa-u {-^=^[((2, 4)2 \u00ae 1)^ + ^1^1,[((2, 2)2 \u00ae 1):] --J=r7(121)[((2, 2)1 \u00ae \u2022 wu \u2022 sv \u2022 s + { - ^ \u00bb 7 ( V ) [ ( ( 2 , 4 ) 2 \u00ae l ) 1 ] - ^ T ? ( , i 1 ) [ ( ( 2 l 2 ) 2 \u00ae l ) 1 ] -1 r? ( 121 )[((2,0)2\u00ae l)1]}o :-zZJir-T7+. 15 {- v 7=^ ( 1 2 1 ) [ ( (2 ) 4)2\u00ae 1),] - 2 ^ T ? ( V ) [ ( ( 2 , 2 ) 2 \u00ae + y|r?(V)[((2, 2)0 \u00ae l) x]}a \u2022 \u2022 S T J \u2022 i + { ~ V ^ O ^ 1 ^ 2 ' 4 ) 2 \u00ae + h ^ \u2122 ^ 2 ' 2 ) 2 0 \" ^r7(V )[((0,2)2[2] 1),] + ^ r 7 ( V ) [ ( (2 ,0 )2 \u00ae 1):] - . ^v \/3r?(V )[((0,0)0[2]\u00ael)1]-^r 7 (V )[((2,2)0(g)l) 1]}<7-7ju;.7j-{\"7!io^ 1 ) [ ( ( 2' 4 ) 2 0 + { Y I 7 1 ^ ^ 2 ) 2 0 + 1 r7(12a)[((2,0)2<8> \\)\\]}8 \u2022 vu \u2022 w + 15 5 ^ f \/ r t l ^ M l V U . \/ l j U l r { 2 [ ( ( 2 ' 4 ) 2 0 + V n77 2 [ ( (2>2)2 \u00ae + -^=r?(121)[((2, 2)1 \u00ae l)!]}a \u2022 su \u2022 wv \u2022 s + { ; J u f ( ' 1 > l ( ( 2 , 4 ) 2 0 + \/n T ? ( \u00bb 1 ) [ ( ( 2 , 2 ) 2 0 + -^r ? (V )[((2 ! 2)1 \u00ae l)!]}a \u2022 su \u2022 vw \u2022 s+ {-y2=^ 12 1 )[((2,4)2 \u00ae ljx] - ^ r ? ( 1 2 ' 1 ) [ ( (2 ) 2)2 \u00ae l)x] + B. Second order terms for IA <-> E P M translation 58 ^-r\/V ' tttO, 2)2 \u00ae l ) ! ]}a \u2022 sw \u2022vu \u2022 s + { V 7 f l 0 ^ 1 ) [ ( ( 2 ' 4 ) 2 \u00ae + {f\/^2' 2 ) 2 \u00ae ^ \" ^ \u00bb r r [ ( ( 2 , 2 ) 1 \u00ae i)i]}<7 \u2022 \u2022 \u2022 s + r? ( 12 1 )[((2, 2)1 l)o]u x u \u2022 sw \u2022 s + ?7 (121)[((2, 2)1 l)0]w xu- sv-s -fc \u201477 ( 121 )[((2, 4)2 \u00ae 1),] ^r7 ( 1 2 ' 1 )[((2,2)l(8i0) 1]t;x s \u2022 a w \u2022 S + (B.4) \"W^KC2, 2 ) 1 \u00ae \u00b0)i]^ x s \u2022 \u2022 s + v 3 { - ^ ( 1 2 1 ) [ ( ( 2 , 2 ) 0 \u00ae 0 ) 0 ] -2 V \/ 3 r 7 ( 1 o 1 ) [ ( ( 0 ) 0 ) 0 [ 2 ] \u00ae 0 ) 0 ] } t ; - u J + Second order terms for I A <-> E P M translation \/fi C. Wavefunction characteristics 60 C. Wavefunction characteristics C . l Channel specifications A l l the 22-channel wavefunctions [25] have the following channel specifica-tions. These 22 channels consist of a l l possible states of the trinucleon system up to and including J = 2, where J is the total angular momentum of the subsystem particle. Channel L I C S s \/ Channel L I C S S I 1 0 0 0 0 1\/2 1 12 2 0 2 1 3\/2 0 2 1 1 0 0 1\/2 0 13 2 2 0 1 1\/2 0 3 1 1 1 0 1\/2 0 14 2 2 1 1 1\/2 0 4 2 2 0 0 1\/2 1 15 2 2 1 1 3\/2 0 5 2 2 1 0 1\/2 1 16 2 2 2 1 3\/2 0 6 1 1 0 1 1\/2 1 17 1 3 2 1 3\/2 1 7 1 1 1 1 1\/2 1 18 3 1 2 1 3\/2 1 8 1 1 1 1 3\/2 1 19 3 3 0 1 1\/2 1 9 1 1 2 1 3\/2 1 20 3 3 1 1 1\/2 1 10 0 0 0 1 1\/2 0 21 3 3 1 1 3\/2 1 11 0 2 2 1 3\/2 0 22 3 3 2 1 3\/2 1 Tab. C.l: Specifications for 22-channel wavefunctions C.2 Model dependent quantities of wavefunctions Table (C.2) lists several important quantities that are dependent on the po-tential of the wavefunction. They are respectively the binding energy given C . Wavefunction characteristics 61 by the wavefunctions, various partial wave probabilities and the numerical normalization (ip | V)num of the wavefunctions. The experimental binding energy Eb of 3 H e is 7.72 M e V and 3 H is 8.48 M e V [34]. Potential E\u201e P(S) P(S') P(P) P(D) {1> 1 ^ ) n u m Bonn A 8.29 92.59% 1.23% 0.030% 6.14% 0.994 Bonn B 8.10 91.61% 1.19% 0.044% 7.16% 0.993 C D Bonn ( 3 He) 7.91 91.61% 1.35% 0.041% 7.01% 0.993 C D Bonn ( 3 H) 7.93 91.63% 1.31% 0.041% 7.01% 0.993 Nijmegen 7.66 90.31% 1.29% 0.065% 8.34% 0.990 Paris 7.38 90.11% 1.40% 0.069% 8.42% 0.988 A V 1 4 7.58 89.86% 1.15% 0.082% 8.90% 0.987 Tab. C.2: Some important quantities of trinucleon wavefunctions. Binding energy Eb in M e V . P(S) denotes the probability of S-wave and so on. D . Nota t ions and conventions 62 D. Notations and conventions Below are some notations and conventions that are often used in this work: 1. Units: Natural units are used. That is, h = c = 1. 2. Four vectors: Let U and V be two four vectors, their four \"dot product\" is defined as Note that U\u00b0 = U0. Sometimes, U is used to denote \\U\\ but the context w i l l always clear the confusion. 3. 7 matrices: The Dirac representation of 7 matrices (the representation of ref. [35]) is used whenever an explicit representation is required. 4. Spherical harmonics: The bipolar spherical harmonic Y^'^3(x, y) is defined as 5. e: e = +\\\/4wa where a is the fine structure constant whenever e is used to denote electric charge quantities. 6. G F : G F is the Fermi coupling constant. G F = 1.16639 x 1 0 - u M e V - 2 . 7. V u d : V u d is the C K M matr ix element connecting up and down quarks. |Vud| = 0.9735 \u00b1 0.0008 [27]. 8. m: m is used to denote the mass of a muon = 105.6583568 \u00b1 5.2 x I O \" 6 M e V [27] throughout. U -V = UaVa = U\u00b0V0 -U-V ( D . l ) YC(mr(y) (D.2) D. No ta t ions and conventions 63 9. m\u00bb: m* is the mass of a charged pion = 139.57018 \u00b1 0 . 0 0 0 3 5 M e V [27]. 10. M t : M t is the average value of masses of 3 H e and 3 H . M t = 2808.66 M e V . The difference between the two masses is ~ 0.6 M e V . 11. Mp: Mp is the average value of masses of a proton and a neutron. B i b l i o g r a p h y 64 Bibliography [1] Y . L i n , The Importance of L i v i n g (Wi l l i am Morrow, New York) 1998. [2] D . H . Wright, T R I U M F research proposal E592, 1992. [3] D . H . Wright et al , Few B o d y Systems Suppl. 0, 1-7 (Springer-Verlag) 1999. [4] L . Kl i eb and H . P. C . Rood , Nuclear Physics A356 (1981) 483-497. [5] L . K l i eb , P h . D . thesis, University of Groningen, The Netherlands, 1982. [6] N . C . Mukhopadhyay, Physics Reports 30 (1977) 1-144. [7] D . F . Measday, Physics Reports (2001), in press. [8] S. L . Adler and Y . Dothan, Phys. Rev. 151 (1966) 1267. [9] C . W . K i m and H . Primakoff, Phys. Rev. 140 (1965) B566. [10] C . W . K i m and H . Primakoff, Phys. Rev. 139 (1965) B1447. [11] H . W . Fearing, Phys. Rev. C21 (1980) 1951. [12] W . - Y . P. Hwang and H . Primakoff, Phys. Rev. C18 (1978) 414. [13] D . H . Beck et al, Phys. Rev. C30 (1984) 1403. [14] C . R . Ottermann et al , Nuclear Physics A436 (1985) 688. [15] F . -P . Juster et al , Phys. Rev. Lett . Vo l . 55, No. 21 (1985) 2261. [16] P. C . Dunn et al, Phys. Rev. C30 (1984) 1403. [17] J . Delorme, Mesons in Nuclei , edited by M . Rho and D . H . Wilk inson (North-Holland, Amsterdam) 1979, V o l I, 107-147. B i b l i o g r a p h y 6 5 [18] J . G . Congleton, P h . D . thesis, University of Br i t i sh Columbia, 1992, unpublished. [19] J . G . Congleton and H . W . Fearing, Nuclear Physics A552 (1993) 534. [20] G . B . Arfken and H . J . Weber, Mathematical Methods for Physicists (Academic Press) 1995. [21] P. Ackerbauer et al, Physics Letters B417 (1998) 224. [22] H . W . Fearing, Phys. Rev. Lett . 35 (1975) 79. [23] W . Schadow, W . Sandhas, J . Haidenbauer and A . Nogga, Few Body Systems 28, 241-258 (2000). [24] W . Schadow, P h . D . thesis, Bonn University, Germany, BONN-IR-97-17 , 1997. [25] W . Schadow, private communication, 2000. [26] W . Glockle, The Quantum Mechanical Few Body Problem, (Springer-Verlag) 1983. [27] D . E . Groom et al . (Particle Da ta Group). Eur . Phys. Jour. C15 (2000) 1 ( U R L : h t t p : \/ \/ p d g . l b l . g o v ) . [28] D . M . Br ink and G . R. Satchler, Angular Momentum (Oxford Science Publications) 1993. [29] A . R . Edmonds, Angular Momentum in Quantum Mechanics (Princeton University Press) 1996. [30] E . A . Peterson, Phys. Rev. 167 (1968) 971-980. [31] J . G . Congleton and E . Truhlik, Phys. Rev. C29 (1996) 956. [32] A . D . Lahiff and I. R . Afnan, Phys. Rev. C56 (1997) 2387-2395. [33] W . H . Press, S. A . Teukolsky, W . T . Vettering and B . P. Flannery, Numerical Recipes in C (Cambridge University Press) 1992. [34] A . H . Wapstra and G . A u d i , Nuclear Physics A432 (1985) 1. Bib l iog raphy 66 [35] J . D . Bjorken and S. D . Drel l , Relativistic Quantum Mechanics (McGraw-Hi l l ) 1964. [36] F . Halzen and A . D . M a r t i n , Quarks and Leptons: an introductory course in modern particle physics (John Wiley & Sons) 1984. [37] H . P. C . Rood and H . A . Tolhoek, Nuclear Physics 70 (1965) 658-692. [38] L . Kl i eb and H . P. C . Rood , Phys. Rev. C29 (1984) 215. [39] L . Kl i eb and H . P. C . Rood , Phys. Rev. C29 (1984) 223. [40] H . Kameyama, M . K a m i m u r a and Y . Fukushima, Phys. Rev. C40 (1989) 974-987. [41] H . W . Fearing and M . S. Welsh, Phys. Rev C46 (1992) 2077. [42] L . Lamport , I ^ T E X : a document preparation system (Addison-Wesley) 1994. [43] E . Truhl ik and F . C . Khanna , nucl-th\/0102006 (2001). [44] M . E . Peskin and D . V . Schroeder, A n introduction to Quantum Fie ld Theory (Addison-Wesley) 1995. [45] A . Messiah, Quantum Mechanics (Dover) 1999. [46] J . L . Friar, Nuclear Physics 87 (1966) 407. [47] P. Chr is t i l l in , M . Rosa-Clot and S. Servadio, Nuclear Physics A345 (1980) 331. [48] J . Govaerts and J - L . Lucio-Mart inez, Nuclear Physics A678 (2000) 110. [49] W . Savitch, Problem Solving with C++ (Addison-Wesley) 1996. [50] M . M a i n and W . Savitch, Data Structures and Other Objects Using C++ (Addison-Wesley) 1997. [51] K . M . Heal, M . L . Hansen and K . M . Rickard, Map leV Release 5 Learn-ing .Guide (Springer) 1998. [52] M . B . Monagan, K . O. Geddes, K . M . Heal, G . Labahn and S. M . Vorkoetter, Map leV Programming Guide (Springer) 1998. B i b l i o g r a p h y [53] T: Wi l l i ams and C . Kelley, gnuplo t Manual 1998. [54] L . Greenfield, The L i n u x Users' Guide 1996. ","type":"literal","lang":"en"}],"http:\/\/www.europeana.eu\/schemas\/edm\/hasType":[{"value":"Thesis\/Dissertation","type":"literal","lang":"en"}],"http:\/\/vivoweb.org\/ontology\/core#dateIssued":[{"value":"2001-05","type":"literal","lang":"en"}],"http:\/\/www.europeana.eu\/schemas\/edm\/isShownAt":[{"value":"10.14288\/1.0085185","type":"literal","lang":"en"}],"http:\/\/purl.org\/dc\/terms\/language":[{"value":"eng","type":"literal","lang":"en"}],"https:\/\/open.library.ubc.ca\/terms#degreeDiscipline":[{"value":"Physics","type":"literal","lang":"en"}],"http:\/\/www.europeana.eu\/schemas\/edm\/provider":[{"value":"Vancouver : University of British Columbia Library","type":"literal","lang":"en"}],"http:\/\/purl.org\/dc\/terms\/rights":[{"value":"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.","type":"literal","lang":"en"}],"https:\/\/open.library.ubc.ca\/terms#scholarLevel":[{"value":"Graduate","type":"literal","lang":"en"}],"http:\/\/purl.org\/dc\/terms\/title":[{"value":"Radiative muon capture by 3He","type":"literal","lang":"en"}],"http:\/\/purl.org\/dc\/terms\/type":[{"value":"Text","type":"literal","lang":"en"}],"https:\/\/open.library.ubc.ca\/terms#identifierURI":[{"value":"http:\/\/hdl.handle.net\/2429\/11546","type":"literal","lang":"en"}]}}~~