pith. machine review for the scientific record. sign in

arxiv: 2605.12946 · v1 · submitted 2026-05-13 · ✦ hep-ph · astro-ph.CO

Recognition: unknown

Minimal Majoron Dark Matter

Authors on Pith no claims yet

Pith reviewed 2026-05-14 18:20 UTC · model grok-4.3

classification ✦ hep-ph astro-ph.CO
keywords Majorondark matterType-I seesawfreeze-inmisalignment mechanismleptogenesisneutrino massesscalar singlet
0
0 comments X

The pith

The Majoron dark matter mass is bounded by about 10 MeV without fine-tuning the initial misalignment angle.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

In the minimal extension of the Standard Model with right-handed neutrinos and a complex scalar singlet, the Majoron emerges as a candidate for dark matter. The paper calculates its relic density by considering production from the thermal plasma via freeze-in and from the oscillating scalar field via the misalignment mechanism. A key result is that matching the observed dark matter abundance without carefully choosing the starting value of the field angle restricts the Majoron mass to no more than order 10 MeV. This constraint arises because higher masses would require either overproduction or underproduction that cannot be adjusted naturally. The analysis also explores how this dark matter scenario can be consistent with generating the matter-antimatter asymmetry through leptogenesis.

Core claim

In the minimal Majoron dark matter model based on the Type-I seesaw with a complex scalar, the combination of freeze-in and misalignment production mechanisms implies that the Majoron mass satisfies m_J ≲ O(10) MeV without fine-tuning of the initial misalignment angle θ_i. For compatibility with thermal leptogenesis involving two right-handed neutrinos, the misalignment-dominated case allows m_J ≲ O(100) eV, while freeze-in dominance requires θ_i ≲ O(0.01).

What carries the argument

The Majoron, the pseudo-Nambu-Goldstone boson from lepton number breaking, whose abundance is set by freeze-in from interactions with the thermal bath and by the misalignment mechanism from the initial field displacement.

If this is right

  • The Majoron can constitute all dark matter only for masses up to O(10) MeV absent fine-tuning of the initial misalignment angle.
  • Thermal leptogenesis with two right-handed neutrinos is compatible either with misalignment-dominated production at m_J ≲ O(100) eV or with mild tuning θ_i ≲ O(0.01) when freeze-in dominates.
  • Observational constraints on dark matter density and stability further restrict the allowed couplings between the Majoron and the Standard Model particles.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Future light dark matter searches sensitive to MeV-scale particles could directly test the upper mass bound derived here.
  • If an independent measurement fixes the Majoron mass above a few MeV, the model would require either additional production channels or acceptance of some initial-condition tuning.
  • The same production mechanisms and mass bounds may apply to other pseudo-Goldstone bosons in minimal extensions that break global symmetries at the seesaw scale.

Load-bearing premise

The Majoron is stable on cosmological timescales, constitutes all of the dark matter, and is produced solely through freeze-in and misalignment in the minimal Type-I seesaw extension.

What would settle it

A confirmed Majoron mass above O(10) MeV that still accounts for the full observed dark matter density without requiring a precisely tuned initial misalignment angle would contradict the central bound.

Figures

Figures reproduced from arXiv: 2605.12946 by Haruto Kitagawa, Kensuke Akita, Koichi Hamaguchi, Tatsuya Yokoyama.

