Multiple cyclotron line-forming regions in GX 301-2

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.Comment: 18 pages, 10 figures, accepted for publication in A&

Cyclotron resonant scattering feature simulations. I. Thermally averaged cyclotron scattering cross sections, mean free photon-path tables, and electron momentum sampling

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.Comment: A&A, in pres

Multiple cyclotron line-forming regions in GX 301−2

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 have been seen in GX 301−2. We find that the CRSFs are very likely independently formed, as their energies are not harmonically related and, if the observed feature were due to 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 a function of phase, between 30 and 38 keV. To explain this variation, we apply a simple model of the accretion column, taking into account the altitude of the line-forming region, the velocity of the in-falling material, and the resulting relativistic effects. We find that in this model the observed energy variation can be explained as being simply due to a variation of the projected velocity and beaming factor of the line-forming region towards us

Cyclotron resonant scattering feature simulations

Context. 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 1012 G) close to the surface of an accreting X-ray pulsar. Due to the complex scattering cross-sections, 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. Aims. 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. Methods. 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 monoenergetic photon injections, producing Green’s functions, which can be used later to generate spectra for arbitrary continua. Results. 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-

A Precessing Be Disc as a Possible Model for Occultation Events in GX 304-1

We report on the RXTE detection of a sudden increase in the absorption column density, NH, during the 2011 May outburst of GX 3041. The NH increased up to 16 10(exp 22) atoms cm(exp 2), which is a factor of 34 larger than what is usually measured during the outbursts of GX 3041 as covered by RXTE. Additionally, an increase in the variability of the hardness ratio as calculated from the energy resolved RXTE-Proportional Counter Array 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 disc around the Be-type companion, we are able to qualitatively explain the NH evolution over time. We are able to constrain the Be disc density to be of the order of 10(exp 11)g cm(exp 3). Our model strengthens the idea of inclined Be discs as origin of double-peaked outbursts as the derived geometry allows accretion twice per orbit under certain conditions