We give an overview of the SImulation of X-ray TElescopes (SIXTE) software package, a generic, mission-independent Monte Carlo simulation toolkit for X-ray astronomical instrumentation. The package is based on a modular approach for the source definition, the description of the optics, and the detector type such that new missions can be easily implemented. The targets to be simulated are stored in a flexible input format called SIMPUT. Based on this source definition, a sample of photons is produced and then propagated through the optics. In order to model the detection process, the software toolkit contains modules for various detector types, ranging from proportional counter and Si-based detectors, to more complex descriptions like transition edge sensor (TES) devices. The implementation of characteristic detector effects and a detailed modeling of the read-out process allow for representative simulations and therefore enable the analysis of characteristic features, such as for example pile-up, and their impact on observations. We present an overview of the implementation of SIXTE from the input source, the imaging, and the detection process, highlighting the modular approach taken by the SIXTE software package. In order to demonstrate the capabilities of the simulation software, we present a selection of representative applications, including the all-sky survey of eROSITA and a study of pile-up effects comparing the currently operating XMM-Newton with the planned Athena-WFI instrument. A simulation of a galaxy cluster with the Athena- X-IFU shows the capability of SIXTE to predict the expected performance of an observation for a complex source with a spatially varying spectrum and our current knowledge of the future instrument.

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 present two observations of the high-mass X-ray binary GX 301-2 with NuSTAR, taken at different orbital phases and different luminosities. We find that the continuum is well described by typical phenomenological models, like a very strongly absorbed NPEX model. However, for a statistically acceptable description of the hard X-ray spectrum we require two cyclotron resonant scattering features (CRSF), one at ~35 keV and the other at ~50 keV. Even though both features strongly overlap, the good resolution and sensitivity of NuSTAR allows us to disentangle them at >=99.9% significance. This is the first time that two CRSFs are seen in GX 301-2. We find that the CRSFs are very likely independently formed, as their energies are not harmonically related and, if it were a single line, the deviation from a Gaussian shape would be very large. We compare our results to archival Suzaku data and find that our model also provides a good fit to those data. We study the behavior of the continuum as well as the CRSF parameters as function of pulse phase in seven phase bins. We find that the energy of the 35 keV CRSF varies smoothly as function of phase, between 30-38 keV. To explain this variation, we apply a simple model of the accretion column, taking the altitude of the line-forming region, the velocity of the in-falling material, and the resulting relativistic effects into account. We find that in this model the observed energy variation can be explained simply due to a variation of the projected velocity and beaming factor of the line forming region towards us.

We present a comprehensive spectral analysis of the BeXRB GRO J1008-57 over a luminosity range of three orders of magnitude using NuSTAR, Suzaku and RXTE data. We find significant evolution of the spectral parameters with luminosity. In particular the photon index hardens with increasing luminosity at intermediate luminosities between $10^{36}$ $-$ $10^{37}$ erg s$^{-1}$. This evolution is stable and repeatedly observed over different outbursts. However, at the extreme ends of the observed luminosity range, we find that the correlation breaks down, with a significance level of at least $3.7\sigma$. We conclude that these changes indicate transitions to different accretion regimes, which are characterized by different deceleration processes, such as Coulomb or radiation breaking. We compare our observed luminosity levels of these transitions to theoretical predications and discuss the variation of those theoretical luminosity values with fundamental neutron star parameters. Finally, we present detailed spectroscopy of the unique "triple peaked" outburst in 2014/15 which does not fit in the general parameter evolution with luminosity. The pulse profile on the other hand is consistent with what is expected at this luminosity level, arguing against a change in accretion geometry. In summary, GRO J1008-57 is an ideal target to study different accretion regimes due to the well constrained evolution of its broad-band spectral continuum over several orders of magnitude in luminosity.

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.

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.

Simultaneous broadband spectral and temporal studies of blazars are an important tool for investigating active galactic nuclei (AGN) jet physics. We study the spectral evolution between quiescent and flaring periods of 22 radio-loud AGN through multi-epoch, quasi-simultaneous broadband spectra. For many of these sources these are the first broadband studies. We use a Bayesian block analysis of \Fermi/LAT light curves in order to determine time ranges of constant flux for constructing quasi-simultaneous SEDs. The shapes of the resulting 81 SEDs are described by two logarithmic parabolas and a blackbody spectrum where needed. For low states the peak frequencies and luminosities agree well with the blazar sequence, higher luminosity implying lower peak frequencies. This is not true for sources in a high state. The $\gamma$-ray photon index in Fermi/LAT correlates with the synchrotron peak frequency in low and intermediate states. No correlation is present in high states. The black hole mass cannot be determined from the SEDs. Surprisingly, the thermal excess often found in FSRQs at optical/UV wavelengths can be described by blackbody emission and not an accretion disk spectrum. The "harder-when-brighter" trend, typically seen in X-ray spectra of flaring blazars, is visible in the blazar sequence. Our results for low and intermediate states, as well as the Compton dominance, are in agreement with previous results. Black hole mass estimates using the parameters from Bonchi (2013) are in agreement with some of the more direct measurements. For two sources, estimates disagree by more than four orders of magnitude, possibly due to boosting effects. The shapes of the thermal excess seen predominantly in flat spectrum radio quasars are inconsistent with a direct accretion disk origin.

We report on the timing and spectral analysis of two Suzaku observations with different flux levels of the high-mass X-ray binary KS 1947+300 during its 2013 outburst. In agreement with simultaneous NuSTAR observations, the continuum is well described by an absorbed power law with a cut-off and an additional black body component. In addition we find fluorescent emission from neutral, He-like, and even H-like iron. We determine a pulse period of ~18.8 s with the source showing a spin-up between the two observations. Both Suzaku observations show a very similar behavior of the pulse profile, which is strongly energy dependent, with an evolution from a profile with one peak at low energies to a profile with two peaks of different widths towards higher energies seen in both, the Suzaku and NuSTAR data. Such an evolution to a more complex profile at higher energies is rarely seen in X-ray pulsars, most cases show the opposite behavior. Pulse phase-resolved spectral analysis shows a variation in the absorbing column density, NH , over pulse phase. Spectra taken during the pulse profile minima are intrinsically softer compared to the pulse phase-averaged spectrum.

We present a timing and spectral analysis of the X-ray pulsar XTE J1946+274 observed with Suzaku during an outburst decline in 2010 October and compare with previous results. XTE J1946+274 is a transient X-ray binary consisting of a Be-type star and a neutron star with a 15.75 s pulse period in a 172 d orbit with 2-3 outbursts per orbit during phases of activity. We improve the orbital solution using data from multiple instruments. The X-ray spectrum can be described by an absorbed Fermi-Dirac cutoff power law model along with a narrow Fe K line at 6.4 keV and a weak Cyclotron Resonance Scattering Feature (CRSF) at ~35 keV. The Suzaku data are consistent with the previously observed continuum flux versus iron line flux correlation expected from fluorescence emission along the line of sight. However, the observed iron line flux is slightly higher, indicating the possibility of a higher iron abundance or the presence of non-uniform material. We argue that the source most likely has only been observed in the subcritical (non-radiation dominated) state since its pulse profile is stable over all observed luminosities and the energy of the CRSF is approximately the same at the highest (~$5 \times 10^{37} $erg s$^{-1}$) and lowest (~$5 \times 10^{36} $erg s$^{-1}$) observed 3-60 keV luminosities.

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.