{"id":"7476911d-c5ff-40dd-a8c3-c1289bd2abd6","arxiv_id":"2608.10123","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A Dirac heavy neutral lepton with a particle-antiparticle asymmetry can evade the hadronic BBN bound by injecting charged pions that restore the standard neutron abundance, opening sub-GeV parameter space accessible to SHiP.","lead":"The paper shows that a special type of heavy neutrino, if it carries a mismatch between particles and antiparticles, can dodge the usual Big Bang nucleosynthesis limits and remain visible to upcoming experiments like SHiP. Read it to see a concrete way that accelerator searches and early-universe observations could jointly measure the lepton asymmetry of the Universe.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Mechanism hinges on exact lepton-number conservation (δm ≲ ℏ/τ_N ≈ 10^-15 eV), an untested model-level premise that the cosmological pipeline cannot validate.","rationale":"I read the paper in good faith and find the mechanism physically plausible and the numerical pipeline extensive, with explicit statements of structural approximations and model dependence. The central logic — that an excess of π+ from an asymmetric Dirac HNL population can push the neutron fraction back to its standard value — is simple and internally consistent. The paper's own caveats about the momentum-averaged flavor treatment and the dependence on the deuterium prediction are real but do not threaten the existence of the mechanism; they affect the quantitative boundaries. The single most load-bearing concern is the requirement of exact lepton-number conservation. The paper acknowledges it directly in 'The laboratory' section, and the reader correctly identified it as the weakest assumption. The cosmological calculation cannot test this premise; it is an input. The model example provides a tree-level realization with exact U(1)_L, but the naturalness of suppressing δm below 10^-15 eV is not established. This justifies the CONDITIONAL verdict: the mechanism is viable if the HNL is sufficiently Dirac, but that condition is not established by the BBN calculation itself. I do not see a stronger internal objection, and I agree with the reader's assessment. No verdict change is needed.","tokens_in":39207,"tokens_out":13538,"duration_ms":142548,"concrete_test":"Compute the radiatively generated Majorana mass for the HNL N in the SuM F model: one-loop contributions from the y_α Yukawa couplings and the dimension-5 Weinberg operator (HH L L)/Λ. Determine the maximum Λ (or the required cancellations) such that δm ≲ ℏ/τ_N ≈ 10^-15 eV for τ_N in the 0.02–1 s range. If Λ must exceed the Planck scale or severe fine-tuning is needed, the exact-lepton-number assumption is not natural, and the scenario's applicability is sharply limited.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim — that a hadronically decaying HNL can evade the BBN bound via charge-asymmetric pion injection — requires the HNL to be Dirac with a conserved lepton number so precise that N→N̄ oscillations do not erase the asymmetry before decay. The paper states this in 'The laboratory' section: any Majorana mass splitting must satisfy δm ≲ ℏ/τ_N ≈ 10^-15 eV. If δm is larger, the injected pion charges become symmetric and the neutron fraction is driven to the isospin plateau, reproducing the standard BBN bound. This is not an internal inconsistency; the paper explicitly frames the scenario for Dirac HNLs. However, it is a load-bearing external premise. The cosmological calculation assumes exact lepton-number conservation and cannot test it. The model example (SuM F) imposes an exactly conserved global U(1)_L, but the naturalness of keeping δm below 10^-15 eV against radiative corrections and Planck-suppressed operators is not demonstrated. If any realistic UV completion generates δm > ℏ/τ_N, the mechanism fails entirely for that completion. Thus the existence claim is conditional on a precisely conserved lepton number, a condition that lies outside the scope of the BBN calculation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a mechanism by which hadronically decaying heavy neutral leptons (HNLs) with lifetimes τ_N ≳ 0.02 s can evade the otherwise near-universal BBN bound. The essential idea is that a Dirac HNL decays into π+ and its antiparticle into π−, so an asymmetry between the N and antiN populations injects an excess of π+, which converts neutrons back to protons and restores the neutron fraction to its standard-BBN trajectory. The same large lepton flavor asymmetries raise N_eff, while HNL decays lower it, and the authors show that the two effects can cancel. The calculation uses a three-stage pipeline: momentum-resolved production of HNLs with exact plasma charge redistribution, a quantum-kinetic flavor-evolution stage, and a momentum-resolved transport stage matched to a BBN network. The authors map the open parameter space in (m_N, τ_N), finding a substantial region within SHiP reach, and identify correlated signatures: an N_eff deficit, an enhanced relic neutrino number density, and a possible first-order QCD transition sourcing gravitational waves. A model example with an exactly conserved global lepton number is provided. The authors also argue that an earlier proposal for relaxing BBN bounds with Majorana HNLs and lepton asymmetries fails once meson-driven p↔n conversion is included.","tokens_in":39470,"tokens_out":9792,"duration_ms":96159,"significance":"If the central claim holds, the result overturns a standard cosmological exclusion for sub-GeV HNLs and turns a future accelerator discovery into a probe of primordial lepton flavor asymmetries. The paper is unusually careful for a scenario paper: the core cancellation is presented transparently in Eqs. (1)–(2); the numerical pipeline passes multiple consistency checks (charge conservation, zero-asymmetry limits reproducing the standard τ_N ≲ 0.02 s bound, convergence in step size and particle number); and the structural approximations of the hybrid QKE/transport treatment are stated explicitly and, in several places, bounded by alternative implementations. The correlated predictions—N_eff deficits, CνB enhancement, and a gravitational-wave target—are falsifiable. The paper also fairly identifies the model-level premise on which the scenario rests: exact or nearly exact lepton-number conservation with δm ≲ 10^-15 eV. If that premise is accepted, the work constitutes a significant and well-validated advance in the phenomenology of heavy neutral leptons.","major_comments":[{"comment":"The central claim that hadronically decaying HNLs can evade the BBN bound is conditional on the Dirac-ness of the HNL: the paper states in 'The laboratory' that any Majorana mass splitting must satisfy δm ≲ ℏ/τ_N ≈ 10^-15 eV for the particle–antiparticle asymmetry to survive until decay. This is a load-bearing external premise: if δm is larger, N↔antiN oscillations erase the asymmetry and the neutron fraction is driven to the isospin plateau, reproducing the standard BBN bound. The model example in SuM F posits an exactly conserved global U(1)_L, but no argument is given for why such a symmetry is natural under radiative corrections or Planck-suppressed operators, and no estimate of the induced δm is provided. I recommend that the authors add a dedicated discussion of the naturalness of δm (e.g., an approximate symmetry with small explicit breaking, a gauged B-L symmetry, or a discrete remnant) or, failing that, explicitly quantify the tuning and state the implications for the scenario if any realistic completion generates δm above the quoted bound.","section":"The laboratory; SuM F"},{"comment":"The final acceptance rule for the domain classification is that a candidate is accepted if it is 'not excluded at one standard error by any window,' i.e., if the 1σ uncertainty interval overlaps the window (SuM C). Since the boundaries of Regions 1–4 in Fig. 1 and the SHiP-reach claims are built from points selected by this rule, it may systematically enlarge the plotted domains relative to a criterion that requires the mean predicted values to lie inside the windows. This is particularly relevant for the δD and ΔN_eff windows, whose widths are comparable to the statistical errors of the benchmarks. Please reclassify the grid with the mean-inside-window criterion, or alternatively report the number and fraction of accepted points whose central values fall outside one or more windows, and state the effect on the final region boundaries.","section":"SuM C (Statistics of the transport)"}],"minor_comments":[{"comment":"The phrase 'the primordial flavor asymmetries are therefore an output of the scenario' is potentially misleading: in the calculation (L_e, L_mu, L_tau) are inputs to the scan, and a hypothetical detection would constrain them along a one-parameter family. Please rephrase to avoid implying that the model itself predicts the asymmetries.","section":"The laboratory"},{"comment":"The estimate r_req ≃ 4–10 in Eq. (2) is usefully contrasted with the actual benchmark q_N values in Table I; I suggest adding a sentence in the main text clarifying that q_N = (r−1)/(r+1) for the quoted benchmarks corresponds to r values below 4, which the full calculation accommodates through the ξ_e shift and additional meson channels.","section":"Eq. (2) and Table I"},{"comment":"The overlap of Regions 2 and 4 in the figure would be easier to read if the hatching and color legend explicitly distinguished the overlap region, since the text states that Regions 2 and 4 overlap.","section":"Fig. 1"},{"comment":"The convergence test at strong degeneracy (0.55 GeV, 0.17 s) shows ΔN_eff rising from 0.09 to 0.24 as the particle number increases from 10^4 to 6.4×10^4 per species, with a residual of ≲0.02 at the adopted 3.2×10^4 setting. This is acceptable, but it would be helpful to state explicitly whether the quoted benchmark ΔN_eff values assume the same sign for that residual at all deficit-window points.","section":"SuM D"},{"comment":"The assertion that 'the two conditions cannot align in any pattern' rests on a qualitative argument, since the general-pattern calculation is explicitly left for future work; I suggest softening the wording to 'are not expected to align' to avoid overclaiming beyond the computed cases.","section":"SuM J"}],"recommendation":"major_revision","confidential_remarks":"This is a well-posed scenario paper with a strong numerical pipeline and an honest presentation of approximations. The principal risk is the exact lepton-number-conservation premise: the authors disclose the δm ≲ 10^-15 eV condition but do not naturalize it, and the model example relies on an exactly conserved global symmetry without discussing quantum-gravity or radiative corrections. This is fixable by a dedicated naturalness discussion or an explicit framing as a tuned but observable benchmark. The acceptance-criterion issue is also fixable and should be addressed before publication. I would not reject: the central mechanism is a legitimate existence proof conditional on clearly stated model assumptions, and the paper makes a significant contribution to the phenomenology of sub-GeV HNLs."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The one thing to know: the central claim is real. A Dirac HNL population with unequal N and anti-N abundances injects more π+ than π−, the strong meson reactions convert neutrons back to protons, and the neutron fraction returns to its standard trajectory. The analytic limit in Eqs. (1)-(2) is transparent (r_req ≃ 4–10, with C_π > 1 from Coulomb attraction in the π−p channel), and the three-stage pipeline confirms compensating solutions exist across an extended (m_N, τ_N) range, including a SHiP-reachable region that standard cosmology closes.\n\nWhat's new: Gelmini et al. 2020 used Majorana HNLs, where N = anti-N so pion charges stay symmetric, and omitted the meson-driven conversion — the dominant HNL effect on BBN in this mass range. The authors' comparison with that reference is direct and fair. They also frame a sub-GeV HNL discovery as a probe of the primordial flavor asymmetries, which is sharp.\n\nThe numerics deserve credit. Charge conservation, the zero-asymmetry limit reproducing τ_N ≲ 0.02 s, the SBBN baseline, step and particle-number convergence, an alternative matching rule between the flavor and transport stages, independent comparisons with Refs. [27] and [63] — the checks are done. They also state plainly what they did not compute: no GW spectrum, no nucleation or wall dynamics, Region 4 only a conditional target. That honesty is in the text.\n\nSoft spots, in proportion. The scenario needs lepton number conserved to δm ≲ ℏ/τ_N ≈ 10⁻¹⁵ eV so N↔anti-N oscillations do not erase the asymmetry before decay. The paper says this openly in \"The laboratory\" section, but it is a load-bearing particle-physics premise the BBN calculation cannot test. The model example imposes exact global U(1)_L; the naturalness of keeping δm that small is not demonstrated, and any UV completion generating a larger Majorana splitting kills the mechanism there. Second, no code or data are released, and a key transport ingredient sits in a companion paper \"to appear\". Third, the quantitative boundaries carry disclosed choices: the adopted deuterium prediction sets the preferred sign of ΔN_eff, acceptance uses a one-sigma overlap, and the QCD criterion has a calibrated factor 3.79, a 20% margin, and assumed slopes. To their credit, they re-run the classification under the alternatives and the domains survive. The mechanism does not hinge on the optics.\n\nWho this is for: anyone in sub-GeV HNL phenomenology, beam-dump forecasts, BBN with decaying relics, or lepton asymmetries. It deserves a serious referee. I'd send it to review; I'd want code and the companion paper before treating the quantitative maps as final.","headline":"A Dirac-HNL mechanism that evades the hadronic BBN bound via charge-asymmetric pion injection; the existence claim is solid, cross-checked, and honestly caveated — send it to referees.","tokens_in":40075,"tokens_out":5105,"would_cite":true,"duration_ms":45571,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Large lepton flavor asymmetries lift the Big Bang nucleosynthesis bound on sub-GeV heavy neutral leptons, opening lifetimes up to about a second to accelerator searches.","keywords":["heavy neutral leptons","Big Bang nucleosynthesis","lepton flavor asymmetry","Dirac neutrinos","effective number of relativistic species","cosmic neutrino background","QCD phase transition","gravitational waves"],"falsifier":"Measure the mass splitting between the two states of a GeV-scale HNL, for instance by searching for lepton-number-violating decays that would signal a non-Dirac component: if the splitting exceeds $\\hbar/\\tau_N\\simeq10^{-15}$ eV for a candidate with $\\tau_N$ in the opened window, the particle–antiparticle asymmetry is erased before decay and the standard BBN bound returns.","tokens_in":38910,"feed_emoji":"🌌","tokens_out":12391,"duration_ms":109248,"temperature":0.7,"pith_summary":"This paper claims that a universe with large lepton flavor asymmetries can evade the standard Big Bang Nucleosynthesis bound on hadronically decaying heavy neutral leptons (HNLs), hypothetical singlet fermions that mix with ordinary neutrinos and are primary targets of upcoming accelerator searches. Normally, mesons from HNL decays convert protons into neutrons, overproducing helium and excluding lifetimes $\\tau_N\\gtrsim0.02$ s almost regardless of abundance. The paper shows that if the HNL is Dirac and inherits an asymmetry between particles and antiparticles, its decays inject more $\\pi^+$ than $\\pi^-$, converting neutrons back to their standard abundance, while the flavor asymmetries and decay heating nearly cancel in the effective number of relativistic neutrino species $N_{\\rm eff}$. This opens the window $m_\\mu+m_\\pi\\lesssim m_N\\lesssim1$ GeV with $\\tau_N\\lesssim1$ s, much of it within reach of the next generation of beam-dump experiments, so an HNL discovery there would signal a non-standard MeV-temperature universe and would pin down the primordial flavor asymmetries up to one direction.","feed_headline":"Lepton asymmetry reopens sub-GeV heavy neutrinos to accelerators","feed_subtitle":"Asymmetric Dirac decays inject extra pions, keeping helium and the neutrino energy density standard.","key_machinery":"The central object is the particle–antiparticle ratio $r\\equiv n_N/n_{\\bar N}$ of the HNL population. While charged pions are present, the neutron fraction tracks the quasi-static equilibrium $X_n \\simeq C_\\pi/(C_\\pi + r)$, with $C_\\pi = \\langle\\sigma v\\rangle_{\\pi^- p\\to n}/\\langle\\sigma v\\rangle_{\\pi^+ n\\to p} = \\mathcal O(1)$; choosing $r_{\\rm req} = C_\\pi(1-X_n^{\\rm SBBN})/X_n^{\\rm SBBN} \\simeq 4$–$10$ pins the neutron fraction to the standard track once the mesons disappear. The mechanism needs the charged-pion channel to track the parent: $N\\to\\pi^+ l^-$ while $\\bar N\\to\\pi^- l^+$, which holds cleanly below about $1$ GeV. Lepton number conservation with $\\delta m\\lesssim\\hbar/\\tau_N$ prevents $N\\leftrightarrow\\bar N$ oscillations from erasing the imbalance before decay.","core_discovery":"The paper's central claim is that the meson-driven BBN bound, which excludes hadronically decaying relics with lifetimes $\\tau\\gtrsim0.02$ s almost independently of their abundance, can be lifted for Dirac heavy neutral leptons that inherit part of the lepton flavor asymmetries of the Universe, without modifying their laboratory decay rates. Because the HNL decays as $N\\to\\pi^+l^-$ and its antiparticle as $\\bar N\\to\\pi^-l^+$, an asymmetry in the population injects an excess of $\\pi^+$, and the strong pion–nucleon conversions return the neutron fraction to its standard value; a cancellation between the asymmetry-driven increase and the decay-driven decrease of $N_{\\rm eff}$ keeps the CMB observable small. The paper demonstrates, with a three-stage momentum-resolved evolution, that these solutions exist over $m_\\mu+m_\\pi\\lesssim m_N\\lesssim1$ GeV and $\\tau_N\\lesssim1$ s, including a substantial region within the reach of upcoming accelerator searches, and that a discovery there would fix the primordial flavor asymmetries up to one remaining direction, which the relic neutrino background and a possible first-order cosmic QCD transition may probe.","pith_inferences":["The paper does not apply the mechanism to other relics, but any hadronically decaying species with a conserved particle–antiparticle asymmetry and charge-correlated meson final states would evade the BBN bound the same way, making the essential condition the asymmetry, not the neutrino nature of the particle.","The paper does not discuss a laboratory readout of the asymmetry, but the charge-tracking relation suggests one: measuring whether $N$ decays preferentially produce $\\pi^+$ or $\\pi^-$ could directly expose the sign of the primordial lepton asymmetry that created the population.","The paper does not claim the current CMB preference for $N_{\\rm eff}<3$ is explained, but the negative $\\Delta N_{\\rm eff}$ realized in part of the opened region gives a concrete target for near-future CMB measurements and a physical route to such a deficit.","The paper leaves supernova constraints as a fixed input; a dedicated core-collapse supernova reanalysis in the sub-GeV mass range could materially shrink or extend the region that accelerator searches could claim."],"forward_implications":["Sub-GeV HNLs with lifetimes 0.02–1 s, including mixing angles that upcoming beam-dump and long-baseline experiments will probe, are not cosmologically excluded provided the Universe carries large lepton flavor asymmetries.","A detection in the opened window would, together with $Y_P$, D/H, and $N_{\\rm eff}$, determine two combinations of the three primordial flavor asymmetries at the percent level, leaving a one-parameter family.","The remaining direction changes the cosmic QCD epoch: some members of the family cross into a first-order transition and could source a stochastic gravitational wave background, while others leave no such signature.","The relic neutrino sea can be enhanced by more than 30% even where $|\\Delta N_{\\rm eff}|$ is tiny, so the cosmic neutrino background becomes a probe separate from the CMB.","The mechanism leaves laboratory decay rates untouched and requires a lepton-number-conserving (Dirac) HNL sector, so lepton-number-violation searches are a direct companion test."],"supporting_citations":[{"why":"supplies the strong meson-driven $p\\leftrightarrow n$ conversion cross sections that produce the BBN bound.","marker":"[3]"},{"why":"derives the pion-driven neutron conversion with Coulomb enhancement, fixing the $C_\\pi$ ratio used in the quasi-static equilibrium.","marker":"[4]"},{"why":"establishes the near-universal $\\tau\\lesssim0.02$ s BBN bound for hadronically decaying relics that the scenario evades.","marker":"[7]"},{"why":"gives a leptoflavorgenesis model that deposits large flavor asymmetries below sphaleron freeze-out.","marker":"[20]"},{"why":"provides the quantum-kinetic treatment of large lepton flavor asymmetries used in the oscillation stage.","marker":"[22]"},{"why":"derives BBN/CMB limits on total lepton number used to anchor the asymmetry windows.","marker":"[23]"},{"why":"locates the critical endpoint at which a large lepton asymmetry makes the cosmic QCD transition first order, defining the gravitational-wave target.","marker":"[28]"},{"why":"supplies the treatment of metastable decay products from long-lived particle decays in the MeV plasma.","marker":"[31]"},{"why":"defines the relic-neutrino capture observable used for the cosmic-neutrino-background enhancement.","marker":"[39]"},{"why":"computes the hadron multiplicities that determine the mass window where pion charge tracks lepton number.","marker":"[64]"}],"fun_headline_variants":["Lepton asymmetry rescues heavy neutrinos from BBN limit","Asymmetric neutrinos slip past BBN to reach SHiP","Flavor asymmetry lifts HNL window for accelerators","Heavy neutrinos escape BBN via lepton asymmetry","Evading BBN: lepton asymmetry opens sub-GeV neutrinos"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The mechanism works only if the heavy neutral lepton is a Dirac particle whose lepton number is conserved so precisely that particle–antiparticle oscillations cannot wash out the asymmetry before decay, which requires a mass splitting below $\\hbar/\\tau_N\\sim10^{-15}$ eV; otherwise the injected pion charges are symmetric and the standard BBN exclusion returns.","fun_headline_variants_meta":{"raw":{"variants":["Lepton asymmetry rescues heavy neutrinos from BBN limit","Asymmetric neutrinos slip past BBN to reach SHiP","Flavor asymmetry lifts HNL window for accelerators","Heavy neutrinos escape BBN via lepton asymmetry","Evading BBN: lepton asymmetry opens sub-GeV neutrinos"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000223,"raw_usage":{"total_tokens":1510,"prompt_tokens":1052,"completion_tokens":458,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":668,"completion_tokens_details":{"reasoning_tokens":372}},"tokens_in":668,"tokens_out":458,"duration_ms":4423,"temperature":1.0,"reasoning_tokens":372,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:11:39.959848+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the mass splitting between the two states of a GeV-scale HNL, for instance by searching for lepton-number-violating decays that would signal a non-Dirac component: if the splitting exceeds $\\hbar/\\tau_N\\simeq10^{-15}$ eV for a candidate with $\\tau_N$ in the opened window, the particle–antiparticle asymmetry is erased before decay and the standard BBN bound returns.","supporting_citations":[],"review_version":1}