For more than 40 years, most astrophysical observations and laboratory studies of two key soft x-ray diagnostic $2p-3d$ transitions, $3C$ and $3D$, in Fe XVII ions found oscillator strength ratios $f(3C)/f(3D)$ disagreeing with theory, but uncertainties had precluded definitive statements on this much studied conundrum. Here, we resonantly excite these lines using synchrotron radiation at PETRA III, and reach, at a millionfold lower photon intensities, a 10 times higher spectral resolution, and 3 times smaller uncertainty than earlier work. Our final result of $f(3C)/f(3D) = 3.09(8)(6)$ supports many of the earlier clean astrophysical and laboratory observations, while departing by five sigmas from our own newest large-scale ab initio calculations, and excluding all proposed explanations, including those invoking nonlinear effects and population transfers.

We demonstrate a widely applicable technique to absolutely calibrate the energy scale of x-ray spectra with experimentally well-known and accurately calculable transitions of highly charged ions, allowing us to measure the K-shell Rydberg spectrum of molecular O$_2$ with 6 meV uncertainty. We reveal a systematic $\sim$450 meV shift from previous literature values, and settle an extraordinary discrepancy between astrophysical and laboratory measurements of neutral atomic oxygen, the latter being calibrated against the aforementioned O$_2$ literature values. Because of the widespread use of such, now deprecated, references, our method impacts on many branches of x-ray absorption spectroscopy. Moreover, it potentially reduces absolute uncertainties there to below the meV level.

We resonantly excite the $K$ series of O$^{5+}$ and O$^{6+}$ up to principal quantum number $n=11$ with monochromatic x rays, producing $K$-shell holes, and observe their relaxation by soft x-ray emission. Some photoabsorption resonances of O$^{5+}$ reveal strong two-electron-one-photon (TEOP) transitions. We find that for the $1s\,2s\,5p\, ^1P_1$ state, TEOP relaxation is by far stronger than the competing, and usually much faster, allowed Auger and radiative decay paths. This enhanced TEOP decay arises from a strong correlation with the near-degenerate upper state $1s\,2p\,4s\, ^1P_1$ of a $K\alpha$ satellite. Even in three-electron systems, TEOP transitions can play a dominant role, and the present results should guide further research on the ubiquitous and abundant many-electron ions where various electronic energy degeneracies and level mixings are far more common.

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 three Chandra High Energy Transmission Gratings observations of the black hole binary Cyg X-1/HDE 226868 at different orbital phases. The stellar wind that is powering the accretion in this system is characterized by temperature and density inhomogeneities including structures, or "clumps", of colder, more dense material embedded in the photoionized gas. As these clumps pass our line of sight, absorption dips appear in the light curve. We characterize the properties of the clumps through spectral changes during various dip stages. Comparing the silicon and sulfur absorption line regions (1.6-2.7 keV $\equiv$ 7.7-4.6 {\AA}) in four levels of varying column depth reveals the presence of lower ionization stages, i.e., colder or denser material, in the deeper dip phases. The Doppler velocities of the lines are roughly consistent within each observation, varying with the respective orbital phase. This is consistent with the picture of a structure that consists of differently ionized material, in which shells of material facing the black hole shield the inner and back shells from the ionizing radiation. The variation of the Doppler velocities compared to a toy model of the stellar wind, however, does not allow us to pin down an exact location of the clump region in the system. This result, as well as the asymmetric shape of the observed lines, point at a picture of a complex wind structure.

The high mass X-ray binary Cep X-4, during its 2014 outburst, showed evidence for an asymmetric cyclotron line in its hard X-ray spectrum. The 2014 spectrum provides one of the clearest cases of an asymmetric line profile among all studied sources with Cyclotron Resonance Scattering Features (CRSF). We present a phase-resolved analysis of NuSTAR and Suzaku data taken at the peak and during the decline phases of this outburst. We find that the pulse-phased resolved spectra are well-fit by a single, symmetric cyclotron feature. The fit parameters vary strongly with pulse phase: most notably the central energy and depth of the cyclotron feature, the slope of the power-law component, and the absorbing column density. We synthesise a phase averaged spectrum using the best-fit parameters for these individual pulse phases, and find that this combined model spectrum has a similar asymmetry in the cyclotron features as discovered in phase-averaged data. We conclude that the pulse phase resolved analysis with simple symmetric line profiles when combined can explain the asymmetry detected in the phase-averaged data.

