{"id":"bda4009c-ffb4-4da6-a680-7907b146dc25","arxiv_id":"2608.02824","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":12,"one_line_summary":"For differentially rotating, strongly magnetized neutron stars on constant angular-momentum sequences, the normalized moment of inertia remains quasi-universal across equations of state when the Breu-Rezzolla fit is extended with magnetic and differential-rotation coefficients.","lead":"This paper tests whether neutron stars still obey a universal moment-of-inertia relation when they rotate differentially and carry very strong magnetic fields. It fits 600 equilibrium models and offers a corrected formula to improve magnetar spin-down and neutron-star merger estimates.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (12) is dimensionally inconsistent: the magnetic factor should be (Bmax R^2/M)^i, not (Bmax R^2/sqrt(GM))^i; as written it does not follow from Eq. (9) and overstates the magnetic correction by M^{i/2}.","rationale":"The central claim is that Eq. (12) is a quasi-universal relation. Before considering the deeper astrophysical questions, I checked the algebra connecting Eq. (12) to the defining fit Eq. (9), and the mismatch is immediate and quantitative. The reader's weakest assumption concerns the additive ansatz, toroidal-only fields, and generalizability to mixed configurations; those are legitimate physical concerns, but they are not the most load-bearing issue. The printed headline formula is internally inconsistent in a way that affects every use of the result, including the values advertised in the abstract and Section 5. The numerical fits behind Fig. 4 may still be correct if the authors actually used Eq. (8) with Eq. (9), but then Eq. (12) must be corrected to read (Bmax R^2/M)^i rather than (Bmax R^2/sqrt(GM))^i. Because this is a fixable algebraic error rather than a collapse of the underlying fitting idea, I keep the reader's CONDITIONAL verdict: the paper should not be accepted until Eq. (12) is corrected and the relation is re-validated, ideally with an out-of-sample check such as leaving out one EoS or one Bmax/A value.","tokens_in":14425,"tokens_out":14482,"duration_ms":127078,"concrete_test":"Re-derive Eq. (12) from Eq. (9): substitute epsilon = Bmax^2 R^4/M^2 and C = M/R into Eq. (9), rewrite the result as a coefficient times C^{-i}, and compare with Eq. (12). Then evaluate both expressions for one configuration from Fig. 4 (e.g., SLy4, A = 8, Bmax = 10^18 G, M = 1.4 M_sun, R = 12 km). If the magnetic contributions differ by the factor M^{i/2}, Eq. (12) is not the relation used to generate Fig. 4; the authors must state which equation was used and correct the other.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.3 defines epsilon = Bmax^2 R^4/M^2 and fits the magnetic part as Delta I~ = sum_i b_i (C/sqrt(epsilon))^{-i} (Eq. 9). Substituting C = M/R gives Delta I~ = sum_i b_i (Bmax R^3/M^2)^i. If this is rewritten as a coefficient multiplying C^{-i}, the magnetic coefficient must be b_i (Bmax R^2/M)^i, not b_i (Bmax R^2/sqrt(GM))^i. With G = c = 1, sqrt(GM) = sqrt(M) has dimension mass^{1/2} in cgs units (or length^{1/2} in geometric units), so the bracket in Eq. (12) is not dimensionless and cannot be added to the dimensionless coefficients a_i and b_i. Numerically, Eq. (12) as printed overestimates the magnetic term by a factor M^{i/2} relative to Eq. (9): for M = 1.4 M_sun, the factors are about 1.18, 1.40, 1.65, and 1.96 for i = 1..4. Therefore Eq. (12) cannot be the relation validated in Fig. 4, which was constructed from Eq. (8) together with Eq. (9). This is an internal algebraic inconsistency, independent of the separate question of whether purely toroidal fields are physically stable.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses the XNS code to construct stationary, axisymmetric equilibrium configurations of neutron stars with 10 hadronic equations of state, differential rotation described by the J-constant law, and purely toroidal magnetic fields, all on fixed-angular-momentum sequences. It proposes a unified phenomenological relation, Eq. (12), between the normalized moment of inertia I~=I/M^3 and compactness C=M/R, with coefficients depending on the differential-rotation parameter A and the maximum magnetic field Bmax. The authors report an EoS-insensitive scatter with average relative error ~8% and maximum relative error ~50% (Fig. 4), and apply the relation to estimate magnetar spin-down luminosities and the rotational kinetic energy of the GW170817 post-merger remnant. The central claim is that quasi-universality of the I~-C relation extends to highly magnetized, differentially rotating neutron stars.","tokens_in":14851,"tokens_out":8154,"duration_ms":68005,"significance":"If the proposed relation were robust, this would be a useful extension of the Breu-Rezzolla quasi-universal relations into a regime relevant to magnetars and binary-neutron-star merger remnants, and the use of 10 EoSs and fixed-J sequences is a reasonable basis for such a study. The paper also demonstrates a concrete phenomenological path for including differential rotation and magnetic fields in observable estimates. However, the central claim is currently undermined by an internal algebraic inconsistency in Eq. (12), by an A-dependence that is fitted through interpolation across only four A values, and by a validation procedure that re-fits the residuals it later reports as agreement. The maximum relative error of ~50% also tempers the term 'quasi-universal.' With corrections and additional validation, the framework could be valuable, but in its present form the main quantitative claims are not established.","major_comments":[{"comment":"Equation (12) as printed is dimensionally inconsistent and does not follow from Eq. (9). With G=c=1, epsilon=Bmax^2 R^4/M^2, so C/sqrt(epsilon) = M^2/(Bmax R^3) and the magnetic contribution in Eq. (9) is Delta I~ = sum_i b_i (C/sqrt(epsilon))^{-i} = sum_i b_i (Bmax R^3/M^2)^i. Rewriting this as coefficients multiplying C^{-i} gives b_i (Bmax R^2/M)^i, not b_i (Bmax R^2/sqrt(GM))^i. Since sqrt(GM)=sqrt(M) has dimensions of length^{1/2} in the G=c=1 units used throughout the paper, the bracket in Eq. (12) is not dimensionless and cannot be added to the dimensionless coefficients a_i and b_i. The printed form overestimates the magnetic term by M^{i/2}; if M is expressed in solar masses, the factors are about 1.18, 1.40, 1.65, and 1.96 for i=1..4. Consequently Eq. (12) cannot be the relation validated in Fig. 4, which was constructed from Eqs. (8)-(9). This error must be corrected and propagated through the applications in Sec. 4.","section":"Sec. 3.2, Eqs. (6)-(7) and Table 2"},{"comment":"The functions a_i(A) and b_i(A) are obtained by fitting four-parameter cubic polynomials in 1/A to only four values of A (infinity, 10, 8, 6). With four parameters and four data points, the fit is an interpolation, so the reported reduced chi-squared values below 0.02-0.03 do not validate the adopted functional form; in fact the reduced chi-squared has zero degrees of freedom for this interpolation. Moreover, the asymptotic limits of Eq. (7) do not reproduce the uniform-rotation coefficients in Table 2 (for example a1(infinity)=2.0895 versus the tabulated 2.1371, and similar discrepancies for a2-a4 and for the b_i limits in Eq. (10) versus Table 3), contradicting the statement that the relation 'naturally recovers' the uniform rotation limit as A tends to infinity. The A-dependence of Eq. (12) is therefore not established beyond the discrete tabulated values; additional A values or an out-of-sample test are needed.","section":"Sec. 3.3, Fig. 4"},{"comment":"The validation in Fig. 4 compares the unified relation with the same data used to construct it. Because Eq. (8) defines the magnetic contribution as the residual after subtracting the rotational fit, and Eq. (9) then fits that residual, the agreement in Fig. 4 is partly guaranteed by construction and does not independently test quasi-universality. The bottom panel of Fig. 4 also reports a maximum relative error of about 50%, which is large for a relation described as quasi-universal; the paper should identify where the worst outliers occur and whether they concentrate at extreme Bmax or at small C. An independent validation, such as leaving out one EoS or one A value and refitting the remaining data, is required to support the claimed 8% average and 50% maximum error estimates.","section":"Footnote 3 and Sec. 5"},{"comment":"The paper restricts the numerical analysis to purely toroidal magnetic fields and acknowledges in footnote 3 that such configurations are unstable. Since the proposed unified relation is intended to cover 'highly magnetized' neutron stars, the absence of any test with mixed or poloidal-dominated field configurations is a significant limitation. The additive decomposition Eq. (8) and the polynomial ansatz Eq. (9) are ad hoc and are tested only for toroidal fields; the paper should either provide evidence that the relation holds for mixed-field equilibria or state explicitly that the claimed universality is limited to toroidal-dominated configurations. The discussion in Sec. 5 should be expanded to reflect this caveat in the central claim.","section":null}],"minor_comments":[{"comment":"The entry 'Relatvistic mean-field' contains a typo; it should read 'Relativistic mean-field.'","section":"Table 1"},{"comment":"The sentence 'The deviation from the uniform rotation relation increases with with increasing differential rotation' contains a duplicated 'with.'","section":"Sec. 5, Discussion"},{"comment":"The phrase 'may reach upto 10^18 G' should be 'may reach up to 10^18 G.'","section":"Sec. 1, Introduction"},{"comment":"The color bar in Fig. 2 shows Bmax in the range 0-6 x 10^17 G, while the text states that Bmax is varied over 10^17-10^18 G; please check the range displayed in the figure and state it consistently.","section":"Fig. 2"},{"comment":"The sentence 'the current catalog of magnetars comprises of a limited number of observations' is awkward; 'comprises a limited number' or 'consists of a limited number' would be clearer.","section":"Sec. 1, Introduction"}],"recommendation":"major_revision","confidential_remarks":"The dimensional inconsistency in Eq. (12) is likely a typographical error that can be fixed, but the interpolation-based A-dependence and the in-sample validation of Fig. 4 require substantive additional work before the central quasi-universality claim is supported. I recommend major revision rather than rejection, as the manuscript's scope and data set are suitable for a revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper extends the Breu-Rezzolla I-C universal relation to toroidally magnetized, differentially rotating neutron stars on constant-J sequences. That is a new and plausible extension. But the headline Eq. (12) has a dimensional inconsistency that, as printed, does not follow from the paper's own Eq. (9). The stress-test note is correct. With that fixed, the paper is a decent numerical fitting study worth a serious referee.\n\nWhat is good: it uses 10 EoSs and the XNS code to map a parameter space (A = 6, 8, 10, infinity; Bmax 10^17 to 10^18 G) that has not been systematically covered before. The writing is clear, and the limitations section is unusually honest about toroidal-only fields, zero temperature, and the phenomenological nature of the fits. The applications to magnetar spin-down and GW170817 are appropriately framed as proof-of-concept.\n\nWhere it gets soft. First, Eq. (12) as printed: substituting epsilon = Bmax^2 R^4/M^2 into Eq. (9) gives a magnetic coefficient (Bmax R^2/M)^i, not (Bmax R^2/sqrt(GM))^i. With G=c=1, sqrt(GM)=sqrt(M) in geometric units, so the bracket is dimensionful and the equation overstates the magnetic term by M^(i/2). For M=1.4 solar masses that is factors ~1.18, 1.40, 1.65, 1.96 for i=1..4. The validation in Fig. 4 presumably used the correct form, so the printed equation does not match the actual fitted relation. It has to be corrected.\n\nSecond, the A-dependence is built from four-parameter cubics through only four A values. That is interpolation, not validation, and the reported reduced chi-squared values near zero are meaningless for testing the functional form. A referee should ask for either more A values or a simpler form with fewer parameters.\n\nThird, the magnetic correction is a residual fit: subtract the rotational fit, then fit the leftover. The maximum relative error of about 50% is large, and the relation is tested only for purely toroidal fields, which the authors themselves note are unstable. That limits claims of universality.\n\nFourth, no error bars are propagated anywhere, and the magnetar application uses a tuning parameter eta to map B_s to B_max, which can swing the spin-down luminosity by orders of magnitude.\n\nBottom line: the core idea is sound and the numerical work appears careful, but the printed central equation is wrong and the A-dependence is overfit. This deserves peer review, but it needs fixing before publication.","headline":"A useful extension of Breu-Rezzolla to magnetized, differentially rotating stars, but the printed central equation is dimensionally inconsistent and the A-dependence is overfit.","tokens_in":15401,"tokens_out":4719,"would_cite":false,"duration_ms":40391,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A quasi-universal relation connects moment of inertia and compactness in neutron stars even when they rotate differentially and carry magnetic fields up to 10^18 G.","keywords":["quasi-universal relations","neutron stars","moment of inertia","compactness","differential rotation","toroidal magnetic fields","magnetars","equation of state"],"falsifier":"Build equilibrium models with poloidal or mixed magnetic field configurations at $B_{\\max}$ up to $10^{18}$ G using a solver that permits them, and compare their $\\tilde I$–$C$ curves to Eq. (12); deviations larger than the paper's roughly 50% maximum relative-error band, or a failure of the additive decomposition, would refute the claimed quasi-universality.","tokens_in":14154,"feed_emoji":"🧲","tokens_out":7486,"duration_ms":58181,"temperature":0.7,"pith_summary":"The paper argues that the normalized moment of inertia $\\tilde I = I/M^3$ of a neutron star is set almost entirely by its compactness $C = M/R$, even for configurations that rotate differentially and carry strong toroidal magnetic fields. This extends a previously known slow-rotation universal relation by allowing the fit coefficients to depend on the differential-rotation parameter $A$ and on magnetic field strength through the dimensionless ellipticity $\\epsilon = B_{\\max}^2 R^4 / M^2$. Using 600 equilibrium models built with ten hadronic equations of state, the paper fits a single phenomenological formula with average relative error of about 8% and maximum relative error of about 50%. If the relation holds, it yields equation-of-state-insensitive estimates of neutron-star structure and shows that ignoring differential rotation and magnetic fields biases magnetar spin-down luminosities and post-merger rotational energies.","feed_headline":"Neutron-star spin formula holds even in strong magnetic fields","feed_subtitle":"Normalized moment of inertia tracks compactness within ~8 percent across ten equations of state, even for magnetars.","key_machinery":"The load-bearing object is the additive decomposition $\\tilde I(C,A,B_{\\max}) = \\tilde I_0(C,A) + \\Delta \\tilde I(C,B_{\\max})$, where the rotational part keeps the inverse-power form in $C$ and the magnetic part is written as a polynomial in $C/\\sqrt{\\epsilon}$ with $\\epsilon = B_{\\max}^2 R^4/M^2$. This scaling variable collapses the magnetic contribution onto a monotonic curve and forces the correction to vanish as $B_{\\max} \\to 0$, so the non-magnetized universal relation is recovered exactly. The coefficients of both polynomials carry all the dependence on differential rotation, parametrized by the length scale $A$ of the $J$-constant rotation law.","core_discovery":"The central claim is that the relation $$\\tilde I(C,A,B_{\\max}) = \\sum_{i=1}^4 \\left[a_i(A) + b_i(A)\\left(\\frac{B_{\\max}$R^{2}$}{\\sqrt{G M}}\\right)^i\\right] $C^{{-i}}$$$ remains valid for uniformly rotating, differentially rotating, and highly magnetized neutron stars alike, provided the sequences have constant angular momentum and lie on the stable branch with $j = J/M^2 \\lesssim 0.35$. The coefficients $a_i(A)$ and $b_i(A)$ are fitted third-order polynomials in inverse powers of $A$, recovering the uniform-rotation limit as $A \\to \\infty$. The paper demonstrates the relation by reproducing its 600 numerical configurations to within an average relative error of about 8% and a maximum of about 50%, with no strong dependence on the underlying equation of state. It then applies the relation to show that magnetic-field corrections can change inferred magnetar spin-down luminosities by orders of magnitude and that the rotational kinetic energy of the GW170817 remnant is constrained near $10^{53}$ erg, with frequency alone unable to break the degeneracy between differential rotation and magnetic field strength.","pith_inferences":["If the same quasi-universality holds for poloidal or mixed magnetic geometries, the relation would let observers infer interior field strengths of magnetars from moment-of-inertia measurements without knowing the equation of state.","The success of the $C/\\sqrt{\\epsilon}$ scaling hints that a virial-type argument, not just a polynomial fit, may underlie the magnetic correction; deriving it would turn the phenomenological formula into a predictive relation.","A natural stress test is to evolve the same constant-$J$ sequences with dynamical magnetic braking: because the relation assumes conserved angular momentum, systems losing $J$ through winds or gravitational waves should deviate in a way the current formula does not capture.","Combining the relation with a second equation-of-state-insensitive observable, such as tidal deformability, may break the degeneracy between differential rotation and magnetic field strength that spin frequency alone cannot resolve."],"forward_implications":["The same fitting formula can be used to estimate $\\tilde I$ for a neutron star when only $M$, $R$, $B_{\\max}$, and the rotation law are known, regardless of which of the ten equations of state is correct.","Magnetar spin-down luminosities inferred from $P$ and $\\dot P$ shift by orders of magnitude once the magnetic-field-corrected moment of inertia replaces the canonical value, so catalog comparisons need a structural correction.","For a fixed $C$, stronger differential rotation lowers $\\tilde I$, meaning constant-angular-momentum sequences with the same compactness have systematically smaller moments of inertia than uniformly rotating stars.","Post-merger remnant analyses that assume uniform rotation and weak magnetization will misestimate rotational kinetic energy; constraining the spin frequency alone leaves the differential-rotation parameter and the magnetic field strength degenerate."],"supporting_citations":[{"why":"Supplies the original inverse-power slow-rotation relation between $\\tilde I$ and $C$ that this paper extends.","marker":"[29]"},{"why":"Provides the $J$-constant differential rotation law used for all equilibrium sequences.","marker":"[38]"},{"why":"Gives the constant-angular-momentum sequence construction followed in the analysis.","marker":"[46]"},{"why":"Underlies the equilibrium solver's extended conformal flatness formulation used to generate models.","marker":"[43]"},{"why":"Defines the magnetic ellipticity $\\epsilon$ used to scale the magnetic contribution.","marker":"[48]"},{"why":"Provides the public magnetar catalog supplying the spin-down data for the luminosity bias estimate.","marker":"[20]"},{"why":"Provides the GW170817 remnant energy estimate against which the paper compares its rotational kinetic energy.","marker":"[49]"}],"fun_headline_variants":["Universal spin-compactness relation holds for magnetars","Magnetar luminosity estimates refined by new universal law","Inertia-compactness universality survives differential rotation","GW170817 remnant energy pinned with quasi-universal rule","Neutron-star spin law extends to sheared, magnetized stars"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The magnetic correction is assumed, not derived, to separate additively from rotation and to follow the chosen polynomial in $C/\\sqrt{\\epsilon}$, and it is calibrated only for purely toroidal fields.","fun_headline_variants_meta":{"raw":{"variants":["Universal spin-compactness relation holds for magnetars","Magnetar luminosity estimates refined by new universal law","Inertia-compactness universality survives differential rotation","GW170817 remnant energy pinned with quasi-universal rule","Neutron-star spin law extends to sheared, magnetized stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000214,"raw_usage":{"total_tokens":1461,"prompt_tokens":1016,"completion_tokens":445,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":632,"completion_tokens_details":{"reasoning_tokens":362}},"tokens_in":632,"tokens_out":445,"duration_ms":4261,"temperature":1.0,"reasoning_tokens":362,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:59:10.185458+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Build equilibrium models with poloidal or mixed magnetic field configurations at $B_{\\max}$ up to $10^{18}$ G using a solver that permits them, and compare their $\\tilde I$–$C$ curves to Eq. (12); deviations larger than the paper's roughly 50% maximum relative-error band, or a failure of the additive decomposition, would refute the claimed quasi-universality.","supporting_citations":[{"cited_title":"Effects of Differential Rotation on the Maximum Mass of Neutron Stars","cited_arxiv_id":"gr-qc/0210012","evidence_quote":"Gives the constant-angular-momentum sequence construction followed in the analysis."},{"cited_title":"Komatsu, Y","cited_arxiv_id":null,"evidence_quote":"Underlies the equilibrium solver's extended conformal flatness formulation used to generate models."},{"cited_title":"Thompson and R","cited_arxiv_id":null,"evidence_quote":"Provides the public magnetar catalog supplying the spin-down data for the luminosity bias estimate."}],"review_version":2}