{"id":"9317cdff-1f8e-4966-9e7e-709dd5385a0c","arxiv_id":"2411.16528","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"GRMHD simulations with relativistic ray tracing show that realistic accretion hotspots are crescent- or band-shaped, not circular, and produce X-ray pulse amplitudes consistent with accreting millisecond pulsars.","lead":"This paper uses 3D general-relativistic magnetohydrodynamic simulations of accreting neutron stars to map the shape of the X-ray hotspots that produce millisecond pulsations. The simulations yield crescent and band-like hotspot shapes and pulse amplitudes in the 1-12% range seen in accreting millisecond pulsars, providing physically motivated inputs for neutron star mass-radius inference.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quantitative pulse-amplitude and harmonic results rest on the isotropic blackbody/no-shock emission assumption that the authors themselves flag; a realistic Comptonized beaming model could shift amplitudes out of the observed range.","rationale":"The reader's weakest-assumption analysis identifies the isotropic blackbody/no-shock surface emission model as the key vulnerability, and the paper's own text supports this: Section 2.2.1 states that the shock and Comptonization component are not modeled, and Section 4 cautions that amplitudes may be overestimated. I agree that this is the single most load-bearing concern because the paper's quantitative bridge to observations—the 1–12% rms range and its consistency with AMXP amplitudes—passes through this emission prescription. The qualitative claim that GRMHD accretion columns produce crescent- and band-shaped footprints is independently supported by the GRMHD simulations and is consistent with earlier non-relativistic MHD results, so I do not see a reason to reject the paper. The cross-code agreement between RAPTOR and X-PSI within roughly 1% is genuine independent support for the ray-tracing pipeline, and the reported convergence check on camera resolution is a positive sign. However, no uncertainty or sensitivity analysis is provided for the reported rms amplitudes, harmonic ratios, or variability coefficients, and no simulation snapshots or analysis scripts are released, which limits independent reproduction. These secondary issues reinforce the conditional verdict but do not by themselves invalidate the central qualitative conclusion. A beamed-emission test would settle whether the quantitative consistency claim survives a more realistic surface model; until then, the conditional verdict is appropriate.","tokens_in":20624,"tokens_out":6020,"duration_ms":72010,"concrete_test":"Recompute the Section 3.2.1 pulse profiles and Figure 5 using RAPTOR with the same time-averaged GRMHD temperature maps, but replace the isotropic Planck intensity with an anisotropic Comptonized beaming model calibrated to AMXPs, for example the Poutanen & Gierliński (2003) intensity pattern or a numerical neutron-star atmosphere/Comptonization model. Compare the fractional rms amplitudes of the fundamental and first harmonic across observer inclinations and magnetic inclinations.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—that simulated fractional rms amplitudes of the fundamental and first harmonic usually range between 1 and 12% and are consistent with observed AMXPs—is computed from surface emission modeled as isotropic blackbody, with the effective temperature set by converting the matter-energy flux via Eq. (12). The radiation-dominated shock and Comptonization component that produce the actual AMXP spectra are explicitly not modeled (Section 2.2.1), and the authors state that pulse amplitudes may be slightly overestimated because of the isotropic blackbody treatment (Section 4). This is load-bearing because pulse fractional amplitudes and harmonic ratios are known to depend sensitively on the angular emission pattern and on the temperature structure of the emitting region. The same non-circular hotspot footprint, when viewed through a Comptonizing shock with anisotropic beaming, can yield significantly different rms amplitudes and a different fundamental-to-harmonic hierarchy. Since the comparison to observed AMXP amplitudes is the main quantitative support for the claim that these realistic hotspot geometries matter, the assumption is not a minor calibration detail. The qualitative conclusion that GRMHD hotspots are non-circular may survive, but the quantitative consistency with observations, and hence the strength of the argument for moving beyond circular spots, is conditional on this emission model.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses global 3D GRMHD simulations (BHAC) of an accreting neutron star with a dipolar magnetic field for three magnetic inclinations (30, 60, and 90 degrees) to extract surface hotspot morphologies, convert the matter-energy flux to an effective temperature via Eq. (12), and compute X-ray pulse profiles with two independent general-relativistic ray-tracing codes (RAPTOR and X-PSI). The main claims are that the hotspots are crescent-shaped at low inclination and elongated bands/bars at higher inclination; that the simulated fractional rms amplitudes of the fundamental and first harmonic usually lie in the 1–12% range, consistent with observed accreting millisecond pulsars; that turbulent accretion flow adds broadband variability on timescales comparable to the spin period; and that electron scattering absorption in the accretion column and disk modifies pulse shapes and introduces additional variability and higher harmonics.","tokens_in":20822,"tokens_out":4830,"duration_ms":50995,"significance":"If the results hold, the paper provides a strong, physically motivated argument that pulse-profile modeling of accretion-powered pulsars should move beyond circular hotspot geometries, a step that is currently missing from standard AMXP pulse-profile analyses. The study is technically careful in several respects: the two ray-tracing codes agree to within about 1%, the camera convergence is checked at the 0.25% level, the harmonic-fit residuals are within 1.2%, and time-dependent ray tracing of individual GRMHD snapshots is used to quantify variability. The authors are also transparent about the main simplifications. However, the quantitative pulse-amplitude predictions, which are the basis for the claimed consistency with observations, rest on an isotropic blackbody emission model with no shock or Comptonization, and the electron-scattering treatment is absorption-only. These are load-bearing assumptions for the numerical values, so the strength of the quantitative conclusions is conditional on them.","major_comments":[{"comment":"The quantitative pulse-amplitude range (1–12% rms) is computed by assuming that the accretion energy is radiated locally as an isotropic blackbody with Teff = (Fm/σSB)^(1/4), with no radiation-dominated shock and no Comptonization. The authors state in Section 4 that the amplitudes may be slightly overestimated because of this assumption. This is load-bearing: pulse rms amplitudes and the fundamental-to-first-harmonic ratio depend sensitively on the angular pattern of surface emission and on the temperature structure of the emitting region. A Comptonized beaming pattern applied to the same non-circular hotspot footprints could shift the predicted amplitudes and harmonic hierarchy relative to the observed AMXP values. I request either a robustness test using a simple anisotropic beaming prescription (for example, a limb-darkening or Comptonized-beaming factor) or a visible softening of the abstract/conclusion statements so that the conditional nature of the amplitude comparison is explicit.","section":"Section 2.2.1 and Eq. (12); discussed in Section 4"},{"comment":"The electron-scattering calculation includes only Thomson absorption opacity (alpha_nu = 0.4 cm^2/g) and ignores line-of-sight scattering and energy exchange. The conclusions that the accretion column and disk introduce significant additional variability and higher harmonics are based on this absorption-only model. Pure absorption removes photons and can overproduce harmonic structure and variability compared with a scattering treatment that redistributes photons in angle and energy. Although the limitation is acknowledged in Section 4, the quantitative values quoted in Section 3.3 (for example, variability increasing from 0.18 to 0.213 at i_obs = 50 degrees and from 0.182 to 0.391 at i_obs = 70 degrees) should be presented more cautiously, for instance as upper limits or as estimates pending a scattering calculation, because they are used to support the harmonic and variability conclusions.","section":"Section 3.3 and Eq. (8)"}],"minor_comments":[{"comment":"Equation (14) fits only the fundamental and first harmonic, but Section 3.3 quotes amplitudes for the second, third, and fourth harmonics. Please specify the extended fit function (and the number of harmonics) used for those cases, since the reported fundamental and first-harmonic amplitudes in the absorption runs could otherwise be biased by unmodeled higher harmonics.","section":"Section 3.3 and Eq. (14)"},{"comment":"The rms values quoted for the fundamental and first harmonic (rms_nu0 and rms_nu1) are obtained from Lorentzian fits to the PSDs, but the fitting procedure is not described in enough detail. Please state the Lorentzian width, the integration range, and how the broadband noise floor is estimated, because the quoted numbers depend on these choices.","section":"Section 3.2.2 and Figure 7"},{"comment":"The coefficients of variation for phi2 are reported as negative values (c_phi2 = -0.279 and -0.344). Since c_x is defined as sigma_x/<x>, a negative value is possible when the mean phase is negative, but the sign is not physically informative. Consider reporting the absolute value or a circular dispersion measure instead.","section":"Section 3.2.2"},{"comment":"The labels 'Middle' and 'Top' in the caption are somewhat unclear; please clarify explicitly that the top and middle rows correspond to the upper and lower magnetic hemispheres, respectively.","section":"Figure 2 caption"},{"comment":"The plasma beta is defined as beta_t = 2p/b^2, but the text later refers simply to 'beta' transitioning smoothly outside the light cylinder. Please use a consistent definition of beta throughout, or clarify that the same definition is implied.","section":"Section 2.1"}],"recommendation":"major_revision","confidential_remarks":"This is a well-executed and transparent simulation study, and the qualitative claim about non-circular hotspots is timely and important for the AMXP pulse-profile modeling community. The main risk is the dependence of the quantitative amplitude and harmonic predictions on the simplified emission model. The requested revisions—adding a beaming-sensitivity test or tempering the quantitative claims, and clarifying the absorption-only scattering treatment—are feasible within the scope of the manuscript. I do not see grounds for rejection, but the central quantitative consistency claim needs to be made conditional or robustified before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: this is a capable computational paper that delivers a genuinely new result—pulse profiles from GRMHD-simulated accretion hotspots for AMXPs—and its qualitative message (hotspots are crescents and bands, not circles) will survive. The quantitative amplitude comparison is shakier, and the authors tell you why.\n\nWhat's actually new: first use of global 3D GRMHD with MRI-driven turbulence to compute hotspot shapes and ray-traced pulses, for inclinations up to 90 degrees. The crescent-to-band transition, the broadband variability from the turbulent flow, and the absorption-induced harmonics from electron scattering are real additions. The code tests are properly done: RAPTOR vs X-PSI agree within ~1%, camera convergence at 0.25%, fit residuals under 1.2%. That is evidence the numerics are under control.\n\nSoft spots, in order: (1) The quantitative pulse amplitudes (1-12% rms) and harmonic ratios come from an emission model that sets the surface to isotropic blackbody with Teff from the matter flux and explicitly drops the shock/Comptonization component. The authors flag this in Sec 2.2.1 and Sec 4, and admit amplitudes may be overestimated. Pulsed fractions and harmonic hierarchy are known to be sensitive to angular beaming, so this assumption is load-bearing for the 'consistent with observed AMXPs' claim. It does not kill the qualitative geometry result, but it means the quantitative consistency deserves an 'in this model' qualifier. (2) The absorption treatment is Thomson opacity only, no scattering into the line of sight; the authors acknowledge this but it affects the higher-harmonic predictions. (3) No uncertainties or sensitivity analyses on the variability coefficients and PSD integrals, and no public data beyond what's in the paper. These are minor-to-moderate.\n\nThe circularity point in the reader's report is fine—they are not fitting, just comparing post hoc, so I don't hold that against them. Self-citation to Das & Porth 2024 is appropriate since the simulations come from there.\n\nAudience: people building pulse profile models for AMXPs, and those planning X-ray timing missions. It deserves a serious referee; the referee should push for a sensitivity study on the emission model or at least a statement of how much the amplitudes would change under Comptonized beaming. I'd bring it to a reading group, mainly to discuss the emission-model caveat.","headline":"Solid GRMHD-era pulse profile paper: the non-circular hotspot results will hold, but the quantitative amplitude claims rest on an emission model the authors themselves flag as approximate.","tokens_in":21421,"tokens_out":2447,"would_cite":true,"duration_ms":23376,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Accreting millisecond pulsar hotspots are crescent- and band-shaped, not circular, and their simulated pulse profiles match observed amplitudes of 1-12% rms.","keywords":["accreting millisecond pulsars","neutron star hotspots","GRMHD simulations","pulse profiles","ray tracing","electron scattering","magnetic inclination","X-ray variability"],"falsifier":"Compare a sample of accreting millisecond pulsars with independently known magnetic and observer inclinations against the simulated crescent and band hotspot family: if observed fractional rms amplitudes consistently fall outside the 1-12% range, or the observed harmonic ratios contradict the simulated trend with magnetic inclination, the central claim would be falsified.","tokens_in":20317,"feed_emoji":"🌟","tokens_out":6854,"duration_ms":62245,"temperature":0.7,"pith_summary":"This paper uses global three-dimensional general-relativistic magnetohydrodynamic simulations of a magnetized neutron star accreting from a disk to ask what shape the X-ray-emitting hotspots on the stellar surface really take, and what pulses those hotspots produce. It finds that the hotspots are not circular: at low magnetic inclination they are crescents around the magnetic axis, and at higher inclination they become elongated bands, with the morphology set entirely by the accretion columns. Ray-traced pulse profiles from these spots have fractional rms amplitudes of the fundamental and first harmonic mostly between 1 and 12 percent, matching the range seen in accreting millisecond pulsars. The paper argues this is strong evidence that pulse-profile modeling of accretion-powered pulsars must abandon the circular-spot assumption. It also shows that the turbulent accretion flow adds broadband variability and that electron-scattering absorption in the accretion column adds higher harmonics and additional variability.","feed_headline":"Accreting pulsar hotspots come as crescents and bands, not circles","feed_subtitle":"Simulated pulse amplitudes of 1-12% rms match real accreting millisecond pulsars—circular-spot models should go.","key_machinery":"The machinery is a chain from simulation to synthetic observation: global three-dimensional general-relativistic magnetohydrodynamic simulations of a dipolar magnetosphere accreting from a turbulent disk produce surface matter-energy flux maps; the flux is converted to an effective blackbody temperature via $T_{\\rm eff} = (F_m/\\sigma_{\\rm SB})^{1/4}$; and two independent general-relativistic ray-tracing codes, one forward in time and one backward, turn the surface temperature maps into bolometric pulse profiles, including stellar oblateness, Doppler boosting, light bending, and, in one code, radiative transfer of Thomson absorption through the accretion column. The spot shapes themselves come from the accretion columns: the star-disk magnetic interaction makes the columns, and thus the spots, non-axisymmetric crescents at low inclination and bands at high inclination. The pulse analysis fits $f(\\phi) = A + B\\sin(\\phi+\\phi_1) + C\\sin(2\\phi+\\phi_2)$ and reports the fractional rms of the fundamental and first harmonic.","core_discovery":"The central claim is that the accretion hotspots of accreting millisecond pulsars, as produced self-consistently by global GRMHD simulations of the star-disk system, form a family of crescent and elongated band shapes that depend on the stellar magnetic inclination, and that these non-circular geometries imprint characteristic pulse profiles. For magnetic inclinations of 30, 60, and 90 degrees, the time-averaged spots range from crescents near the magnetic axis to an equatorial belt, and the fractional rms amplitudes of the fundamental and first harmonic of the resulting pulses usually lie between 1 and 12 percent, consistent with observed accreting millisecond pulsars. The paper further claims that the turbulent accretion flow produces pulse-to-pulse variability on timescales comparable to the rotation period, and that electron scattering in the accretion column and disk obscures the antipodal hotspot, strengthens the fundamental, introduces higher harmonics, and increases overall variability. The conclusion drawn is that pulse-profile modeling for accretion-powered pulsars needs to move beyond circular spot shapes.","pith_inferences":["If the simulated spot family is representative, current circular-spot pulse-profile fits to accreting millisecond pulsars may be underestimating geometric systematics in inferred neutron-star masses and radii; a natural next step is to fit the crescent and band family to archival data.","The predicted trend that the first harmonic strengthens and the fundamental weakens as magnetic inclination increases toward 90 degrees gives a population-level test: sources with independently inferred high magnetic obliquity should show stronger harmonic content.","The absorption-driven growth of the fundamental with rising mass accretion rate implies that an individual pulsar should show systematic pulse-shape evolution across an outburst, with higher harmonics appearing at the bright end; this could be checked with existing long monitoring campaigns."],"forward_implications":["Pulse-profile models for accretion-powered pulsars should use non-circular hotspot shapes, because fitting circular spots may bias inferred mass, radius, and geometry.","Simulated fractional rms amplitudes of 1-12% for the fundamental and first harmonic match observed accreting millisecond pulsars, supporting accretion-column-determined hotspots as the origin of their pulsations.","The turbulent accretion flow produces significant broadband variability that exceeds the pulsed variability by a factor of 1.3-3.3, so observed light-curve variability should be treated as partly accretion noise.","Electron scattering in the accretion column and disk obscures the antipodal hotspot at higher accretion rates, increasing the fundamental amplitude and generating higher harmonics, so pulse shape should evolve with mass accretion rate.","Time-variable spot shapes on rotation-period timescales cause pulse-to-pulse variability with larger scatter in the fundamental amplitude and phase than in the first harmonic."],"supporting_citations":[{"why":"Supplies the global GRMHD simulations of neutron stars with inclined magnetospheres, including the accretion column dynamics and surface flux maps analyzed here.","marker":"Das & Porth 2024"},{"why":"Earlier MHD simulations of MRI-driven accretion onto stars at low inclinations produce near ring-like hotspots, providing the comparison benchmark for the crescent shapes.","marker":"Romanova et al. 2012"},{"why":"Pulse-profile modeling of accreting millisecond pulsars that adopts circular hotspot geometries, the assumption this paper argues against.","marker":"Salmi et al. 2018"},{"why":"Recent AMXP pulse-profile model also using circular spots, relevant to the conclusion that modeling must move beyond circular spots.","marker":"Dorsman et al. 2025"},{"why":"Analytical model showing electron scattering in the accretion column affects pulse evolution, motivating and providing comparison for the absorption treatment.","marker":"Ahlberg et al. 2024"},{"why":"Describes the Comptonized power-law component and anisotropic beaming of AMXP spectra, used to justify ignoring Comptonization and to note possible amplitude overestimation.","marker":"Poutanen & Gierliński 2003"},{"why":"Review defining the observational context of accreting millisecond pulsars and the hotspot origin of their pulsations.","marker":"Patruno & Watts 2021"},{"why":"One of the observational references establishing the typical 3-10% rms fractional amplitudes of AMXP pulsations used for comparison.","marker":"Wijnands & van der Klis 1998"}],"fun_headline_variants":["Accreting pulsars: hotspots are crescents, not circles, simulation shows","GRMHD: Neutron star hotspots shaped like crescents and bands","Pulsar hotspots defy circular models, new simulation finds","Non-circular hotspots explain accreting pulsar pulses","Simulated pulsar pulses match real amplitudes (1-12% rms)"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The simulations assume that the hotspot radiates as an isotropic blackbody with temperature set by the matter-energy flux, ignoring the radiation-dominated shock, Comptonization, and anisotropic beaming that real accreting millisecond pulsars show, so the quantitative pulse amplitudes and shapes could shift if a more complete emission model were used.","fun_headline_variants_meta":{"raw":{"variants":["Accreting pulsars: hotspots are crescents, not circles, simulation shows","GRMHD: Neutron star hotspots shaped like crescents and bands","Pulsar hotspots defy circular models, new simulation finds","Non-circular hotspots explain accreting pulsar pulses","Simulated pulsar pulses match real amplitudes (1-12% rms)"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000806,"raw_usage":{"total_tokens":3592,"prompt_tokens":1048,"completion_tokens":2544,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":664,"completion_tokens_details":{"reasoning_tokens":2451}},"tokens_in":664,"tokens_out":2544,"duration_ms":18214,"temperature":1.0,"reasoning_tokens":2451,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:02:51.830077+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare a sample of accreting millisecond pulsars with independently known magnetic and observer inclinations against the simulated crescent and band hotspot family: if observed fractional rms amplitudes consistently fall outside the 1-12% range, or the observed harmonic ratios contradict the simulated trend with magnetic inclination, the central claim would be falsified.","supporting_citations":[],"review_version":1}