How the spin-phase variability of cyclotron lines shapes the pulsed fraction spectra: insights from 4U 1538-52

TL;DR

This study investigates how the phase-dependent variability of cyclotron lines influences the energy-dependent pulse profiles of 4U 1538-52, revealing a broad bump near the cyclotron energy driven by spectral changes.

Contribution

The paper introduces a detailed physical model linking phase-dependent spectral variability, especially cyclotron line depth, to pulse profile features in 4U 1538-52, supported by simulations matching observations.

Findings

Broad bump near cyclotron energy in pulse profiles

Phase-dependent spectral variability drives pulse profile features

Modeling indicates high inclination and magnetic obliquity

Abstract

We study the energy-dependent pulse profile of 4U 1538-52 and its phase-dependent spectral variability, with emphasis on the behavior around the cyclotron resonant scattering feature at around 21 keV. We analyze all available NuSTAR observations of 4U 1538-52. We decompose energy-resolved pulse profiles into Fourier harmonics to study their energy dependence. Specifically, we compute pulsed fraction spectra, cross-correlation and lag spectra, identifying discontinuities and linking them to features in the phase-averaged spectra. We perform phase-averaged and phase-resolved spectral analyses to probe spectral variability and its relation to pulse profile changes. Finally, we interpret our findings via physical modeling of energy- and angle-dependent pulse profile emission, performing radiative transfer in a homogeneous slab-like atmosphere under conditions relevant to 4U 1538-52. The emission is projected onto the observer's sky plane to derive expected observables. In contrast to the dips in pulsed fraction spectra observed in other sources (e.g., Her X-1), we find a broad bump near the cyclotron resonance energy in 4U 1538-52. This increase is driven primarily by phase-dependent spectral variability, especially by strong variations in cyclotron line depth across different phase intervals. We interpret the observed contrast between dips and bumps in various sources as arising from phase-dependent variations of cyclotron line depth relative to the phase-modulated flux. We model the X-ray emission from an accreting neutron star and find that our simulations indicate high values of both the observer's inclination and the magnetic obliquity, along with a 10-15 degrees asymmetry between the locations of the magnetic poles. Assuming this geometry, we satisfactorily reproduce the observed pulse profiles and introduce general trends in the observables resulting from the system's geometry.