{"id":"f398d6ac-384b-452a-9bc6-0e0ddf5cbc3e","arxiv_id":"2504.14918","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Future neutrino detectors DUNE and Hyper-Kamiokande could detect magnetic monopoles through Callan-Rubakov processes, with projected flux sensitivities reaching an order of magnitude below the Parker bound for relativistic monopoles.","lead":"DUNE and Hyper-Kamiokande could detect magnetic monopoles by looking for very high-energy antiprotons created when monopoles hit electrons, and by looking for monopole-catalyzed proton decay. If the expected Callan-Rubakov reactions happen at their maximal rates, these detectors could see monopole fluxes below the Parker bound, the standard astrophysical limit.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (3.2) counts every electron in the 1 m rock shell as a detected target, yet the paper's own FLUKA study reports only 0.41 antiprotons (or 2.02 charged hadrons) per primary reach the detector, so the shell enhancement and flux limits are not propagated through an acceptance factor.","rationale":"The reader's weakest assumption identifies exactly the same gap: the shell target electrons are counted at 100% signal efficiency in Eq. (3.2), despite the FLUKA result in Section 3.2 showing that only a fraction of produced antiprotons (or a small number of charged hadrons per primary) actually reach the detector. This is the most load-bearing concern for the antiproton-production half of the central claim because it directly enters the only event-rate formula and the quoted flux limits. It is also an internal inconsistency rather than a disagreement with external consensus: the paper itself supplies the Monte Carlo numbers that should be used as an acceptance factor. The non-relativistic proton-decay recast is independent and appears sound given the benchmark cross-section and stated lifetime sensitivities. The background-free assumption is another simplification, but the paper's 30 GeV muon-energy cut and normalised background estimate make that issue at least arguable; the shell-acceptance gap is concrete and quantitative. Because the correction is likely to weaken but not erase the order-of-magnitude claim below the Parker bound, the appropriate disposition is to keep the reader's CONDITIONAL verdict rather than rejecting the paper.","tokens_in":16014,"tokens_out":21389,"duration_ms":211094,"concrete_test":"Recompute the DUNE and Hyper-K event rates and the 90% C.L. flux limits of Fig. 5 with Eq. (3.2) modified to N = N_det σ Φ t_run + ε_shell N_shell σ Φ t_run, where ε_shell is taken first as 0.41 (antiproton survival from FLUKA) and second as the FLUKA-based probability that at least one charged hadron above 30 GeV enters the detector, e.g. 1 − exp(−2.02) under a Poisson assumption. Then compare the resulting limits with the Parker bound and with the paper's Fig. 5. As a further check, rerun FLUKA with antiproton production vertices uniformly distributed through the 1 m shell rather than a beam incident on the outer face, and use the average survival fraction as ε_shell.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The relativistic-monopole flux limits are driven by the event-rate formula N = N_target t_run Φ ∫ (dσ/dt) dt in Eq. (3.2). In Section 3.2, N_target is taken to include the 1 m crust shell, adding 2.58e34 electron targets for DUNE and 1.62e34 for Hyper-K. This treats every shell interaction as producing a detected signal with efficiency 1. However, the FLUKA simulation in the same section, for 861.4 GeV antiprotons traversing 1 m of silicon, 2 m of air, and then the LAr detector, reports only 0.41 antiprotons per primary arriving with E ≥ 30 GeV, and 2.02 charged hadrons per primary in total. No acceptance, survival probability, or per-interaction detection efficiency multiplies the shell contribution to N_target in Eq. (3.2). For DUNE the shell is roughly 70% of the total N_target, so a per-interaction efficiency in the plausible range 0.4–0.9 shifts the claimed 238% enhancement and the resulting exclusion bands by factors of order 1.7–3.4. Hyper-K is less affected because its shell is only about 26% of N_target. The paper's own simulation therefore contains the correction needed, but that correction is not propagated into the rate calculation; the 'perfect detector' assumption stated in Section 3.3 cannot recover particles that never enter the fiducial volume.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that future neutrino detectors DUNE and Hyper-Kamiokande can constrain magnetic monopole fluxes through two Callan-Rubakov signatures: high-energy antiproton production by relativistic monopoles (energies near 900 GeV) and monopole-catalysed proton decay by non-relativistic monopoles. Section 2 presents the cross-sections, kinematics, and the result that the outgoing antiproton carries at least 861.4 GeV in the laboratory frame for any allowed scattering angle and energy. Section 3 estimates atmospheric-neutrino backgrounds, introduces a 1 m rock-shell expansion of the target volume, and defines the signal-event formula N = N_target t_run Φ ∫ (dσ/dt) dt. Section 4 derives exclusion regions in the monopole mass-energy plane and recasts proton-decay lifetime sensitivities into flux limits, Φ ≲ 2.3×10^-23 cm^-2 s^-1 sr^-1 for Hyper-K and 1.1×10^-22 for DUNE. The central claim is that these experiments can reach fluxes about an order of magnitude below the Parker bound in the relativistic mass range 10^5–10^11 GeV.","tokens_in":16360,"tokens_out":9965,"duration_ms":99236,"significance":"If the analysis withstands scrutiny, the paper provides a novel and interesting application of Callan-Rubakov processes to large future detectors, with a distinctive high-energy antiproton signature and an explicit kinematic floor at 861.4 GeV. The cross-section input is not purely ad hoc: Eq. (2.3) is consistent with the earlier Kazama-Yang-Goldhaber result and with the relativistic pairwise-helicity formalism, so the central physics is independently grounded. The proton-decay catalysis limits are a simple recast of projected lifetime sensitivities, and the forward-going kinematic difference from GUT proton decay is a useful experimental handle. The paper is generally well structured and the simulations use publicly available tools, which aids reproducibility.","major_comments":[{"comment":"The shell contribution to N_target is entered with unit weight, but the paper's own FLUKA simulation in Section 3.2 reports that only 0.41 antiprotons per primary (or 2.02 charged hadrons per primary, including secondaries) reach the LAr detector after traversing 1 m of silicon and 2 m of air. Since the 1 m shell supplies about 70% of the DUNE target electrons (2.58×10^34 out of 3.67×10^34), counting all shell targets with weight 1 is not equivalent to assuming a perfect detector: it assumes that every antiproton produced anywhere in the shell enters the active volume with full efficiency. If the signal is a single high-energy antiproton, the DUNE enhancement drops from the claimed 238% to roughly 100% once the 0.41 survival probability is applied; if instead every charged hadron above 30 GeV is counted as a separate signal, the shell weight should be 2.02, and then the background estimate must additionally account for muon-induced hadronic secondaries. In either case Eq. (3.2) must be modified to fold the propagation probability into the shell term, and the exclusion bands and flux limits in Fig. 5 and the abstract need to be recomputed. This is load-bearing for the paper's headline sensitivity claim.","section":"Secs. 3.2–3.3, Eq. (3.2)"}],"minor_comments":[{"comment":"The Super-Kamiokande comparison contains an exponent typo: with τ/B > 1.6×10^34 yr, Eq. (4.2) gives a flux limit of order 1.1×10^-22 cm^-2 s^-1 sr^-1, not 1.1×10^-23 as printed; the sentence is also internally inconsistent with the statement that this is one order of magnitude larger than the Hyper-K limit.","section":"Sec. 4.2, after Eq. (4.2)"},{"comment":"The flux convention is used inconsistently: the abstract quotes limits in cm^-2 s^-1 sr^-1, while Section 4.1 fixes Φ = 4π×10^-16 cm^-2 s^-1; the authors should specify clearly in each place whether Φ is per steradian or integrated over the sphere.","section":"Abstract and Sec. 4.1"},{"comment":"The statement that the Hyper-K case follows from the DUNE FLUKA simulation because 'the same will happen for Hyper-K' is not self-evident, since water, LAr, and the surrounding geometry differ; a separate propagation study for water (or a clear argument for equivalence) would be more convincing.","section":"Sec. 3.2"},{"comment":"The sentence reporting '2.02 charged hadrons/primary' should state explicitly whether all of these charged hadrons satisfy the 30 GeV energy cut, since the preceding sentence appears to say they do but the wording is ambiguous.","section":"Sec. 3.2"},{"comment":"The caption describes bars representing fluxes that produce between 10 and 100 events, but the text does not specify what the vertical extent of each bar means in statistical terms or whether the endpoints are computed with the same target model used in the exclusion contours.","section":"Fig. 5, bottom-right panel"}],"recommendation":"major_revision","confidential_remarks":"The acceptance-factor problem is fixable within the paper's scope: the authors should either apply the FLUKA-derived survival probability to the shell contribution in Eq. (3.2) or present the shell and detector contributions separately. The proton-decay recast is straightforward and likely correct, so the overall direction is sound; however, the DUNE relativistic-monopole limit as stated is not supported until the shell propagation is folded in. I would not reject the paper on this basis."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper has a real new idea — using Callan-Rubakov antiproton synthesis as a monopole search at DUNE and Hyper-K — and a nice kinematic observation: the outgoing antiproton always carries at least ~861 GeV regardless of scattering angle or monopole energy. That's worth knowing. But the event-rate calculation has a specific internal inconsistency. Section 3.2 uses FLUKA to argue that the 1 m rock shell is a reasonable target expansion, and the simulation finds only 0.41 antiprotons per primary arrive at the detector with E ≥ 30 GeV. Then Eq. (3.2) counts every electron in the shell with 100% acceptance, no survival factor. For DUNE the shell is ~70% of the targets, so the claimed 238% enhancement and the corresponding flux limits are optimistic by roughly a factor of 1.7 if you apply their own 0.41 number. The actual average survival for antiprotons produced uniformly in the shell is probably a bit higher than 0.41, but the qualitative point stands. This doesn't kill the paper — the projected sensitivity would still sit around or below the Parker bound — but the exact limits and the 'enhanced volume' claims need correction.\n\nWhat's genuinely good: the cross-section they use is from Ref. [23] (Khoze, a co-author), but they show Eq. (2.3) is consistent with the old KYG result and with Csáki et al., so the central physics isn't circular. The ~861 GeV plateau is a clean, new result. The proton-decay catalysis recast is straightforward and reproducible from the quoted proton-lifetime sensitivities. The paper is clearly written and the derivation is transparent.\n\nSoft spots: (1) the missing acceptance factor above; (2) the background-free assumption with 100% efficiency is optimistic, though for a first sensitivity estimate it's a defensible simplifying choice; (3) there are small internal inconsistencies in the text: the intro quotes 1.5e-23 for Hyper-K's proton-decay flux while the abstract and summary say 2.3e-23, and the parenthetical about Super-K's limit (1.1e-23) doesn't match a direct calculation (should be ~1.1e-22). These are minor but should be fixed. The comparison to existing IceCube limits for beta < 0.995 is also worth a sentence, to position the new channel properly.\n\nBottom line: the central argument holds up qualitatively. The paper deserves a proper referee — it's a serious, well-grounded contribution to monopole searches with a real new channel and a real kinematic insight, but the rate calculation needs to be revised to propagate the survival probability.","headline":"A genuinely new monopole signature with a clean kinematic result, but the shell target counting ignores the paper's own FLUKA survival fraction and overstates the DUNE sensitivity by a factor of ~1.7.","tokens_in":16903,"tokens_out":8779,"would_cite":true,"duration_ms":78114,"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":"DUNE and Hyper-Kamiokande can probe magnetic monopoles through 900 GeV antiprotons and catalysed proton decay.","keywords":["magnetic monopoles","Callan-Rubakov effect","proton decay catalysis","antiproton production","DUNE","Hyper-Kamiokande","Parker bound","monopole flux limits"],"falsifier":"Recompute the event rate with the shell's targets weighted by a survival or acceptance factor of 0.41 antiprotons per primary (or 2.02 charged hadrons per primary) rather than 100%, and check whether the 90% C.L. flux limits still sit below the Parker bound. A dedicated simulation that generates 861.4 GeV antiprotons uniformly throughout the 1 m shell and counts those entering the fiducial volume above 30 GeV would settle the point directly.","tokens_in":15798,"feed_emoji":"🧲","tokens_out":10491,"duration_ms":87765,"temperature":0.7,"pith_summary":"This paper argues that two planned neutrino observatories, DUNE and Hyper-Kamiokande, can double as magnetic monopole detectors. The mechanism is the Callan-Rubakov effect, in which a monopole scattering off an electron converts it into an antiproton with at least 861.4 GeV of energy in the laboratory frame, or alternatively catalyses proton decay when the monopole is slow. For relativistic monopoles the authors derive projected flux limits around $\\Phi \\lesssim 10^{-16}\\,\\mathrm{cm^{-2}\\,s^{-1}\\,sr^{-1}}$, an order of magnitude below the Parker bound, after 15 years of running. For non-relativistic monopoles they recast proton-decay lifetime sensitivities into flux limits of $\\Phi \\lesssim 2.3\\times10^{-23}$ (Hyper-K) and $\\Phi \\lesssim 1.1\\times10^{-22}$ (DUNE) $\\mathrm{cm^{-2}\\,s^{-1}\\,sr^{-1}}$. The kinematics restrict the relativistic sensitivity to monopole masses between $10^5$ and $10^{11}$ GeV, and the fixed antiproton energy lies far above the atmospheric-neutrino-induced muon background, making the search effectively background-free.","feed_headline":"DUNE and Hyper-K could catch monopoles below the Parker bound","feed_subtitle":"The antiproton channel alone reaches ~10^-16 cm^-2 s^-1 sr^-1, an order of magnitude deeper than the Parker bound.","key_machinery":"The load-bearing machinery is the Callan-Rubakov cross-section for antiproton synthesis, Eq. (2.3): $\\frac{d\\sigma}{d\\Omega} = \\frac12 \\frac{|\\mathbf{p}^{\\rm cm}_p|}{|\\mathbf{p}^{\\rm cm}_e|} \\frac{q_J^2}{|\\mathbf{p}^{\\rm cm}_e|^2} \\left[\\sin(\\theta^{\\rm cm}/2)\\right]^{4|q_J|-2}$, where $q_J$ is half the monopole magnetic charge in units of $2\\pi\\hbar/e$. It is IR-dominated, with size set by the GeV strong-interaction scale rather than the GUT scale, and it combines with the kinematic threshold $E_{\\rm th}\\simeq (m_p/m_e)m_M\\simeq 1.8\\times10^3 m_M$ to enforce the relativistic condition $\\gamma \\gtrsim 2000$ and the $10^5$–$10^{11}$ GeV mass window. For slow monopoles, the companion catalysed proton-decay cross-section Eq. (2.5) governs the $M+p\\to M+e^++\\pi^0$ channel used in the proton-lifetime recast.","core_discovery":"The central discovery is that a monopole–electron Callan-Rubakov scattering always produces an antiproton of at least 861.4 GeV in the laboratory frame once the reaction is kinematically allowed, regardless of the scattering angle or the monopole mass and energy. This fixed high energy places the signal far above the dominant atmospheric-neutrino muon background, whose energies are at most about 20 GeV. The cross-section, Eq. (2.3), scales as $q_J^2/(\\mathrm{GeV}^2)$ and is not suppressed by the GUT scale, so the event rate is governed by the strong-interaction scale. With a 1 m rock shell included as target, DUNE gains 238% in effective electron targets and Hyper-K gains 25%; the projected 90% C.L. sensitivities are $\\Phi \\lesssim 10^{-16}\\,\\mathrm{cm^{-2}\\,s^{-1}\\,sr^{-1}}$ for relativistic monopoles after 15 years, and $\\Phi \\lesssim 2.3\\times10^{-23}$ (Hyper-K) or $1.1\\times10^{-22}$ (DUNE) $\\mathrm{cm^{-2}\\,s^{-1}\\,sr^{-1}}$ for non-relativistic monopoles through catalysed proton decay.","pith_inferences":["The paper's event-rate formula counts every electron in the 1 m shell as a 100%-efficient target, even though its own simulation reports only 0.41 antiprotons (or 2.02 charged hadrons) reaching the detector per primary; applying that survival factor would shrink the shell contribution and push the flux limits upward, by roughly a factor of two for DUNE.","The same simulation method could optimise the shell thickness: a thicker shell adds targets but reduces survival probability, so the effective detection volume is maximised at some intermediate thickness rather than at exactly 1 m.","If the antiproton plateau at 861.4 GeV is confirmed, the energy sideband above about 30 GeV in any large underground detector becomes a monopole search window, extending this analysis to detectors not discussed in the paper.","The recast proton-decay limits scale linearly with the assumed catalysis cross-section; if a more realistic cross-section is smaller than $\\sigma = 1.4\\times10^{-21}\\,\\mathrm{cm^2}$, the flux limits degrade proportionally."],"forward_implications":["A null observation of the antiproton signal at DUNE or Hyper-K after 15 years would exclude relativistic monopole fluxes at the level of $\\Phi \\sim 10^{-16}\\,\\mathrm{cm^{-2}\\,s^{-1}\\,sr^{-1}}$, an order of magnitude below the Parker bound.","The fixed 861.4 GeV antiproton energy separates the signal kinematically from the dominant atmospheric-neutrino background, whose muons peak below about 20 GeV, so a 30 GeV energy cut leaves the search essentially background-free.","For non-relativistic monopoles, non-observation of $p\\to e^+\\pi^0$ at Hyper-K would constrain fluxes to $\\Phi \\lesssim 2.3\\times10^{-23}\\,\\mathrm{cm^{-2}\\,s^{-1}\\,sr^{-1}}$, comparable to existing Super-K limits, while DUNE would reach $\\Phi \\lesssim 1.1\\times10^{-22}$.","Monopole-catalysed proton decay, if seen, would present forward-going decay products rather than the back-to-back topology of GUT gauge-boson-mediated decay, providing a kinematic discrimination handle.","The 1 m rock-shell expansion raises the effective target-electron count by 238% at DUNE and 25% at Hyper-K, effectively turning surrounding material into part of the detector for this signature."],"supporting_citations":[{"why":"Supplies the Callan-Rubakov cross-section formulae (Eqs. 2.3 and 2.5) on which all event-rate calculations rest.","marker":"[23]"},{"why":"Provides the QED monopole-electron scattering result that the branching-ratio-stripped version of the cross-section must reproduce.","marker":"[24]"},{"why":"Provides Hyper-Kamiokande's fiducial volume and the 90% C.L. proton-lifetime sensitivity used for the catalysis recast.","marker":"[17]"},{"why":"Provides DUNE's detector parameters and its proton-lifetime sensitivity used in the recast.","marker":"[40]"},{"why":"Supplies the isotropic-flux treatment and the comparison point for relativistic monopole limits from IceCube.","marker":"[9]"},{"why":"Supplies the stopping powers and energy-loss rates in rock used to determine which monopole energies survive to the detector.","marker":"[33]"},{"why":"Provides the atmospheric neutrino flux used to estimate the muon background energy distribution.","marker":"[34]"},{"why":"Used to simulate charged-current atmospheric neutrino events for the muon background in both detectors.","marker":"[35]"},{"why":"Used to simulate antiproton propagation through the shell and detector in support of the shell-expansion concept.","marker":"[36]"}],"fun_headline_variants":["DUNE and Hyper-K could spot monopoles via 900 GeV antiprotons","Monopole scattering yields >861 GeV antiprotons at DUNE","DUNE, Hyper-K probe monopole flux below Parker bound","Monopole hunt at future neutrino detectors: antiproton signal"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that every electron in the surrounding 1 m rock shell contributes to the signal with 100% efficiency: the event-rate formula Eq. (3.2) counts the full shell as target, even though the paper's own simulation finds that only 0.41 antiprotons (or 2.02 charged hadrons) per primary actually reach the detector.","fun_headline_variants_meta":{"raw":{"variants":["DUNE and Hyper-K could spot monopoles via 900 GeV antiprotons","Monopole scattering yields >861 GeV antiprotons at DUNE","DUNE, Hyper-K probe monopole flux below Parker bound","Monopole hunt at future neutrino detectors: antiproton signal"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000593,"raw_usage":{"total_tokens":2815,"prompt_tokens":1017,"completion_tokens":1798,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":633,"completion_tokens_details":{"reasoning_tokens":1732}},"tokens_in":633,"tokens_out":1798,"duration_ms":12687,"temperature":1.0,"reasoning_tokens":1732,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:37:08.993112+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the event rate with the shell's targets weighted by a survival or acceptance factor of 0.41 antiprotons per primary (or 2.02 charged hadrons per primary) rather than 100%, and check whether the 90% C.L. flux limits still sit below the Parker bound. A dedicated simulation that generates 861.4 GeV antiprotons uniformly throughout the 1 m shell and counts those entering the fiducial volume above 30 GeV would settle the point directly.","supporting_citations":[{"cited_title":"Kazama, C.N","cited_arxiv_id":null,"evidence_quote":"Provides the QED monopole-electron scattering result that the branching-ratio-stripped version of the cross-section must reproduce."},{"cited_title":"Derkaoui, G","cited_arxiv_id":null,"evidence_quote":"Supplies the stopping powers and energy-loss rates in rock used to determine which monopole energies survive to the detector."},{"cited_title":"FLUKA CERN","cited_arxiv_id":null,"evidence_quote":"Used to simulate antiproton propagation through the shell and detector in support of the shell-expansion concept."}],"review_version":1}