{"id":"0a44295c-af04-47ca-848b-1c394b372f23","arxiv_id":"2608.12992","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Finite temperature makes theoretical white dwarf models in scalar-tensor gravity larger at the same mass, creating a degeneracy with modified gravity signatures.","lead":"This paper adds finite temperature to white dwarf models in three modified gravity theories and compares their mass-radius curves. It finds that heat puffs the stars up, and that this puffiness can look like modified gravity.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The high-temperature radius increase may be dominated by an arbitrary pressure cutoff in non-degenerate outer layers, making the central radius claim partly numerical.","rationale":"The reader's weakest-assumption identification is that the finite-temperature Chandrasekhar EoS is pushed beyond validity in the non-degenerate outer layers and at T = 10^8 K. My stress-test sharpens this into a concrete numerical mechanism: the outer layers behave like an isothermal ideal gas, so the stellar radius is not defined by P = 0 but by the arbitrary cutoff P = 10^-10 P_c, and the high-T radius excess may be dominated by that cutoff rather than by genuine stellar structure. This is the most load-bearing place because the paper's novelty is the temperature-induced radius increase, and the largest such increase is reported precisely in the regime where the EoS and the surface definition are weakest. The paper is transparent about EoS validity limits, which is why I do not recommend moving to REJECT or UNVERDICTED; the concern is addressable with a targeted recomputation. The dilaton parameter conversion issue and absence of released code are real weaknesses, but they do not threaten the central temperature effect as directly as the surface/cutoff problem. The recommended verdict is unchanged: CONDITIONAL, with the additional explicit condition that the T = 10^8 K radius predictions must be shown to be insensitive to the pressure cutoff and to the non-degenerate outer-layer treatment.","tokens_in":17257,"tokens_out":6870,"duration_ms":82146,"concrete_test":"Recompute the GR and massive Brans-Dicke mass-radius branches at T = 10^8 K with the surface pressure cutoff changed from P(r_s) = 10^-10 P_c to P(r_s) = 10^-14 P_c, and also with a matched outer-layer EoS that includes arbitrary-degeneracy electron pressure, ion pressure, and a non-isothermal temperature profile. If the radii shift by more than a few percent, especially at low masses where the T = 10^8 K curves converge, then the central high-temperature radius claim is an artifact of the cutoff and the CONDITIONAL verdict should remain with that caveat made explicit; if radii are stable, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central observable claim is that finite temperature expands the white dwarf radius at fixed mass, most strongly at T = 10^8 K, where the paper also reports convergence of all modified-gravity curves to a common large radius. Section 4 defines the stellar surface as P(r_s) = 10^-10 P_c and asserts that lowering the cutoff does not change radii, while Secs. 5.1 and 6 concede that the finite-temperature Chandrasekhar EoS loses validity as the gas becomes less degenerate and that realistic models need non-degenerate outer layers. The problem is that in exactly the low-density outer layers relevant at T = 10^8 K, the finite-temperature EoS approaches an isothermal ideal-gas form with P proportional to rho. In hydrostatic equilibrium such an atmosphere has no natural surface: the radius at which P/P_c = 10^-10 moves with the arbitrary pressure cutoff, roughly logarithmically in P_c/P_cut. The T = 10^8 K branches, including the claimed common large-radius limit, occupy precisely this low-density, weakly degenerate regime, so the reported radii may be substantially an artifact of the P = 10^-10 P_c convention rather than a robust physical prediction. Because the mass is almost unchanged by construction under the rest-mass-only energy-density approximation, the radius is the only diagnostic carrying the temperature signal, which makes this sensitivity especially concerning.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript studies the equilibrium structure of white dwarfs at finite temperature in three scalar-tensor gravity theories: massive Brans-Dicke theory, symmetron screening, and dilaton screening. The authors derive a finite-temperature Chandrasekhar equation of state, integrate the spherically symmetric Einstein-frame field equations with a shooting method, and present mass-radius relations, effective gravitational constants, and radial profiles for temperatures up to 10^8 K. The central claim is that a non-zero temperature increases the white dwarf radius while leaving the total mass essentially unchanged, and that this thermal effect can partially mimic or mask modified-gravity signatures. The paper also reports that the T = 10^8 K branches converge to a common large radius across all theories, while repeatedly cautioning that the equation of state loses validity in precisely that regime.","tokens_in":17554,"tokens_out":6948,"duration_ms":75100,"significance":"If the central claim is robust, the paper provides a useful investigation of an under-explored degeneracy between thermal effects and scalar-tensor gravity modifications in white dwarfs, with potential implications for white dwarf cooling models and for interpreting mass-radius observations. The manuscript is transparent about its numerical setup, benchmarks the GR limit against Ref. [19], and explicitly flags the EoS validity limits and an unresolved dilaton parameter conversion. The main reservation is that the radius, the only temperature-sensitive diagnostic in the presented results, is defined by a pressure cutoff in a regime where the EoS is admitted to become invalid; the robustness of the central claim therefore depends on a quantitative treatment of the outer layers and of the thermal energy density that is not currently provided.","major_comments":[{"comment":"The stellar radius is defined numerically by the pressure cutoff P(r_s) = 10^-10 P_c, with the assertion that demanding more orders of magnitude does not significantly change the radii. At T = 10^8 K the low-density outer layers become non-degenerate, as the authors themselves state in Secs. 5.1 and 6; in this regime the finite-temperature EoS approaches an isothermal ideal-gas form P ∝ ρ, for which hydrostatic equilibrium has no natural surface and r_s shifts with the arbitrary cutoff. Since the T = 10^8 K branches and their convergence to a common large radius are presented as a main result, and since the mass is nearly unchanged by construction, the radius is the only temperature-sensitive diagnostic. Please quantify the cutoff dependence of r_s at T = 10^8 K and, if possible, at T = 10^7 K, or restrict the claimed radius increase to densities and temperatures where the Chandrasekhar EoS remains valid.","section":"Section 4 and Section 5.1"},{"comment":"The energy density is approximated as rest-mass only, described as 'safely approximated to consist only of this term', but no estimate of the thermal energy contribution is given. The full energy density in Eq. (1b) contains a temperature-dependent part that becomes relatively important at low degeneracy and high temperature. Using only the rest-mass density in the hydrostatic equilibrium equation (15) removes a priori any direct temperature dependence of the energy density that enters the mass integral, so the conclusion that the mass is 'essentially unchanged' is partly an artifact of this approximation. Please quantify the thermal energy density relative to the rest-mass energy density over the full range of densities and temperatures used (especially ρ_c ~ 10^5-10^7 g cm^-3 at T = 10^8 K) and show the corresponding change in the inferred stellar mass.","section":"Section 2.1 and Eq. (15)"},{"comment":"The dilaton results rely on parameter values that were adopted after an unresolved factor-of-10^-3 discrepancy with the theoretical conversion from Ref. [42]. The authors state that the reported mass-radius relations are recovered only for parameters smaller by this factor than those obtained from the theoretical conversion, and that the source of the discrepancy was not identified. Since the dilaton mass-radius curves and their temperature shifts are part of the quantitative results, this unresolved conversion should be resolved or its effect on the conclusions quantified; as written, the dilaton predictions rest on an empirically adjusted parameter set.","section":"Section 5.1, footnote 1"}],"minor_comments":[{"comment":"The abstract and conclusions state the finite-temperature radius increase without the important caveat, emphasized in Secs. 5.1 and 6, that at T = 10^8 K the equation of state is reaching its validity limit and that realistic models require non-degenerate outer layers. Please qualify the main claim accordingly.","section":"Abstract and Section 6"},{"comment":"The statement that lowering the pressure cutoff does not significantly change the radii is not supported by a convergence test. A small figure or table showing r_s as a function of P_cut/P_c for representative models, including at least one T = 10^8 K case, would make the numerical boundary treatment reproducible.","section":"Section 4"},{"comment":"There are several typographical errors: Sec. 5.2 has 'it's value' instead of 'its value'; Sec. 6 has 'grativational' instead of 'gravitational'; Appendix B has 'mehtod' instead of 'method'; Sec. 5.1 has 'looses' instead of 'loses'. A careful proofreading pass is needed.","section":"Throughout"},{"comment":"The symmetron parameters are adopted from Ref. [42] and are outside the observational bounds quoted from Ref. [44]. The paper states this is intentional for comparison and mechanism illustration, which is acceptable, but the conclusions should explicitly remind the reader that the quantitative mass-radius shifts for the symmetron are not predictions for observationally allowed parameters.","section":"Section 5.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and the authors are transparent about their limitations, including the unresolved dilaton parameter conversion and the EoS validity boundary. The main load-bearing issue is the pressure-cutoff definition of the stellar surface at high temperature, which directly affects the paper's primary radius claim; this needs either a quantitative demonstration of cutoff independence or a restriction of the conclusions to the degenerate regime. The thermal energy density approximation should also be quantified before the 'mass essentially unchanged' statement can be trusted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does something useful and honest: it takes the known finite-temperature Chandrasekhar EoS and runs it through three scalar-tensor gravity models that were previously only treated at zero temperature. The new content is the specific mass-radius relations, effective gravitational constant curves, and radial profiles at T up to 1e8 K, and the paper openly discusses the degeneracy this creates between thermal and modified-gravity signatures. The GR benchmark matches [19], the EoS formulas are standard, and the authors flag where the EoS loses validity. That is solid forward modeling, not a hidden fit.\n\nThe main soft spot is the surface definition. The radius is set by P(r_s) = 1e-10 P_c, and the paper asserts that lowering the cutoff does not change radii. But at T = 1e8 K, the outer layers are weakly degenerate and the EoS approaches an isothermal ideal-gas form with P ~ rho. In hydrostatic equilibrium, that means the radius at which P/P_c = 1e-10 shifts with the cutoff, roughly logarithmically in P_c/P_cut. The T = 1e8 K branches, including the claimed common large-radius limit, live exactly in this regime. So the central radius claim at high T is partly convention-dependent, and the convergence of all modified-gravity curves to one large radius may be more numerical than physical. The authors do warn that the EoS is losing validity there, but they still present those radii as the main temperature effect, so this is not a minor caveat.\n\nOther issues: no convergence tests or error bars, the rest-mass-only energy-density approximation is never quantified, the dilaton parameter-conversion discrepancy is reported but left unresolved, and no code or data is released. These are addressable and the paper is transparent about the dilaton issue.\n\nThe stress-test note about the surface cutoff is the one that matters. If the authors rerun the high-T cases with the pressure cutoff varied, or use a realistic non-degenerate envelope, I would expect the T = 1e8 K radii to shift noticeably. The lower-temperature results, where the effect is small anyway, are probably robust.\n\nBottom line: the paper is worth a serious referee, but the referee should push on the surface-convention dependence before accepting the headline radius increase at high T. For temperatures up to ~1e6 K, the result is solid and the degeneracy point is well taken. For 1e8 K, the claims need a clearer numerical check. I would want to see that check before citing the large-radius branches.","headline":"A transparent finite-temperature WD study in scalar-tensor theories with a plausible central result but a surface-definition soft spot that needs addressing before the high-T radius claims can be trusted.","tokens_in":18038,"tokens_out":631,"would_cite":false,"duration_ms":7876,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.50.Kd","04.40.Dg","97.20.Rp"],"model":"deepseek-v4-flash","headline":"The paper shows that a white dwarf at finite temperature is larger at the same mass in three modified-gravity theories, and that at 10^8 K thermal effects can mimic or mask changes in the gravity parameters.","keywords":["white dwarfs","finite temperature equation of state","mass-radius relation","massive Brans-Dicke theory","symmetron screening","dilaton screening","scalar-tensor gravity","effective gravitational constant"],"falsifier":"Rebuild one of the paper's $10^{8}$ K models, for example the low-mass, large-radius branch near 0.75 solar masses, using a more complete finite-temperature treatment that includes non-degenerate electrons, ions, and a radiative outer envelope; if the radius at the same mass is no longer larger than the $10^{4}$ K model, or the branch disappears, the central claim is falsified.","tokens_in":17058,"feed_emoji":"🌡️","tokens_out":6444,"duration_ms":65920,"temperature":0.7,"pith_summary":"The paper asks whether a white dwarf's internal temperature changes the equilibrium size and mass predicted by modified-gravity theories, and whether temperature could be mistaken for a genuine gravitational effect. It computes white-dwarf models with a finite-temperature version of the Chandrasekhar equation of state in three scalar-tensor theories: massive Brans-Dicke, symmetron, and dilaton. The central result is that heating the star leaves its total mass almost unchanged but inflates its radius; below roughly $10^{6}$ K the effect is negligible, and above that it grows sharply. Because the high-temperature mass-radius branches of different theories overlap, the paper concludes that finite temperature creates partial degeneracies that must be accounted for when using white-dwarf observations to constrain gravity.","feed_headline":"Heat inflates white dwarfs in modified gravity","feed_subtitle":"Finite-temperature models give larger radii at fixed mass and blur differences among gravity theories.","key_machinery":"The central object is the finite-temperature Chandrasekhar equation of state, built from relativistic Fermi-Dirac integrals $F_k(\\eta,\\beta)$ with $\\beta=k_B T/(m_e c^2)$ and degeneracy parameter $\\eta=\\mu/(k_B T)$, which replaces the zero-temperature step-function occupancy of the degenerate electron gas. This EOS is coupled to the Einstein-frame scalar-tensor field equations and the hydrostatic equilibrium equation $P' = -(\\rho c^2 + P)(\\tilde{\\varphi}' + \\alpha \\tilde{\\phi}')$, with the scalar field governed by an effective potential set by each theory's conformal factor $A(\\tilde{\\phi})$ and self-interaction $V(\\tilde{\\phi})$. The stellar surface is fixed by $P=0$, and the effective gravitational constant is $G_{\\mathrm{eff}} = A^2(\\tilde{\\phi}) G_N$. This machinery carries the argument because it is the temperature-dependent pressure that expands the star while the scalar field modifies the force balance.","core_discovery":"On its own terms, the paper establishes that in all three theories a nonzero temperature shifts the mass-radius relation: at fixed total mass the equilibrium radius increases with temperature while the mass stays essentially unchanged, because thermal pressure adds to degeneracy pressure. Up to about $10^{6}$ K, deviations from zero temperature are very small; at $10^{7}$ K and $10^{8}$ K they become significant. In every theory the effective gravitational constant is weaker at the stellar center than at infinity, with the largest deviations in massive Brans-Dicke and much smaller deviations in the screened theories. At $10^{8}$ K the mass-radius curves of all theories converge at large radii and low masses, which the paper reads as a hint that thermal effects dominate modified-gravity effects in that regime, while noting that the equation of state is starting to lose validity there. The interior radial profiles of the scalar field, pressure, and metric are mostly temperature independent except near $10^{8}$ K.","pith_inferences":["Extending the paper's logic, combining mass-radius data with independent temperature or cooling-rate measurements could lift the degeneracy the paper identifies, because temperature and modified-gravity parameters affect cooling ages differently.","The high-temperature convergence of all theories at low masses suggests that any modified-gravity interpretation of an unusually large or small white dwarf should first check whether a hot standard model can reproduce the same radius; that check is not performed in the paper.","A natural next step would be to feed the same finite-temperature EOS into rotating or magnetized white-dwarf models, where thermal expansion could shift the stability boundary relevant to super-Chandrasekhar candidates.","The paper's own validity caveat implies the 10^8 K results are provisional: replacing the outer layers with a non-degenerate atmosphere and adding neutrino emission or shell burning could either sharpen or erase the apparent degeneracies."],"forward_implications":["For all three theories, a white dwarf of fixed mass has a larger radius at higher temperature, with negligible change below about 10^6 K and significant change above.","White-dwarf cooling tracks from 10^8 K down to 10^6 K should not assume a constant radius, because the thermal expansion affects the relation between cooling time and observed size.","Finite temperature produces partial overlaps between different theory-parameter curves in the mass-radius plane, so parameter constraints drawn from observations must include temperature or risk being biased.","At 10^8 K the mass-radius curves for massive Brans-Dicke, symmetron, and dilaton converge at large radii and low masses, suggesting that temperature effects dominate modified-gravity effects there, subject to the paper's own caveat about equation-of-state validity.","In all three theories the central effective gravitational constant is smaller than its value at infinity, and temperature leaves this mostly unchanged except for a saturation effect at 10^8 K that the paper attributes to a possible numerical or EOS artifact."],"supporting_citations":[{"why":"Supplies the zero-temperature Chandrasekhar equation of state and the maximum white-dwarf mass that the finite-temperature version generalizes.","marker":"[5]"},{"why":"Provides the finite-temperature white-dwarf equilibrium structure in general relativity that the paper uses as its GR baseline.","marker":"[19]"},{"why":"Supplies the symmetron and dilaton parameter values and white-dwarf screening setup that the paper reproduces and extends.","marker":"[42]"},{"why":"Provides the numerical shooting-method setup and prior scalar-tensor stellar models that the code is built on.","marker":"[43]"},{"why":"Justifies the chosen massive Brans-Dicke scalar-field mass range through observational constraints.","marker":"[41]"},{"why":"Defines the symmetron screening mechanism that fixes the conformal factor and potential used in the model.","marker":"[39]"},{"why":"Defines the dilaton screening mechanism that fixes the conformal factor and runaway potential used in the model.","marker":"[40]"}],"fun_headline_variants":["Hot white dwarfs swell in modified gravity","Thermal pressure inflates white dwarf radii in modified gravity","Finite temperature widens white dwarf radii in modified gravity","Heat can dominate modified gravity in hot white dwarfs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the finite-temperature Chandrasekhar equation of state, with energy density dominated by rest mass and electrons treated as a degenerate Fermi gas, remains valid through the entire star up to $10^{8}$ K, including the outer layers where the paper concedes the gas becomes less degenerate and the EOS starts losing validity.","fun_headline_variants_meta":{"raw":{"variants":["Hot white dwarfs swell in modified gravity","Thermal pressure inflates white dwarf radii in modified gravity","Finite temperature widens white dwarf radii in modified gravity","Heat can dominate modified gravity in hot white dwarfs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000649,"raw_usage":{"total_tokens":2913,"prompt_tokens":811,"completion_tokens":2102,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":427,"completion_tokens_details":{"reasoning_tokens":2051}},"tokens_in":427,"tokens_out":2102,"duration_ms":15505,"temperature":1.0,"reasoning_tokens":2051,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:01:45.600509+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Rebuild one of the paper's $10^{8}$ K models, for example the low-mass, large-radius branch near 0.75 solar masses, using a more complete finite-temperature treatment that includes non-degenerate electrons, ions, and a radiative outer envelope; if the radius at the same mass is no longer larger than the $10^{4}$ K model, or the branch disappears, the central claim is falsified.","supporting_citations":[{"cited_title":"massive Brans-Dicke theory, symmetron and dilaton screening, at different finite temperatures","cited_arxiv_id":null,"evidence_quote":"Supplies the zero-temperature Chandrasekhar equation of state and the maximum white-dwarf mass that the finite-temperature version generalizes."},{"cited_title":"While the field is driven to zero in high density re- gions, symmetry breaking causes it to converge towards the potential minimum at low densities","cited_arxiv_id":null,"evidence_quote":"Provides the finite-temperature white-dwarf equilibrium structure in general relativity that the paper uses as its GR baseline."},{"cited_title":"Equilibrium structure of white dwarfs at finite temperatures","cited_arxiv_id":"1510.02024","evidence_quote":"Supplies the symmetron and dilaton parameter values and white-dwarf screening setup that the paper reproduces and extends."},{"cited_title":"Mach’s principle and a relativistic theory of gravitation,","cited_arxiv_id":null,"evidence_quote":"Provides the numerical shooting-method setup and prior scalar-tensor stellar models that the code is built on."},{"cited_title":"Stability criterion for white dwarfs in Palatini $f(R)$ gravity","cited_arxiv_id":"2111.08029","evidence_quote":"Defines the symmetron screening mechanism that fixes the conformal factor and potential used in the model."}],"review_version":1}