Axion-assisted resonance oscillation rescues the Dodelson–Widrow mechanism

The keV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textrm{keV}$$\end{document} scale sterile neutrino was a qualified candidate for dark matter particles in the Dodelson-Widrow mechanism. But the mixing angle, needed to provide enough amount of dark matter, is in contradiction with the astrophysical observations. To alleviate such tension, we introduce an effective interaction, i.e. ga(ϕ/Λ)∂μaνα¯γμγ5να\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g_a (\phi /\Lambda )\partial _{\mu }a \overline{\nu _\alpha }\gamma ^{\mu } \gamma _5 \nu _\alpha $$\end{document}, among Standard Model neutrino να\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu _\alpha $$\end{document}, axion a, and singlet ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document}. The axial-vector interaction form is determined by the axion shift symmetry, and the singlet ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document} with dynamically varied vacuum expectation value is introduced to reinforce the axial-vector coupling strength and evade the stringent neutrino oscillation constraints. The effective potential generated by the new interaction could cancel the SM counterpart, resulting in an enhanced converting probability between SM neutrino and sterile neutrino. Hence, the production rate of sterile neutrinos can be substantially enlarged with smaller mixing compared to the DW mechanism.


INTRODUCTION
The mystery surrounding dark matter (DM) production and neutrino mass generation is intriguing.The seesaw mechanism, which introduces an extra neutrino ν s , is the minimal and natural framework for solving these two issues.After the electroweak symmetry breaking, the Dirac mass terms yv να ν s + h.c. will inevitably result in a small mixing between the active and sterile neutrinos, ν 4 = cos θ ν s + sin θ ν α , which makes the sterile neutrino an ideal dark matter candidate in the Dodelson-Widrow (DW) mechanism [1].This mass mixing θ is too tiny in the early universe to thermalize the sterile neutrinos with the Standard Model (SM) bath.Thus there were no initial distributions for sterile neutrino.However, by the collisions of active neutrinos with other SM particles, even a slight amount of mixing can produce a sizable population of sterile neutrinos.The production process will end once the SM neutrinos decoupled from the thermal bath and the relic abundance of sterile neutrinos is attained with θ ∼ 10 −6 and m ν s ∼ 100keV [2].Unfortunately the active-sterile neutrino mixing is already in tension with searches for DM decaying into monochromatic X-rays in this minimal setup [3][4][5][6][7].Together with the bounds from the observation on the dwarf galaxies [2,[8][9][10], the parameter space has been almost ruled out.
Given these challenges, less minimal hypotheses that would increase the rate of sterile neutrinos and reduce the restrictions imposed by astrophysical observations have been proposed.The existence of a significant lepton-number asymmetry in the early universe to enlarge the sterile neutrino production rate resonantly is a common new physics approach [11] to alleviate the above tension.The problem with this approach is that it can not explain why the lepton asymmetry is significantly larger than the baryon asymmetry, and get a strong constraint from Big Bang Nucleosynthesis.Including additional self-interaction between active neutrinos is another el-egant approach [12][13][14] to erasing the tension.It helps efficiently produce sterile neutrino DM and keeps neutrinos in thermal equilibrium with themselves for a longer time.The neutrinophilic scalar or vector is responsible for mediating the self-interaction, which results in a connection between production rate and dark radiation ∆N eff .A similar approach to generating the appropriate relic abundance of ν s is the freezein production via Higgs singlet decay or vector boson decay processes [15][16][17][18][19]. Ref. [20] proposed a solution by introducing interactions with axion-like particles, which breaks the connection of the mixing angles that responsible for sterile neutrino production in the early universe and for the decay of sterile neutrino at late times.The size of the mixing decreases as the energy density of axion-like particles diluted over the cosmological expansion, thus avoiding the Xray bounds.Even though all above proposals can explain the sterile neutrino relic abundance, they go far beyond the minimal configuration.The resonant oscillation seems to be the best option if we insist on preserving the grace of the DW mechanism.
It is demonstrated that the presence of a primordial condensate of the dark photon can alter the dispersion relation of the neutrinos in the early universe [21].A greater active neutrino self-energy allows for a level crossover with its heavier sterile counterpart, enabling resonant oscillations between the two states.We notice that the underlying principle of the approach is that a dark photon originates from the misalignment mechanism where it has to couple to the curvature R with extra substantial, highly tuned couplings.But these couplings trigger perturbative unitarity violations in longitudinal photon-graviton scattering at low energies [22].This straightforward and elegant approach also faces problems.
We point out that the axion can be used to generate resonant oscillation, as it naturally originates from the misalignment mechanism [23,24].The minimal interaction between axion and active neutrino g a ∂ µ aν α γ µ γ 5 ν α suffers from stringent constraints, where the magnitude of g a required for resonance to exist contradicts the neutrino oscillation in the current universe.In this letter, we propose a brand-new method to generate resonant oscillations in the early universe and remove the neutrino oscillation constraints in the late universe through the axion assisted DW mechanism (denoted as aDW hereafter), with interaction g a (φ /Λ)∂ µ aν α γ µ γ 5 ν α .The main idea is that the vacuum expectation value of the singlet φ is not a constant but a function of temperature, where it becomes significant in the early universe and vanishes in the late universe.In contrast to the minimal axion assisted DW mechanism, the ensuing resonant oscillations enable an observable relic abundance associated with significantly reduced portal couplings g a , evading the astrophysical constraints.

PRODUCTION OF STERILE NEUTRINO
Sterile neutrino is among the simplest solutions which generate the nonzero neutrino masses confirmed by oscillation experiments.The mixing between it and SM neutrino will then naturally give birth to the initial abundance of sterile neutrino.An accurate calculation of sterile neutrino production requires solving the evolution equation in the matrix density formalism, while in the collision-dominated period (which it's indeed the case of the production peak of sterile neutrino, i.e.T ∼ 0.1 GeV), the quasi-classical Boltzmann equation is a good approximation [25], here f ν s and f ν α stand for the distribution functions of sterile neutrino and SM neutrino, Γ α is the reaction rate where SM neutrino takes part in.⟨P m ⟩ is the averaged converting probability between sterile neutrino and SM neutrino, expressed as [26,27] ⟨P m (ν α ↔ ν s ; p,t)⟩ ≃ 1 2 with ∆ = m 2 s /2p, the effective potential is shown as V , which characterizes the effects of SM neutrinos passing through the thermal background.The effective potential consists of contributions from SM electroweak interactions and any other new physics.For DW mechanism, there is only the SM electroweak contribution, which is approximated as [21,28] for temperature between the mass of corresponding charged lepton and electroweak phase transition.The α stands for the fine structure constant, s w is the sine of Weinberg angle, and G F is the Fermi constant.
To obtain the relic abundance of sterile neutrinos, one can derive the equation for r ≡ n s /n α , with n i denoting the number density of sterile (SM) neutrino when i = s (α).The resulting equation reads Here entropy density conservation has been considered, and g * s is the degrees of freedom with respect to entropy density.γ on the right hand side is defined as For constant g * s , γ/H stands for the relative number of sterile neutrino, compared to SM neutrino, in each log-interval of T .Thus in the simple DW mechanism there are two free parameters, i.e., the mixing angle θ and the sterile neutrino mass m s .The relic abundance of sterile neutrino will be determined once these two parameters are fixed.In Figure 1, we show the variation of γ/H with temperature during a period of [10 −3 , 10] GeV in the DW mechanism (the red line).We choose the mixing angle as sin 2 2θ = 10 −12 , and the sterile neutrino mass is selected as m s = 10 keV.The production peak reaches γ/H ∼ 10 −6 at temperature ∼ 0.2 GeV.
Such a framework is simple enough to produce sterile neutrino, just from SM neutrino oscillation.However this mechanism is facing crises from various experiments.As the sterile neutrino are much heavier than SM neutrino, decay process of Hence the observation of X-ray [3][4][5][6][7] can put strong bounds on the mixing angle.Moreover, for fermionic DM, there is a lower mass limit due to the requirement that the phase-space density of DM does not exceed that of the degenerate Fermi gas [8].A bound of m s > 2 keV is set from the analysis of DM phase space distribution in dwarf galaxies [2,9,10].As can be seen in Figure 2, the X-ray observation has excluded almost the whole parameter space in DW mechanism.Considering the constraints of dwarf galaxies, there is no room for the DW mechanism to survive.

THE AXION ASSISTED DW MECHANISM
We consider a new interaction of SM neutrino in an effective operator form including the axion field a with mass m a , and a singlet φ .The axial-vector interaction form is fixed by the axion shift symmetry.The inclusion of φ is to avoid severe constraints of neutrino oscillation observations.Scattering between axion and SM neutrino will modify the vacuum dispersion relation of neutrino, hence generating a Mikheyev-Smirnov-Wolfenstein matter effect [30][31][32].In later discussions we can see that, without the inclusion of φ , the coupling strength g a needed to give the correct DM relic density is in strong contradiction with the oscillation observations.The coupling strength matrix g a has a negative mass dimension [M] −1 .Λ is the ultraviolate cutoff scale.Due to the large occupation number nature, the axion fields can be approximated in a classical wave form of a(t, x) = a 0 cos(m a t − ⃗ p •⃗ x), with a 0 = √ 2ρ a /m a and ρ a being the energy density of axions.Equation 6indicates that axions could be generated during the early stages of the universe.They could be produced thermally and maintained equilibrium with the SM bath.However, the resulting density is calculated to be significantly lower than that of cold dark matter [33].Consequently, they cannot be considered the primary component of dark matter.In our consideration, we treat sterile neutrino as a dominant component of dark matter, while for the axion we set ρ a as 1% of the dark matter energy density.
Containing the singlet φ in the effective interaction is an important ingredient in our setup.φ is a morphon field [34] which develops a dynamically varied vacuum expectation value by interacting with Ricci scalar R at the early stage of the universe, The coupling constant ξ quantifies the interaction strength between the morphon field and curvature scalar.One common choice is ξ = 1/6, which comes from the requirement of conformal invariance [35,36].A comprehensive analysis of the temperature dependence of Ricci-scalar is carried out in [37] for the time range from inflationary epoch to the present day.
In our calculation, we consider the evolution of Ricci scalar from T ∼ 1 MeV to T ∼ 10 GeV.Thus, at high temperatures the scalar field develops a nonzero vacuum expectation value, The scalar vacuum expectation value decreases with T until it vanishes at low temperatures.This is the main feature that it helps generate resonant oscillations in the early universe and evade constraints in the late universe, such as the neutrino oscillation bounds.Now the additional contribution to the effective potential of SM neutrino can be obtained from the new introduced effective interaction in Equation 6, in the first-order approximation as The derivative on space part could be neglected as the velocity of axion in the Milky Way is of O(10 −3 )c.Then the effective potential is further approximated as where the plane-wave form of the axion field has been taken.We have mentioned above that the SM weak interactions contribute a negative thermal potential V SM = −C 2 W pT 4 .While the new contribution could be either positive or negative, depending on the axion mass and time (or temperature).Therefore, in principle there exists cancellation between the contributions from SM and the new interaction.It would substantially enlarge the averaged converting probability ⟨P m ⟩, and resonantly produce the sterile neutrino dark matter.As a result, the mixing angle could be much smaller than that in the original DW mechanism.This can help to release the strong tension against the X-ray observations.The cancellation condition is explicitly written as We firstly give a rough estimation to show the incapability of the minimal interaction.Actually the minimal interaction would arise when φ develops a vacuum expectation value, under the condition ξ R(T ) > m 2 φ .Specifically we label the minimal interaction as g ′ a ∂ µ aν α γ µ γ 5 ν α , here g ′ a relates to g a as g ′ a = g a v φ (T )/Λ.The generated potential is in a form of g ′ a √ 2ρ a sin(m a t).The former two terms in Equation 11have different dependencies on temperature, i.e. ∆ cos 2θ ∝ T −1 and V SM ∝ T 5 with setting p = T for simplicity.Thus the lowest value of the sum is estimated as ∼ 4.1(m s /10 keV) 5/3 eV.We find that the coupling strength g ′ a needed for the effective potential, to have a considerable cancellation with the former two terms in Equation 11, is much above the upper limits from the oscillation experiments.For instance, g ′ a should be no less than ∼ 10 4 eV −1 for the configuration we considered about dark matter components.However it's much larger than the upper bounds from oscillation experiments [38].Taking DUNE as an example, the upper bounds of g ′ a lie in [3 × 10 −12 , 3 × 10 −9 ] eV −1 for m a in [10 −22 , 10 −9 ] eV [38].Thus if the minimal interaction were to persist until the present epoch, it would significantly impact neutrino oscillations, imposing restrictions both on the parameter g ′ a and, consequently, on g a .However, the oscillation limitations would not appear in our case, since the introduced φ has a temperature-dependent vacuum expectation value(see Equation 8).The selected m φ , i.e. m φ = 10 −15 eV, would make the vacuum expectation value of φ to be vanished at around 5 MeV.Thus, at the present time, both the minimal interaction and the effective potential in 9, 10 are absent, as v φ is vanished.The impact on oscillations is also eliminated.As a result, we can get rid of the oscillation constraints.Another simple interaction form between SM neutrino and scalar is φ νν.From Ref. [39] we know it only results in a neutrino mass correction, and is irrelevant to the resonance oscillation.Neutrino coalescence into axion and φ could take place from the effective interaction 6, it may affect the cooling process in supernovae.In order to avoid any discrepancies with the observation of SN1987a, it is necessary for the coupling strength to guarantee that the amount of energy carried away by the new particles does not exceed the value permitted by the neutrino-only scenario.This condition is referred to as the "Raffelt criterion" [40].In essence, this requirement could be met if g a /Λ is much less than the strength of weak interaction, i.e.G F .The resonant production of sterile neutrino and constraints from supernovae cooling necessitate an exceedingly tiny value for λ φ .For instance, with typical pa-rameters as shown in Figure 1, the ratio g a /( λ φ Λ) is set to 10 −14 eV −2 .To satisfy the condition that g a /Λ is much smaller than G F , it implies that λ φ must be much less than 10 −74 .Such a tiny λ φ might seem unnatural, but one potential solution is the Clockwork mechanism [41].By introducing a series of scalar gear fields φ i (i = 1...N), a scalar potential can be constructed as [41] The iterative relationships between these gear fields, determined by the Euler-Lagrange equation of motion, result in φ 1 being exponentially suppressed as φ 1 ∼ q −N ϕ (0) , where ϕ 0 is defined as This allows us to express the quartic interaction in terms of an effective coupling λ φ , which can be naturally generated from an O(1) coupling λ via λ φ ∼ q −4N λ , for q > 1.For instance, choose q = 5, N = 28 and λ at O(1), the resultant λ φ ∼ 10 −80 .To reveal the resonance effects, we depict the variation of the two potentials, as functions of temperature in Figure 1.In practice, we normalize these two terms by ∆ cos 2θ , and momentum is set to p = T .The purple line indicates the relative magnitude of SM potential, and the blue line stands for the relative magnitude of the new potential V a .The black dashed line represents the normalization factor ∆ cos 2θ .
Without considering large deviations between the sterile neutrino momentum and temperature, it can be seen from Figure 1 that ∆ cos 2θ is dominant in the region of T ≲ 0.2 GeV, while the SM potential |V SM | is dominant in the region of T ≳ 0.8 GeV.The new potential V a oscillates with temperature and terminates at T ∼ 5 MeV because the vacuum expectation value vanishes at lower temperatures.As the temperature increases, the oscillation frequency becomes smaller since the relation between time and temperature is t ∝ T −2 .The maximum value that can be reached by the oscillated V a gradually approaches and exceeds ∆ cos 2θ .Therefore, as the temperature increases, the cancellation between the three terms is expected to be more pronounced, until |V SM | increases so that SM contribution becomes dominant.
In Figure 1, we also show the dependence of γ/H on the temperature for both DW (red line) and aDW (green line) mechanisms.The aDW scenario will go back to DW for the high and low temperatures.For the intermediate temperature, the cancellation leads to an enhancement on γ/H in the aDW mechanism, whose magnitude is up to ∼ 6 orders higher in the comparison of the two scenarios.Thus the production of sterile neutrino dark matter is expected to be greatly enhanced for aDW.
Base on the resonance effects, the sterile neutrino could be produced with a much smaller mixing angle.In Figure 2, we show the correct relic density in (m s , sin 2 2θ ) plane.The green line stands for the available parameter space in aDW scenario, and the red line represents the traditional DW mechanism.The most stringent constraints come from the X-ray observation and phase space distribution of dwarf galaxies, which are shown by the brown and purple excluded regions.The most parameter space that predicted by DW mechanism have been excluded.In our scenario containing the new interactions, the values of the mixing angle can be much lower comparatively, meanwhile providing the correct amount of relic abundance.For the parameters we have selected, sterile neutrinos with mass in [2,12] keV can satisfy all the astrophysical constraints.In Figure 2, we also show the detecting capabilities of two forthcoming experiments, i.e.Athena [42] and eXTP [43,44].We see that they could detect some parameter region that is more stringent than the X-ray constraints, and can be used to test the feasibility of our scenario.

CONCLUSION
Sterile neutrino dark matter, generated via the mixing with SM neutrino in DW mechanism, suffers severe conflict between its relic density and the astrophysical observations.Xray constraints excluded the most parameter space of this scenario.To release this tension, we have introduced an effective axial-vector interaction between axion, SM neutrino, and a singlet in this work.The new interactions bring a new contribution to SM neutrino effective potential, then the cancellation occurs between this new contribution and SM one.The sterile neutrino thus can be resonantly produced with a significant enhancement factor.Around the typical DW-peak temperature, the resonance enhancement could produce sterile neutrino with an amount of several orders higher than DW mechanism.Such that the much smaller mixing angles can satisfy the stringent astrophysical constraints and provide the right abundance of the relic DM at the same time.An inevitable consequence of new interaction with neutrino is the modification of the neutrino oscillation.Since the vacuum expectation value of the new scalar dynamically varies in temperature, the new effective potential would vanish at the current temperature.So it's safe to ignore it in our case.

Figure 1 :
Figure 1: Production of sterile neutrino in the DW and aDW mechanisms.The red and green lines stand for the variation of γ/H with temperature T in DW and aDW mechanisms, respectively.The parameters are set as sin 2 2θ = 10 −12 , m s = 10 keV, m φ = 10 −15 eV, ξ = 1/6, m a = 10 −10 eV, and a combination of g a /( λ φ Λ) is set to 10 14 eV −2 .The blue and purple lines show relative size of the SM potential and the new potential to ∆ cos 2θ .

Figure 2 :
Figure 2: Correct relic density of sterile neutrino DM as functions of mass and mixing angle.The red line represents the proper parameter space in the traditional Dodelson-Widrow mechanism (DW), while the green line stands for that in the aDW mechanism (aDW).The choices of model parameters are m φ = 10 −15 eV, ξ = 1/6, m a = 10 −10 eV, and g a /( λ φ Λ) = 10 14 eV −2 .Constraints from X-ray observation are depicted in brown region, and lower bound on mass of sterile neutrino from dwarf galaxies is shown by purple region.The detecting capabilities of two forthcoming experiments, i.e.Athena and eXTP, are showed in orange and blue dashed lines respectively.