REVIEW 3 major objections 5 minor 6 references
Exploring constraints on the core radius and density jumps inside Earth using atmospheric neutrino oscillations
T0 review · 3 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read Simulated ICAL data show that atmospheric neutrinos would measure Earth's core-mantle density jump to about 16% precision at 1σ.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [§II and §III] 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.
- [§II] 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.
- [§III] 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.'
minor comments (5)
- [Fig. 1 caption] The label 'Radial Dis' in the figure appears to be an incomplete phrase; it should read 'Radial distance' or similar.
- [Fig. 2 and Fig. 3] 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.
- [Introduction and Conclusions] 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.
- [References] 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.
- [§III] 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.
Circularity Check
No significant circularity: the paper is an explicit sensitivity forecast using simulated data and the same five-layer model for generation and fitting, so its constraints are closure results rather than disguised inputs.
full rationale
The paper does not present a derivation of Earth structure from real data; it presents an expected median sensitivity study. It explicitly states that the prospective data are simulated assuming the standard five-layered density profile and then fitted with modified density jumps, so the quoted intervals such as [5.1, 7.0] g/cm3 are statistical forecasts for how well ICAL could recover assumed input parameters, not predictions derived from independent measurements. The two free parameters (R_CMB and Δρ_CMB) are constrained by the simulated oscillation data, not fixed by the model constraints; the constraints only reduce the parameter space from ten to two. The use of the same model for pseudo-data and fit is a standard closure test for sensitivity, and it does not make the result equivalent to the input by construction. The paper relies on the authors' prior works (Refs. [5] and [6]) for simulation details, binning, and the definition of Δχ2, but this is normal self-citation of the underlying methodology rather than a load-bearing circular premise. No equation is shown to reduce a predicted quantity to a fitted input, and no uniqueness theorem or ansatz is smuggled in via self-citation. The 'hydrostatic equilibrium condition' is only implemented as radial monotonicity of layer densities, which is a modeling approximation and a possible correctness concern, but it is not circularity. Overall, no specific circular step can be quoted, so the score is 0.
Assumptions & free parameters
free parameters (2)
- R_CMB (core-mantle boundary radius) =
scanned 2500-4500 km; standard value 3480 km
- Delta_rho_CMB (density jump at the core-mantle boundary) =
standard value 6.0 g/cm3; constrained to [5.1, 7.0] g/cm3 at 1 sigma (NO, with CID)
assumptions (4)
- domain assumption The five-layered uniform-density Earth model adequately represents the electron density profile encountered by neutrinos.
- standard math Neutrino oscillation probabilities in matter are governed by the standard MSW effect depending on electron density.
- domain assumption Earth's total mass and moment of inertia are known and fixed to PREM values, and the inner-core/outer-core density ratio equals the PREM ratio.
- ad hoc to paper The 'hydrostatic equilibrium condition' is equivalent to monotonic density (rho_inner > rho_outer).
Cite this review
Pith. "Pith review of Exploring constraints on the core radius and density jumps inside Earth using atmospheric neutrino oscillations." pith.science (2026). https://pith.science/paper/CAPB6OPO
@misc{pith2026250107621,
author = {Pith},
title = {Pith review of: Exploring constraints on the core radius and density jumps inside Earth using atmospheric neutrino oscillations},
year = {2026},
howpublished = {\url{https://pith.science/paper/CAPB6OPO}},
note = {Machine review of arXiv:2501.07621}
}
read the original abstract
Atmospheric neutrinos, through their weak interactions, can serve as an independent tool for exploring the internal structure of Earth. The information obtained would be complementary to that provided by seismic and gravitational measurements. The Earth matter effects in neutrino oscillations depend upon the energy of neutrinos and the electron density distribution that they encounter during their journey through Earth, and hence, can be used to probe the inner structure of Earth. In this contribution, we demonstrate how well an atmospheric neutrino experiment, such as an iron calorimeter detector (ICAL), would simultaneously constrain the density jumps inside Earth and determine the location of the core-mantle boundary. In this work, we employ a five-layered density model of Earth, where the layer densities and core radius are modified to explore the parameter space, ensuring that the mass and moment of inertia of Earth remain constant while satisfying the hydrostatic equilibrium condition. We further demonstrate that the charge identification capability of an ICAL-like detector would play a crucial role in obtaining these correlated constraints.
Figures
Reference graph
Works this paper leans on
-
[5]
A. K. Upadhyay, A. Kumar, S. K. Agarwalla, and A. Dighe, (2024), arXiv:2405.04986 [hep-ph]
arXiv 2024
-
[1]
A. M. Dziewonski and D. L. Anderson, Phys. Earth Planet. Interiors 25, 297 (1981)
1981
- [2]
-
[3]
S. P. Mikheev and A. Y. Smirnov, Sov. J. Nucl. Phys. 42, 913 (1985), [Yad.Fiz.42:1441-1448,1985]
work page 1985
-
[4]
S. Ahmed et al. (ICAL), Pramana 88, 79 (2017), arXiv:1505.07380 [physics.ins-det]
arXiv 2017
-
[6]
A. K. Upadhyay, A. Kumar, S. K. Agarwalla, and A. Dighe, JHEP 04, 068 (2023), arXiv:2211.08688 [hep- ph]
arXiv 2023
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.