Pith. sign in

REVIEW 4 major objections 5 minor 72 references

In right-handed neutrino portal freeze-in, ignoring dark-sector internal dynamics underestimates the final dark matter abundance by 30% or more than 95%.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-03 17:01 UTC pith:HIKUQZ3Y

load-bearing objection The 30%/95% freeze-in error claim is plausible but unverifiable as written; the stress-test's number-conservation objection misses the N1N1→χχ channel, so the worry is not fatal, but the paper needs code or explicit matrix elements to be trustworthy. the 4 major comments →

arxiv 2512.10762 v2 pith:HIKUQZ3Y submitted 2025-12-11 hep-ph

Electroweak right-handed neutrino portal dark matter

classification hep-ph
keywords dark matterright-handed neutrino portalfreeze-infreeze-outType-I seesawrelic abundanceBoltzmann equationsheavy neutral leptons
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper establishes that in freeze-in dark matter production through an electroweak right-handed neutrino portal, the internal interactions among dark-sector particles materially change the predicted relic density. Solving the full coupled Boltzmann equations — including a separate dark-sector temperature — gives a final dark matter abundance that is roughly 30% higher than the common shortcut of treating each dark species as an independent freeze-in component and then adding late decays, and more than 95% higher when the three-point decay N↔χφ is forbidden. The authors construct three concrete Type-I seesaw realizations (regular, structure cancellation, and split seesaw) whose parameters are fixed by neutrino oscillation data and experimental constraints, and they show that a fermion–scalar dark sector can reproduce the observed relic density through either secluded freeze-out or freeze-in. The upshot is that relic-density calculations in portal freeze-in models should not decouple the dark species; the coupled system is required.

Core claim

The paper claims that in the right-handed neutrino portal freeze-in framework, internal dark-sector interactions substantially alter the final DM relic abundance. When the dark sector reaches internal thermal equilibrium, one must introduce a hidden-sector temperature T_h distinct from the visible-sector temperature T_v and evolve it together with the comoving number densities. The authors compare their full multi-temperature treatment with the approximate two-step procedure standard in parts of the literature — first compute freeze-in for each species separately, then add late decays — and find that the approximation underestimates the final abundance by about 30% when the three-point proce

What carries the argument

The central machinery is the set of coupled Boltzmann equations for the comoving number densities Y_χ, Y_φ, and Y_Ni, supplemented by the temperature ratio η = T_v/T_h when the dark sector self-thermalizes. The portal is the renormalizable Yukawa y_x χφN, which after seesaw mixing becomes a sum over neutrino mass eigenstates. The three seesaw benchmarks (regular, structure cancellation, and split) provide numerically fixed portal couplings derived from neutrino data, and the full collision terms include all relevant two-to-two, one-to-two, and internal dark-sector processes with real-intermediate-state subtraction.

Load-bearing premise

The dark sector is assumed to couple to the Standard Model only through the right-handed neutrino portal Yukawa y_x χφN; if any additional coupling exists — Higgs portal, kinetic mixing, a light dark photon, or extra Yukawas — the freeze-in/freeze-out yields and the quoted 30%/95% error figures would change.

What would settle it

For Model d (Benchmark-SS1) and Model e (Benchmark-SS2), compute the relic density with an independent Boltzmann solver that includes all internal dark-sector interactions and compare to the paper's full values (Ωχh²=0.12 in both cases). If the independent-species treatment with late decays also gives 0.12 within a few percent for Model d, or if removing the internal interactions changes the abundance by less than a few percent, the paper's 30%/95% discrepancy claim would be falsified.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Freeze-in relic densities in right-handed neutrino portal models must be computed from the coupled system; the independent-species approximation is not reliable when the dark sector self-thermalizes.
  • The paper's benchmarks show that both freeze-out and freeze-in can yield Ωχh²≈0.12 while satisfying current neutrino and collider constraints, giving concrete targets for LHC and future heavy-neutral-lepton searches.
  • In Case-SC, the O(10^-2) portal couplings make the scenario particularly accessible to collider searches and neutrino experiments, so a positive signal would tie dark matter directly to seesaw neutrino mass generation.
  • When the three-point channel N↔χφ is absent, the error from the approximate treatment exceeds 95%, meaning that any study using the shortcut in such models should be revisited.
  • The multi-temperature treatment used here provides a template for other hidden-sector models where the dark sector has its own temperature and internal equilibrium.

