The µ d capture rate in effective field theory

Transcription

1 Physics Letters B The µ d capture rate in effective field theory S. Ando, T.-S. Park, K. Kubodera, F. Myhrer Department of Physics and Astronomy, University of South Carolina, Columbia, SC 98, USA Received 5 September 1; received in revised form 11 March ; accepted March Editor: W. Haxton Abstract Muon capture on the deuteron is studied in a framework that essentially incorporates heavy baryon chiral perturbation theory HBχPT. It is found that by far the dominant contribution to capture comes from a region of the final three-body phasespace in which the energy of the two neutrons is sufficiently small for HBχPT to be applicable. The single unknown low-energy constant having been fixed from the tritium beta decay rate, our calculation contains no free parameter. Our estimate of the capture rate is consistent with the existing data. The relation between capture and the νd reactions, which are important for the SNO experiments, is briefly discussed. Elsevier Science B.V. All rights reserved. PACS: 1.39.Fe; 3.4.-s 1. Introduction Electroweak processes in the two-nucleon systems invite detailed studies for multiple reasons. From the nuclear physics point of view, these processes offer a valuable testing ground of the basic inputs of nuclear physics. In astrophysics, the precise knowledge of the pp fusion cross section is of crucial importance for building a reliable model for stellar evolution [1]. Furthermore, experiments at the Sudbury Neutrino Observatory SNO [,3] to observe solar neutrinos with a heavy-water Cerenkov counter have made it extremely important to estimate the νd reaction cross sections with high precision. In this note we study capture: µ + d ν µ + n + n, in a formalism motivated by effective field theory EFT. Our work is connected to the above-mentioned urgent need of accurate estimates of the νd cross sections, σ νd. To expound this connection, we will first explain the standard nuclear physics approach SNPA, see, e.g., Ref. [4]. This is a highly successful method for describing nuclear responses to electroweak probes. In this approach we consider one-body 1B impulse approximation terms and two-body B exchange-current terms acting on non-relativistic nuclear wave functions, with the exchange currents derived from a one-boson exchange model. The vertices in the relevant Feynman diagrams are obtained from a Lagrangian constructed to satisfy the lowenergy theorems and current algebra [5], while the nuclear wave functions are generated by solving the A-body addresses: sando@nuc3.psc.sc.edu S. Ando, tspark@nuc3.psc.sc.edu T.-S. Park, kubodera@sc.edu K. Kubodera, myhrer@sc.edu F. Myhrer //$ see front matter Elsevier Science B.V. All rights reserved. PII: S

