{"id":"55b59d2a-57e3-456f-af86-d93b68e6f1b7","arxiv_id":"2411.12344","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Hyper-Kamiokande could measure the average core density with an uncertainty around -14%/+40% at nominal resolution, improving to about -8%/+10% with better resolution, after 6500 days.","lead":"This paper calculates how well the future Hyper-Kamiokande neutrino detector could measure the density of the Earth's core and mantle layers by watching how atmospheric neutrinos oscillate as they pass through the planet. If the estimate holds, neutrino oscillation tomography could become an independent check on geophysical models like PREM.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The nominal +39.5% upper edge sits outside the paper's own M⊕/I⊕/EHE-allowed range (κC<0.10), so the abstract's headline 2σ interval is not a physical posterior interval.","rationale":"The reader's strongest claim identifies the same numerical conflict, but her weakest_assumption emphasizes missing systematics. I agree with the reader's conditional verdict but think the constraint conflict is more directly load-bearing because it affects the literal headline number and is verifiable from the paper's own equations. The favorable-case interval with sin²θ23=0.50 (upper +9.8%) is consistent with the constraints, so the paper's method is not wholly undermined; the nominal-case claim and the sin²θ23=0.45 favorable case need re-expression with the prior truncation. A corrected abstract would report the nominal 95% upper edge as ≈+10%, with the unconstrained +39.5% labeled as a statistical-only sensitivity. Verdict unchanged: conditional acceptance remains appropriate pending this correction.","tokens_in":21448,"tokens_out":7039,"duration_ms":72802,"concrete_test":"Take the nominal sin²θ23=0.45 sensitivity scan and impose the EHE inequality Eq. (9) as a hard cut while keeping the M⊕/I⊕ compensation Eq. (10). Recompute the largest κC with Δχ²≤5.99; if the scan is cut off before the Δχ² crossing, the 95% upper edge is the constraint value ≈+10%, not +39.5%. As a quick arithmetic check, evaluate Eq. (9) at κC=0.395: 6.54 g/cm³ < 2.13 g/cm³ fails, proving that the model used for the +39.5% point is unphysical.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2.1 derives hard limits on the core-density deviation: Eq. (12) gives κC < 0.089 from the EHE condition ¯ρuman(1+κuman) < ¯ρlman(1+κlman), and Eqs. (13)-(14) quote −0.33 < κC < 0.10. Despite this, the abstract and Section 4 state that in the nominal case (Eres=30%, θzres=20°, sin²θ23=0.45) HK determines ¯ρC at 2σ with κC up to +0.395. At κC=0.395, the compensating relations Eq. (10) give κuman=0.818 and κlman=−0.565, so the upper mantle density becomes 3.60×(1.818)=6.54 g/cm³ while the lower mantle becomes 4.90×(0.435)=2.13 g/cm³, violating Eq. (9) by a wide margin. Thus the positive branch of the nominal sensitivity curve is evaluated on Earth models that fail the hydrostatic-equilibrium constraint the paper claims to implement. Table 2 acknowledges this by replacing the 95% upper limit with a dash (the constraint limit), and the text concedes the discrepancy in Section 4. The abstract's quoted +39.5% is therefore not the constrained 2σ uncertainty; the actual posterior upper edge is ≈+10%. This is an internal inconsistency in the paper's central quantitative claim, not merely a conservative or optimistic choice.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a forward-model sensitivity study of neutrino-oscillation Earth tomography with the planned Hyper-Kamiokande detector. It uses the public Super-Kamiokande atmospheric-neutrino simulation release as a proxy for HK, scales it to HK's fiducial volume and 6500 days of exposure, and evaluates sensitivity to constant fractional deviations of the average core, lower-mantle, and upper-mantle densities from PREM. The analysis imposes total Earth mass, moment of inertia, and hydrostatic-equilibrium constraints, and studies how the expected 68% and 95% intervals depend on neutrino energy and zenith-angle resolution and on sin^2(theta_23), for normal mass ordering only. The headline result is that in a nominal resolution scenario HK could determine the core density at 2 sigma with an asymmetric uncertainty of about (-14.5%)/+39.5%, while a more favorable resolution scenario gives roughly (-8.3%)/+9.8%.","tokens_in":21747,"tokens_out":5249,"duration_ms":52723,"significance":"If the quoted sensitivities are correct after accounting for the Earth constraints, the paper would provide a useful quantitative projection of HK's Earth-tomography capability and would complement earlier studies for DUNE, IceCube, and ORCA. The work has several concrete strengths: it builds on a public, realistic SK simulation release; it uses quantile-weighted oscillation probabilities rather than simple bin-center estimates; it implements M_earth, I_earth, and hydrostatic-equilibrium constraints analytically and numerically; and it documents the dependence of the sensitivity on detector resolution and on sin^2(theta_23) in a compact table. The paper also honestly notes in several places that the geophysical bounds limit positive core-density deviations to about 10%. However, the central quantitative claim is presently not stated consistently: the abstract and Section 4 quote a +39.5% upper uncertainty for the nominal case that is incompatible with the constraints derived in Section 2.1, while Table 2 silently truncates the same interval.","major_comments":[{"comment":"The headline nominal-case interval (-14.5%, +39.5%) is not a physically allowed posterior interval under the constraints the paper itself imposes. Equation (14) gives -0.30 < kappa_C < 0.10 and Eq. (12) gives kappa_C < 0.089; at kappa_C = 0.395 the compensation relations in Eq. (10) give kappa_uman = 0.818 and kappa_lman = -0.565, which violate the hydrostatic-equilibrium inequality in Eq. (9). Table 2 correctly replaces the 95% upper limit with a dash for the HK Nominal rows, and Section 4 admits the constraint, but the Abstract and the main text of Section 4 still quote +39.5% as the 2-sigma uncertainty. The abstract and conclusions should quote the constrained upper edge (approximately +10%) or must explicitly label the unconstrained +39.5% value as a separate, unphysical sensitivity estimator.","section":"Abstract; Sec. 2.1, Eqs. (9)-(14); Sec. 4; Table 2"},{"comment":"The sensitivity statistic is computed as Delta chi^2 = sum_i (O_i - E_i)^2 / E_i from expectation values only, with no nuisance parameters for flux normalization, cross-section uncertainties, detector efficiencies, or energy-scale biases. The 'nominal' scenario E30%&20 deg is implemented in Sec. 3.6 only as a Gaussian smearing of reconstructed energy and angle; it is not a model of systematic uncertainties. Therefore the quoted intervals are statistical-plus-idealized-resolution projections, and the statement in Section 4 that the measurement sensitivity is 'dominated by systematic uncertainties' is not supported by the analysis as written, because no systematic uncertainties are actually varied or profiled. The authors should either introduce a nuisance-parameter treatment or explicitly restrict all claims to 'statistical and resolution-only' sensitivity.","section":"Sec. 3.7; Sec. 3.6; Sec. 4"}],"minor_comments":[{"comment":"The caption states that the dark blue bands correspond to varying kappa_C in the interval '0.30 < kappa_C < 0.10', which is internally inconsistent; it should presumably read '-0.30 < kappa_C < 0.10'.","section":"Fig. 6 caption"},{"comment":"There are two unit/typo issues in the density ranges: '7.36 g/cm2' should be 'g/cm3', and the interval written as '1.12 g/cm3 < rho_lman < 4.23 g/cm3' appears to refer to the upper mantle density rho_uman rather than the lower mantle.","section":"Sec. 2.1, text after Eq. (14)"},{"comment":"The abstract and conclusions contain typographical errors in the quoted numbers: '(-9.3%/+31/7%)' should be '(-9.3%/+31.7%)', and the 'more favorable' case is described with 'zenith angle resolution theta_zres of 10%' where 10 degrees is clearly intended.","section":"Abstract and Sec. 5"},{"comment":"The quantile-weighted probability in Eq. (16) uses a product of independent energy and angle weights w_i w_j, which implicitly assumes that the true energy and true angle distributions are uncorrelated within each bin. Since the public SK release provides only 1D quantiles, this assumption should be stated explicitly as a limitation of the method.","section":"Sec. 3.4, Eq. (16)"}],"recommendation":"major_revision","confidential_remarks":"The paper is an external sensitivity study rather than an HK collaboration result, which is properly acknowledged in the acknowledgments. The main obstacle to publication is the mismatch between the abstract's headline numbers and the constraints derived and partially implemented in the paper itself; this is fixable with a careful revision of the abstract, Section 4, and Table 2. The lack of a nuisance-parameter treatment should also be addressed or explicitly qualified, since otherwise the absolute precision claims are easy to over-interpret."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe one thing to know: this paper does something genuinely useful and then undercuts it with a headline it can't actually claim. It takes the public Super-Kamiokande atmospheric-neutrino release, scales it to Hyper-Kamiokande, and projects sensitivity to the average densities of the core, lower mantle, and upper mantle, while imposing M⊕, I⊕, and hydrostatic-equilibrium constraints. Those numbers are new, and the input is reproducible: the SK prefit release is public, and the constraint coefficients come from a separate parameter-free derivation. The authors also deserve credit for being honest in Table 2, where the 95% upper limits that would violate the constraints are replaced by dashes, and for flagging an error in the earlier Ref. [60].\n\nThe soft spot is the abstract, and it is not cosmetic. The nominal case (Eres=30%, θzres=20°, sin²θ23=0.45) is quoted as determining core density at 2σ with uncertainty (-14.5%)/+39.5%. But Section 2.1 derives κC<0.10 from the M⊕/I⊕/EHE constraints, and the text itself concedes this in Section 4. At κC=0.395, the compensating relations Eq. (10) force the upper mantle density to 6.54 g/cm³ and the lower mantle to 2.13 g/cm³, which is exactly the violation of Eq. (9) the paper says it enforces. So the positive branch of the nominal sensitivity curve is evaluated on Earth models the paper has already ruled out. Table 2 knows this; the abstract and conclusions don't. The actual constrained posterior upper edge is roughly +10%, and that should be the quoted number.\n\nThe second issue is less severe but real. The χ² is computed from expectation values only, with no nuisance parameters for flux normalization or cross sections. The \"nominal\" case is essentially resolution smearing, not a full systematic model. The paper openly notes the gap to the statistical-only case, so the quoted intervals are optimistic, not fraudulent—but the word \"systematics\" is doing more work than the analysis actually performs.\n\nProportionally: the central idea is sound. The favorable-resolution case (E20%, θzres=10°) gives (-8.3%)/+9.8%, which is inside the constraint and looks like a plausible HK capability. The nominal case, once truncated to the constraint, is still a meaningful bound: roughly -14.5%/+10% at 2σ. That is worth knowing.\n\nWho is this for? Anyone planning atmospheric-neutrino Earth tomography with water Cherenkov detectors. It deserves peer review, but the authors should be required to reconcile the abstract with their own constraints. A revised version that quotes the truncated interval and drops the unphysical upper edge is a solid contribution.","headline":"Genuinely useful HK sensitivity projection undermined by an abstract that quotes an interval inconsistent with the paper's own geophysical constraints.","tokens_in":22267,"tokens_out":3207,"would_cite":false,"duration_ms":32640,"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":"Hyper-Kamiokande can measure the Earth's average core density to within about ±10 percent by watching atmospheric neutrinos oscillate through the planet.","keywords":["neutrino oscillation tomography","Earth core density","atmospheric neutrinos","Hyper-Kamiokande","PREM reference model","Earth mass and moment of inertia constraints","matter effects","mantle-core density jump"],"falsifier":"Run the first 6500 days of real Hyper-Kamiokande data through the paper's own 12-sample binning and compute $\\Delta\\chi^2$ as a function of $\\kappa_C$ against the PREM prediction. If the resulting $2\\sigma$ interval for $\\kappa_C$ spans the entire physically allowed range $[-0.30, +0.10]$ without excluding either edge, the claimed core-density determination is refuted. A faster in-situ check: measure HK's actual reconstructed energy and angular resolutions; if they are worse than $E_{\\mathrm{res}} \\approx 20\\%$ and $\\theta_{z\\mathrm{res}} \\approx 10^\\circ$, the favorable-case $(-8.3\\%)/+9.8\\%$ interval is unreachable under the paper's own scaling.","tokens_in":21263,"feed_emoji":"🌍","tokens_out":20727,"duration_ms":190497,"temperature":0.7,"pith_summary":"The paper asks what the future Hyper-Kamiokande detector can learn about the Earth's deep interior from atmospheric neutrinos that pass through the planet. It argues that matter effects on neutrino oscillations leave a record of the average densities of the core, lower mantle, and upper mantle, and that after 6500 days HK can determine the core average density $\\bar{\\rho}_C$ at $2\\sigma$ confidence with uncertainty $(-14.5\\%)/+39.5\\%$ under nominal resolutions, or $(-8.3\\%)/+9.8\\%$ if energy and angular resolution improve as hoped. It further shows that independent geophysical knowledge — the Earth's measured mass $M_\\oplus$, moment of inertia $I_\\oplus$, and hydrostatic equilibrium — caps any positive core-density deviation at about $+10\\%$, so neutrinos and geophysics together bound the core density inside a narrow band. If the projection holds, neutrino oscillations become a non-seismic probe of the planet's layered interior, independently checking the sharp density jump at the mantle-core boundary.","feed_headline":"6500 days of Hyper-Kamiokande data could weigh Earth's core to ±10%","feed_subtitle":"Oscillations of neutrinos crossing the planet could confirm the mantle-core density jump without seismology.","key_machinery":"The central mechanism is the Mikheyev--Smirnov--Wolfenstein matter potential $V = \\sqrt{2}\\,G_F N_e$, which makes neutrino oscillation probabilities depend on the electron density along the trajectory; binning atmospheric neutrinos in energy and zenith angle lets each bin sample a different chord through the Earth, and core-crossing chords ($\\cos\\theta_z < -0.84$) are where the core's imprint is largest, aided by the resonance-like mantle-core interference effect in the 2--10 GeV range. Deviations from the PREM reference profile are parametrized by constant scale factors $\\rho'_i = (1+\\kappa_i)\\rho_i$ for the core, lower mantle, and upper mantle, and the geophysical side conditions are imposed analytically: conserving $M_\\oplus$ and $I_\\oplus$ fixes $\\kappa_{\\mathrm{lman}} = -1.43\\,\\kappa_C$ and $\\kappa_{\\mathrm{uman}} = 2.070\\,\\kappa_C$, while hydrostatic equilibrium restricts $\\kappa_C$ to roughly $[-0.33, +0.09]$. Sensitivities then come from a $\\Delta\\chi^2 = \\sum_i (O_i - E_i)^2/E_i$ comparison between PREM and modified-Earth expected counts, computed from the public Super-Kamiokande pre-fit simulation (12 multi-GeV and partially contained samples, with the collaboration's flux and interaction modeling) scaled to HK's volume and re-binned for each assumed resolution.","core_discovery":"The paper's claim is that Hyper-Kamiokande, a water Cherenkov detector eight times larger than Super-Kamiokande, can after 6500 days of live time determine the Earth's average core density $\\bar{\\rho}_C$ at $2\\sigma$ confidence from atmospheric neutrino oscillations alone, with relative uncertainty $(-14.5\\%)/+39.5\\%$ in the nominal configuration ($E_{\\mathrm{res}} = 30\\%$, $\\theta_{z\\mathrm{res}} = 20^\\circ$, $\\sin^2\\theta_{23} = 0.45$) and $(-8.3\\%)/+9.8\\%$ in the favorable configuration ($E_{\\mathrm{res}} = 20\\%$, $\\theta_{z\\mathrm{res}} = 10^\\circ$, $\\sin^2\\theta_{23} = 0.58$). The authors stress that the positive side of these intervals is not the whole story: any modified density profile must conserve the measured Earth mass $M_\\oplus$ and moment of inertia $I_\\oplus$ and obey hydrostatic equilibrium ($\\bar{\\rho}'_{\\mathrm{uman}} < \\bar{\\rho}'_{\\mathrm{lman}} < \\bar{\\rho}'_C$), which limits positive core-density deviations to about $+10\\%$ — so the $+39.5\\%$ edge of the nominal interval is physically excluded, and the tables mark such cases with constraint-supplied limits. The same analysis yields correlated bounds on the lower and upper mantle average densities, and the authors conclude that with $\\sin^2\\theta_{23}=0.50$ and favorable resolution the $2\\sigma$ core interval of $10.08\\ \\mathrm{g/cm^3} \\le \\bar{\\rho}_C \\le 12.06\\ \\mathrm{g/cm^3}$ around the PREM value $10.99\\ \\mathrm{g/cm^3}$ would constitute independent, non-seismic evidence for the existence of at least three major density layers inside the Earth.","pith_inferences":["The paper's own gap between statistical-only and nominal sensitivity implies the measurement is systematics-limited, so the favorable-case interval will be reached only if HK's real reconstruction and calibration outperform Super-Kamiokande's; a conservative reading places the realized core-density uncertainty between the nominal and favorable columns.","The same three-layer machinery applies directly to other future neutrino detectors (DUNE, ORCA, INO): because the geophysical constraints link the three densities, combining detectors could break the degeneracy between core density, mantle density, and the electron fraction $Y_e$, which this paper holds fixed.","A near-term validation is possible: the favorable-resolution scenario predicts a specific zenith-angle-dependent distortion of up-going multi-GeV rates that grows with core-crossing path length, so early HK data can test the projection long before 6500 days accumulate."],"forward_implications":["A $2\\sigma$ core-density determination near $10.08$--$12.06\\ \\mathrm{g/cm^3}$ (favorable resolution, $\\sin^2\\theta_{23} = 0.50$) would independently confirm at least three major density layers inside the Earth, including a large density jump between mantle and core — a result currently resting almost entirely on seismology.","The true value of $\\sin^2\\theta_{23}$ materially changes the reach: nominal-case uncertainty shifts from $(-14.5\\%)/+39.5\\%$ at $\\sin^2\\theta_{23} = 0.45$ to $(-9.3\\%)/+31.7\\%$ at $0.58$, so pinning down $\\theta_{23}$ sharpens the tomographic measurement.","Improving the reconstructed zenith-angle resolution buys more sensitivity than a comparable improvement in energy resolution, so reconstruction R&D directly translates into geophysical reach.","Because the statistical-only sensitivity is far better than the nominal one, the measurement is dominated by systematic effects; reducing them is what separates the nominal from the favorable-case outcome.","The geophysical constraints cap positive core deviations near $+10\\%$, so any HK preference for $\\kappa_C > 0.10$ would signal a conflict with established Earth-mass and hydrostatic-equilibrium assumptions rather than a new density measurement."],"supporting_citations":[{"why":"The Super-Kamiokande atmospheric oscillation analysis whose simulated event rates, sample selections, and 6500-day normalization supply the baseline for the sensitivity projection.","marker":"[72]"},{"why":"The public Super-Kamiokande simulation release that the paper uses directly as its input event-rate distributions, scaled to Hyper-Kamiokande's volume.","marker":"[74]"},{"why":"The atmospheric neutrino flux model used to generate the simulated rates that the sensitivity calculation relies on.","marker":"[75]"},{"why":"PREM, the reference Earth density model against which all core and mantle density deviations are defined.","marker":"[10]"},{"why":"The derivation of the Earth mass, moment-of-inertia, and hydrostatic-equilibrium constraints that yields the compensation factors $\\kappa_{\\mathrm{lman}} = -1.43\\,\\kappa_C$ and $\\kappa_{\\mathrm{uman}} = 2.070\\,\\kappa_C$.","marker":"[60]"},{"why":"The mantle-core interference (neutrino oscillation length resonance-like) effect that makes core-crossing oscillations sensitive to the mantle-core density jump.","marker":"[61]"},{"why":"The global neutrino oscillation fit that fixes the PMNS parameters and defines the $3\\sigma$ range of $\\sin^2\\theta_{23}$ scanned in the sensitivity study.","marker":"[73]"},{"why":"The neutrino interaction model used inside the Super-Kamiokande simulation that produces the event rates the analysis reweights.","marker":"[76]"}],"fun_headline_variants":["Neutrino oscillations to weigh Earth's core within 10%","HK's neutrino data could map Earth's core in 6500 days","Neutrino tomography of Earth: HK to measure core density","Atmospheric neutrinos could weigh Earth's core to 10%","6500 days of neutrino data: HK to probe Earth's core"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The projections assume that Super-Kamiokande's public pre-fit simulation, scaled to Hyper-Kamiokande's eightfold larger volume and re-smeared with assumed Gaussian resolutions, faithfully represents how HK will record 6500 days of atmospheric neutrinos, and that systematic uncertainties such as flux normalization and cross-section errors will not degrade the quoted precision; the paper's own large gap between its statistical-only and nominal benchmarks shows that real systematics decide which column of numbers comes true.","fun_headline_variants_meta":{"raw":{"variants":["Neutrino oscillations to weigh Earth's core within 10%","HK's neutrino data could map Earth's core in 6500 days","Neutrino tomography of Earth: HK to measure core density","Atmospheric neutrinos could weigh Earth's core to 10%","6500 days of neutrino data: HK to probe Earth's core"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001453,"raw_usage":{"total_tokens":6092,"prompt_tokens":1426,"completion_tokens":4666,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":1042,"completion_tokens_details":{"reasoning_tokens":4575}},"tokens_in":1042,"tokens_out":4666,"duration_ms":35929,"temperature":1.0,"reasoning_tokens":4575,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T17:38:00.451814+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the first 6500 days of real Hyper-Kamiokande data through the paper's own 12-sample binning and compute $\\Delta\\chi^2$ as a function of $\\kappa_C$ against the PREM prediction. If the resulting $2\\sigma$ interval for $\\kappa_C$ spans the entire physically allowed range $[-0.30, +0.10]$ without excluding either edge, the claimed core-density determination is refuted. A faster in-situ check: measure HK's actual reconstructed energy and angular resolutions; if they are worse than $E_{\\mathrm{res}} \\approx 20\\%$ and $\\theta_{z\\mathrm{res}} \\approx 10^\\circ$, the favorable-case $(-8.3\\%)/+9.8\\%$ interval is unreachable under the paper's own scaling.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"PREM, the reference Earth density model against which all core and mantle density deviations are defined."},{"cited_title":"Neutrino Tomography of the Earth: the Earth Total Mass, Moment of Inertia and Hydrostatic Equilibrium Constraints","cited_arxiv_id":"2406.13727","evidence_quote":"The derivation of the Earth mass, moment-of-inertia, and hydrostatic-equilibrium constraints that yields the compensation factors $\\kappa_{\\mathrm{lman}} = -1.43\\,\\kappa_C$ and $\\kappa_{\\mathrm{uman}} = 2.070\\,\\kappa_C$."},{"cited_title":"Hayato, ”A neutrino interaction simulation program library NEUT,” Acta Phys","cited_arxiv_id":null,"evidence_quote":"The neutrino interaction model used inside the Super-Kamiokande simulation that produces the event rates the analysis reweights."}],"review_version":1}