{"id":"3e46cf71-0272-4823-9ca0-942e67dc84a9","arxiv_id":"1908.02638","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"If neutron matter becomes spin-polarized inside neutron stars, the maximum possible mass would be capped at 2.6 to 2.9 solar masses.","lead":"This paper argues that a transition to spin-polarized neutron matter would cap the maximum neutron star mass at 2.6 to 2.9 solar masses. It combines microscopic nuclear force calculations with GW170817 observations to show that large spin-polarized cores in neutron stars are disfavored.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Upper mass bound rests on unvalidated extrapolation of the spin-polarized EOS beyond 2 nsat, where the phase transition occurs for stiff EOSs.","rationale":"The paper is a careful, honest theoretical study: it constructs a broad EOS family, includes three independent SPM calculations, tests a Gibbs-like smearing of the transition, and states its caveats. The reader's conditional verdict is appropriate. The single most load-bearing technical point is the reliability of the SPM EOS at the transition density for stiff unpolarized EOSs. Because the maximum-mass bound is obtained by replacing the unpolarized EOS with the SPM EOS once the latter becomes favorable, the crossing density controls Mmax; for EOSs that would otherwise reach 3–4 M⊙, the crossing is expected above the 2 nsat range of the microscopic SPM points. Eq. (1) is a plausible phenomenological form, and the AFDMC check over 1–2 nsat is a real supporting test, but it cannot certify the extrapolation to the densities that set the 2.9 M⊙ endpoint. A beta-equilibrated treatment of the SPM phase would be a further caveat, but the paper's Gibbs-smearing sensitivity analysis supports the direction that including it would lower Mmax, so the extrapolation uncertainty is the sharper concern. The proposed re-computation with alternative high-density continuations directly tests whether the central bound is an artifact of the fit.","tokens_in":11425,"tokens_out":7746,"duration_ms":91347,"concrete_test":"Compute Mmax for the stiffest unpolarized EOS in the speed-of-sound band using three different high-density continuations of the AFDMC VE1 SPM band beyond 2 nsat: (i) the paper's Eq. (1) fit; (ii) a one-parameter polytrope anchored to the last two Table 1 points; (iii) a linear extrapolation of E/N in density. Record the phase-transition density and Mmax for each case. If the transition density lies above 2 nsat and Mmax stays below 2.9 M⊙ in all three continuations, the extrapolation concern is resolved; if Mmax varies by more than about 0.1–0.2 M⊙ or exceeds 2.9 M⊙, the upper bound is an artifact of the chosen fit and the claim should be reported as conditional on the extrapolation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is the high-density SPM EOS used in the Maxwell construction. The SPM energy is a two-term power-law fit, Eq. (1), anchored to microscopic calculations that reach only about 2 nsat. For the stiff unpolarized EOSs that would otherwise give the largest Mmax, the unpolarized and SPM branches cross above 2 nsat, so the bound Mmax ≤ 2.6–2.9 M⊙ is set by the extrapolated SPM energy rather than by the microscopic input. The upper AFDMC bound, which produces the 2.9 M⊙ end of the range, comes from the upper edge of the chiral-EFT uncertainty band and is precisely the quantity least constrained by the calculations. The in-paper check described in §3 (fit to n ≤ nsat, compare with data up to 2 nsat) validates the extrapolation only within the fitted density range; it does not certify the form beyond 2 nsat. If the true SPM energy beyond 2 nsat is slightly stiffer than the extrapolated upper bound, the transition density moves upward and Mmax can exceed 2.9 M⊙. The authors remark that the stiffest AFDMC case is likely overestimated by regulator artifacts, which is plausible but not a quantitative bound on the extrapolation error. The paper is internally consistent and explicitly conditional, but the headline limit as stated is not robust to the functional form of Eq. (1) outside the microscopic range.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript argues that a first-order phase transition from unpolarized neutron-star matter to fully spin-polarized neutron matter (SPM) generically cuts off the high-pressure branch of the neutron-star equation of state, limiting the maximum mass Mmax to approximately 2.6-2.9 Msun. The authors combine a general ensemble of low-density EOSs from chiral EFT/Quantum Monte Carlo with three microscopic calculations of the SPM energy (AFDMC with local chiral N2LO interactions, MBPT with chiral N3LO interactions, and BHF with phenomenological potentials), fit the SPM energies with the two-power-law form of Eq. (1), and construct hybrid EOSs via a Maxwell construction. They then solve the Tolman-Oppenheimer-Volkoff equations and compare the resulting mass-radius relations and Mmax values with constraints inferred from GW170817. The central, observation-independent conclusion is that an EOS that would otherwise support Mmax above about 2.6-2.9 Msun is viable only if the SPM transition is absent; once the transition is included, the stable branch ends earlier and the maximum mass is lowered.","tokens_in":11741,"tokens_out":7869,"duration_ms":92293,"significance":"If the calculation is correct, this is an important observation-independent theoretical bound: it is more restrictive than generic speed-of-sound extensions of the EOS and complementary to the empirical GW170817 limits. The paper's strengths are that the SPM EOS is computed from several independent many-body methods with explicit uncertainty bands; the Maxwell construction is varied, and smearing the transition only lowers Mmax, so the bound is conservative with respect to the transition width; and the GW170817 comparison is used after the fact rather than as input. The data in Table 1 and the enumeration of EOS variants make the calculation transparent. The main weakness is the reliance on Eq. (1) outside the density range where the microscopic SPM calculations are actually constrained.","major_comments":[{"comment":"The headline upper bound (Mmax about 2.9 Msun) is set by the extrapolated SPM energy above 2 nsat. For the stiff unpolarized EOSs that would otherwise give the largest masses, the crossing between the unpolarized and SPM branches can occur above 2 nsat, i.e., beyond the density range of the AFDMC and MBPT calculations. The fit-quality check described in Section 3 (fit to n <= nsat and compare with data up to 2 nsat) validates the functional form only inside the fitted density range; it does not certify the two-power-law behavior beyond 2 nsat. If the true SPM energy rises slightly more steeply than the extrapolated upper band, the transition density moves upward and Mmax can exceed the quoted limit. Please either restrict the claim to a value robust against this extrapolation (about 2.6 Msun from the MBPT/BHF results) or provide a quantitative estimate of the extrapolation error above 2 nsat.","section":"Section 3, Eq. (1)"},{"comment":"The upper end of the reported range comes from the AFDMC upper uncertainty band, which the authors themselves describe as 'most likely overestimated' by local regulator artifacts. Quoting 2.6-2.9 Msun as the limit therefore mixes a robust result with a source explicitly identified as unreliable. The authors should either quantify the regulator contribution to that band or report the 2.9 Msun value only as an unvalidated extreme, with the firm theoretical bound given by the MBPT/BHF results (around 2.6 Msun).","section":"Section 3, AFDMC upper bound"}],"minor_comments":[{"comment":"The rendering '2.6$-$$2.9' should be corrected to a proper en-dash range.","section":"Abstract"},{"comment":"The text contains 'many many follow-up observations'; this should be edited to a single 'many'.","section":"Section 1"},{"comment":"The statement that the mass of the SPM domain is '<= 0.005 Msun, largely a result of numerical discretization artifacts' is confusing; please clarify whether these are stars with essentially no physical SPM core and explain how the discretization affects the quoted bound.","section":"Section 3"},{"comment":"For reproducibility, please provide the fitted parameter values (a, alpha, b, beta) for each SPM EOS variant, or at least specify the fit ranges and weighting used for the upper/lower bounds.","section":"Section 3, Eq. (1)"},{"comment":"The caption should state explicitly that the gray areas are the baseline EOS band without the SPM transition, while the hatched areas are the SPM-inclusive bands; the meaning of the solid red line in panel (a) (the centroid) is clear but should be stated for all panels.","section":"Figure 1 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is transparent about its caveats and the reference list is appropriate. The main issue for the editor is whether the upper end of the claimed bound (2.9 Msun) should appear in the abstract, given that it rests on an extrapolation beyond the microscopic range and on an uncertainty band the authors themselves regard as likely overestimated. I see no circularity or novelty concerns; the required changes are to harden or soften the headline claim, which is within the scope of a major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe short version: Tews and Schwenk argue that a phase transition to spin-polarized neutron matter (SPM) caps the maximum neutron-star mass at roughly 2.6–2.9 solar masses, and they back it with new AFDMC calculations of SPM using chiral interactions. The idea is genuinely new—previous SPM work didn't connect to the neutron-star mass-radius relation this way—and the paper is careful to separate the microphysical input from the observational comparison. GW170817 enters only after the bound is derived, so there's no circularity.\n\nWhat's good: The AFDMC calculation of SPM is a real addition, and they validate the density extrapolation of Eq. (1) against the calculated points between nsat and 2 nsat. The Maxwell construction and the ensemble of unpolarized EOSs are handled cleanly; the Gibbs-smearing check shows the Maxwell construction gives the most conservative (highest) Mmax, so the bound isn't an artifact of the transition construction. The paper is also honest about the proton fraction corrections (~10%) and the magnetic-field caveat.\n\nThe soft spot is the one the authors themselves flag: the SPM EOS is only computed to about 2 nsat, and for the stiffest unpolarized EOSs the phase crossing sits above that. The upper end of the claimed range, 2.9 M⊙, comes from the upper edge of the AFDMC uncertainty band—exactly the part least constrained by the underlying calculation. The in-paper extrapolation test only certifies the fit within the fitted density range, not beyond 2 nsat. So the headline 'limit' is a conditional statement, not a theorem. If the true SPM EOS turns out somewhat stiffer above 2 nsat, stars above 2.9 M⊙ are not excluded. That doesn't kill the paper—the authors say as much—but it means the bound should be read as a physics-motivated conjecture with quantified uncertainty, not a hard microphysical cap.\n\nBottom line: This is a solid, serious paper that deserves a real referee. It will be cited because it gives the dense-matter community a concrete mechanism and an explicit (if conditional) mass ceiling that agrees with astrophysical inferences. I'd want the extrapolation sensitivity discussed further, but I would not desk-reject it.","headline":"A physically motivated cap on the neutron-star maximum mass from spin-polarized matter, with new AFDMC results, though the bound is only as strong as the extrapolated polarized EOS.","tokens_in":12243,"tokens_out":2613,"would_cite":true,"duration_ms":26978,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["26.60.-c","97.60.Jd"],"model":"deepseek-v4-flash","headline":"A phase transition to spin-polarized neutron matter would cap neutron-star masses below 2.6-2.9 solar masses.","keywords":["neutron stars","equation of state","spin-polarized neutron matter","maximum mass","GW170817","chiral effective field theory","quantum Monte Carlo","phase transition"],"falsifier":"A clean falsification would be a neutron star with a precisely measured mass above about 2.9 solar masses, or a measurement showing that the fully polarized phase is stiffer at high density than the near-free-Fermi-gas behavior assumed here.","tokens_in":11248,"feed_emoji":"⭐","tokens_out":11264,"duration_ms":104977,"temperature":0.7,"pith_summary":"The paper sets out to show that dense neutron-star matter, if it switches into a spin-polarized phase at high density, cannot support neutron stars heavier than about 2.6-2.9 solar masses. This matters because most general extrapolations of the nuclear equation of state permit masses up to 3-4 solar masses, while the new transition removes those extreme stars and lands naturally in the mass range inferred from the gravitational-wave and electromagnetic observations of GW170817. The argument combines microscopic many-body calculations of spin-polarized neutron matter with a Maxwell construction that joins the ordinary and polarized phases at equal pressure and chemical potential. The authors find that the transition softens the equation of state so strongly that the maximum neutron-star mass is set by the onset of the polarized phase, not by other extension schemes.","feed_headline":"Spin-polarized cores cap neutron stars below 2.9 solar masses","feed_subtitle":"A phase transition to fully aligned neutron spins would cap the heaviest neutron stars, matching GW170817.","key_machinery":"The load-bearing object is the equation of state of spin-polarized neutron matter (SPM), neutron matter in which all spins are aligned so that a single spin state is occupied and Pauli blocking makes interactions weak. The paper computes SPM with three independent methods: auxiliary-field diffusion Monte Carlo with local chiral forces, many-body perturbation theory with chiral $N^3$LO interactions, and Brueckner-Hartree-Fock with phenomenological potentials. Each result is extended to higher density with the functional form $E_{\\rm pol}/N = a(n/n_{\\rm sat})^\\alpha + b(n/n_{\\rm sat})^\\beta$, and for every allowed unpolarized neutron-star equation of state from a causal speed-of-sound extension, a Maxwell construction transitions to SPM at the crossing pressure. The crossing of the two equations of state is what sets the end of the stable branch and therefore the maximum mass.","core_discovery":"The central claim is that a phase transition from unpolarized neutron-star matter to fully spin-polarized neutron matter in the core limits the maximum mass of neutron stars to $M_{\\rm max}\\lesssim 2.6$-$2.9\\,M_\\odot$, depending on which microscopic nuclear force is used in the spin-polarized calculation. Without the transition, the same general speed-of-sound extension of the equation of state allows $M_{\\rm max}$ up to about $3$-$4\\,M_\\odot$. The mechanism is that interactions in spin-polarized neutron matter are weak, close to a free Fermi gas, so the pressure rises only slowly once that phase is thermodynamically preferred, and the stable branch of the mass-radius relation ends near the crossing point. The paper further finds that the masses thus obtained agree with the upper limits inferred from the kilonova and gravitational-wave signal of GW170817, and that neutron stars containing a large spin-polarized core are ruled out by the radius constraint from the same event.","pith_inferences":["A testable extension: a precise measurement of a neutron star above roughly $2.9\\,M_\\odot$ would contradict the predicted ceiling, while a star near $2.5\\,M_\\odot$ with a small radius would sharpen the crossing density.","The same Maxwell-construction machinery could be applied to other candidate high-density phases; the paper notes quark matter cannot yet give a firm bound, but a future first-principles quark-matter equation of state would provide a natural point of comparison.","If the spin-polarized phase is realized, radius measurements near the maximum mass should show a kink or flattening as mass approaches the ceiling, distinguishing this transition from smoother equations of state."],"forward_implications":["If the transition exists, no isolated neutron star can exceed about $2.9\\,M_\\odot$; the heaviest measured neutron stars must sit near or below that ceiling.","General equation-of-state extension schemes that allow arbitrarily stiff pressure at high density overestimate the maximum mass unless the spin-polarized phase is absent or much stiffer than current calculations indicate.","Neutron stars with a substantial spin-polarized core are effectively excluded: the mass in the polarized phase is at most about $0.02\\,M_\\odot$ once GW170817 radius constraints are applied.","The theoretical upper bound and the independent upper limits from GW170817 reinforce each other, because the theoretical bound is derived without using any merger information."],"supporting_citations":[{"why":"Supplies the speed-of-sound extension and the family of unpolarized equations of state that define the allowed band, and provides the GW170817 radius limit used to rule out stars with spin-polarized cores.","marker":"Tews et al. 2018b"},{"why":"Provides the chiral many-body perturbation theory equation of state for spin-polarized neutron matter.","marker":"Krüger et al. 2015"},{"why":"Provides the Brueckner-Hartree-Fock spin-polarized equations of state and the earlier conclusion that polarized matter does not appear in neutron-star cores.","marker":"Vidana et al. 2002"},{"why":"Supplies the local chiral Hamiltonians and Quantum Monte Carlo framework on which the auxiliary-field diffusion Monte Carlo calculations are based.","marker":"Lynn et al. 2016"},{"why":"Introduces the functional form used to extrapolate spin-polarized energies per particle to higher density.","marker":"Gandolfi et al. 2009"},{"why":"Reports the GW170817 gravitational-wave signal that underlies the radius and remnant constraints used for comparison.","marker":"Abbott et al. 2017a"},{"why":"Gives the kilonova-based upper bound on maximum mass that the paper's results are compared against.","marker":"Margalit & Metzger 2017"},{"why":"Gives a merger-based maximum-mass estimate around 2.2 solar masses used as a comparison.","marker":"Shibata et al. 2017"},{"why":"Gives empirical relations between uniformly and differentially rotating maximum masses used to bound Mmax from GW170817.","marker":"Rezzolla et al. 2018"}],"fun_headline_variants":["Spin-polarized phase caps neutron star mass below 2.9 suns","Spin-aligned cores limit neutron star mass to 2.6-2.9 suns","Neutron star mass capped by spin-polarized phase, GW170817 agrees","Spin-polarized matter sets neutron star max mass near 2.9 suns"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The bound depends on the spin-polarized neutron-matter equation of state being reliable at the density where the phase transition happens, and for stiff unpolarized equations of state that density lies beyond the range of the microscopic calculations and is reached by a simple two-term power-law extrapolation.","fun_headline_variants_meta":{"raw":{"variants":["Spin-polarized phase caps neutron star mass below 2.9 suns","Spin-aligned cores limit neutron star mass to 2.6-2.9 suns","Neutron star mass capped by spin-polarized phase, GW170817 agrees","Spin-polarized matter sets neutron star max mass near 2.9 suns"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000638,"raw_usage":{"total_tokens":2919,"prompt_tokens":903,"completion_tokens":2016,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":519,"completion_tokens_details":{"reasoning_tokens":1927}},"tokens_in":519,"tokens_out":2016,"duration_ms":13655,"temperature":1.0,"reasoning_tokens":1927,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:38:52.789422+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A clean falsification would be a neutron star with a precisely measured mass above about 2.9 solar masses, or a measurement showing that the fully polarized phase is stiffer at high density than the near-free-Fermi-gas behavior assumed here.","supporting_citations":[{"cited_title":"2002, Phys","cited_arxiv_id":null,"evidence_quote":"Provides the Brueckner-Hartree-Fock spin-polarized equations of state and the earlier conclusion that polarized matter does not appear in neutron-star cores."}],"review_version":1}