Engineering Papers⌕ Search

Engineering topics

Ralf Ballhausen

Publications and source records attributed to Ralf Ballhausen.

Xstar Atomic Database: The PyXstar Package

We present a progress report on the development of PyXstar, a Python package to manage the data (input, output, intermediate, atomic database, and model-grids) associated with the XSTAR code for treating photoionized and collisionally ionized plasmas. The PyXstar modular structure and database retrieval scheme are described, and its functionality is illustrated with Python functions and classes for performing database searches. We briefly compare PyXstar with two other Python spectrum modeling tools: PyNeb and PyAtomDB.

Claudio Mendoza↗

Fitting strategies of accretion column models and application to the broadband spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT dataset of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies.

X-rays: binaries↗

NuSTAR Observation of GROJ1744−28 at Low Mass Accretion Rate

Context.Neutron stars in low-mass X-ray binaries (LMXBs) are important systems to study the physics of accretion onto compact ob-jects. The system GRO J1744−28 is particularly interesting, as it usually shows clear pulsations as well as X-ray bursts. Additionally,there are claims for a magnetic field of 5×1011G through the detection of a Cyclotron Resonant Scattering Feature (CRSF).Aims.We present the spectral analysis of GRO J1744−28 using∼29 ks of NuSTAR data taken in 2017 February at a low luminosityof 3.2×1036erg s−1(3–50 keV). Our goal is to study the variability of the source spectrum with pulse phase and to search for the claimed CRSF. Methods.The continuum spectrum is modeled with an absorbed power-law with exponential cut-off, and an additional iron line component. We find no obvious indications for a CRSF and therefore perform a detailed cyclotron line search using statistical methods.We perform this search on pulse phase-averaged as well as phase-resolved spectra. Results.GRO J1744−28 was observed at one of the lowest luminosities studied so far. The previously detected Type II X-ray bursts re absent. Clear pulsations at a period of 2.141124(9) Hz are detected. The pulse profile shows an indication of a secondary peak,which was not seen at higher flux. The 4σupper limit for the strength of a CRSF in the 3–20 keV band is 0.07 keV, lower than the strength of the line found at higher luminosity. Conclusions.The detection of pulsations shows that the source did not enter the “propeller” regime, even though the source flux of4.15×10−10erg cm−2s−1was almost one order of magnitude below the threshold for the propeller regime claimed in previous studies on this source. The transition into the propeller regime in GRO J1744−28 must therefore be below a luminosity of 3.2×1036erg s−1(3–50 keV), which implies a surface magnetic field.2.9×1011G and mass accretion rate.1.7×1016g s−1. A change of the CRSF depth as function of luminosity is not unexpected and has been observed in other sources. This result possibly implies a change in emission geometry as function of mass accretion rate to reduce the depth of the line below our detection limit.

Ole König↗

Dust and gas absorption in the High Mass X-ray Binary IGR J16318−4848

Context. With an absorption column density on the order of 10(exp 24)per sq. cm, IGR J16318-4848 is one of the most extreme cases of a highly obscured High Mass X-ray Binary. Besides the overall continuum absorption, the source spectrum exhibits a strong iron and nickel fluorescence line complex at 6.4 keV. Previous empirical modeling of these features and comparison with radiative transfer simulations raised questions about the structure and covering fraction of the absorber and the profile of the fluorescence lines. Aims. We aim at a self-consistent description of the continuum absorption, the absorption edges and the fluorescent lines to constrain properties of the absorbing material, such as ionization structure and geometry. We further investigate the effects of dust absorption on the observed spectra and the possibility of fluorescent emission from dust grains. Methods. We use XMM-Newton and NuSTAR spectra to first constrain empirically the incident continuum and fluorescence lines. Next, we use XSTAR to construct a customized photoionization model where we vary the ionization parameter, column density, and covering fraction. In a third step, we model the absorption and fluorescence in a dusty olivine absorber and employ both, a simple analytical model for the fluorescent line emission and Monte Carlo radiative transfer, spectral shapes and line fluxes that are very close to the data are generated. Results. Our empirical spectral modeling is in agreement with previous works. Our second model, the single gas absorber does not describe the observational data. In particular, irrespective of the ionization state or column density of the absorber, a much higher covering fraction than previously estimated is needed to produce the strong fluorescence lines and the large continuum absorption. A dusty, spherical absorber (modeled as consisting of olivine dust, although the nature of dust cannot be constrained) is well able to produce the observed continuum absorption and edges. Conclusions. A dense, dusty absorber in the direct vicinity of the source consisting of dust offers a consistent description of both the strong continuum absorption and the strong emission features in the X-ray spectrum of IGR J16318􀀀4848. In particular, for low optical depth of individual grains the dust will contribute significantly to the fluorescent emission, which is the case for typical densities and grain size distribution models.

Ralf Ballhausen↗

High Resolution Photoexcitation Measurements Exacerbate the Long-Standing Fe XVII Oscillator Strength Problem

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 rations 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.

Electron-ion collisions↗

The Giant Outburst of 4U 0115+634 in 2011 with Suzaku and RXTE

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 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 keV. Apart from these features, we find for the first time evidence for a He-like Fe XXV emission line at keV and weak H-like Fe XXVI emission close to keV.

Matthias Bissinger né Kühnel↗