{"id":"6c1fb33c-201a-45f1-8283-1f563fb90a88","arxiv_id":"2507.11379","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"CDM3Y equations of state with incompressibility 230 to 330 MeV reproduce observed neutron-star tidal deformability constraints and yield a nearly universal linear relation between log tidal deformability and compactness.","lead":"This paper computes how neutron stars deform under tidal forces using a family of nuclear equations of state based on the M3Y nucleon-nucleon interaction. It finds that the stiffer members of this family, with incompressibility 230 to 330 MeV, match current gravitational-wave and X-ray constraints and give a simple linear relation between tidal deformability and stellar compactness.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Headline 'narrower ranges' lacks a defined statistical comparison; Table I mixes mass-posterior and K0 spreads, so the central claim is ambiguous and contradicted for GW190425 m2 on some readings.","rationale":"The most load-bearing element of the paper is not the microphysical composition assumption (which the authors explicitly qualify) but the advertised quantitative result: the CDM3Y-230..330 EOS band is narrower than the experimental ranges. If this claim is not well-defined, the abstract and conclusions overstate the paper's contribution. The ambiguity is visible in Table I: the '1.4 Msun' row and the 'GW190425 m2' row have essentially the same central mass (1.4 vs 1.405) but report Lambda = 311 ± 98 vs 978 ± 649, indicating that the latter includes propagation of the mass posterior while the former does not. The observed intervals are 90% credible intervals from nonparametric EOS inference [82]; comparing them with an ad hoc K0-envelope is apples-to-oranges. Under a fixed-mass reading, the GW190425 m2 Reid interval is as wide as the observed one; under an envelope reading, the CDM3Y union extends to higher Lambda than observed. No unique conclusion follows. This concern is testable and internal; it does not depend on whether phase transitions occur. The reader flagged the K0=230-330 selection and Eq. (22) overclaim as red flags; I agree with those but find the undefined 'range' statistic more central. The Eq. (22) universal relation is also under-supported (only 10 EOSs from one interaction family, fitted rather than predicted, and compared with other fits that differ by up to 40% in Fig. 3(d)), but it is secondary to the headline. Because the paper can be fixed by adding a clear statistical analysis, the verdict stays conditional acceptance, i.e., unchanged.","tokens_in":21934,"tokens_out":17045,"duration_ms":214283,"concrete_test":"Recompute Table I under two explicit definitions. (1) Fixed mass: for each event, set M to the posterior median and sample K0 uniformly from the CDM3Y family (first 230-330, then the full 200-330), computing the central 90% interval of Lambda and R. (2) Full posterior: draw M from the observed posterior and K0 from the same prior, and report the posterior predictive 90% interval. For each definition, compute the width ratio and overlap fraction with the observed 90% interval for GW170817 m1/m2 and GW190425 m1/m2. If the CDM3Y interval is not narrower under a clear definition (e.g., width ratio < 0.8) for at least one component, the headline claim should be revised to a per-object statement or removed. Also repeat with K0=200 included to test selection sensitivity.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline claim that the CDM3Y-230..330 EOSs yield 'more limited ranges of tidal deformability and radius ... than their experimentally inferred ranges' is not backed by a well-defined statistical comparison. Table I reports CDM3Y R and Lambda with uncertainties that appear to mix two very different sources: the spread over K0 at a fixed mass (e.g., the 1.4 Msun row: Lambda = 311 ± 98) and the spread induced by the component-mass posterior (e.g., GW190425 m2, nominally 1.405 Msun but Lambda = 978 ± 649). The experimental intervals, by contrast, are 90% credible intervals from EOS-agnostic inference that already marginalize over mass and waveform systematics. Because 'range' is never defined (envelope over K0, central 90% interval, or union over all objects), different natural definitions yield opposite answers: for GW190425 m2, the Reid CDM3Y interval (47-1453) is nearly identical in width to the observed interval (60-1433), and the Paris interval (329-1627) extends beyond it. Thus the central quantitative claim is ambiguous and, on at least one reading, contradicted by the paper's own table. The claim's force also depends on having preselected K0=230-330; including the softer K0=200 EOS would widen the band for low masses where it supports a star. This is not a mere style issue: it is the main advertised result.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the tidal deformability and compactness of neutron stars using the CDM3Y-Paris and CDM3Y-Reid equations of state with saturation incompressibility K0 between 200 and 330 MeV. The authors solve the TOV equations together with the tidal Love-number differential equation for beta-equilibrated npeμ matter with a unified crust, compare their predictions for radius and dimensionless tidal deformability against GW170817, GW190425, and NICER pulsar constraints, and conclude that the K0 = 230–330 MeV subset yields narrower ranges of radius and tidal deformability than the experimentally inferred ranges. They also examine the mass dependence of the compactness and the tidal Love number, and propose a linear relation, ln Lambda = 11.327 - 32.8 C (Eq. 22), for neutron stars with M >= M_sun, fitted to ten CDM3Y models. The paper is a standard application of well-established machinery to a specific family of finite-range interactions.","tokens_in":22266,"tokens_out":7860,"duration_ms":94431,"significance":"If the narrowing claim were rigorously established, the paper would provide a useful constraint on the CDM3Y interaction family and on the nuclear incompressibility from astrophysical data. The tabulated predictions for the GW170817 and GW190425 components, the comparison with NICER constraints, and the comparison with recent Skyrme and RMF calculations are informative reference material. The numerical scheme is standard and the manuscript is honest about the model dependence of its EOS. However, the central narrowing claim is not yet well defined, the proposed universal relation is tested only within one model family, and there is an internal inconsistency in the light-neutron-star behavior. These issues are fixable and the underlying calculations appear sound, so the paper is a reasonable candidate for revision rather than rejection.","major_comments":[{"comment":"The headline claim that the CDM3Y-230 to CDM3Y-330 EOSs give 'more limited ranges of tidal deformability and radius ... than their experimentally inferred ranges' is not backed by a defined statistical comparison. In Table I, the quoted uncertainties for the present calculations appear to mix two distinct sources: the spread over K0 at fixed mass and the spread induced by the component-mass posterior, while the experimental entries are 90% credible intervals that already marginalize over mass and waveform systematics. On the data as presented, the claim is false for at least one row: for GW190425 m2 the experimental Lambda interval is 60-1433 (521 +912/-461), whereas the Paris CDM3Y result is 978 ± 649, i.e. 329-1627, which extends beyond the observed interval; the Reid result, 47-1453, has essentially the same width as the observed interval. The authors should define the 'range' precisely (envelope over K0, central 90% interval, or union over objects), separate the K0 and mass-posterior contributions, and recompute the comparison for every object. The claim may survive for the 1.4 M_sun row, where the CDM3Y intervals are indeed narrower, but it cannot be stated as a general result based on the current table.","section":"Table I and Sec. III"},{"comment":"The narrowing claim depends on a post-hoc selection of the K0 = 230-330 MeV subset. The paper includes K0 = 200 MeV in the figures but states that 'only the soft EOS (K0 = 200 MeV) fails to reproduce the indicated Lambda for the NS objects of masses larger than 1.36 M_sun' and that its maximum mass does not exceed 1.4-1.5 M_sun. Excluding an EOS because it fails the very constraints that are later used to claim agreement makes the resulting 'narrower ranges' partly a consequence of the selection, not an independent prediction. The authors should report the full K0 = 200-330 MeV envelope and provide an explicit statement of whether the narrowing claim is meant as an in-sample description of the selected EOSs or as an out-of-sample prediction. In its present form, the central quantitative result is ambiguous and partially circular.","section":"Sec. III, K0-band selection"},{"comment":"The claim that Eq. (22), ln(Lambda) = 11.327 - 32.8 C with R^2 = 0.992, provides a 'unified description ... independent of the details of the employed EOS' is not supported by the evidence presented. The fit uses ten EOSs from a single family (CDM3Y-Paris and CDM3Y-Reid with K0 = 230-330 MeV), all based on the same functional form, the same npeμ composition, and the same assumption of no phase transitions. R^2 is an in-sample goodness-of-fit measure, not a test of universality across independent EOS families. Moreover, Fig. 3(d) itself shows that the various universal relations differ by up to 40% at C = 0.12 and about 24% at C = 0.27, with the present relation deviating from its own linear fit by 32% at C = 0.12. The authors should test Eq. (22) against at least a few independent EOS families and report the residuals and mass range of validity before calling it EOS-independent.","section":"Eq. (22) and Fig. 3(d)"},{"comment":"There is an internal inconsistency concerning the behavior of light neutron stars. The abstract states: 'For light NS (M < M_sun), both k2 and C decreases upon decreasing the NS mass, which enhances its tidal deformability.' In Sec. III, however, the authors report that upon decreasing the NS mass from 0.6 to 0.2 M_sun, the compactness increases from about 0.02 to 0.074 and the tidal Love number increases from about 0.02 to 0.082, and that these increasing rates strongly decrease the tidal deformability by about two orders of magnitude. These two statements cannot both be correct. Since Eq. (11) shows Lambda = (2/3) k2 / C^5, an increase in C with decreasing mass would dominate and reduce Lambda, not enhance it. The abstract and the discussion in Sec. III must be reconciled, and the direction of the effect should be stated correctly with the mass range clearly specified.","section":"Abstract and Sec. III"}],"minor_comments":[{"comment":"The caption contains the duplicated phrase 'MeV MeV'; this should be corrected to 'MeV'.","section":"Fig. 3 caption"},{"comment":"Several labels in Fig. 1 and Fig. 2 appear as corrupted text, e.g., 'npe/s109 NS matter'; these should be replaced with the intended notation, such as 'npeμ NS matter'.","section":"Fig. 1 labels"},{"comment":"Several references use 'el al.' instead of 'et al.' (e.g., Refs. [83], [85]-[90], [110]); these should be corrected for consistency.","section":"References"},{"comment":"The text attributes the inner-crust EOS to 'Douchen et al.' but the cited Ref. [79] is by Douchin and Haensel; the name should be spelled consistently with the reference list.","section":"Sec. III, crust EOS attribution"},{"comment":"The positivity condition for the core-crust transition is stated as K_mu_i > 0, but the text refers to 'the positivity condition of compressibility'; it would be clearer to state that the thermodynamic stability condition requires the appropriate incompressibility to be positive.","section":"Sec. II, Eq. (20)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a competent application of standard tools to a specific EOS family, but the advertised narrowing result needs a precise statistical formulation, and the universal-relation claim needs out-of-sample testing. The internal inconsistency between the abstract and Sec. III on light-neutron-star behavior should be corrected before the paper can be considered for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the two things to know. First, the paper is a careful, useful computation: it runs TOV plus the standard Love-number equation on the CDM3Y-Paris and CDM3Y-Reid EOS family across K0=200-330 MeV, and tabulates radii, tidal deformabilities, and Love numbers for canonical NS and for the GW170817/GW190425 components. That is a clean, reproducible addition to the literature, especially the Reid results and the K0>270 range, which the earlier Tan et al. scan didn't cover. The comparison with NICER and GW constraints is honest, and the qualitative findings (Love number rises with stiffness for M>Msun, falls for light stars) come out clearly.\n\nThe soft spot is the advertised headline. The abstract and Sec. IV claim the K0=230-330 EOSs give 'more limited ranges' of R and Lambda than experiment. The stress-test note is right: 'range' is never defined. Table I mixes two different spreads - the K0 spread at fixed mass and the component-mass posterior spread. For GW190425 m2 the CDM3Y-Paris band (329-1627) is actually wider than the observed 90% interval (60-1433), so on a natural reading the claim is false. The narrower result only holds cleanly for the 1.4 Msun canonical star, where CDM3Y gives Lambda=311±98 vs observed ~451(+241,-279). Also, the K0=230-330 window is selected after seeing what matches; including K0=200 widens the low-mass band. So the main conclusion as written is not robust. The authors need to either define the range comparison carefully (e.g., central 90% intervals at fixed mass) or soften the claim.\n\nThe Lambda-C relation in Eq. (22) is fine as a family-specific fit, but calling it 'EOS-independent' is too strong. It is a linear fit to ten EOSs from one interaction family, with the usual caveat that such relations are approximate. The comparison with Yagi-Yunes and other fits is a nice sanity check, though it adds little beyond existing results.\n\nThe math is standard and there are no apparent slips in the TOV/Love-number integration. The no-phase-transition assumption (uniform npe-mu matter) is standard for this kind of hadronic EOS paper; the authors even mention hybrid stars, so it's not hidden. The tidal constraints are external, so the main validation is not circular. There is some self-citation, but the cited prior work is part of the same program and the extension is real.\n\nVerdict: this deserves peer review, but with a major revision request on the statistical framing. I would not put much weight on the 'narrower ranges' claim until it is properly defined. The tables and qualitative trends are worth publishing. I'd bring it to reading group only if we want to talk about how easily 'range' can be ambiguous; otherwise it's a standard EOS-applications paper.","headline":"Solid CDM3Y tidal-deformability scan, but the headline 'narrower ranges' claim needs a defined statistical basis before it can be believed.","tokens_in":22846,"tokens_out":4546,"would_cite":false,"duration_ms":48985,"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 compact family of neutron-star equations of state passes tidal and radius tests.","keywords":["tidal deformability","neutron stars","equation of state","compactness","CDM3Y interaction","M3Y interaction","tidal Love number","gravitational waves"],"falsifier":"A secure measurement of the 1.4-solar-mass radius outside roughly 11.4–12.8 km, or of its tidal deformability outside roughly 200–430, would contradict the paper's central consistency claim; a confirmed neutron star above about 2.4 solar masses would exceed the maximum masses that the $K_0 = 230\\text{–}330$ MeV family allows.","tokens_in":21715,"feed_emoji":"🌌","tokens_out":6784,"duration_ms":76704,"temperature":0.7,"pith_summary":"This paper sets out to show that a single family of semi-microscopic nuclear equations of state, built from the M3Y finite-range nucleon-nucleon interaction in its density-dependent CDM3Y form, can simultaneously account for the tidal deformability inferred from the binary neutron-star mergers GW170817 and GW190425 and the radius constraints from X-ray pulse-profile observations. If true, the nuclear saturation incompressibility is pinned to 230–330 MeV, and the allowed radius and tidal deformability of neutron stars are narrower than current empirical bounds. The paper also claims an equation-of-state-independent linear relation between the logarithm of the dimensionless tidal deformability and compactness for stars heavier than one solar mass. A sympathetic reader would care because this links laboratory nuclear-matter parameters directly to what gravitational-wave detectors and X-ray telescopes measure.","feed_headline":"Neutron-star tidal tests favor one family of equations of state","feed_subtitle":"A single family of nuclear models matches merger and X-ray data, narrowing the allowed radius band.","key_machinery":"The engine of the calculation is the CDM3Y density-dependent form of the finite-range M3Y nucleon-nucleon interaction, parameterized to give equations of state with saturation incompressibility 200–330 MeV. These equations of state feed the Tolman-Oppenheimer-Volkoff structure equations together with the differential equation for the tidal Love number $y(r)$, whose surface value $y_R$ enters the compactness-dependent formula for $k_2$ and hence for the dimensionless tidal deformability $\\Lambda = 2k_2/(3C^5)$. The equation-of-state-independent claim is carried by a linear fit to $\\ln\\Lambda$ versus compactness for $M \\ge M_\\odot$, Eq. (22), which collapses the ten CDM3Y-Paris and CDM3Y-Reid curves onto one line.","core_discovery":"The central claim is that the CDM3Y-Paris and CDM3Y-Reid equations of state with saturation incompressibility $230 \\le K_0 \\le 330$ MeV reproduce the tidal deformability estimates from GW170817 and GW190425 and the X-ray pulse-profile radius constraints for massive pulsars, while producing a narrower band of radius and deformability than observation alone allows. The paper derives the relation $\\ln(\\Lambda) = 11.327 - 32.8\\, C$, with coefficient of determination $R^2 = 0.992$, for neutron stars with $M \\ge M_\\odot$, and interprets it as independent of the details of the equation of state. It further shows that for stars below one solar mass the tidal deformability is largely insensitive to the stiffness of the core, whereas above one solar mass stiffer matter enhances deformation because it lowers compactness while raising the tidal Love number.","pith_inferences":["Editorial extension: if the $\\ln\\Lambda$–$C$ relation holds beyond the CDM3Y family, a single precise radius or compactness measurement could be used to predict the deformability seen in the next binary neutron-star inspiral, sharpening the interpretation of gravitational-wave data.","Editorial extension: the paper's restriction to hadronic npe$\\mu$ cores sets up a testable extension—repeat the calculation with a quark or hyperon phase added; a hybrid-star branch that still fits GW170817 would show whether the 230–330 MeV band is genuinely preferred or merely one hadronic realization.","Editorial extension: the same machinery could be applied to the moment of inertia, which the paper mentions but does not correlate in the same equation-of-state-independent way, potentially yielding another compactness-based relation testable by pulsar timing."],"forward_implications":["The 230–330 MeV incompressibility window becomes a sharp prediction: radii of canonical neutron stars lie near 11.4–12.8 km, tighter than current empirical ranges.","The linear $\\ln\\Lambda$–$C$ relation gives a simple way to convert a compactness measurement into a tidal deformability prediction for any star heavier than one solar mass.","Stiffer equations of state produce larger radii and larger tidal deformabilities for massive neutron stars, so a future precise measurement of either quantity would discriminate among the family members.","The paper's maximum masses (about 2.1 solar masses for these EOSs) imply that any confirmed neutron star above roughly 2.4 solar masses would fall outside this hadronic family.","The results support the lower end of the inferred tidal-deformability constraints for the heavier merger components and the upper end for lighter components."],"supporting_citations":[{"why":"Supplies the CDM3Y-Paris and CDM3Y-Reid EOS family and its mass-radius profiles that the paper extends to tidal deformability.","marker":"[23]"},{"why":"Provides the GW170817 and GW190425 tidal-deformability constraints with which the computed Lambda values are compared.","marker":"[39]"},{"why":"Provides combined pulsar, gravitational-wave, and X-ray constraints on radii and tidal deformability used as empirical benchmarks.","marker":"[82]"},{"why":"Provides an X-ray pulse-profile constraint on radius and tidal deformability for a massive pulsar used in the comparisons.","marker":"[83]"},{"why":"Introduces the CDM3Y density-dependent form of the M3Y interaction on which all the EOSs in this paper are built.","marker":"[68]"},{"why":"Fixes the isovector density-dependence parameters and symmetry-energy properties of the CDM3Y EOSs.","marker":"[53]"},{"why":"Fixes the isoscalar density-dependence parameterization that determines saturation properties and the incompressibility range.","marker":"[54]"},{"why":"Reports the GW170817 detection whose observed signal motivates the tidal-deformability constraints used here.","marker":"[17]"},{"why":"Reports the GW190425 event whose component-mass tidal limits are compared with the EOS predictions.","marker":"[81]"}],"fun_headline_variants":["M3Y equations of state tighten neutron star radius and tidal limits","Stiffness matters only above one solar mass for neutron star tides","Semi-microscopic EOS narrow neutron star tidal and radius band","Universal compactness–deformability relation fits neutron star tide data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that neutron star cores contain only ordinary hadronic matter (protons, neutrons, electrons, and muons) in beta equilibrium, with no phase transition to quark or hyperon matter; if such a transition happens inside a real star, agreement with observed tidal constraints would not by itself validate the equation of state.","fun_headline_variants_meta":{"raw":{"variants":["M3Y equations of state tighten neutron star radius and tidal limits","Stiffness matters only above one solar mass for neutron star tides","Semi-microscopic EOS narrow neutron star tidal and radius band","Universal compactness–deformability relation fits neutron star tide data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000742,"raw_usage":{"total_tokens":3366,"prompt_tokens":1054,"completion_tokens":2312,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":670,"completion_tokens_details":{"reasoning_tokens":2236}},"tokens_in":670,"tokens_out":2312,"duration_ms":21956,"temperature":1.0,"reasoning_tokens":2236,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:09:26.375493+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A secure measurement of the 1.4-solar-mass radius outside roughly 11.4–12.8 km, or of its tidal deformability outside roughly 200–430, would contradict the paper's central consistency claim; a confirmed neutron star above about 2.4 solar masses would exceed the maximum masses that the $K_0 = 230\\text{–}330$ MeV family allows.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides combined pulsar, gravitational-wave, and X-ray constraints on radii and tidal deformability used as empirical benchmarks."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the CDM3Y density-dependent form of the M3Y interaction on which all the EOSs in this paper are built."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Fixes the isovector density-dependence parameters and symmetry-energy properties of the CDM3Y EOSs."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Fixes the isoscalar density-dependence parameterization that determines saturation properties and the incompressibility range."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the GW190425 event whose component-mass tidal limits are compared with the EOS predictions."}],"review_version":1}