2 6 S. Ando et al. / Physics Letters B Schrödinger equation, H Ψ A =E Ψ A, where the Hamiltonian H contains realistic phenomenological nuclear potentials. The most elaborate study of capture based on SNPA was carried out by Tatara et al. TKK [6] and by Adam et al. [7]. Now, the best available estimation of σ νd based on SNPA is due to Nakamura et al. NSGK [8], while that based on EFT is due to Butler, Chen, and Kong BCK [9]. Since EFT is a general framework [1], it can give model-independent results, provided all the low-energy coefficients LEC in the effective Lagrangian, L eff, are predetermined. L eff considered by BCK, however, does contain one unknown parameter L 1A, which they adjusted to reproduce σ νd obtained by NSGK. After this adjustment, the results of BCK are found to be in perfect agreement with those of NSGK. The fact that an ab initio calculation modulo one free parameter based on EFT reproduces σ νd of NSGK extremely well offers strong support to the calculation based on SNPA. At the same time, it stresses the importance of carrying out an EFT calculation free from an adjustable parameter. An interesting possibility is to use capture data as input to control the unknown LEC. An immediate question, however, is whether this process is gentle enough to be amenable to EFT. The substantial energy transfer accompanying the disappearance of a muon can lead to a region of the final three-particle phase space in which the intrinsic state of the two neutrons receives such a large momentum that the applicability of EFT becomes a delicate issue. Let this unfavorable kinematical region be called the dangerous region. The problem of the dangerous region is reminiscent of the difficulty one encounters in applying EFT to threshold pion production in N + N N + N + π [11 14]. It will turn out see below, however, that, unlike the pion production case, capture receives only a tiny fraction of contribution from the dangerous region, and therefore the theoretical uncertainty caused by the dangerous region is practically negligible. The EFT calculation in [9] used the power divergence subtraction scheme PDS [15]. We employ here a formalism in which the transition operators are derived from irreducible diagrams in heavy-baryon chiral perturbation theory HBχPT, while the transition matrix elements are obtained with the use of the initial and final nuclear wave functions obtained in SNPA. For convenience, we refer to this approach as EFT*. The use of the SNPA wave functions causes some degree of deviation from genuine EFT but, as discussed in [16], EFT* can nevertheless reduce the model-dependence of SNPA drastically. According to [16], a next-to-next-to-next-toleading order N 3 LO calculation in EFT* contains one unknown LEC, denoted by ˆd R.LikeL 1A discussed by BCK, the parameter ˆd R controls the strength of a short-range exchange-current term and, once ˆd R is fixed from data, we can make a definite prediction for σ νd. We, therefore, investigate here the relation between ˆd R and the capture rate, Γ. Our study is essentially of exploratory nature, given the present limited accuracy see below of the experimental value of Γ. Another important point concerning ˆd R is that, as emphasized in [16], the strength of ˆd R can be reliably related to the tritium β-decay rate, Γβ t. Thus, using the experimental value of Γ β t, which is known with high precision, one can determine ˆd R and then proceed to make predictions on various two-nucleon weak-interaction processes, including the capture rate, pp fusion rate, and νd cross sections. We will present here the first estimate of the capture rate obtained in this approach; the pp fusion rate has already been discussed in [16], and the νd cross sections will be reported elsewhere [17].. The capture rate Although capture can in principle occur from the two hyperfine states S = 1/ ands = 3/, the capture is known to take place practically uniquely from the hyperfine doublet state. Therefore, concentrating on this dominant capture, we refer to hyperfine-doublet capture simply as capture and denote the hyperfinedoublet capture rate by Γ. The measured value of Γ is Γ exp = 49 ± 4 s 1 [18] and Γ exp = 47 ± 9 s 1 [19]. We remark that a high-precision measurement of Γ is being contemplated at PSI [].

3 S. Ando et al. / Physics Letters B Denoting by L the orbital angular momentum of the two-neutron relative motion in the final state, we can write Γ = Γ L, 1 L=,1,... where Γ L L= is the rate of capture leading to the L state. Here we shall be primarily concerned with Γ since it is this quantity that contains information about ˆd R. The contributions of Γ L L 1 are significant,1 but their calculation is not expected to involve any major EFT-related issues. In general, due to the centrifugal force, the L 1 contributions cannot be too sensitive to short-range physics, which implies that chiral expansion for them should converge rapidly. Specifically, the L = 1 contributions are dominated by the axial-charge AC and E1 transitions, whose one-body operators which are NLO in chiral counting are well known. The lowest order meson-exchange corrections MEC to the one-body operators come from soft one-pion-exchange OPE, which is N LO in chiral counting. These soft-ope terms, dictated by chiral symmetry, are well known, and they are model-independent. For L states, within the accuracy of our evaluation only one-body contributions [6] have to be included. In our exploratory study, therefore, we concentrate on a detailed evaluation of Γ L=,andforΓ L L 1 we simply use the results obtained by Ref. [6]. Muon capture by the deuteron is effectively described by the current current Hamiltonian of weak interactions H W = G V d 3 xl α xj α x + h.c., where the leptonic and the hadronic charged currents are L α x= ψ ν xγ α 1 γ 5 ψ µ x and J α x= V α A α a=1 x i V α A α a= x, respectively, and G V = GeV [1]; α a is the Lorentz isospin index. In the center-of-mass system of the initial µ d atom from which capture occurs, we can safely assume p µ = p d =. Consequently, the four-momentum transfer to the leptonic system, q α p ν p µ α, reads q, q = E ν m µ, p ν.the capture amplitude is then given by 3 f H W i = G V Ψ µ d l α Ψ nn q, p; s 1 s j α q Ψ d s d, where the Fourier-transformed currents are j α q d 3 x e i q x J α x, 4 5 l α e i x q νp ν,s ν L α x µ p µ,s µ =ū ν p ν,s ν γ α 1 γ 5 u µ p µ,s µ. 6 In Eq. 4, Ψ µ d = 1/ πa 3/ is the 1S wave function of the µ d atom, where a m µ + m d /m µ m d α, with α 1/ the fine structure constant. Ψ d s d in Eq. 4 represents the deuteron wave function with the z-component of its spin s d ; Ψ nn q, p; s 1 s represents the final nn wave function, with total nn momentum q, relative nn momentum p = p 1 p /, and the z-components of the neutron spins, s 1 and s. It is easy to obtain Γ = G V Ψ µ d 3 p d 3 p ν 4J + 1 π 3 π 3 πδ E J S = J s 1 s s µ s d Ψnn j α q Ψ d l α 1,s µ ; 1,s d J,S, 7 1 For example, Γ L=1 1 3 Γ according to Ref. [6]. The small correction due to the finite size of the deuteron is taken into account in the actual calculation.

