We present an analysis of the spectral shape and pulse profile of the accretion-powered pulsar 4U 1626-67 observed with Suzaku and NuSTAR during a spin-up state. The pulsar, which experienced a torque reversal to spin-up in 2008, has a spin period of 7.7 s. Comparing the phase-averaged spectra obtained with Suzaku in 2010 and with NuSTAR in 2015, we find that the spectral shape changed between the two observations: the 3-10 keV flux increased by 5% while the 30-60 keV flux decreased significantly by 35%. Phase-averaged and phase-resolved spectral analysis shows that the continuum spectrum observed by NuSTAR is well described by an empirical NPEX continuum with an added broad Gaussian emission component around the spectral peak at 20 keV. Taken together with the observed Pdot value obtained from Fermi/GBM, we conclude that the spectral change between the Suzaku and NuSTAR observations was likely caused by an increase of the accretion rate. We also report the possible detection of asymmetry in the profile of the fundamental cyclotron line. Furthermore, we present a study of the energy-resolved pulse profiles using a new relativistic ray tracing code, where we perform a simultaneous fit to the pulse profiles assuming a two-column geometry with a mixed pencil- and fan-beam emission pattern. The resulting pulse profile decompositions enable us to obtain geometrical parameters of accretion columns (inclination, azimuthal and polar angles) and a fiducial set of beam patterns. This information is important to validate the theoretical predictions from radiation transfer in a strong magnetic field.

Accreting X-ray pulsars are among the best observed objects of X-ray astronomy with a rich data set of observational phenomena in the spectral and timing domain. While the general picture for these sources is well established, the detailed physics behind the observed phenomena are often subject of debate. We present recent observational, theoretical and modeling results for the structure and dynamics of the accretion column in these systems. Our results indicate the presence of different accretion regimes and possible explanations for observed variations of spectral features with luminosity.

Accretion-powered X-ray pulsars exhibit significant variability of the Cyclotron Resonance Scattering Feature (CRSF) centroid energy on pulse-to-pulse timescales, and also on much longer timescales. Two types of spectral variability are observed. For sources in group 1, the CRSF energy is negatively correlated with the variable source luminosity, and for sources in group 2, the opposite behavior is observed. The physical basis for this bimodal behavior is currently not understood. We explore the hypothesis that the accretion dynamics in the group 1 sources is dominated by radiation pressure near the stellar surface, and that Coulomb interactions decelerate the gas to rest in the group 2 sources. We derive a new expression for the critical luminosity such that radiation pressure decelerates the matter to rest in the supercritical sources. The formula for the critical luminosity is evaluated for 5 sources, using the maximum value of the CRSF centroid energy to estimate the surface magnetic field strength. The results confirm that the group 1 sources are supercritical and the group 2 sources are subcritical, although the situation is less clear for those highly variable sources that cross over the critical line. We also explain the variation of the CRSF energy with luminosity as a consequence of the variation of the characteristic emission height. The sign of the height variation is opposite in the supercritical and subcritical cases, hence creating the observed bimodal behavior.