REVIEW 4 major objections 3 minor 26 references
Heavy Neutral Lepton Decay Searches using Solar Neutrinos
T0 review · 4 major / 3 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read Solar boron-8 neutrinos could probe heavy neutral lepton mixing far below current bounds.
desk verdict Genuinely new differential width and interesting detector projections, but the central rate formula is internally inconsistent in its |UeN|^2 scaling and no code is provided, so the numbers cannot be trusted as written. 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 load-bearing object is the fully differential decay width $d^2\Gamma/(dl_0\,d\cos\theta)$ for $N \to \nu_e e^+ e^-$, including neutral-current/charged-current interference. The coefficients $X = [(g_V+1)^2 + (g_A+1)^2]$ and $Z = [(g_V+1)^2 - (g_A+1)^2]$ appear in the width; $X$ multiplies the structure already present in the charged-current-only result, while $Z$ multiplies new interference terms proportional to $m_e^2$. This width feeds the event-rate formula that includes the solar HNL flux, the survival factor from the Sun to Earth, and the decay-in-detector probability, and it is complemented by a kinematic calculation showing that the $e^+e^-$ opening angle stays above about 50 degrees for HNL energies from 2 to 16 MeV.
What would settle it
A detector-level simulation of the gamma-ray background in a 400-ton argon or 50-ton xenon detector that finds even a handful of $e^+e^-$-like coincidence events in the design exposure would invalidate the quoted 2.3-event sensitivity, and would directly test the tracking projections.
Extended reading notes
Core claim
The central claim is that the decay $N \to \nu_e e^+ e^-$, fed by solar boron-8 neutrinos, is a viable and potentially background-free channel for searching heavy neutral leptons with masses between the electron-positron threshold and about 16 MeV. The paper's new closed-form differential decay width reproduces the earlier charged-current result and adds the neutral-current and interference terms, organized through the coefficients $X = (g_V+1)^2 + (g_A+1)^2$ and $Z = (g_V+1)^2 - (g_A+1)^2$. With this width, the Borexino 0.122 kt-yr exposure gives a 90% confidence limit that is somewhat weaker than the collaboration's published limit, demonstrating that the interference is non-negligible. The paper also finds that a 100 kt-yr low-energy DUNE could improve the limit by about a factor of 30 using a calorimetric analysis, while track-resolving detectors XLZD and Argo could improve it by 30–80 for 1–5 kt-yr exposures, and a background-free kiloton-scale LiquidO detector could reach more than two orders of magnitude better. The lab-frame opening angle between the $e^+$ and $e^-$ is always at least about 50 degrees across the parameter space, which is what makes the two-track tagging signature resolvable.
Load-bearing premise
The projected multi-order-of-magnitude improvements assume that kiloton-scale detectors can tag electron-positron pairs with high efficiency and zero background using position and timing cuts alone, a claim not yet backed by a detector simulation.
Editorial extensions
If this is right
- A low-energy upgraded DUNE with 100 kton-year exposure could reach a mixing-angle limit about 30 times stronger than Borexino through a calorimetric analysis of the $e^+e^-$ pair plus annihilation gammas.
- Tracking detectors XLZD (50 tons of xenon) and Argo (400 tons of argon) could improve the limit by factors of 30–80 with 1–5 kton-year exposures, using the two tracks and the two back-to-back 0.511 MeV gammas as the signature.
- If LiquidO can tag $e^+e^-$ pairs with zero background in kiloton-scale detectors, a 25 kton-year exposure would push the mixing-angle reach more than two orders of magnitude beyond current bounds.
- Because the $e^+e^-$ opening angle never drops below about 50 degrees, the two-track signature should be resolvable throughout the entire 2–16 MeV HNL mass range.
- The published Borexino limit is somewhat weaker once the NC/CC interference is included, so future solar-neutrino analyses should use the new width rather than the charged-current-only result.
Reading between the lines
- The tracking projections rest on the assumption that kiloton-scale detectors can tag $e^+e^-$ pairs with zero background; if argon-39 beta decays or external gamma pair conversion cannot be fully rejected, the realistic sensitivity likely falls back toward the calorimetric DUNE level, and a detector simulation is the missing test.
- The same $e^+e^-$ final state and tracking signature could be applied to other high-energy neutrino sources, such as diffuse supernova neutrinos, to extend HNL searches beyond the 16 MeV reach of solar boron-8 neutrinos.
- Because the interference term changes both the rate and the spectral shape of the $e^+e^-$ energy distribution, a shape-based fit to existing or future solar-neutrino data could independently validate the new decay width.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper proposes using the solar boron-8 neutrino flux as a source of heavy neutral leptons (HNLs) and searches for their decay N -> nu_e e+ e- in large detectors. It derives a fully differential decay width including charged- and neutral-current interference, reanalyzes the Borexino bound, and projects 90% C.L. sensitivities for DUNE, XLZD, Argo, and LiquidO. The headline quantitative claim is that future detectors can improve the |U_eN|^2 reach by one to two orders of magnitude relative to existing bounds.
Significance. The analytic expression for the NC/CC interference in the differential width is a potentially useful contribution, and the idea of using the solar 8B flux as an HNL source is interesting. If the rates were correct, the paper would provide new, falsifiable search channels for MeV-scale HNLs. However, as written, the event-rate normalization is not reproducible, the treatment of the initial HNL spin is internally inconsistent, and the tracking-based projections rest on an unquantified background-free assumption. These issues block acceptance of the numerical results.
major comments (4)
- [Integrated Decay Rates, Eqs. (4), (5), (9)] Equation (4) has the wrong normalization for the number of decays and is dimensionally inconsistent. In natural units, V T times the integrated flux leaves a negative mass dimension, and the factor (m_N/p_N)[Gamma_e |U_eN|^2/Gamma_tot] does not supply the required factor of an actual decay width. A dimensionally correct expression must contain Gamma_tot through the inverse decay length, or equivalently Gamma_e^{actual}. If Gamma_e and Gamma_tot are the actual partial and total widths, then Gamma_e/Gamma_tot is the branching ratio and the explicit |U_eN|^2 in Eq. (4) is spurious. If they are reduced widths with the mixing angle removed, the text after Eq. (4) and the branching-ratio interpretation must be changed. Equation (9) repeats the same structure, and because no numerical implementation is provided, the limits in Figs. 1 and 2 cannot be checked against the formulas. I do not find the |U_eN|^6 power counting to be the cleanest statement of the problem (once Gamma_tot is also proportional to |U_eN|^2, the ratio Gamma_e/Gamma_tot is mixing-independent), but the normalization defect is real and load-bearing. The authors must correct the rate formula and regenerate all limits and sensitivity curves.
- [SM Decays of HNL and Appendix, Eq. (3)] The treatment of the HNL spin is contradictory. The main text says that the differential width is obtained by summing over final-state spins and averaging over the two possible spin states of N, which would set |s_vec|=0 and remove the |s_vec| cos(theta) terms in Eq. (3). The appendix, however, retains those terms and says that because the calculation is interested in the polarization of the initial state, only the final-state spins are summed. If the polarization terms are physical, the solar HNL production mechanism must provide the polarization vector s_vec as a function of energy; no such distribution is given. If the terms are to vanish, Eq. (3), Eq. (23), and the surrounding discussion must be revised consistently. Since Eq. (9) is built from this differential distribution, the inconsistency directly affects the predicted event spectra.
- [Future Sensitivities and Fig. 2] The tracking-based projections assume that e+e- pairs can be tagged with zero background and high efficiency. The text itself states that a detailed detector study is required for argon-39 and that the LiquidO projection is based on an assumed background-free tagging capability. Because the headline improvement of more than two orders of magnitude is carried by these curves, especially the 25 kton-yr LiquidO line, the paper should either provide a quantitative background model and an efficiency estimate (external gamma pair conversion, Compton pile-up, annihilation-gamma misidentification) or explicitly label those curves as idealized ceilings rather than sensitivities.
- [Existing Constraints] The Borexino reanalysis is not reproducible. The authors say that they digitize the electron recoil data from Fig. 4 of Ref. [1] and perform a binned likelihood fit with that data as background, but they do not specify the binning, the elastic-scattering background normalization and its uncertainty, the treatment of systematic errors, the signal efficiency, or the method used to derive the 90% C.L. contour. This matters because the comparison in Fig. 1 is used to demonstrate that the NC/CC interference is non-negligible.
minor comments (3)
- [Throughout] There are many typographical errors and inconsistent notations, for example 'predicitions', 'sensivity', and |U_eH|^2 where |U_eN|^2 is clearly intended; the manuscript needs a careful proofread.
- [Eq. (3)] The variables in Eq. (3) are not fully defined: Q is introduced as a momentum transfer but appears in dimensionless ratios without an explicit scaling by m_N, and the definitions of g_V and g_A and their relation to the Standard Model couplings should be stated explicitly.
- [Fig. 2] Figure 2 is very dense; the constraint labels and sensitivity curves should be enlarged or placed in a readable legend so that each curve can be identified without close scrutiny.
Circularity Check
No circular reasoning detected; the derivation chain is self-contained against external benchmarks.
full rationale
The paper's central calculation is a first-principles decay-width derivation (Eq. 3 and the appendix) combined with externally supplied inputs: the standard solar boron-8 neutrino flux, the published Borexino data, and independent detector specifications and exposures. The production flux in Eq. (5) and the decay-rate expression in Eq. (9) are not defined in terms of the sensitivity bounds they produce; the bounds are outputs of a likelihood fit, not inputs. No fitted parameter is renamed as a prediction, and no load-bearing uniqueness theorem or ansatz is imported from the authors' prior work. The explicit assumptions, such as background-free e+e- tagging for LiquidO, are stated as assumptions and are not disguised as derived results. The possible over-counting of |UeN|^2 between Eqs. (5) and (9) is an internal normalization and correctness concern, not a circularity, because it does not make any predicted limit true by construction; it is a mechanical error that would make the numerical projections unreliable but does not reduce the derivation to its own inputs.
Assumptions & free parameters
free parameters (2)
- DUNE energy resolution =
7%, energy-independent above 5 MeV
- Tracking detection efficiency =
100% (implicit)
assumptions (4)
- domain assumption The solar boron-8 flux at Earth, without flavor conversion, is the production source for HNLs via the mixing |UeN|^2.
- ad hoc to paper The HNL spin polarization is either zero or averaged in the rate calculation.
- domain assumption Detector backgrounds for the tracking mode are negligible after spatial and temporal cuts.
- standard math The effective Hamiltonian in Eq. (10) with SM couplings gV and gA describes the NC and CC contributions to N to nu e+ e-.
Cite this review
Pith. "Pith review of Heavy Neutral Lepton Decay Searches using Solar Neutrinos." pith.science (2026). https://pith.science/paper/6GY4LTCP
@misc{pith2026250604099,
author = {Pith},
title = {Pith review of: Heavy Neutral Lepton Decay Searches using Solar Neutrinos},
year = {2026},
howpublished = {\url{https://pith.science/paper/6GY4LTCP}},
note = {Machine review of arXiv:2506.04099}
}
abstract
We study the sensitivity to the decay of a heavy neutral lepton into $e^+e^-$-pairs using the solar boron-8 neutrino flux as source. We provide a fully differential cross section for this process including the interference of neutral and charged current amplitudes. We revisit a previous bound from Borexino and make predicitions for the expected sensitivity in future large liquid noble gas detectors, like XLZD, Argo and DUNE, as well as high-resolution scintillator detectors based on the LiquidO technology. We find that more than two orders of magnitude improvement in mixing angle reach is possible relative to existing bounds.
Figures
Reference graph
Works this paper leans on
-
[1]
G. Bellini et al. (Borexino), Phys. Rev. D 88, 072010 (2013), arXiv:1311.5347 [hep-ex]
arXiv 2013
-
[2]
F. Capozzi, S. W. Li, G. Zhu, and J. F. Beacom, Phys. Rev. Lett. 123, 131803 (2019), arXiv:1808.08232 [hep- ph]
arXiv 2019
-
[3]
J. Aalbers et al. , J. Phys. G 50, 013001 (2023), arXiv:2203.02309 [physics.ins-det]
arXiv 2023
-
[4]
A. B. McDonald, Nucl. Phys. B 1003, 116436 (2024)
work page 2024
-
[5]
A. Cabrera et al. (LiquidO), Commun. Phys. 4, 273 (2021), arXiv:1908.02859 [physics.ins-det]
arXiv 2021
-
[6]
D. Gorbunov and M. Shaposhnikov, JHEP 10, 015 (2007), [Erratum: JHEP 11, 101 (2013)], arXiv:0705.1729 [hep-ph]
arXiv 2007
-
[7]
R. E. Shrock, Phys. Rev. D 24, 1275 (1981)
work page 1981
- [8]
Show all 26 references
-
[9]
Hagner, M
C. Hagner, M. Altmann, F. von Feilitzsch, L. Oberauer, Y. Declais, and E. Kajfasz, Phys. Rev. D 52, 1343 (1995)
1995
-
[10]
Oberauer, F
L. Oberauer, F. Von Feilitzsch, and R. L. Mossbauer, Phys. Lett. B 198, 113 (1987)
1987
-
[11]
V. I. Kopeikin, L. A. Mikaelyan, and S. A. Fayans, JETP Lett. 51, 86 (1990)
1990
-
[12]
D. A. Bryman and R. Shrock, Phys. Rev. D 100, 073011 (2019), arXiv:1909.11198 [hep-ph]
2019 arXiv
-
[13]
de Gouvˆ ea and A
A. de Gouvˆ ea and A. Kobach, Phys. Rev. D 93, 033005 (2016), arXiv:1511.00683 [hep-ph]
2016 arXiv
-
[14]
Grossheim et al
A. Grossheim et al. (TWIST), Phys. Rev. D 80, 052012 (2009), arXiv:0908.4270 [hep-ex]
2009 arXiv
-
[15]
Bayes (TWIST), J
R. Bayes (TWIST), J. Phys. Conf. Ser. 408, 012071 (2013)
2013
-
[16]
Drewes et al., JCAP 01, 025 (2017), arXiv:1602.04816 [hep-ph]
M. Drewes et al., JCAP 01, 025 (2017), arXiv:1602.04816 [hep-ph]
2017 arXiv
-
[17]
Boyarsky, M
A. Boyarsky, M. Drewes, T. Lasserre, S. Mertens, and O. Ruchayskiy, Prog. Part. Nucl. Phys. 104, 1 (2019), arXiv:1807.07938 [hep-ph]
2019 arXiv
-
[18]
Sabti, A
N. Sabti, A. Magalich, and A. Filimonova, JCAP 11, 056 (2020), arXiv:2006.07387 [hep-ph]
2020 arXiv
-
[19]
S. A. Baranov et al., Phys. Lett. B302, 336 (1993)
1993
-
[20]
Bernardi et al., Phys
G. Bernardi et al., Phys. Lett. B 203, 332 (1988)
1988
-
[21]
F. P. Calaprice and D. J. Millener, Phys. Rev. C 27, 1175 (1983)
1983
-
[22]
A. I. Derbin, A. V. Chernyi, L. A. Popeko, V. N. Mura- tova, G. A. Shishkina, and S. I. Bakhlanov, JETP Lett. 57, 768 (1993)
1993
-
[23]
Aguilar-Arevalo et al
A. Aguilar-Arevalo et al. (PIENU), Phys. Rev. D 97, 072012 (2018), arXiv:1712.03275 [hep-ex]
2018 arXiv
-
[24]
P. D. Bolton, F. F. Deppisch, and P. S. Bhupal Dev, JHEP 03, 170 (2020), arXiv:1912.03058 [hep-ph]
2020 arXiv
-
[25]
Apilluelo et al
J. Apilluelo et al. (LiquidO), Nucl. Instrum. Meth. A 1071, 170075 (2025), arXiv:2406.13054 [physics.ins-det]
2025 arXiv
-
[26]
Apilluelo et al
J. Apilluelo et al. (LiquidO), (2025), arXiv:2503.02541 [physics.ins-det]. 6 Cross Section The charged current result was previously computed by Shrock [7] and we confirm this earlier result. In the following we will present the result including the neutral current and resulti...
2025 arXiv
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.