4 8 S. Ando et al. / Physics Letters B where J = 1/, E is the energy difference between the final and initial states, Ψ d Ψ d s d and Ψ nn Ψ nn q, p; s 1 s. 3. HBχPT Lagrangian and the hadronic currents To derive the transition operators, we adopt Weinberg s counting rule. In HBχPT the leading order LO Lagrangian is given by with L = B[iv D + ig A S ]B 1 A C A BΓ A B + f π Tr i µ i µ + f π 4 Trχ + D µ µ + 1 [ ξ, µ ξ ] i ξ R µ ξ i ξl µξ, µ = 1 ξ {, µ ξ } + i ξ R µ ξ i ξl µξ 8 and χ + = ξ χξ + ξχ ξ,wherer µ = τ a Va µ + Aa ν and L µ = τ a Va µ Aa ν denote the external gauge fields. In the absence of the external scalar and pseudo-scalar fields χ = m π, and we define the pion field as ξ = exp i f τ π π. It is convenient to choose the four-velocity v µ and the spin operator S µ as v µ = 1, and S µ =, σ/. The next-to-leading-order NLO Lagrangian including the 1/m N terms in the one-nucleon sector is given in [] while that in the two-nucleon sector is given in [1]. 3 Combining them, we can write the NLO Lagrangian relevant to our case as v L 1 = B µ v ν g µν D µ D ν + c 1 Tr χ + + m N + c m N 4c g A m N v i + 4c 3 i i B [S µ,s ν] [i µ,i ν ] i 1 + c 6 m N [ S µ,s ν] f + µν 4id 1 BS B BB + id ɛ abc ɛ µνλδ v µ ν,a BS λ τ b B BS δ τ c B +, where ɛ 13 = 1, µ = τ a a µ, f µν + = ξ µ L ν ν L µ i[l µ, L ν ] ξ + ξ µ R ν ν R µ i[r µ, R ν ] ξ. 9 We find it convenient to use the dimensionless low-energy constants ĉ s and ˆd s defined by c 1,,3,4 = 1 m N ĉ 1,,3,4, d 1, = g A m N f π ˆd 1,. 1 The values of these low energy constants, ĉ 1,,3,4, are taken from Ref. []: ĉ 1 =.6 ±.13, ĉ = 1.67 ±.9, ĉ 3 = 3.66 ±.8, ĉ 4 =.11 ±.8 and c 6 = κ V = 3.7. These values were determined at tree level or NLO in the one-nucleon sector, which correspond to N 3 LO in our two-nucleon calculation Our definition of the pion field differs from that used in Ref. [1] by a minus sign.