Where Pith is reading between the lines

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

  • The same failure of the independent-species approximation likely occurs in other freeze-in portal constructions (Higgs portal, dark photon portal) with self-interacting dark sectors; the magnitude of the bias should scale with the strength of internal number-changing versus number-conserving interactions.
  • A dark sector that reaches internal equilibrium during freeze-in will have a non-thermal or multi-temperature momentum distribution for the dark matter, which could leave imprints in structure formation (e.g., Lyman-α forest) or in the effective number of relativistic species; these are testable with upcoming surveys.
  • The >95% error when N↔χφ is forbidden suggests that the underestimate is driven by the absence of a channel that directly transfers abundance between dark species; mapping this error as a function of the mass splitting m_φ−m_χ and the coupling y_x would give a practical diagnostic for when the shortcut can be used.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The paper studies dark matter in a right-handed neutrino portal framework, with a dark fermion χ and complex scalar φ coupled to the SM only through Yukawa interactions with right-handed Majorana neutrinos. Three seesaw realizations (regular, structure-cancellation, split seesaw) are constructed, and Particle Swarm Optimization is used to fit seesaw parameters to neutrino data. The authors solve coupled Boltzmann equations for the dark and neutrino sectors for both freeze-out and freeze-in production, and claim that a common treatment of freeze-in—treating each dark species as an independent freeze-in component and adding late decays afterward—underestimates the final relic abundance by ~30% when N↔χφ is open and by >95% when it is closed. The paper discusses experimental constraints, benchmark models, and the multi-temperature hidden sector formalism.

Significance. If the 30%/95% claim is correct, it provides a useful cautionary result for freeze-in calculations in neutrino portal models, where internal dark sector dynamics can significantly modify the relic density. The systematic construction of viable seesaw benchmarks via PSO is a contribution, and the multi-temperature formalism is published previously. However, the central quantitative claim is not fully reproducible from the manuscript: the approximate procedure is not specified algorithmically, the matrix elements and cross-sections are not given in closed form, and no code is provided. The paper also helpfully identifies that dark-number-changing processes involving the portal neutrino must be included, which resolves a potential paradox in Model d.

major comments (4)
  1. [§4.3.2, Models d and e] The 'common approximate procedure' is described only in prose and never written as an algorithm or equations. The central 30%/95% claim is defined as the difference between this procedure and the full solution, so the procedure must be specified precisely: the list of production processes for each species, how LH→χϕ is apportioned between χ and ϕ, when decays are applied, and which reverse/annihilation processes (e.g., χχ→N1N1) are omitted. Without this, the comparison is not well-defined and cannot be checked.
  2. [§4.2, Eqs. (41)-(52)] The Boltzmann equations are written in terms of thermally averaged cross-sections ⟨σv⟩, decay widths ⟨Γ⟩, and J-functionals whose explicit forms are not given; only references to Feynman diagrams are provided. Since the 30%/95% result is a numerical output of a private solver, the authors should provide the code, the full matrix elements, or at least a validation against a simpler limit. As it stands, the readers cannot independently reproduce the central numerical claim.
  3. [§4.3.2, Model d] In Model d, M1 < mχ + mϕ, so N1→χϕ is closed, but the text states that χχ→N1N1 and ϕϕ→N1N1 occur and N1 subsequently decays to SM particles. These processes do not conserve the total dark number, so the final χ yield is not simply the total dark number produced by LH→χϕ. This should be stated explicitly when comparing with the approximate procedure, and the contribution of N1-mediated channels to the >95% error should be quantified. Otherwise readers may incorrectly infer a number-conservation paradox.
  4. [Table 2 and Fig. 8] The 30%/95% error is demonstrated only for two tuned benchmark points where yx is adjusted so the full calculation gives Ωχh²=0.12. The abstract and conclusion present these numbers as general statements. Please show the error as a function of yx over the viable parameter range, and ideally for more freeze-in benchmarks, to establish that the effect is not an artifact of the selected points. Also clarify that the 'reproduction' of Ω=0.12 is a fit to yx, not a parameter-free prediction.
minor comments (5)
  1. [Throughout] Typos and grammar: 'Majorara' (§2.3.1), 'diaganolization' (App. B), 'šCsterile' (Fig. 3 caption), 'the eventually decrease' (§4.3.2), missing space in 'electroweakright-handedneutrinoportaldarkmatterprovidesarobustframework' (Conclusion).
  2. [Eq. (35)] The definition of yx'_i is ambiguous, as the index contraction with the seesaw matrix U is not written explicitly. Please define yx'_j = Σ_i yx_i U_{i+3,j} (or the equivalent) to clarify which mass eigenstates couple to the dark sector.
  3. [References] Refs. [37] and [63] have placeholder-looking DOIs ('10.1103/vshn-gbdp' and '10.1103/df3g-32t9') and Ref. [38] is 'to appear'. These should be corrected or updated before publication.
  4. [Appendix H, Eq. (128)] The equation contains a garbled expression 'ρ dρh/dTh' that is likely a typo for 'dρh/dTh'. The surrounding derivation should be proofread.
  5. [§4.3.2, Model d paragraph] The sentence 'since the right-handed neutrino N1 cannot decay into dark particles, the freeze-in contribution from N1→χϕ is absent. As a result, the approximate treatment deviates strongly...' is logically unclear. The deviation arises from neglecting internal dark-sector interactions and the treatment of N1 dynamics, not merely from the absence of N1→χϕ. Rephrase for clarity.

