We present an analysis of X-ray spectra of the high mass X-ray binary 4U 0115+634 as observed with Suzaku and RXTE in 2011 July, during the fading phase of a giant X-ray outburst. We used a continuum model consisting of an absorbed cutoff power-law and an ad-hoc Gaussian emission feature centered around 8.5 keV, which we discuss to be due to cyclotron emission. Our results are consistent with a fundamental cyclotron absorption line centered at ${\sim}10.2$ keV for all observed flux ranges. At the same time we rule out significant influence of the 8.5 kev Gaussian on the CRSF parameters, which are not consistent with the cyclotron line energies and depths of previously reported flux-dependent descriptions. We also show that some continuum models can lead to artificial line-like residuals in the analyzed spectra, which are then misinterpreted as unphysically strong cyclotron lines. Specifically, our results do not support the existence of a previously claimed additional cyclotron feature at ${\sim}15$ keV. Apart from these features, we find for the first time evidence for a He-like Fe XXV emission line at ${\sim}6.7$ keV and weak H-like Fe XXVI emission close to ${\sim}7.0$ keV.

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.

We report on two NuSTAR observations of the HMXB A 0535+26 taken toward the end of its normal 2015 outburst at very low $3-50$ keV luminosities of ${\sim}1.4\times10^{36}$ erg/s and ${\sim}5\times10^{35}$ erg/s which are complemented by 9 Swift observations. The data clearly confirm indications seen in earlier data that the source's spectral shape softens as it becomes fainter. The smooth, exponential rollover at high energies present in the first observation evolves to a much more abrupt steepening of the spectrum at $20-30$ keV. The continuum evolution can be well described with emission from a magnetized accretion column, modeled using the compmag model modified by an additional Gaussian emission component for the fainter observation. Between the two observations, the optical depth changes from $0.75\pm0.04$ to $0.56^{+0.01}_{-0.04}$, the electron temperature remains constant, and there is an indication that the column decreases in radius. Since the energy resolved pulse profiles remain virtually unchanged in shape between the two observations, the emission properties of the accretion column, however, reflect the same accretion regime. This conclusion is also confirmed by our result that the energy of the cyclotron resonant scattering feature (CRSF) at ${\sim}45$ keV is independent of the luminosity, implying that the magnetic field in the region in which the observed radiation is produced is the same in both observations. Finally, we also constrain the evolution of the continuum parameters with rotational phase of the neutron star. The width of the CRSF could only be constrained for the brighter observation. Based on Monte-Carlo simulations of CRSF formation in single accretion columns, its pulse phase dependence supports a simplified fan beam emission pattern. The evolution of the CRSF width is very similar to that of the CRSF depth, which is in disagreement with expectations.

Cyclotron resonant scattering features (CRSFs) are formed by scattering of X-ray photons off quantized plasma electrons in the strong magnetic field (of the order 10^12 G) close to the surface of an accreting X-ray pulsar. The line profiles of CRSFs cannot be described by an analytic expression. Numerical methods such as Monte Carlo (MC) simulations of the scattering processes are required in order to predict precise line shapes for a given physical setup, which can be compared to observations to gain information about the underlying physics in these systems. A versatile simulation code is needed for the generation of synthetic cyclotron lines. Sophisticated geometries should be investigatable by making their simulation possible for the first time. The simulation utilizes the mean free path tables described in the first paper of this series for the fast interpolation of propagation lengths. The code is parallelized to make the very time consuming simulations possible on convenient time scales. Furthermore, it can generate responses to mono-energetic photon injections, producing Green's functions, which can be used later to generate spectra for arbitrary continua. We develop a new simulation code to generate synthetic cyclotron lines for complex scenarios, allowing for unprecedented physical interpretation of the observed data. An associated XSPEC model implementation is used to fit synthetic line profiles to NuSTAR data of Cep X-4. The code has been developed with the main goal of overcoming previous geometrical constraints in MC simulations of CRSFs. By applying this code also to more simple, classic geometries used in previous works, we furthermore address issues of code verification and cross-comparison of various models. The XSPEC model and the Green's function tables are available online at http://www.sternwarte.uni-erlangen.de/research/cyclo .

Electron cyclotron resonant scattering features (CRSFs) are observed as absorption-like lines in the spectra of X-ray pulsars. A significant fraction of the computing time for Monte Carlo simulations of these quantum mechanical features is spent on the calculation of the mean free path for each individual photon before scattering, since it involves a complex numerical integration over the scattering cross section and the (thermal) velocity distribution of the scattering electrons. We aim to numerically calculate interpolation tables which can be used in CRSF simulations to sample the mean free path of the scattering photon and the momentum of the scattering electron. The tables also contain all the information required for sampling the scattering electron's final spin. The tables were calculated using an adaptive Simpson integration scheme. The energy and angle grids were refined until a prescribed accuracy is reached. The tables are used by our simulation code to produce artificial CRSF spectra. The electron momenta sampled during these simulations were analyzed and justified using theoretically determined boundaries. We present a complete set of tables suited for mean free path calculations of Monte Carlo simulations of the cyclotron scattering process for conditions expected in typical X-ray pulsar accretion columns (0.01<B/B_{crit}<=0.12, where B_{crit}=4.413x10^{13} G and 3keV<=kT<15keV). The sampling of the tables is chosen such that the results have an estimated relative error of at most 1/15 for all points in the grid. The tables are available online at http://www.sternwarte.uni-erlangen.de/research/cyclo.