5 S. Ando et al. / Physics Letters B The one-body currents The one-body currents can be obtained either by an explicit HBχPT calculation or by the Foldy Wouthuysen FW reduction of the well-known relativistic expressions. The former method requires one-loop diagrams which consist of vertices from L for N LO contributions and those which contain one vertex from L 1 for N 3 LO contributions; also needed are the corresponding counter-terms from L and L 3. We adopt here the FW reduction method for convenience. Since the range of t q = m µ m µ E ν for capture is small m µ t m µ, the t-dependences in the standard form factors, F1, V t and G At, give only less than % effects, which can be reliably taken into account by expansion in t, F1 V t = 1 + t 6 r V + Ot, F V t = κ V + Ot, G A t = g A 1 + t 6 r A + Ot. Keeping the terms linear in t is consistent with HBχ PT to the order we calculate. Some caution is requiredfor the G P t term, which contains the pion-pole contributions, G P t m N βt m NG A t m π t, where βt is a slowly varying function. Comparison with the explicit HBχPT calculation up to N LO [3] leads to βt = f π g πnn 1 Q 3 13 g A m N 6 r A m π + O Λ 3 = 1 + [. 1.5 ±.3 ] %. χ The authors of Ref. [6] found that Γ is reduced only by % when β increases by 1%. Thus, limiting ourselves to the β = 1 case entails at most.4% error. 4 The resulting one-body vector 1B current components are V, 1B q= [ τi e i q r i 1 + t 6 r V q 8m κ V i q σ ] i p i q i N 4m κ V N 4m, N V 1B q= i τ i e i q r i [ p i + i 1 + κ V q σ i m N κ V i σ i p i 1 q where τi τi x iτ y i /and p i = p i + p i/. The one-body axial-vector current components can be written for convenience as A α 1B = Âα 1B + qα m π t q βâβ 1B, ω 4m N ], which defines  1B,where Â, 1B q= i ˆ A 1B q= i [ g A τi e i q r σi p i i g A τ i e i q r i m N [ σ i 1 + t 6 r A ω σ ] i q, 8m N + p i σ i p i σ i p ] i 1 q σ i q + i q p i. 4m N 16 4 If the deviation of β from 1 were important, we would have to include the one-body pseudoscalar term, P 1B in Eq. 15, as is necessary for the two-body P B = P in Eq. 18.

6 3 S. Ando et al. / Physics Letters B The above equations correspond to N LO in HBχPT [4]. Apart from the mentioned t-dependence of the form factors, the N LO contribution is found to be negligible,.1%, indicating a rapid convergence. We therefore limit ourselves to N LO for the one-body currents although, to be completely consistent, we should in principle include N 3 LO. 3.. Two-body exchange currents Since the evaluation of the two-body exchange current is the focus of this work we discuss the various parts of this exchange current in detail. We write the two-body vector and axial-vector currents as V µ B k 1, k and A µ B k 1, k,where k i = p i p i is the momentum transferred to the ith nucleon. For the two-body vector-charge current we know V B k 1, k = OQ 4 5 [5], which therefore can be ignored. The spatial components of the vector current have a one-pion-exchange contribution of order Q [6]: V B k 1, k = iτ1 τ g A 4f π [ σ1 σ k m π k σ σ 1 k 1 m π k 1 + σ 1 k 1 m π k 1 σ k m π k k k 1 ] + O Q 4, where f π 93 MeV is the pion decay constant, and we make use of the notation τ 1 τ τ 1 τ x iτ 1 τ y. Analogously to Eq. 15, we write the axial-vector current as [16,7], 17 with A α B = Âα B + qα m π t qβ Â β B + P, Â B k 1, k = i τ1 τ g A g A ˆ A B k 1, k = m N fπ P g A m k 1, k { = π m N fπ 8ĉ 1 τ ĉ +ĉ 3 g A 8 σ 1 k 1 4fπ m π g A k m 1 N fπ {[ i τ 1 τ p 1 + 4ĉ 3 τ k + ĉ τ 1 k1 σ 1 k 1 m π + 1, k 1 τ 1 τ σ 1 k c ] 6 τ 1 τ σ k σ 1 q 4 m π k + [ ˆd 1 τ1 σ 1 + τ σ + ˆd τ 1 τ σ 1 σ ] } + 1, σ k m π k } Only one combination of the LEC, ˆd 1 and ˆd, is relevant for the capture process, ˆd R ˆd 1 + ˆd + 1 3ĉ3 + 3ĉ Exactly the same combination of LEC s appears in triton β-decay, pp-fusion and the solar hep process [16]. Adopting the same strategy as in Ref. [16], we fix ˆd R from Γβ t exp, the experimental value of the tritium β-decay rate. 5 The terms of OQ ν correspond to those of N ν LO.

