Peter Kretschmar

Context. Galactic high-mass X-ray binaries (HMXBs) are important systems for studying accretion mechanisms onto compact objects and for investigating the complex stellar winds of massive stars. In particular, HMXBs that host a neutron star allow us to reveal the structure of the accreted material in X-ray pulsars and consequently to investigate how matter behaves under extreme conditions of pressure and density. These are major scientific goals for XRISM and NewAthena . Aims. Here we report on the first out-of-eclipse XMM-Newton observation of the HMXB 4U 1538−52, complemented by NuSTAR coverage. Our campaign aimed to investigate stellar-wind variability and continuum changes with high-resolution spectroscopy at a critical orbital phase: when the neutron star is in inferior conjunction. Methods. Thanks to simultaneous observations covering both soft and hard X-rays, we obtain the most detailed X-ray view of the accreted material in 4U 1538−52 to date. In particular, we perform time-resolved spectroscopy down to the pulse period of the neutron star to highlight wind clumping properties and accretion structures. Results. In this dataset, we observe a bright flare reaching ∼10 ³⁷ erg s ⁻¹ probably induced by the accretion of a 10 ²⁰ g clump, followed by a luminosity dip forming a heartbeat-like episode. This event is followed by three local absorption peaks with local variability of the order of the pulse period and a gradual hardening of the underlying spectrum throughout the observation. This could indicate the presence of both small-scale and large-scale overdense structures in the vicinity of the neutron star, which can be attributed to clumps and filamentary structures embedded in the accretion wake. These observational evidences are further supported and reproduced by 3D hydrodynamic simulations.

Context. Pulse profiles probe the emission geometry of accreting X-ray pulsars, but their observed shapes may depend on instrumental response and observational setup. The pulsed fraction spectrum provides a compact spectro-timing observable that can both trace localized spectral features and serve as a quantitative cross-calibration diagnostic. Aims. We assess the consistency of energy-resolved pulse profiles obtained with simultaneous XMM-Newton/EPIC-pn and NuSTAR/FPM observations of Vela X-1 and investigate the broadband pulsed fraction spectrum as a diagnostic of spectral features from 1 to 70 keV. Methods. We constructed energy–phase matrices for both instruments and derived pulsed fraction spectra after carefully accounting for instrumental and observational effects. We quantified the residual systematics in the overlapping 3–10 keV band. We then modeled the broadband pulsed fraction spectra phenomenologically and searched for timing signatures of spectral features. Results. After correcting for instrumental effects, the pulsed fraction spectra derived strictly over the common exposure intervals of the two instruments agree within 5% in their overlapping 3–10 keV range. Remaining discrepancies larger than 5% are confined to the iron-line region and can be attributed to the different energy resolutions of the two instruments. The broadband pulsed fraction spectrum reveals significant localized features corresponding to known emission lines in the soft band and to the fundamental and harmonic cyclotron resonant scattering features at ∼25 and ∼55 keV. An orbital-phase-resolved modeling of the pulsed fraction spectrum of EPIC-pn shows that the soft band features strongly depend on the value of the equivalent absorption column, with emission line signatures becoming progressively suppressed during highly absorbed intervals. Conclusions. The pulsed fraction spectrum serves both as a quantitative cross-calibration diagnostic and as a powerful spectro-timing diagnostic. Its modeling provides independent constraints on spectral features, complementing traditional phase-averaged spectroscopy.

We present a comprehensive analysis of the X-ray spectrum of the highly obscured high-mass X-ray binary IGR J16318−4848 based on ∼250 ks Chandra/HETGS observations covering different orbital phases. We detected a strong narrow Fe K α emission line and interpret our finding as the iron line originating in the stellar wind close to the compact object. We tentatively identify a Compton shoulder associated with the Fe K α line. The ratio of the Compton-shoulder flux to that of the Fe K α core exhibits orbital-phase variability, reaching a maximum of ∼10% around orbital phase ϕ _orb ≃ 0.33. In addition, we do not detect ionized lines, implying that the plasma in the system is predominantly cold. We preliminarily report that the Fe K α line centroid exhibits an orbital modulation with an amplitude of 167 ± 32 km s ⁻¹ , adopting the best-known circular orbit in the source. This study presents the first detailed view of the high-resolution X-ray spectra of IGR J16318−4848 across different orbital phases, providing timely context for the newly taken XRISM observations of the source.

