{"id":"a43adbd4-75e6-4f71-a0fd-a19f64275548","arxiv_id":"2608.13236","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"White dwarfs in minimal dilatonic gravity are sub-Chandrasekhar for all scalar Compton lengths, with maximum mass dropping from 1.425 to 1.09 solar masses at λΦ = 500 km.","lead":"Researchers compute white dwarf models in a modified gravity theory with a massive scalar field, finding they are always lighter than the standard Chandrasekhar limit. The same coupling would create detectable altitude differences in Earth-orbiting satellites, so the theory can be tested with existing data.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the sub-Chandrasekhar claim is robust, though one supporting verification statement is inaccurate.","rationale":"The reader's weakest assumption was the linearized-exterior Robin condition, worrying that the mass eigenvalues, m_*/m_tot split, and deep-screening results depend on it. After reviewing the paper, this concern is real but not load-bearing for the central claim. The maximum mass as a function of λ_Φ is determined by configurations with R_*/λ_Φ ≲ 5, where the shooting method is reliable and independent of the Robin condition; the two methods agree to high precision there. The deep-screening regime is only needed for low-density stars, whose masses are far below the Chandrasekhar limit anyway, so an error there would not create a super-Chandrasekhar branch. The analytic λ_Φ→∞ limit provides a rigorous anchor for the strongest suppression. The paper's Sec. 5.1 statement that |Φ_c−1| ≤ 10⁻⁴ 'throughout' is demonstrably wrong for high-density models — the compactness alone gives |Φ_c−1| ~ 10⁻³ — but this is a verification error, not a soundness error, because the linearization error is quadratic in the small field and thus ~10⁻⁶ of the masses. The code is not yet public, which hampers reproducibility but is not a technical flaw. I therefore see no change to the ACCEPT verdict; a minor correction to the verification claim would be in order.","tokens_in":16266,"tokens_out":27330,"duration_ms":259493,"concrete_test":"Recompute the high-density maximum-mass sequence for λ_Φ = 1500 km using a full nonlinear exterior BVP instead of the linearized Robin condition, e.g., by extending the shooting method to the exterior with a matched asymptotic expansion and bisecting on the absence of the growing mode; if m* or m_tot shifts by more than 0.1%, the sub-Chandrasekhar claim and the disphere splits would need re-examination.","verdict_should_be":"UNCHANGED","load_bearing_attack":"No load-bearing objection found. The central claim that MDG white dwarfs are strictly sub-Chandrasekhar for every λ_Φ is supported by two independent solvers, by exact polytropic checks, and by the analytic λ_Φ→∞ limit (G_eff = 4/3 G_N, mass ratio 0.6495) which the numerics match to 0.06%. The collocation method's linearized-exterior Robin condition (Eq. 9) is not load-bearing for the maximum-mass sequence: those configurations have R_*/λ_Φ ≲ 5, where the shooting method (which does not use the Robin condition) is valid and agrees with collocation to 10⁻⁴–10⁻⁶. The deep-screening regime (R_*/λ_Φ > 20) where only collocation works concerns low-density stars, not the maximum mass. One internal inconsistency: Sec. 5.1 claims |Φ_c−1| ≲ 10⁻⁴ 'throughout', but high-density models (Table 4, λ_Φ = 1500 km, ρ_c = 3.2×10¹⁰ g cm⁻³) have compactness GM/(Rc²) ~ 10⁻³, implying |Φ_c−1| ~ 10⁻³. The linearization remains valid at the 10⁻³ level (nonlinear corrections ~10⁻⁶), so this does not threaten the central result, but the verification statement should be corrected.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies static, spherically symmetric white dwarfs in minimal dilatonic gravity (MDG), a Brans–Dicke theory with ω=0, fixed coupling α^2=1/3, and a free Compton length λΦ. It solves the relativistic structure equations with two independent numerical methods (a shooting method and a new free-boundary collocation method with a linearized-exterior Robin condition), using polytropic, Chandrasekhar, and Salpeter-corrected equations of state. The central finding is that MDG white dwarfs are strictly sub-Chandrasekhar for all λΦ, with the maximum mass decreasing from 1.425 M⊙ in GR to 1.27 M⊙ at λΦ=200 km and 1.09 M⊙ at 500 km, and with the suppression ratio nearly independent of the equation of state. Including rigid rotation near mass shedding, the most massive observed white dwarfs bound λΦ ≲300 km; a gravitational-redshift mass–radius test gives λΦ ≲720 km. The paper also predicts an altitude-dependent Kepler-inferred GM⊕ from the Earth-orbit dilaton field and introduces the distinction between the stellar mass m* and the total gravitating mass m_tot that includes the exterior dilaton 'disphere.'","tokens_in":16603,"tokens_out":5483,"duration_ms":53524,"significance":"If the central result holds, it is significant: it rules out super-Chandrasekhar white-dwarf branches in this one-parameter theory and turns the observed massive white-dwarf population into a direct constraint on λΦ. The paper's strengths are substantial: two independent solvers agree to 1e-6–1e-4 along mass–radius sequences; the GR limit recovers exact Lane–Emden polytropes; the analytic λΦ→∞ limit (G_eff=4G/3, mass ratio 0.6495) is matched to 0.06%; and the EOS dependence of the suppression ratio is small (at most 0.26 percentage points between Chandrasekhar and Salpeter). The free-boundary collocation method is a useful technical contribution that removes the exponential stiffness of shooting methods. The potential weakness concerning the linearized-exterior Robin condition in Eq. (9) does not, on close reading, endanger the maximum-mass claim, because the maximum-mass configurations lie in the regime where the shooting method is valid and the two methods agree; the deep-screening regime where only collocation operates concerns low-density stars, not the mass limit.","major_comments":[],"minor_comments":[{"comment":"The verification statement that \"the central field remains about 1 (|Φ_c−1| ≲ 10^-4) throughout\" is contradicted by the high-density rows of Table 4: for λΦ=1500 km and ρ_c=3.2×10^10 g cm^-3, the compactness GM/(Rc^2) is about 10^-3, implying |Φ_c−1| ∼ 10^-3. The linearization underlying the Robin condition remains valid at the 10^-3 level because nonlinear corrections are much smaller, but the sentence should be corrected or explicitly restricted to the low-density models.","section":"Sec. 5.1, paragraph (ii)"},{"comment":"The caption of the middle panel states that I(λΦ) is shown for ρ_c=10^7 g cm^-3, while the text reports the sequence I/I_GR = 0.959, 0.865, 0.666, 0.584 for \"the ρ_c=10^9 g cm^-3 star.\" One of these density values should be corrected so that the figure and text agree.","section":"Fig. 5 caption and Sec. 5.4"},{"comment":"The sentence \"the light-field limit is therefore the strong-coupling limit of the theory\" is confusing as written: λΦ→∞ makes the dilaton field unscreened and hence effectively strongly coupled to matter, but calling this the \"light-field limit\" while also calling it the \"strong-coupling limit\" in the same sentence could mislead readers. Please rephrase to distinguish the field mass (small, hence light) from the effective gravitational coupling (large, 4G/3).","section":"Sec. 5.2"},{"comment":"The text contains a garbled rendering of \"Eöt-Wash\" (\"E¨ ot-Wash\"). In addition, the quoted inverse-square-law bound λΦ ≲ 10^-4 m from Ref. [32] should be stated more precisely, since it depends on the assumed coupling strength α^2=1/3 used throughout the paper.","section":"Sec. 6"}],"recommendation":"minor_revision","confidential_remarks":"The reader's report and the stress-test note align with my own reading: the central sub-Chandrasekhar claim is robust, and the concern about the linearized-exterior Robin condition is not load-bearing for the maximum-mass sequences. The manuscript needs only local corrections, notably to the verification statement in Sec. 5.1 and the Fig. 5 caption inconsistency."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis one deserves a real referee slot. It is the first white-dwarf calculation in minimal dilatonic gravity, and the headline result--strictly sub-Chandrasekhar maximum mass for every dilaton Compton length--holds up under scrutiny. The numerics are unusually careful: two independent solvers agree to 1e-6-1e-4, the GR and Lane-Emden limits are recovered, and the analytic lambda -> infinity limit (G_eff = 4G/3, mass ratio 0.6495) is reproduced to 0.06%. The EOS sensitivity of the suppression ratio is genuinely small, which makes the lambda bounds credible.\n\nThe authors are also honest about the main observational weakness: the ZTF J1901+1458-based bound depends on composition and on an assumed mass-radius relation, so they wisely present the gravitational-redshift bound as the better-controlled one, even though it is weaker. They disclose that the code will be released after publication.\n\nSoft spots are minor. The big one is in Sec. 5.1: the text says |Phi_c - 1| is 1e-4 \"throughout,\" but the high-density lambda = 1500 km models in Table 4 have compactness ~1e-3, implying central dilaton deviations around 1e-3. The statement should be corrected. It is not load-bearing: the maximum-mass sequence lives at R*/lambda <= 5, where the shooting solver (which does not use the Robin condition) is valid and agrees with the collocation method; the deep-screening regime where the linearized-exterior assumption matters involves low-density stars, not the interesting maximum masses. The Earth-orbit bound is approximate, and the authors themselves note that Eot-Wash already excludes these lambda values, so the astrophysical window is likely closed from the laboratory side. That caps the significance but does not undermine the central result.\n\nCitation pattern looks clean: relevant modified-gravity white-dwarf literature is cited, the Fiziev foundation is acknowledged, and there is no obvious self-citation inflation. The free-boundary collocation method is a useful contribution in its own right, since it removes the exponential stiffness that has limited scalar-tensor stellar-structure calculations.\n\nRecommendation: send to peer review. Ask the author to fix the 1e-4 wording, to clarify the validity range of the Robin condition, and to make the code public before or at publication. The central claim is solid and the paper will be useful to anyone working on scalar-tensor constraints from compact objects.","headline":"Careful numerical work that closes the white-dwarf window for minimal dilatonic gravity; the main sub-Chandrasekhar claim holds, with one overstated verification check.","tokens_in":17124,"tokens_out":2004,"would_cite":true,"duration_ms":22882,"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":"Every white dwarf in minimal dilatonic gravity stays below the Chandrasekhar limit, regardless of the dilaton's Compton length.","keywords":["white dwarfs","minimal dilatonic gravity","Brans-Dicke theory","Chandrasekhar limit","scalar-tensor gravity","dilaton Compton length","mass-radius relation","free-boundary collocation"],"falsifier":"Measure a white dwarf's dynamical mass above the maximum that the same equation of state allows in general relativity, about 1.43 solar masses for the ideal Chandrasekhar EOS; unless differential rotation, which the paper does not model, is responsible, the strict sub-Chandrasekhar prediction would be falsified. A cleaner laboratory check is to compare two independent satellite geodesy determinations of $GM_\\oplus$ at different altitudes: if they agree to better than $10^{-3}$ for $\\lambda_\\Phi \\gtrsim 300$ km, the predicted percent-level altitude dependence is excluded.","tokens_in":16074,"feed_emoji":"⭐","tokens_out":10646,"duration_ms":94703,"temperature":0.7,"pith_summary":"Minimal dilatonic gravity (MDG) is a one-parameter scalar-tensor theory: a Brans-Dicke-like dilaton with fixed coupling $\\alpha^2=1/3$ and a free Compton length $\\lambda_\\Phi$. This paper solves the full relativistic stellar-structure equations for white dwarfs in MDG and finds that the maximum mass is always lower than in general relativity, strictly sub-Chandrasekhar for every $\\lambda_\\Phi$. The maximum drops from $1.425\\,M_\\odot$ in GR to $1.27\\,M_\\odot$ at $\\lambda_\\Phi=200$ km and $1.09\\,M_\\odot$ at 500 km, and the suppression ratio is almost independent of the equation of state. Because the most massive observed white dwarfs sit near the GR limit, the paper derives upper bounds $\\lambda_\\Phi \\lesssim 300$ km including rotation, and $\\lambda_\\Phi \\lesssim 720$ km from the gravitational-redshift mass-radius relation. A new free-boundary collocation method with a linearized-exterior Robin condition makes the screened regime computable in double precision; the same parameter predicts a percent-level altitude dependence of the Kepler-inferred $GM_\\oplus$, connecting white dwarfs to Earth-orbit tests.","feed_headline":"White dwarfs stay sub-Chandrasekhar in minimal dilatonic gravity","feed_subtitle":"Maximum mass falls from 1.425 to 1.09 solar masses, pinning the dilaton's range below roughly 300 km.","key_machinery":"The load-bearing machinery is the set of relativistic structure equations for a static, spherically symmetric star in MDG, with the Brans-Dicke coupling fixed to $\\omega=0$ (scalar coupling $\\alpha^2=1/3$) and a one-parameter 'withholding' potential whose curvature at its minimum sets the dilaton mass $m_\\Phi=\\lambda_\\Phi^{-1}$. The second essential piece is the new free-boundary collocation scheme: the exterior dilaton solution is linearized to a single decaying Yukawa mode and condensed into a Robin boundary condition $p_\\Phi(r_*)=\\Delta_*/(4\\pi r_*^2)(1/\\lambda_\\Phi+1/r_*)(\\Phi_*-1)$ at the stellar surface. This eliminates the $e^{R_*/\\lambda_\\Phi}$ amplification that makes shooting methods fail beyond $R_*/\\lambda_\\Phi\\simeq 20$, so the screened regime, where deviations from GR are exponentially small, can be computed in double precision up to $R_*/\\lambda_\\Phi\\simeq 35$ and continued analytically beyond. An analytic control is the in-matter equilibrium $\\Phi_{\\rm eq}=[1-q]^{-1/2}$ with $q=16\\pi(\\epsilon-3p)\\lambda_\\Phi^2/3$, which the deep-screening solutions approach as expected.","core_discovery":"The central claim is that no choice of the dilaton Compton length $\\lambda_\\Phi$ yields a super-Chandrasekhar white dwarf in MDG. Over the full central-density range, the mass-radius curve sits below the general-relativistic one, with the suppression growing from 3.1% at $\\lambda_\\Phi=500$ km to 18.9% at 1500 km, and the maximum stellar mass falling monotonically from $1.425\\,M_\\odot$ in GR to $1.268\\,M_\\odot$ at 200 km and $1.090\\,M_\\odot$ at 500 km. Part of the deficit is stored in the exterior dilaton field, the 'disphere': $m_{\\rm tot}$ includes this gravitating but non-baryonic component, and for large $\\lambda_\\Phi$ the disphere can hold about 20% of the total mass, though only a sub-percent fraction near the observationally allowed window. The suppression factor is nearly independent of the equation of state (a 0.2% difference between the ideal Chandrasekhar and Coulomb-corrected cases), so the derived bounds on $\\lambda_\\Phi$ are not microphysics artifacts. With rigid rotation near mass shedding adding only 4% to 6%, the paper concludes that consistency with the most massive known white dwarf requires $\\lambda_\\Phi \\lesssim 300$ km, and the gravitational-redshift test requires $\\lambda_\\Phi \\lesssim 720$ km.","pith_inferences":["I infer that the white-dwarf bound and the Earth-orbit signal are two windows on the same scale: a future null result in cross-mission $GM_\\oplus$ comparisons at the $10^{-3}$ level would close the astrophysically interesting $\\lambda_\\Phi$ window, while a positive signal would predict a specific 3% to 11% suppression of the most massive white dwarfs.","A consequence the author leaves implicit is that the linearized-exterior collocation method, if stable as described, should transfer to any massive scalar-tensor theory whose exterior is approximately Yukawa; a natural test is to apply it to chameleon or $f(R)$ white dwarfs and check where the single-mode Robin condition breaks.","Because the dilaton mass is density-independent, the theory's astrophysical window is in tension with local inverse-square-law experiments, which already require $\\lambda_\\Phi\\lesssim 10^{-4}$ m; the paper suggests that a density-dependent scalar mass would resolve this, and I infer that such a resolution would remove the strict sub-Chandrasekhar prediction as a sharp test."],"forward_implications":["No value of $\\lambda_\\Phi$ produces a super-Chandrasekhar branch; the maximum mass is monotonically suppressed as $\\lambda_\\Phi$ grows, so MDG cannot explain super-luminous Type Ia supernovae through massive progenitors.","The most massive observed white dwarfs, which sit near the GR limit, restrict $\\lambda_\\Phi$ to roughly 100 to 170 km in the non-rotating case and to $\\lesssim 300$ km once rigid rotation near mass shedding is allowed.","The population-level gravitational-redshift test, which constrains $m_*/R_*$, yields the independent bound $\\lambda_\\Phi \\lesssim 720$ km, corresponding to $m_\\Phi \\gtrsim 2.7\\times 10^{-13}$ eV/$c^2$, and at $2\\sigma$ the bound is $\\lambda_\\Phi \\lesssim 1100$ km.","At the same Compton lengths the Kepler-inferred $GM_\\oplus$ becomes altitude-dependent: two satellite geodesy determinations at different altitudes would disagree by $3.5\\times 10^{-5}$ at $\\lambda_\\Phi=100$ km, $2.9\\times 10^{-3}$ at 300 km, and $1.8\\times 10^{-2}$ at 700 km.","The suppression ratio is insensitive to the equation of state, with a 0.2% difference between ideal and Coulomb-corrected models, so the $\\lambda_\\Phi$ bounds transfer to white-dwarf compositions beyond those explicitly modeled."],"supporting_citations":[{"why":"It supplies the general-relativistic baseline maximum mass that MDG must suppress.","marker":"[7]"},{"why":"It defines the MDG action and the withholding potential whose curvature sets the dilaton mass.","marker":"[16]"},{"why":"It supplies the relativistic structure equations and dimensionless units used throughout the calculation.","marker":"[17]"},{"why":"It identifies the dilaton mass $m_\\Phi=\\lambda_\\Phi^{-1}$ and the Klein-Gordon behavior around the GR value.","marker":"[18]"},{"why":"It provides the observed white-dwarf mass and radius distribution used for the population context and maximum-mass comparison.","marker":"[2]"},{"why":"It provides the gravitational-redshift mass-radius test that yields the $\\lambda_\\Phi \\lesssim 720$ km bound.","marker":"[3]"},{"why":"It provides the observed mass of the most massive known white dwarf, used to derive the $\\lambda_\\Phi \\lesssim 300$ km constraint.","marker":"[4]"},{"why":"It establishes the Newtonian screening result that MDG extends to full relativity and that fixes the expected sign of the mass suppression.","marker":"[12]"},{"why":"It supplies the Coulomb-lattice correction used to test the equation-of-state sensitivity of the suppression ratio.","marker":"[21]"},{"why":"It supplies the local inverse-square-law bounds that create the tension with the astrophysical $\\lambda_\\Phi$ window.","marker":"[32]"}],"fun_headline_variants":["No super-Chandrasekhar white dwarfs in dilatonic gravity","Dilaton gravity caps white dwarf mass below Chandrasekhar","White dwarfs lose mass in minimal dilatonic gravity","Dilaton range constrained by white dwarf mass limits","Sub-Chandrasekhar limit holds in minimal dilatonic gravity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that just outside the star's surface the dilaton field is always close enough to its general-relativistic value that the exterior can be replaced by a single decaying Yukawa mode and a linearized Robin condition; the computed masses, the stellar-versus-disphere split, and the deep-screening results all depend on that approximation.","fun_headline_variants_meta":{"raw":{"variants":["No super-Chandrasekhar white dwarfs in dilatonic gravity","Dilaton gravity caps white dwarf mass below Chandrasekhar","White dwarfs lose mass in minimal dilatonic gravity","Dilaton range constrained by white dwarf mass limits","Sub-Chandrasekhar limit holds in minimal dilatonic gravity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00027,"raw_usage":{"total_tokens":1701,"prompt_tokens":1098,"completion_tokens":603,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":714,"completion_tokens_details":{"reasoning_tokens":517}},"tokens_in":714,"tokens_out":603,"duration_ms":5310,"temperature":1.0,"reasoning_tokens":517,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:41:35.880381+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure a white dwarf's dynamical mass above the maximum that the same equation of state allows in general relativity, about 1.43 solar masses for the ideal Chandrasekhar EOS; unless differential rotation, which the paper does not model, is responsible, the strict sub-Chandrasekhar prediction would be falsified. A cleaner laboratory check is to compare two independent satellite geodesy determinations of $GM_\\oplus$ at different altitudes: if they agree to better than $10^{-3}$ for $\\lambda_\\Phi \\gtrsim 300$ km, the predicted percent-level altitude dependence is excluded.","supporting_citations":[{"cited_title":"Chandrasekhar,The Maximum Mass of Ideal White Dwarfs, Astrophys","cited_arxiv_id":null,"evidence_quote":"It supplies the general-relativistic baseline maximum mass that MDG must suppress."},{"cited_title":"Withholding Potentials, Absence of Ghosts and Relationship between Minimal Dilatonic Gravity and f(R) Theories","cited_arxiv_id":"1209.2695","evidence_quote":"It defines the MDG action and the withholding potential whose curvature sets the dilaton mass."},{"cited_title":"Compact static stars in minimal dilatonic gravity","cited_arxiv_id":"1402.2813","evidence_quote":"It supplies the relativistic structure equations and dimensionless units used throughout the calculation."},{"cited_title":"The mass of dark scalar and phase space analysis of realistic models of static spherically symmetric objects","cited_arxiv_id":"1512.03931","evidence_quote":"It identifies the dilaton mass $m_\\Phi=\\lambda_\\Phi^{-1}$ and the Klein-Gordon behavior around the GR value."},{"cited_title":"Measuring The Mass-Radius Relation of White Dwarfs Using Wide Binaries","cited_arxiv_id":"2310.19866","evidence_quote":"It provides the gravitational-redshift mass-radius test that yields the $\\lambda_\\Phi \\lesssim 720$ km bound."},{"cited_title":"Screening Mechanisms on White Dwarfs: Symmetron & Dilaton","cited_arxiv_id":"2505.05871","evidence_quote":"It establishes the Newtonian screening result that MDG extends to full relativity and that fixes the expected sign of the mass suppression."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the Coulomb-lattice correction used to test the equation-of-state sensitivity of the suppression ratio."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the local inverse-square-law bounds that create the tension with the astrophysical $\\lambda_\\Phi$ window."}],"review_version":1}