{"id":"1d2a7e01-76ab-4088-a70b-76644f5b21e8","arxiv_id":"2412.13867","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"For neutron stars near the critical compactness, scalar-tensor gravity can change computed X-ray pulse fluxes by up to 80% and inferred radii by about 10%.","lead":"This paper computes X-ray pulse shapes from hot spots on very compact neutron stars in scalar-tensor gravity, including light bending, time delays, and multiple images. It finds up to 80% flux differences from general relativity for stars near the critical compactness, suggesting pulse profiles could test modified gravity.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 80% flux claim rests on a two-image prescription whose geometric validity at the critical compactness is asserted rather than verified from the exact integral for the secondary image.","rationale":"The reader's weakest_assumption focuses on atmospheric beaming, rotation, oblateness, and the idealized spherical setup. Those are real limitations for astrophysical applicability, and the paper itself is transparent about them in Sections 3 and 5, so they are appropriate for a CONDITIONAL verdict. However, my reading of the manuscript places the more load-bearing concern one step earlier: before asking whether the idealized model describes a real pulsar, one should ask whether the idealized model's light curve is correctly computed in the high-compactness regime. The two-image prescription (center psi_sec = 2*pi - psi_0, turned on by the criterion cos(psi_0 + Delta_psi) <= cos(psi_c), using the same flux kernel for both images) is the mechanism that generates the 80% difference claimed in the abstract and the strongest_claim. The paper does not show a derivation of this prescription from the Just-metric geodesics, nor a validation against direct ray tracing; it cites the analogous GR treatment (Sotani and Miyamoto 2018 for the threshold and Sotani 2020 for spot-size effects), but in the STT case psi_c depends on Q and A_s in a nonlinear way (Eq. 26, Fig. 1), and the Jacobian dcos(alpha)/dcos(psi) diverges near psi=pi, making the numerical handling of the secondary image delicate. This is an internal-correctness risk, not merely a disagreement with consensus, and it is the type of concern that a numerical cross-check would settle. That said, I do not see a reason to reject: the benchmark below the critical compactness agrees with Hu et al. (Appendix B) and the GR limit is reproduced. So the appropriate verdict remains CONDITIONAL, but conditioned not only on astrophysical realism (atmosphere, rotation, oblateness) but also on verification of the multiple-image flux construction. This is why I record partial agreement with the reader: we both land on CONDITIONAL, but my load-bearing concern is shifted from model realism to the geometric/numerical correctness of the secondary image in the very regime where the headline effect lives.","tokens_in":23033,"tokens_out":1929,"duration_ms":16028,"concrete_test":"Recompute the secondary-image contribution by direct geodesic ray tracing in the Just metric: for the Table 1 models, sample photon initial directions from the full stellar surface, integrate the null geodesic equations (Eqs. 21-24) to infinity, and accumulate flux by impact-parameter binning without assuming psi_sec = 2*pi - psi_0 or using the dalpha/dpsi Jacobian. Compare the resulting light-curve peak height and phase with Figs. 4 and B9.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central result, the up-to-80% flux difference in Figs. 4, 5, and B9, is produced entirely by the secondary-image brightening when psi_c exceeds pi: the spot crossing the region behind the star is assumed to contribute a second image centered at psi_sec = 2*pi - psi_0, with the same flux formula (Eq. 29) applied to that image, starting when cos(psi_0 + Delta_psi) <= cos(psi_c) (Section 3, after Eq. 31). This is analogous to how a Schwarzschild critical-compactness star is treated in Sotani and Miyamoto and in de Lima et al., but those works verified the two-image construction against direct photon geodesic counts. Here, the identification of the secondary image center and the applicability of the same differential-flux formula (which involves dcos(alpha)/dcos(psi) and a Jacobian that is singular at psi=pi) are asserted, not derived from the Just metric geodesics. At the near-critical values used (M/R 0.275-0.305, with psi_c near pi; Fig. 1), a small error in where or how the secondary image turns on would directly set the height and phase of the brightening peak that is the quantitative headline. The paper's own text concedes the setup is too simplified for real sources and flags atmosphere, oblateness, and rotation as needed additions, so the concern is about the demonstration of the 80% magnitude itself, not about realism alone.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"Ottoni et al. model X-ray pulsed light curves from a single finite hot spot on a neutron star in massless scalar-tensor theory, using the analytic Just exterior metric and the Silva-Yunes differential flux formula. They focus on compactness near the GR critical value M/R=0.284, where the whole stellar surface becomes visible and multiple images of the spot can appear. For doppelgänger models with the same mass and radius but scalar charge Q=0.5 and surface conformal factor As=0.95 or 1.05, they report bolometric flux differences relative to GR of up to 80%, driven by a secondary-image brightening when the spot crosses the lensed region behind the star, and a nonlinear shift of the critical view angle with compactness. They also estimate that STT can change the inferred neutron-star radius by up to about 10% and argue that improved radius measurements from NICER and future missions could eventually constrain scalar charges.","tokens_in":23313,"tokens_out":18317,"duration_ms":165807,"significance":"The paper identifies a promising qualitative effect: near the critical compactness, the pulse profile becomes sensitive to the scalar charge through the appearance and brightness of a secondary image. Its strengths are the use of an exact analytic exterior solution, the recovery of the GR limit, the agreement with Hu et al. below the critical compactness within 5% (Fig. B7), and a clear statement of the model's limitations. The nonlinear dependence of the critical view angle on compactness in Fig. 1 is also interesting. However, the quantitative 80% claim rests on a two-image construction whose validity for the Just metric is asserted rather than demonstrated, and on parameter combinations whose physical realizability is not established. The paper is therefore a useful exploratory study, but the headline effect is not yet anchored to a consistently derived observable.","major_comments":[{"comment":"The secondary-image prescription is asserted rather than derived for the Just metric. The paper places the secondary image center at ψ_sec = 2π − ψ_0, activates it when cos(ψ_0 + Δψ) ≤ cos(ψ_c), and applies the same integrand dα/dψ to it. This is the established Schwarzschild treatment in Refs. [17,23,62], but for the Just metric with scalar charge the equivalence of this prescription with the true lensed secondary image, including the correct Jacobian on the secondary branch, is not shown. Since the up-to-80% brightening in Figs. 4 and B9 is produced entirely by this secondary-image contribution near ψ_c, an error in the turn-on location or in the branch of dα/dψ would directly set the headline number. I request either a validation by backward ray tracing in the Just metric for representative parameters, or an analytic derivation of the secondary-branch Jacobian, or an explicit statement that the two-image construction is an assumption that has not been verified against direct photon geodesic counts.","section":"§3, Eq. (29) and the turn-on condition after Eq. (31)"},{"comment":"The doppelgänger models treat Q and As as independent inputs for fixed M and R. In a consistent STT, Q and As are not free parameters; they are fixed by the interior stellar solution and the coupling function. The paper does not show that (Q=0.5, As=0.95) or (Q=0.5, As=1.05) with M=2.1 M⊙ and R=10.918 km corresponds to any scalarized stellar model. Figure 2 demonstrates that high-compactness branches exist for the ENG and MPA1 equations of state with ξ=−3 and 25, but it does not report the values of Q and As along those branches at M/R=0.284. Without such a consistency check, the 80% flux difference is an upper envelope over a parameter space that may be unphysical. I recommend computing Q and As self-consistently for representative couplings and re-evaluating the lightcurve comparison with those values, or at least restricting the lightcurve survey to parameter combinations that are known to occur on a solution sequence.","section":"§3, Table 1"},{"comment":"The 80% maximum difference is produced by the secondary-image brightening, which is dominated by near-tangential emission (α ≈ π/2). The paper correctly notes in the conclusion that atmospheric effects can attenuate tangentially emitted photons, but because this is the mechanism of the headline effect, the observational claim should be presented as an upper limit under isotropic, atmosphere-free emission. A simple sensitivity test with a limb-darkening or absorption model would show how much of the 80% survives. Without it, the title-level claim that X-ray pulse profiles can 'probe' STT remains unquantified, even though the paper is upfront about the need for atmospheres in future work.","section":"§4, Figs. 4 and B9; §5"}],"minor_comments":[{"comment":"Please specify the integration domain over the spot in terms of ψ and φ, and define the spot's local azimuthal coordinate; as written, dψ dφ leaves the orientation of the spot relative to the rotation axis ambiguous for a finite circular spot.","section":"§3, Eq. (29)"},{"comment":"Define ¯a_s and ¯b_s in the caption (as a/ρ_s and b/ρ_s) and state explicitly that ρ_s is obtained by inverting Eq. (15) with the Jordan-frame radius R; otherwise the reader cannot verify the entries.","section":"Table 1"},{"comment":"State clearly that the left boundary of the shaded band corresponds to A_s=1.05 and the right boundary to A_s=0.95; the current wording is easy to misread.","section":"Figure 1 caption"},{"comment":"The phrase 'half the rational phase' should read 'half the rotational phase'.","section":"§4, second paragraph"},{"comment":"Clarify that the effective compactness values 0.234 and 0.258 are fitted to make the approximate curve follow the numerical one, and that the reported 7% error applies at α = π/2.","section":"Appendix A"},{"comment":"The paper does not report a numerical convergence test for the spot integration, which is relevant near the critical angle where dα/dψ diverges; please add a brief convergence statement, including the number of grid points in ψ and φ.","section":"General numerical methods"}],"recommendation":"major_revision","confidential_remarks":"This is an honest exploratory paper, but the headline 80% effect needs to be anchored by two things: a validation of the secondary-image construction for the Just metric, and a demonstration that the chosen (Q, A_s) pairs are realizable in some consistent STT model. The atmosphere issue is a limitation the authors already acknowledge, so it is a matter of framing and sensitivity rather than correctness. If the validation is added, the paper would be a reasonable contribution to the pulse-profile modeling literature. The scope fits a theory-oriented journal, though the observational claims should be tempered until the model includes atmospheres and self-consistent stellar structure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look if you work on pulse-profile tests of gravity. The genuinely new thing here is the first finite-spot X-ray pulse-profile calculation in scalar-tensor theory at compactnesses at and above the GR critical value (M/R ≈ 0.284), where the spot's secondary image becomes visible. Previous STT pulse papers used infinitesimal spots or stayed below critical compactness; previous high-compactness multiple-image work was done in GR. This paper sits at the intersection and it's a natural, useful step. The execution is mostly careful. They use the published Silva-Yunes flux formulas and the Just exterior metric as inputs, check the GR limit, and their below-critical comparison with Hu et al. (2021) agrees within 5% (Fig. B7). They are also refreshingly candid about the limitations: no atmosphere, no oblateness, spherical spacetime even at 700 Hz, and they say outright the model is too simplified for real sources. The qualitative claim—that STT shifts the critical angle so a spot crossing the far side of the star gets eclipsed or produces a lensed brightening, with flux differences that can be large—is plausible and worth taking seriously. The soft spot is the secondary-image construction, and it is not a minor detail. The paper places the secondary image at ψ_sec = 2π − ψ_0, switches it on when cos(ψ_0+Δψ) ≤ cos(ψ_c), and applies the same differential flux formula, which involves a Jacobian that diverges at ψ=π. That recipe is exactly the one Sotani and Miyamoto used for GR, but they verified it against direct photon geodesic counting. Here the Just-metric geodesics are never used to validate the secondary image position or flux. Since the 80% brightening in Figs. 4 and B9 is entirely produced by this secondary image, the quantitative headline is only as strong as that assertion. A referee should ask for a geodesic ray-tracing check or an independent numerical integration in the Just metric. The remaining concerns are milder. Q=0.5 and A_s=0.95/1.05 are chosen by hand rather than derived from a scalarized stellar model, so the size of the effect is parameter-dependent; there is no code, no error budget, and no statistical forecast against NICER data. Those are addressable and don't undermine the qualitative conclusion. The paper also correctly notes the degeneracy between compactness and scalar charge, so an independent radius measurement would be needed. Bottom line: this is a solid proof-of-concept that points to a promising observable for high-mass pulsars. The authors know the model is simplified. I'd send it to a serious referee, with the explicit request to verify the secondary-image prescription. I'd cite it if I worked in this area, and I'd bring it to a reading group to debate how confident we should be in the 80%.","headline":"First finite-spot STT pulse profiles at critical compactness; the 80% flux claim is a plausible but unverified secondary-image effect that deserves a hard look.","tokens_in":23904,"tokens_out":3798,"would_cite":true,"duration_ms":34917,"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":"The paper argues that near the critical general-relativistic compactness $M/R = 0.284$, a single hot spot's X-ray pulse becomes a sensitive probe of scalar-tensor gravity because the critical deflection angle shifts with the scalar charge…","keywords":["neutron stars","scalar-tensor gravity","X-ray pulse profiles","light bending","hot spots","spontaneous scalarization","Just metric","compactness"],"falsifier":"One concrete check is to recompute the STT pulse profile with a realistic atmospheric radiative-transfer model and an oblate rotating star: if the secondary-image contribution at half rotation phase drops below about 10% of the direct flux, the claimed 80% contrast cannot appear in real data. Observationally, a high-precision phase-resolved observation of a high-mass pulsar with an independently measured mass and radius that is fit well by a GR model at or above $M/R = 0.284$ without any extra brightening would rule out the paper's fiducial scalar charge.","tokens_in":22791,"feed_emoji":"🌟","tokens_out":12813,"duration_ms":109799,"temperature":0.7,"pith_summary":"The paper tries to establish that X-ray pulse profiles of neutron stars near the maximum-compactness regime—about $GM/(Rc^2) = 0.284$, where general relativity predicts the whole stellar surface becomes visible—are unusually sensitive to scalar-tensor gravity. Computing the flux of a single finite hot spot in the exact scalar-tensor exterior spacetime (the Just metric), with light bending, Shapiro time delay, and Doppler effects included, the authors find that a scalar charge $Q=0.5$ with a surface conformal factor $A_s=1.05$ can make the bolometric flux differ by up to 80% from the GR prediction when the spot crosses the far side of the star. The physical reason is that the critical deflection angle $\\psi_c$, which decides whether an invisible shadow zone or a second lensed image forms behind the star, shifts nonlinearly with compactness in scalar-tensor theory. If this is right, accurate X-ray pulse shapes of massive pulsars could constrain the scalar charge and scalar-tensor parameters, and radius inferences from light curves could be biased by up to about 10% if scalarization is ignored. The paper also argues that current X-ray radius uncertainties are too large to distinguish the theories, but next-generation X-ray and gravitational-wave measurements could close the gap.","feed_headline":"Scalar gravity can alter neutron-star X-ray pulses by 80%","feed_subtitle":"For the densest pulsars, a lensed second image makes pulse shapes a practical test of general relativity.","key_machinery":"The machinery is the Einstein-frame Just metric, the exact spherical exterior solution of scalar-tensor gravity, parametrized by the ADM mass $b=2M$, a length scale $a$, and the scalar charge $q$, with $a/b = \\sqrt{1+Q^2}$; the surface conformal factor $A_s$ converts the Jordan-frame stellar radius used by observers into the Einstein-frame radial coordinate $\\rho_s$ used in ray tracing. Photon paths are null geodesics of this metric, and the key derived object is the integral for the cumulative deflection angle $\\psi(\\alpha; \\bar a_s, Q)$ (Eq. 26) together with its value at tangential emission, the critical angle $\\psi_c = \\psi(\\pi/2)$. The flux is obtained by integrating the scalar-tensor differential flux formula over a finite circular spot, with the integration switching to the secondary image at $\\cos(\\psi_0+\\Delta\\psi) \\le \\cos\\psi_c$. The critical angle is the switch that turns a small STT-GR difference in geometry into a large difference in observable flux, which is why the paper identifies high compactness as the regime where pulse profiles can probe scalarization.","core_discovery":"The central claim is that the pulse profile of a highly compact neutron star is a sharp scalar-charge detector because the photon deflection integral (Eq. 26) makes the critical view angle $\\psi_c$ depend nonlinearly on both compactness and the scalar charge $Q$, through the ratio $a/b = \\sqrt{1+Q^2}$ of the Just exterior metric. At the GR threshold $M/R = 0.284$ the tangent ray bends through exactly $\\pi$, so the whole surface is visible; in scalar-tensor theory the same compactness can give $\\psi_c < \\pi$, leaving an invisible zone that eclipses the spot, or $\\psi_c > \\pi$, creating a second, lensed image whose flux adds a brightening near half the rotation period. For the paper's representative models ($Q=0.5$ and $A_s = 0.95$ or 1.05, same mass and Jordan-frame radius as a $2.1\\,M_\\odot$ GR star), that secondary-image brightening produces bolometric flux differences up to 80% relative to GR, and the relative change grows nonlinearly as compactness approaches the critical value while staying below roughly 5% for compactness well below it. The paper presents these results as an upper-limit estimate for massless scalar-tensor theory, with massive scalar fields suppressing the effect.","pith_inferences":["Editorial extension: the paper assumes isotropic bolometric surface emission, but real neutron-star atmospheres tend to suppress tangentially emitted photons, so the 80% contrast is likely an upper bound and the practical sensitivity depends on how much of the secondary-image flux survives.","Editorial extension: applying the same ray-tracing pipeline to two-spot or multipolar spot geometries, as many millisecond pulsars require, could either amplify or wash out the secondary-image signature; this is a concrete modelling step the paper leaves undone.","Editorial extension: the paper works at fixed $Q=0.5$; a natural follow-up is a Bayesian fit of GR and STT templates to one high-mass pulsar with an independent radius prior, which would convert the claimed sensitivity into posterior bounds on the scalar charge."],"forward_implications":["For a neutron star near $M/R = 0.284$ with an edge-on spot crossing the far side, the STT light curve can show either a full eclipse or a lensed brightening where GR shows a smooth non-zero flux, so a single pulse shape plus an independent mass can separate the two theories.","Because the critical compactness shifts by up to about 10% between GR and STT for the models considered, fitting STT light curves changes the inferred radius by a similar amount, so ignoring scalarization is a systematic bias in radius measurements.","Even in geometries where the spot never crosses the special region behind the star, including rotation at 700 Hz still produces STT-GR differences from time delay and redshift integrated over a finite spot, so the test is not limited to the edge-on configuration.","Below the critical-compactness regime, STT-GR flux differences stay under about 5%, concentrating the probe in the high-mass, high-compactness pulsars.","Smaller hot spots give more pronounced secondary-image brightening peaks, so future searches should target sources with small, near-equatorial spots."],"supporting_citations":[{"why":"introduced the exponential coupling model that produces spontaneous scalarization, the physical effect the paper aims to make observable.","marker":"[6]"},{"why":"supplied the scalar-tensor differential flux formula and deflection/time-delay integrals that the paper integrates over a finite spot.","marker":"[20]"},{"why":"identified the GR compactness 0.284 at which the entire neutron-star surface becomes visible, setting the multiple-image threshold used throughout.","marker":"[17]"},{"why":"showed how spot size shapes the brightening peak for highly compact stars, which the paper uses to interpret its secondary-image flux.","marker":"[23]"},{"why":"computed finite-spot scalar-tensor light curves with a massive scalar field, providing the comparison regime where STT-GR differences are small.","marker":"[22]"},{"why":"provided the exact Just exterior metric in the Einstein frame that the ray tracing is built on.","marker":"[34]"},{"why":"gave a model-independent bound on the scalar charge for canonical-mass stars, supporting the paper's choice of $Q=0.5$ as still viable at high masses.","marker":"[41]"},{"why":"showed that a massive scalar field suppresses scalarization, justifying the paper's treatment of its massless results as an upper limit.","marker":"[42]"}],"fun_headline_variants":["X-ray pulses from ultradense neutron stars test scalar-tensor gravity","Scalar gravity can twist neutron-star pulse shapes by up to 80%","Ultracompact pulsars: a probe for deviations from general relativity","Pulsar X-ray light curves might expose scalar-tensor gravity","Second image in neutron-star pulse reveals scalar charge"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation rests on treating the spinning, mass-loaded star as a spherical, non-rotating Just spacetime with a fixed scalar charge and isotropic surface emission; if a real atmosphere suppresses the tangentially emitted rays that form the second image, or if rotation deforms the star's exterior, the predicted brightening could shrink or disappear.","fun_headline_variants_meta":{"raw":{"variants":["X-ray pulses from ultradense neutron stars test scalar-tensor gravity","Scalar gravity can twist neutron-star pulse shapes by up to 80%","Ultracompact pulsars: a probe for deviations from general relativity","Pulsar X-ray light curves might expose scalar-tensor gravity","Second image in neutron-star pulse reveals scalar charge"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000778,"raw_usage":{"total_tokens":3531,"prompt_tokens":1130,"completion_tokens":2401,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":746,"completion_tokens_details":{"reasoning_tokens":2311}},"tokens_in":746,"tokens_out":2401,"duration_ms":16641,"temperature":1.0,"reasoning_tokens":2311,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:42:52.071199+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"One concrete check is to recompute the STT pulse profile with a realistic atmospheric radiative-transfer model and an oblate rotating star: if the secondary-image contribution at half rotation phase drops below about 10% of the direct flux, the claimed 80% contrast cannot appear in real data. Observationally, a high-precision phase-resolved observation of a high-mass pulsar with an independently measured mass and radius that is fit well by a GR model at or above $M/R = 0.284$ without any extra brightening would rule out the paper's fiducial scalar charge.","supporting_citations":[{"cited_title":"Silva, N","cited_arxiv_id":null,"evidence_quote":"supplied the scalar-tensor differential flux formula and deflection/time-delay integrals that the paper integrates over a finite spot."},{"cited_title":"Sotani, U","cited_arxiv_id":null,"evidence_quote":"identified the GR compactness 0.284 at which the entire neutron-star surface becomes visible, setting the multiple-image threshold used throughout."},{"cited_title":"Sotani, Light curves from highly com- pact neutron stars with spot size effect","cited_arxiv_id":null,"evidence_quote":"showed how spot size shapes the brightening peak for highly compact stars, which the paper uses to interpret its secondary-image flux."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"computed finite-spot scalar-tensor light curves with a massive scalar field, providing the comparison regime where STT-GR differences are small."},{"cited_title":"Just, Notizen: The motion of mercury according to the theory of thiry and lich- nerowicz","cited_arxiv_id":null,"evidence_quote":"provided the exact Just exterior metric in the Einstein frame that the ray tracing is built on."},{"cited_title":"Model-Independent Comparisons of Pulsar Timings to Scalar-Tensor Gravity","cited_arxiv_id":"1107.3585","evidence_quote":"gave a model-independent bound on the scalar charge for canonical-mass stars, supporting the paper's choice of $Q=0.5$ as still viable at high masses."},{"cited_title":"Ramazano˘ glu, F","cited_arxiv_id":null,"evidence_quote":"showed that a massive scalar field suppresses scalarization, justifying the paper's treatment of its massless results as an upper limit."}],"review_version":1}