Modeling the light curves of ultraluminous X-ray sources as precession

We present a freely available xspec model for the modulations seen in the long-term light curves of multiple ultraluminous X-ray sources (ULXs). By incorporating the physics of multiple electron scatterings (ray traced with a Monte-Carlo routine), we go beyond analytical predictions and show that the geometrical beaming of radiation in the conical outflow can be more than a factor of 100 for opening angles smaller than $10^\circ$. We apply our new model to the long-term, well sampled Swift light curve of the recently confirmed ULX pulsar NGC 5907 X-1 with an established period of 78 days. Our results suggest that geometrical beaming together with a slight precession of the conical wind can describe the light curve with a consistent set of parameters for the wind. The small opening angle of roughly $10\mathrm{-}13^\circ$ implies a highly super-critical flow and boosting factors at the order of $\mathcal{B}=60\mathrm{-}90$ that would yield a fairly low surface magnetic field strength of $2\times 10^{10}\,$Gauss.Comment: accepted by MNRAS for publication (7 pages, 6 figures

X-ray reflected spectra from accretion disk models. III. A complete grid of ionized reflection calculations

We present a new and complete library of synthetic spectra for modeling the component of emission that is reflected from an illuminated accretion disk. The spectra were computed using an updated version of our code XILLVER that incorporates new routines and a richer atomic data base. We offer in the form of a table model an extensive grid of reflection models that cover a wide range of parameters. Each individual model is characterized by the photon index \Gamma of the illuminating radiation, the ionization parameter \xi at the surface of the disk (i.e., the ratio of the X-ray flux to the gas density), and the iron abundance A_{Fe} relative to the solar value. The ranges of the parameters covered are: 1.2 \leq \Gamma \leq 3.4, 1 \leq \xi \leq 10^4, and 0.5 \leq A_{Fe} \leq 10. These ranges capture the physical conditions typically inferred from observations of active galactic nuclei, and also stellar-mass black holes in the hard state. This library is intended for use when the thermal disk flux is faint compared to the incident power-law flux. The models are expected to provide an accurate description of the Fe K emission line, which is the crucial spectral feature used to measure black hole spin. A total of 720 reflection spectra are provided in a single FITS file{\url{http://hea-www.cfa.harvard.edu/~javier/xillver/}} suitable for the analysis of X-ray observations via the atable model in XSPEC. Detailed comparisons with previous reflection models illustrate the improvements incorporated in this version of XILLVER.Comment: 70 pages, 21 figures, submitted to Ap

A Suzaku, NuSTAR, and XMM-Newton view on variable absorption and relativistic reflection in NGC 4151

We disentangle X-ray disk reflection from complex line-of-sight absorption in the nearby Seyfert NGC 4151, using a suite of Suzaku, NuSTAR, and XMM-Newton observations. Extending upon earlier published work, we pursue a physically motivated model using the latest angle-resolved version of the lamp-post geometry reflection model relxillCp_lp together with a Comptonization continuum. We use the long-look simultaneous Suzaku/NuSTAR observation to develop a baseline model wherein we model reflected emission as a combination of lamp-post components at the heights of 1.2 and 15.0 gravitational radii. We argue for a vertically extended corona as opposed to two compact and distinct primary sources. We find two neutral absorbers (one full-covering and one partial-covering), an ionized absorber ($\log \xi = 2.8$), and a highly-ionized ultra-fast outflow, which have all been reported previously. All analyzed spectra are well described by this baseline model. The bulk of the spectral variability between 1 keV and 6 keV can be accounted for by changes in the column density of both neutral absorbers, which appear to be degenerate and inversely correlated with the variable hard continuum component flux. We track variability in absorption on both short (2 d) and long ($\sim$1 yr) timescales; the observed evolution is either consistent with changes in the absorber structure (clumpy absorber at distances ranging from the broad line region (BLR) to the inner torus or a dusty radiatively driven wind) or a geometrically stable neutral absorber that becomes increasingly ionized at a rising flux level. The soft X-rays below 1 keV are dominated by photoionized emission from extended gas that may act as a warm mirror for the nuclear radiation.Comment: 21 pages, 19 figures, 8 tables, accepted for publication by A&

Normalizing a relativistic model of X-ray reflection: Definition of the reflection fraction and its implementation in relxill

Aims. The only relativistic reflection model that implements a parameter relating the intensity incident on an accretion disk to the observed intensity is relxill. The parameter used in earlier versions of this model, referred to as the reflection strength, is unsatisfactory; it has been superseded by a parameter that provides insight into the accretion geometry, namely the reflection fraction. The reflection fraction is defined as the ratio of the coronal intensity illuminating the disk to the coronal intensity that reaches the observer. Methods. The relxill model combines a general relativistic ray-tracing code and a photoionization code to compute the component of radiation reflected from an accretion that is illuminated by an external source. The reflection fraction is a particularly important parameter for relativistic models with well-defined geometry, such as the lamp post model, which is a focus of this paper. Results. Relativistic spectra are compared for three inclinations and for four values of the key parameter of the lamp post model, namely the height above the black hole of the illuminating, on-axis point source. In all cases, the strongest reflection is produced for low source heights and high spin. A low-spin black hole is shown to be incapable of producing enhanced relativistic reflection. Results for the relxill model are compared to those obtained with other models and a Monte Carlo simulation. Conclusions. Fitting data by using the relxill model and the recently implemented reflection fraction, the geometry of a system can be constrained. The reflection fraction is independent of system parameters such as inclination and black hole spin. The reflection-fraction parameter was implemented with the name refl_frac in all flavours of the relxill model, and the non-relativistic reflection model xillver, in v0.4a (18 January 2016)

Optical-NIR spectroscopy of the puzzling gamma-ray source 3FGL 1603.9-4903/PMN J1603-4904 with X-shooter

The Fermi/LAT instrument has detected about two thousands Extragalactic High Energy (E > 100 MeV) gamma-ray sources. One of the brightest is 3FGL 1603.9-4903, associated to the radio source PMN J1603-4904. Its nature is not yet clear, it could be either a very peculiar BL Lac or a CSO (Compact Symmetric Object) radio source, considered as the early stage of a radio galaxy. The latter, if confirmed, would be the first detection in gamma-rays for this class of objects. Recently a redshift z=0.18 +/- 0.01 has been claimed on the basis of the detection of a single X-ray line at 5.44 +/- 0.05 keV interpreted as a 6.4 keV (rest frame) fluorescent line. We aim to investigate the nature of 3FGL 1603.9-4903/PMN J1603-4904 using optical to NIR spectroscopy. We observed PMN J1603-4904 with the UV-NIR VLT/X-shooter spectrograph for two hours. We extracted spectra in the VIS and NIR range that we calibrated in flux and corrected for telluric absorption and we systematically searched for absorption and emission features. The source was detected starting from ~6300 Ang down to 24000 Ang with an intensity comparable to the one of its 2MASS counterpart and a mostly featureless spectrum. The continuum lacks absorption features and thus is non-stellar in origin and likely non-thermal. On top of this spectrum we detected three emission lines that we interpret as the Halpha-[NII] complex, the [SII] 6716,6731 doublet and the [SIII] 9530 line, obtaining a redshift estimate of z= 0.2321 +/- 0.0004. The equivalent width of the Halpha-[NII] complex implies that PMN J1603-4904 does not follow the observational definition of BL Lac, the line ratios suggest that a LINER/Seyfert nucleus is powering the emission. This new redshift measurement implies that the X-ray line previously detected should be interpreted as a 6.7 keV line which is very peculiar.Comment: Published in Astronomy and Astrophysic

A redshifted Fe Kα\alpha line from the unusual gamma-ray source PMN J1603-4904

Multiwavelength observations have revealed the highly unusual properties of the gamma-ray source PMN J1603-4904, which are difficult to reconcile with any other well established gamma-ray source class. The object is either a very atypical blazar or compact jet source seen at a larger angle to the line of sight. In order to determine the physical origin of the high-energy emission processes in PMN J1603-4904, we study the X-ray spectrum in detail. We performed quasi-simultaneous X-ray observations with XMM-Newton and Suzaku in 2013 September, resulting in the first high signal-to-noise X-ray spectrum of this source. The 2-10 keV X-ray spectrum can be well described by an absorbed power law with an emission line at 5.44$\pm$0.05 keV (observed frame). Interpreting this feature as a K{\alpha} line from neutral iron, we determine the redshift of PMN J1603-4904 to be z=0.18$\pm$0.01, corresponding to a luminosity distance of 872$\pm$54 Mpc. The detection of a redshifted X-ray emission line further challenges the original BL Lac classification of PMN J1603-4904. This result suggests that the source is observed at a larger angle to the line of sight than expected for blazars, and thus the source would add to the elusive class of gamma-ray loud misaligned-jet objects, possibly a {\gamma}-ray bright young radio galaxy.Comment: 5 pages, 1 figure, A&A accepte

Investigating source confusion in PMN J1603−-4904

PMN J1603$-$4904 is a likely member of the rare class of $\gamma$-ray emitting young radio galaxies. Only one other source, PKS 1718$-$649, has been confirmed so far. These objects, which may transition into larger radio galaxies, are a stepping stone to understanding AGN evolution. It is not completely clear how these young galaxies, seen edge-on, can produce high-energy $\gamma$-rays. PMN J1603$-$4904 has been detected by TANAMI Very Long Baseline Interferometry (VLBI) observations and has been followed-up with multiwavelength observations. A Fermi/LAT $\gamma$-ray source has been associated with it in the LAT catalogs. We have obtained Chandra observations of the source in order to consider the possibility of source confusion, due to the relatively large positional uncertainty of Fermi/LAT. The goal was to investigate the possibility of other X-ray bright sources in the vicinity of PMN J1603$-$4904 that could be counterparts to the $\gamma$-ray emission. With Chandra/ACIS, we find no other sources in the uncertainty ellipse of Fermi/LAT data, which includes an improved localization analysis of 8 years of data. We further study the X-ray fluxes and spectra. We conclude that PMN J1603$-$4904 is indeed the second confirmed $\gamma$-ray bright young radio galaxy.Comment: 4 pages, 3 figures, accepted for publication in A&

Formation of phase lags at the cyclotron energies in the pulse profiles of magnetized, accreting neutron stars

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.Comment: 4 pages, 4 figures; updated reference lis