We report new spectral modeling of the accreting X-ray pulsar Hercules X- 1. Our radiation-dominated radiative shock model is an implementation of the analytic work of Becker & Wolff on Comptonized accretion flows onto magnetic neutron stars. We obtain a good fit to the spin-phase averaged 4 to 78 keV X-ray spectrum observed by the Nuclear Spectroscopic Telescope Array during a main- on phase of the Her X-1 35-day accretion disk precession period. This model allows us to estimate the accretion rate, the Comptonizing temperature of the radiating plasma, the radius of the magnetic polar cap, and the average scattering opacity parameters in the accretion column. This is in contrast to previous phenomenological models that characterized the shape of the X-ray spectrum but could not determine the physical parameters of the accretion flow. We describe the spectral fitting details and discuss the interpretation of the accretion flow physical parameters.

In a previous work, we introduced a tool for analyzing multiple datasets simultaneously, which has been implemented into ISIS. This tool was used to fit many spectra of X-ray binaries. However, the large number of degrees of freedom and individual datasets raise an issue about a good measure for a simultaneous fit quality. We present three ways to check the goodness of these fits: we investigate the goodness of each fit in all datasets, we define a combined goodness exploiting the logical structure of a simultaneous fit, and we stack the fit residuals of all datasets to detect weak features. These tools are applied to all RXTE-spectra from GRO 1008-57, revealing calibration features that are not detected significantly in any single spectrum. Stacking the residuals from the best-fit model for the Vela X-1 and XTE J1859+083 data evidences fluorescent emission lines that would have gone undetected otherwise.

Parallel computing and steadily increasing computation speed have led to a new tool for analyzing multiple datasets and datatypes: fitting several datasets simultaneously. With this technique, physically connected parameters of individual data can be treated as a single parameter by implementing this connection into the fit directly. We discuss the terminology, implementation, and possible issues of simultaneous fits based on the X-ray data analysis tool Interactive Spectral Interpretation System (ISIS). While all data modeling tools in X-ray astronomy allow in principle fitting data from multiple data sets individually, the syntax used in these tools is not often well suited for this task. Applying simultaneous fits to the transient X-ray binary GRO J1008-57, we find that the spectral shape is only dependent on X-ray flux. We determine time independent parameters such as, e.g., the folding energy E_fold, with unprecedented precision.

Context: Accretion-powered X-ray pulsars show highly energy-dependent and complex pulse-profile morphologies. Significant deviations from the average pulse profile can appear, in particular close to the cyclotron line energies. These deviations can be described as energy-dependent phase lags, that is, as energy-dependent shifts of main features in the pulse profile. Aims: Using a numerical study we explore the effect of cyclotron resonant scattering on observable, energy-resolved pulse profiles. Methods: We generated the observable emission as a function of spin phase, using Monte Carlo simulations for cyclotron resonant scattering and a numerical ray-tracing routine accounting for general relativistic light-bending effects on the intrinsic emission from the accretion columns. Results: We find strong changes in the pulse profile coincident with the cyclotron line energies. Features in the pulse profile vary strongly with respect to the average pulse profile with the observing geometry and shift and smear out in energy additionally when assuming a non-static plasma. Conclusions: We demonstrate how phase lags at the cyclotron energies arise as a consequence of the effects of angular redistribution of X-rays by cyclotron resonance scattering in a strong magnetic field combined with relativistic effects. We also show that phase lags are strongly dependent on the accretion geometry. These intrinsic effects will in principle allow us to constrain a system's accretion geometry.

We report on an outburst of the high mass X-ray binary 4U 0115+63 with a pulse period of 3.6s in 2008 March/April as observed with RXTE and INTEGRAL. During the outburst the neutron star's luminosity varied by a factor of 10 in the 3--50\,keV band. In agreement with earlier work we find evidence for five cyclotron resonance scattering features at ~10.7, 21.8, 35.5, 46.7, and 59.7keV. Previous work had found an anti-correlation between the fundamental cyclotron line energy and the X-ray flux. We show that this apparent anti-correlation is probably due to the unphysical interplay of parameters of the cyclotron line with the continuum models used previously, e.g., the negative and positive exponent power law (NPEX). For this model, we show that cyclotron line modeling erroneously leads to describing part of the exponential cutoff and the continuum variability, and not the cyclotron lines. When the X-ray continuum is modeled with a simple exponentially cutoff power law modified by a Gaussian emission feature around 10keV, the correlation between the line energy and the flux vanishes and the line parameters remain virtually constant over the outburst. We therefore conclude that the previously reported anti-correlation is an artifact of the assumptions adopted in the modeling of the continuum.

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.

Seventeen accreting neutron star pulsars, mostly high mass X-ray binaries with half of them Be-type transients, are known to exhibit Cyclotron Resonance Scattering Features (CRSFs) in their X-ray spectra, with characteristic line energies from 10 to 60 keV. To date about two thirds of them, plus a few similar systems without known CRSFs, have been observed with Suzaku. We present an overview of results from these observations, including the discovery of a CRSF in the transient 1A 1118-61 and pulse phase resolved spectroscopy of GX 301-2. These observations allow for the determination of cyclotron line parameters to an unprecedented degree of accuracy within a moderate amount of observing time. This is important since these parameters vary - e.g., with orbital phase, pulse phase, or luminosity - depending on the geometry of the magnetic field of the pulsar and the properties of the accretion column at the magnetic poles. We briefly introduce a spectral model for CRSFs that is currently being developed and that for the first time is based on these physical properties. In addition to cyclotron line measurements, selected highlights from the Suzaku analyses include dip and flare studies, e.g., of 4U 1907+09 and Vela X-1, which show clumpy wind effects (like partial absorption and/or a decrease in the mass accretion rate supplied by the wind) and may also display magnetospheric gating effects.