Modeling the observed emission from highly magnetized accreting neutron stars is essential for interpreting their X-ray spectra and variability. This task is challenging due to the dynamic, multidimensional nature of the accretion columns and the strong gravitational field. Radiation propagation outside the column, within the neutron star's gravitational field, significantly shapes the observed emission and its variability with rotational phase. Geometry -- the location of the magnetic poles and the observer’s inclination -- is one of the key factors influencing the observables due to the anisotropy of the problem. In this work, we study how visibility effects in the gravitational field of the neutron star influence the direct emission from an accretion column, highlighting their importance when interpreting plasma parameters from observed spectra. Building on the physical model from the previous paper in this series, we investigate the directly observable X-ray flux from one and two columns with fan-beam wall emission, down-boosted by the bulk flow, focusing on the impact of geometry. We find that the observed flux strongly depends on the geometrical setup, which we illustrate across various observables. Particularly notable are special geometries in which a column on the far side of the neutron star aligns with the line of sight. In these cases, shadowing and strong light bending produce a narrow intense peak in pulse profiles at soft energies and a pronounced dip at higher energies (≳ 20, ). The combination of shadowing and height dependence of emission also leads to spectral softening. Flux-derived luminosities can be overestimated by up to a factor of 10. The fundamental cyclotron resonant scattering feature (CRSF) appears predominantly in emission across geometries, and we discuss possible causes. keV Our results demonstrate that geometry has a major impact on pulse profiles, phase-averaged spectra, phase–energy maps, and the anisotropy factor relevant for luminosity estimates. Future studies should include the reflected component in addition to the direct emission considered here, and account for a more accurate treatment of resonant redistribution in the CRSF-forming region.

Accretion-powered X-ray pulsars show complex spectral and timing behavior that reflects the extreme conditions near the neutron star. Modeling these observables remains challenging, particularly at high mass-accretion rates, dot M 17 , , ^ g s -1 , where the extended accretion column plays a role in shaping anisotropic emission seen by a remote observer. Here, we present a modular framework that links radiation processes inside the column to the emission observed at infinity. The framework combines existing models for continuum and cyclotron resonant scattering feature (CRSF) formation in sidewall column emission with ray tracing, thereby bridging the gap between the rest-frame emission and the observable flux. We model the radiation from the accretion column using continuum emission from the column walls, calculated in the regime of saturated Comptonization in the presence of a radiation-dominated shock. This emission then passes through a thin outer layer with a fast-moving bulk flow, where CRSFs, derived from Monte Carlo simulations, are imprinted. Finally, a geometry-independent ray-tracing code accounts for light bending in the Schwarzschild metric, incorporating the shape and location of the emitting regions and neutron star rotation to compute the observed phase- and energy-dependent flux. We present the height-dependent anisotropic emission in the column rest frame. For the chosen model, approximately 70% of the radiation is emitted within 1, of the base of the column, predominantly at large angles to the magnetic field. Relativistic boosting in the bulk flow affects the appearance of CRSFs in the spectra, leading to a strong dependence of the line locations on the angle relative to the magnetic field. The proposed theoretical framework allows a straightforward combination of internal emission models with ray tracing to study the observed flux. Here, we focus on the joint treatment of the two-dimensional column structure and CRSF formation in the neutron star rest frame. The observed phase and spectral variations of the flux are highly sensitive to the location of the emitting regions and the observer's viewing angle, which is explored in an accompanying paper.