7 S. Ando et al. / Physics Letters B To facilitate the calculations, we perform a Fourier transformation FT of the above two-body currents. To control short-range physics in performing FT, we introduce a Gaussian cut-off regulator S Λ k k = exp 3 Λ. where Λ is a cut-off parameter. It is to be emphasized that, although our calculation without regularization involves no infinities, we still need a regulator since EFT, by definition, is valid only up to a certain momentum scale. The regulated delta and Yukawa functions read δ 3 Λ r d 3 q π 3 S Λ q e i q r = Λ3 4π d 3 q y Λ m, r π 3 S Λ q e i q r 1 q + m. 3/ exp Λ r We remark that this is exactly the same regularization method as used in Ref. [16]. In performing FT, we need to specify the time components of the momentum transferred to the nucleons. Energy conservation imposes the constraint: k1 + k = q = m µ E ν. In our calculation we will adopt the so-called fixed-kinematics assumption FKA [11], where the energy transfer is assumed to be shared equally between the two nucleons, i.e., k1 = k = m µ E ν /, which naturally brings in the quantity m π m π m µ E ν /4. The uncertainty related to FKA becomes large as q grows. The contribution from the large q region, however, will turn out to be so tiny that the assumptions related to ki cause little uncertainty in our calculation. 4, 4 4. The capture rate for the transition to the 1 S nn state The deuteron and the 1 S wave function may be written as ψ d r; s d = 1 [ u d r + S ] 1ˆr w d r χ 1,sd ξ,, ψ r = 1 u rχ, ξ 1, 1 4πr 8 4πr 5 with dr [ u d r + w d r] = 1 and lim u r = sin δ r p [cos pr + cot δ sin pr]. Here S 1 ˆr = 3 σ 1 ˆr σ ˆr σ 1 σ, χ ξ is the Pauli spinor isospinor, and δ is the nn 1 S phase shift. To facilitate numerical work, we approximate E as E = E ν + m + p + E ν 6 4 M pν p + O 4m 3, where m m n = MeV is the neutron mass, M m µ + m d = MeV. In our calculation we will neglect the O p ν p 4m term since, as we shall show, the major contributions comes from the low p p region. 3 Choosing the z-axis along p ν, we write q = p ν = E ν ẑ. This simplifies the structure of the transition amplitudes as ψ j q Ψ d s d =δ sd,m t, ψ ê λ j q Ψ d s d =δ sd,λm λ, 7

