Quantifying the Information Gain from Future High-Precision Radius Measurements for Identifying Twin Neutron Stars

Bao-An Li; Xavier Grundler

REVIEW 3 major objections 3 minor 1 cited by

Quantifying the Information Gain from Future High-Precision Radius Measurements for Identifying Twin Neutron Stars

T0 review · 3 major / 3 minor · reviewed 2026-08-01 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper claims that measuring neutron-star radii to about 0.2 km precision extracts most of the information available for identifying twin neutron stars, with further improvements yielding diminishing returns.

desk verdict The 0.2-km saturation benchmark is contradicted by the paper's own entropy and distinguishability numbers, but the framework is a useful step toward planning radius-precision requirements. read the letter →

arxiv 2607.18124 v2 pith:EKHESDH7 submitted 2026-07-20 astro-ph.HE astro-ph.GAastro-ph.SRnucl-exnucl-th

Bao-An Li , Xavier Grundler This is my paper · ORCID

classification astro-ph.HEastro-ph.GAastro-ph.SRnucl-exnucl-th

keywords twinneutronstarsradiusmeasurementsBayesianinferenceinformationgainfirst-orderphasetransitionmass-radiusrelationShannonentropyobservationalprecision

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper asks how precise neutron-star radius measurements must be to identify twin stars—pairs of stars with the same mass but different radii, a signature of a strong first-order phase transition inside supradense matter. Using mock radius data with uncertainties from 0.9 km down to 0.1 km, the authors run Bayesian inference and find three regimes: above 0.6 km the posterior is prior-dominated, between 0.2 and 0.6 km information is gained most rapidly, and below 0.2 km the posterior saturates. They conclude that about 0.2 km precision is sufficient to extract most of the available information, a quantitative benchmark for designing future X-ray and gravitational-wave missions.

What carries the argument

The central object is the ratio ΔR/(2σ_R), the intrinsic twin-star radius separation divided by twice the observational uncertainty. The analytical distinguishability probability P_dis = Φ(ΔR/(2σ_R))—the Gaussian cumulative probability that an observed radius is assigned to the correct branch—governs branch resolvability and saturates near unity when the ratio exceeds about 2. The Bayesian inference uses a nine-parameter meta-model equation of state with a uniform prior on the hadron-quark transition density over (1–6)ρ₀, and the posterior distribution of ΔR is combined with P_dis into an observational efficiency integral and the Shannon entropy of the posterior.

What would settle it

A concrete test: re-run the identical Bayesian inference using a different equation-of-state meta-model (for example, a speed-of-sound parameterization) or a wider prior on the transition density, and check whether the posterior mean of ΔR still flattens at σ_R ≈ 0.2 km; a continued steep decline below 0.2 km would falsify the claimed saturation threshold as a general result.

Watch

Extended reading notes

Core claim

Within the paper's Bayesian framework, the posterior distribution of the maximum twin-star radius separation ΔR narrows most rapidly as the assumed radius uncertainty σ_R decreases from 0.6 km to 0.2 km, and essentially stops narrowing below 0.2 km. This three-regime structure—prior-dominated, rapid-gain, and information-saturation—is reproduced by an analytical model where the probability of statistically resolving the two mass-radius branches is Φ(ΔR/(2σ_R)), and by complementary measures based on branch observational efficiency and Shannon entropy. The paper therefore claims that radius measurements with precision around 0.2 km are sufficient to identify twin neutron stars, and that furth

Load-bearing premise

The 0.2-km benchmark is computed inside the paper's own Bayesian framework—a nine-parameter meta-model equation of state and a uniform prior on the hadron-quark transition density over (1–6)ρ₀—so if the true equation of state or prior lies outside this family, the saturation threshold could shift.

Editorial extensions

If this is right

If the benchmark holds, future missions aiming for ~0.2 km radius precision will capture most of the twin-star information; designs targeting 0.1 km add little.
The three-regime decomposition gives a rule of thumb for when a measurement campaign transitions from prior-dominated to likelihood-dominated.
The analytical distinguishability formula can be used to quickly estimate twin-star observability from any predicted ΔR and assumed σ_R.
The framework extends to other questions where progressively precise observations constrain a parameter, identifying the point of diminishing returns.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

