N ¯ N interaction from chiral effective ﬁeld theory and its application to neutron-antineutron oscillations

. An ¯ NN potential is introduced which is derived within chiral effective ﬁeld theory and ﬁtted to up-to-date ¯ NN phase shifts and inelasticities, provided by a proper phase-shift analysis of available ¯ pp scattering data. As an application of this interaction neutron-antineutron oscillations in the deuteron are considered. In particular, results for the deuteron lifetime are presented, evaluated in terms of the free-space n − ¯ n oscillation time, utilizing that ¯ NN potential together with an NN interaction likewise derived within chiral e ↵ ective ﬁeld theory.


Introduction
Recently the Jülich-Bonn group has established a high-precisionNN potential within chiral e↵ective field theory (EFT) [1]. Starting point for that work was Ref. [2] where a new generation of NN potentials derived in the framework of chiral EFT was presented. In that publication a novel local regularization scheme was introduced and applied to the pion-exchange contributions of the NN force. Furthermore, an alternative scheme for estimating the theoretical uncertainty was proposed that no longer depends on a variation of the cuto↵s [3,4]. Those concepts were adopted and implemented in the study of theNN interaction [1]. Specifically, anNN potential was derived up to next-to-next-to-next-to-leading order (N 3 LO) in the perturbative expansion, thereby extending a previous work by the Jülich-Bonn group that had considered theNN force up to N 2 LO [5]. In both cases the strength parameters of the contact terms that arise in the EFT framework [4,6] (i.e. the low-energy constants or simply LECs) have been fixed by a fit to the phase shifts and inelasticities provided by a proper phase-shift analysis ofpp scattering data [7].
As an application of this interaction neutron-antineutron oscillations in the deuteron have been considered [8]. Neutron-antineutron (n −n) oscillations involve a change of the baryon number (B) by two units (|∆B| = 2). An experimental observation would allow a glimpse on physics beyond the standard model, see e.g. [9][10][11]. Since in such oscillations B is violated the process satisfies one of the Sakharov conditions [12] that have been formulated in order to explain the observation that there is more matter than anti-matter in the universe [13].
The key quanitity in this subject is the free n −n oscillation time, ⌧ n−n . The presently best experimental limit on it is ⌧ n−n > 0.86 ⇥ 10 8 s ⇡ 2.7 yr (with 90 % C.L.) [14]. Additional information can be deduced from studies of n −n oscillations in a nuclear environment. Corresponding experiments have been performed, e.g., for 56 Fe [15], 16 O [16], and for the deuteron ( 2 H) [17], while others are planned [18]. In such a case the oscillation process is suppressed as compared to the free situation. The pertinent lifetime ⌧ nuc is commonly expressed in terms of the one in free space as [10] where R is an intranuclear suppression factor, also called reduced lifetime, that depends on the specific nucleus. It can be calculated from nuclear theory and then can be used to relate the measured lifetimes of those nuclei with the free n −n oscillation time [10], see, e.g., Refs. [19][20][21][22][23]. For a long time the suppression factors published in 1983 [20] have been used as standard by experimentalists in the interpretation of their measurements [15,17]. For example, in case of the deuteron the corresponding value is R ⇠ (2.40 − 2.56) ⇥ 10 22 s −1 , a prediction based on the phenomenological antinucleon-nucleon (NN) potentials by Dover and Richard [24,25]. Recently, however, those values have been challenged in a work by Oosterhof et al. [26]. In that study an e↵ective field theory for the |∆B| = 2 interaction is constructed and the quantity R is evaluated within the power counting scheme proposed by Kaplan, Savage, and Wise (KSW) [27,28] for the nucleon-nucleon (NN) andNN interactions. The value of R for the deuteron obtained in that approach is (1.1 ± 0.3) ⇥ 10 22 s −1 , about a factor 2 smaller than the one by Dover et al. [20].
In the light of this controversal situation a new calculation of the suppression factor for the deuteron has been performed by us [8]. It was prompted by the aim to utilize the discussed modern chiral interactions for the involved NN [2] andNN [1] systems. Specifically, in case of theNN interaction most of the available precisepp scattering data (mostly from the LEAR facility at CERN [29]) have appeared only after the publication of the potentials used in Ref. [20]. Therefore, an update is long overdue. Of course, the main motivation was the aforementiond discrepancy reported in Ref. [26] and the prospect to find a plausible explanation for that di↵erence.
The paper is structured in the following way: In Sect. 2 the employedNN potential is introduced and some results fornp, relevant for the calculation of n −n oscillations in the deuteron, are presented. In Sect. 3 a basic description of the formalism for evaluating the n −n oscillations in the deuteron is provided. Results for the oscillations, specifically for the suppression factor R, are presented in Sect. 4. The paper closes with a brief summary.