8 3 S. Ando et al. / Physics Letters B where ê ± = ˆx ± iŷ/, ê =ẑ,andλ =±1,. We decompose the matrix elements into vector and axial vector current contributions, M t,λ = M t,λ [V ] M t,λ [A], and arrive at Γ L= = G V Ψ µ π p max dp p E ν 1 E ν M 1 + M M t. M 3 Note that M 1 = M +1 [V ]+M +1 [A], M = M [A] and, to the order under consideration, M t = M t [A]. The matrix elements of the vector current are M λ [V ]=λ dr {qu u d j j w d µv ωu 1 m u d + w } d µ V 1 j 1 N + λ4 g A 8fπ q dr u 1/ 1/ dx 4m N 8 [ j x u d j xw d y L 3 yl 1 xqrj1 x u d + w d y L + 1 j x 3 u d j x + j x ] w d y1 L, 9 where jn x j nqrx are the spherical Bessel functions, yn L y nλ m π + 1 4x 4 q,r, y 1Λ m, r r r y Λm, r, y Λ m, r 1 m r 1 r r r y Λm, r. Using Eq. 18, we obtain for the axial current { } Mt [A] = M λ [A] M t [Â] = ga { [Â] Mt [Â] M λ ga f π M λ [Â] = ga 4 dr } + 1 m π t { ω δ λ, q } ωmt [Â] q M [Â] + M [ P ], { 1 u 1 m ud w d j1 qω N 8m u ud j + } w d j N [ 1 ĉ +ĉ 3 g A mµ E ν 8 t6 p dr {[1 + ra g A m N f π 6m u N [ y1 1 dr 3m N [ m N δ λ, q ] 8m N dr u ud w d j1 y 1Λ r, ]u u d j w d j λ u d w d j λ ] } w d j O kin 1 + c 6 r 4 m π 3 y Λ ĉ 3 + ĉ δ λ, q u u d + w d j 1 u u d j w d j λ 3 31

9 where M [ P ] = 4 S. Ando et al. / Physics Letters B m π 3 y Λ ĉ 3 ĉ 4 1 u wd j u d w ] d j λ + ˆd R δ Λ ru u d, 4 g A m N f π dr O kin q = δ λ, 8 u ud w d j j j λ u w d, 1 4 [ ĉ 1 m y 1Λ π r u ud ] w d j1, j + j λ [ u u d w d u ud ] w d j λ 1 3δ λ,j, j n = j n 1 qr, u 1 r = u r u r, u r = u r r 3 u r + 3 u r r r. In the above expressions the curly brackets denote 1B contributions, and for clarity we have suppressed the dependence on r in some equations. 5. Results Table 1 shows Γ L= as a function of the cut-off parameter, Λ. As discussed, the short-range exchange current contribution depends on the single low-energy constant ˆd R, see Eqs. 3,, and ˆd R determined from Γ t β exp is a function of Λ see Ref. [16]. We observe that the variation of Γ L= over the range of Λ under consideration is less than.7s 1.The ˆd R -dependence in the table indicates the importance of the contribution of the short-distance exchange current. Without the ˆd R term, Γ L= would change as much as 16 s 1 for Λ = 5 8 MeV. Thus, renormalizing the ˆd R -term using Γ t β exp reduces the variation of Γ with respect to Λ by a factor, leading to the practically Λ-independent behavior of Γ. Considering this stability we will hereafter only discuss the case corresponding to Λ = 6 MeV and ˆd R = The capture rate contains several interference terms, which are listed in Table in a cumulative manner. We note that the axial charge AC plays only a minor role; its destructive interference with GT decreases the capture rate by 1s 1. Meanwhile, the M1 contribution interferes constructively with GT, increasing Γ by 59 s 1. Furthermore, the two-body MEC in the L = channel increases the capture rate by 13 s 1. Table 1 L = capture rate in s 1 calculated as a function of the cutoff Λ. Also listed are the corresponding values of ˆd R determined from Γβ t exp [16] Λ MeV ˆd R Γ L= [s 1 ] 5 1. ± ˆd R ˆd R = 45. ± ± ˆd R +.13 ˆd R = 45.3 ± ± ˆd R +.7 ˆd R = 45.7 ±.6

