The X-ray Spectra of Accreting Pulsars: Studies of Three Sources Using Empirical and Phenomenological Models

Accreting X-ray pulsars are a class of astrophysical objects consisting of a neutron star in a binary system with a stellar companion. Matter expelled by the companion star is captured by the neutron star's gravity; as this matter falls towards the neutron star's surface, is is compressed and heated, giving off X-rays. As the matter falls the last few miles above the neutron star surface, a number of physical processes compete for dominance, resulting in a highly complex environment governed by the interplay of magnetic, hydrodynamical, and radiative processes. The resulting spectrum often shows broad absorption-like features called cyclotron lines, which provide the only direct measurement of the magnetic field of a neutron star and act as probes of the properties of the accretion column, and their behavior with respect to changes in the accretion rate onto the neutron star has been of interest in recent years. My work in this dissertation brings together nearly 20 years of data from three X-ray satellites to study the X-ray emission from accreting pulsars, with a focus on the hard X-ray continuum and cyclotron lines. I present results for the accreting pulsars 4U 1538-522 and 4U 1907+09, examining the behavior of their cyclotron lines with respect to their luminosity, finding evidence for a positive correlation between the line energy and luminosity in 4U 1907+09. A combined analysis of most of the available X-ray data for the accreting pulsar 4U 1538-522 shows no such correlation in this source, either positive or negative. However, I do present evidence that the cyclotron line energy in 4U 1538-522 has shifted upwards by ~5% in recent years compared to measurements from 10-20 years ago. I additionally carry out an extensive analysis of the environment around 4U 1538-522 using the soft X-ray detectors aboard the satellite Suzaku. I finally present a set of new results from the transient X-ray pulsar V 0332+53, which I fit with a new physics-based model for the accretion column. These fits allow me to constrain the size and temperature of the accretion column, as well as the relative contributions of different processes in the column to the overall observed spectrum