REVIEW 3 major objections 6 minor 40 references
Freeze-in leptogenesis with sterile neutrino self-interactions
T0 review · 3 major / 6 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read Scalar-mediated self-interactions among sterile neutrinos, with mass and Yukawa matrices misaligned, can boost the produced baryon asymmetry by orders of magnitude and greatly relax the fine-tuned mass degeneracy required in vanilla ARS…
desk verdict A serious model-building paper with a new structural result — off-diagonal scalar-mediated self-interactions generate the baryon asymmetry at O(F^4) and relax the ARS mass degeneracy — but the scalar-SM thermal equilibrium is assumed without a specified coupling. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the scalar-induced thermal potential $V_\phi = (\pi^2 a_R/(432\,\zeta(3)\,T_{ws}))\,Y\cdot Y^T$, the real part of the sterile-neutrino self-energy from forward scattering on the thermal $\phi$ background. In the quantum kinetic equation it contributes a commutator term $-i[V_\phi,\delta n_N]/2$ that couples off-diagonal sterile correlations to the diagonal charges. Its key effect is a new source term $2\,V_{\phi,ij}\,\mathrm{Im}[\delta n_{N,ij}^{\mathrm{odd}}]$ in the evolution of sterile charges, which generates charges already at first nontrivial order in $F^2$ instead of requiring the $O(F^6)$ washout-driven mechanism of vanilla ARS. The paper also includes the scalar decay and inverse-decay rate $\Gamma_\phi$ and shows it is subdominant in the parameter region where the analytical solution applies.
What would settle it
A direct calculation of the equilibration rate of $\phi$ with the SM plasma would settle the claim: if the only couplings of $\phi$ are to sterile neutrinos, $\phi$ cannot be in equilibrium and the enhancement does not occur as described, while a portal strong enough to maintain equilibrium may itself erase the asymmetry or violate existing constraints on light scalars.
Extended reading notes
Core claim
Working in the mass basis of the two sterile neutrinos, the authors add a real scalar singlet $\phi$ with a real symmetric Yukawa interaction $Y_{ij}\bar N_i N_j \phi$. The scalar's forward scattering produces a thermal potential $V_\phi \propto Y\cdot Y^T$ that acts as an additional contribution to the oscillation Hamiltonian. When $Y$ and the Majorana mass matrix $M$ are diagonal in different bases, this potential mixes the diagonal and off-diagonal components of the sterile density matrix. In the oscillatory regime the authors solve the quantum kinetic equations perturbatively in $F$ and obtain sterile charges at order $O(F^2)$ and a baryon asymmetry at order $O(F^4)$, replacing the $O(F^6)$ suppression of the vanilla ARS mechanism. For a benchmark with $M_1\simeq M_2=1$ GeV, $\Delta M^2=10^{-4}$ GeV$^2$, $Y_{ii}=10^{-6}$, $Y_{ij}=1.1\times 10^{-7}$ and $m_\phi=3$ GeV, the produced asymmetry is $|B|/s=8.76\times 10^{-11}$, matching the observed value while vanilla ARS with the same mass splitting falls short by at least an order of magnitude.
Load-bearing premise
The load-bearing premise is that the new scalar $\phi$ is in thermal equilibrium with the Standard Model plasma, since the thermal potential, the production rate, and the scalar number density all rely on that equilibrium; the paper does not specify the coupling that would maintain it.
Editorial extensions
If this is right
- Sterile-neutrino mass splittings orders of magnitude larger than in vanilla ARS become compatible with the observed baryon asymmetry; the paper shows this for $\Delta M^2$ up to $10^{-2}$ GeV$^2$.
- The baryon asymmetry no longer requires the tiny $O(F^6)$ combination of Standard Model Yukawa couplings; two powers of $F$ are replaced by two powers of the new scalar coupling $Y$, so the scenario can work with smaller active-sterile mixing.
- For large enough $Y_{ij}$, the decay rate brings the sterile neutrinos into equilibrium before sphaleron freeze-out and the asymmetry is washed out, recovering the suppression seen in earlier scalar-extension studies.
- The analytical approximation matches the full numerical result within a factor of about two over its range of validity, giving a closed-form estimate of the produced baryon asymmetry in terms of a generalized hypergeometric function.
- The required mass degeneracy is relaxed to the point that the benchmark point with $\Delta M^2=10^{-4}$ GeV$^2$ and $M\simeq 1$ GeV reproduces $|B|/s\approx 8.76\times 10^{-11}$.
Reading between the lines
- The paper assumes $\phi$ is in thermal equilibrium with the SM plasma but specifies no coupling between $\phi$ and SM particles; if $\phi$ interacts only with sterile neutrinos, its abundance would be produced by freeze-in rather than equilibrium, and the enhancement would need re-evaluation.
- A testable extension would be to compute the minimal scalar-SM portal (for example, Higgs mixing) needed to maintain $\phi$ equilibrium and check against constraints from BBN, the CMB, and collider searches for light scalars; this would map the viable parameter region.
- The generalized case $Y_{ii}\neq Y_{jj}$, which the paper leaves for future work, could produce a qualitatively different thermal-mass splitting and possibly extend the enhancement into regimes where the present approximation breaks down.
- The same scalar-induced potential may affect sterile neutrino dark matter production, so the mechanism could link baryogenesis and dark matter in a minimal two-sterile-neutrino setup.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper investigates whether a scalar singlet φ coupled to two sterile neutrinos via a real symmetric Yukawa matrix can relax the mass-degeneracy condition of ARS freeze-in leptogenesis. The authors add the scalar-mediated forward-scattering potential Vφ and the decay/inverse-decay rates Γφ to the quantum kinetic equations (Eq. (3.6)), solve the system numerically, and derive an analytic approximation in the oscillatory regime. The central claim is that when the scalar Yukawa matrix Y and the sterile mass matrix M are misaligned, the thermal potential generates sterile charges at order F^2, so the final baryon asymmetry scales as F^4 instead of the vanilla F^6, boosting the BAU by more than two orders of magnitude and allowing mass splittings up to ΔM^2 ~ 10^-2 GeV^2 for M ~ 1 GeV.
Significance. If the model is completed, the mechanism is an interesting and original way to alleviate the ARS fine-tuning. The analytic derivation is transparent and yields a distinct parametric scaling (F^4 vs F^6) that can be tested in future numerical studies. The authors provide an explicit benchmark reproducing the observed BAU and a parameter-space region (Fig. 4). However, the quantitative claim is conditional on the unstated assumption that φ is thermalized with the SM plasma, and the numerical solvability of the results is limited by the lack of any description of the numerical method.
major comments (3)
- [Sec. 3 (first paragraph) and Eq. (3.1)] The thermal-equilibrium assumption for φ is not supported by the Lagrangian. The only φ interaction is L_int = Y_ij \bar N_i N_j φ, with Y_ii ~ 10^-6 and Y_ij ~ 10^-7, so φ is coupled only to sterile neutrinos that are themselves populated by freeze-in; no φ-SM operator is present. Consequently, φ would not be in thermal equilibrium with the SM, and the thermal potential Vφ, the production rate Γφ, and n_eq^φ used in Eqs. (3.2), (3.3), (3.6), and (3.7) would be suppressed. The O(F^4) enhancement in Eq. (3.12) and the benchmark in Table 2 rely on this equilibrium bath. The authors should specify a concrete φ-SM coupling that maintains equilibrium, demonstrate that it does not introduce additional washout or modify the QKE, and re-evaluate the parameter space under that model completion.
- [Sec. 3.1 and Table 2] The quantitative results (the factor ≲2 agreement in Fig. 6 and the benchmark |B|/s = 8.76×10^-11) are obtained from numerical solutions of Eq. (3.6), but no details of the numerical method are given: no solver, step-size control, error tolerances, or validation against the Y→0 limit. The numerical results are therefore not reproducible as presented. The authors should describe the integration scheme and provide convergence tests or release the code.
- [Appendix C, Eq. (C.19)] The derivation states that the sterile-neutrino width vanishes because the decay into a heavier scalar is kinematically forbidden. Since mφ = 3 GeV and M1 ≈ M2 = 1 GeV, the process φ → N_i N_j is kinematically allowed and is exactly the production process used in Eq. (3.3). The distinction between the vanishing Γs and the non-vanishing production term Ξ should be clarified, or the derivation corrected; otherwise the new collision terms in Eq. (3.6) are not fully justified.
minor comments (6)
- [Eq. (3.6)] The notation for the new interaction terms appears to be missing division symbols; the terms should read -(z^2/n_eq){Γφ, δn} and -(z^2/n_eq^2) δn Γφ δn. Please clarify.
- [Sec. 3.1, Fig. 2 caption] The dashed line for vanilla ARS is described as "light blue"; in the figure it appears dashed purple. Please check consistency.
- [Eq. (3.2) and Ref. [34]] The thermal potential is taken from the authors' previous work without a derivation; since it is central to the claimed enhancement, a brief derivation or explicit expression in this paper would be helpful.
- [Table 1] The labels "m2_2" and "m2_3" are ambiguous; they should be the solar and atmospheric mass-squared differences.
- [Sec. 4 (Conclusions)] The conclusions mention the strongly overdamped regime of [23] but do not discuss whether the new self-interactions could open other regimes where the approximations fail; a sentence on this would be useful.
- [Eq. (3.7)] The definition of Γφ has different powers of T and n_eq than the thermally averaged production rate in Eq. (3.3); please verify the consistency of the integrated version.
Circularity Check
No significant circularity: the O(F^4) boost is derived from the input Lagrangian's thermal potential and is not fitted to the observed baryon asymmetry; the only self-citation ([34]) supplies a non-load-bearing computation of that potential.
full rationale
The paper's central claim is that the scalar-mediated thermal potential, rather than the vacuum mass splitting, drives sterile neutrino oscillations and generates sterile charges at O(F^2), leading to a baryon asymmetry at O(F^4). This scaling follows from the analytic solution of the quantum kinetic equations: Eq. (3.12) gives q_Ni = V_phi Im[F*F^T] G times a hypergeometric factor, with V_phi defined in Eq. (3.7) from the input Yukawa matrix Y. The statement that 'the sterile potential is able to generate charges already at order O(F^2)' is a derived consequence of solving Eqs. (3.10) and (3.11), not an assumption. No free parameter is fitted to the observed baryon-to-entropy ratio: the benchmark point in Table 2 is one point on a parameter scan whose output |B|/s is computed and then compared with the observed band; the observed value appears only as a reference line in Figs. 2-4. The paper's assumption that 'phi is in thermal equilibrium with the SM plasma' is an explicit model idealization, and the Lagrangian (3.1) indeed contains no phi-SM coupling, so the model is incomplete as stated. That is a physical robustness limitation, not a circularity: the output asymmetry is not built into the input assumption. The only self-citation is Ref. [34] for the computation of the thermal potential H_th_phi in Eq. (3.2); that is a standard self-energy computation from the same interaction, and Appendix C independently derives the QKE structure of the new terms. The self-citation is therefore not load-bearing for the central O(F^4) result, and the derivation is self-contained against external benchmarks. Thus the appropriate circularity score is low, reflecting one minor non-load-bearing self-citation rather than any reduction of the prediction to its inputs.
Assumptions & free parameters
free parameters (6)
- Yii (diagonal scalar Yukawa coupling) =
10^-6 or 10^-7
- Yij (off-diagonal scalar Yukawa coupling) =
1.1 x 10^-7 benchmark; varied from 10^-10 to 10^-5
- m_phi (scalar mass) =
3 GeV
- Delta M^2 (sterile neutrino mass-squared difference) =
10^-4 GeV^2 benchmark; varied
- M1, M2 (sterile neutrino masses) =
1 GeV
- Casas-Ibarra parameters (omega, delta, alpha1, alpha2) =
omega = 3 pi/4 + 2.16 i, delta = 3 pi/2, alpha1 = 0, alpha2 = -2 pi
assumptions (6)
- domain assumption The standard ARS quantum kinetic equation (2.7) with the simplifying assumptions leading to Eq. (2.11) is an adequate description of the system.
- domain assumption The scalar phi remains in thermal equilibrium with the SM plasma throughout the relevant epoch.
- domain assumption The scalar is sufficiently light, m_phi << T, so the thermal potential takes the form in Eq. (3.2).
- domain assumption Perturbative expansion in the small SM Yukawa coupling F, with the replacement (delta n_N+/-)_ij -> -neq delta_ij on the right-hand side, captures the leading asymmetry.
- ad hoc to paper The scalar Yukawa matrix Y is real and symmetric and is not diagonal in the same basis as the sterile mass matrix M.
- domain assumption The Majorana width of the sterile neutrinos from the new interaction vanishes because N_i -> N_j phi is kinematically forbidden.
invented entities (1)
-
Scalar singlet phi
Cite this review
Pith. "Pith review of Freeze-in leptogenesis with sterile neutrino self-interactions." pith.science (2026). https://pith.science/paper/YKRNKYNW
@misc{pith2026241214044,
author = {Pith},
title = {Pith review of: Freeze-in leptogenesis with sterile neutrino self-interactions},
year = {2026},
howpublished = {\url{https://pith.science/paper/YKRNKYNW}},
note = {Machine review of arXiv:2412.14044}
}
read the original abstract
Sterile neutrinos are a simple yet compelling addition to the Standard Model. For right-handed neutrinos with masses below the electroweak scale, leptogenesis can proceed through CP-violating oscillations of the sterile neutrinos. This is known as ARS or freeze-in leptogenesis. However, the ARS scenario requires the right-handed neutrinos to have a high degree of mass degeneracy. In this work, we study an extension of the SM that introduces a scalar singlet in addition to the two sterile neutrinos required to generate the baryon asymmetry. The new scalar interacts with the sterile neutrinos via a Yukawa interaction. This leads to an additional rate for the production and destruction of the sterile neutrinos and to a novel contribution to the effective potential. For the case in which the mass and the new Yukawa matrices are not diagonal in the same basis, we find that the effective potential can boost the baryon asymmetry of the universe by several orders of magnitude. This significantly alleviates the fine-tuned mass condition required in vanilla ARS leptogenesis.
Reference graph
Works this paper leans on
-
[1]
Aghanim et al.,Planck 2018 results
Planck Collaboration, N. Aghanim et al.,Planck 2018 results. VI. Cosmological parameters, Astron. Astrophys.641 (2020) A6, [arXiv:1807.06209]. [Erratum: Astron.Astrophys. 652, C4 (2021)]. – 20 –
arXiv 2020
-
[2]
Navas et al.,Review of particle physics, Phys
Particle Data Group Collaboration, S. Navas et al.,Review of particle physics, Phys. Rev. D 110 (2024), no. 3 030001
work page 2024
-
[3]
T.-H. Yeh, J. Shelton, K. A. Olive, and B. D. Fields,Probing physics beyond the standard model: limits from BBN and the CMB independently and combined, JCAP 10 (2022) 046, [arXiv:2207.13133]
arXiv 2022
-
[4]
M. Fukugita and T. Yanagida,Baryogenesis Without Grand Unification, Phys. Lett. B174 (1986) 45–47
work page 1986
-
[5]
V. A. Kuzmin, V. A. Rubakov, and M. E. Shaposhnikov,On the Anomalous Electroweak Baryon Number Nonconservation in the Early Universe, Phys. Lett. B155 (1985) 36
work page 1985
-
[6]
T. Asaka and M. Shaposhnikov,The νMSM, dark matter and baryon asymmetry of the universe, Phys. Lett. B620 (2005) 17–26, [hep-ph/0505013]
arXiv 2005
- [7]
-
[8]
The nuMSM, leptonic asymmetries, and properties of singlet fermions
M. Shaposhnikov,The nuMSM, leptonic asymmetries, and properties of singlet fermions, JHEP 08 (2008) 008, [arXiv:0804.4542]
work page Pith review arXiv 2008
Show all 40 references
-
[9]
Asaka and H
T. Asaka and H. Ishida,Flavour Mixing of Neutrinos and Baryon Asymmetry of the Universe, Phys. Lett. B692 (2010) 105–113, [arXiv:1004.5491]
2010 arXiv
-
[10]
Canetti and M
L. Canetti and M. Shaposhnikov,Baryon Asymmetry of the Universe in the NuMSM, JCAP 09 (2010) 001, [arXiv:1006.0133]
2010 arXiv
-
[11]
Pilaftsis,CP violation and baryogenesis due to heavy Majorana neutrinos, Phys
A. Pilaftsis,CP violation and baryogenesis due to heavy Majorana neutrinos, Phys. Rev. D56 (1997) 5431–5451, [hep-ph/9707235]
1997 arXiv
-
[12]
Pilaftsis and T
A. Pilaftsis and T. E. J. Underwood,Resonant leptogenesis, Nucl. Phys. B692 (2004) 303–345, [hep-ph/0309342]
2004 arXiv
-
[13]
E. K. Akhmedov, V. A. Rubakov, and A. Y. Smirnov,Baryogenesis via neutrino oscillations, Phys. Rev. Lett.81 (1998) 1359–1362, [hep-ph/9803255]
1998 arXiv
-
[14]
Drewes, B
M. Drewes, B. Garbrecht, P. Hernandez, M. Kekic, J. Lopez-Pavon, J. Racker, N. Rius, J. Salvado, and D. Teresi,ARS Leptogenesis, Int. J. Mod. Phys. A33 (2018), no. 05n06 1842002, [arXiv:1711.02862]
2018 arXiv
-
[15]
Beacham et al.,Physics Beyond Colliders at CERN: Beyond the Standard Model Working Group Report, J
J. Beacham et al.,Physics Beyond Colliders at CERN: Beyond the Standard Model Working Group Report, J. Phys. G47 (2020), no. 1 010501, [arXiv:1901.09966]
2020 arXiv
-
[16]
Flood, R
I. Flood, R. Porto, J. Schlesinger, B. Shuve, and M. Thum,Hidden-sector neutrinos and freeze-in leptogenesis, Phys. Rev. D105 (2022), no. 9 095025, [arXiv:2109.10908]
2022 arXiv
-
[17]
Fischer, M
O. Fischer, M. Lindner, and S. van der Woude,Robustness of ARS leptogenesis in scalar extensions, JHEP 05 (2022) 149, [arXiv:2110.14499]
2022 arXiv
-
[18]
Escudero and S
M. Escudero and S. J. Witte,The hubble tension as a hint of leptogenesis and neutrino mass generation, Eur. Phys. J. C81 (2021), no. 6 515, [arXiv:2103.03249]
2021 arXiv
-
[19]
Caputo, P
A. Caputo, P. Hernandez, and N. Rius,Leptogenesis from oscillations and dark matter, Eur. Phys. J. C79 (2019), no. 7 574, [arXiv:1807.03309]
2019 arXiv
-
[20]
Alanne, T
T. Alanne, T. Hugle, M. Platscher, and K. Schmitz,Low-scale leptogenesis assisted by a real scalar singlet, JCAP 03 (2019) 037, [arXiv:1812.04421]
2019 arXiv
-
[21]
A. D. Sakharov,Violation of CP Invariance, C asymmetry, and baryon asymmetry of the universe, Pisma Zh. Eksp. Teor. Fiz.5 (1967) 32–35
1967
-
[22]
Hambye and D
T. Hambye and D. Teresi,Baryogenesis from L-violating Higgs-doublet decay in the density-matrix formalism, Phys. Rev. D96 (2017), no. 1 015031, [arXiv:1705.00016]. – 21 –
2017 arXiv
-
[23]
Drewes, B
M. Drewes, B. Garbrecht, D. Gueter, and J. Klaric,Leptogenesis from Oscillations of Heavy Neutrinos with Large Mixing Angles, JHEP 12 (2016) 150, [arXiv:1606.06690]
2016 arXiv
-
[24]
Hernández, M
P. Hernández, M. Kekic, J. López-Pavón, J. Racker, and N. Rius,Leptogenesis in GeV scale seesaw models, JHEP 10 (2015) 067, [arXiv:1508.03676]
2015 arXiv
-
[25]
Hernández, M
P. Hernández, M. Kekic, J. López-Pavón, J. Racker, and J. Salvado,Testable Baryogenesis in Seesaw Models, JHEP 08 (2016) 157, [arXiv:1606.06719]
2016 arXiv
-
[26]
Eijima, M
S. Eijima, M. Shaposhnikov, and I. Timiryasov,Parameter space of baryogenesis in theνMSM, JHEP 07 (2019) 077, [arXiv:1808.10833]
2019 arXiv
-
[27]
Drewes and B
M. Drewes and B. Garbrecht,Leptogenesis from a GeV Seesaw without Mass Degeneracy, JHEP 03 (2013) 096, [arXiv:1206.5537]
2013 arXiv
-
[28]
A. D. Dolgov,Neutrinos in the Early Universe, Sov. J. Nucl. Phys.33 (1981) 700–706
1981
-
[29]
Sigl and G
G. Sigl and G. Raffelt,General kinetic description of relativistic mixed neutrinos, Nucl. Phys. B 406 (1993) 423–451
1993
-
[30]
Barbieri and A
R. Barbieri and A. Dolgov,Neutrino oscillations in the early universe, Nucl. Phys. B349 (1991) 743–753
1991
-
[31]
Asaka, S
T. Asaka, S. Eijima, and H. Ishida,Kinetic Equations for Baryogenesis via Sterile Neutrino Oscillation, JCAP 02 (2012) 021, [arXiv:1112.5565]
2012 arXiv
-
[32]
Notzold and G
D. Notzold and G. Raffelt,Neutrino Dispersion at Finite Temperature and Density, Nucl. Phys. B 307 (1988) 924–936
1988
-
[33]
Quimbay and S
C. Quimbay and S. Vargas-Castrillon,Fermionic dispersion relations in the standard model at finite temperature, Nucl. Phys. B451 (1995) 265–304, [hep-ph/9504410]
1995 arXiv
-
[34]
M. D. Astros and S. Vogl,Boosting the production of sterile neutrino dark matter with self-interactions, JHEP 03 (2024) 032, [arXiv:2307.15565]
2024 arXiv
-
[35]
J. A. Casas and A. Ibarra,Oscillating neutrinos and µ → e, γ, Nucl. Phys. B618 (2001) 171–204, [hep-ph/0103065]
2001 arXiv
-
[36]
M. C. Gonzalez-Garcia, M. Maltoni, and T. Schwetz,NuFIT: Three-Flavour Global Analyses of Neutrino Oscillation Experiments, Universe 7 (2021), no. 12 459, [arXiv:2111.03086]
2021 arXiv
-
[37]
Drewes and B
M. Drewes and B. Garbrecht,Combining experimental and cosmological constraints on heavy neutrinos, Nucl. Phys. B921 (2017) 250–315, [arXiv:1502.00477]
2017 arXiv
-
[38]
Abada, G
A. Abada, G. Arcadi, V. Domcke, and M. Lucente,Lepton number violation as a key to low-scale leptogenesis, JCAP 11 (2015) 041, [arXiv:1507.06215]
2015 arXiv
-
[39]
Brambilla, M
N. Brambilla, M. A. Escobedo, J. Soto, and A. Vairo,Heavy quarkonium suppression in a fireball, Phys. Rev. D97 (2018), no. 7 074009, [arXiv:1711.04515]
2018 arXiv
-
[40]
M. L. Bellac,Thermal Field Theory. Cambridge Monographs on Mathematical Physics. Cambridge University Press, 3, 2011. – 22 –
2011
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.