10 34 S. Ando et al. / Physics Letters B Table Cumulative contributions to Γ calculated for Λ = 6 MeV and ˆd R = The row labeled 1B corresponds to the case that contains one-body contributions only, while the row labeled 1B + B to the case that includes both the one-body and MEC contributions. The three columns labeled L = show contributions from the L = channel, with the contributions of the different transition operators displayed in a cumulative manner. The fifth column gives contribution from the L 1 channels, as evaluated in TKK, Ref. [6], and the last column shows the sum of the L = andl 1 contributions Γ [s 1 ] L = L 1 Total GT GT + AC GT + AC + M1 1B B + B Table 3 Matrix elements calculated for representative values of E nn MeV and the cumulative L = capture rate for the case: Λ = 6 MeV and ˆd R = In each entry for the matrix element, the first number preceding a + or sign gives the one-body contribution, while the second number gives the two-body contribution E nn M +1 [A] M +1 [V ] M [A] M t [A] Γ L= [s 1 ] Enn max Our final result for Γ L= = 45 s 1 in Table should be compared with TKK s result, Γ L= TKK = 59 s 1. 6 By adding the 1 L 5 contribution, Γ L 1 = 141 s 1, calculated by TKK, we arrive at the total capture rate Γ = 386 s 1, to be compared with TKK s result Γ TKK = s 1. As mentioned earlier, a primary question is whether capture process is gentle enough for applying HBχ PT with reasonable confidence. As noted the dangerous region for HBχPT occurs when the two neutrons carry most of the final energy. To address this issue, it is useful to consider the differential capture rate, dγ /de nn,where E nn m n + p m n is the energy of the final two-neutron relative motion. An equally informative quantity is the cumulative capture rate Γ E nn E nn dγ de nn de nn. From these quantities we can assess to what extent capture is free from the dangerous kinematic region. We show in Table 3 the matrix elements, M +1 [A], M +1 [V ], M [A] and M t [A], calculated for representative values of E nn,andforλ = 6 MeV and ˆd R = Table 3 also gives Γ L= E nn. The graphical representation of Γ E nn can be found in Fig. 1. We learn from Table 3 that the matrix elements decrease quite fast as E nn increases, a feature that can be easily understood as follows. The 1 S nn radial wave function is proportional to sin δ /p =±[p cot δ + p ] 1/. Since the nn scattering length is very large, p cot δ diminishes rapidly when the nn relative momentum p gets small. The examination of Table 3 also reveals that the one-body amplitudes decrease more quickly than the two-body amplitudes. This is a consequence of the softness of the deuteron wave We have re-run the code of TKK using g A = TKK s original result corresponding to g A = 1.6 was Γ L= TKK = 57 s 1.

11 S. Ando et al. / Physics Letters B Fig. 1. Cumulative µ d capture rate in s 1 calculated for Λ = 6 MeV and ˆd R = The dashed line gives the L = contribution, Γ L= E nn, while the solid line shows the total contribution, Γ E nn Γ L= E nn + Γ L 1 E nn. The empty and solid circles for the values at E nn = Enn max 1 MeV, for L = andl, respectively. function, which cannot supply high momentum transfers needed for producing large values of p. As a result, the contributions from high E nn where the applicability of EFT is questionable is negligible. For instance, the contribution to Γ L= from E nn > 3 MeV is just 1.1s 1, and that from E nn > 5MeVislessthan.1s 1. We now can make a rough estimate of the theoretical error associated with this calculation. Uncertainty related to the G P term, Eq. 1 or β is 1s 1, while uncertainty reflecting the Λ-dependence is less than 1 s 1 ; uncertainty in Γβ t exp or that in ˆd R for a given Λ can affect Γ at the level of 1 s 1. Furthermore, owing to the above-discussed gentleness of the capture kinematics, the higher-order corrections to the B MEC should converge rapidly in powers of m µ /m N.1; the uncertainly due to the higher-order contributions in expected to be 1s 1. If we assign a rather conservative error, s 1,totheL 1 contributions obtained in Ref. [6], the overall uncertainty in our estimate becomes 5 s 1 or 1% in the total capture rate. As mentioned, there is a serious disagreement between the two measured values of Γ. Our theoretical result is consistent with Γ exp in Ref. [18]. In the present exploratory study we have not considered radiative corrections [8], which are expected to be smaller than the existing uncertainty in Γ exp. When the planned precision measurement of the Γ at PSI [] is realized, the issue of radiative corrections should certainly be addressed. The EFT approach as described here will provide a useful tool for this purpose as well. Once the accuracy in Γ exp is significantly improved, we will be able to use capture to determine the low energy constant ˆd R, a quantity critically important for the accurate evaluation of the νd cross sections used in the analysis of the SNO experiments. At present the tritium β-decay is a much more accurate source of information on ˆd R than capture, but it is hoped that in the near future Γ will provide an independent constraint on ˆd R. We consider this redundancy extremely important. Acknowledgements We thank P. Kammel and J.-W. Chen for useful communications. This work is supported in part by NSF Grant No. PHY and No. INT

