{"id":"0ca803fb-d6da-494c-8fea-58fd30d1051e","arxiv_id":"2501.07621","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":2.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Simulations show ICAL could constrain Earth's core-mantle density jump to about 16-18% precision and simultaneously locate the core radius.","lead":"This paper forecasts how well the proposed ICAL neutrino detector could measure Earth's core-mantle boundary radius and density jumps using atmospheric neutrinos. It finds that with one megaton-year of data, ICAL could pin the core-mantle density jump to roughly 16-18% precision.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The forecast's five-layer uniform-density model and unstated mass-to-electron density conversion (Z/A) could bias the quoted constraints; a PREM-closure test is needed.","rationale":"The reader's weakest assumption correctly identifies the five-layer model and the unvaried Z/A conversion as the key fragility. My stress-test focuses on the same link: neutrino oscillations are sensitive to electron density, yet the paper's parameter space and constraints are defined in terms of mass density, with no explicit mapping or uncertainty on composition. The absence of a closure test against the real PREM profile means the quoted 1σ intervals are conditional on the model being a faithful proxy. This is load-bearing because even a perfect detector cannot correct a systematic offset from using the wrong Earth model. I agree with the reader's CONDITIONAL verdict; the central claim is plausible but not fully supported by the standalone proceedings text. The concrete test would settle whether the model misspecification actually shifts the forecasts enough to matter. Other concerns (e.g., omitted systematics) are secondary because they would degrade precision, while this concern could bias the central values themselves.","tokens_in":4920,"tokens_out":7260,"duration_ms":76179,"concrete_test":"Simulate 1 Mt·yr ICAL data with the full PREM electron density profile (using layer-dependent Z/A: ~0.467 in the core, ~0.495 in the mantle) and the same detector response, flux, and binning as in Ref. [6]; then fit these data using the paper's five-layer model as the theory. If the recovered best-fit (∆ρ_CMB, R_CMB) shifts by more than the 1σ contour width in Fig. 3, or if the effective PREM-equivalent values fall outside the 1σ contours, the forecasted constraints are biased by model misspecification and the quoted precision is not robust to the Earth model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central sensitivity claim in §III rests on the five-layer density model introduced in §II, but the paper never specifies how mass density ρ is converted to the electron density n_e that enters the oscillation Hamiltonian. If a single global Z/A (e.g., 0.5) is used, the core (Z/A ≈ 0.467) and mantle (≈ 0.495) are mis-modeled; the resulting ~6% difference in n_e at fixed ρ is comparable to the quoted 16% precision and can shift the inferred ∆ρ_CMB and R_CMB. Furthermore, the simulated data are generated with the same five-layer model as the fit (a closure test), so the quoted [5.1, 7.0] g/cm³ interval is a self-fit that excludes model misspecification. The 'hydrostatic equilibrium condition' in §II is enforced only as monotonic density (ρ_inner > ρ_outer), not as full pressure balance; this weak prior could either admit or exclude profiles that a more physical Earth model would rule out, though the effect on the contours is not tested. If the real Earth's electron density profile differs from the five-layer uniform proxy, the best-fit parameters and the coverage of the 1σ contours could be biased, undermining the claim that ICAL would determine the CMB density jump and radius.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This proceedings contribution reports a sensitivity forecast for the proposed 50 kt INO-ICAL detector to constrain Earth's internal density structure using atmospheric neutrino oscillations. The authors adopt a five-layer uniform-density model of Earth, with the core-mantle boundary radius R_CMB and the density jump at that boundary as two free parameters, while fixing total mass, moment of inertia, outer-mantle density, other layer radii, and the inner-core/outer-core density ratio. Pseudo-data are generated from the standard five-layer profile, and fits are performed to modified profiles. The main quantitative results are that with 1 Mt·yr exposure and charge identification, the 1σ allowed range for the CMB density jump is [5.1, 7.0] g/cm³ for normal mass ordering, corresponding to about 16% precision, and that two-dimensional contours can simultaneously constrain R_CMB and the density jumps at the IC-OC, CMB, IM-MM, and MM-OM boundaries. The paper concludes that an ICAL-like detector would provide information complementary to seismic and gravitational probes of Earth's interior.","tokens_in":5205,"tokens_out":3669,"duration_ms":40671,"significance":"If the forecast is robust, the result is significant: it would establish that a weak-interaction probe can independently measure the density jump and location of the core-mantle boundary, complementing PREM-based information. The analysis uses a standard simulation-and-fit procedure, and the emphasis on the role of charge identification is a useful, explicit demonstration. However, the paper is a compact proceedings summary that delegates most technical detail to Ref. [5], and the quoted precision is derived from a self-consistent closure test in which the same five-layer model generates the pseudo-data and is used in the fit. As a forecast, the [5.1, 7.0] g/cm³ interval is internally credible, but the manuscript does not currently establish that the interval survives realistic model-misspecification and systematic effects. The central claim is therefore plausible but not yet fully supported in this standalone contribution.","major_comments":[{"comment":"The manuscript never specifies how the mass densities ρ_i of the five-layer model are converted to the electron density n_e that enters the matter potential in the oscillation Hamiltonian. A global conversion with a single Z/A value (for example 0.5) would mis-model the core (Z/A ≈ 0.467) relative to the mantle (Z/A ≈ 0.495). The resulting difference of roughly 6% in n_e at fixed mass density is comparable to the quoted 16% precision on Δρ_CMB, so this omission is potentially load-bearing for the [5.1, 7.0] g/cm³ interval. The authors should state the conversion used, and preferably show that a PREM-based electron-density profile or a layer-dependent Z/A prescription does not shift the contours by more than a small fraction of the quoted precision.","section":"§II and §III"},{"comment":"The 'hydrostatic equilibrium condition' is imposed only as monotonic density, ρ_inner > ρ_outer, and not as a full pressure-balance or equation-of-state condition. This weak prior defines the gray 'unphysical' regions in Figs. 2 and 3, but the paper does not test how much of the allowed parameter space is actually a consequence of this weak condition rather than of the neutrino data. A more realistic hydrostatic constraint (e.g., requiring approximate pressure balance or consistency with a known equation of state) could shrink or shift the allowed region. The authors should either justify the monotonic-density condition as a deliberate conservative choice or show that the main contours are insensitive to strengthening it.","section":"§II"},{"comment":"The predicted intervals and contours come from fitting pseudo-data that were generated from the same five-layer model used in the fit. This is a valid closure test for an idealized sensitivity, but it excludes model misspecification. The paper should add a PREM-closure test: generate pseudo-data from the full continuous PREM profile (or from a four- or six-layer model with different transition radii) and fit with the five-layer model, to estimate the bias in R_CMB and Δρ_CMB that would arise if the real electron-density profile deviates from the assumed piecewise-uniform shape. Relatedly, no systematic uncertainties (flux normalization, cross sections, detector energy and angular resolution, charge-ID efficiency) are included; even a brief statement of the dominant systematics and their estimated impact on the quoted precision is needed for a standalone claim about what ICAL 'would constrain.'","section":"§III"}],"minor_comments":[{"comment":"The label 'Radial Dis' in the figure appears to be an incomplete phrase; it should read 'Radial distance' or similar.","section":"Fig. 1 caption"},{"comment":"The gray regions are labeled 'Hydro. Cond. Violated' in a way that could be compressed or clarified; consider one unified legend entry explaining that these regions violate the monotonic-density condition described in §II.","section":"Fig. 2 and Fig. 3"},{"comment":"The abstract and body say 'demonstrate how well ... would constrain,' but the analysis is explicitly a median-sensitivity forecast over simulated data. The wording should be made consistently prospective to avoid implying that real data have been analyzed.","section":"Introduction and Conclusions"},{"comment":"The paper relies heavily on Ref. [5] for the definition of Δχ² and for the underlying scenarios, and on Ref. [4] for the detector simulation. The text should state explicitly which results are new to this contribution and which are reproduced from [5], especially since the figures are already taken from there.","section":"References"},{"comment":"The sentence 'Without the CID capability, the sensitivity gets reduced significantly if the true neutrino mass ordering is IO' is vague; a quantitative statement (e.g., the 1σ interval or its width in the IO, no-CID case) would be more informative.","section":"§III"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is essentially a proceedings-style summary of Ref. [5], and most of the technical content is delegated to that reference. For a full journal attachment, the referee's major comments are appropriate requests for completeness rather than fundamental objections. The electron-density conversion and the hydrostatic-condition issue are particularly important because they affect whether the quoted [5.1, 7.0] g/cm³ interval is a statement about the real Earth or a statement about the five-layer toy model. I recommend major revision rather than rejection because the underlying forecasting methodology is standard and the missing items can be supplied within the scope of the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this is a NuFact proceedings writeup whose three figures and numerical results are explicitly taken from the authors' earlier arXiv paper (Ref [5]). As a standalone contribution it adds no new physics, no new data, and no new analysis. Second, the sensitivity forecast itself is standard closure-test forecasting: pseudo-data are generated from a five-layer Earth model and then fit with modified profiles. That is exactly how one should estimate an experiment's reach, and the numbers are credible.\n\nWhat the paper does well: it clearly describes the five-layer parameterization and the eight constraints (fixed mass, moment of inertia, radii, density ratios) that reduce it to two free parameters. The point that charge identification matters is made with a clean comparison of w/ and w/o CID curves. For a reader who wants the bottom line from Ref [5] without wading through the full paper, this is a serviceable summary.\n\nSoft spots. All three figures are reproductions, so the only new text is descriptive. The paper does not specify how mass density ρ is converted to the electron density n_e that appears in the oscillation Hamiltonian. If the conversion uses a single global Z/A, the core-mantle difference (about 0.467 vs 0.495) is ignored, a ~6% effect on n_e at fixed ρ. That is roughly a third of the quoted 16–18% precision and could shift the inferred parameters. Because the same conversion is used in generating the pseudo-data and in the fit, the closure test does not expose this; it is a modeling bias, not a statistical one. The hydrostatic equilibrium condition is also enforced only as monotonic density, not full pressure balance. These are real caveats, but they are minor for a forecast: they do not invalidate the quoted sensitivities, they just make the forecast conditional on the five-layer electron-density ansatz.\n\nWho this is for: people tracking neutrino-geophysics forecasts, especially for ICAL/INO. It deserves a serious referee only as a proceedings note; for a journal it is a duplicate of Ref [5]. I would not cite this paper over the original, but I would keep Ref [5] in view. If this came to a journal as a regular article, I would recommend desk rejection on grounds of lack of novelty; for a proceedings volume, it is acceptable as a summary.","headline":"A clean, thin proceedings summary of a prior sensitivity study; the forecast is credible but adds no new results and leaves the electron-density conversion implicit.","tokens_in":5713,"tokens_out":3176,"would_cite":false,"duration_ms":33952,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Simulated ICAL data show that atmospheric neutrinos would measure Earth's core-mantle density jump to about 16% precision at 1σ.","keywords":["atmospheric neutrinos","Earth matter effects","core-mantle boundary","density jumps","neutrino oscillations","ICAL detector","charge identification","Earth interior structure"],"falsifier":"Re-run the same simulation with the full continuous PREM density profile, and additionally allow the electron-density conversion ($Z/A$) to vary with depth, as the true model; if the reconstructed $R_{\\rm CMB}$ and $\\Delta\\rho_{\\rm CMB}$ shift by more than the quoted $1\\sigma$ widths, the five-layer uniform-density simplification is responsible for the reported constraints.","tokens_in":4730,"feed_emoji":"🌍","tokens_out":5458,"duration_ms":50902,"temperature":0.7,"pith_summary":"Atmospheric neutrinos that cross Earth feel matter effects that depend on the electron density along their path, so their oscillation patterns carry information about Earth's interior. This paper asks how well a future iron calorimeter detector (ICAL), with 1 Mt·yr of data, could measure the density jump at the core-mantle boundary and the boundary's radius. The paper works with a five-layer Earth model in which total mass, moment of inertia, and most radii are held fixed; only the core-mantle radius and a few density jumps are free. Its central forecast is that ICAL would constrain the core-mantle density jump to roughly 16–18% precision, and that charge identification is needed to achieve this. If true, neutrino oscillations become an independent, weak-interaction-based probe of Earth's deep structure, complementary to seismology and gravimetry.","feed_headline":"Neutrino oscillations could measure the core-mantle density jump to 16%","feed_subtitle":"A 1 Mt·yr ICAL exposure would pin down the core-mantle boundary and its density jump—if the five-layer Earth model holds.","key_machinery":"The central mechanism is Earth-matter-induced MSW neutrino oscillations, quantified by comparing simulated ICAL data against modified five-layer density profiles through $Δ\\chi^2$ statistics. The two free parameters are the core-mantle radius $R_{\\rm CMB}$ and the density jump at the core-mantle boundary $\\Delta\\rho_{\\rm CMB}$, while the other eight layer parameters are tied by fixed total mass, fixed moment of inertia, fixed outer-mantle density, fixed radii for most layers, and the PREM-based density ratio of the inner to outer core. Because the oscillation probability depends on the electron density profile along the neutrino's chord through Earth, moving the boundary or changing the jump alters the resonance crossing pattern, which is what makes the reconstructed parameters sensitive to the density structure.","core_discovery":"The paper claims that simulated data from a 50 kt ICAL detector operating for 20 years, equivalent to 1 Mt·yr exposure, would yield a median sensitivity such that the density jump at the core-mantle boundary is constrained to [5.1, 7.0] g/cm³ at 1σ when the true mass ordering is normal and charge identification is used; with inverted ordering the interval is [4.9, 7.0] g/cm³. In the two-dimensional analysis, the same data produce 1σ (2 d.o.f.) contours in the planes of each density jump versus the core-mantle radius, showing that the CMB location and all four boundary jumps can be constrained simultaneously. The charge identification capability of the detector is shown to be crucial, especially for inverted ordering, where without it the sensitivity is significantly reduced. The paper frames this as an independent cross-check on the Preliminary Reference Earth Model, obtained through weak interactions rather than through seismic or gravitational measurements.","pith_inferences":["The paper fixes the composition through a single electron-density conversion, so a testable extension would let the average $Z/A$ vary by layer; the forecasted constraints would likely widen or shift, especially if core composition differs from the PREM assumption.","A real measurement would combine neutrino-derived density constraints with seismic and gravitational priors, and the paper's orthogonal sensitivity suggests a joint inversion could break degeneracies that each probe alone leaves open.","The same five-layer formalism could be applied to other long-baseline atmospheric detectors or neutrino telescopes, where longer baselines through the core would amplify the same matter-effect resonances and sharpen the constraints."],"forward_implications":["With 1 Mt·yr exposure and charge identification, ICAL would bound the core-mantle density jump to $[5.1, 7.0]\\,\\mathrm{g/cm^3}$ at $1\\sigma$ under normal mass ordering, about 16% precision.","Without charge identification, the $1\\sigma$ precision degrades to about 20% for normal ordering and much more severely for inverted ordering, so particle/antiparticle discrimination is load-bearing for the measurement.","Two-dimensional contours in the density-jump versus $R_{\\rm CMB}$ planes show that neutrino data could simultaneously shrink the allowed ranges for all four boundary jumps and the core-mantle radius.","An inverted true mass ordering gives a weaker $1\\sigma$ bound of $[4.9, 7.0]\\,\\mathrm{g/cm^3}$, about 18% precision.","The same framework can be applied to the other density jumps (inner-core/outer-core, inner-mantle/middle-mantle, and middle-mantle/outer-mantle) because fixing Earth's mass and moment of inertia links those jumps to the CMB parameters."],"supporting_citations":[{"why":"Provides the PREM density profile that defines the standard Earth model and the baseline values for layer densities and the core-mantle boundary.","marker":"[1]"},{"why":"Introduces the matter effect in neutrino propagation through Earth, the physical mechanism the analysis relies on.","marker":"[2]"},{"why":"Develops the MSW resonance formalism that makes oscillation probabilities sensitive to density jumps along the neutrino path.","marker":"[3]"},{"why":"Supplies the ICAL detector simulation and event-response details used to generate the prospective data.","marker":"[4]"},{"why":"This is the parent study whose methods, density profiles, and definitions of $\\Delta\\chi^2_{1\\mathrm{D-DJ}}$ and $\\Delta\\chi^2_{\\mathrm{DJ-CMB}}$ are followed here; all quantitative results in this contribution come from it.","marker":"[5]"},{"why":"Provides the event binning scheme and simulation setup used to compute the expected median sensitivities.","marker":"[6]"}],"fun_headline_variants":["Neutrino oscillations could map Earth's core-mantle boundary","Atmospheric neutrinos pinpoint density jumps inside Earth","Neutrino charge ID key to Earth interior tomography","Simulated neutrinos constrain Earth's core radius and jumps","Atmospheric neutrinos as a new probe of Earth's interior"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The forecast assumes that a five-layer uniform-density Earth, with fixed total mass and moment of inertia and PREM-based density ratios, captures every density variation that affects neutrino oscillations.","fun_headline_variants_meta":{"raw":{"variants":["Neutrino oscillations could map Earth's core-mantle boundary","Atmospheric neutrinos pinpoint density jumps inside Earth","Neutrino charge ID key to Earth interior tomography","Simulated neutrinos constrain Earth's core radius and jumps","Atmospheric neutrinos as a new probe of Earth's interior"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000493,"raw_usage":{"total_tokens":2416,"prompt_tokens":931,"completion_tokens":1485,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":547,"completion_tokens_details":{"reasoning_tokens":1404}},"tokens_in":547,"tokens_out":1485,"duration_ms":12304,"temperature":1.0,"reasoning_tokens":1404,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:38:50.225477+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the same simulation with the full continuous PREM density profile, and additionally allow the electron-density conversion ($Z/A$) to vary with depth, as the true model; if the reconstructed $R_{\\rm CMB}$ and $\\Delta\\rho_{\\rm CMB}$ shift by more than the quoted $1\\sigma$ widths, the five-layer uniform-density simplification is responsible for the reported constraints.","supporting_citations":[{"cited_title":"Wolfenstein, Phys","cited_arxiv_id":null,"evidence_quote":"Introduces the matter effect in neutrino propagation through Earth, the physical mechanism the analysis relies on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Develops the MSW resonance formalism that makes oscillation probabilities sensitive to density jumps along the neutrino path."}],"review_version":1}