Editorial inference: the 0.2-km threshold is likely sensitive to the chosen equation-of-state parameterization and the prior on the transition density; a different generative model could shift the saturation point, so the benchmark should be re-checked with alternative meta-models.
Editorial inference: the same distinguishability logic applies to other binary-branch problems in astrophysics, such as separating hadronic and hybrid branches in mass-tidal-deformability space, where the dimensionless ratio of separation to uncertainty plays the same role.
Editorial inference: a testable extension is to run the same analysis with a prior that includes transition densities above 6ρ₀ or with a speed-of-sound parameterization, and check whether the 0.2-km saturation persists.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Dear X,

The headline is that this paper is a clearly written extension of the authors' Bayesian twin-star program, but its central 0.2-km saturation benchmark is contradicted by its own numbers. The three-regime description of the posterior mean is a reasonable read of their Fig. 3, but the claim that 'most of the information is acquired before 0.2 km' does not survive the entropy and distinguishability curves in the same paper.

What is actually new: they apply Shannon entropy and an efficiency integral to the twin-star identification problem and give a quantitative threshold (0.2 km) for diminishing returns, with a three-regime decomposition (prior-dominated, rapid gain, saturation). No one has stated this benchmark before, so the question they address is timely for mission planning. The paper is also honest in scope: it says explicitly that the result is within the adopted Bayesian framework and EOS meta-model.

What it does well: the analytical branch-distinguishability formula (Eq. 1) is textbook but correctly applied, and they correctly distinguish the two roles of σ_R: selecting EOSs in the likelihood versus setting the resolution scale for branch separation. The posterior plots and table indicate a real computational effort.

The soft spots are substantial. First, the internal contradiction: Section 2.3 says the Shannon entropy 'increases most rapidly between σ_R=0.1 and 0.2 km.' Since entropy measures residual uncertainty, a steep rise over that interval means the biggest information gain per km of precision happens exactly between 0.1 and 0.2 km. Similarly, Eq. (1) gives P_dis ≈ 0.89 at ΔR=0.5 km, σ=0.2 km versus 0.994 at 0.1 km; the residual misclassification rate drops by more than an order of magnitude. That is not saturation. The authors do hedge by saying the different functionals 'need not display identical features,' but then they still claim the analyses 'consistently indicate' the 0.2-km benchmark. That step is not supported.

Second, the entire quantitative result is computed inside the authors' own posterior from their earlier meta-model EOS and prior on transition density (1–6)ρ_0. That is legitimate as a framework-dependent statement, but the saturation threshold is a property of that prior, not of nature. The paper itself concedes the threshold 'may depend somewhat' on the adopted EOS parameterization. Third, no code, data, or MCMC diagnostics are provided, and the run lengths differ across σ_R values (300,000 steps for 0.9 and 0.5 km, 600,000 for 0.2 and 0.1 km), which could affect the comparison.

Who is this for: someone planning future radius-measurement campaigns wants to know whether pushing from 0.2 to 0.1 km pays off. The paper gives a useful framework for asking that question, but the answer it advertises is likely too pessimistic — 0.2 km is not the saturation point even on their own metrics. A revised version moving the benchmark to ~0.1 km or restricting the claim to the posterior mean of ΔR would be much harder to attack.

I would send it to a referee: the question is important, the framework is reusable, and the internal contradiction is fixable. My own verdict would be major revision rather than acceptance.

Referee Report

3 major / 3 minor

Summary. This Letter uses Bayesian inference with mock radius measurements of a canonical 1.4 solar-mass neutron star to quantify how the precision of future radius observations would improve our ability to identify twin neutron stars. The radius uncertainty is varied from 0.9 km to 0.1 km. The paper analyzes the posterior distribution of the maximum radius separation ΔR, an analytical Gaussian branch-distinguishability probability, a branch observational efficiency, and the Shannon entropy of the posterior. It identifies three regimes: a prior-dominated regime for σ_R ≳ 0.6 km, a rapid information-gain regime for 0.2 ≲ σ_R ≲ 0.6 km, and an information-saturation regime for σ_R ≲ 0.2 km, concluding that 0.2-km precision already extracts most of the available information for identifying twin NSs within the adopted Bayesian framework and EOS meta-model.