Figure 1
Figure 1. Figure 1: Evolution of the Majoron abundance YJ = nJ /s from freeze-in production for M1 = M2 = 108 GeV, with g = 10−2 (f ≃ 1.4 × 1010 GeV) in the left panel and g = 10−4 (f ≃ 1.4 × 1012 GeV) in the right panel. The individual contributions from the different freeze-in production channels, as well as their sum, are shown separately. where xi = s/M2 i and βi = (1 − 4M2 i /s) 1/2 with s being the Mandelstam variable. … view at source ↗
Figure 2
Figure 2. Figure 2: Values of g required to reproduce the observed DM abundance through freeze-in production as a function of the Majoron mass mJ , for M1 = M2 ≡ M and z = 0. The blue, orange, and green curves correspond to M = 106 , 108 , and 1010 GeV, respectively. T ∼ M, one finds YJ ∼ σn2 NiH−1/s ∼ g 4/M, and hence Ω (FI) J h 2 ∝ mJ YJ ∝ mJ g 4/M. In this regime, we numerically find Ω (FI) J h 2 ≃ 0.26  mJ 10−5 GeV   M… view at source ↗
Figure 3
Figure 3. Figure 3: Parameter space in the (mJ , f) plane for θi = 0.1 (left) and θi = 0.01 (right), with M1 = M2 = M and z = 0. The colored contours show the parameter regions reproducing the observed DM abundance, Ω (FI) J h 2 + Ω(mis) J h 2 = ΩDMh 2 ≃ 0.12 for M = 106 GeV (blue), 108 GeV (orange), and 1010 GeV (green). For each M, the contour has two branches: the upper branch is dominated by misalignment production and th… view at source ↗
Figure 4
Figure 4. Figure 4: Values of various combinations of K required to reproduce the observed DM abundance as a function of mJ , with M1 = M2 = M and z = 0, for θi = 0.1 (left) and θi = 0.01 (right). The blue and orange curves correspond to M = 106 GeV and 108 GeV, respectively, following the same convention as in [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Values of f required to reproduce the observed DM abundance via the misalignment mech￾anism (navy curves, for θi = 0.1, 0.01, and 0.001 from bottom to top) and via freeze-in production (colored curves) as a function of mJ , with M2 = 10M1 and TR = M2. The blue, orange, and green curves correspond to the three benchmark points with M1 = 1012 , 1011, and 2.7 × 1010 GeV, respec￾tively, for which successful le… view at source ↗
Figure 6
Figure 6. Figure 6: Dependence of tr K and | − Kee + Kµµ + Kττ | on the Casas-Ibarra parameter z = a + ib and the Majorana phase α, for M2/M1 = 1 (left) and M2/M1 = 10 (right). The bands show the range obtained by varying a ∈ [0, π] and α ∈ [0, 2π] for each value of b, normalized to the benchmark case z = α = 0. Using Eq. (7), one finds vfK = mDm † D = U   0 0 0 0 ˆmν 2 [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Dependence of K′ on the Casas-Ibarra parameter z = a + ib and the Majorana phase α as a function of mJ , for M2/M1 = 1 (top) and M2/M1 = 10 (bottom), with cosh 2b = 1 (left) and cosh 2b = 10 (right). The bands show the range obtained by varying a and α for each value of mJ , normalized to the benchmark case z = α = 0. values of mJ , for M2/M1 = 1, 10 (top and bottom) and cosh 2b = 1, 10 (left and right). H… view at source ↗
read the original abstract

We study Majoron dark matter (DM) in its minimal realization, based on the Type-I seesaw framework extended by a SM-singlet complex scalar. Remaining agnostic about the origin and value of the Majoron mass, we evaluate the DM abundance from both the freeze-in and misalignment mechanisms, and identify the viable parameter space consistent with observational constraints. Without fine-tuning of the initial misalignment angle, we find that the Majoron mass is bounded by $m_J \lesssim \mathcal{O}(10)~\mathrm{MeV}$. We also discuss compatibility with thermal leptogenesis. Successful leptogenesis with two right-handed neutrinos favors misalignment-dominated production with the Majoron mass $m_J \lesssim \mathcal{O}(100)~\mathrm{eV}$, while freeze-in dominated production is compatible with leptogenesis only with a mild fine-tuning of the initial misalignment angle, $\theta_i \lesssim \mathcal{O}(0.01)$.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript examines Majoron dark matter in the minimal Type-I seesaw framework extended by a SM-singlet complex scalar. It calculates the dark matter relic abundance from freeze-in and misalignment mechanisms while remaining agnostic about the origin of the Majoron mass. The central result is that, without fine-tuning of the initial misalignment angle, the Majoron mass satisfies m_J ≲ O(10) MeV, with additional analysis of compatibility with thermal leptogenesis for two right-handed neutrinos (misalignment-dominated production allowing m_J ≲ O(100) eV, and freeze-in requiring mild θ_i fine-tuning).

Significance. If the production calculations hold, the work supplies useful bounds on a minimal Majoron DM candidate and its interplay with leptogenesis. The agnostic treatment of the mass origin broadens the applicability of the results. The explicit linkage between DM production channels and successful leptogenesis is a strength that could inform model-building in neutrino and dark-matter phenomenology.

major comments (2)
  1. [§3.2] §3.2 (Misalignment Production): The relic-density calculation assumes standard axion-like temperature-dependent mass activation and oscillation onset. Because the paper remains agnostic about the Majoron mass origin, the mass term could arise from explicit breaking or operators that activate near the electroweak scale rather than a high-scale phase transition; this would shift the Hubble crossing and change the misalignment abundance by orders of magnitude, directly affecting the headline bound m_J ≲ O(10) MeV for θ_i ~ O(1).
  2. [§5] §5 (Leptogenesis compatibility): The statement that freeze-in-dominated production is compatible with leptogenesis only for mild fine-tuning θ_i ≲ O(0.01) lacks an explicit quantitative measure of fine-tuning (e.g., Barbieri-Giudice measure) and a clear mapping of the overlapping parameter space with the two-RH-neutrino leptogenesis constraints; this is load-bearing for the compatibility claim.
minor comments (2)
  1. [Abstract] Abstract: The phrase 'observational constraints' is used without enumerating the specific bounds applied (e.g., CMB, BBN, Lyman-α); adding a short list would improve readability.
  2. [Notation] Notation throughout: The approximate symbols O(10) MeV and O(100) eV should be accompanied by the precise numerical upper limits obtained from the scans or analytic expressions for transparency.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment point by point below, indicating the revisions made.

read point-by-point responses
  1. Referee: §3.2 (Misalignment Production): The relic-density calculation assumes standard axion-like temperature-dependent mass activation and oscillation onset. Because the paper remains agnostic about the Majoron mass origin, the mass term could arise from explicit breaking or operators that activate near the electroweak scale rather than a high-scale phase transition; this would shift the Hubble crossing and change the misalignment abundance by orders of magnitude, directly affecting the headline bound m_J ≲ O(10) MeV for θ_i ~ O(1).

    Authors: We agree that the calculation relies on the standard high-scale mass activation assumption common to axion-like misalignment. Given our agnostic stance on the mass origin, we have revised §3.2 to explicitly state this assumption and to note that low-scale activation (e.g., near the electroweak scale) would alter the Hubble crossing and relic abundance, potentially relaxing the bound. The headline result m_J ≲ O(10) MeV is presented under the standard high-scale case, which is the typical scenario in the minimal model. revision: yes

  2. Referee: §5 (Leptogenesis compatibility): The statement that freeze-in-dominated production is compatible with leptogenesis only for mild fine-tuning θ_i ≲ O(0.01) lacks an explicit quantitative measure of fine-tuning (e.g., Barbieri-Giudice measure) and a clear mapping of the overlapping parameter space with the two-RH-neutrino leptogenesis constraints; this is load-bearing for the compatibility claim.

    Authors: We acknowledge the need for greater precision. In the revised manuscript we have added the Barbieri-Giudice fine-tuning measure applied to θ_i, confirming that θ_i ≲ O(0.01) corresponds to mild tuning of order 10^{-2}. We have also included an explicit mapping of the overlapping parameter space (with an accompanying figure) between the DM production channels and the viable region for successful thermal leptogenesis with two right-handed neutrinos, thereby strengthening the compatibility statement. revision: yes

Circularity Check

0 steps flagged

No circularity: bounds derived from external cosmology and observations

full rationale

The paper evaluates Majoron DM abundance via standard freeze-in and misalignment formulas applied to the Type-I seesaw extension, remaining explicitly agnostic about mass origin. The central bound m_J ≲ O(10) MeV follows from requiring the observed relic density without θ_i fine-tuning, using external Hubble evolution and observational inputs rather than any internal fit, self-definition, or self-citation chain. No step reduces the target result to the paper's own parameters by construction.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 1 invented entities

The claim rests on the Type-I seesaw extension, standard cosmological production channels, and the assumption that the Majoron is the sole dark-matter component; no new free parameters are introduced beyond the Majoron mass and initial angle already discussed.

free parameters (2)
  • Majoron mass m_J
    Remains agnostic about its origin; bounds are derived rather than fitted.
  • initial misalignment angle theta_i
    Discussed with and without fine-tuning; values such as theta_i ≲ O(0.01) are stated as conditions.
axioms (2)
  • domain assumption Type-I seesaw framework extended by SM-singlet complex scalar
    Basis for the minimal Majoron model stated in the abstract.
  • domain assumption Standard freeze-in and misalignment production mechanisms
    Used to evaluate DM abundance without derivation in the abstract.
invented entities (1)
  • Majoron no independent evidence
    purpose: Dark matter candidate arising from spontaneous symmetry breaking
    Postulated in the minimal extension; no independent evidence supplied beyond the model itself.

pith-pipeline@v0.9.0 · 5455 in / 1377 out tokens · 22423 ms · 2026-05-14T18:20:22.695107+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

50 extracted references · 40 canonical work pages · 23 internal anchors

  1. [1]

    Minkowski,µ→eγat a Rate of One Out of10 9 Muon Decays?,Phys

    P. Minkowski,µ→eγat a Rate of One Out of10 9 Muon Decays?,Phys. Lett. B67(1977) 421

  2. [2]

    Yanagida,Horizontal gauge symmetry and masses of neutrinos,Conf

    T. Yanagida,Horizontal gauge symmetry and masses of neutrinos,Conf. Proc. C7902131 (1979) 95

  3. [3]

    Complex Spinors and Unified Theories

    M. Gell-Mann, P. Ramond and R. Slansky,Complex Spinors and Unified Theories,Conf. Proc. C790927(1979) 315 [1306.4669]

  4. [4]

    R. N. Mohapatra and G. Senjanovic,Neutrino Mass and Spontaneous Parity Nonconservation, Phys. Rev. Lett.44(1980) 912

  5. [5]

    Schechter and J

    J. Schechter and J. W. F. Valle,Neutrino Masses in SU(2) x U(1) Theories,Phys. Rev. D22 (1980) 2227. 15

  6. [6]

    Fukugita and T

    M. Fukugita and T. Yanagida,Baryogenesis Without Grand Unification,Phys. Lett. B174 (1986) 45

  7. [7]

    Chikashige, R

    Y. Chikashige, R. N. Mohapatra and R. D. Peccei,Are There Real Goldstone Bosons Associated with Broken Lepton Number?,Phys. Lett. B98(1981) 265

  8. [8]

    G. B. Gelmini and M. Roncadelli,Left-Handed Neutrino Mass Scale and Spontaneously Broken Lepton Number,Phys. Lett. B99(1981) 411

  9. [9]

    Schechter and J

    J. Schechter and J. W. F. Valle,Neutrino decay and spontaneous violation of lepton number, Phys. Rev. D25(1982) 774

  10. [10]

    S. B. Giddings and A. Strominger,STRING WORMHOLES,Phys. Lett. B230(1989) 46

  11. [11]

    Abbott and M

    L. Abbott and M. B. Wise,Wormholes and global symmetries,Nucl. Phys. B325(1989) 687

  12. [12]

    L. J. Hall, K. Jedamzik, J. March-Russell and S. M. West,Freeze-In Production of FIMP Dark Matter,JHEP03(2010) 080 [0911.1120]

  13. [13]

    I. Z. Rothstein, K. S. Babu and D. Seckel,Planck scale symmetry breaking and majoron physics, Nucl. Phys. B403(1993) 725 [hep-ph/9301213]

  14. [14]

    The KeV Majoron as Dark Matter Particle

    V. Berezinsky and J. W. F. Valle,The KeV majoron as a dark matter particle,Phys. Lett. B 318(1993) 360 [hep-ph/9309214]

  15. [15]

    Decaying warm dark matter and neutrino masses

    M. Lattanzi and J. W. F. Valle,Decaying warm dark matter and neutrino masses,Phys. Rev. Lett.99(2007) 121301 [0705.2406]

  16. [16]

    X-ray photons from late-decaying majoron dark matter

    F. Bazzocchi, M. Lattanzi, S. Riemer-Sørensen and J. W. F. Valle,X-ray photons from late-decaying majoron dark matter,JCAP08(2008) 013 [0805.2372]

  17. [17]

    P.-H. Gu, E. Ma and U. Sarkar,Pseudo-Majoron as Dark Matter,Phys. Lett. B690(2010) 145 [1004.1919]

  18. [18]

    Updated CMB, X- and gamma-ray constraints on majoron dark matter

    M. Lattanzi, S. Riemer-Sorensen, M. Tortola and J. W. F. Valle,Updated CMB and x- andγ-ray constraints on Majoron dark matter,Phys. Rev. D88(2013) 063528 [1303.4685]

  19. [19]

    F. S. Queiroz and K. Sinha,The Poker Face of the Majoron Dark Matter Model: LUX to keV Line,Phys. Lett. B735(2014) 69 [1404.1400]

  20. [20]

    Sub-GeV dark matter as pseudo-Goldstone from the seesaw scale

    M. Frigerio, T. Hambye and E. Masso,Sub-GeV dark matter as pseudo-Goldstone from the seesaw scale,Phys. Rev. X1(2011) 021026 [1107.4564]

  21. [21]

    Cold light dark matter in extended seesaw models

    S. Boulebnane, J. Heeck, A. Nguyen and D. Teresi,Cold light dark matter in extended seesaw models,JCAP04(2018) 006 [1709.07283]

  22. [22]

    Cold keV dark matter from decays and scatterings

    J. Heeck and D. Teresi,Cold keV dark matter from decays and scatterings,Phys. Rev. D96 (2017) 035018 [1706.09909]

  23. [23]

    Y. Abe, Y. Hamada, T. Ohata, K. Suzuki and K. Yoshioka,TeV-scale Majorogenesis,JHEP07 (2020) 105 [2004.00599]

  24. [24]

    S. K. Manna and A. Sil,Majorons revisited: Light dark matter as a FIMP,Phys. Rev. D108 (2023) 075026 [2212.08404]

  25. [25]

    X.-J. Xu, S. Zhou and J. Zhu,TheνR-philic scalar dark matter,JCAP04(2024) 012 [2310.16346]. 16

  26. [26]

    Chiang, S.-Y

    C.-W. Chiang, S.-Y. Ho and V. Q. Tran,TeV-scale unification of light dark matter and neutrino mass,2603.05200

  27. [27]

    Neutrino Lines from Majoron Dark Matter

    C. Garcia-Cely and J. Heeck,Neutrino Lines from Majoron Dark Matter,JHEP05(2017) 102 [1701.07209]

  28. [28]

    Akita and M

    K. Akita and M. Niibo,Updated constraints and future prospects on majoron dark matter,JHEP 07(2023) 132 [2304.04430]

  29. [29]

    Liang, X

    Q. Liang, X. Ponce Díaz and T. T. Yanagida,Axion Detection Experiments Can Probe Majoron Models,Phys. Rev. Lett.134(2025) 151803 [2406.19083]

  30. [30]

    Probing Majoron Dark Matter with Gravitational Wave Detectors

    I. Obata and T. T. Yanagida,Probing Majoron Dark Matter with Gravitational Wave Detectors, 2604.08193

  31. [31]

    J. A. Casas and A. Ibarra,Oscillating neutrinos andµ→e, γ,Nucl. Phys. B618(2001) 171 [hep-ph/0103065]. [32]Particle Data Groupcollaboration,Review of particle physics,Phys. Rev. D110(2024) 030001

  32. [32]

    NuFit-6.0: Updated global analysis of three-flavor neutrino oscillations

    I. Esteban, M. C. Gonzalez-Garcia, M. Maltoni, I. Martinez-Soler, J. P. Pinheiro and T. Schwetz,NuFit-6.0: updated global analysis of three-flavor neutrino oscillations,JHEP12 (2024) 216 [2410.05380]

  33. [33]

    Leptogenesis for Pedestrians

    W. Buchmuller, P. Di Bari and M. Plumacher,Leptogenesis for pedestrians,Annals Phys.315 (2005) 305 [hep-ph/0401240]

  34. [34]

    Villasenor, B

    B. Villasenor, B. Robertson, P. Madau and E. Schneider,New constraints on warm dark matter from the Lyman-αforest power spectrum,Phys. Rev. D108(2023) 023502 [2209.14220]

  35. [35]

    N. F. Bell, M. J. Dolan and S. Robles,Searching for Sub-GeV Dark Matter in the Galactic Centre using Hyper-Kamiokande,JCAP09(2020) 019 [2005.01950]

  36. [36]

    Akita, G

    K. Akita, G. Lambiase, M. Niibo and M. Yamaguchi,Neutrino lines from MeV dark matter annihilation and decay in JUNO,JCAP10(2022) 097 [2206.06755]

  37. [37]

    Astrophysical and Terrestrial Constraints on Singlet Majoron Models

    A. Pilaftsis,Astrophysical and terrestrial constraints on singlet Majoron models,Phys. Rev. D 49(1994) 2398 [hep-ph/9308258]

  38. [38]

    Heeck and H

    J. Heeck and H. H. Patel,Majoron at two loops,Phys. Rev. D100(2019) 095015 [1909.02029]. [40]Fermi-LATcollaboration,Updated search for spectral lines from Galactic dark matter interactions with pass 8 data from the Fermi Large Area Telescope,Phys. Rev. D91(2015) 122002 [1506.00013]

  39. [39]

    Fischer, D

    S. Fischer, D. Malyshev, L. Ducci and A. Santangelo,New constraints on decaying dark matter from INTEGRAL/SPI,Mon. Not. Roy. Astron. Soc.520(2023) 6322 [2211.06200]

  40. [40]

    Constraining Light Dark Matter with Diffuse X-Ray and Gamma-Ray Observations

    R. Essig, E. Kuflik, S. D. McDermott, T. Volansky and K. M. Zurek,Constraining Light Dark Matter with Diffuse X-Ray and Gamma-Ray Observations,JHEP11(2013) 193 [1309.4091]

  41. [41]

    Caputo, M

    A. Caputo, M. Negro, M. Regis and M. Taoso,Dark matter prospects with COSI: ALPs, PBHs and sub-GeV dark matter,JCAP02(2023) 006 [2210.09310]. 17

  42. [42]

    J. A. Tomsick et al.,The Compton Spectrometer and Imager,PoSICRC2023(2023) 745 [2308.12362]

  43. [43]

    K. E. O’Donnell and T. R. Slatyer,Constraints on dark matter with future MeV gamma-ray telescopes,Phys. Rev. D111(2025) 083037 [2411.00087]

  44. [44]

    C. A. Argüelles, D. Delgado, A. Friedlander, A. Kheirandish, I. Safa, A. C. Vincent et al.,Dark matter decay to neutrinos,Phys. Rev. D108(2023) 123021 [2210.01303]

  45. [45]

    A lower bound on the right-handed neutrino mass from leptogenesis

    S. Davidson and A. Ibarra,A Lower bound on the right-handed neutrino mass from leptogenesis, Phys. Lett. B535(2002) 25 [hep-ph/0202239]

  46. [46]

    P. H. Chankowski and K. Turzynski,Limits on T(reh) for thermal leptogenesis with hierarchical neutrino masses,Phys. Lett. B570(2003) 198 [hep-ph/0306059]

  47. [47]

    Neutrino Masses, Lepton Flavor Mixing and Leptogenesis in the Minimal Seesaw Model

    W.-l. Guo, Z.-z. Xing and S. Zhou,Neutrino Masses, Lepton Flavor Mixing and Leptogenesis in the Minimal Seesaw Model,Int. J. Mod. Phys. E16(2007) 1 [hep-ph/0612033]

  48. [48]

    New aspects of leptogenesis bounds

    S. Blanchet and P. Di Bari,New aspects of leptogenesis bounds,Nucl. Phys. B807(2009) 155 [0807.0743]

  49. [49]

    Granelli, K

    A. Granelli, K. Moffat, Y. F. Perez-Gonzalez, H. Schulz and J. Turner,ULYSSES: Universal LeptogeneSiS Equation Solver,Comput. Phys. Commun.262(2021) 107813 [2007.09150]

  50. [50]

    Granelli, C

    A. Granelli, C. Leslie, Y. F. Perez-Gonzalez, H. Schulz, B. Shuve, J. Turner et al.,ULYSSES, universal LeptogeneSiS equation solver: Version 2,Comput. Phys. Commun.291(2023) 108834 [2301.05722]. 18