{"id":"41028238-4de4-476f-a1ae-c15b6434483d","arxiv_id":"2504.21313","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"During the 2023 outburst of GK Persei, the spin-pulse amplitude grew and the accretion column faded as optical brightness rose, consistent with curtain absorption rather than a change in intrinsic X-ray emission.","lead":"X-ray and optical observations of the 2023 outburst of GK Persei show that the 351-second X-ray pulse of the white dwarf is best explained by rotating accreting gas, the curtain, blocking the hot accretion column. The same data track how the inner disk and the accretion column evolve as the outburst brightens and fades.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Phase-resolved fits in §3.5 do not statistically exclude a pure normalization decrease; without a Δχ² comparison against alternative models, the curtain-absorption origin of the 351.3-s pulse is not established.","rationale":"The reader's weakest assumption is precisely the load-bearing point: the off-pulse dimming is attributed to a higher absorbing column while the intrinsic accretion-column spectrum is assumed unchanged. My review agrees and sharpens it: the phase-resolved fits in §3.5 do not report a statistical comparison against a model in which only the power-law normalization changes, which is the natural alternative (occultation or intrinsic emissivity change). With a fixed photon index and a limited NICER bandpass, NH,max and normalization are degenerate, so the qualitative statement that the column was 'much higher' is not sufficient. The high-energy NuSTAR pulse with almost no energy dependence further illustrates that photoelectric absorption alone cannot explain all pulse components, and the paper itself concedes auxiliary mechanisms. The independent support includes public NICER/NuSTAR data and a standard spectral model, so the concern is testable rather than fatal. The CONDITIONAL verdict from the reader remains appropriate; the condition should be sharpened to require a quantitative model-selection test before the curtain-absorption origin is treated as settled.","tokens_in":20319,"tokens_out":6044,"duration_ms":69336,"concrete_test":"Re-run the on/off-pulse spectral fits of §3.5 for T1–T4 (public NICER data) with two nested models: (A) tie the PL normalization and photon index between on- and off-pulse phases, free only pwab; (B) tie pwab to the time-averaged values, free only the PL normalization. Report Δχ²/ΔAIC and the 90% confidence range of NH,max for each time zone. If model B is not significantly worse, or if the NH,max difference is consistent with zero, the curtain-absorption origin of the pulse is not established. Optionally add a third fit letting both vary and test whether NH,max still changes significantly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim in §4.3 is that the 351.3-s X-ray pulse is produced by spin-phase-dependent photoelectric absorption by the curtain, not by a change in the intrinsic accretion-column emission. The evidentiary basis is the on/off-pulse spectral modeling in §3.5 (Figure 10), where the power-law component drops and the partial-covering column is reported to be 'much higher' off-pulse. The paper does not provide the phase-resolved best-fit parameters, confidence intervals, or a model-selection statistic for these fits. This matters because with the photon index fixed at Γ=1.48 and only the 0.3–8 keV NICER bandpass, the partial-covering column density NH,max and the power-law normalization are strongly degenerate: a lower normalization at constant absorption produces nearly the same spectrum as a higher column at constant normalization, especially when the covering fraction is also free. The prose conflates the observed decrease in absorbed PL flux with an inferred increase in absorber column; that inference is exactly the assumption needing testing. An occultation of the lower part of the column by the WD/disk, or a genuine decline in the column emissivity, would also reduce the PL flux while leaving the BB and line components nearly unchanged, matching Figure 10. The high-energy NuSTAR pulse (Figure 9) with almost no energy dependence is an especially clear region where photoelectric absorption cannot produce the pulse, and the paper itself invokes electron scattering or geometry; this shows the absorption interpretation already needs auxiliary mechanisms. The load-bearing condition for the headline claim is therefore that an absorption-only model is statistically preferred over an intrinsic-flux-change model in the phase-resolved fits. That condition is not demonstrated by the reported analysis.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"Using simultaneous NICER and NuSTAR X-ray observations, optical photometry from Tomo-e Gozen, and AAVSO archive data, the authors study the 2023 outburst of the intermediate polar GK Per. They decompose the X-ray spectrum into a soft blackbody, a set of emission lines, and a multi-temperature bremsstrahlung (accretion-column) component, and track the evolution of these components over the outburst. They detect the 351.3-s white-dwarf spin pulse in X-rays, find that its amplitude and energy dependence evolve with optical brightness, and interpret the pulse as primarily due to spin-phase-dependent photoelectric absorption of the accretion column by the accretion curtain. They also detect a ~5700-s optical QPO and propose an origin linked to an irradiated vertically thick structure near the disk edge, and interpret the outburst light curve as an inside-out disk-instability outburst.","tokens_in":20503,"tokens_out":7691,"duration_ms":73245,"significance":"The observational dataset is rich and the paper covers a broad range of phenomena: spectral evolution, spin-pulse evolution, QPOs, and the outburst mechanism. If the central interpretation holds, it would strengthen the accretion-curtain model for GK Per and demonstrate how pulse properties respond to changes in mass accretion rate. The paper explicitly reports the spectral model, shows the data and fits, and connects the pulse evolution to the spectral evolution. However, the central claim concerning the curtain-absorption origin of the spin pulse is not yet supported by a fully quantitative statistical analysis: the phase-resolved spectral fit parameters and model-comparison statistics are not reported. The paper would be strengthened by providing those details and by clarifying the role of alternative mechanisms at high energies.","major_comments":[{"comment":"The paper states that the column density of the absorber was much higher in the off-pulse phase than in the on-pulse phase, but it does not quote the phase-resolved best-fit values of NH,max, the covering fraction, the power-law normalization, or their confidence intervals, nor any Δχ² between the on- and off-pulse fits. Because the 0.3–8 keV band with Γ fixed at 1.48 leaves NH,max and the power-law normalization partially degenerate, the conclusion in Section 4.3 requires a statistical comparison against an alternative model in which only the power-law normalization changes or the column is held fixed. Please provide these numbers, or explicitly restrict the claim to a qualitative statement.","section":"Section 3.5 / Figure 10"},{"comment":"The pulse amplitudes are quoted (for example, 'consistently lower than 20%' at <1 keV and 'highest around the outburst maximum') without formal uncertainties. The claimed energy dependence and its evolution (weaker energy dependence when the system is brighter) are central to the curtain-density interpretation; without amplitude errors, the significance of this evolution cannot be evaluated. Add error bars or confidence intervals for the pulse amplitudes in each energy band and time zone.","section":"Section 3.5 / Figure 8"},{"comment":"Fixing the photon index at Γ=1.48 for all time zones and for both pulse phases may bias the inferred pwab column density if the intrinsic spectrum varies. The authors should either free Γ in the phase-resolved fits and quote its best-fit value and error, or explicitly test and report that the difference in NH,max between on-pulse and off-pulse fits survives when Γ is allowed to vary.","section":"Sections 3.3 and 3.5"},{"comment":"The NuSTAR pulse in the 3–50 keV band shows almost no energy dependence, which the paper attributes to electron scattering or a geometrical change in the visible area of the accretion column. This means photoelectric absorption cannot be the sole, or even the dominant, cause of the pulse at high energies. To make the headline claim in the abstract and Section 4.3 more precise, the paper should quantify the relative contributions of absorption, scattering, and geometry, or explicitly limit the curtain-absorption interpretation to below ~8 keV.","section":"Section 4.3 / Figure 9"}],"minor_comments":[{"comment":"The phrase 'while the the BB and line fluxes' contains a duplicated article.","section":"Section 3.5"},{"comment":"The definition of time zone T4 as 'BJD 2459900–2460000' appears to be a typo; from the context it should be BJD 2459990–2460000.","section":"Section 3.5"},{"comment":"The significance of the PDM detections (351.32 s and 5699 s) should be quantified with false-alarm probabilities (for example, via bootstrap or Monte Carlo), not only the θ statistic.","section":"Section 3.4"},{"comment":"The phrase '/greaterorsimilar10 keV' is a LaTeX error; it should read '≳10 keV'.","section":"Section 4.2"},{"comment":"The cross-calibration factor of 0.98 between the NuSTAR and NICER normalizations is mentioned in the Figure 3 caption but should also be stated in the text of Section 3.2.","section":"Section 3.2 / Figure 3"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a valuable dataset and a mostly coherent interpretation. The main reservation is the lack of quantitative detail for the phase-resolved spectral fits that underpin the central claim; this should be addressed before publication. The authors should also be careful not to overstate the 'major cause' claim given the high-energy pulse behavior."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a straightforward observational study of the 2023 GK Per outburst: NICER monitoring, one simultaneous NuSTAR spectrum, Tomo-e Gozen high-speed optical photometry, and AAVSO coverage. The new data are real and the paper does a good job of describing the spectral evolution: the broadband fit with blackbody, line emitter, and cooling-flow components is internally consistent, and the on-pulse/off-pulse contrast—where the power-law flux drops while the blackbody and lines stay roughly constant—is a nice, direct handle on the pulse mechanism. The light-curve and spin-pulse evolution over the outburst is cleanly presented. The paper is an honest extension of earlier work (Zemko et al. 2017; Pei et al. 2024), not an overclaiming novelty.\n\nThe soft spot is exactly where the stress-test note points. Section 4.3 says the column density of the absorber was 'much higher' off-pulse, and uses that to conclude the spin pulse is caused by curtain absorption. But the phase-resolved fits in §3.5 are not shown with errors or a model-selection statistic. With the photon index fixed and only the NICER bandpass, a decrease in the power-law normalization at constant absorption can mimic an increase in partial-covering column. The paper does not report Δχ² against an intrinsic-flux-change model, so the central inference—that the absorber column, not the column emissivity, changes—is not actually established by statistical evidence. That does not make the interpretation wrong; it makes it under-supported. The paper itself admits auxiliary mechanisms (electron scattering, geometry) are needed for the energy-independent hard pulse, which is worth reading as a signal of how far the absorption model goes.\n\nOther issues are minor: pulse amplitudes are quoted without formal uncertainties; the T4 BJD interval in §3.5 is a typo (2459900–2460000 should be 2459990–2460000); no code or parameter files are shipped. The paper's own discussion of unresolved line-emitter and QPO origins is honest.\n\nBottom line: this is a paper for people who work on intermediate polars and dwarf-nova outbursts. It deserves a serious referee—the data are useful, the analysis is mostly standard, and the central claim is plausible but needs a quantitative model comparison before it can be taken as settled. I would send it to review with the expectation that the authors test the absorption-only interpretation against an intrinsic-flux-change model.","headline":"A solid, useful observational paper on the 2023 GK Per outburst, but the curtain-absorption interpretation of the pulse rests on a model comparison the paper does not actually report.","tokens_in":21285,"tokens_out":3064,"would_cite":true,"duration_ms":28383,"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":"During GK Persei's 2023 outburst, the white dwarf's 351.3-second X-ray spin pulse was produced by spin-phase-dependent photoelectric absorption of the accretion column by the curtain, not by intrinsic variability of the column.","keywords":["cataclysmic variable stars","dwarf novae","DQ Herculis stars","intermediate polars","X-ray sources","stellar accretion disks","accretion curtain","white-dwarf spin pulse"],"falsifier":"Compare the fully deabsorbed accretion-column flux between on-pulse and off-pulse phases by removing the best-fit partial-covering column: if the model is right, the intrinsic fluxes should match. If the intrinsic flux still differs significantly after deabsorption, or if the 6.4 keV iron line and the continuum drop together without an absorption-edge increase at off-pulse phase, the pulse is not caused by curtain absorption.","tokens_in":19989,"feed_emoji":"🛰️","tokens_out":7223,"duration_ms":70372,"temperature":0.7,"pith_summary":"This paper argues that the 351.3-second X-ray spin pulse of the white dwarf in GK Persei, seen throughout the 2023 outburst, is produced by photoelectric absorption of the accretion column by the surrounding accreting gas, the accretion curtain, rather than by the column intrinsically brightening and fading. The evidence comes from phase-resolved X-ray spectra: at off-pulse phase the continuum from the accretion column drops sharply while the fitted absorbing column rises, and the other spectral components—the soft blackbody from the irradiated white-dwarf surface and the keV emission lines—hardly change. This absorption picture also explains why the pulse deepens and its energy dependence weakens as the optical brightness and accretion rate rise. If the claim is correct, the pulse depth becomes a direct, quantitative probe of the geometry and density of the inner accretion flow during outbursts, with implications for how matter is channeled onto magnetic white dwarfs.","feed_headline":"Spin pulse in GK Persei traced to absorbing accretion curtain","feed_subtitle":"NICER and NuSTAR phase-resolved spectra tie the 351-second pulse to absorption by the accreting curtain.","key_machinery":"The central object is the accretion curtain: the magnetically channeled gas flow from the truncated disk onto the white dwarf, modeled as a partial photoelectric absorber (the pwab component in XSPEC) whose column density along the line of sight changes with the white dwarf's 351.3-second spin phase. Phase-resolved fits of on-pulse and off-pulse spectra, using a fixed intrinsic power-law continuum, measure this phase-dependent column; the constancy of the blackbody and line components across phase isolates the column as the only varying absorber, carrying the argument that the pulse is an absorption effect.","core_discovery":"The paper claims that during the 2023 outburst of GK Persei the 351.3-second X-ray spin pulse of the white dwarf is produced by spin-phase-dependent photoelectric absorption: at the off-pulse phase the partial-covering absorbing column rises to about $10.7\\times10^{22}$ cm$^{-2}$ while the intrinsic multi-temperature bremsstrahlung spectrum of the accretion column remains essentially unchanged. The blackbody component from the irradiated white-dwarf surface and the keV-scale emission lines vary little between on- and off-pulse phases, so the pulse is not intrinsic to the column's emission. The same absorption picture accounts for the observed correlations: as the optical brightness and inferred mass accretion rate rise, the curtain column grows, the observed column flux drops, the pulse amplitude increases, and the energy dependence of the pulse below 8 keV weakens. At the outburst end, the emergence of a double-peaked high-energy pulse and the near energy-independence above 10 keV are attributed to a receding inner disk edge that exposes the second pole, and to electron scattering or geometric visibility, respectively.","pith_inferences":["If the curtain-absorption interpretation is right, the spin pulse amplitude can serve as a close-to-real-time gauge of the mass accretion rate onto the white dwarf, since the curtain column density tracks the accretion rate; this would make intermediate polars in outburst useful clocks for disk-instability dynamics.","A direct way to separate absorption from occultation observationally is to track the iron K$\\alpha$ fluorescence line and the photoelectric absorption edge across spin phase: an absorption origin predicts the edge depth and the line-to-continuum ratio change together, whereas occultation would dim the continuum and line together without an edge change.","The same phase-resolved spectral decomposition could be applied to the 2018 outburst data of GK Persei and to other intermediate polars with dwarf-nova outbursts to test whether the pulse-depth versus accretion-rate scaling found here is a general property of the accretion-curtain mechanism rather than specific to this outburst.","The claim that the ~5700 s optical quasi-periodicity comes from an irradiated vertically thick structure rotating at its Keplerian radius predicts that the period should stay constant across the outburst; a future long-cadence optical campaign during an outburst could check that directly."],"forward_implications":["In brighter states, the pulse amplitude at >2 keV grows while the observed accretion-column flux drops, so the two can be used together to track the curtain column density and hence the mass accretion rate.","Pulse profile shape and the low-energy pulse peak lag are tied to the twisted geometry of the accretion column and curtain, meaning pulse timing and phase lags carry structural information about the inner disk edge.","The pulse is weak in the 1-2 keV band until mid-outburst because the line-emitting region outside the curtain dilutes the modulation; the pulse appears in that band once the line flux fades.","At the end of the outburst, the double-peaked high-energy pulse indicates that the inner disk edge expanded and the second magnetic pole became visible, providing a geometric probe of the inner radius.","The outburst itself is an inside-out disk-instability event with a stagnation phase, consistent with the disk-instability model applied to magnetically truncated disks."],"supporting_citations":[{"why":"Introduced the accretion-curtain model that connects spin-phase-dependent absorption to the X-ray pulse in intermediate polars.","marker":"Rosen et al. 1988"},{"why":"Developed the curtain interpretation and its observational expectations, which the paper's spectral fits directly test.","marker":"Hellier et al. 1991"},{"why":"Discovered the 351-s white-dwarf spin pulse in GK Persei, the signal whose origin this paper establishes.","marker":"Watson et al. 1985"},{"why":"Provided the abundance and column-density framework and the quiescent spectral parameters that the 2023 fits extend.","marker":"Zemko et al. 2017"},{"why":"Reported the pulse-amplitude growth with accretion rate in the 2018 outburst, the trend continued here with phase-resolved spectra.","marker":"Pei et al. 2024"},{"why":"Supplied the mass-attenuation coefficients that show electron scattering dominates above ~10 keV, invoked for the weak high-energy energy dependence.","marker":"Hayashi et al. 2018"}],"fun_headline_variants":["Accretion curtain, not column, shapes GK Persei's X-ray pulse","GK Persei's spin pulse is a curtain effect, not an emission change","Curtain absorption explains white dwarf's 351-s X-ray pulse","X-ray pulse in GK Persei: absorption, not emission, drives it","Spin pulse in GK Persei: absorbing curtain is the cause"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The off-pulse dimming is interpreted as extra photoelectric absorption of an otherwise unchanged accretion-column spectrum; if the dimming were instead caused by the white dwarf or disk hiding the column, or by the column itself emitting less at that phase, the curtain-absorption conclusion fails.","fun_headline_variants_meta":{"raw":{"variants":["Accretion curtain, not column, shapes GK Persei's X-ray pulse","GK Persei's spin pulse is a curtain effect, not an emission change","Curtain absorption explains white dwarf's 351-s X-ray pulse","X-ray pulse in GK Persei: absorption, not emission, drives it","Spin pulse in GK Persei: absorbing curtain is the cause"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00099,"raw_usage":{"total_tokens":4243,"prompt_tokens":1038,"completion_tokens":3205,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":654,"completion_tokens_details":{"reasoning_tokens":3114}},"tokens_in":654,"tokens_out":3205,"duration_ms":21638,"temperature":1.0,"reasoning_tokens":3114,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:06:42.384205+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare the fully deabsorbed accretion-column flux between on-pulse and off-pulse phases by removing the best-fit partial-covering column: if the model is right, the intrinsic fluxes should match. If the intrinsic flux still differs significantly after deabsorption, or if the 6.4 keV iron line and the continuum drop together without an absorption-edge increase at off-pulse phase, the pulse is not caused by curtain absorption.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Developed the curtain interpretation and its observational expectations, which the paper's spectral fits directly test."},{"cited_title":"G., King, A","cited_arxiv_id":null,"evidence_quote":"Discovered the 351-s white-dwarf spin pulse in GK Persei, the signal whose origin this paper establishes."},{"cited_title":"2024, MNRAS, 529, 1463, doi: 10.1093/mnras/stae650","cited_arxiv_id":null,"evidence_quote":"Reported the pulse-amplitude growth with accretion rate in the 2018 outburst, the trend continued here with phase-resolved spectra."}],"review_version":1}