Circularity Check

0 steps flagged

No significant circularity: the 30%/95% claim is a controlled Boltzmann-scheme comparison at fixed couplings; relic-density matching is an acknowledged benchmark fit, not a prediction.

full rationale

The central quantitative claim is a comparison between two solution schemes at fixed benchmark parameters: the full coupled multi-temperature Boltzmann system (Eqs. 48-51, Appendix H) and an approximate independent-component plus late-decay procedure. The approximate treatment is an external approximation rather than a rearrangement of the full equations, so the 30%/95% discrepancy is a computed numerical result, not a relation true by construction. The benchmark models do contain a fitted element: yx is chosen so that the full calculation returns Ωχh^2 ≈ 0.12, as acknowledged in the freeze-out table caption ('Parameters are chosen to satisfy ... the observed abundance of the dark matter') and visible in the Ωχh^2 column of Table 2. This makes 'reproducing the observed relic abundance' a benchmark-matching exercise rather than an independent prediction, but it does not infect the 30%/95% comparison, because both numbers are evaluated with the same chosen couplings. The hidden-temperature formalism is drawn substantially from the authors' prior work (Refs [33-37]), but Appendix H re-derives the relevant equations, and the method does not assume or contain the 30%/95% result. No equation in the paper is set equal to its target by construction; no uniqueness theorem is imported; and no fitted parameter is disguised as a prediction. The approximate procedure is described in prose rather than as an explicit algorithm, which is a reproducibility concern but not evidence of circularity. No significant circularity found.

Axiom & Free-Parameter Ledger

4 free parameters · 5 axioms · 2 invented entities

The calculation rests on standard seesaw physics plus a self-cited multi-temperature formalism and the pure-portal assumption. The DM abundance is fit by tuning yx and masses; the seesaw parameters are fit to neutrino data. The paper's genuinely new quantitative claim (30%/95% error) is computed but not independently verifiable from the text alone.

free parameters (4)
  • Dark Yukawa coupling yx = e.g., 0.433 (FO Model a), 2.2e-6 (FI Model a), 0.29 (SS1), 1.7e-12 (SS3)
    Adjusted in Tables 1 and 2 so that the calculated Ωχh² equals the observed 0.12 (or the combined 0.12 in Model f). This is a fit to the relic-density target, not a prediction.
  • Dark sector masses mχ, mφ = e.g., (210,350), (240,360), (440,450) GeV
    Benchmark mass points chosen by hand; they set kinematics for decays and annihilations.
  • Seesaw input parameters (yν_ij, Mi) for Benchmarks RS/SC/SS = Appendix I; e.g., Benchmark-RS: yν entries 1e-7–4.6e-7, Mi = 200,220,260 GeV
    Fitted by Particle Swarm Optimization to reproduce |U_PMNS| (Eq. 30), mass-squared differences (Eqs. 32–33), mass-sum bound (Eq. 31), and for Case-SC non-unitarity bounds (Appendix D). Fitted to independent data.
  • Case-SC fine-tuning parameter ε = 6.34546e-10
    Tiny parameter in the structure-cancellation matrix Eq. (19), chosen so that an otherwise rank-1 Yukawa structure generates the observed tiny neutrino masses. Ad hoc small parameter.
axioms (5)
  • standard math Type-I seesaw mixing: light neutrino mass matrix mν ≈ -MD MN^{-1} MD^T via Schur-complement diagonalization (Eq. 24)
    Standard seesaw mechanism, used as background throughout Section 2.
  • domain assumption Multi-temperature hidden-sector Boltzmann formalism (Eqs. 48–52, Appendix H) from Ref [33] is valid
    The coupled evolution with hidden temperature Th is taken from the authors' own previous work (Ref [33]) and applied without re-derivation. Published, peer-reviewed, but not independently reproduced here.
  • ad hoc to paper Dark sector has no additional couplings to the SM besides the right-handed neutrino portal
    Section 3 assumes the dark U(1) gauge boson is heavy/decoupled and the evolutions of χ, φ are purely portal-driven. If other portals exist, the calculated relic densities and the 30%/95% error figures change.
  • domain assumption Classical Maxwell-Boltzmann statistics and standard thermally-averaged cross sections are used in all collision terms
    Standard in DM Boltzmann solvers; quantum statistics and statistical factors are neglected, which can shift freeze-in yields at the O(1–10%) level.
  • domain assumption The 'common approximate procedure' (independent freeze-in components + late decays) is representative of the literature
    The paper critiques this approximation but does not cite a specific reference that uses it, making the exact target of the 30%/95% comparison difficult to verify.
