Dark matter implications of the neutron anomaly

Motivated by the neutron decay anomaly, we reconsider the neutron decay model n→χχχ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n \rightarrow \chi \chi \chi $$\end{document}, where the new species χ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi $$\end{document} plays the role of dark matter. We precisely compute the χ¯n→χχ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{\chi } n \rightarrow \chi \chi $$\end{document} rate finding that fitting the anomaly compatibly with all bounds needs at least two ”generations” of the χ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi $$\end{document} particle.


Introduction
The current measuraments of the neutron life-time present an inconsistency.Two methods are employed and give different values.One called beam method measures the β − decay rate, Γ β n = 1/(888 ± 2 s), by counting the protons produced from a beam of cold neutrons [1,2].The other called bottle method measures the total neutron decay width, Γ tot n = 1/(878.3± 0.3 s), by storing ultra-cold neutrons in a magnetic bottle [3,4,5].These measuraments are in disagreement at ∼ 4.6σ with the Standard Model, that predicts Γ tot n = Γ β .It has been proposed that ∆Γ = Γ tot n − Γ β ≈ 1.2 × 10 −5 (1 ± 0.21) 1 sec can be due to an extra decay channel of the neutron.The new decay needs to be nearly invisible with a branching ratio around 1%. Various authors proposed a n → χγ [7] decay into a new neutral fermion χ with mass M slightly below the neutron mass so that E γ ≈ m n − M is small.This new particle χ does not decay back to SM charged particles if its mass 1 arXiv:2212.04998v1[hep-ph] 9 Dec 2022 is M < m p + m e and can thus be considered as a dark matter candidate, a possibility disfavoured by tests [8,9] which don't see the final photon with the predicted energy.Furthermore, the χ's thermalize inside neutron stars (NS) [10,11] softening the equation of state1 too much, reducing the NS mass (M N S ) below 0.7M while M N S ∼ 2M is observed.In this perspective we consider a different model that was proposed in [6].The new decay is n → χχχ.In the more minimal model χ carries baryon number B χ = 1/3.The chemical potential of χ in thermal equilibrium in a neutron star is fixed in terms of the chemical potential of the conserved charges as µ χ = B χ µ n .This ensures a substantial reduction of the impact on neutron stars.In section 2, we discuss the main features of the dark species.In section 3, we show that all the bounds are satisfied if more than one χ species are considered.In section 4, we investigate direct detections of DM.In section 5, we show the equations of state of the neutron star resulting from calculations for different number of generation N of the particle χ compared with predictions from the SM.The conclusions are presented in the last section.

Model
The particle χ must be light enough to allow the decay n → χχχ, so M m n /3.The DM mass M is strongly constrained by the following thresholds.Proton stability gives M > (m p − m e )/3 = 312.59MeV.Hydrogen decay H → χχχν e is kinematically open for M < (m p + m e )/3 ≈ 312.93 MeV.The bound from Beryllium ( 8 Be) nuclear decay implies M < (m n − E Be )/3 ≈ 312.63 MeV where E Be = 1.664MeV [7].The number of generations N for this new particle will be important when making comparisons with the experimental bounds.If one assumes more generations of χ i the NS mass remains close enough to the SM limit [6].In particular, considering the behaviour of NS mass as function of NS radius, M N S undergoes only a small softening, thus remaining compatible with the data.Different SM computations lead to maximal NS masses between 1.8M and 2.6M and minimal radii between 10 km and 14 km, in apparent agreement with the data.We compute the N = 1 and N = 3 cases.The theory employed is the following, with n the neutron field and Ψ the dark particle field: The new interaction term is unusual because we want three final particles and no anti-particles.Following [6] we consider the vectorial-axial and scalar interactions.We study in detail the cases with N = 3 generations which will prove to be in agreement with the data.