Two \textit{NuSTAR} observations of the luminous X-ray pulsar LMC X-4 in October and November 2015 captured several bright accretion flares from this source, which has a long history of stable pulse and superorbital behavior. We present a timing analysis of these data in which we detect a rapid pulse turn-on in association with the accretion flares, during which the source reaches super-Eddington luminosities. Pulsations, which are normally seen from this source, are found to only occur for approximately one hour before and during the bright flares. Beyond one hour before and after the flares, we find pulsations to be weak or nonexistent, with fractional RMS amplitudes of less than 0.05. At the onset of the flare, the pulse profiles exhibit a phase shift of 0.25 cycles that could be associated with a change in the emission geometry. This increase in pulse strength occurring well before the flare cannot be explained by the propeller effect, and potentially offers a connection between the magnetic properties of pulsars that accrete close to their Eddington limits and ULX pulsars.

We present an X-ray spectral and timing analysis of two $NuSTAR$ observations of the transient Be X-ray binary SAX J2103.5+4545 during its April 2016 outburst, which was characterized by the highest flux since $NuSTAR$'s launch. These observations provide detailed hard X-ray spectra of this source during its bright precursor flare and subsequent fainter regular outburst for the first time. In this work, we model the phase-averaged spectra for these observations with a negative and positive power law with an exponential cut-off (NPEX) model and compare the pulse profiles at different flux states. We found that the broad-band pulse profile changes from a three peaked pulse in the first observation to a two peaked pulse in the second observation, and that each of the pulse peaks has some energy dependence. We also perform pulse-phase spectroscopy and fit phase-resolved spectra with NPEX to evaluate how spectral parameters change with pulse phase. We find that while the continuum parameters are mostly constant with pulse phase, a weak absorption feature at ~12 keV that might, with further study, be classified as a cyclotron line, does show strong pulse phase dependence.

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.

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.

We report on the RXTE detection of a sudden increase in the absorption column density, $N_\mathrm{H}$, during the 2011 May outburst of GX 304-1. The $N_\mathrm{H}$ increased up to ${\sim}16\times 10^{22}$ atoms cm$^{-2}$, which is a factor of 3-4 larger than what is usually measured during the outbursts of GX 304-1 as covered by RXTE. Additionally, an increase in the variability of the hardness ratio as calculated from the energy resolved RXTE-PCA light curves is measured during this time range. We interpret these facts as an occultation event of the neutron star by material in the line of sight. Using a simple 3D model of an inclined and precessing Be disk around the Be type companion, we are able to qualitatively explain the $N_\mathrm{H}$ evolution over time. We are able to constrain the Be-disk density to be on the order of $10^{-11}$ g cm$^{-3}$. Our model strengthens the idea of inclined Be disks as origin of double-peaked outbursts as the derived geometry allows accretion twice per orbit under certain conditions.

We report the identification of a bright hard X-ray source dominating the M31 bulge above 25 keV from a simultaneous NuSTAR-Swift observation. We find that this source is the counterpart to Swift J0042.6+4112, which was previously detected in the Swift BAT All-sky Hard X-ray Survey. This Swift BAT source had been suggested to be the combined emission from a number of point sources; our new observations have identified a single X-ray source from 0.5 to 50 keV as the counterpart for the first time. In the 0.5-10 keV band, the source had been classified as an X-ray binary candidate in various Chandra and XMM studies; however, since it was not clearly associated with Swift J0042.6+4112, the previous E < 10 keV observations did not generate much attention. This source has a spectrum with a soft X-ray excess (kT~ 0.2 keV) plus a hard spectrum with a power law of Gamma ~ 1 and a cutoff around 15-20 keV, typical of the spectral characteristics of accreting pulsars. Unfortunately, any potential pulsation was undetected in the NuSTAR data, possibly due to insufficient photon statistics. The existing deep HST images exclude high-mass (>3 Msun) donors at the location of this source. The best interpretation for the nature of this source is an X-ray pulsar with an intermediate-mass (<3 Msun) companion or a symbiotic X-ray binary. We discuss other possibilities in more detail.

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 .

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.

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.