{"id":"fb843369-7125-461e-82e9-db979ac585b4","arxiv_id":"1908.05958","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A JUNO-like reactor neutrino detector could reach a 90% confidence sensitivity of about 2.4 antineutrino events in a 1000-second gravitational-wave coincidence window, implying fluence limits of order 10^8 cm^-2 and detectable distances of 1 to 3 Mpc.","lead":"Future liquid-scintillator neutrino detectors like JUNO could see electron antineutrinos from merging neutron stars and black holes out to about 1 to 3 million light-years, according to this simulation study. The paper estimates how sensitive these detectors would be to neutrinos arriving together with gravitational-wave signals, and how adding a second detector would improve the reach.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"BBH distance reach rests on an ad hoc 3 s emission duration; since D_UL scales as sqrt(ΔT), a 0.1 s burst would shrink the claimed 1-3 Mpc reach by about a factor of 5.","rationale":"The paper's statistical framework appears sound: the Feldman-Cousins sensitivity μ90 = 2.44 for n_obs = n_bkg in 1000 s is the standard quantity, and the 38.5% improvement from an identical extra detector follows from combining Poisson likelihoods with the same μ_C per effective proton. The MC multiplicity formula (Eq. 3.2) is derived correctly in Appendix A. The weakest link is the astronomical reach, because it is the only headline number obtained by inserting external model luminosities and an assumed emission duration into Eq. (4.4). The paper itself flags the absence of a BBH emission-time estimate and adopts 3 s for consistency, applying a hypermassive-neutron-star cooling timescale to prompt accretion-disk emission. This is not a matter of consensus but of internal applicability: the same equation is used with a duration that is not supported by the cited source models. A single numerical test, recomputing the BBH/BH-NS rows with the simulation's own emission window, would determine whether the claimed 1-3 Mpc reach survives. Until then, the conditional verdict is appropriate; the concern does not invalidate μ90 or the relative improvement from an extra detector, but it does weaken the absolute distance claim for a large fraction of GW sources.","tokens_in":16684,"tokens_out":6438,"duration_ms":59885,"concrete_test":"Recompute the BBH/BH-NS rows of Table 4 using the actual neutrino-emission duration from the cited accretion-torus simulations (Caballero et al. 2016, arXiv:1510.06011) instead of ΔT = 3 s; if the simulation reports a time-dependent luminosity, integrate L(t) over the emission window to obtain total energy rather than multiplying peak L_s by 3 s. Then recalculate D_UL via Eq. (4.4). If the BBH/BH-NS distances drop below 1 Mpc or change by more than 50% relative to Table 4, the '1-3 Mpc' reach in the abstract must be qualified as BNS-only and the BBH distance claim revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing premise is the conversion of the 90% C.L. fluence upper limit into a detectable distance in Section 4.2, Eq. (4.4): D_UL = sqrt(L_s ΔT / (4π F_UL <E>)). Table 3 lists neutrino luminosities L_s from merger simulations, and the paper multiplies these by a single assumed emission duration ΔT = 3 s for all source types. The text states explicitly that for BBH mergers 'there are no estimations on the emission time of neutrinos yet' and that 3 s is adopted 'for consistency', based on hypermassive neutron star cooling times of 2-3 s [33,34]. That cooling argument does not apply to prompt accretion-torus emission in BBH/BH-NS models, where the neutrino-emitting transient is typically of order tens to hundreds of milliseconds. Because D_UL scales as the square root of ΔT, using 3 s when the true duration is 0.1 s overstates the distance reach by a factor of sqrt(30) ≈ 5.5; with a 10 ms burst, the overestimate is a factor of about 17. The same concern applies to multiplying a peak luminosity by a duration at all, since the cited simulations provide time-dependent luminosities and the total emitted energy may be far less than L_s × 3 s. Since BBH events dominate the GW catalog, the headline 'detectable distance of about 1-3 Mpc for gravitational wave sources' is not supported for BBH/BH-NS and the stated reach should be restricted to BNS models or explicitly parametrized by ΔT.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes using the future JUNO reactor neutrino experiment (and RENO-50 as a combined or alternative detector) to search for low-energy electron antineutrinos from gravitational wave sources, motivated by the fact that no coincident neutrino candidates have been found to date. The authors simulate background in the nominal JUNO detector, study the time correlation of adjacent background events via a cut-off-time bundling method, and then compute a 90% C.L. sensitivity upper limit μ90 = 2.44 for a 1000 s coincidence window using a Feldman-Cousins style Poisson likelihood. They translate this limit into a neutrino fluence sensitivity under monochromatic and Fermi-Dirac spectrum assumptions, and then, using published merger model luminosities for BNS, BH-NS, and BBH sources, convert the fluence limits into detectable distances of about 1-3 Mpc. They also evaluate the improvement from an extra identical JUNO-like detector, claiming about 38% sensitivity improvement and about 28% distance improvement, and from a JUNO + RENO-50 combination, finding similar gains.","tokens_in":16955,"tokens_out":15084,"duration_ms":143151,"significance":"If the method is sound, it would provide a quantitative estimate of the multi-messenger reach of large liquid-scintillator detectors for low-energy neutrinos from compact binary mergers, a topic of continuing interest. The paper has some clear strengths: the analytic derivation of the multiplicity distribution of background bundles in Appendix A is explicit; the background table is concrete; and the use of the Feldman-Cousins approach for small signals is appropriate. However, the headline distance reach depends on several model-dependent choices, most notably an ad hoc 3 s neutrino emission duration for BBH and BH-NS mergers, and the manuscript contains numerical inconsistencies in the quoted fluence endpoints. The core Poisson sensitivity calculation is reproducible in principle, but the model-dependent distance claims need substantial qualification.","major_comments":[{"comment":"The distance reach D90_UL is computed from D = sqrt(L_s ΔT / (4π F_UL <E>)) using a single emission duration ΔT = 3 s for BNS, BH-NS, and BBH sources. The text acknowledges that for BBH mergers there are no estimates of the neutrino emission time and sets 3 s 'for consistency', based on hypermassive neutron star cooling times. Since D90_UL scales as sqrt(ΔT), an actual prompt emission duration of ~0.1 s (typical for accretion-torus emission in BH-NS/BBH models) reduces the claimed distance reach by a factor of about 5.5, and a 10 ms burst by a factor of about 17. In addition, the luminosities in Table 3 are peak values around 10 ms after merger; multiplying a peak luminosity by 3 s may substantially overestimate the total emitted energy. The headline '1-3 Mpc' is therefore not supported for BBH/BH-NS sources and should be reported as a function of ΔT or restricted to BNS models.","section":"Sec. 4.2, Eq. (4.4), Table 3"},{"comment":"The reported monochromatic fluence sensitivity is internally inconsistent by four orders of magnitude: the abstract and the conclusion state the range 6×10^10 to 4×10^10 cm^-2, while the paragraph after Eq. (4.5) in Sec. 4.2 gives '6×10^10 cm^-2 to about 4×10^6 cm^-2'. For a monochromatic fluence limit the energy dependence is F_UL ∝ 1/σ(E_ν) ≈ 1/E_ν^2, so the ratio of the 1.8 MeV to 120 MeV endpoints should be roughly (120/1.8)^2 ≈ 4.4×10^3; neither stated pair satisfies this. The authors must correct the endpoint, provide a sample calculation with the values of μ90, N_T, ε, and σ(E), and make the abstract consistent with the body.","section":"Abstract, Sec. 4.2, Sec. 6"},{"comment":"The paper introduces a time-correlation method (bundling adjacent background events with a cut-off time) as its new analysis technique, but the sensitivity calculation in Sec. 4.1 is a simple counting experiment in a fixed 1000 s window, not a clustering analysis. The two approaches give very different sensitivities: Sec. 3 concludes from Table 2 that a bundle of 24 events is needed for unambiguous identification, while Sec. 4.1 quotes μ90 = 2.44. This discrepancy is not explained. The authors should either use the time-correlation information in the Feldman-Cousins likelihood or clearly state that Sec. 3 is only a background-characterization study and that the quoted sensitivity comes from the counting analysis.","section":"Sec. 3 and Sec. 4.1"},{"comment":"The 90% C.L. sensitivity μ90 = 2.44 is a purely statistical limit that treats the background rate nbkg as exactly known. The background rates in Table 1 are estimates without uncertainties, and the detection efficiency ε and reactor power are used without error propagation. Systematic uncertainties on these inputs can shift the upper limit and, via Eqs. (4.3)-(4.5), the fluence and distance constraints. The authors should quantify the effect of plausible systematic variations (e.g., ±10% on the reactor neutrino rate) or justify that they are negligible for the quoted precision.","section":"Sec. 4.1, Eqs. (4.1)-(4.2), Table 1"},{"comment":"Equation (4.3) explicitly assumes no neutrino oscillation. Since the luminosities in Table 3 are source properties and the fluence F is evaluated at Earth, the relation between source ν̄e luminosity and detected ν̄e fluence should account for flavor conversion over astrophysical baselines. The oscillation factors are of order unity but can shift D90_UL by a factor of a few in either direction; this should be stated and, if a no-oscillation assumption is retained, justified as conservative.","section":"Sec. 4.2, Eq. (4.3)"}],"minor_comments":[{"comment":"There are numerous typographical errors (e.g., 'backgroud' for 'background', 'varing' for 'varying', 'revelas' for 'reveals', 'beseline' for 'baseline', 'Finially' for 'Finally', 'conincident' for 'coincident'). A careful proofread is needed.","section":"Throughout"},{"comment":"In Eqs. (A.1)-(A.4), the infinitesimal interval is written as P(1; R dt1) dt1; since P(1; R dt1) is already a probability, the extra dt1 makes the expression dimensionally inconsistent. It should be R dt1 (or P(1;R dt1) without the trailing dt1). The final result (3.2) is correct.","section":"Appendix A, Eqs. (A.1)-(A.4)"},{"comment":"The sentence 'the accumulated ν̄e signals emitted from gravitational wave sources can’t be identified from backgrounds until over 23 events within the time window of between 1000 s and 24000 s' is unclear because a signal is expected at a known GW time, so the relevant background is the rate in a much shorter interval; this should be rephrased.","section":"Sec. 3, after Table 2"},{"comment":"The best-fit value μ_best is not defined; for a Poisson likelihood with known background it should be max(0, nobs - nbkg), and this definition should be stated.","section":"Sec. 4.1, Eq. (4.2)"},{"comment":"The caption says the red curve is 'the most conservative detection capability,' but it is not explained why the monochromatic assumption is the most conservative; a sentence in the text would help.","section":"Fig. 4 caption"},{"comment":"The sentence 'Varing the baseline away from Yangjiang and Taishan NPP and target mass' is garbled; it should read 'Varying the baseline and target mass...'.","section":"Sec. 5, paragraph on extra detector scan"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a 2019 preprint and has not been updated to reflect later JUNO design information. The main technical weakness is the ad hoc 3 s emission duration applied to BBH/BH-NS sources, which directly undermines the headline distance reach; the internal inconsistency between the abstract and Sec. 4.2 fluence endpoints also needs to be fixed. If the authors restrict the distance claims to BNS cases or parametrize them by ΔT, correct the numerical errors, and clarify the relationship between the Sec. 3 clustering method and the Sec. 4 counting sensitivity, the paper could be publishable as a methods-oriented contribution. In its current form, the internal contradictions prevent a positive recommendation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. The paper is a clean, narrow sensitivity study: it gives the first quantitative projections I know of for JUNO and RENO-50 in the search for low-energy antineutrinos coincident with gravitational-wave events, including a scan of what an extra identical detector buys. The bundle-multiplicity formula (3.2) is derived in the appendix and checks out; the Monte Carlo matches it. The μ90 = 2.44, the fluence limits, and the 38%/28% gains from a second JUNO-like detector are all plausible and internally consistent, assuming the standard Feldman-Cousins/Poisson machinery and the nobs=nbkg convention.\n\nThe real soft spot is the distance conversion, not the counting statistics. The authors honestly say there are no estimates for BBH neutrino emission time and adopt ΔT=3 s \"for consistency.\" Since D_UL scales as sqrt(ΔT), that single choice controls the headline 1-3 Mpc reach for BBH/BH-NS. A 0.1 s burst would move the reach down by about a factor of 5; a few milliseconds would be worse. For BNS, the 2-3 s cooling argument is more defensible, but multiplying a peak luminosity by a duration is still a crude proxy for the time-integrated emission. The fix is easy: parametrize by ΔT, or quote total emitted energy, and restrict the unqualified headline to BNS.\n\nThere is also a four-order-of-magnitude mismatch in the fluence range: the abstract says 6×10^10 to 4×10^10 cm^-2, Section 4.2 says 6×10^10 to 4×10^6 cm^-2. Simple E^-2 scaling from the IBD cross section gives about 1×10^7 cm^-2 at 120 MeV, so both endpoints need checking. Neutrino oscillation is ignored in Eq. (4.3); at Mpc distances that is a tens-of-percent effect and should be acknowledged. No code or data artifacts are provided, which limits reproducibility, and there are no systematic uncertainties. Those are limitations, not fatal flaws for a projection paper.\n\nThe citation pattern looks normal—published JUNO background rates, standard merger simulation luminosities, earlier KamLAND/SK/Borexino searches. I would not desk-reject this. It deserves a serious referee, but I would ask the referee to require a corrected fluence range, a ΔT-parametrized distance claim, and a BBH claim that does not overreach. Then it would be a useful reference.","headline":"A useful first sensitivity projection for JUNO/RENO-50 to GW neutrinos, but the advertised 1-3 Mpc BBH reach hangs on an admitted placeholder of 3 s emission duration and a fluence-range typo.","tokens_in":17580,"tokens_out":5319,"would_cite":true,"duration_ms":53520,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A JUNO-class reactor neutrino detector can set a 90% C.L. limit of 2.44 antineutrino events per 1000 s, reaching 1-3 Mpc for merger models.","keywords":["gravitational-wave multi-messenger astronomy","electron antineutrinos","JUNO","reactor neutrino experiment","inverse beta decay","Feldman-Cousins sensitivity","neutrino fluence upper limit","binary neutron star mergers"],"falsifier":"A numerical-relativity simulation of a binary black hole merger with an accretion disk that computes the time-integrated $\\bar\\nu_e$ energy would settle the most fragile part of the distance claim: if that energy is below roughly $10^{54}$ erg, the value implied by the Table 3 luminosities times $\\Delta T=3$ s, then the claimed black-hole-merger distances overstate the reach. For the headline sensitivity itself, an independent implementation of the Feldman-Cousins construction with $n_{\\rm obs}=n_{\\rm bkg}=0.745$ events per 1000 s should reproduce $\\mu_{90}=2.44$; a materially different value would falsify the claimed limit.","tokens_in":16408,"feed_emoji":"🔭","tokens_out":10443,"duration_ms":95451,"temperature":0.7,"pith_summary":"The paper argues that future medium-baseline reactor neutrino detectors such as JUNO can serve as sensitive monitors of the low-energy electron antineutrinos that might accompany gravitational-wave transients. Using the background rate at the nominal JUNO detector, it shows that a non-detection in a 1000-second window would set a 90% confidence-level sensitivity of $\\mu_{90} = 2.44$ signal events, a fluence limit of roughly $6\\times10^{10}$ to $4\\times10^{10}$ cm$^{-2}$ for monoenergetic antineutrinos between 1.8 and 120 MeV, and about 1-3 Mpc distance reach when model luminosities for binary neutron star, black hole-neutron star, and binary black hole mergers are assumed. It then shows that a second identical detector improves the fluence sensitivity by about 38% and the distance reach by about 28%, while a JUNO plus RENO-50 combination gives nearly the same improvement. The point is that reactor neutrino experiments, built for oscillation physics, can double as multi-messenger observatories for nearby compact-object mergers.","feed_headline":"Reactor neutrino detector could trace GW mergers to 3 Mpc","feed_subtitle":"A quiet 1000-second window still sets the tightest low-energy antineutrino fluence limits for nearby mergers.","key_machinery":"The method's two devices are a Poisson time-bundling statistic and the Feldman-Cousins confidence-belt construction. With a cutoff time $\\Delta t$, the expected number of background bundles of multiplicity $m$ is $N(m)=RT\\,e^{-2R\\Delta t}(1-e^{-R\\Delta t})^{m-1}$, so a cluster of $m$ events within $\\Delta t$ can be recognized as an excess. The Feldman-Cousins construction on a Poisson likelihood ratio turns the expected background and observed candidates into the 90% C.L. upper limit $\\mu_{90}$. That limit is then converted to a fluence via the inverse $\\beta$-decay cross section and the number of target protons, and to a distance via $L=F\\,4\\pi D^{2}\\langle E\\rangle$ using model luminosities and either a monochromatic or a Fermi-Dirac spectrum.","core_discovery":"The paper's central claim is that a non-detection of $\\bar\\nu_e$ in a 1000 s window at the nominal JUNO experiment yields a 90% C.L. sensitivity upper limit of $\\mu_{90}=2.44$ events, corresponding to a fluence limit of roughly $6\\times10^{10}$ to $4\\times10^{10}$ cm$^{-2}$ for monoenergetic $\\bar\\nu_e$ between 1.8 and 120 MeV and to $F_{UL}^{90}\\sim 1$-$3\\times10^{8}$ cm$^{-2}$ for the model spectra, and that this translates to a detectable distance $D_{UL}^{90}\\sim 1$-$3$ Mpc under the BNS, BH-NS, and BBH luminosity models in Table 3. It further claims that adding a second identical JUNO-like detector improves the fluence sensitivity by 38.5% and the distance reach by 27.5%, and that combining JUNO with RENO-50 gives a 37.7% improvement, nearly the same as adding the identical detector.","pith_inferences":["Because $D_{UL}$ scales as $\\sqrt{L_s\\,\\Delta T}$, the assumed 3 s emission duration is a large lever: if neutrino emission from a merger lasted 30 s, the same detector would reach about three times farther, and if it lasted 0.3 s, about one-third as far.","The paper's figure of merit shows diminishing returns from extra target mass; beyond a few tens of kilotons, reducing backgrounds through a farther baseline or other means may buy more sensitivity than adding more mass.","If a nearby binary neutron star merger is ever observed in neutrinos, the same fluence-to-distance conversion could be inverted to measure the source's neutrino luminosity directly, testing the model predictions rather than assuming them."],"forward_implications":["During a future gravitational-wave trigger, a null result from JUNO would set the tightest low-energy $\\bar\\nu_e$ fluence limits yet from a reactor neutrino detector, about twenty times better than KamLAND.","Only a few $\\bar\\nu_e$ events arriving within tens of seconds would stand out from background, so the experiment can flag nearby mergers without needing large event statistics.","The 1-3 Mpc reach means neutrino-coincident gravitational-wave detections are limited to very nearby mergers, making each successful observation rare and the constraints from non-detection valuable.","An extra identical detector, or the planned RENO-50 detector analyzed jointly, improves the fluence sensitivity by roughly 38%, and the paper quantifies how the gain varies with target mass and baseline.","The same multiplicity-bundling test can be applied to any short-timescale transient search at future reactor neutrino experiments, not only gravitational-wave triggers."],"supporting_citations":[{"why":"Supplies the nominal JUNO detector configuration, target mass, detection efficiency, and background rates used throughout.","marker":"[18]"},{"why":"Provides the KamLAND search whose energy range, time window, and sensitivity this paper extends and compares against.","marker":"[13]"},{"why":"Supplies the Feldman-Cousins unified confidence-belt method used to compute the 90% C.L. sensitivity upper limit.","marker":"[29]"},{"why":"Supplies the inverse beta-decay cross section used to convert the signal-event upper limit into a fluence limit.","marker":"[31]"},{"why":"Provides binary-neutron-star merger neutrino luminosities and average energies used in Table 3.","marker":"[23]"},{"why":"Provides black-hole-neutron-star merger neutrino luminosity and average energy values used in Table 3.","marker":"[26]"},{"why":"Provides the accretion-disk model used for the black-hole merger luminosity entries in Table 3.","marker":"[27]"}],"fun_headline_variants":["JUNO's silence maps neutrino fluence from GWs","Extra detector boosts GW neutrino sensitivity by 38%","Non-detection sets tightest limits on GW neutrinos","Neutrino fluence limits from GW sources with JUNO","JUNO + RENO-50 improves GW neutrino reach"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The distance reach of 1-3 Mpc assumes every gravitational-wave source class, including binary black hole mergers, emits $\\bar\\nu_e$ with the model luminosities in Table 3 and for a common emission duration $\\Delta T=3$ s; the paper itself states that no estimate for the BBH neutrino emission duration exists and adopts 3 s \"for consistency.\"","fun_headline_variants_meta":{"raw":{"variants":["JUNO's silence maps neutrino fluence from GWs","Extra detector boosts GW neutrino sensitivity by 38%","Non-detection sets tightest limits on GW neutrinos","Neutrino fluence limits from GW sources with JUNO","JUNO + RENO-50 improves GW neutrino reach"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000643,"raw_usage":{"total_tokens":3077,"prompt_tokens":1185,"completion_tokens":1892,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":801,"completion_tokens_details":{"reasoning_tokens":1807}},"tokens_in":801,"tokens_out":1892,"duration_ms":13196,"temperature":1.0,"reasoning_tokens":1807,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:01:15.164355+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A numerical-relativity simulation of a binary black hole merger with an accretion disk that computes the time-integrated $\\bar\\nu_e$ energy would settle the most fragile part of the distance claim: if that energy is below roughly $10^{54}$ erg, the value implied by the Table 3 luminosities times $\\Delta T=3$ s, then the claimed black-hole-merger distances overstate the reach. For the headline sensitivity itself, an independent implementation of the Feldman-Cousins construction with $n_{\\rm obs}=n_{\\rm bkg}=0.745$ events per 1000 s should reproduce $\\mu_{90}=2.44$; a materially different value would falsify the claimed limit.","supporting_citations":[{"cited_title":"Search for electron antineutrinos associated with gravitational wave events GW150914 and GW151226 using KamLAND","cited_arxiv_id":"1606.07155","evidence_quote":"Provides the KamLAND search whose energy range, time window, and sensitivity this paper extends and compares against."},{"cited_title":"Low mass binary neutron star mergers : gravitational waves and neutrino emission","cited_arxiv_id":"1510.06398","evidence_quote":"Provides binary-neutron-star merger neutrino luminosities and average energies used in Table 3."},{"cited_title":"The black hole spin influence on accretion disk neutrino detection","cited_arxiv_id":"1510.06011","evidence_quote":"Provides the accretion-disk model used for the black-hole merger luminosity entries in Table 3."}],"review_version":1}