invented entities (2)
  • Dark fermion χ no independent evidence
    purpose: Dark matter candidate; stable under a dark Z2/U(1); couples to right-handed neutrinos via yx χφN.
    Mass and coupling yx are benchmark parameters fitted to the relic density; no independent collider or astrophysical signature is computed for χ itself.
  • Dark complex scalar φ no independent evidence
    purpose: Partner for χ in the portal Yukawa; heavier than χ; decays to χ+ν before BBN or contributes to DM if long-lived.
    Its mass, coupling, and lifetime are chosen to satisfy benchmark conditions; no independent evidence is provided.

pith-pipeline@v1.3.0-alltime-deepseek · 37472 in / 14235 out tokens · 145085 ms · 2026-08-03T17:01:33.361363+00:00 · methodology

0 comments
read the original abstract

We study dark matter coupled to the Standard Model via electroweak scale right-handed neutrinos in a Type-I seesaw framework. We consider a minimal dark sector containing a fermion $\chi$ and a complex scalar $\phi$ whose only connection to the Standard Model is through renormalizable Yukawa interactions with right-handed Majorana neutrinos, thus realizing a neutrino portal after seesaw mixing. We discuss three representative realizations of electroweak right-handed neutrinos arising from the Type-I seesaw mechanism, spanning small, tiny, and ultraweak couplings to the Standard Model sector, so that the dark particles can either undergo secluded freeze-out or be produced via freeze-in. Instead of merely estimating the order of magnitude of the seesaw couplings, we use the Particle Swarm Optimization algorithm to obtain viable seesaw parameter sets consistent with neutrino data and other constraints, and then compute the coupled evolution of the dark particles and right-handed neutrinos, reproducing the observed dark matter relic abundance in representative benchmark scenarios. For freeze-out, dark matter depletion is controlled by coupled dark sector dynamics, requiring a full Boltzmann treatment for a reliable relic abundance. For freeze-in, internal dark interactions also alter the relic density: treating hidden particles as independent components with late decays added afterward can misestimate the abundance by $30\%$ or even $95\%$, depending on the interaction structure. Electroweak right-handed neutrino portal dark matter thus provides a robust and predictive framework that tightly connects neutrino physics, heavy neutral lepton phenomenology, and the cosmological dark matter relic density, offering a well-motivated benchmark for complementary collider, neutrino, and cosmological probes at the high energy frontier.

Figures

Figures reproduced from arXiv: 2512.10762 by Ao Li, Wan-Zhe Feng, Zi-Hui Zhang, Zong-Huan Ye.