Significance. If the 0.2-km saturation claim were sound, the paper would provide a useful quantitative benchmark for planning next-generation X-ray timing and gravitational-wave observatories, and the proposed framework for measuring information gain as a function of observational precision would be of general methodological value. The paper is clearly written, uses a published Bayesian EOS meta-model, and supplies an analytical distinguishability formula (Eq. 1) that is independently checkable. However, the central quantitative claim is contradicted by the paper's own reported information-theoretic and distinguishability metrics, so the significance of the benchmark is currently not established.

major comments (3)

[Section 2.3, Fig. 5, and Eq. (1)] The paper's own metrics contradict the claimed 'information-saturation regime for σ_R ≲ 0.2 km'. Section 2.3 explicitly states that the Shannon entropy 'increases most rapidly between σ_R = 0.1 and 0.2 km', meaning that the entropy drop—and hence the information gain—when improving from 0.2 km to 0.1 km is the largest in the plotted range. This is the opposite of saturation. Similarly, Eq. (1) gives P_dis = Φ(1.25) ≈ 0.89 for ΔR = 0.5 km at σ_R = 0.2 km, but P_dis = Φ(2.5) ≈ 0.994 at σ_R = 0.1 km; the remaining misclassification probability drops by more than an order of magnitude. Thus the abstract's 'most of the information ... is acquired before ... 0.2 km' is not supported by the paper's own quantitative measures.
[Section 2 opening and Summary] The 0.2-km threshold is computed within a single EOS meta-model (Zhang et al. 2018; Xie & Li 2020, 2021) with a uniform prior on the hadron–quark transition density over (1–6)ρ_0. The paper acknowledges this caveat ('within the present Bayesian framework'), but then generalizes: 'the existence of an information-saturation regime is expected to be a generic consequence of Bayesian inference.' No robustness test is provided—e.g., varying the meta-model parameterization, the transition-density prior, or the mock central radius. Because the central claim is intended as an observational benchmark, the model-dependence of the threshold should be quantified or the generalization softened.
[Section 2.1] The claim that 'most of the posterior narrowing occurs when the observational precision increases from σ_R = 0.9 km to about 0.2 km, whereas further reductions to 0.1 km produce only modest changes' is not supported by any plotted or tabulated measure of posterior width or variance. Moreover, the Shannon entropy behavior reported in Section 2.3 suggests that the posterior continues to narrow substantially below 0.2 km. If the authors intend a narrower meaning of 'narrowing' (e.g., of the posterior mean only), this should be stated explicitly and quantified; otherwise the statement appears inconsistent with the paper's own entropy analysis.

minor comments (3)

[Abstract and Section 3] The phrase 'These complementary analyses consistently indicate' is inaccurate given the entropy behavior described in Section 2.3; at minimum, the abstract should be reworded to reflect the actual consistency (or lack thereof) among the adopted information measures.
[Section 2.2] The text says that for ΔR/(2σ_R) ≳ 2 the branches are 'nearly completely' resolved, and that P_dis = 0.89 at σ_R = 0.2 km is 'close to the saturation region'. Since the saturation region is defined by ratios ≳2, a ratio of 1.25 (P_dis = 0.89) is not close to saturation under the paper's own definition. This wording should be revised.
[Table 1] The fraction of twin-producing EOSs among accepted EOSs is roughly constant across σ_R (14.75%, 13.30%, 10.86%, 14.04%). The paper should clarify whether the 'information gain' refers only to the conditional distribution of ΔR given twin-producing EOSs, rather than to the posterior probability of the twin hypothesis itself; otherwise the reader may expect this fraction to respond to measurement precision.

Circularity Check

1 steps flagged · score 4.0 of 10