TheNN interaction in chiral EFT
A detailed description of the derivation of the chiralNN potential using the Weinberg power counting can be found in Ref. [1]. For information on the NN potential employed in the present calculation see Ref. [2]. As already indicated above, the chiral potential contains pion exchanges and a series of contact interactions with an increasing number of derivatives. Up to N 3 LO there are contributions from one-, two-and three-pion exchanges. Those are identical to the ones that appear in the NN potential [2]. However, there is a sign change in case of an odd number of exchanged pions due to its negative G parity, i.e.
being the corresponding contributions to the NN force. On the other hand, the contact interaction, VN N cont , cannot be taken over simply from the NN case. Those contact terms represent e↵ectively the short-range part of the interaction and, therefore, the G parity of the indivudal contributions remains unresolved. Thus, the strength parameters associated with the arising contact terms, the LECs, need to be determined in a fit toNN data. We fix them by fitting to the phase shifts and inelasticity parameters of the PWA of Zhou and Timmermans [7]. How this is done is described in detail in Ref. [1]. Note that there are more independent LECs in theNN case than in NN because in the former system there is no restriction from the Pauli principle. In addition, in our approachNN annihilation (into multi-meson channels) is likewise parameterized by contact terms, see Refs. [1,5] for explicit expressions.
Once the potential is established, the reaction amplitude is obtained from the solution of a relativistic Lippmann-Schwinger (LS) equation. It reads in partial-wave projected form Here p 00 and p 0 are the (moduli of the) center-of-massNN momenta in the final and initial states, respectively, and E k = where k is the on-shell momentum. We adopt a relativistic scattering equation so that our amplitudes fulfill the relativistic unitarity condition at any order, as done also in the NN sector [3,6]. On the other hand, relativistic corrections to the potential are calculated order by order. They appear first at N 3 LO in the Weinberg scheme, see Appendix A in Ref. [1]. L, L 0 , etc., specifies the orbital angular momentum, considering that theNN system can be in uncoupled (spin-singlet or triplet) states where L 00 = L 0 = L = J or in coupled partial waves with L 00 , L 0 , L = J − 1, J + 1.     Since the integral in the LS equation (2) is divergent for the chiral potentials [3,6] a regularization needs to be introduced. For that the regularization scheme of Ref. [2] is utilized, where a local regulator is used for the pion-exchange contributions and a nonlocal regulator for the contact terms: The explict form of the regulation functions is f R (r) = h 1 − exp(−r 2 /R 2 ) i 6 and f ⇤ (p, p 0 ) = exp(−(p 02 + p 2 )/⇤ 2 ). The cuto↵ radius R is varied in the range R = 0.7-1.2 fm [2] where ⇤ = 2/R is used for relating the momentum-space cuto↵ parameter with the cuto↵ radius.
A complete overview of our results forNN scattering up to N 3 LO in chiral EFT can be found in Ref. [1] while possibleNN bound states are discussed in Ref. [30]. Here we focus on quantities that are relevant for the study of n −n oscillations in the deuteron. Specifically, we look at the 3 S 1 -3 D 1 phase shifts in thenp (isospin I = 1) channel (Fig. 1) and thenp cross section, see Fig. 2. Note that data for the latter reaction were not included in the PWA [7] because they are less precise than the ones forpp. The bands in Figs. 1 and 2 represent the estimated uncertainty of our results. Here we follow the suggestion of Ref. [2] and use as measure the expected size of the higher-order corrections together with the actual size of the higher-order corrections, see Ref. [1] for details.
In the study of n −n oscillations below we employ also an earlier (N 2 LO)NN potential published by our group. A description of that interaction can be found in Ref. [5]. We consider these twoNN interactions because they are based on rather di↵erent regularization schemes. In the earlier potential [5] a non-local exponential exponential regulator was employed for the whole potential while, as outlined above, in the N 3 LO interaction [1] a local regulator was adopted for the evaluation of the one-and two-pion contributions. Comparing the pertinent results allows us to shed light on the question in how far the choice of the regulator influences the predictions. For exploring the sensitivity of the results to the deuteron wave function we employ also those from two meson-exchange potentials [31,32], besides the ones calculated consistently within chiral EFT.
The evaluation of n −n oscillations in the deuteron is done along the formalism presented in Refs. [19,20]. However, since our study is performed in momentum space we describe the main steps below. The starting point is the eigenvalue (Schrödinger) equation [19] Here, V np and Vn p are the potentials in the np andnp systems and | np i and | np i are the corresponding wave functions. The systems are coupled via V n−n which is given by the n −n transition matrix element δm n−n where the latter is proportional to the inverse of the n −n oscillation time, i.e. V n−n = δm n−n =~/⌧ n−n [10].
To leading order thenp component | np i obeys the equation where E d is the unperturbed energy of the deuteron and | d i is the corresponding deuteron wave function. The decay width of the deuteron, Γ d , is then [19] We solve Eq. (5) in momentum space. Performing a partial wave decomposition and taking into account the coupling of the 3 S 1 and 3 D 1 channels, the above integral equation reads with L, L 0 = 0, 2. Note that E d is the total energy corresponding to the deuteron, i.e. E d − and the width is provided by The deuteron lifetime ⌧ d is given by ⌧ d =~/Γ d . The interesting quantity is the so-called reduced lifetime R [19,20,23] which relates the free n −n oscillation lifetime with that of the deuteron,