We report on a pair of X-ray Multi-Mirror Mission (XMM-Newton) observations of the Supergiant X-ray binary 4U 1909+07, which were performed on 2021 October 3 and 8, respectively. We measure the neutron star rotation period in each observation to be ∼602.62 s. This continues a long spin-up trend that has persisted since 2001 where the neutron star spin period was found to be ∼604.66 s. In our timing analysis, we observe strong variations in the amplitude of the 1–10 keV pulse profile as a function of time, and for the first time, we find a low flux interval extending for a single pulse period in which pulsations are no longer detected. We interpret this low flux interval as a pulse dropout similar to those observed in Vela X-1 and GX 301-2, which were each explained by a low-density cavity in the wind driving the propeller effect. In our time-resolved spectral analysis, we observed the spectral continuum, which can be described as an absorbed power law modified by a high-energy cutoff, to significantly soften during the pulse-dropout phase. No evidence of an increasing absorption column density was found. The observed softening in 4U 1909+07 also supports an interpretation that the observed pulse dropout may be driven by the propeller effect, but the quasi-spherical settling accretion regime cannot be ruled out.

Context. X-ray spectra of accretion-powered X-ray pulsars can often be described using a power-law continuum with a high-energy cutoff, which might be further modified by additional spectral components. The Be X-ray binary system 4U 0115+63 is well known for having one of the highest numbers of detected harmonics of its cyclotron resonant scattering features (CRSFs), a pronounced spectral component known as the “10 keV feature”, and quasiperiodic oscillations (QPOs) with a period of about 500 s during outbursts. Aims. The changes in count rate by a factor of two during the ∼500 s QPOs allow us to probe the variation in the spectral components with flux. We study the “10 keV feature” in emission, aiming to disentangle it from the broadband continuum and CRSFs and investigate its origin. Methods. We focus on the flux-dependent behavior of the CRSF and its harmonics, and particularly the contribution of the “10 keV feature”, as seen in the flux-resolved analysis of two NuSTAR observations of the 2015 outburst. Results. Comparing the flux-resolved spectra of a given observation with the respective total dataset revealed a distinct change in overall spectral shape at the position of the “10 keV feature” but no comparable deviation at the energies of the harmonic CRSFs. The change associated with the “10 keV feature” does not seem to involve its centroid energy, which remains constant within a given observation. We find indications for an anticorrelation between the continuum flux and the ratio of the “10 keV feature” flux to the continuum flux within each observation. Conclusions. The analysis strengthens previous claims that the “10 keV feature” shows some independence from the remaining features. This result supports the interpretation that the “10 keV feature” has a different formation mechanism than the continuum emission, although its origin lies within the same physical environment.

Context. Accreting X-ray pulsars (XRPs) undergo different physical regimes depending on different mass accretion rates, with related changes in the emitted X-ray spectral properties. Recent observations have shown a dramatic change in the emission properties of this class of sources observed at low luminosity. Aims. We explore the timing and spectral properties of the XRP MAXI J0655-013 observed in the low-luminosity regime (about 5 × 10 ³³ erg s ⁻¹ ) to witness the corresponding spectral shape and pulse profiles. Methods. We employed recent XMM-Newton and NuSTAR pointed observations of the MAXI J0655-013 X-ray activity during the low-luminosity stage. We explored several spectral models to fit the data and test theoretical expectations of the dramatic transition of the spectral shape compared to the higher luminosity regime. We studied the pulsating nature of the source and found a precise timing solution. We explored the energy-resolved pulse profiles and the derived energy dependence of different pulsed fraction estimators (PF _minmax and PF _rms ). We also obtained NuSTAR pulsed fraction spectra (PFS) at different luminosity regimes. Results. The MAXI J0655-013 spectrum is well fit by a double Comptonization model, in agreement with recent observational results and theoretical expectations that explain the observed spectrum as being composed of two distinct bumps, each dominated by different polarization modes. We measured a spin period of 1081.86 ± 0.02 s, consistent with the source spinning up compared to previous observations, yielding an upper limit for the magnetic field strength of B ≲ 9 × 10 ¹³ G. The pulse profiles show a single broad peak interrupted by a sharp dip that coincides with an increase in the hardness ratio, and thus likely due to absorption. For the low-luminosity observation, the PF _minmax increases with energy up to ∼100% in the 10–30 keV band, while the PF _rms remains steady at ∼60%. The PFS obtained at high luminosity shows evidence of an iron K α emission line but no indications of a cyclotron line.