Abstract Cyclotron resonant scattering features in accreting X-ray pulsars are often difficult to detect, especially when shallow or variable. Recent studies have shown that combining spectral and timing analyses enhances their detectability. We investigated the evolution of energy-resolved pulse profiles of the X-ray pulsar 4U 1901+03 during its 2019 giant outburst, focusing on the 30–40 keV range where there have been disputed claims of a cyclotron line detection and its properties.We analysed four NuSTAR observations of 4U 1901+03 at different luminosities. We studied energy-resolved pulse profiles using harmonic decomposition, cross-correlation analysis, energy–phase maps, and pulsed-fraction spectra, and we used Bayesian spectral modelling to assess the presence and properties of a cyclotron line. We detected significant spectral–timing variability in the 30–40 keV range, becoming stronger at lower luminosities. We found a pronounced drop in the pulsed fraction near ≈ 35 keV only in the lowest accretion state and in the first harmonic of one intermediate-luminosity observation. Adopting a Bayesian informative approach on the spectral parameters, we find evidence for acyclotron line in all the examined energy spectra, with an average centroid energy of Ecyc ≈ 32 keV, varying by only ≈ 1.6%, with anti-correlation between line depth and luminosity.We show that a spectral-timing combined approach is more sensitive than phase-averaged spectroscopy to shallow cyclotron features. The luminosity-dependent evolution of pulse profiles and cyclotron line depth point to a drastic change in the emission geometry and accretion flow structure. Additional information The current repository contains cornerplots of the posterior distribution of the selected model through bayesian spectral fitting and the raw NuSTAR data analyzed, which can be downloaded also via HEASARC.
Ambrosi et al. (Tue,) studied this question.