Partial circularity: the 0.2-km saturation benchmark leans on same-group citations; the paper's own entropy and P_dis curves also tension the claim.

self citation load bearing [Section 2.2, final paragraph]
"once σ_R ≲0.2 km, not only the posterior distribution of ΔR but also those of all nine EOS parameters exhibit only weak additional evolution as the observational precision is further improved Li et al. (2024a); Grundler and Li (2025); Li et al. (2026)."

The central saturation threshold is argued to follow from the saturation of all nine EOS parameters, and that saturation is justified by citing three prior papers from the same group (Li, Grundler, Zhang, Xie). Those papers use the same meta-model EOS and the same Bayesian priors, so they are not independent, machine-checked, or externally falsified support for the present conclusion. The current paper's own Fig. 3 provides some independent evidence, so the circularity is partial rather than total; still, a load-bearing part of the derivation chain reduces to a self-citation chain.

full rationale

The analytical distinguishability model (Eqs. 1-5) is textbook Gaussian classification and is externally grounded; no circularity is found there. The Bayesian posterior analysis is an explicitly stated mock-data exercise with the prior and meta-model disclosed, and the 0.2-km conclusion is repeatedly qualified as 'within the present Bayesian framework'. The main circularity concern is the self-citation chain supporting the saturation: the paper invokes Li et al. (2024a), Grundler & Li (2025), and Li et al. (2026) to assert that the EOS-parameter posteriors have saturated below 0.2 km, and uses that to explain the central benchmark. These are same-author, same-framework papers, not independent external evidence. Because Fig. 3 of the present paper independently shows a flattening of the posterior mean, the self-citation is load-bearing but not the sole basis, so a score of 4 is appropriate rather than 6 or higher. Separately, and as a correctness concern rather than a circularity, Section 2.3 states that the Shannon entropy 'increases most rapidly between σ_R=0.1 and 0.2 km'. Since entropy quantifies residual posterior uncertainty, this means the largest information gain occurs when precision improves from 0.2 to 0.1 km, which contradicts the Summary's 'information-saturation regime for σ_R≲0.2 km'. Similarly, Eq. (1) with ΔR≈0.5 km gives P_dis≈0.89 at σ_R=0.2 km versus ≈0.994 at 0.1 km, so the misclassification rate drops by more than an order of magnitude. Those internal inconsistencies weaken the quantitative benchmark but are not a reduction of a prediction to its inputs, so they are flagged here rather than counted as an additional circular step.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The paper introduces no new physical entities. The free parameters are the mock data central value, the transition-density prior, and the inherited EOS meta-model parameters. The axioms are mostly domain assumptions about the EOS model family and Gaussian measurement errors, plus a hand-set prior. The most fragile premise is the inherited self-cited saturation claim that anchors the 0.2-km threshold.

free parameters (3)

Mock NS radius central value R_1.4 = 11.9 km
Chosen by hand as the true radius for the mock 1.4 M_sun measurement; the inferred ΔR posteriors and saturation threshold depend on this central value.
Transition density prior range = 1-6 rho_0 (uniform)
Prior on the hadron-quark transition density adopted from earlier work; the information content is measured relative to this prior.
Nine EOS meta-model parameters = Not specified in this paper
The nine parameters of the meta-model EOS are constrained by earlier Bayesian analyses; this paper reuses those posteriors as its effective prior, so the central result inherits their fitted values.

assumptions (5)

