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 .

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.

A study of archival RXTE, Swift, and Suzaku pointed observations of the transient high mass X-ray binary GRO J1008-57 is presented. A new orbital ephemeris based on pulse arrival timing shows times of maximum luminosities during outbursts of GRO J1008-57 to be close to periastron at orbital phase -0.03. This makes the source one of a few for which outburst dates can be predicted with very high precision. Spectra of the source in 2005, 2007, and 2011 can be well described by a simple power law with high energy cutoff and an additional black body at lower energies. The photon index of the power law and the black body flux only depend on the 15-50 keV source flux. No apparent hysteresis effects are seen. These correlations allow to predict the evolution of the pulsar's X-ray spectral shape over all outbursts as a function of just one parameter, the source's flux. If modified by an additional soft component, this prediction even holds during GRO J1008-57's 2012 type II outburst.

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.

Reliable timing calibration is essential for the accurate comparison of XMM-Newton light curves with those from other observatories, to ultimately use them to derive precise physical quantities. The XMM-Newton timing calibration is based on pulsar analysis. However, as pulsars show both timing noise and glitches, it is essential to monitor these calibration sources regularly. To this end, the XMM-Newton observatory performs observations twice a year of the Crab pulsar to monitor the absolute timing accuracy of the EPIC-pn camera in the fast Timing and Burst modes. We present the results of this monitoring campaign, comparing XMM-Newton data from the Crab pulsar (PSR B0531+21) with radio measurements. In addition, we use five pulsars (PSR J0537-69, PSR B0540-69, PSR B0833-45, PSR B1509-58 and PSR B1055-52) with periods ranging from 16 ms to 197 ms to verify the relative timing accuracy. We analysed 38 XMM-Newton observations (0.2-12.0 keV) of the Crab taken over the first ten years of the mission and 13 observations from the five complementary pulsars. All the data were processed with the SAS, the XMM-Newton Scientific Analysis Software, version 9.0. Epoch folding techniques coupled with \chi^{2} tests were used to derive relative timing accuracies. The absolute timing accuracy was determined using the Crab data and comparing the time shift between the main X-ray and radio peaks in the phase folded light curves. The relative timing accuracy of XMM-Newton is found to be better than 10^{-8}. The strongest X-ray pulse peak precedes the corresponding radio peak by 306\pm9 \mus, which is in agreement with other high energy observatories such as Chandra, INTEGRAL and RXTE. The derived absolute timing accuracy from our analysis is \pm48 \mus.

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.

We present a detailed analysis of observations of the high mass X-ray binary Cen X-3 spanning two consecutive binary orbits performed with the RXTE satellite in early March 1997. The PCA and HEXTE light curves both show a clear reduction in count rate after mid-orbit for both binary revolutions. We therefore analyze two broad band spectra for each orbit, before and after mid-orbit. Consistent with earlier observations these four joint PCA and HEXTE spectra can be well described using a phenomenological pulsar continuum model, including an iron emission line and a cyclotron resonance scattering feature. While no strong spectral variations were detected, the second half of orbit 2 shows a tendency toward being softer and more strongly absorbed. In order to follow the orbital phase-dependent evolution of the spectrum in greater detail, we model spectra for shorter exposures, confirming that most spectral parameters show either a gradual or sudden change for the second half of the second orbit. A comparison with a simple wind model indicates the existence of an accretion wake in this system. We also present and discuss high resolution pulse profiles for several different energy bands, as well as their hardness ratios. PCA and HEXTE spectra were created for 24 phase bins and fitted using the same model as in the phase averaged case. Systematic pulse phase-dependent variations of several continuum and cyclotron line parameters were detected, most notably a significant increase of the cyclotron line energy during the early rise of the main peak, followed by a gradual decrease. We show that applying a simple dipole model for the magnetic field is not sufficient to describe our data.

(abbreviated version of the abstract) We study the physics of cyclotron line formation in the high-energy spectra of accreting X-ray pulsars using Monte Carlo methods, assuming that the line-forming region is a low-density electron plasma in a sub-critical magnetic field. We investigate the dependence of the shape of the fundamental line on angle, geometry, optical depth and temperature. We also discuss variations of the line ratios for non-uniform magnetic fields. These numerical predictions for the line profiles are linked to results from observational data analysis using an XSPEC model based on the Monte Carlo simulations. We apply this model to observational data from RXTE and INTEGRAL. The predicted strong emission wings of the fundamental cyclotron feature are not found in observational data, hinting at a bottom illuminated slab geometry for line formation.

The Be/X-ray binary 3A 0535+262 has the highest magnetic field determined by cyclotron line studies of all accreting X-ray pulsars, despite an open debate if the fundamental line was rather at ~50 or above 100 keV as observed by different instruments in past outbursts. The source went into quiescence for more than ten years since its last outbursts in 1994. Observing during a `normal' outburst August/September 2005 with Integral and RXTE we find a strong cyclotron line feature at ~45 keV and have for the first time since 1975 determined the low energy pulse profile.