Pionic final state interactions and the hypertriton lifetime

We analyze the contribution of pionic final state interactions (FSI) in the weak decay of the hypertriton. Focusing on the 3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^3$$\end{document}He channel, we find a contribution of the pionic FSI of the order of 18%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$18\%$$\end{document}. Assuming a fixed value for the branching ratio R3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_3$$\end{document} for the decay width into 3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^3$$\end{document}He over the decay width into 3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^3$$\end{document}He and pd final states, we find values for the hypertriton lifetime that are consistent with the world average as well as recent measurements by the ALICE Collaboration.


I. INTRODUCTION
The properties of hypernuclei provide an ideal playground to test the strong interaction in the strange sector.For example, hyperons open a window onto the interior of nuclei since they are not constrained by the Pauli principle.Moreover, they are expected to play an important role in neutron stars [1].Hyperon-nucleon interactions based on chiral effective field theory (EFT) have been derived in Refs.[2][3][4].For light hypernuclei, a description in pionless EFT is also possible [5][6][7][8].An overview of theoretical and experimental advances in hypernuclear physics was given by Gal and collaborators [9].A comprehensive collection of hypernuclear data can be found at [10].Here we focus on the few-body physics of hypernuclei.Since there are no bound hyperon-nucleon and hyperon-hyperon systems, the simplest hypernucleus is the hypertriton, a three-body bound state of a proton, a neutron, and a Λ hyperon.It plays the same role as the "deuteron" in standard nuclear physics and provides a gateway to understand nuclear physics in the strangeness sector.
The properties of the hypertriton have been calculated using phenomenological interaction models as well as effective field theories (see, e.g., Refs.[5,8,[11][12][13][14][15][16]).Furthermore, first lattice QCD calculations of light hypernuclei have become available for unphysical pion masses [17].Since the Λ separation energy of the hypertriton, B Λ , is small compared to the binding energy of the deuteron, B d ≈ 2.2 MeV, the hypertriton can be viewed as a Λd bound state at low resolution.The most frequently cited value for this separation energy is B Λ = (0.13 ± 0.05) MeV [18], resulting in a large separation of the Λ from the deuteron of about 10 fm [8].The recent value from ALICE [19], B Λ = (0.102±0.063 (stat.)±0.067(syst.))MeV, is slightly smaller but fully compatible.
Since the branching ratio R 3 = Γ3 He /(Γ3 He + Γ pd ) for the decay width into 3 He over the decay width into 3 He and pd is known to be in the range R 3 = 0.3−0.4[20,23,25,29,32] and the contribution from the break up into three nucleons is small [33], one can directly relate decay channels to the complete decay width of the hypertriton.In contrast to the direct calculation of both channels as done, e.g., in Ref. [16], however, one loses the predictive power for the branching ratio.In addition, new preliminary results of the STAR collaboration suggest a slightly smaller value, R 3 = 0.272 ± 0.030(stat.)± 0.042(syst.),which is still consistent with the identification of the hypertriton ground state as J = 1/2 + [34].
In the recent theoretical and experimental investigations, the influence of final state interactions of the produced pion with other decay products on the hypertriton lifetime has been discussed.The purely kinematical argument that the large pion momentum should lead to a small final state interaction seems valid.Therefore pionic final state interactions were often neglected in earlier calculations.An explicit treatment of Pérez-Obiol et al. [15] by distorting the outgoing pion waves leads to a contribution of the order of 10%.Moreover, a treatment of the pions using optical potentials for the hypertriton was done in Ref. [35].Both calculations come to the result that there is a small but relevant contribution to the hypertriton lifetime.
In our previous work on the hypertriton lifetime [16], which neglected the substructure of the deuteron, we relied on the kinematical suppression argument to neglect the pionic contribution, as a 10% effect would come with other higher-order corrections from the hypertriton structure.Given the importance of the hypertriton lifetime puzzle, a reevaluation of this assumption in the power counting of the effective theory is in order.In the present work, we explicitly calculate the pionic final state interactions within the framework of [16] and discuss the impact on the hypertriton lifetime.
The FSI calculation of Pérez-Obiol [15] et al. is purely done for the trinucleon-pion final state, and we will follow their strategy.We limit our FSI calculation to the 3 He π − channel and use the known branching ratio R 3 as well as the ∆I = 1/2 rule to connect to the 3 H π 0 as well as the p d π − and n d π 0 breakup channels.Pionic final state interactions can be treated in this channel in the most straightforward way since only the trinucleon-pion interaction contributes.For the breakup into a deuteron and a nucleon a direct calculation of the pionic final state interaction is much more challenging since one has to include the multiple scattering series for the contributions of nucleon-pion and deuteron-pion scattering.
On the experimental side, the data for 3 He−π − scattering is sparse.Experimental data as well as calculations for the scattering length are available (see, e.g., Refs.[36][37][38][39][40]).While there are also cross section data at finite energy [41][42][43], an actual phaseshift analysis is not available.This is problematic since the outgoing pion momentum in the decay is fixed at a rather large momentum of k ∼ 114 MeV and hence a parametrization of the scattering amplitude by the scattering length only is not sufficient.In addition, the off-shell T-matrix would be necessary for a complete treatment.The problems and implications of this situation will be discussed below.
The paper is structured as follows: We start with a short recap of the most important features of the decay of the hypertriton into a trinucleon ( 3 He/ 3 H) state and a pion in Sec.II.In Sec.II, we also investigate two different types of FSI approximations and their calculation, before discussing our result and its implications on the hypertriton lifetime in Sec.III.

II. TRINUCLEON CHANNEL WITH FINAL STATE INTERACTIONS
We start be reviewing the key parts of the calculation for the trinucleon channel carried out in Ref. [16].The contributions in the triton channel can be related to the helium channel using the ∆I = 1/2 rule.In the following, we thus focus on the helium channel.Because we only have two outgoing particles, energy and momentum conservation leads to a fixed pion momentum k.The final decay amplitude reads [16]: where G F is the Fermi coupling constant, I q is the result of a loop integral, and A π , B π are the baryonic decay coupling constants.The wave function normalizations of 3 He and the hypertriton are given by Z3 He and Z3 Λ H , respectively.Moreover, the energy of the outgoing pion is ω k = M2 π + k 2 with k a fixed momentum determined by mass differences and the Λ separation energy B Λ .For further details see Ref. [16]. 1he two contributions to the hypertriton decay are depicted in Fig. 1.Our treatment of FSI assumes that possible contributions from the three-body cut, which arises when the pion, nucleon and deuteron are on-shell, are small.A full calculation of the FSI contribution (right panel of Fig. 1) is beyond our reach because of the lack of a phase shift analysis for 3 He-pion scattering.We therefore establish different approximations for the loop momentum in the rescattering contribution and discuss their implications for the hypertriton lifetime.The evaluation of the spin structure of the final state interaction is straightforward.Since the pion does not carry spin no additional spin structure is obtained, i.e., the spin structure generated by the weak vertex is inherited from the problem without pionic FSI.The two amplitudes for direct decay, M / FSI , and decay with pionic FSI, M FSI , shown in Fig. 1 must be added coherently, The first diagram simply gives the direct decay amplitude M / FSI = F (k). Defining E = k 2 /(2µ) as the kinetic energy of the 3 He and pion in the final state with µ = M π M3 He /(M π + M3 He ), the contribution of the second diagram in Fig. 1 can be written as: Here is the elastic π − − 3 He scattering amplitude, which describes the final state interactions. 2 Note that the dependence of the amplitudes F and T on the energy E has been suppressed.
The principal value integral now determines a momentum k.This momentum is the characteristic momentum of the decay amplitude F .Defining where a factor of π/2 has been pulled out of k for convenience, we obtain The standard on-shell approximation for the π − − 3 He intermediate state can be obtained by simply discarding the contribution of the principal value integral, i.e. setting k = 0. (See, e.g., Ref. [44] for an application to η − 3 He scattering.)In this case, we obtain for the total contribution of the two diagrams in Fig. 1: where we have assumed that the decay amplitude F (k) is real.From the form of P / E it is obvious that this approximation can only lead to a reduction of the total amplitude.
Including the contribution of the principal value integral, we obtain which depends on the characteristic momentum k.
In order to estimate k we evaluate the integrand leading to Eq. ( 1) without delta distribution enforcing energy conservation, cf.Eq. ( 17) in Ref. [16].The result of this calculation is depicted in Fig. 2. It shows a strong peak around k ∼ 100 MeV, resulting in a value k = 2k/π ∼ 64 MeV, which is a good estimate for the characteristic momentum k.The dependence on the Λ separation energy B Λ is shown in the inset.As expected the dependence is very weak, since the width goes to zero as √ B Λ [16].For B Λ ≤ 0.5 MeV, it is of the order of 3% or less and hence can be neglected.In principle, similar characteristic momenta can also be defined for the three-and four-body decays, as done by Kamada et al. [33].In the next section, we use Eqs.(7,8) together with the available experimental information on π − − 3 He scattering to estimate the contribution of pionic FSI to the hypertriton lifetime.

III. RESULTS AND DISCUSSION
The scattering data for pion scattering on the trinucleon is very sparse.Cross section data for charged pion scattering on 3 He is available at higher energies (see, for example, Ref. [42]), while data for lower energies can be found in Ref. [43].However, we are not aware of a phase shift analysis for 3 He−π ± scattering.The scattering length for charged pion scattering off 3 He, in contrast, is well explored in X-ray experiments [36][37][38] as well as in theoretical calculations [39,40].The resulting value is small, a ≈ 50 • 10 −3 M −1 π .However, replacing the phase shift with its threshold value is not expected to be a good approximation for the hypertriton decay as the FSI happens far away from threshold.The black lines represent FSI factors calculated by the on-shell method while the red ones take loop effects into consideration.The dashed lines represent a calculation where the π−He T-Matrix is replaced by its threshold value.The cyan line is the value of the phase shift obtained by an analysis of the data in Ref. [43].
Before considering the available 3 He−π ± scattering data, we estimate the effect of pionic FSI based on general considerations.Results for the correction factors P / E and P E are shown in Fig. 2 as a function of the phase shift δ.The on-shell FSI correction factor satisfies P / E (δ) ∼ |cos(δ)| 2 and thus can only lead to a reduction of the decay width.If the phase shift is approximated by the threshold value determined by the scattering length, denoted by P / E (a) and P E (a) in Fig. 2, the correction is very small.Therefore, we focus now on the more realistic FSI factor P E (δ), which includes both loop effects and information about the phase shift at finite energy.Even small phases, such as 5 • allow correction of the order 10% in the hypertriton width.This would make it of the same order as, e.g., next-to-leading order hypernuclear contributions which can be estimated from the Λd effective range and scattering length as r Λd /a Λd ≈ 0.1 at B Λ = 0.13 MeV [8].Since our previous result for the hypertriton lifetime [16] is consistent with the free Λ lifetime, such a correction would lead to lifetimes larger than 75% of the free one, τ Λ .A 10% correction as obtained by Perez-Obiol [15], would lead to a lifetime of around 0.9τ Λ .The maximum enhancement of about 30% is reached for δ ≈ 30 • .
In the next step, we explicitly use the differential cross section data by Fournier et al. at T lab π = (45.1 ± 1.0) MeV [43].This pion kinetic energy corresponds to a center-of-mass momentum k ≈ 117 MeV which is close to the fixed momentum k ≈ 114 MeV of the outgoing pion in the hypertriton decay and the peak momentum k ≈ 100 MeV from the left panel of Fig. 2. A fit of the theoretical expression for the differential cross section to the experimental data leads to a scattering phase shift δ 0 ≈ 11 • .(Further details of this analysis are given in Appendix A.) This value is indicated by the cyan line in Fig. 2. Using this estimate, we obtain an enhancement of 18% using the correction factor P E (δ).
We now focus on the consequences for the hypertriton lifetime τ as a function of the Λ separation energy B Λ .We again consider two scenarios to obtain the full decay width from the width for decay into a trinucleon and a pion: (i) we use the branching ratio R 3 as calculated in Ref. [16] without pionic FSI and (ii) we use the experimental value for R 3 as input.Using scenario (i) we obtain a lifetime τ ≈ 0.84τ Λ which is almost independent of the Λ separation energy B Λ .Note, however, that this result is not valid in the limit B Λ → 0 where τ should approach the free Λ lifetime τ Λ .In scenario (ii), we use the experimental branching ratio R 3 to acquire the full lifetime of the hypertriton.The result of this procedure is shown in Fig. 3.We apply the average of the branching ratios R 3 obtained in Refs.[20,23,25,29,32] as well as the new proposed value by STAR [34] with and without the inclusion of pionic FSI.For comparison, we show the calculation by Pérez-Obiol [15] (yellow region), the world averages for B Λ and τ , and the result of our previous calculation [16].In order to not clutter the plot any further, we only draw the EFT uncertainty band for the original result.The grey band represents the result obtained in Ref. [16].In addition, we show the results using a fixed branching ratio R 3 to calculate the total width with and without pionic final state interactions, similar to the approach in Ref. [15], which is given in yellow for reference.The green areas give the world average for the hypertriton lifetime as well as the Λ separation energy according to Ref. [10].In order to keep the plot readable we show EFT uncertainty bands only for the results of Ref. [16].
Note that the lifetime goes to infinity as B Λ → 0 since Γ3 He → 0. This is clearly an artefact of assuming a constant value of R 3 .It is not present in the calculation of Ref. [16] where R 3 also vanishes in this limit.Using the experimental value of R 3 , the results of Refs.[15,16] are in agreement with each other within uncertainties.Moreover, the inclusion of pionic final state interactions in both scenarios (i) and (ii) is consistent with the world averages of B Λ as well as the lifetime.Smaller binding energies or lifetimes closer to the free Λ lifetime are also compatible.Higher binding energies seem to have more tension with the current average lifetime.
In this paper, we discussed different methods to include pionic final state interactions in our calculation of the hypertriton lifetime in the EFT approach of Ref. [16].Following the strategy of Ref. [15], we calculate the decay into 3 He and a π − explicitly and recover the full decay width by using the known branching ratio R 3 for the decay width into 3 He over the decay width into 3 He and pd final states.The contribution of the neutral channels is obtained from the phenomenological ∆I = 1/2 rule.Using the on-shell approximation for the FSI correction factor leads to small corrections and an increased lifetime.Lifting this restriction and including loop effects yields contributions similar in size as found in Ref. [15] in a pion distorted wave approach.We obtain an ehancement of the width by 18%, which is consistent with the original power counting of the EFT in Ref. [16].Thus the contribution of pionic FSI must be included at next-to-leading order.Applying a fixed branching ratio allows predictions for the total width as a function of the Λ separation energy.The results obtained in this way are consistent with the current experimental averages of the hypertriton lifetime and binding energy.process is elastic (η i = 1).Least square fits for the dominant contributions stabilise at l = 3.Higher angular momenta do not improve the reduced χ 2 of the fit.The inclusion of inelastic parts in the fit does yield further insights.This is due to the high number of parameters, even for a low order fit, compared to the 13 data points available.The result of the fit with l = 3 is shown in Fig. 4 degrees.Adding inelastic parts to the fit (η i ̸ = 1) does not decrease the reduced χ 2 .The S-Wave phase shift is stable against adding higher partial waves and inelastic parts, however, the data will be clearly overfitted in this cases.Explicit Coulomb effects are expected to be small at high energies and thus neglected.

Figure 1 .
Figure 1.Amplitude for the weak decay of the hypertriton into 3 He and a pion: direct decay (left panel) and decay with pionic final state interactions (right panel).The amplitude F is the weak decay amplitude, while amplitude T describes the strong final state interaction

Figure 2 .
Figure 2. Left panel: Square of the decay amplitude before enforcing energy conservation for a Λ separation energy B Λ = 0.13 MeV.The inset shows the weak dependence of the peak position k on the value of B Λ .Right panel: results for different approximations for the FSI factor P (δ).The black lines represent FSI factors calculated by the on-shell method while the red ones take loop effects into consideration.The dashed lines represent a calculation where the π−He T-Matrix is replaced by its threshold value.The cyan line is the value of the phase shift obtained by an analysis of the data in Ref.[43].

ΛFigure 3 .
Figure3.Lifetime of the hypertriton τ relative to the free Λ lifetime τ Λ .The grey band represents the result obtained in Ref.[16].In addition, we show the results using a fixed branching ratio R 3 to calculate the total width with and without pionic final state interactions, similar to the approach in Ref.[15], which is given in yellow for reference.The green areas give the world average for the hypertriton lifetime as well as the Λ separation energy according to Ref.[10].In order to keep the plot readable we show EFT uncertainty bands only for the results of Ref.[16].