The Neutron Decay Anomaly
In this section we study the lifetime of the free neutron in the case of N = 3 generations of the χ particle.
This will be important in the context of direct detection later.In the previous papers, e.g.[6], were given only estimatations of N > 1 cases.First we calculate the decay widths of the free neutron and subsequently we calculate the lifetimes for various processes involving the neutron.The n → χχχ decay with N = 1 arises from 4-fermion effective operators in Eq. (2.1).The decay widths for scalar left-right and vector-axial couplings were calculated in [6]: The Lagrangian2 for the N = 3 generation process is: The free neutron decay widths using the scalar and vectorial3 couplings are: The decay width in the N = 1 case is proportional to 3 , with = 1 − 3M mn .In the case N > 1 the width is ∝ 2 .The bounds on Γ n→χχχγ , Γ H→χχχνe and Γ H→χχχνeγ are satisfied in both cases [6].

Direct Detection of Dark Matter and Cosmology
The neutron anomaly is reproduced for Λ χn v so the interaction rates of χ particles at T ∼ m n are below electroweak rates.The χ species can be produced by freeze-in [6].The χn ↔ χχ process leads to unusual DM direct detection signals.
It is kinematically open today and the non-relativistic χ and n in its initial state produce relativistic χ with energy E = 2m n /3 [6].In this section we compute its rate and compare with bounds.

Signals from χn → χχ process
We consider the scalar (left-right) coupling, Γ = g L P L + g R P R .Later, calculations with the vector-axial coupling will also be considered.First, we set N = 1 generations for χ and give a calculation of cross section.In particular, for this process, the initial states are non-relativistic and the products of scattering are relativistic.The prediction for σv rel is: Fixing ∆Γ from the free neutron decay rate gives the cross section as a function of the χ mass.The cross section is well below current bounds from direct detection experiments for the value of Λ χn motivated by the neutron decay anomaly.However different experiments are more sensitive.Dark matter scattering makes ordinary matter radioactive since a neutron that disappears within a nucleus leaves a hole, triggering nuclear deexcitations and decays 4 .The expression for effective lifetime is Bounds on the neutron effective lifetime are found in various experiments: 4.9 10 26 yr from KAMIOKANDE [14] 2.5 10 29 yr from SNO [15] 5.8 10 29 yr from KAMLAND [16] (4. 3) The DM-induced neutron lifetime satisfies the bounds in Eq. ( 4.3), if the number of χ generations are larger than one.So, in this perspective we fix N = 3.We take into account the lagrangian in Eq. (3.2).The expressions for the cross sections per relative velocity are:  The plot shows that the N = 3 case satisfies all the bounds.The predictions for N = 1 are excluded by SNO and by KamLand.The main reason for this difference is the proportionality to different powers of the = 1 − 3M mn factor.Γ(n → χχχ) in the N = 1 case is ∝ 3 whereas in the N = 3 case it is ∝ 2 .Since Γ(n → χχχ) for N > 1 are ∝ 2 all these cases satisfy bounds.

Neutron Star Physics
Neutron stars are described by the Tolmann-Oppenheimer-Volkoff (TOV) equations.
These equations can be solved for any equation of state that gives the density ρ in terms of the pressure ℘, starting with an arbitrary pressure at the center, r = 0, where M (0) = 0, and by evolving outwards up to the neutron star radius r = R at which ρ(R) = 0.In the same way one can predict the relation between the radius R and the total mass M .We add new particles χ to neutron and to sub-dominant SM particles, in order to have ρ = ρ n + ρ χ and ℘ = ℘ n + ℘ χ .The energy density and the pressure are Thermal equilibrium of n ↔ χχχ relates the chemical potential of χ, µ χ = M 2 + p 2 χ to the chemical potential of neutrons, µ χ = µ n /3.In our calculations we adopt the BSk24 from [18].