Results and discussion
Our results for the reduced lifetime R for the deuteron are summarized in Table 1. They are based on our N 3 LO interaction with cuto↵ R 0 = 0.9 fm [1]  in Ref. [20] where theNN potentials DR 1 and DR 2 by Dover-Richard [24,25] have been utilized. Furthermore we include results from the calculation of Oosterhof et al. performed directly within EFT on the basis of the KSW approach. In this case R can be represented in a compact analytical form which reads up to NLO [26] Obviously, the only parameter here is thenp 3 S 1 scattering length. All other quantities that enter are well established NN observables, cf. Ref. [26] for details. Note that in that paper, the scattering length Re an p was taken from Ref. [1]. Table 1. Reduced lifetime R calculated for the χEFTNN potentials from Refs. [1,5], together with information on the pertinentnp 3 S 1 scattering length. Results for the Dover-Richard potentials DR 1 and DR 2 are taken from Ref. [20]. The corresponding scattering lengths are from Ref. [33]. Predictions based on Eq. (11), i.e. on the KSW approach applied in Ref. [26], are indicated too.
χEFT N 2 LO [5] χEFT N Table 1 reveals that the values for R predicted by the chiralNN interactions are fairly similar to those obtained for the DR potentials in the past. By contrast, the results based on the framework employed by Oosterhof et al. [26] are rather di↵erent. Since the scattering length from the N 3 LO chiralNN interaction [1] is utilized in that work, the large discrepancy observed in Ref. [26] is certainly not due to di↵erences in Im an p but must be primarily a consequence of the di↵erent approaches.
For investigating the sensitivity of our results to the used ingredients we performed various exploratory calculations. Specifically, we employed the NLO and N 2 LO variants of the consideredNN (and NN) interactions. The corresponding predictions for R were found to lie within a range of (2.48 − 2.65) ⇥ 10 22 s −1 . If one takes this variation as measure for the uncertainty due to the nuclear structure, i.e. the NN andNN interactions (wave functions), a value of roughly R = (2.6 ± 0.1) ⇥ 10 22 s −1 can be deduced. Application of the method proposed in Ref. [2] for estimating the uncertainty to the calculation based on theNN interaction from 2017 [1], say, leads to a slightly smaller uncertainty. We have also varied the deuteron wave functions alone. As an extreme case we even took wave functions from phenomenological NN potentials derived in an entirely di↵erent framework, namely in the meson-exchange picture [31,32]. Also here the obtained values for R remained within the range given above.  [20], provided that these potentials describe thepp data at low energies.
Admittedly, we do not have a ready explanation for the di↵erence of our results (and those of Ref. [20]) to the ones of Oosterhof et al. [26]. However, we believe that it is due to the fact that in the latter work the width Γ d is evaluated following the perturbative scheme  [28]. In that scheme there is no proper deuteron wave function. Rather one works with an e↵ectively constructed wave function that is represented in terms of an irreducible two-point function [26,28]. This seems to work well for some electromagnetic form factors of the deuteron, at least at low momentum transfer [28,34]. On the other hand, the quadrupole moment of the deuteron is overestimated by 40 % [28], which hints that the properties of the wave function at large distances (small momenta) are not that well represented in this scheme. Clearly, this should have an impact on the quantity studied in the present work as well. Note that a comparable agreement (mismatch) with regard to the KSW scheme has been also observed in studies of the electric dipole moment (magnetic quadrupole moment) of the deuteron [35][36][37]. In any case, one should not forget that there is convergence problem of the KSW approach for NN partial waves where the tensor force from pion exchange is present [38]. It a↵ects specifically the 3 S 1 -3 D 1 channel where difficulties appear already for momenta around 100 MeV/c, see [38] and also the discussions in Refs. [4,39].

Summary
In this presentation we have introduced anNN potential [1], established within chiral e↵ective field theory, which has been fitted to up-to-dateNN phase shifts and inelasticities provided by a phase-shift analysis of availablepp scattering data [7].
As an application of this interaction neutron-antineutron oscillations in the deuteron have been considered. In particular, results for the deuteron lifetime have been presented, evaluated in terms of the free-space n −n oscillation time, utilizing thatNN potential in combination with an NN interaction likewise derived within chiral e↵ective field theory. The value obtained for the so-called reduced lifetime R which relates the free-space n −n oscillation time ⌧ n−n with the deuteron lifetime is found to be R = (2.6 ± 0.1) ⇥ 10 22 s −1 , where the quoted uncertainty is due to the NN andNN interactions (wave functions).
Our prediction for R agrees with the value obtained by Dover and collaborators almost four decades ago [20] but deviates from recent EFT calculations, based on the perturbative scheme proposed by Kaplan, Savage, and Wise [26], by about a factor of 2. A possible explanation for the di↵erence could be that the KSW scheme does not involve a proper deuteron wave function. Rather this ingredient is represented e↵ectively in terms of an irreducible twopoint function. It is known from past studies that the KSW approach fails to describe quantities that depend more sensitively on the wave function like, for example, the quadrupole moment of the deuteron [28].