A relativistic model for the non-mesonic weak decay of the 12C hypernucleus Λ Francesco Conti, Andrea Meucci, Carlotta Giusti, and Franco Davide Pacati Dipartimento di Fisica Nucleare e Teorica, Universita` degli Studi di Pavia and Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, I-27100 Pavia, Italy (Dated: January 21, 2011) A fully relativistic finite nucleus wave-function approach to the non-mesonic weak decay of the 12C hypernucleus is presented. The model is based on the calculation of the amplitudes of the Λ tree-level Feynman diagrams for the ΛN → NN process and includes one-pion exchange and one- kaon exchange diagrams. The pseudo-scalar and pseudo-vector choices for the vertex structure are compared. Final-state interactions between each one of the outgoing nucleons and the residual nucleusareaccountedforbyacomplexphenomenologicalopticalpotential. InitialΛN andfinalNN short-rangecorrelationsareincludedbymeansofphenomenologicalcorrelationfunctions. Numerical 1 resultsarepresentedanddiscussedforthetotalnon-mesonicdecaywidthΓnm =Γn+Γp,theΓn/Γp 1 ratio, theaΛ intrinsic asymmetry parameter, and thekinetic energy and angular spectra. 0 2 PACSnumbers: 21.80.+aHypernuclei;24.10.JvRelativisticModels n a J I. INTRODUCTION 0 2 The birth of hypernuclear physics dates back to 1952 [1] when the first hypernuclear fragment originated from the collision of a high-energy cosmic proton and a nucleus of the photographic emulsion exposed to cosmic rays was ] observed through its weak decays, revealing the presence of an unstable particle: this was interpreted as due to h t the formation of a nucleus in which a neutron is replaced by the Λ hyperon, i.e., the lightest strange baryon. A - hypernucleus is a bound system of neutrons, protons, and one or more hyperons. Only the lightest hyperon, the Λ, l c is stable with respect to esoenergetic strong and electromagnetic processes in nuclear systems. Therefore, the most u stable hypernuclei are those made up of nucleons and a Λ particle. We denote with AX a hypernucleus with Z n Λ protons, (A−Z−1) neutrons, and a Λ (Λ-hypernucleus). [ Hypernuclei represent a unique laboratory for the study of strong and weak interactions of hyperons and nucleons 2 through the investigation of hypernuclear structure and decay. The Λ particle, which does not have to obey the v Pauli principle, is an ideal low-energy probe of the nuclear environment which allows a deepening of classicalnuclear 0 physicssubjects,suchastheroleofnuclearshellmodelsandthedynamicaloriginofthenuclearspin-orbitinteraction. 3 6 Hypernuclearphysicsalsoestablishesa bridgebetweennuclearandhadronicphysics,sincemanyrelatedissuescanin 3 principleunravelthe roleplayedbyquarksandgluonspartonicdegreesoffreedominsidenuclei. Inthis direction,the . study of hybrid theories combining meson-exchange mechanisms with direct quark interactions have the potentiality 2 to teach us something on the confinement phenomenon, an issue still far from being satisfactorily understood. 1 9 In Λ-hypernuclei the Λ can decay via either a mesonic or a non-mesonic strangeness-changing weak interaction 0 process. In the nuclear medium the mesonic decay, Λ → Nπ, which is the same decay of a free Λ, is strongly : suppressed, but in the lightest hypernuclei, by the effect of the Pauli principle on the produced nucleon, whose v i momentum (∼ 100 MeV/c) is well below the Fermi momentum. In the non-mesonic weak decay (NMWD) the pion X produced in the weak Λ → Nπ transition is virtual and gets absorbed by neighbor nucleons. Then, two or three r nucleons with high momenta (∼ 400 MeV/c) are emitted. We can distinguish between one and two-nucleon induced a decays, according to whether the Λ interacts with a single nucleon, either a proton, Λp→np (decay width Γ ), or a p neutron, Λn → nn (Γ ), or with a pair of correlated nucleons, ΛNN → nNN (Γ ). Mesons heavier than the pion n 2 can also mediate these transitions. The NMWD process is only possible in the nuclear environment and represents the dominant decay channel in hypernuclei beyond the s-shell. The total weak decay rate is given by the sum of the mesonic (Γ ) and non-mesonic (Γ ) contributions: m nm Γ =Γ +Γ , (1) tot m nm with Γ =Γ +Γ , Γ =Γ +Γ . (2) nm 1 2 1 p n The fundamental interest in the NMWD mode is that it provides a unique tool to study the weak strangeness changing(|∆S =1|)baryon-baryoninteractionΛN →nN,inparticularitsparityconservingpart,thatismuchmore difficult to study with the weak NN → NN transition, that is overwhelmed by the parity-conserving strong NN interaction. Since no stable hyperon beams are available, the weak process ΛN →nN can be investigated only with bound strange systems. The study of the inverse process pn→pΛ would however be useful. 2 Although the relevance of the NMWD channelwas recognizedsince the early days ofhypernuclear physics,only in recent yearsthe field has experienced greatadvances due to the conception and realizationof innovative experiments and to the development of elaborated theoretical models [2–8]. FormanyyearsthemainopenprobleminthedecayofhypernucleihasbeentheΓ /Γ puzzle,i.e.,thedisagreement n p between theoreticalpredictions and experimentalresults of the ratio between the neutron- and proton-induceddecay widths: forallthe consideredhypernucleitheexperimentalratio,intherange∼0.5÷2,wasstronglyunderestimated (by about one order of magnitude) by the theoretical results. The Γ /Γ ratio directly depends on the isospin n p structure of the weak process driving the hypernuclear decay. The analysis of the ratio is a complicated task, due to difficulties in the experimentalextractions, which require the detection of the decay products, especially neutrons, and to the presence of additional competing effects, such as final-state interactions (FSI) of the outgoing nucleons and two-nucleon induced decays, which could in part mask and modify the original information. In the first theoretical calculations the one-pion-exchange (OPE) nonrelativistic picture was adopted as a natural starting point in the description of the ΛN →nN process, mainly on the basis of its success in predicting the basic features of the strong NN interaction. The first OPE models were able to reproduce the non-mesonic decay width Γ = Γ +Γ but predicted too small Γ /Γ ratios [9–14]. It thus seemed that the theoretical approaches tend nm n p n p to underestimate Γ and overestimate Γ . A solution of the puzzle then requires devising dynamical effects able to n p increase the n-induced channel and decrease the p-induced one. In the following years the theoretical framework was improved including the exchange of all the pseudo-scalar and vectormesons,intheformofafullone-meson-exchange(OME)model,orproperlysimulatingadditionaleffects,above all initial short-range correlations (SRC) and FSI, by means of direct quark mechanisms and many-body techniques [9, 10, 12, 13, 15–24]. In particular, the inclusion of K exchanges seems essential to improve the agreement between theory and experiments. Only a few of these calculations have been able to predict a sizeable increase of the Γ /Γ n p ratio [13–15, 24], but no fundamental progresshas been achievedconcerning the deep dynamical originof the puzzle. The situation has been considerably clarified during the very last years, thanks to considerable progress in both experimental techniques [25–31] and theoretical treatments [13–15, 24, 32–40]. From the experimental point of view, the new generationof KEK experiments has been able to measure the fundamental observables for the 5He and 12C Λ Λ hypernucleiwithmuchmoreprecisionascomparedwiththe“old”data,alsoprovidingthefirstresultsofsimultaneous one-proton and one-neutron energy spectra, which can be directly compared with model calculations. Very recently, it has also been possible to obtain for the first time coincidence measurements of the nucleon pairs emitted in the non-mesonic decay,with valuable informationonthe correspondingangularand energycorrelations. These new data further refine our experimental knowledge of the hypernuclear decay rates, also allowing a cleaner and more reliable extraction of the Γ /Γ ratio. From the theoretical point of view, crucial steps towards the solution of the puzzle n p have been carriedout, mainly through a non-trivial reanalysis of the pure experimental results by means of a proper consideration of FSI and rescattering mechanisms, inside the nuclear medium, for the outgoing nucleons, as well as of the two-nucleon induced channel. This strict interplay between theory and experiments is at the basis of the present belief that the Γ /Γ puzzle has been solved. In particular, this is due to the study of nucleon coincidence n p observables,recently measuredatKEK[29,30], whose weak-decay-modelindependent analysiscarriedoutin [32,33] yields values of Γ /Γ around 0.3÷0.4 for the 5He and 12C hypernuclei, in satisfactory agreement with the most n p Λ Λ recent theoretical evaluations [13–15, 24]. New, more precise results are expected from forthcoming experiments at DAΦNE [41] and J-PARC [42]. Anotherintriguingissueisrepresentedbytheasymmetryoftheangularemissionofnon-mesonicdecayprotonsfrom polarizedhypernuclei. Thelargemomentumtransferinvolvedinthen(π+,K+)Λreactioncanbeexploitedtoproduce finalhypernuclearstates characterizedbya relevantamountofspin-polarization,preferentiallyalignedalongthe axis normal to the reaction plane [43, 44]. The hypernuclear polarization mainly descends from a non-negligible spin-flip term in the elementary nπ+ → ΛK+ scattering process, which in turn interferes with the spin-non-flip contribution [45]: in free space, and for |p |=1.05 GeV and θ ≃15o, the final hyperon polarization is about 75%. π K Polarization observables represent a natural playground to test the present knowledge of the NMWD reaction mechanism, being strictly related to the spin-parity structure of the elementary ΛN → nN interaction. Indeed, by focusing on the p-induced channel, experiments with polarized hypernuclei revealed the existence of an asymmetry in the angular distributions of the emitted protons with respect to the hypernuclear polarization direction. Such an asymmetry originates from an interference effect between parity-violating and parity-conserving amplitudes for the −→ Λp→npelementaryprocess,andcanthuscomplementtheexperimentalinformationontheΓ andΓ partialdecay n p rates,whichareinsteadmainlydeterminedbythe parity-conservingcontributions. Asforthe Γ /Γ ratio,FSIcould n p play a crucial role in determining the measured value of this observable. The asymmetry puzzle concerns the strong disagreementbetween theoretical predictions and experimental extrac- tions of the so-called intrinsic asymmetry parameter a . The first asymmetry measurements [43, 44] with limited Λ statistics gave large uncertainties and even inconsistent results. The very recent and more accurate data from KEK- E508 [5, 47, 48] favour small values of a , compatible with a vanishing value. Moreover, the observed asymmetry Λ 3 parameters are negative for 12C and positive (and smaller, in absolute value) for 5He. Theoretical models generally Λ Λ predict negative and larger values of a . FSI effects do not improve the agreement with data [49]. The inclusion, Λ within the usualframeworkof nonrelativisticOME models, ofthe exchangeofcorrelatedanduncorrelatedpionpairs [50]greatlyimprovesthesituation. Indeed,itonlyslightlymodifiesthenon-mesonicdecayratesandtheΓ /Γ ratio, n p but the modification in the strength and sign of some relevant decay amplitudes is crucial and yields asymmetry parameters which lie well within the experimental observations. In particular, a small and positive value is now −→ predicted for a in 5He. This important achievement justifies the claim that also the asymmetry puzzle has finally Λ Λ found a solution. Recent experimental and theoretical studies have led to a deeper understanding of some fundamental aspects of the NMWD of Λ-hypernuclei. Froma theoreticalpoint ofview, the standardapproachtowardsthese topics has been strictlynonrelativistic,withbothnuclearmatterandfinite nucleicalculationsconvergingtowardssimilarconclusions: nonrelativisticfullone-meson-exchangeplustwo-pion-exchangemodels,basedonthe polarization-propagatormethod (PPM)[51,52]oronthewave-functionmethod(WFM)[11,13,14],seemabletoreproducealltherelevantobservables for the 5He and 12C light-medium hypernuclei. A crucial contribution to this achievement is however due to a non- Λ Λ trivial theoretical analysis of KEK most recent coincidence data, based on the proper consideration and simulation of nuclear FSI and two-nucleon induced decays [33, 50, 53, 54]. In those nonrelativistic models many theoretical ingredients are included with unavoidable approximations. Initial-state interactions and strong ΛN SRC are treated in a phenomenologicalway,though based on microscopic models calculations. The inclusion of the full pseudo-scalar andvectormesonsspectra,inparticularthestrangeK andK∗mesons,inthecontextofcompleteOMEnonrelativistic calculations, somewhat clashes with the still poor knowledge of the weak Λ-N-meson coupling constants for mesons heavier than the pion. Their evaluation requires model calculations which unavoidably introduce a certain degree of uncertainty in the correspondingconclusionsabout the Γ /Γ ratio andother observables. In order to reproduce few n p observables, i.e., the decay rates Γ , Γ , and Γ , and the asymmetry parameter a for 5He and 12C, these models nm n p Λ Λ Λ needtoinclude manydynamicaleffects,suchasthe exchangeofall the possiblemesons,plus two-pionexchange,plus phenomenological σ mesons, plus the corresponding interferences, and, moreover, strong nuclear medium effects in theformofnon-trivialFSI.Althoughtheinclusionofmanytheoreticalingredientscanbe consideredasanaturaland desirable refinement of the simple OPE models, all the improvements do not seem to provide a significantly deeper insight into the decay dynamics. Despite all the theoretical efforts, the solution of the Γ /Γ and a parameter n p Λ puzzlesseemsdue to effects,suchasdistortion,scatteringandabsorptionofthe primarynucleonsby the surrounding nuclear medium, rather than to the weak-strong interactions driving the elementary ΛN → nN or ΛNN → nNN process. The re-analysis of the recent KEK experimental data [25, 29, 30, 55] and the corresponding extraction of Γ /Γ are indeed completely independent of the weak-decay-mechanisms[33, 38], but depend strongly on the model n p adopted to describe FSI and on somewhat arbitrary assumptions, e.g. on the ratio Γ /Γ between two-nucleon and 2 1 one-nucleon induced non-mesonic decay rates. Thepresentlyavailableexperimentalinformationonhypernucleardecayisstilllimitedandaffectedbyuncertainties of both experimental and theoretical nature. Moreover,single-nucleonspectra seem to point at a possible systematic protonsunderestimation[38]. Thenewgenerationofexperimentsplannedinvariouslaboratoriesworldwideisexpected to produce more precise data on the already studied observables as well as new valuable information in the form of differential energy and angular decay particles spectra. In spite of the recent important achievements, the NMWD of hypernuclei deserves further experimental and the- oretical investigation. From the theoretical point of view, the role of relativity is almost unexplored. But for a few calculations in [46, 56, 57, 64] no fully relativistic model has been exploited to draw definite conclusions about the role of relativity in the description of the weak decay dynamics. In this paper we present a fully relativistic model for the NMWD of 12C [58]. The adopted framework consists of Λ a finite nucleus WFM approach based on Dirac phenomenology. As a first step the model includes only OPE and one-kaon-exchange (OKE) diagrams, and is limited to one-nucleon induced decay. We are aware that the neglected contributions could play an important role in the decays. Our aim is to explain all the at least qualitative features of the hypernuclear NMWD with a conceptually simple model, in terms of a few physical mechanisms and free parameters. We stress that, dealing with a fully relativistic treatment of the weak dynamics based on the calculation of Feynman diagrams within a covariant formalism, it is quite difficult to directly compare such an approach and its results to standard nonrelativistic OME calculations. We will thus rather focus on the internal coherence and on the theoretical motivations of the model. The model is presented in Sec. 2. Numerical results for the total non-mesonic decay width Γ , the Γ /Γ ratio, nm n p the a intrinsic asymmetry parameter, as well as for kinetic energy and angular spectra are presented and discussed Λ in Sec. 3. The sensitivity to the choiceof the main theoreticalingredients is investigated. The theoreticalpredictions of the model are compared with the most recent experimental results. Some conclusions are drawn in Sec. 4. 4 II. MODEL In this Sectionwe presenta fully relativistic finite nucleus wave-functionapproachto study the NMWD of the 12C Λ hypernucleus. Our model is based on a fully relativistic evaluation of the elementary amplitude for the ΛN → NN process, which, at least in the impulse approximation, is the fundamental interaction responsible for the NMWD. Covariant,complex amplitudes are calculated in terms of proper Feynman diagrams. The tree-leveldiagraminvolves a weak and a strong current, connected by the exchange of a single virtual meson. Integrations over the spatial positions of the two vertices as well as over the transferred 3-momentum are performed. As a first approximation, only OPE and OKE diagrams are considered. Possible two-nucleoninduced contributions are neglected, even if they could play an important role in the hypernuclear decay phenomenology. InterestedreaderscanfindfurtherdetailsaboutthepresentmodelinthePhDthesisofRef. [58],whereanextensive analysis of the adopted formalism as well as of the involved theoretical ingredients is provided. In the calculation of the hypernuclear decay rate the Feynman amplitude must then be properly included into a many-body treatment for nuclear structure. The amplitude is therefore only a part of the complete calculation, but it is the basic ingredientof the model and involvesall the relevantinformation onthe dynamical mechanisms driving the decay process. Shortrangecorrelationsarealsoincludedinthemodel,coherentlywithwhatcommonlydoneinmostnonrelativistic calculations, since the relatively high nucleon energies involved in the hypernuclear NMWD can in principle probe quite small baryon-baryondistances, where strong interactions may be active and play an important role. Following a phenomenological approach, we have chosen to include initial SRC effects by means of a multiplicative local and energy-independent function, whose general form [46] provides an excellent parametrization of a realistic ΛN cor- relation function obtained from a G-matrix nonrelativistic calculation [59, 60]. The problem of ensuring a correct implementation of such a nonrelativistic SRC function within a relativistic, covariant formalism has been addressed in Ref. [61] and shown to be tightly connected with the choice of the interaction vertices. For full generality,we also choose to account for possible strong NN short range interactions acting on the two final emitted nucleons, again adopting a simple phenomenological average correlation function [13] which provides a good description of nucleon pairs in 4He [62] as calculated with the Reid soft-core interaction [63]; such final-state SRC could in principle play an important role, and they complement the final-state interactions between each of one of the two emitted nucleons and the residual nulceus, that is accounted for in our model by a relativistic complex optical potential. A. Coupling ambiguities In order to devise a relativistic treatment of the elementary ΛN → NN process, great care must be devoted to the choice of the Dirac-Lorentz structure for the strong and weak parity-conserving vertices. The pseudo-scalar (PS) prescription, that consists in a iγ Dirac structure, and the pseudo-vector (PV) one, that contains a γ γµ∂ 5 5 µ axial-vector structure, are in principle equivalent, at least for positive energy on-shell states, because they descend from equivalent Lagrangians. However, ambiguities arise when one tries to take into account SRC in terms of a multiplicative local and energy-independent function f(r). Such ambiguities are not of dynamical origin and should not be mistaken as relativistic effects: they are simply bound to the phenomenological way of including (initial and final) short range correlations, by matching a nonrelativistic correlation function within a relativistic Feynman diagram approach. The crucial observation, in this regard, is that it is possible to give theoretical reasons [61] to prefer the PV coupling in its modified version where the 4-derivative operates on the propagator (PV′), over the PS coupling and also over the standard PV one, where the 4-derivative acts on the matrix element. On the one hand, the PV′ choice permits to recover,in the nonrelativistic limit, the standardOPEpotential,multiplied by f(r), which is commonly used as the starting point in nonrelativistic calculations,whereas the PS and PV couplings yield a simple Yukawa function in the same limit: this allows, at least in principle, a comparison between relativistic and nonrelativistic results. On the other hand, a microscopic model of (initial) SRC effects, adopting standard PS NNπ vertices and introducing an additional ω-exchange mechanism simultaneous to the OPE dominant one, produces a result analogous to what can be derived in a phenomenological tree-level approach contemplating the inclusion of a SRC function, provided in this case the modified derivative PV′ coupling, rather than the PS one, is employed. The main feature is the development of an explicit dependence of the interaction matrix elements on the exchanged three-momentum (through the momentum involved in the corresponding loop integrals,in the microscopic model, or the derivative effect of the PV′ coupling, in the tree-level phenomenological approach). When dealing with a more complexmodelfornuclearstructure,wedonotgenerallyusepositiveenergyon-shellstates. Stillthe generalmessage keeps its validity, though the details of the explicit calculations may be different. In order to correctly treat SRC, nuclear currents showing a dependence on the 3-momentum transfer q are needed, which in the simple model above correspond to matrix elements between external spinors and intermediate spinors carrying the momentum of the 5 intermediate stateexcitedby the heavymeson. This canbe achievedusing the PV′ coupling actingonthe pionfield, while the use of the PS or standard PV couplings, as done for instance in Ref. [46], would generate nuclear currents independent of q, corresponding, within the considered simple SRC model, to matrix elements between spinors all carrying the external momenta. B. Pseudo-scalar couplings As a first example, we employ a PS coupling for the strong vertex and for the parity-conserving part of the weak interaction. The ΛN →nN fundamental process can then be decomposed into a weak ΛNπ vertex, governed by the weak Hamiltonian H(w) =iG m2 Ψ¯(s)(A +B γ )τ ·φ Ψ(b), (3) ΛNπ F π N π π 5 π Λ and a strong NNπ vertex, driven by the Hamiltonian H(s) =ig Ψ¯(s)γ τ ·φ Ψ(b). (4) NNπ NNπ N 5 π N The Dirac spinorsΨ(b) andΨ(b) are the wavefunctions ofthe bound Λ hyperonandnucleoninside the hypernucleus, Λ N Ψ¯(s) isthe Diracspinorrepresentingthescatteringwavefunctionofeachoneofthe twofinalnucleons,τ isthevector N formed by the three Pauli matrices, and φ is the isovector pion field. The Fermi weak constant G and the pion π F mass m give G m2 ≃2.21×10−7. The empirical constants A =1.05 and B =−7.15 are adjusted to the free Λ π F π π π decayanddeterminethe strengthsofthe parity-violatingandparity-conservingnon-mesonicweakrates,respectively. Finally, g =13.16 is the strong NNπ coupling. The initial Λ and final nucleon fields, Ψ(b) and Ψ¯(s), are defined NNπ Λ N inspace-spin-isospinspaceandtheyaredescribedbya space-spinparttimes atwocomponentisospinor. Inaddition, the Λ field is represented as a pure m =−1/2 state to enforce the empirical ∆I =1/2 selection rule. tΛ The relativistic Feynman amplitude for the two-body matrix element describing the ΛN →NN transition, driven by the exchange of a virtual pion, can be written as T(PS) = iG m2 g d4x d4y fini(|x−y|) fi,π F π NNπ ΛN Z Z × Ψ¯(s) (x)(A +B γ )τaΨ(b) (x) k1,ms1,mt1 π π 5 1 αΛ,µΛ,−1/2 h i × δab∆ (x−y) ffin(|x−y|) π NN × Ψ¯(s) (y)γ τbΨ(b) (y) , (5) k2,ms2,mt2 5 2 αN,µN,mtN h i where Ψ(b) and Ψ(b) are the bound Λ and nucleon wave functions, with α = {nlj} quantum αΛ,µΛ,−1/2 αN,µN,mtN N,Λ N,Λ numbers and total spin (isospin) projections µ (m , m = −1/2), and Ψ(s) , with i = 1,2, are the N,Λ tN tΛ ki,msi,mti scattering wave functions for the two final nucleons emitted in the hypernuclear NMWD, with asymptotic momenta k and spin (isospin) projections m (m ). In both the initial and the final baryon wave functions it is possible to i si ti factor out the isospin 2-spinors as well as the energy-dependent exponentials: Ψ(x) ≡ ψ(x)e−iEx0χmt , where E t=1/2 is the total energy of the considered baryon. The factor fini(|x−y|) represents a short-range two-body correlation ΛN functionactingonthe initialΛ andN baryons,andsimilarlyffin(|x−y|)describespossibleshortrangeinteractions NN betweenthetwofinalnucleonsemergingfromtheinteractionvertex. ∆ (x−y)istheFouriertransformoftheproduct π of the pion propagator with the vertex form factors (supposed to be equal for the strong and weak vertices), i.e., d4q eiq·(x−y) ∆ (x−y)= F2(q2) . (6) π (2π)4 q2−m2 +iε π Z π After performing time integrations in Eq. (5) and taking advantage of the q0 part of the integral in Eq. (6), we get 6 for the relativistic amplitude the expression T(PS) = iG m2 g I d3x d3y fini(|x−y|) fi,π F π NNπ ΛN Z Z × ψ¯(s) (x)(A +B γ )ψ(b) (x) k1,ms1 π π 5 αΛ,µΛ h i × ∆ (|x−y|) ffin(|x−y|) ψ¯(s) (y)γ ψ(b) (y) π NN k2,ms2 5 αN,µN h i × (2π) δ(E +E −E −E ) , (7) 1 2 Λ N where d3q e−iq·(x−y) ∆ (|x−y|)≡ F2 (q0)2−q2 , (8) π (2π)3(q0)2−q2−m2 +iε π Z π (cid:12)q0=qe0 (cid:0) (cid:1) (cid:12) (cid:12) with q0 =E −E =E −E and I is an isospin factor that depends on the conside(cid:12)red decay channel (either p- or Λ 1 2 N n- induced), i.e., e I ≡ χmt1 †τ χ−1/2 · χmt2 †τ χmtN . (9) 1/2 1 1/2 1/2 2 1/2 (cid:20)(cid:16) (cid:17) (cid:21) (cid:20)(cid:16) (cid:17) (cid:21) It is easy to check that I is different from zero only for the charge-conservingprocesses Λp→np and Λn→nn. For the calculation of the integral over the 3-momentum transfer q in Eq. (7) we choose a monopolar form factor, i.e., Λ2 −m2 F (q2)≡ π π , (10) π Λ2 −q2 π where m ≃ 140 MeV is the pion mass and Λ ≃ 1.3 GeV is the cut-off parameter [46]. The initial ΛN correlation π π function adopted in our calculations is [46] fini(r)= 1−e−r2/a2 n+br2e−r2/c2, (11) ΛN (cid:16) (cid:17) with n=2, a=0.5, b=0.25, c=1.28, while the final NN correlation function is chosen as [13] ffin(r)=1−j (q r), (12) NN 0 c where j (x) = sinx is the first spherical Bessel function, and q = 3.93 fm−1. The two correlation functions of Eqs. 0 x c (11) and (12) are plotted in Fig. 1. As a consequance of the approximations adopted in the present calculation, we could then, for practical purposes, treat initial and final SRC as a whole, in terms of an overall correlation functions defined as f(|x−y|)=fini(|x−y|)ffin(|x−y|). (13) ΛN NN C. Pseudo-vector couplings When we use derivative PV′ couplings, the equivalent of Eq. (5) is T(PV ′) = iG m2 g d4x d4y fi,π F π NNπ Z Z iB × Ψ¯(s) (x) A − π γ γµ∂x τaΨ(b) (x) h k1,ms1,mt1 (cid:18) π 2M¯ 5 µ(cid:19) 1 αΛ,µΛ,−1/2 i × δab∆ (x−y)f(|x−y|) π i × Ψ¯(s) (y) − γ γν∂y τbΨ(b) (y) , (14) h k2,ms2,mt2 (cid:18) 2MN 5 ν(cid:19) 2 αN,µN,mtN i 7 fHrL 1.2 1 0.8 0.6 0.4 0.2 0 0 1 2 3 4 5 6 r@fmD FIG. 1: Initial ΛN fini(r) (Eq. (11), solid line) and final NN ffin(r) (Eq. (12), dashed line) correlation functions, with r ΛN NN representing thebaryon-baryonrelative distance. where M¯ = (M +M )/2, with M = 1.1156 GeV, and now f(|x−y|) is defined as in Eq. (13). In Eq. (14) the N Λ Λ space-time derivatives act just on the pion propagator ∆ (x−y) (given in Eq. (8)) and not on the SRC function π f(|x−y|), which is considered as a phenomenological ingredient entering Eq. (14) in a factorized form. Thus, the involved derivatives translate into multiplicative terms under the integral over q. After the time integrations in Eq. (14), we obtain the final expression of the relativistic amplitude T(PV ′) = iG m2 g I d3x d3y fi,π F π NNπ Z Z d3q e−iq·(x−y) × F2(q2) f(|x−y|) (2π)3q2−m2 +iε π "Z π B × ψ¯(s) (x) A + π γ /q ψ(b) (x) h k1,ms1 (cid:18) π 2M¯ 5 (cid:19) αΛ,µΛ i 1 × ψ¯(s) (y) − γ /q ψ(b) (y) h k2,ms2 (cid:18) 2MN 5 (cid:19) αN,µN i × (2π) δ(E +E −E −E ) , (15) 1 2 Λ N # where I, F(q2), and f(|x−y|) are defined in Eqs. (9), (10), and (13), respectively. Eq. (15) must be evaluated at q0 =q0 ≡E −E =E −E . Λ 1 2 N ThecrucialdifferencebetweenEq. (7)andEq. (15)isthatinEq. (15),obtainedadoptingderivativePV′ couplings at the vertices,the matrix elements betweenthe initial bound andthe final scattering states explicitly depend on the e 3-momentum transfer q, while using PS couplings the matrix elements in Eq. (7) are independent of q. Though representinga computational complication, the q-dependence of the matrix elements is a desirable feature in connection with the problem of correctly including short range correlations in a fully relativistic formalism as the onedevelopedhere. TheuseofPS vertices,whichproducesbaryonicmatrixelementsonlydependingontheexternal variables, is not coherent with a simple but significant model of the physical mechanism behind SRC, based on the simultaneousexchangeofapionplus oneormoreheavymesons. The correlatedFeynmanamplitudes involvebox(or more complex) diagrams and one expects that the interaction matrix elements explicitly depend on the momentum 8 involved in the corresponding loop integrals. This in turn represents a strong motivation to consider PV′ couplings as the most appropriate ones for a fully relativistic approach to the NMWD of Λ-hypernuclei, since they prove able to mimic such a physical effect. D. Initial- and final-state wave functions The main theoretical ingredients entering the relativistic amplitudes of Eqs. (7) and (15) are the vertices oper- ators and the initial and final baryon wave functions. Since we adopt a covariant description for the strong and weak interaction operators, the involved wave functions are required to be 4-spinors. Their explicit expressions are obtained within the framework of Dirac phenomenology in presence of scalar and vector relativistic potentials. In the calculations presented in this work the bound nucleon states are taken as self-consistent Dirac-Hartree solutions derived within a relativistic mean field approach, employing a relativistic Lagrangian containing σ, ω, and ρ mesons contributions [68–72]. Slight modifications also permit to adapt such an approach to the determination of the initial Λ wave function and binding energy. The explicit form of the bound-state wave functions reads g (|r|) Yµ(Ω) nκ κ ψµ (r)= , (16) nκ if (|r|) Yµ (Ω)! nκ −κ where the 2-components spin-orbital is written as Yµ(Ω)≡ l µ 1/2 µ j µ Yµl(Ω)χµs , (17) κ l s l 1/2 µl,µXs=±1/2(cid:16) (cid:12) (cid:17) (cid:12) (cid:12) with 1 l =κ if κ>0 j =|κ|− and ; (18) 2 (l =−κ−1 if κ<0 nistheradialquantumnumberandκdeterminesboththetotalandtheorbitalangularmomentumquantumnumbers. The normalization of the radial wave functions is given by d|r||r|2 f (|r|) 2+ g (|r|) 2 =1. (19) nκ nκ Z (cid:16)(cid:12) (cid:12) (cid:12) (cid:12) (cid:17) The outgoingnucleonswavefunctions arecalcu(cid:12)latedby(cid:12)mean(cid:12)softhe (cid:12)relativisticenergy-dependentcomplex optical potentials of Ref. [73], which fits protonelastic-scattering data on severalnuclei in an energy range up to 1040 MeV. In the explicit construction of the ejectile states, the direct Pauli reduction method is followed. It is well known that a Dirac 4-spinor, commonly represented in terms of its two Pauli 2-spinor components ψk,ms(r)=(cid:18)χφkk,,mmss((rr))(cid:19), (20) can be written in terms of its positive energy component φ(r) as φk,ms(r) ψk,ms(r)= M +E+(σS·(|kr)|)−V(|r|) φk,ms(r), (21) (cid:20) N (cid:21)   whereS(|r|)andV(|r|)arethescalarandvectorpotentialsforthefinalnucleonwithenergyE. Theuppercomponent φ(r) can be related to a Schr¨odinger-like wave function φ(r) by the Darwin factor D(|r|), i.e., φ(r)≡ eD(|r|) φ(r), (22) p with e M +E+S(|r|)−V(|r|) D(|r|)≡ N . (23) M +E N 9 The two-component wave function φ(r) is solution of a Schr¨odinger equation containing equivalent central and spin- orbit potentials, which are functions of the energy-dependent relativistic scalar and vector potentials S and V. Its general form is given by e M +E φk,ms(r)= N2E 4πil ulj(|r|) Ylµj;|k|(Ωr) r lmXljµ h i e × lml 1/2ms j µ Ylml∗(Ωk). (24) (cid:16) (cid:12) (cid:17) (cid:12) (cid:12) E. Decay rates Inthecompletecalculationofthetotalandpartialdecayrates,aswellasofpolarizationobservables,thedynamical information on the elementary ΛN →NN process, given by the amplitudes in Eq. (7) or Eq. (15), are included in a many-body calculation for nuclear structure. The weak non-mesonic total decay rate is defined as [46, 65] d3k d3k 1 2 Γ = (2π)δ(M −E −E −E ) nm (2π)3 (2π)3 H R 1 2 Z Z M2 1 × N |M |2. (25) fi E E 2J +1 1 2 H MJXH{R} ms1ms2 mt1mt2 The energy-conservingdelta function connects the sum of the asymptotic energiesE of the two outgoing nucleons, 1,2 coming from the underlying ΛN → NN microscopic process, with the difference between the initial hypernucleus mass M and the total energy E of the residual (A−2)-particle system after the decay. A sum over E is also H R R usually understood. Integration over the phase spaces of the two final nucleons is needed, since the decay rate is a fully inclusive observable. Moreover, the sums in Eq. (25) encode an average over the initial hypernucleus spin projections M , where J is the hypernucleus total spin, a sum over all the spin and isospin quantum numbers of JH H the residual (A−2)-system, {R}≡{J ,M ,T ,M }, as well as a sum over the spin and isospin projections of the R R R TR two outgoing nucleons m and m , respectively. If we choose a reference frame in which, for instance, the zˆ-axis s1,2 t1,2 is aligned along the momentum k , and exploiting the energy-conservationin the delta function, the six-dimensional 1 integral in Eq. 25 can be reduced to a two-dimensional integral, one over the energy of one of the two final nucleons and the other one over the relative angle between the momenta of the two nucleons (due to azimuthal symmetry), which can be performed numerically. The expression for the NMWD rate Γ can be decomposed into a sum over n- and p-induced decay processes nm without any interference effects, i.e., Γ = Γ [m ]=Γ(p) +Γ(n) , (26) nm nm tN nm nm mXtN where Γ [m ] is defined as in Eq. (25) and |M |2 is evaluated with a fixed value of the initial-nucleon nm tN fi isospin projection, m = 1/2 for p-induced and m = −1/2 for n-induced channels. Actually, in each term tN tN of |M |2 the m quantum numbers are fixed, so that |M |2 would involve products of the kind mt1mt2 fi t1,2 fi TA∗ [m ,m ,m′ ]TA [m ,m ,m ],where,inprinciple,alsointerferenceeffects, m′ 6=m ,areallowed. How- fi,Pπ t1 t2 tN fi,π t1 t2 tN tN tN ever, the non diagonal products with m′ 6=m are necessarily zero, since if one of the two amplitudes is non-zero tN tN the other one must vanish as a consequence of the charge-conservation isospin factor I of Eq. (9) (same final state but different initial states, pΛ or nΛ). Therefore, only the diagonal terms, m′ = m , contribute and without tN tN interferences the coherent sum over m becomes an incoherent one. tN Thenucleartransitionamplitude,fromtheinitialhypernuclearstatetothefinalstateofan(A−2)residualnucleus and the two outgoing nucleons, is defined as M =hf|Mˆ |ii (27) fi ΛN→NN and can be represented in terms of the elementary two-body ΛN → NN relativistic Feynman amplitude, T , fi,π which contains all the relevant information about the weak-strong dynamics driving the global decay process. The final A-particle state |fi must be further specified and decomposed into products of antisymmetric two-nucleonand 10 residual (A−2)-nucleon wave functions. An explicit decomposition for the initial hypernuclear wave function |ii can be developed following the approach introduced in Ref. [46], which is based on a weak-coupling scheme, i.e., the isoscalarΛ hyperon is assumed to be in the 1s ground state and it only couples to the ground-state wave function 1/2 of the (A−1)-nucleon core. As discussed in Ref. [46], this weak-coupling approximationhas been able to yield quite good results in hypernuclear shell-model calculations [66]. The final expression for M is fi 2 |M |2[m ] = T M 1/2m T M fi tN R TR tN H TH (cid:16) (cid:12) (cid:17) (cid:12) 2 × A J T J(cid:12) T ,j m c H R R N tN XjN D n(cid:12) E (cid:12) (cid:12) 2 × J M j m J M R R N N c c " XmN (cid:16) (cid:12) (cid:17) (cid:12) × J M 1/2m J(cid:12) M 2 TA 2 , (28) c c Λ H JH fi,π # where m = +1/2 for Γ(p) and −1/2 for Γ((cid:16)n), {J ,M ,T(cid:12)(cid:12)(cid:12) ,M }(cid:17)ar(cid:12)(cid:12)e the(cid:12)(cid:12)spin-isospin quantum numbers for tN nm nm H JH H TH the initial hypernucleus, {j ,m } are the initial-nucleon total spin and its third component, {J ,M } are the same N N c c quantum numbers for the (A−1)-nucleon core, and, finally, m is the initial Λ total spin projection. In Eq. (28) Λ hJ T {|J T ,j i are the real coefficients of fractional parentage (c.f.p.), which allow the decomposition of the c H R R N initial (A−1)-nucleoncore wavefunctions in terms of states involvinga single nucleon coupled to a residual(A−2)- nucleon steatee. The factor A is produced by the combination of initial- and final-state antisymmetrization factors with the number of ΛN pairs contributing to the total decay rate. Eq. (28) neglects possible quantum interference effects between different values of j (and m ), namely we are ruling out interferences between different shells (s N N 1/2 and p ) for the initial nucleon. Thus the calculation does not require the c.f.p., but only the spectroscopic factors 3/2 S =A(c.f.p)2, that can be taken, e.g. from Ref. [46]. F. Antisymmetrization and isospin factors A crucial role in determining the Γ /Γ ratio is played by the isospin content of the model, namely the I factors n p defined in Eq. (9) in terms of the SU(2) isospin operators (generally represented by the 2×2 Pauli matrices) and of the corresponding isospin 2-spinors for the initial Λ and N as well as for the two final nucleons. Takingadvantageofthe∆I =1/2isospinselectionrule,fromtheisospinpointofview,theΛbehaveslikeaneutron state. We can then explicitly represent the isospin spinors for the p, n and Λ baryons as χ ≡χmtN=1/2 = 1 , χ ≡χmtN=−1/2 = 0 , χ =χ . (29) p 1/2 0 n 1/2 1 Λ n (cid:18) (cid:19) (cid:18) (cid:19) With these definitions, the I isospin factors can be evaluated for all the possible combinations of the m , m , and tN t1 m isospin projection quantum numbers. They are non-zero only for those processes in which charge is conserved, t2 namely Λp→np and Λn→nn. We obtain I [m =1/2, m =−1/2, m =1/2] ≡ I(d) = −1, (30) tN t1 t2 Λp→np I [m =1/2, m =1/2, m =−1/2] ≡ I(e) = 2, (31) tN t1 t2 Λp→np I [m =−1/2, m =−1/2, m =−1/2] ≡ I(d) = I(e) = 1. (32) tN t1 t2 Λn→nn Λn→nn Alltheotherspossibilitiesimplychargeviolationandgivezero. The(d)and(e)apicesrefertothedirectandexchange diagrams of the relativistic Feynman amplitudes for the elementary processes. The use of the isospin formalism means that we are treating the neutron and the proton as two indistinguishable particles; therefore the final state is composed of two identical particles and this requires the antisymmetrization of the TA amplitude. The antisymmetrization acts on the two final nucleons, exchanging their spin-isospin quantum fi,π numbers, m , m , and their momenta k within the matrix elements defining the T complex amplitude. We can si ti i fi,π thus define TA ≡T(d) −T(e) , (33) fi,π fi,π fi,π