Figure 1
Figure 1. Figure 1: Exclusion limits at the 95% confidence level on the active-sterile mixing parameter [PITH_FULL_IMAGE:figures/full_fig_p015_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Exclusion limits at the 95% confidence level on the active sterile mixing parameter [PITH_FULL_IMAGE:figures/full_fig_p016_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Exclusion limits at the 95% confidence level on the activešCsterile mixing parameter [PITH_FULL_IMAGE:figures/full_fig_p016_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: A schematic illustration of the dark sector evolution via freeze-out. At high [PITH_FULL_IMAGE:figures/full_fig_p021_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: An exhibition of the dark sector freeze-out evolution for Models [PITH_FULL_IMAGE:figures/full_fig_p023_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: For Case-RS and SC, the right-handed neutrinos reach thermal equilibrium during [PITH_FULL_IMAGE:figures/full_fig_p024_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: For Case-SS, the portal right-handed neutrino [PITH_FULL_IMAGE:figures/full_fig_p025_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: An exhibition of the freeze-in evolution for Model [PITH_FULL_IMAGE:figures/full_fig_p028_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Relevant Feynman diagrams for 2 → 2 scattering processes are shown. The first row depicts dark particle χ, ϕ annihilation into light or heavy neutrinos after seesaw mixing; the inverse processes lead to freeze-in production of the dark sector particles. The second row shows internal dark sector interactions mediated by light or heavy neutrinos. The third row illustrates processes arising from χϕ → LH after… view at source ↗
Figure 10
Figure 10. Figure 10: Feynman diagrams of the heavy neutrino N1,2,3 decaying to dark sector particles χ, ϕ, and to SM particles. 42 [PITH_FULL_IMAGE:figures/full_fig_p042_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Feynman diagram of ϕ decay to χ and light (always present) or heavy (if mϕ > mχ + mNi ) neutrinos. where g v eff (g h eff) is the visible (hidden) effective degrees of freedom. Further one can deduct the derivative of ρv, ρh w.r.t. Th ρ dρh dTh = [PITH_FULL_IMAGE:figures/full_fig_p043_11.png] view at source ↗

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

72 extracted references · 2 canonical work pages

  1. [1]

    Khachatryanet al.[CMS], Phys

    V. Khachatryanet al.[CMS], Phys. Lett. B748, 144-166 (2015) doi:10.1016/j.physletb.2015.06.070 [arXiv:1501.05566 [hep-ex]]

  2. [2]

    A. M. Sirunyanet al.[CMS], Phys. Rev. Lett.120, no.22, 221801 (2018) doi:10.1103/PhysRevLett.120.221801 [arXiv:1802.02965 [hep-ex]]

  3. [3]

    A. M. Sirunyanet al.[CMS], JHEP01, 122 (2019) doi:10.1007/JHEP01(2019)122 [arXiv:1806.10905 [hep-ex]]

  4. [4]

    Tumasyanet al.[CMS], Phys

    A. Tumasyanet al.[CMS], Phys. Rev. Lett.131, no.1, 011803 (2023) doi:10.1103/PhysRevLett.131.011803 [arXiv:2206.08956 [hep-ex]]

  5. [5]

    Aadet al.[ATLAS], Eur

    G. Aadet al.[ATLAS], Eur. Phys. J. C83, no.9, 824 (2023) doi:10.1140/epjc/s10052- 023-11915-y [arXiv:2305.14931 [hep-ex]]

  6. [6]

    Hayrapetyanet al.[CMS], JHEP06, 123 (2024) doi:10.1007/JHEP06(2024)123 [arXiv:2403.00100 [hep-ex]]

    A. Hayrapetyanet al.[CMS], JHEP06, 123 (2024) doi:10.1007/JHEP06(2024)123 [arXiv:2403.00100 [hep-ex]]

  7. [7]

    Aadet al.[ATLAS], Phys

    G. Aadet al.[ATLAS], Phys. Lett. B856, 138865 (2024) doi:10.1016/j.physletb.2024.138865 [arXiv:2403.15016 [hep-ex]]

  8. [8]

    Kozynets, P

    T. Kozynets, P. Eller, A. Zander, M. Ettengruber and D. J. Koskinen, JHEP05, 130 (2025) doi:10.1007/JHEP05(2025)130 [arXiv:2407.20388 [hep-ph]]

  9. [9]

    A. D. Dolgov, S. H. Hansen, G. Raffelt and D. V. Semikoz, Nucl. Phys. B590, 562-574 (2000) doi:10.1016/S0550-3213(00)00566-6 [arXiv:hep-ph/0008138 [hep-ph]]

  10. [10]

    I. A. Zelko, T. Treu, K. N. Abazajian, D. Gilman, A. J. Benson, S. Birrer, A. M. Nierenberg and A. Kusenko, Phys. Rev. Lett.129, no.19, 191301 (2022) doi:10.1103/PhysRevLett.129.191301 [arXiv:2205.09777 [hep-ph]]

  11. [11]

    B. M. Roach, S. Rossland, K. C. Y. Ng, K. Perez, J. F. Beacom, B. W. Grefenstette, S. Horiuchi, R. Krivonos and D. R. Wik, Phys. Rev. D107, no.2, 023009 (2023) doi:10.1103/PhysRevD.107.023009 [arXiv:2207.04572 [astro-ph.HE]]

  12. [12]

    G. H. Du, T. N. Li, P. J. Wu, L. Feng, S. H. Zhou, J. F. Zhang and X. Zhang, [arXiv:2501.10785 [astro-ph.CO]]

  13. [13]

    Cortina Gilet al.[NA62], Phys

    E. Cortina Gilet al.[NA62], Phys. Lett. B807, 135599 (2020) doi:10.1016/j.physletb.2020.135599 [arXiv:2005.09575 [hep-ex]]

  14. [14]

    Cortina Gilet al.[NA62], Phys

    E. Cortina Gilet al.[NA62], Phys. Lett. B816, 136259 (2021) doi:10.1016/j.physletb.2021.136259 [arXiv:2101.12304 [hep-ex]]

  15. [15]

    Bloch-Devauxet al.[NA62], [arXiv:2507.07345 [hep-ex]]

    B. Bloch-Devauxet al.[NA62], [arXiv:2507.07345 [hep-ex]]

  16. [16]

    Köllenberger [KATRIN], PoSDISCRETE2022, 078 (2024) doi:10.22323/1.431.0078 47

    L. Köllenberger [KATRIN], PoSDISCRETE2022, 078 (2024) doi:10.22323/1.431.0078 47

  17. [17]

    Acharyaet al.[KATRIN], Nature648, no.8092, 70-75 (2025) doi:10.1038/s41586- 025-09739-9 [arXiv:2503.18667 [hep-ex]]

    H. Acharyaet al.[KATRIN], Nature648, no.8092, 70-75 (2025) doi:10.1038/s41586- 025-09739-9 [arXiv:2503.18667 [hep-ex]]

  18. [18]

    S. P. Li and X. J. Xu, JCAP06, 047 (2023) doi:10.1088/1475-7516/2023/06/047 [arXiv:2212.09109 [hep-ph]]

  19. [19]

    Escudero, N

    M. Escudero, N. Rius and V. Sanz, JHEP02, 045 (2017) doi:10.1007/JHEP02(2017)045 [arXiv:1606.01258 [hep-ph]]

  20. [20]

    Escudero, N

    M. Escudero, N. Rius and V. Sanz, Eur. Phys. J. C77, no.6, 397 (2017) doi:10.1140/epjc/s10052-017-4963-x [arXiv:1607.02373 [hep-ph]]

  21. [21]

    Y. L. Tang and S. h. Zhu, JHEP01, 025 (2017) doi:10.1007/JHEP01(2017)025 [arXiv:1609.07841 [hep-ph]]

  22. [22]

    Bandyopadhyay, E

    P. Bandyopadhyay, E. J. Chun, R. Mandal and F. S. Queiroz, Phys. Lett. B788, 530-534 (2019) doi:10.1016/j.physletb.2018.12.003 [arXiv:1807.05122 [hep-ph]]

  23. [23]

    Borah, M

    D. Borah, M. Dutta, S. Mahapatra and N. Sahu, Phys. Rev. D105, no.1, 015004 (2022) doi:10.1103/PhysRevD.105.015004 [arXiv:2110.00021 [hep-ph]]

  24. [24]

    Coito, C

    L. Coito, C. Faubel, J. Herrero-García, A. Santamaria and A. Titov, JHEP08, 085 (2022) doi:10.1007/JHEP08(2022)085 [arXiv:2203.01946 [hep-ph]]

  25. [25]

    Becker, Eur

    M. Becker, Eur. Phys. J. C79, no.7, 611 (2019) doi:10.1140/epjc/s10052-019-7095-7 [arXiv:1806.08579 [hep-ph]]

  26. [26]

    Bian and Y

    L. Bian and Y. L. Tang, JHEP12, 006 (2018) doi:10.1007/JHEP12(2018)006 [arXiv:1810.03172 [hep-ph]]

  27. [27]

    Cosme, M

    C. Cosme, M. Dutra, T. Ma, Y. Wu and L. Yang, JHEP03, 026 (2021) doi:10.1007/JHEP03(2021)026 [arXiv:2003.01723 [hep-ph]]

  28. [28]

    Y. Du, F. Huang, H. L. Li and J. H. Yu, JHEP12, 207 (2020) doi:10.1007/JHEP12(2020)207 [arXiv:2005.01717 [hep-ph]]

  29. [29]

    Chianese, B

    M. Chianese, B. Fu and S. F. King, JCAP01, 034 (2021) doi:10.1088/1475- 7516/2021/01/034 [arXiv:2009.01847 [hep-ph]]

  30. [30]

    A. Liu, F. L. Shao, Z. L. Han, Y. Jin and H. Li, Eur. Phys. J. C83, no.5, 423 (2023) doi:10.1140/epjc/s10052-023-11609-5 [arXiv:2205.11846 [hep-ph]]

  31. [31]

    Barman, P

    B. Barman, P. S. Bhupal Dev and A. Ghoshal, Phys. Rev. D108, no.3, 035037 (2023) doi:10.1103/PhysRevD.108.035037 [arXiv:2210.07739 [hep-ph]]

  32. [32]

    A. Liu, F. L. Shao, Z. L. Han, Y. Jin and H. Li, Phys. Rev. D109, no.5, 055027 (2024) doi:10.1103/PhysRevD.109.055027 [arXiv:2308.12588 [hep-ph]]

  33. [33]

    Aboubrahim, W

    A. Aboubrahim, W. Z. Feng, P. Nath and Z. Y. Wang, Phys. Rev. D103, no.7, 075014 (2021) doi:10.1103/PhysRevD.103.075014 [arXiv:2008.00529 [hep-ph]]. 48

  34. [34]

    Aboubrahim, W

    A. Aboubrahim, W. Z. Feng, P. Nath and Z. Y. Wang, JHEP06, 086 (2021) doi:10.1007/JHEP06(2021)086 [arXiv:2103.15769 [hep-ph]]

  35. [35]

    W. Z. Feng, Z. H. Zhang and K. Y. Zhang, JCAP05, 112 (2024) doi:10.1088/1475- 7516/2024/05/112 [arXiv:2312.03837 [hep-ph]]

  36. [36]

    W. Z. Feng, J. Li and P. Nath, Phys. Rev. D110, no.1, 015020 (2024) doi:10.1103/PhysRevD.110.015020 [arXiv:2403.09558 [hep-ph]]

  37. [37]

    W. Z. Feng and Z. H. Zhang, Phys. Rev. D112, no.3, 035004 (2025) doi:10.1103/vshn- gbdp [arXiv:2405.19431 [hep-ph]]

  38. [38]

    W. Z. Feng, Y. N. Wang and Y. Zhang, to appear

  39. [39]

    L. E. Ibanez and A. M. Uranga, JHEP03, 052 (2007) doi:10.1088/1126- 6708/2007/03/052 [arXiv:hep-th/0609213 [hep-th]]

  40. [40]

    Cvetic, R

    M. Cvetic, R. Richter and T. Weigand, Phys. Rev. D76, 086002 (2007) doi:10.1103/PhysRevD.76.086002 [arXiv:hep-th/0703028 [hep-th]]

  41. [41]

    Antusch, L

    S. Antusch, L. E. Ibanez and T. Macri, JHEP09, 087 (2007) doi:10.1088/1126- 6708/2007/09/087 [arXiv:0706.2132 [hep-ph]]

  42. [42]

    Cvetic and P

    M. Cvetic and P. Langacker, Phys. Rev. D78, 066012 (2008) doi:10.1103/PhysRevD.78.066012 [arXiv:0803.2876 [hep-th]]

  43. [43]

    L. E. Ibanez and R. Richter, JHEP03, 090 (2009) doi:10.1088/1126-6708/2009/03/090 [arXiv:0811.1583 [hep-th]]

  44. [44]

    Cvetic, J

    M. Cvetic, J. Halverson, P. Langacker and R. Richter, JHEP10, 094 (2010) doi:10.1007/JHEP10(2010)094 [arXiv:1001.3148 [hep-th]]

  45. [45]

    W. l. Guo and Z. z. Xing, Phys. Lett. B583, 163-172 (2004) doi:10.1016/j.physletb.2003.12.043 [arXiv:hep-ph/0310326 [hep-ph]]

  46. [46]

    J. w. Mei and Z. z. Xing, Phys. Rev. D69, 073003 (2004) doi:10.1103/PhysRevD.69.073003 [arXiv:hep-ph/0312167 [hep-ph]]

  47. [47]

    W. l. Guo, Z. z. Xing and S. Zhou, Int. J. Mod. Phys. E16, 1-50 (2007) doi:10.1142/S0218301307004898 [arXiv:hep-ph/0612033 [hep-ph]]

  48. [48]

    N. Nath, Z. z. Xing and J. Zhang, Eur. Phys. J. C78, no.4, 289 (2018) doi:10.1140/epjc/s10052-018-5751-y [arXiv:1801.09931 [hep-ph]]

  49. [49]

    Chianese and S

    M. Chianese and S. F. King, JCAP09, 027 (2018) doi:10.1088/1475-7516/2018/09/027 [arXiv:1806.10606 [hep-ph]]

  50. [50]

    Chianese, B

    M. Chianese, B. Fu and S. F. King, JCAP03, 030 (2020) doi:10.1088/1475- 7516/2020/03/030 [arXiv:1910.12916 [hep-ph]]. 49

  51. [51]

    Z. z. Xing and Z. h. Zhao, Rept. Prog. Phys.84, no.6, 066201 (2021) doi:10.1088/1361- 6633/abf086 [arXiv:2008.12090 [hep-ph]]

  52. [52]

    Ingelman and J

    G. Ingelman and J. Rathsman, Z. Phys. C60, 243-254 (1993) doi:10.1007/BF01474620

  53. [53]

    Kersten and A

    J. Kersten and A. Y. Smirnov, Phys. Rev. D76, 073005 (2007) doi:10.1103/PhysRevD.76.073005 [arXiv:0705.3221 [hep-ph]]

  54. [54]

    Barry, W

    J. Barry, W. Rodejohann and H. Zhang, JHEP07, 091 (2011) doi:10.1007/JHEP07(2011)091 [arXiv:1105.3911 [hep-ph]]

  55. [55]

    Zhang, Phys

    H. Zhang, Phys. Lett. B714, 262-266 (2012) doi:10.1016/j.physletb.2012.06.074 [arXiv:1110.6838 [hep-ph]]

  56. [56]

    K. Y. Jung and K. Siyeon, J. Korean Phys. Soc.81, no.12, 1211-1224 (2022) doi:10.1007/s40042-022-00674-w [arXiv:2205.13860 [hep-ph]]

  57. [57]

    Rauch, [arXiv:1610.08420 [hep-ph]]

    M. Rauch, [arXiv:1610.08420 [hep-ph]]

  58. [58]

    Achardet al.[L3], Phys

    P. Achardet al.[L3], Phys. Lett. B517, 67-74 (2001) doi:10.1016/S0370-2693(01)00993- 5 [arXiv:hep-ex/0107014 [hep-ex]]

  59. [59]

    de Blas, EPJ Web Conf.60, 19008 (2013) doi:10.1051/epjconf/20136019008 [arXiv:1307.6173 [hep-ph]]

    J. de Blas, EPJ Web Conf.60, 19008 (2013) doi:10.1051/epjconf/20136019008 [arXiv:1307.6173 [hep-ph]]

  60. [60]

    Cheung, Y

    K. Cheung, Y. L. Chung, H. Ishida and C. T. Lu, Phys. Rev. D102, no.7, 075038 (2020) doi:10.1103/PhysRevD.102.075038 [arXiv:2004.11537 [hep-ph]]

  61. [61]

    A. Das, S. Jana, S. Mandal and S. Nandi, Phys. Rev. D99, no.5, 055030 (2019) doi:10.1103/PhysRevD.99.055030 [arXiv:1811.04291 [hep-ph]]

  62. [62]

    H. Yang, B. Long and C. F. Qiao, Nucl. Phys. B1004, 116576 (2024) doi:10.1016/j.nuclphysb.2024.116576 [arXiv:2309.16233 [hep-ph]]

  63. [63]

    A. Das, J. Li, S. Mandal, T. Nomura and R. Zhang, Phys. Rev. D112, no.3, 035008 (2025) doi:10.1103/df3g-32t9 [arXiv:2410.21956 [hep-ph]]

  64. [64]

    Blennow, P

    M. Blennow, P. Coloma, E. Fernandez-Martinez, J. Hernandez-Garcia and J. Lopez- Pavon, JHEP04, 153(2017)doi:10.1007/JHEP04(2017)153[arXiv:1609.08637[hep-ph]]

  65. [65]

    Astieret al.[NOMAD], Phys

    P. Astieret al.[NOMAD], Phys. Lett. B570, 19-31 (2003) doi:10.1016/j.physletb.2003.07.029 [arXiv:hep-ex/0306037 [hep-ex]]

  66. [66]

    Avvakumovet al.[NuTeV], Phys

    S. Avvakumovet al.[NuTeV], Phys. Rev. Lett.89, 011804 (2002) doi:10.1103/PhysRevLett.89.011804 [arXiv:hep-ex/0203018 [hep-ex]]

  67. [67]

    D. V. Forero, C. Giunti, C. A. Ternes and M. Tortola, Phys. Rev. D104, no.7, 075030 (2021) doi:10.1103/PhysRevD.104.075030 [arXiv:2103.01998 [hep-ph]]

  68. [68]

    Abeet al.[T2K], Phys

    K. Abeet al.[T2K], Phys. Rev. D103, no.11, 112008 (2021) doi:10.1103/PhysRevD.103.112008 [arXiv:2101.03779 [hep-ex]]. 50

  69. [69]

    M. A. Aceroet al.[NOvA], Phys. Rev. D98, 032012 (2018) doi:10.1103/PhysRevD.98.032012 [arXiv:1806.00096 [hep-ex]]

  70. [70]

    Adamsonet al.[MINOS+], Phys

    P. Adamsonet al.[MINOS+], Phys. Rev. Lett.122, no.9, 091803 (2019) doi:10.1103/PhysRevLett.122.091803 [arXiv:1710.06488 [hep-ex]]

  71. [71]

    E. W. Kolb and S. Wolfram, Nucl. Phys. B172, 224 (1980) [erratum: Nucl. Phys. B 195, 542 (1982)] doi:10.1016/0550-3213(82)90012-8

  72. [72]

    G. F. Giudice, A. Notari, M. Raidal, A. Riotto and A. Strumia, Nucl. Phys. B685, 89-149 (2004) doi:10.1016/j.nuclphysb.2004.02.019 [arXiv:hep-ph/0310123 [hep-ph]]. 51