Equations of State
We investigate the case in which the number of generations is N = 2 and the decay involves two different dark particles, χ 1 and χ 2 .The decay is n ↔ χ 1 χ 1 χ 2 .The kinematic threshold imposed by Beryllium is taken as the lower limit for the sum of the masses of the dark particles.Thus, we obtain where E Be = 1.664MeV.The plot in Fig. 5.1 shows the allowed mass range.[18].Left.For the SM (black curve); n ↔ χ 1 χ 1 χ 2 for N=2 case (blue curves); n ↔ χχχ for N=2 case in the limit of equal masses of χs (red curves).Right.The corresponding relationship between the radius and mass of neutron stars.
In Fig. 5.3 we show the results of the TOV equations for different models.The red curve in the plot on the left shows how the equations of state calculated in the SM soften too much if n ↔ χγ is in thermal equilibrium (i.e.µ n = µ χ ) with M = m n .The blue curves in the plot on the left show how thermal equilibrium of n ↔ χχχ (i.e.µ χ = µ n /3), considering the limit of equal masses M = m n /3 and for different number of generations N = 1, 2, 3, leads to a milder change in the equation of state.As the number of generations increases, the curve undergoes an increasing softening compared to the result in the SM.The curve scales approximately linearly with N .Consequently, in the plot on the right we show that the relation between the neutron star mass M and radius R in presence of n ↔ χχχ remains close to the SM limit (in particular taking into account the observed neutron stars with mass M ≈ 2M ) even when the number of generations increases, in contrast with what happens if n ↔ χγ is in thermal equilibrium.In particular, while n ↔ χγ reduces the maximum mass of neutron stars in contradiction to the data, n ↔ χχχ leads to a slight reduction comparable to the current SM uncertainties.The radius of the neutron star is reduced equally slightly and is compatible with the data.Various calculations in the SM lead to maximum neutron star masses between 1.8M and 2.6M and minimum radii between 10 and 14 km [19], in apparent agreement with the data.Therefore, if the SM is able to explain the observed neutron stars, so is its n → χχχ extension.More precise future calculations and observations may be able to verify this slight difference.[18].In the SM (black curve); n ↔ χχχ for different number of generations N (different blue curves); n ↔ χγ in N=1 (red curves).Right.The corresponding relationship between the radius and mass of neutron stars.This shows that the observed NS with masses around two solar masses are compatible with n ↔ χχχ, but not with n ↔ χγ.The solutions below the peaks at smaller radii are unstable.

Conclusions
In this work we investigated the decay model n → χ i χ j χ k for different number of generations of the species χ to solve the neutron lifetime anomaly.The results are in agreement with the experimental bounds from neutron physics.In Section 3 we have shown the neutron decay rates for different couplings in the case N = 3.This implies a fully invisible Hydrogen decay rate compatible with the inverse universe age and with the bound from the Borexino experiment.In Section 4 we calculated the cross section and effective lifetimes for the process χn → χχ, showing that the case N = 3 satisfies the bounds from Kamiokande, SNO and KamLand.In Section 5 we studied the equations of state for different number of generations of χ species.
We proposed a model for the N = 2 case considering two different species χ 1 and χ 2 with different masses, and we showed that it produces an enhancement of the EoS respect to the prediction in the case with equal masses.The equation of state ℘(ρ) of neutron stars shown in Fig. 5.3 are, for each N ≤ 3, comparable with the uncertainties in the SM.Results for the N = 2, 3 cases in the limit of equal masses for the χ i are in agreement with both χn → χχ and neutron stars.These analyses shed more light on the neutron anomaly and we conclude that the n → χχχ decay model is predictive, satisfying the experimental bounds on DM and neutron physics for N > 1.

. 4 )
In previous expressions the density of DM is ρ = 0.4 GeV/cm 3 .The N = 1 case calculations gives equal expressions for the vectorial coupling and the scalar coupling.Γ(n → χχχ) is fixed from the anomaly.We show our results in Fig.4.1.

Figure 4 . 1 :
Figure 4.1: Plot of effective τ n as a function of the mass of the dark particle χ in the allowed kinematic range.We show the two different cases: for N=1, excluded (SNO and KamLand bounds) and for N=3 above the bounds

Figure 5 . 3 :
Figure 5.3: Left.Equations of state for neutron stars considering the BSk24 equation of state from[18].In the SM (black curve); n ↔ χχχ for different number of generations N (different blue curves); n ↔ χγ in N=1 (red curves).Right.The corresponding relationship between the radius and mass of neutron stars.This shows that the observed NS with masses around two solar masses are compatible with n ↔ χχχ, but not with n ↔ χγ.The solutions below the peaks at smaller radii are unstable.