12 36 S. Ando et al. / Physics Letters B References [1] J.N. Bahcall, R. May, Astrophys. J [] The SNO Collaboration, Phys. Lett. B , and references therein. [3] Q.R. Ahmad et al., Phys. Rev. Lett [4] J. Carlson, R. Schiavilla, Rev. Mod. Phys [5] M. Chemtob, M. Rho, Nucl. Phys. A ; E. Ivanov, E. Truhlík, Nucl. Phys. A ; E. Ivanov, E. Truhlík, Nucl. Phys. A [6] N. Tatara, Y. Kohyama, K. Kubodera, Phys. Rev. C [7] J. Adam, E. Truhlík, S. Ciechanowicz, K.-M. Schmitt, Nucl. Phys. A [8] S. Nakamura, T. Sato, V. Gudkov, K. Kubodera, Phys. Rev. C [9] M. Butler, J.-W. Chen, X. Kong, Phys. Rev. C [1] For reviews, see: U. van Kolck, Prog. Part. Nucl. Phys ; S.R. Beane et al., in: M. Shifman Ed., At the Frontier of Particle Physics Handbook of QCD, Vol. 1, World Scientific, Singapore, 1, p. 133; G.E. Brown, M. Rho, hep-ph/131, Phys. Rep., accepted. [11] B.-Y. Park, F. Myhrer, J.R. Morones, T. Meissner, K. Kubodera, Phys. Rev. C [1] T.D. Cohen, J.L. Friar, G.A. Miller, U. van Kolck, Phys. Rev. C [13] T. Sato, T.-S.H. Lee, F. Myhrer, K. Kubodera, Phys. Rev. C [14] S. Ando, T.-S. Park, D.-P. Min, Phys. Lett. B [15] D.B. Kaplan, M.J. Savage, M.B. Wise, Phys. Lett. B ; D.B. Kaplan, M.J. Savage, M.B. Wise, Nucl. Phys. B [16] T.-S. Park, L.E. Marcucci, R. Schiavilla, M. Viviani, A. Kievsky, S. Rosati, K. Kubodera, D.-P. Min, M. Rho, nucl-th/166; T.-S. Park, L.E. Marcucci, R. Schiavilla, M. Viviani, A. Kievsky, S. Rosati, K. Kubodera, D.-P. Min, M. Rho, nucl-th/171. [17] S. Ando et al., in progress. [18] M. Cargelli et al., in: M. Morita, H. Ejiri, H. Ohtsubo, T. Sato Eds., Proceedings of the XXIII Yamada Conference on Nuclear Weak Processes and Nuclear Structure, Osaka, 1989, World Scientific, Singapore, 1989, p [19] G. Bardin et al., Nucl. Phys. A [] P. Kammel et al., nucl-ex/11, to appear in Hyperfine Interactions; P. Kammel, private communication. [1] J.C. Hardy, I.S. Towner, V.T. Koslowsky, E. Hagberg, H. Schmeing, Nucl. Phys. A [] V. Bernard, N. Kaiser, U.-G. Meißner, Nucl. Phys. B [3] V. Bernard, H.W. Fearing, T.R. Hemmert, U.-G. Meißner, Nucl. Phys. A [4] S. Ando, F. Myhrer, K. Kubodera, Phys. Rev. C [5] T.-S. Park, K. Kubodera, D.-P. Min, M. Rho, Phys. Lett. B [6] T.-S. Park, D.-P. Min, M. Rho, Phys. Rev. Lett ; T.-S. Park, D.-P. Min, M. Rho, Nucl. Phys. A [7] T.-S. Park, D.-P. Min, M. Rho, Phys. Rep ; T.-S. Park, I.S. Towner, K. Kubodera, Nucl. Phys. A ; T.-S. Park, H. Jung, D.-P. Min, Phys. Lett. B ; T.-S. Park, K. Kubodera, D.-P. Min, M. Rho, Astrophys. J [8] J.F. Beacom, S.J. Parke, Phys. Rev. D