standard math Standard probability and information-theory identities (Bayes' rule, Gaussian CDF, Shannon entropy) are valid.
Used throughout Sections 2.2-2.4 without proof; standard textbook results.
domain assumption Radius measurement errors are Gaussian with known variance σ_R.
Section 2.2 assumes R_1 ~ N(μ_1, σ_R), R_2 ~ N(μ_2, σ_R) and the Bayesian likelihood uses Gaussian mock data.
domain assumption The EOS meta-model of Zhang et al. (2018) and Xie & Li (2020, 2021) spans the relevant space of supradense matter.
The paper's central conclusion is conditional on this model family; no comparison to alternative parameterizations is made.
ad hoc to paper Uniform prior on the hadron-quark transition density over (1-6)ρ_0.
A hand-set prior inherited from earlier same-group work; it is not derived from independent data and directly shapes the posterior saturation behavior.
ad hoc to paper Posterior saturation below σ_R ≈ 0.2 km, cited from Li et al. (2024a), Grundler & Li (2025), and Li et al. (2026), is correct and transferable to the ΔR posterior.
The paper relies on earlier self-cited results to explain why the ΔR posterior saturates; this is a load-bearing premise rather than a result demonstrated with new external data.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Quantifying the Information Gain from Future High-Precision Radius Measurements for Identifying Twin Neutron Stars." pith.science (2026). https://pith.science/paper/EKHESDH7

@misc{pith2026260718124,
  author       = {Pith},
  title        = {Pith review of: Quantifying the Information Gain from Future High-Precision Radius Measurements for Identifying Twin Neutron Stars},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/EKHESDH7}},
  note         = {Machine review of arXiv:2607.18124}
}

abstract

Twin neutron stars (NSs), characterized by identical gravitational masses but different radii, are among the most promising astrophysical signatures of a strong first-order hadron--quark phase transition in supradense matter. We investigate how increasingly precise NS radius measurements improve the Bayesian inference of twin-star observability using mock radius data for a canonical $1.4\,M_\odot$ NS. Radius uncertainties are varied from the current level of about $0.9$ km to the $\approx 0.1$ km precision anticipated from future X-ray and gravitational-wave observations. We quantify the information gained using the posterior distribution of the maximum twin-star radius separation $\Delta R$ together with an analytical model of branch distinguishability and complementary information-theoretic measures based on the branch observational efficiency and the Shannon entropy. The combined analyses reveal three inference regimes: a prior-dominated regime for $\sigma_R \gtrsim 0.6$ km, a rapid information-gain regime for $0.2 \lesssim \sigma_R \lesssim 0.6$ km, and an information-saturation regime for $\sigma_R \lesssim 0.2$ km. These complementary analyses consistently indicate that radius measurements with a precision of about $0.2$ km already extract most of the information available for identifying twin NSs within the present Bayesian framework. Beyond establishing a quantitative observational benchmark for future high-precision radius measurements, this work provides a general Bayesian framework for quantifying the information gain from progressively more precise observations and identifying the point of diminishing scientific returns.

Figures

Figures reproduced from arXiv: 2607.18124 by the authors.

**Figure 2.** Posterior probability distributions of ∆R and ∆M for Both (left) and Disconnected (right) categories from Bayesian inference using the mock radius data of R1.4 = 11.9 ± σR km data with σR = 0.1, 0.2, 0.5, 0.9 km, respectively. 0. 0 0. 1 0. 2 0. 3 0. 4 0. 5 0. 6 0. 7 0. 8 0. 9 1 . 0 1 . 1 2 1 . 1 4 1 . 1 6 1 . 1 8 1 . 20 1 . 22 1 . 24 Both category M e a n o f ∆ R ( k m ) σR (km) Maximum twi n star separati on [PITH… view at source ↗

**Figure 3.** Posterior mean of the maximum radius separation [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗

Figures from the paper (2 more)

**Figure 4.** Figure 4: The branch distinguishability probability [PITH_FULL_IMAGE:figures/full_fig_p004_4.png]

**Figure 5.** Figure 5: Normalized Shannon entropy, S/S (σR = 0.1) (black, left axis), and normalized branch observational efficiency, η/η(σR = 0.1) (blue, right axis), as functions of the assumed radius uncertainty σR for the Both category [PITH_FULL_IMAGE:figures/full_fig_p005_5.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

How Neutron Star Radii Encode the Dense-Matter Equation of State and Hadron-Quark Transition
astro-ph.HE 2026-08 conditional novelty 6.0 of 10

A Bayesian analysis of mock neutron star radius measurements shows that precision radii mainly constrain the symmetry-energy parameters L and Ksym and the quark transition density, while transition strength and quark ...

Pith tools

Reviewed August 1, 2026 · model on record in the stance chip above.