Showing 1–50 of 80 results for author: Middleton, M J

The High Energy X-ray Probe (HEX-P): Studying Extreme Accretion with Ultraluminous X-ray Sources

Authors: Matteo Bachetti, Matthew J. Middleton, Ciro Pinto, Andrés Gúrpide, Dominic J. Walton, Murray Brightman, Bret Lehmer, Timothy P. Roberts, Georgios Vasilopoulos, Jason Alford, Roberta Amato, Elena Ambrosi, Lixin Dai, Hannah P. Earnshaw, Hamza El Byad, Javier A. García, Gian Luca Israel, Amruta Jaodand, Kristin Madsen, Chandreyee Maitra, Shifra Mandel, Kaya Mori, Fabio Pintore, Ken Ohsuga, Maura Pilia , et al. (3 additional authors not shown)

Abstract: Ultraluminous X-ray sources (ULXs) represent an extreme class of accreting compact objects: from the identification of some of the accretors as neutron stars to the detection of powerful winds travelling at 0.1-0.2 c, the increasing evidence points towards ULXs harbouring stellar-mass compact objects undergoing highly super-Eddington accretion. Measuring their intrinsic properties, such as the acc… ▽ More Ultraluminous X-ray sources (ULXs) represent an extreme class of accreting compact objects: from the identification of some of the accretors as neutron stars to the detection of powerful winds travelling at 0.1-0.2 c, the increasing evidence points towards ULXs harbouring stellar-mass compact objects undergoing highly super-Eddington accretion. Measuring their intrinsic properties, such as the accretion rate onto the compact object, the outflow rate, the masses of accretor/companion -- hence their progenitors, lifetimes, and future evolution -- is challenging due to ULXs being mostly extragalactic and in crowded fields. Yet ULXs represent our best opportunity to understand super-Eddington accretion physics and the paths through binary evolution to eventual double compact object binaries and gravitational wave sources. Through a combination of end-to-end and single-source simulations, we investigate the ability of HEX-P to study ULXs in the context of their host galaxies and compare it to XMM-Newton and NuSTAR, the current instruments with the most similar capabilities. HEX-P's higher sensitivity, which is driven by its narrow point-spread function and low background, allows it to detect pulsations and broad spectral features from ULXs better than XMM-Newton and NuSTAR. We describe the value of HEX-P in understanding ULXs and their associated key physics, through a combination of broadband sensitivity, timing resolution, and angular resolution, which make the mission ideal for pulsation detection and low-background, broadband spectral studies. △ Less

Submitted 10 November, 2023; v1 submitted 8 November, 2023; originally announced November 2023.

Comments: 24 pages, 9 figures. Accepted for publication in Frontiers of Astronomy and Space Science. Minor corrections included

Long-Term X-Ray/UV Variability in ULXs

Authors: Norman Khan, Matthew. J. Middleton

Abstract: The focus of NASA's Swift telescope has been transients and target-of-opportunity observing, resulting in many observations of ultraluminous X-ray sources (ULXs) over the last ~20 years. For the vast majority of these observations, simultaneous data has been obtained using both the X-ray telescope (XRT) and the ultraviolet and optical telescope (UVOT), providing a unique opportunity to study coupl… ▽ More The focus of NASA's Swift telescope has been transients and target-of-opportunity observing, resulting in many observations of ultraluminous X-ray sources (ULXs) over the last ~20 years. For the vast majority of these observations, simultaneous data has been obtained using both the X-ray telescope (XRT) and the ultraviolet and optical telescope (UVOT), providing a unique opportunity to study coupled variability between these bands. Using a sample of ~40 ULXs with numerous repeat observations, we extract stacked images to characterise the spatial extent of the UV-Optical emission and extract long-term light curves to search for first-order linear correlations between the UV and X-ray emission. We find that a small subset may show weakly correlated joint variability, while other sources appear to display non-linear relationships between the bands. We discuss these observations in the context of several theoretical models: precession, irradiation of the outer accretion disc and irradiation of the companion star. We conclude that more complicated analysis or higher quality data may be required to accurately constrain the nature of the joint X-ray and UV/optical emission in these sources. △ Less

Submitted 15 August, 2023; originally announced August 2023.

Comments: 14 pages, 13 Figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 524, Issue 3, pp.4302-4314, September 2023

Drop in the hard pulsed fraction and a candidate cyclotron line in IGR J16320-4751 seen by NuSTAR

Authors: Arash Bodaghee, Alan J. -L. Chiu, John A. Tomsick, Varun Bhalerao, Eugenio Bottacini, Maica Clavel, Cody Cox, Felix Fürst, Matthew J. Middleton, Farid Rahoui, Jerome Rodriguez, Pat Romano, Joern Wilms

Abstract: We report on a timing and spectral analysis of a 50-ks NuSTAR observation of IGR J16320-4751 (= AX J1631.9-4752); a high-mass X-ray binary hosting a slowly-rotating neutron star. In this observation from 2015, the spin period was 1,308.8+/-0.4 s giving a period derivative dP/dt ~ 2E-8 s s-1 when compared with the period measured in 2004. In addition, the pulsed fraction decreased as a function of… ▽ More We report on a timing and spectral analysis of a 50-ks NuSTAR observation of IGR J16320-4751 (= AX J1631.9-4752); a high-mass X-ray binary hosting a slowly-rotating neutron star. In this observation from 2015, the spin period was 1,308.8+/-0.4 s giving a period derivative dP/dt ~ 2E-8 s s-1 when compared with the period measured in 2004. In addition, the pulsed fraction decreased as a function of energy, as opposed to the constant trend that was seen previously. This suggests a change in the accretion geometry of the system during the intervening 11 years. The phase-averaged spectra were fit with the typical model for accreting pulsars: a power law with an exponential cutoff. This left positive residuals at 6.4 keV attributable to the known iron K-alpha line, as well as negative residuals around 14 keV from a candidate cyclotron line detected at a significance of 5-sigma. We found no significant differences in the spectral parameters across the spin period, other than the expected changes in flux and component normalizations. A flare lasting around 5 ks was captured during the first half of the observation where the X-ray emission hardened and the local column density decreased. Finally, the binary orbital period was refined to 8.9912+/-0.0078 d thanks to Swift/BAT monitoring data from 2005-2022. △ Less

Submitted 11 May, 2023; originally announced May 2023.

Comments: 17 pages, 11 figures, Referee-revised version accepted for publication in the Astrophysical Journal

A new sample of transient ultraluminous X-ray sources serendipitously discovered by Swift/XRT

Authors: Murray Brightman, Jean-Marie Hameury, Jean-Pierre Lasota, Ranieri D. Baldi, Gabriele Bruni, Jenna M. Cann, Hannah Earnshaw, Felix Fürst, Marianne Heida, Amruta Jaodand, Margaret Lazzarini, Matthew J. Middleton, Dominic J. Walton, Kimberly A. Weaver

Abstract: Ultraluminous X-ray sources (ULXs) are our best laboratories for studying extreme super-Eddington accretion. Most studies of these objects are of relatively persistent sources, however there is growing evidence to suggest a large fraction of these sources are transient. Here we present a sample of five newly reported transient ULXs in the galaxies NGC 4945, NGC 7793 and M81 serendipitously discove… ▽ More Ultraluminous X-ray sources (ULXs) are our best laboratories for studying extreme super-Eddington accretion. Most studies of these objects are of relatively persistent sources, however there is growing evidence to suggest a large fraction of these sources are transient. Here we present a sample of five newly reported transient ULXs in the galaxies NGC 4945, NGC 7793 and M81 serendipitously discovered in Swift/XRT observations. Swift monitoring of these sources have provided well sampled lightcurves, allowing for us to model the lightcurves with the disk instability model of Hameury & Lasota (2020) which implies durations of 60-400 days and that the mass accretion rate through the disk is close to or greater than the Eddington rate. Of the three source regions with prior HST imaging, color magnitude diagrams of the potential stellar counterparts show varying ages of the possible stellar counterparts. Our estimation of the rates of these sources in these three galaxies is 0.4-1.3 year$^{-1}$. We find that while persistent ULXs dominate the high end of galaxy luminosity functions, the number of systems that produce ULX luminosities are likely dominated by transient sources. △ Less

Submitted 2 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ

Spectral calculations of 3D RMHD simulations of super-Eddington accretion onto a stellar-mass black hole

Authors: Brianna S. Mills, Shane W. Davis, Yan-Fei Jiang, Matthew J. Middleton

Abstract: We use the Athena++ Monte Carlo (MC) radiation transfer module to post-process simulation snapshots from non-relativistic Athena++ radiation magnetohydrodynamic (RMHD) simulations. These simulations were run using a gray (frequency integrated) approach but were also restarted and ran with a multi-group approach that accounts for Compton scattering with a Kompaneets operator. These simulations prod… ▽ More We use the Athena++ Monte Carlo (MC) radiation transfer module to post-process simulation snapshots from non-relativistic Athena++ radiation magnetohydrodynamic (RMHD) simulations. These simulations were run using a gray (frequency integrated) approach but were also restarted and ran with a multi-group approach that accounts for Compton scattering with a Kompaneets operator. These simulations produced moderately super-Eddington accretion rates onto a 6.62 $M_\odot$ black hole. Since we only achieve inflow equilibrium out to 20-25 gravitational radii, we focus on the hard X-ray emission. We provide a comparison between the MC and RMHD simulations showing that the treatment of Compton scattering in the gray RMHD simulations underestimates the gas temperature in the regions above and below the accretion disk. In contrast, the restarted multi-group snapshots provides a treatment for the radiation field that is more consistent with the MC calculations, and result in post-processed spectra with harder X-ray emission compared to their gray snapshot counterparts. We characterize these MC post-processed spectra using commonly employed phenomenological models used for spectral fitting. We also attempt to fit our MC spectra directly to observations of the ultraluminous X-ray source (ULX) NGC 1313 X-1, finding best fit values that are competitive to phenomenological model fits, indicating that first principle models of super-Eddington accretion may adequately explain the observed hard X-ray spectra in some ULX sources. △ Less

Submitted 16 April, 2023; originally announced April 2023.

Comments: Submitted to ApJ; 20 pages, 15 figures,

Probing the nature of the low state in the extreme ultraluminous X-ray pulsar NGC 5907 ULX1

Authors: F. Fuerst, D. J. Walton, G. L. Israel, M. Bachetti, D. Barret, M. Brightman, H. P. Earnshaw, A. Fabian, M. Heida, M. Imbrogno, M. J. Middleton, C. Pinto, R. Salvaterra, T. P. Roberts, G. A. Rodríguez Castillo, N. Webb

Abstract: NGC 5907 ULX1 is the most luminous ultra-luminous X-ray pulsar (ULXP) known to date, reaching luminosities in excess of 1e41 erg/s. The pulsar is known for its fast spin-up during the on-state. Here, we present a long-term monitoring of the X-ray flux and the pulse period between 2003-2022. We find that the source was in an off- or low-state between mid-2017 to mid-2020. During this state, our pul… ▽ More NGC 5907 ULX1 is the most luminous ultra-luminous X-ray pulsar (ULXP) known to date, reaching luminosities in excess of 1e41 erg/s. The pulsar is known for its fast spin-up during the on-state. Here, we present a long-term monitoring of the X-ray flux and the pulse period between 2003-2022. We find that the source was in an off- or low-state between mid-2017 to mid-2020. During this state, our pulse period monitoring shows that the source had spun down considerably. We interpret this spin-down as likely being due to the propeller effect, whereby accretion onto the neutron star surface is inhibited. Using state-of-the-art accretion and torque models, we use the spin-up and spin-down episodes to constrain the magnetic field. For the spin-up episode, we find solutions for magnetic field strengths of either around 1e12G or 1e13G, however, the strong spin-down during the off-state seems only to be consistent with a very high magnetic field, namely, >1e13G. This is the first time a strong spin-down is seen during a low flux state in a ULXP. Based on the assumption that the source entered the propeller regime, this gives us the best estimate so far for the magnetic field of NGC 5907 ULX1. △ Less

Submitted 7 February, 2023; originally announced February 2023.

Comments: 10 pages, 1 figure, accepted for publication in A&A

Journal ref: A&A 672, A140 (2023)

Global Three-Dimensional Radiation Magnetohydrodynamic Simulations of Accretion onto a Stellar Mass Black Hole at Sub- and Near-critical Accretion Rates

Authors: Jiahui Huang, Yan-Fei Jiang, Hua Feng, Shane W. Davis, James M. Stone, Matthew J. Middleton

Abstract: We present global 3D radiation magnetohydrodynamical simulations of accretion onto a 6.62 solar mass black hole with quasi-steady state accretion rates reaching 0.016 to 0.9 times the critical accretion rate, which is defined as the accretion rate to power the Eddington luminosity assuming a 10% radiative efficiency, in different runs. The simulations show no sign of thermal instability over hundr… ▽ More We present global 3D radiation magnetohydrodynamical simulations of accretion onto a 6.62 solar mass black hole with quasi-steady state accretion rates reaching 0.016 to 0.9 times the critical accretion rate, which is defined as the accretion rate to power the Eddington luminosity assuming a 10% radiative efficiency, in different runs. The simulations show no sign of thermal instability over hundreds of thermal timescales at 10 $r_{\rm g}$. The energy dissipation happens close to the mid-plane in the near-critical runs and near the disk surface in the low accretion rate run. The total radiative luminosity inside $\sim$20 $r_{\rm g}$ is about 1% to 30% the Eddington limit, with a radiative efficiency of about 6% and 3%, respectively, in the sub- and near-critical accretion regimes. In both cases, self-consistent turbulence generated by the magnetorotational instability (MRI) leads to angular momentum transfer, and the disk is supported by magnetic pressure. Outflows from the central low-density funnel with a terminal velocity of $\sim$0.1$c$ are seen only in the near-critical runs. We conclude that these magnetic pressure dominated disks are thermally stable and thicker than the $α$ disk, and the effective temperature profiles are much flatter than that in the $α$ disks. The magnetic pressure of these disks are comparable within an order of magnitude with the previous analytical magnetic pressure dominated disk model. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 17 pages, 13 figures, 3 tables, accepted for publication in ApJ

The ultraluminous X-ray source M81 X-6: a weakly magnetised neutron star with a precessing accretion disc?

Authors: Roberta Amato, Andres Gúrpide, Natalie A. Webb, Olivier Godet, Matthew J. Middleton

Abstract: We investigate the nature of the ULX M81 X-6, which has been suggested to harbour a neutron star (NS), by studying its long-term X-ray spectral and temporal evolution, using the rich set of available archival data from XMM-Newton, Chandra, NuSTAR, and Swift/XRT. We tracked the evolution of the source on the hardness-intensity diagram and find that the source oscillates between two main states: one… ▽ More We investigate the nature of the ULX M81 X-6, which has been suggested to harbour a neutron star (NS), by studying its long-term X-ray spectral and temporal evolution, using the rich set of available archival data from XMM-Newton, Chandra, NuSTAR, and Swift/XRT. We tracked the evolution of the source on the hardness-intensity diagram and find that the source oscillates between two main states: one characterised by a hard and luminous spectrum and the other at low hardness and luminosity. The properties of the soft component remain constant between the two states, suggesting that changes in the mass-transfer rate are not driving the spectral transitions. Instead, the bi-modal behaviour of the source and the known super-orbital period would point to the precession of the accretion disc. Here, we tested two theoretical models: (1) Lense-Thirring precession, which can explain the super-orbital period if the NS has a magnetic field $B$ $\lesssim10^{10}$ G, supporting the idea of M81 X-6 as a weakly magnetised NS, and (2) precession due to the torque of the NS magnetic field, which leads to $B \gtrsim$ 10$^{11}$ G. However, the latter scenario, assuming M81 X-6 shares similar properties with other NS-ULXs, is disfavoured because it would require magnetic field strengths ($B>10^{15}$ G) much higher than those known for other pulsating ULXs. We further show that the contribution from the hard component attributed to the putative accretion column sits just below the typical values found in pulsating ULXs, which, together with the low value of the pulsed fraction ($\leq10$\%) found for one XMM-Newton/pn observation, could explain the source's lack of pulsations. The spectral properties and variability of M81 X-6 can be accounted for if the accretor is a NS with a low magnetic field. Under the hypothesis of Lense-Thirring precession, we predict a spin period of the NS of a few seconds. △ Less

Submitted 14 December, 2022; originally announced December 2022.

Comments: 14 pages, 10 figures. Accepted for publication on Astronomy & Astrophysics

Journal ref: A&A 669, A130 (2023)

The variability behavior of NGC 925 ULX-3

Authors: Hannah P. Earnshaw, Murray Brightman, Fiona A. Harrison, Marianne Heida, Amruta Jaodand, Matthew J. Middleton, Timothy P. Roberts, Dominic J. Walton

Abstract: We report the results of a 2019-2021 monitoring campaign with Swift and associated target-of-opportunity observations with XMM-Newton and NuSTAR, examining the spectral and timing behavior of the highly variable ultraluminous X-ray source (ULX) NGC 925 ULX-3. We find that the source exhibits a 127-128 day periodicity, with fluxes typically ranging from 1e-13 to 8e-13 ergs/s/cm2. We do not find str… ▽ More We report the results of a 2019-2021 monitoring campaign with Swift and associated target-of-opportunity observations with XMM-Newton and NuSTAR, examining the spectral and timing behavior of the highly variable ultraluminous X-ray source (ULX) NGC 925 ULX-3. We find that the source exhibits a 127-128 day periodicity, with fluxes typically ranging from 1e-13 to 8e-13 ergs/s/cm2. We do not find strong evidence for a change in period over the time that NGC 925 ULX-3 has been observed, although the source may have been in a much lower flux state when first observed with Chandra in 2005. We do not detect pulsations, and we place an upper limit on the pulsed fraction of ~40% in the XMM-Newton band, consistent with some previous pulsation detections at low energies in other ULXs. The source exhibits a typical ULX spectrum that turns over in the NuSTAR band and can be fitted using two thermal components. These components have a high temperature ratio that may indicate the lack of extreme inner disk truncation by a magnetar-level magnetic field. We examine the implications for a number of different models for superorbital periods in ULXs, finding that a neutron star with a magnetic field of ~10^12 G may be plausible for this source. The future detection of pulsations from this source would allow for the further testing and constraining of such models. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: 13 pages, 6 figures, 2 tables. Accepted for publication in ApJ

Constraining Energy-dependent emissivity profiles of AGN inflows

Authors: D. I. Ashton, M. J. Middleton

Abstract: The emissivity of the accretion flow is a key parameter affecting the shape of both the energy and variability power spectrum of AGN. We explore the energy-dependence of the power spectrum for five AGN, across the XMM-Newton bandpass, and across the 0.01-1 mHz frequency range, finding a ubiquitous flattening of the power spectrum towards higher energies. We develop a framework to explore this beha… ▽ More The emissivity of the accretion flow is a key parameter affecting the shape of both the energy and variability power spectrum of AGN. We explore the energy-dependence of the power spectrum for five AGN, across the XMM-Newton bandpass, and across the 0.01-1 mHz frequency range, finding a ubiquitous flattening of the power spectrum towards higher energies. We develop a framework to explore this behaviour and thereby extract the energy dependence of the emissivity assuming a simple disc-like geometry for the inflow. We find that the emissivity ranges from R^-2 at energies around the soft excess and increases to R^-4 or steeper above ~4-6 keV. We describe the changing emissivity index with a linear function in energy, finding the best-fitting slopes to vary between AGN. We attempt to correlate the slope of the linear function against key AGN parameters but, as yet, the sample size is too small to confirm hints of a correlation with Eddington ratio. △ Less

Submitted 21 April, 2022; originally announced April 2022.

Comments: 12 pages, 6 figures, 2 tables. Accepted for publication in MNRAS

The X-ray Disk/Wind Degeneracy in AGN

Authors: M. L. Parker, G. A. Matzeu, J. H. Matthews, M. J. Middleton, T. Dauser, J. Jiang, A. M. Joyce

Abstract: Relativistic Fe K emission lines from accretion disks and from disk winds encode key information about black holes, and their accretion and feedback mechanisms. We show that these two processes can in principle produce indistinguishable line profiles, such that they cannot be disentangled spectrally. We argue that it is likely that in many cases both processes contribute to the net line profile, a… ▽ More Relativistic Fe K emission lines from accretion disks and from disk winds encode key information about black holes, and their accretion and feedback mechanisms. We show that these two processes can in principle produce indistinguishable line profiles, such that they cannot be disentangled spectrally. We argue that it is likely that in many cases both processes contribute to the net line profile, and their relative contributions cannot be constrained purely by Fe K spectroscopy. In almost all studies of Fe K emission to date, a single process (either disk reflection or wind Compton scattering) is assumed to dominate the total line profile. We demonstrate that fitting a single process emission model (pure reflection or pure wind) to a hybrid line profile results in large systematic biases in the estimates of key parameters, such as mass outflow rate and spin. We discuss various strategies to mitigate this effect, such as including high energy data covering the Compton hump, and the implications for future X-ray missions. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: 22 pages, 17 figures, accepted for publication in MNRAS

An 8.56 keV absorption line in the hyperluminous X-ray source in NGC 4045: ultra-fast outflow or cyclotron line?

Authors: Murray Brightman, Peter Kosec, Felix Fürst, Hannah Earnshaw, Marianne Heida, Matthew J Middleton, Daniel Stern, Dominic J Walton

Abstract: We report on the discovery of an absorption line at $E=8.56^{+0.05}_{-0.11}$ keV detected with a significance of $>3.3σ$ in the NuSTAR and XMM-Newton spectra of a newly discovered hyperluminous X-ray source (HLX, $L_{\rm X}>10^{41}$ ergs$^{-1}$) in the galaxy NGC 4045 at a distance of 32 Mpc. The source was first discovered serendipitously in a Swift/XRT observation of the galaxy, and Swift monito… ▽ More We report on the discovery of an absorption line at $E=8.56^{+0.05}_{-0.11}$ keV detected with a significance of $>3.3σ$ in the NuSTAR and XMM-Newton spectra of a newly discovered hyperluminous X-ray source (HLX, $L_{\rm X}>10^{41}$ ergs$^{-1}$) in the galaxy NGC 4045 at a distance of 32 Mpc. The source was first discovered serendipitously in a Swift/XRT observation of the galaxy, and Swift monitoring reveals a highly variable source changing by over an order of magnitude from maximum to minimum. The origin of the absorption line appears likely to be by highly ionized iron with a blue shift of 0.19$c$, indicating an ultrafast outflow (UFO). However, the large equivalent width of the line (EW$=-0.22^{+0.08}_{-0.09}$ keV) paired with the lack of other absorption lines detected are difficult to reconcile with models. An alternative explanation is that the line is due to a cyclotron resonance scattering feature (CRSF), produced by the interaction of X-ray photons with the powerful magnetic field of a neutron star. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: Accepted for publication in ApJ

Detection of periodic signals in AGN red noise light curves: empirical tests on the Auto-Correlation Function and Phase Dispersion Minimization

Authors: S. Krishnan, A. G. Markowitz, A. Schwarzenberg-Czerny, M. J. Middleton

Abstract: Active galactic nucleus (AGN) emission is dominated by stochastic, aperiodic variability which overwhelms any periodic/quasi-periodic signal (QPO) if one is present. The Auto Correlation Function (ACF) and Phase Dispersion Minimization (PDM) techniques have been used previously to claim detections of QPOs in AGN light curves. In this paper we perform Monte Carlo simulations to empirically test QPO… ▽ More Active galactic nucleus (AGN) emission is dominated by stochastic, aperiodic variability which overwhelms any periodic/quasi-periodic signal (QPO) if one is present. The Auto Correlation Function (ACF) and Phase Dispersion Minimization (PDM) techniques have been used previously to claim detections of QPOs in AGN light curves. In this paper we perform Monte Carlo simulations to empirically test QPO detection feasibility in the presence of red noise. Given the community's access to large databases of monitoring light curves via large-area monitoring programmes, our goal is to provide guidance to those searching for QPOs via data trawls. We simulate evenly-sampled pure red noise light curves to estimate false alarm probabilities; false positives in both tools tend to occur towards timescales longer than (very roughly) one-third of the light curve duration. We simulate QPOs mixed with pure red noise and determine the true-positive detection sensitivity; in both tools, it depends strongly on the relative strength of the QPO against the red noise and on the steepness of the red noise PSD slope. We find that extremely large values of peak QPO power relative to red noise (typically $\sim 10^{4-5}$) are needed for a 99.7 per cent true-positive detection rate. Given that the true-positive detections using the ACF or PDM are generally rare to obtain, we conclude that period searches based on the ACF or PDM must be treated with extreme caution when the data quality is not good. We consider the feasibility of QPO detection in the context of highly-inclined, periodically self-lensing supermassive black hole binaries. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: 21 pages, 14 figures, Published in MNRAS

Journal ref: MNRAS, Volume 508, Issue 3, pp.3975-3994 (2021)

The impact of precession on the observed population of ULXs

Authors: Norman Khan, Matthew. J. Middleton, Grzegorz Wiktorowicz, Thomas Dauser, Timothy P. Roberts, Joern Wilms

Abstract: The discovery of neutron stars powering several ultraluminous X-ray sources (ULXs) raises important questions about the nature of the underlying population. In this paper we build on previous work studying simulated populations by incorporating a model where the emission originates from a precessing, geometrically beamed wind-cone, created by a super-critical inflow. We obtain estimates -- indepen… ▽ More The discovery of neutron stars powering several ultraluminous X-ray sources (ULXs) raises important questions about the nature of the underlying population. In this paper we build on previous work studying simulated populations by incorporating a model where the emission originates from a precessing, geometrically beamed wind-cone, created by a super-critical inflow. We obtain estimates -- independent of the prescription for the precession period of the wind -- for the relative number of ULXs that are potentially visible (persistent or transient) for a range of underlying factors such as the relative abundance of black holes or neutron stars within the population, maximum precessional angle, and LMXB duty cycle. We make initial comparisons to existing data using a catalogue compiled from XMM-Newton. Finally, based on estimates for the precession period, we determine how the eROSITA all-sky survey (eRASS) will be able to constrain the underlying demographic. △ Less

Submitted 21 October, 2021; originally announced October 2021.

Comments: 14 pages, 8 Figures, Accepted submission MNRAS

Long-term pulse period evolution of the ultra-luminous X-ray pulsar NGC 7793 P13

Authors: F. Fuerst, D. J. Walton, M. Heida, M. Bachetti, C. Pinto, M. J. Middleton, M. Brightman, H. P. Earnshaw, D. Barret, A. C. Fabian, P. Kretschmar, K. Pottschmidt, A. Ptak, T. Roberts, D. Stern, N. Webb, J. Wilms

Abstract: Ultra-luminous X-ray pulsars (ULXPs) provide a unique opportunity to study super-Eddington accretion. We present the results of a monitoring campaign of ULXP NGC 7793 P13. Over our four-year monitoring campaign with Swift, XMM-Newton, and NuSTAR, we measured a continuous spin-up with $\dot P$ ~ -3.8e-11 s/s. The strength of the spin-up is independent of the observed X-ray flux, indicating that des… ▽ More Ultra-luminous X-ray pulsars (ULXPs) provide a unique opportunity to study super-Eddington accretion. We present the results of a monitoring campaign of ULXP NGC 7793 P13. Over our four-year monitoring campaign with Swift, XMM-Newton, and NuSTAR, we measured a continuous spin-up with $\dot P$ ~ -3.8e-11 s/s. The strength of the spin-up is independent of the observed X-ray flux, indicating that despite a drop in observed flux in 2019, accretion onto the source has continued at largely similar rates. The source entered an apparent off-state in early 2020, which might have resulted in a change in the accretion geometry as no pulsations were found in observations in July and August 2020. We used the long-term monitoring to update the orbital ephemeris and the periodicities seen in both the observed optical/UV and X-ray fluxes. We find that the optical/UV period is very stable over the years, with $P_\text{UV}$ = 63.75 (+0.17, -0.12) d. The best-fit orbital period determined from our X-ray timing results is 64.86 +/- 0.19 d, which is almost a day longer than previously implied, and the X-ray flux period is 65.21+/- 0.15 d, which is slightly shorter than previously measured. The physical origin of these different flux periods is currently unknown. We study the hardness ratio to search for indications of spectral changes. We find that the hardness ratios at high energies are very stable and not directly correlated with the observed flux. At lower energies we observe a small hardening with increased flux, which might indicate increased obscuration through outflows at higher luminosities. We find that the pulsed fraction is significantly higher at low fluxes. This seems to imply that the accretion geometry already changed before the source entered the deep off-state. We discuss possible scenarios to explain this behavior, which is likely driven by a precessing accretion disk. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: 10 pages, 6 figures, accepted for publication in A&A

Journal ref: A&A 651, A75 (2021)

XMM-Newton campaign on the ultraluminous X-ray source NGC 247 ULX-1: outflows

Authors: C. Pinto, R. Soria, D. Walton, A. D'Ai, F. Pintore, P. Kosec, W. N. Alston, F. Fuerst, M. J. Middleton, T. P. Roberts, M. Del Santo, D. Barret, E. Ambrosi, A. Robba, H. Earnshaw, A. Fabian

Abstract: Most ULXs are believed to be powered by super-Eddington accreting neutron stars and, perhaps, black holes. Above the Eddington rate the disc is expected to thicken and to launch powerful winds through radiation pressure. Winds have been recently discovered in several ULXs. However, it is yet unclear whether the thickening of the disc or the wind variability causes the switch between the classical… ▽ More Most ULXs are believed to be powered by super-Eddington accreting neutron stars and, perhaps, black holes. Above the Eddington rate the disc is expected to thicken and to launch powerful winds through radiation pressure. Winds have been recently discovered in several ULXs. However, it is yet unclear whether the thickening of the disc or the wind variability causes the switch between the classical soft and supersoft states observed in some ULXs. In order to understand such phenomenology and the overall super-Eddington mechanism, we undertook a large (800 ks) observing campaign with XMM-Newton to study NGC 247 ULX-1, which shifts between a supersoft and classical soft ULX state. The new observations show unambiguous evidence of a wind in the form of emission and absorption lines from highly-ionised ionic species, with the latter indicating a mildly-relativistic outflow (-0.17c) in line with the detections in other ULXs. Strong dipping activity is observed in the lightcurve and primarily during the brightest observations, which is typical among soft ULXs, and indicates a close relationship between the accretion rate and the appearance of the dips. The latter is likely due to a thickening of the disc scale-height and the wind as shown by a progressively increasing blueshift in the spectral lines. △ Less

Submitted 4 June, 2021; v1 submitted 22 April, 2021; originally announced April 2021.

Comments: 18 pages, 16 figures, accepted for publication in MNRAS

Quasi-periodic dipping in the ultraluminous X-ray source, NGC 247 ULX-1

Authors: W. N. Alston, C. Pinto, D. Barret, A. D'Ai, M. Del Santo, H. Earnshaw, A. C. Fabian, F. Fuerst, E. Kara, P. Kosec, M. J. Middleton, M. L. Parker, F. Pintore, A. Robba, T. P. Roberts, R. Sathyaprakash, D. Walton, E. Ambrosi

Abstract: Most ultraluminous X-ray sources (ULXs) are believed to be stellar mass black holes or neutron stars accreting beyond the Eddington limit. Determining the nature of the compact object and the accretion mode from broadband spectroscopy is currently a challenge, but the observed timing properties provide insight into the compact object and details of the geometry and accretion processes. Here we rep… ▽ More Most ultraluminous X-ray sources (ULXs) are believed to be stellar mass black holes or neutron stars accreting beyond the Eddington limit. Determining the nature of the compact object and the accretion mode from broadband spectroscopy is currently a challenge, but the observed timing properties provide insight into the compact object and details of the geometry and accretion processes. Here we report a timing analysis for an 800 ks XMM-Newton campaign on the supersoft ultraluminous X-ray source, NGC 247 ULX-1. Deep and frequent dips occur in the X-ray light curve, with the amplitude increasing with increasing energy band. Power spectra and coherence analysis reveals the dipping preferentially occurs on $\sim 5$ ks and $\sim 10$ ks timescales. The dips can be caused by either the occultation of the central X-ray source by an optically thick structure, such as warping of the accretion disc, or from obscuration by a wind launched from the accretion disc, or both. This behaviour supports the idea that supersoft ULXs are viewed close to edge-on to the accretion disc. △ Less

Submitted 19 May, 2021; v1 submitted 22 April, 2021; originally announced April 2021.

Comments: Accepted MNRAS version following minor comments. 8 pages, 6 figures

The black hole spin in GRS 1915+105, revisited

Authors: Brianna S. Mills, Shane W. Davis, Matthew J. Middleton

Abstract: We estimate the black hole spin parameter in GRS 1915+105 using the continuum-fitting method with revised mass and inclination constraints based on the very long baseline interferometric parallax measurement of the distance to this source. We fit Rossi X-ray Timing Explorer observations selected to be accretion disk-dominated spectral states as described in McClinotck et al. (2006) and Middleton e… ▽ More We estimate the black hole spin parameter in GRS 1915+105 using the continuum-fitting method with revised mass and inclination constraints based on the very long baseline interferometric parallax measurement of the distance to this source. We fit Rossi X-ray Timing Explorer observations selected to be accretion disk-dominated spectral states as described in McClinotck et al. (2006) and Middleton et al. (2006), which previously gave discrepant spin estimates with this method. We find that, using the new system parameters, the spin in both datasets increased, providing a best-fit spin of $a_*=0.86$ for the Middleton et al. data and a poor fit for the McClintock et al. dataset, which becomes pegged at the BHSPEC model limit of $a_*=0.99$. We explore the impact of the uncertainties in the system parameters, showing that the best-fit spin ranges from $a_*= 0.4$ to 0.99 for the Middleton et al. dataset and allows reasonable fits to the McClintock et al. dataset with near maximal spin for system distances greater than $\sim 10$ kpc. We discuss the uncertainties and implications of these estimates. △ Less

Submitted 27 January, 2021; originally announced January 2021.

Comments: 12 pages, 5 figures, Submitted to ApJ

Searching for energy-resolved quasi-periodic oscillations in AGN

Authors: D. I. Ashton, M. J. Middleton

Abstract: X-ray quasi-periodic oscillations (QPOs) in AGN allow us to probe and understand the nature of accretion in highly curved space-time, yet the most robust form of detection (i.e. repeat detections over multiple observations) has been limited to a single source to-date, with only tentative claims of single observation detections in several others. The association of those established AGN QPOs with a… ▽ More X-ray quasi-periodic oscillations (QPOs) in AGN allow us to probe and understand the nature of accretion in highly curved space-time, yet the most robust form of detection (i.e. repeat detections over multiple observations) has been limited to a single source to-date, with only tentative claims of single observation detections in several others. The association of those established AGN QPOs with a specific spectral component has motivated us to search the XMM-Newton archive and analyse the energy-resolved lightcurves of 38 bright AGN. We apply a conservative false alarm testing routine folding in the uncertainty and covariance of the underlying broad-band noise. We also explore the impact of red-noise leak and the assumption of various different forms (power-law, broken power-law and lorentzians) for the underlying broad-band noise. In this initial study, we report QPO candidates in 6 AGN (7 including one tentative detection in MRK~766) from our sample of 38, which tend to be found at characteristic energies and, in four cases, at the same frequency across at least two observations, indicating they are highly unlikely to be spurious in nature. △ Less

Submitted 4 January, 2021; originally announced January 2021.

Comments: 23 pages, 19 figures. Accepted for publication in MNRAS

A New Transient Ultraluminous X-ray Source in NGC 7090

Authors: D. J. Walton, M. Heida, M. Bachetti, F. Furst, M. Brightman, H. Earnshaw, P. A. Evans, A. C. Fabian, B. W. Grefenstette, F. A. Harrison, G. L. Israel, G. B. Lansbury, M. J. Middleton, S. Pike, V. Rana, T. P. Roberts, G. A. Rodriguez Castillo, R. Salvaterra, X. Song, D. Stern

Abstract: We report on the discovery of a new, transient ultraluminous X-ray source (ULX) in the galaxy NGC 7090. This new ULX, which we refer to as NGC 7090 ULX3, was discovered via monitoring with $Swift$ during 2019-20, and to date has exhibited a peak luminosity of $L_{\rm{X}} \sim 6 \times 10^{39}$ erg s$^{-1}$. Archival searches show that, prior to its recent transition into the ULX regime, ULX3 appea… ▽ More We report on the discovery of a new, transient ultraluminous X-ray source (ULX) in the galaxy NGC 7090. This new ULX, which we refer to as NGC 7090 ULX3, was discovered via monitoring with $Swift$ during 2019-20, and to date has exhibited a peak luminosity of $L_{\rm{X}} \sim 6 \times 10^{39}$ erg s$^{-1}$. Archival searches show that, prior to its recent transition into the ULX regime, ULX3 appeared to exhibit a fairly stable luminosity of $L_{\rm{X}} \sim 10^{38}$ erg s$^{-1}$. Such strong long-timescale variability may be reminiscent of the small population of known ULX pulsars, although deep follow-up observations with $XMM$-$Newton$ and $NuSTAR$ do not reveal any robust X-ray pulsation signals. Pulsations similar to those seen from known ULX pulsars cannot be completely excluded, however, as the limit on the pulsed fraction of any signal that remains undetected in these data is $\lesssim$20\%. The broadband spectrum from these observations is well modelled with a simple thin disc model, consistent with sub-Eddington accretion, which may instead imply a moderately large black hole accretor ($M_{\rm{BH}} \sim 40 ~ M_{\odot}$). Similarly, though, more complex models consistent with the super-Eddington spectra seen in other ULXs (and the known ULX pulsars) cannot be excluded given the limited signal-to-noise of the available broadband data. The nature of the accretor powering this new ULX therefore remains uncertain. △ Less

Submitted 19 November, 2020; v1 submitted 17 November, 2020; originally announced November 2020.

Comments: 11 pages, 7 figures; accepted for publication in MNRAS

Blueshifted absorption lines from X-ray reflection in IRAS 13224-3809

Authors: A. C. Fabian, C. S. Reynolds, J. Jiang, C. Pinto, L. C. Gallo, M. L. Parker, A. N. Lasenby, W. N. Alston, D. J. K Buisson, E. M. Cackett, B. De Marco, J. Garcia, E. Kara, P. Kosec, M. J. Middleton, J. M. Miller, G. Miniutti, D. J. Walton, D. R. Wilkins, A. J. Young

Abstract: We explore a disc origin for the highly-blueshifted, variable absorption lines seen in the X-ray spectrum of the Narrow Line Seyfert 1 galaxy IRAS13224-3809. The blueshift corresponds to a velocity of about 0.25c. Such features in other Active Galactic Nuclei are often interpreted as UltraFast Outflows (UFOs). The velocity is of course present in the orbital motions of the inner disk. The absorpti… ▽ More We explore a disc origin for the highly-blueshifted, variable absorption lines seen in the X-ray spectrum of the Narrow Line Seyfert 1 galaxy IRAS13224-3809. The blueshift corresponds to a velocity of about 0.25c. Such features in other Active Galactic Nuclei are often interpreted as UltraFast Outflows (UFOs). The velocity is of course present in the orbital motions of the inner disk. The absorption lines in IRAS13224-3809 are best seen when the flux is low and the reflection component of the disk is strong relative to the power-law continuum. The spectra are consistent with a model in which the reflection component passes through a thin, highly-ionized absorbing layer at the surface of the inner disc, the blue-shifted side of which dominates the flux due to relativistic aberration (the disc inclination is about 70 deg). No fast outflow need occur beyond the disc. △ Less

Submitted 15 February, 2020; originally announced February 2020.

Comments: 6 pages, 5 figures, accepted for publication in MNRAS

A dynamic black hole corona in an active galaxy through X-ray reverberation mapping

Authors: William N. Alston, Andrew C. Fabian, Erin Kara, Michael L. Parker, Michal Dovciak, Ciro Pinto, Jiachen Jiang, Matthew J. Middleton, Giovanni Miniutti, Dominic J. Walton, Dan R. Wilkins, Douglas J. K. Buisson, Maria D. Caballero-Garcia, Edward M. Cackett, Barbara De Marco, Luigi C. Gallo, Anne M. Lohfink, Chris S. Reynolds, Phil Uttley, Andrew J. Young, Abderahmen Zogbhi

Abstract: X-ray reverberation echoes are assumed to be produced in the strongly distorted spacetime around accreting supermassive black holes. This signal allows us to spatially map the geometry of the inner accretion flow - a region which cannot yet be spatially resolved by any telescope - and provides a direct measure of the black hole mass and spin. The reverberation timescale is set by the light travel… ▽ More X-ray reverberation echoes are assumed to be produced in the strongly distorted spacetime around accreting supermassive black holes. This signal allows us to spatially map the geometry of the inner accretion flow - a region which cannot yet be spatially resolved by any telescope - and provides a direct measure of the black hole mass and spin. The reverberation timescale is set by the light travel path between the direct emission from a hot X-ray corona and the reprocessed emission from the inner edge of the accretion disc. However, there is an inherent degeneracy in the reverberation signal between black hole mass, inner disc radius and height of the illuminating corona above the disc. Here, we use a long X-ray observation of the highly-variable active galaxy, IRAS 13224-3809, to track the reverberation signal as the system evolves on timescales of a day. With the inclusion of all the relativistic effects, modelling reveals that the height of the X-ray corona increases with increasing luminosity, providing a dynamic view of the inner accretion region. This simultaneous modelling allows us to break the inherent degeneracies and obtain an independent timing-based estimate for the mass and spin of the black hole. The uncertainty on black hole mass is comparable to the leading optical reverberation method, making X-ray reverberation a powerful technique, particularly for sources with low optical variability. △ Less

Submitted 21 October, 2020; v1 submitted 17 January, 2020; originally announced January 2020.

Comments: Accepted Nat. Ast. version

Journal ref: NatAs 4, 597 (2020)

The (re)appearance of NGC 925 ULX-3, a new transient ULX

Authors: Hannah P. Earnshaw, Marianne Heida, Murray Brightman, Felix Fürst, Fiona A. Harrison, Amruta Jaodand, Matthew J. Middleton, Timothy P. Roberts, Rajath Sathyaprakash, Daniel Stern, Dominic J. Walton

Abstract: We report the discovery of a third ULX in NGC 925 (ULX-3), detected in November 2017 by Chandra at a luminosity of $L_{\rm X} = (7.8\pm0.8)\times10^{39}$ erg s$^{-1}$. Examination of archival data for NGC 925 reveals that ULX-3 was detected by Swift at a similarly high luminosity in 2011, as well as by XMM-Newton in January 2017 at a much lower luminosity of $L_{\rm X} = (3.8\pm0.5)\times10^{38}$… ▽ More We report the discovery of a third ULX in NGC 925 (ULX-3), detected in November 2017 by Chandra at a luminosity of $L_{\rm X} = (7.8\pm0.8)\times10^{39}$ erg s$^{-1}$. Examination of archival data for NGC 925 reveals that ULX-3 was detected by Swift at a similarly high luminosity in 2011, as well as by XMM-Newton in January 2017 at a much lower luminosity of $L_{\rm X} = (3.8\pm0.5)\times10^{38}$ erg s$^{-1}$. With an additional Chandra non-detection in 2005, this object demonstrates a high dynamic range of flux of factor >26. In its high-luminosity detections, ULX-3 exhibits a hard power-law spectrum with $Γ=1.6\pm0.1$, whereas the XMM-Newton detection is slightly softer, with $Γ=1.8^{+0.2}_{-0.1}$ and also well-fitted with a broadened disc model. The long-term light curve is sparsely covered and could be consistent either with the propeller effect or with a large-amplitude superorbital period, both of which are seen in ULXs, in particular those with neutron star accretors. Further systematic monitoring of ULX-3 will allow us to determine the mechanism by which ULX-3 undergoes its extreme variability and to better understand the accretion processes of ULXs. △ Less

Submitted 19 February, 2020; v1 submitted 2 January, 2020; originally announced January 2020.

Comments: 9 pages, 3 figures, 1 table. Accepted for publication in ApJ; updated to accepted version with minor revisions

The Unusual Broadband X-ray Spectral Variability of NGC 1313 X-1 seen with XMM$-$Newton, Chandra and NuSTAR

Authors: D. J. Walton, C. Pinto, M. Nowak, M. Bachetti, R. Sathyaprakash, E. Kara, T. P. Roberts, R. Soria, M. Brightman, C. R. Canizares, H. P. Earnshaw, F. Furst, M. Heida, M. J. Middleton, D. Stern, L. Tao, N. Webb, W. N. Alston, D. Barret, A. C. Fabian, F. A. Harrison, P. Kosec

Abstract: We present results from the major coordinated X-ray observing program on the ULX NGC 1313 X-1 performed in 2017, combining $XMM$-$Newton$, $Chandra$ and $NuSTAR$, focusing on the evolution of the broadband ($\sim$0.3-30.0 keV) continuum emission. Clear and unusual spectral variability is observed, but this is markedly suppressed above $\sim$10-15 keV, qualitatively similar to the ULX Holmberg IX X… ▽ More We present results from the major coordinated X-ray observing program on the ULX NGC 1313 X-1 performed in 2017, combining $XMM$-$Newton$, $Chandra$ and $NuSTAR$, focusing on the evolution of the broadband ($\sim$0.3-30.0 keV) continuum emission. Clear and unusual spectral variability is observed, but this is markedly suppressed above $\sim$10-15 keV, qualitatively similar to the ULX Holmberg IX X-1. We model the multi-epoch data with two-component accretion disc models designed to approximate super-Eddington accretion, allowing for both a black hole and a neutron star accretor. With regards to the hotter disc component, the data trace out two distinct tracks in the luminosity-temperature plane, with larger emitting radii and lower temperatures seen at higher observed fluxes. Despite this apparent anti-correlation, each of these tracks individually shows a positive luminosity-temperature relation. Both are broadly consistent with $L\propto{T}^{4}$, as expected for blackbody emission with a constant area, and also with $L\propto{T}^{2}$, as may be expected for an advection-dominated disc around a black hole. We consider a variety of possibilities for this unusual behaviour. Scenarios in which the innermost flow is suddenly blocked from view by outer regions of the super-Eddington disc/wind can explain the luminosity-temperature behaviour, but are difficult to reconcile with the lack of strong variability at higher energies, assuming this emission arises from the most compact regions. Instead, we may be seeing evidence for further radial stratification of the accretion flow than is included in the simple models considered, with a combination of winds and advection resulting in the suppressed high-energy variability. △ Less

Submitted 21 April, 2020; v1 submitted 21 November, 2019; originally announced November 2019.

Comments: 18 pages, 7 figures, accepted for publication in MNRAS

Discovery of a soft X-ray lag in the Ultraluminous X-ray Source NGC 1313 X-1

Authors: E. Kara, C. Pinto, D. J. Walton, W. N. Alston, M. Bachetti, D. Barret, M. Brightman, C. R. Canizares, H. P. Earnshaw, A. C. Fabian, F. Furst, P. Kosec, M. J. Middleton, T. P. Roberts, R. Soria, L. Tao, N. A. Webb

Abstract: Ultraluminous X-ray Sources (ULXs) provide a unique opportunities to probe the geometry and energetics of super-Eddington accretion. The radiative processes involved in super-Eddington accretion are not well understood, and so studying correlated variability between different energy bands can provide insights into the causal connection between different emitting regions. We present a spectral-timi… ▽ More Ultraluminous X-ray Sources (ULXs) provide a unique opportunities to probe the geometry and energetics of super-Eddington accretion. The radiative processes involved in super-Eddington accretion are not well understood, and so studying correlated variability between different energy bands can provide insights into the causal connection between different emitting regions. We present a spectral-timing analysis of NGC 1313 X-1 from a recent XMM-Newton campaign. The spectra can be decomposed into two thermal-like components, the hotter of which may originate from the inner accretion disc, and the cooler from an optically thick outflow. We find correlated variability between hard (2-10 keV) and soft (0.3-2 keV) bands on kilosecond timescales, and find a soft lag of ~150 seconds. The covariance spectrum suggests that emission contributing to the lags is largely associated with the hotter of the two thermal-like components, likely originating from the inner accretion flow. This is only the third ULX to exhibit soft lags. The lags range over three orders of magnitude in amplitude, but all three are ~5 to ~20 percent of the corresponding characteristic variability timescales. If these soft lags can be understood in the context of a unified picture of ULXs, then lag timescales may provide constraints on the density and extent of radiatively-driven outflows. △ Less

Submitted 21 November, 2019; originally announced November 2019.

Comments: Re-submitted to MNRAS after minor revisions, 7 pages, 5 figures

XMM-Newton Campaign On Ultraluminous X-ray Source NGC 1313 X-1: Wind vs. State Variability

Authors: C. Pinto, D. J. Walton, E. Kara, M. L. Parker, R. Soria, P. Kosec, M. J. Middleton, W. N. Alston, A. C. Fabian, M. Guainazzi, T. P. Roberts, F. Fuerst, H. P. Earnshaw, R. Sathyaprakash, D. Barret

Abstract: Most ultraluminous X-ray sources (ULXs) are thought to be powered by neutron stars and black holes accreting beyond the Eddington limit. If the compact object is a black hole or a neutron star with a magnetic field $\lesssim10^{12}$ G, the accretion disc is expected to thicken and launch powerful winds driven by radiation pressure. Evidence of such winds has been found in ULXs through the high-res… ▽ More Most ultraluminous X-ray sources (ULXs) are thought to be powered by neutron stars and black holes accreting beyond the Eddington limit. If the compact object is a black hole or a neutron star with a magnetic field $\lesssim10^{12}$ G, the accretion disc is expected to thicken and launch powerful winds driven by radiation pressure. Evidence of such winds has been found in ULXs through the high-resolution spectrometers onboard XMM-Newton, but several unknowns remain, such as the geometry and launching mechanism of these winds. In order to better understand ULX winds and their link to the accretion regime, we have undertaken a major campaign with XMM-Newton to study the ULX NGC 1313 X-1, which is known to exhibit strong emission and absorption features from a mildly-relativistic wind. The new observations show clear changes in the wind with a significantly weakened fast component (0.2c) and the rise of a new wind phase which is cooler and slower (0.06-0.08c). We also detect for the first time variability in the emission lines which indicates an origin within the accretion disc or in the wind. We describe the variability of the wind in the framework of variable super-Eddington accretion rate and discuss a possible geometry for the accretion disc. △ Less

Submitted 13 January, 2020; v1 submitted 21 November, 2019; originally announced November 2019.

Comments: 20 pages, 16 figures, 4 tables, accepted for publication in MNRAS

A rapidly-changing jet orientation in the stellar-mass black hole V404 Cygni

Authors: James C. A. Miller-Jones, Alexandra J. Tetarenko, Gregory R. Sivakoff, Matthew J. Middleton, Diego Altamirano, Gemma E. Anderson, Tomaso M. Belloni, Rob P. Fender, Peter G. Jonker, Elmar G. Körding, Hans A. Krimm, Dipankar Maitra, Sera Markoff, Simone Migliari, Kunal P. Mooley, Michael P. Rupen, David M. Russell, Thomas D. Russell, Craig L. Sarazin, Roberto Soria, Valeriu Tudose

Abstract: Powerful relativistic jets are one of the main ways in which accreting black holes provide kinetic feedback to their surroundings. Jets launched from or redirected by the accretion flow that powers them should be affected by the dynamics of the flow, which in accreting stellar-mass black holes has shown increasing evidence for precession due to frame dragging effects that occur when the black hole… ▽ More Powerful relativistic jets are one of the main ways in which accreting black holes provide kinetic feedback to their surroundings. Jets launched from or redirected by the accretion flow that powers them should be affected by the dynamics of the flow, which in accreting stellar-mass black holes has shown increasing evidence for precession due to frame dragging effects that occur when the black hole spin axis is misaligned with the orbital plane of its companion star. Recently, theoretical simulations have suggested that the jets can exert an additional torque on the accretion flow, although the full interplay between the dynamics of the accretion flow and the launching of the jets is not yet understood. Here we report a rapidly changing jet orientation on a timescale of minutes to hours in the black hole X-ray binary V404 Cygni, detected with very long baseline interferometry during the peak of its 2015 outburst. We show that this can be modelled as Lense-Thirring precession of a vertically-extended slim disk that arises from the super-Eddington accretion rate. Our findings suggest that the dynamics of the precessing inner accretion disk could play a role in either directly launching or redirecting the jets within the inner few hundred gravitational radii. Similar dynamics should be expected in any strongly-accreting black hole whose spin is misaligned with the inflowing gas, both affecting the observational characteristics of the jets, and distributing the black hole feedback more uniformly over the surrounding environment. △ Less

Submitted 12 June, 2019; originally announced June 2019.

Comments: This is a pre-print of an article published in Nature. The final published version is available online at: http://dx.doi.org/10.1038/s41586-019-1152-0 . 42 pages, 14 figures (5 as tables)

Journal ref: Nature (2019), 569, 374-377

The discovery of weak coherent pulsations in the ultraluminous X-ray source NGC 1313 X-2

Authors: R. Sathyaprakash, T. P. Roberts, D. J. Walton, F. Fuerst, M. Bachetti, C. Pinto, W. N. Alston, H. P. Earnshaw, A. C. Fabian, M. J. Middleton, R. Soria

Abstract: We report the detection of weak pulsations from the archetypal ultraluminous X-ray source (ULX) NGC 1313 X-2. Acceleration searches reveal sinusoidal pulsations in segments of two out of six new deep observations of this object, with a period of $\sim$ 1.5 s and a pulsed fraction of $\sim$ 5%. We use Monte Carlo simulations to demonstrate that the individual detections are unlikely to originate in… ▽ More We report the detection of weak pulsations from the archetypal ultraluminous X-ray source (ULX) NGC 1313 X-2. Acceleration searches reveal sinusoidal pulsations in segments of two out of six new deep observations of this object, with a period of $\sim$ 1.5 s and a pulsed fraction of $\sim$ 5%. We use Monte Carlo simulations to demonstrate that the individual detections are unlikely to originate in false Poisson noise detections given their very close frequencies; their strong similarity to other pulsations detected from ULXs also argues they are real. The presence of a large bubble nebula surrounding NGC 1313 X-2 implies an age of order 1 Myr for the accreting phase of the ULX, which implies that the neutron star's magnetic field has not been suppressed over time by accreted material, nor has the neutron star collapsed into a black hole, despite an average energy output into the nebula two orders of magnitude above Eddington. This argues that most of the accreted material has been expelled over the lifetime of the ULX, favouring physical models including strong winds and/or jets for neutron star ULXs. △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: 6 pages, 4 figures, accepted for publication by MNRAS. Comments are welcomed

Puzzling Blue Dips in the Black Hole Candidate Swift J1357.2-0933, from ULTRACAM, SALT, ATCA, Swift and NuSTAR

Authors: J. A. Paice, P. Gandhi, P. A. Charles, V. S. Dhillon, T. R. Marsh, D. A. H. Buckley, M. M. Kotze, A. Beri, D. Altamirano, M. J. Middleton, R. M. Plotkin, J. C. A. Miller-Jones, D. Russell, J. Tomsick, W. Díaz-Merced, R. Misra

Abstract: We present rapid, multiwavelength photometry of the low-mass X-ray binary Swift J1357.2-0933 during its 2017 outburst. Using several sets of quasi-simultaneous ULTRACAM/NTT (optical), NuSTAR (X-ray), XRT/Swift (X-ray), SALT (optical) and ATCA (radio) observations taken during outburst decline, we confirm the frequent optical dipping that has previously been noted both in outburst and in quiescence… ▽ More We present rapid, multiwavelength photometry of the low-mass X-ray binary Swift J1357.2-0933 during its 2017 outburst. Using several sets of quasi-simultaneous ULTRACAM/NTT (optical), NuSTAR (X-ray), XRT/Swift (X-ray), SALT (optical) and ATCA (radio) observations taken during outburst decline, we confirm the frequent optical dipping that has previously been noted both in outburst and in quiescence. We also find: 1) that the dip frequency decreases as the outburst decays, similar to what was seen in the previous outburst, 2) that the dips produce a shape similar to that in binary systems with partial disc occultations, 3) that the source becomes significantly bluer during these dips, indicating an unusual geometry compared to other LMXB dippers, and 4) that dip superposition analysis confirms the lack of an X-ray response to the optical dips. These very unusual properties appear to be unique to Swift J1357.2-0933, and are likely the result of a high binary inclination, as inferred from features such as its very low outburst X-ray luminosity. From this analysis as well as X-ray/optical timing correlations, we suggest a model with multi-component emission/absorption features with differing colours. This could include the possible presence of a sporadically occulted jet base and a recessed disc. This source still hosts many puzzling features, with consequences for the very faint X-ray transients population. △ Less

Submitted 10 June, 2019; v1 submitted 23 May, 2019; originally announced May 2019.

Comments: 15 pages, 13 figures. Revised and resubmitted to MNRAS

All at once: transient pulsations, spin down and a glitch from the Pulsating Ultraluminous X-ray Source M82 X-2

Authors: Matteo Bachetti, Thomas J. Maccarone, Murray Brightman, McKinley C. Brumback, Felix Fürst, Fiona A. Harrison, Marianne Heida, Gian Luca Israel, Matthew J. Middleton, John A. Tomsick, Natalie A. Webb, Dominic J. Walton

Abstract: M82 X-2 is the first pulsating ultraluminous X-ray source (PULX)to be identified. Since the discovery in 2014, NuSTAR has observed the M82 field 15 times throughout 2015 and 2016. In this paper, we report the results of pulsation searches in all these datasets, and find only one new detection. This new detection allows us to refine the orbital period of the source and measure an average spin down… ▽ More M82 X-2 is the first pulsating ultraluminous X-ray source (PULX)to be identified. Since the discovery in 2014, NuSTAR has observed the M82 field 15 times throughout 2015 and 2016. In this paper, we report the results of pulsation searches in all these datasets, and find only one new detection. This new detection allows us to refine the orbital period of the source and measure an average spin down rate between 2014 and 2016 of 5x10^-11 Hz/s, which is in contrast to the strong spin up seen during the 2014 observations and represent the first detection of spin down in a PULX system. Thanks to the improved orbital solution allowed by this new detection, we are also able to detect pulsations in additional segments of the original 2014 dataset. We find a glitch superimposed on the very strong and variable spin-up already reported, the first positive glitch identified in a PULX system. We discuss the new findings in the context of current leading models for PULXs. △ Less

Submitted 19 November, 2019; v1 submitted 15 May, 2019; originally announced May 2019.

Comments: 15 pages, 7 figures, 5 tables. Revised version

Discovery of a radio transient in M81

Authors: G. E. Anderson, J. C. A. Miller-Jones, M. J. Middleton, R. Soria, D. A. Swartz, R. Urquhart, N. Hurley-Walker, P. J. Hancock, R. P. Fender, P. Gandhi, S. Marko, T. P. Roberts

Abstract: We report the discovery of a radio transient in the spiral galaxy M81. The transient was detected in early 2015 as part of a two-year survey of M81 made up of 12 epochs using the Karl G. Jansky Very Large Array. While undetected on 2014 September 12, the source was first detected on 2015 January 2, from which point it remained visible at an approximately constant luminosity of… ▽ More We report the discovery of a radio transient in the spiral galaxy M81. The transient was detected in early 2015 as part of a two-year survey of M81 made up of 12 epochs using the Karl G. Jansky Very Large Array. While undetected on 2014 September 12, the source was first detected on 2015 January 2, from which point it remained visible at an approximately constant luminosity of $L_{R} = 1.5 \pm 0.1 \times 10^{24}$erg s$^{-1}$ Hz$^{-1}$ at the observing frequency of 6 GHz for at least 2 months. Assuming this is a synchrotron event with a rise-time between 2.6 and 112 days, the peak luminosity (at equipartition) corresponds to a minimum energy of $10^{44} \lesssim E_{min} \lesssim 10^{46}$erg and jet power of $P_{min} \sim 10^{39}$erg s$^{-1}$, which are higher than most known X-ray binaries. Given its longevity, lack of short-term radio variability, and the absence of any multi-wavelength counterpart (X-ray luminosity $L_{X} \lesssim 10^{36}$erg s$^{-1}$), it does not behave like known Galactic or extragalactic X-ray binaries. The M81 transient radio properties more closely resemble the unidentified radio transient 43.78+59.3 discovered in M82, which has been suggested to be a radio nebula associated with an accreting source similar to SS 433. One possibility is that both the new M81 transient and the M82 transient may be the birth of a short-lived radio bubble associated with a discrete accretion event similar to those observed from the ULX Holmberg II X-1. However, it is not possible to rule out other identifications including long-term supernova shockwave interactions with the surrounding medium from a faint supernova or a background active galaxy. △ Less

Submitted 16 September, 2019; v1 submitted 9 May, 2019; originally announced May 2019.

Comments: 18 pages, 5 figures, 6 tables, accepted for publication in MNRAS

Breaking the limit: Super-Eddington accretion onto black holes and neutron stars

Authors: M. Brightman, M. Bachetti, H. P. Earnshaw, F. Fürst, J. García, B. Grefenstette, M. Heida, E. Kara, K. K. Madsen, M. J. Middleton, D. Stern, F. Tombesi, D. J. Walton

Abstract: With the recent discoveries of massive and highly luminous quasars at high redshifts ($z\sim7$; e.g. Mortlock et al. 2011), the question of how black holes (BHs) grow in the early Universe has been cast in a new light. In order to grow BHs with $M_{\rm BH} > 10^9$ M$_{\odot}$ by less than a billion years after the Big Bang, mass accretion onto the low-mass seed BHs needs to have been very rapid (V… ▽ More With the recent discoveries of massive and highly luminous quasars at high redshifts ($z\sim7$; e.g. Mortlock et al. 2011), the question of how black holes (BHs) grow in the early Universe has been cast in a new light. In order to grow BHs with $M_{\rm BH} > 10^9$ M$_{\odot}$ by less than a billion years after the Big Bang, mass accretion onto the low-mass seed BHs needs to have been very rapid (Volonteri & Rees, 2005). Indeed, for any stellar remnant seed, the rate required would need to exceed the Eddington limit. This is the point at which the outward force produced by radiation pressure is equal to the gravitational attraction experienced by the in-falling matter. In principle, this implies that there is a maximum luminosity an object of mass $M$ can emit; assuming spherical accretion and that the opacity is dominated by Thompson scattering, this Eddington luminosity is $L_{\rm{E}} = 1.38 \times 10^{38} (M/M_{\odot})$ erg s$^{-1}$. In reality, it is known that this limit can be violated, due to non-spherical geometry or various kinds of instabilities. Nevertheless, the Eddington limit remains an important reference point, and many of the details of how accretion proceeds above this limit remain unclear. Understanding how this so-called super-Eddington accretion occurs is of clear cosmological importance, since it potentially governs the growth of the first supermassive black holes (SMBHs) and the impact this growth would have had on their host galaxies (`feedback') and the epoch of reionization, as well as improving our understanding of accretion physics more generally. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: A Science White Paper submitted to the Astro2020 Decadal Survey

Thermal stability of winds driven by radiation pressure in super-Eddington accretion disks

Authors: Ciro Pinto, Missagh Mehdipour, Dom J. Walton, Matthew J. Middleton, Tim P. Roberts, Andrew C. Fabian, Matteo Guainazzi, Roberto Soria, Peter Kosec, Jan-Uwe Ness

Abstract: Ultraluminous X-ray sources (ULXs) are mainly powered by accretion in neutron stars or stellar-mass black holes. Accreting at rates exceeding the Eddington limit by factors of a few up to hundreds, radiation pressure is expected to inflate the accretion disc, and drive fast winds that have in fact been observed at significant fractions of the speed of light. Given the super-Eddington luminosity, t… ▽ More Ultraluminous X-ray sources (ULXs) are mainly powered by accretion in neutron stars or stellar-mass black holes. Accreting at rates exceeding the Eddington limit by factors of a few up to hundreds, radiation pressure is expected to inflate the accretion disc, and drive fast winds that have in fact been observed at significant fractions of the speed of light. Given the super-Eddington luminosity, the accretion disc will be thicker than in sub-Eddington accretors such as common active galactic nuclei and X-ray binaries, leading to a different spectral energy distribution and, possibly, a different thermal status of the wind. Here we show the first attempt to calculate the photoionization balance of the winds driven by strong radiation pressure in thick discs with a focus on ULXs hosting black holes or non-magnetic neutron stars. We find that the winds are generally in thermally stable equilibrium, but long-term variations in the accretion rate and the inclination due to precession may have significant effects on the wind appearance and stability. Our model trends can explain the observed correlation between the spectral residuals around 1 keV and the ULX spectral state. We also find a possible correlation between the spectral hardness of the ULX, the wind velocity and the ionization parameter in support of the general scenario. △ Less

Submitted 5 December, 2019; v1 submitted 14 March, 2019; originally announced March 2019.

Comments: 17 pages, 9 figures, 5 tables, accepted for publication in MNRAS

The black hole X-ray transient Swift J1357.2-0933 as seen with Swift and NuSTAR during its 2017 outburst

Authors: Aru Beri, B. E. Tetarenko, A. Bahramian, Diego Altamirano, Poshak Gandhi, G. R. Sivakoff, N. Degenaar, M. J. Middleton, R. Wijnands, J. V. Hernándz Santisteban, John A. Paice

Abstract: We report on observations of black hole Swift J1357.2-0933, a member of the modest population of very faint X-ray transients. This source has previously shown intense dips in the optical lightcurve, a phenomena that has been linked to the existence of a "unique toroidal structure" in the inner region of the disc, seen at a high inclination. Our observations, carried out by the Neil Gehrels Swift a… ▽ More We report on observations of black hole Swift J1357.2-0933, a member of the modest population of very faint X-ray transients. This source has previously shown intense dips in the optical lightcurve, a phenomena that has been linked to the existence of a "unique toroidal structure" in the inner region of the disc, seen at a high inclination. Our observations, carried out by the Neil Gehrels Swift and NuSTAR X-ray observatories, do not show the presence of intense dips in the optical light curves. We find that the X-ray light curves do not show any features that would straightforwardly support an edge-on configuration or high inclination configuration of the orbit. This is similar to what was seen in the X-ray observations of the source during its 2011 outburst. Moreover, the broadband spectra were well described with an absorbed power-law model without any signatures of the cut-off at energies above 10 keV, or any reflection from the disc or the putative torus. Thus, the X-ray data do not support the unique obscuring torus scenario proposed for J1357. We also performed a multi-wavelength study using the data of X-ray telescope and Ultraviolet/Optical Telescope aboard Swift, taken during the 4.5 months duration of the 2017 outburst. This is consistent with what was previously inferred for this source. We found a correlation between the simultaneous X-ray and ultraviolet/optical data and our study suggests that most of the reprocessed flux must be coming out in the ultraviolet. △ Less

Submitted 28 February, 2019; originally announced March 2019.

Comments: 13 pages, Accepted for publication in MNRAS

A new, clean catalogue of extragalactic non-nuclear X-ray sources in nearby galaxies

Authors: H. P. Earnshaw, T. P. Roberts, M. J. Middleton, D. J. Walton, S. Mateos

Abstract: We have created a new, clean catalogue of extragalactic non-nuclear X-ray sources by correlating the 3XMM-DR4 data release of the XMM-Newton Serendipitous Source Catalogue with the Third Reference Catalogue of Bright Galaxies and the Catalogue of Neighbouring Galaxies, using an improved version of the method presented in Walton et al. (2011). Our catalogue contains 1,314 sources, of which 384 are… ▽ More We have created a new, clean catalogue of extragalactic non-nuclear X-ray sources by correlating the 3XMM-DR4 data release of the XMM-Newton Serendipitous Source Catalogue with the Third Reference Catalogue of Bright Galaxies and the Catalogue of Neighbouring Galaxies, using an improved version of the method presented in Walton et al. (2011). Our catalogue contains 1,314 sources, of which 384 are candidate ultraluminous X-ray sources (ULXs). The resulting catalogue improves upon previous catalogues in its handling of spurious detections by taking into account XMM-Newton quality flags. We estimate the contamination of ULXs by background sources to be 24 per cent. We define a 'complete' subsample as those ULXs in galaxies for which the sensitivity limit is below $10^{39}$ erg/s and use it to examine the hardness ratio properties between ULX and non-ULX sources, and ULXs in different classes of host galaxy. We find that ULXs have a similar hardness ratio distribution to lower-luminosity sources, consistent with previous studies. We also find that ULXs in spiral and elliptical host galaxies have similar distributions to each other independent of host galaxy morphology, however our results do support previous indications that the population of ULXs is more luminous in star-forming host galaxies than in non-star-forming galaxies. Our catalogue contains further interesting subpopulations for future study, including Eddington Threshold sources and highly variable ULXs. We also examine the highest-luminosity (L$_X$ > $5 \times 10^{40}$ erg/s) ULXs in our catalogue in search of intermediate-mass black hole candidates, and find nine new possible candidates. △ Less

Submitted 11 December, 2018; originally announced December 2018.

Comments: 20 pages, 14 figures, 7 tables. Accepted for publication by MNRAS

NuSTAR reveals the hidden nature of SS433

Authors: M. J. Middleton, D. J. Walton, W. Alston, T. Dauser, S. Eikenberry, Y-F Jiang, A. C. Fabian, F. Fuerst, M. Brightman, H. Marshall, M. Parker, C. Pinto, F. A. Harrison, M. Bachetti, D. Altamirano, A. J. Bird, G. Perez, J. Miller-Jones, P. A. Charles, S. Boggs, F. Christensen, W. Craig, K. Forster, B. Grefenstette, C. Hailey , et al. (3 additional authors not shown)

Abstract: SS433 is the only Galactic binary system known to persistently accrete at highly super-critical (or hyper-critical) rates, similar to those in tidal disruption events, and likely needed to explain the rapid growth of those very high redshift quasars containing massive SMBHs. Probing the inner regions of SS433 in the X-rays is crucial to understanding this system, and super-critical accretion in ge… ▽ More SS433 is the only Galactic binary system known to persistently accrete at highly super-critical (or hyper-critical) rates, similar to those in tidal disruption events, and likely needed to explain the rapid growth of those very high redshift quasars containing massive SMBHs. Probing the inner regions of SS433 in the X-rays is crucial to understanding this system, and super-critical accretion in general, but is highly challenging due to obscuration by the surrounding wind, driven from the accretion flow. NuSTAR observed SS433 in the hard X-ray band across multiple phases of its 162 day super-orbital precession period. Spectral-timing tools allow us to infer that the hard X-ray emission from the inner regions is likely being scattered towards us by the walls of the wind-cone. By comparing to numerical models, we determine an intrinsic X-ray luminosity of $\ge$ 2$\times$10$^{37}$ erg/s and that, if viewed face on, we would infer an apparent luminosity of $>$ 1$\times$10$^{39}$ erg/s, confirming SS433's long-suspected nature as an ultraluminous X-ray source (ULX). We present the discovery of a narrow, $\sim 100$s lag due to atomic processes occurring in outflowing material travelling at least 0.14-0.29c, which matches absorption lines seen in ULXs and -- in the future -- will allow us to map a super-critical outflow for the first time. △ Less

Submitted 4 May, 2021; v1 submitted 24 October, 2018; originally announced October 2018.

Comments: 14 pages, 10 figures, accepted for publication in MNRAS

A tale of two periods: determination of the orbital ephemeris of the super-Eddington pulsar NGC 7793 P13

Authors: F. Fuerst, D. J. Walton, M. Heida, F. A. Harrison, D. Barret, M. Brightman, A. C. Fabian, M. J. Middleton, C. Pinto, V. Rana, F. Tramper, N. Webb, P. Kretschmar

Abstract: We present a timing analysis of multiple XMM-Newton and NuSTAR observations of the ultra-luminous pulsar NGC 7793 P13 spread over its 65d variability period. We use the measured pulse periods to determine the orbital ephemeris, confirm a long orbital period with P_orb = 63.9 (+0.5,-0.6) d, and find an eccentricity of e <= 0.15. The orbital signature is imprinted on top of a secular spin-up, which… ▽ More We present a timing analysis of multiple XMM-Newton and NuSTAR observations of the ultra-luminous pulsar NGC 7793 P13 spread over its 65d variability period. We use the measured pulse periods to determine the orbital ephemeris, confirm a long orbital period with P_orb = 63.9 (+0.5,-0.6) d, and find an eccentricity of e <= 0.15. The orbital signature is imprinted on top of a secular spin-up, which seems to get faster as the source becomes brighter. We also analyse data from dense monitoring of the source with Swift and find an optical photometric period of 63.9 +/- 0.5 d and an X-ray flux period of 66.8 +/- 0.4 d. The optical period is consistent with the orbital period, while the X-ray flux period is significantly longer. We discuss possible reasons for this discrepancy, which could be due to a super-orbital period caused by a precessing accretion disk or an orbital resonance. We put the orbital period of P13 into context with the orbital periods implied for two other ultra-luminous pulsars, M82 X-2 and NGC 5907 ULX and discuss possible implications for the system parameters. △ Less

Submitted 14 May, 2018; originally announced May 2018.

Comments: 10 pages, 6 figures, accepted for publication in A&A

Journal ref: A&A 616, A186 (2018)

The remarkable X-ray variability of IRAS 13224-3809 I: The variability process

Authors: W. N. Alston, A. C. Fabian, D. J. K. Buisson, E. Kara, M. L. Parker, A. M. Lohfink, P. Uttley, D. R. Wilkins, C. Pinto, B. De Marco, E. M. Cackett, M. J. Middleton, D. J. Walton, C. S. Reynolds, J. Jiang, L. C. Gallo, A. Zogbhi, G. Miniutti, M. Dovciak, A. J. Young

Abstract: We present a detailed X-ray timing analysis of the highly variable NLS1 galaxy, IRAS 13224-3809. The source was recently monitored for 1.5 Ms with XMM-Newton which, combined with 500 ks archival data, makes this the best studied NLS1 galaxy in X-rays to date. We apply standard time- and Fourier-domain in order to understand the underlying variability process. The source flux is not distributed log… ▽ More We present a detailed X-ray timing analysis of the highly variable NLS1 galaxy, IRAS 13224-3809. The source was recently monitored for 1.5 Ms with XMM-Newton which, combined with 500 ks archival data, makes this the best studied NLS1 galaxy in X-rays to date. We apply standard time- and Fourier-domain in order to understand the underlying variability process. The source flux is not distributed lognormally, as would be expected for accreting sources. The first non-linear rms-flux relation for any accreting source in any waveband is found, with $\mathrm{rms} \propto \mathrm{flux}^{2/3}$. The light curves exhibit significant strong non-stationarity, in addition to that caused by the rms-flux relation, and are fractionally more variable at lower source flux. The power spectrum is estimated down to $\sim 10^{-7}$ Hz and consists of multiple peaked components: a low-frequency break at $\sim 10^{-5}$ Hz, with slope $α< 1$ down to low frequencies; an additional component breaking at $\sim 10^{-3}$ Hz. Using the high-frequency break we estimate the black hole mass $M_\mathrm{BH} = [0.5-2] \times 10^{6} M_{\odot}$, and mass accretion rate in Eddington units, $\dot m_{\rm Edd} \gtrsim 1$. The non-stationarity is manifest in the PSD with the normalisation of the peaked components increasing with decreasing source flux, as well as the low-frequency peak moving to higher frequencies. We also detect a narrow coherent feature in the soft band PSD at $0.7$ mHz, modelled with a Lorentzian the feature has $Q \sim 8$ and an $\mathrm{rms} \sim 3$ %. We discuss the implication of these results for accretion of matter onto black holes. △ Less

Submitted 17 September, 2018; v1 submitted 28 March, 2018; originally announced March 2018.

Comments: Accepted to MNRAS. 19 pages, 17 figures

A Potential Cyclotron Resonant Scattering Feature in the ULX Pulsar NGC 300 ULX1 seen by NuSTAR and XMM-Newton

Authors: D. J. Walton, M. Bachetti, F. Fuerst, D. Barret, M. Brightman, A. C. Fabian, B. W. Grefenstette, F. A. Harrison, M. Heida, J. Kennea, P. Kosec, R. M. Lau, K. K. Madsen, M. J. Middleton, C. Pinto, J. F. Steiner, N. Webb

Abstract: Based on phase-resolved broadband spectroscopy using $XMM$-$Newton$ and $NuSTAR$, we report on a potential cyclotron resonant scattering feature at $E \sim 13$ keV in the pulsed spectrum of the recently discoverd ULX pulsar NGC 300 ULX1. If this interpretation is correct, the implied magnetic field of the central neutron star is $B \sim 10^{12}$ G (assuming scattering off electrons), similar to th… ▽ More Based on phase-resolved broadband spectroscopy using $XMM$-$Newton$ and $NuSTAR$, we report on a potential cyclotron resonant scattering feature at $E \sim 13$ keV in the pulsed spectrum of the recently discoverd ULX pulsar NGC 300 ULX1. If this interpretation is correct, the implied magnetic field of the central neutron star is $B \sim 10^{12}$ G (assuming scattering off electrons), similar to that estimated from the observed spin-up of the star, and also similar to known Galactic X-ray pulsars. We discuss the implications of this result for the connection between NGC 300 ULX1 and the other known ULX pulsars, particularly in light of the recent discovery of a likely proton Cyclotron line in another ULX, M51 ULX-8. △ Less

Submitted 4 April, 2018; v1 submitted 20 March, 2018; originally announced March 2018.

Comments: 6 pages, 3 figures. Accepted for publication in ApJL

Evidence for Pulsar-like Emission Components in the Broadband ULX Sample

Authors: D. J. Walton, F. Fuerst, M. Heida, F. A. Harrison, D. Barret, D. Stern, M. Bachetti, M. Brightman, A. C. Fabian, M. J. Middleton

Abstract: We present broadband X-ray analyses of a sample of bright ultraluminous X-ray sources with the goal of investigating the spectral similarity of this population to the known ULX pulsars, M82 X-2, NGC7793 P13 and NGC5907 ULX. We perform a phase-resolved analysis of the broadband XMM-Newton+NuSTAR dataset of NGC5907 ULX, finding that the pulsed emission from the accretion column in this source exhibi… ▽ More We present broadband X-ray analyses of a sample of bright ultraluminous X-ray sources with the goal of investigating the spectral similarity of this population to the known ULX pulsars, M82 X-2, NGC7793 P13 and NGC5907 ULX. We perform a phase-resolved analysis of the broadband XMM-Newton+NuSTAR dataset of NGC5907 ULX, finding that the pulsed emission from the accretion column in this source exhibits a similar spectral shape to that seen in both M82 X-2 and NGC7793 P13, and that this is responsible for the excess emission observed at the highest energies when the spectra are fit with accretion disk models. We then demonstrate that similar 'hard' excesses are seen in all the ULXs in the broadband sample. Finally, for the ULXs where the nature of the accretor is currently unknown, we test whether the hard excesses are all consistent with being produced by an accretion column similar to those present in M82 X-2, NGC7793 P13 and NGC5907 ULX. Based on the average shape of the pulsed emission, we find that in all cases a similar accretion column can successfully reproduce the observed data, consistent with the hypothesis that this ULX sample may be dominated by neutron star accretors. Compared to the known pulsar ULXs, our spectral fits for the remaining ULXs suggest that the non-pulsed emission from the accretion flow beyond the magnetosphere makes a stronger relative contribution than the component associated with the accretion column. If these sources do also contain neutron star accretors, this may help to explain the lack of detected pulsations. △ Less

Submitted 12 March, 2018; originally announced March 2018.

Comments: 13 pages, 6 figures; accepted for publication in ApJ

Magnetic field strength of a neutron-star-powered ultraluminous X-ray source

Authors: Murray Brightman, Fiona A. Harrison, Felix Fuerst, Matthew J. Middleton, Dominic J. Walton, Daniel Stern, Andrew C. Fabian, Marianne Heida, Didier Barret, Matteo Bachetti

Abstract: Ultraluminous X-ray sources (ULXs) are bright X-ray sources in nearby galaxies not associated with the central supermassive black hole. Their luminosities imply they are powered by either an extreme accretion rate onto a compact stellar remnant, or an intermediate mass ($\sim100-10^5$ M$_{\odot}$) black hole. The recent detection of coherent pulsations coming from three bright ULXs demonstrates th… ▽ More Ultraluminous X-ray sources (ULXs) are bright X-ray sources in nearby galaxies not associated with the central supermassive black hole. Their luminosities imply they are powered by either an extreme accretion rate onto a compact stellar remnant, or an intermediate mass ($\sim100-10^5$ M$_{\odot}$) black hole. The recent detection of coherent pulsations coming from three bright ULXs demonstrates that some of these sources are powered by accretion onto a neutron star, implying accretion rates significantly in excess of the Eddington limit, a high degree of geometric beaming, or both. The physical challenges associated with the high implied accretion rates can be mitigated if the neutron star surface field is high - in the magnetar regime ($10^{14}$ G), since this suppresses the electron scattering cross section, reducing the radiation pressure that chokes off accretion for high luminosities. One of the few ways to determine surface magnetic fields is through the detection of cyclotron resonance scattering features (CRSFs) produced by the transition of charged particles between quantized Landau levels. To date, CRSFs have only been detected in Galactic accreting pulsars. Here we present the detection at 3.8-σ significance of a strong absorption line at a rest-frame energy of 4.5 keV in the Chandra spectrum of a ULX in M51. We find that this feature is likely to be a CRSF produced by the strong magnetic field of a neutron star. Assuming scattering off electrons, the magnetic field strength is implied to be $\sim10^{11}$ G, however the line is narrower than any electron CRSFs previously observed, and assuming thermal broadening, the implied temperature is significantly cooler than the accretion disk or column. The line shape is, however, consistent with a proton resonance scattering feature, implying that the neutron star has a magnetic field near the surface of B$\sim10^{15}$ G. △ Less

Submitted 6 March, 2018; originally announced March 2018.

Comments: Author version of the paper published in Nature Astronomy on 26 February 2018. See free-to-view published version here: http://rdcu.be/HQpR

Is there a UV/X-ray connection in IRAS 13224-3809?

Authors: D. J. K. Buisson, A. M. Lohfink, W. N. Alston, E. M. Cackett, C. -Y. Chiang, T. Dauser, B. De Marco, A. C. Fabian, L. C. Gallo, J. A. García, J. Jiang, E. Kara, M. J. Middleton, G. Miniutti, M. L. Parker, C. Pinto, P. Uttley, D. J. Walton, D. R. Wilkins

Abstract: We present results from the optical, ultraviolet and X-ray monitoring of the NLS1 galaxy IRAS 13224-3809 taken with Swift and XMM-Newton during 2016. IRAS 13224-3809 is the most variable bright AGN in the X-ray sky and shows strong X-ray reflection, implying that the X-rays strongly illuminate the inner disc. Therefore, it is a good candidate to study the relationship between coronal X-ray and dis… ▽ More We present results from the optical, ultraviolet and X-ray monitoring of the NLS1 galaxy IRAS 13224-3809 taken with Swift and XMM-Newton during 2016. IRAS 13224-3809 is the most variable bright AGN in the X-ray sky and shows strong X-ray reflection, implying that the X-rays strongly illuminate the inner disc. Therefore, it is a good candidate to study the relationship between coronal X-ray and disc UV emission. However, we find no correlation between the X-ray and UV flux over the available ~40 day monitoring, despite the presence of strong X-ray variability and the variable part of the UV spectrum being consistent with irradiation of a standard thin disc. This means either that the X-ray flux which irradiates the UV emitting outer disc does not correlate with the X-ray flux in our line of sight and/or that another process drives the majority of the UV variability. The former case may be due to changes in coronal geometry, absorption or scattering between the corona and the disc. △ Less

Submitted 18 December, 2017; originally announced December 2017.

Comments: 9 pages, 8 figures, MNRAS submitted: this version in response to reviewers comments; in original form 2017 October 2

Lense-Thirring precession in ULXs as a possible means to constrain the neutron star equation-of-state

Authors: M. J. Middleton, P. C. Fragile, M. Bachetti, M. Brightman, Y-F. Jiang, W. C. G. Ho, T. P. Roberts, A. R. Ingram, T. Dauser, C. Pinto, D. J. Walton, F. Fuerst, A. C. Fabian, N. Gehrels

Abstract: The presence of neutron stars in at least three ultraluminous X-ray sources is now firmly established and offers an unambiguous view of super-critical accretion. All three systems show long-timescale periods (60-80 days) in the X-rays and/or optical, two of which are known to be super-orbital in nature. Should the flow be classically super critical, i.e. the Eddington limit is reached locally in t… ▽ More The presence of neutron stars in at least three ultraluminous X-ray sources is now firmly established and offers an unambiguous view of super-critical accretion. All three systems show long-timescale periods (60-80 days) in the X-rays and/or optical, two of which are known to be super-orbital in nature. Should the flow be classically super critical, i.e. the Eddington limit is reached locally in the disc (implying surface dipole fields that are sub-magnetar in strength), then the large scale-height flow can precess through the Lense-Thirring effect which could provide an explanation for the observed super-orbital periods. By connecting the details of the Lense-Thirring effect with the observed pulsar spin period, we are able to infer the moment-of-inertia and therefore equation-of-state of the neutron star without relying on the inclination of, or distance to the system. We apply our technique to the case of NGC 7793 P13 and demonstrate that stronger magnetic fields imply stiffer equations of state. We discuss the caveats and uncertainties, many of which can be addressed through forthcoming radiative magnetohydrodynamic (RMHD) simulations and their connection to observation. △ Less

Submitted 20 November, 2017; originally announced November 2017.

Comments: 14 pages, 9 figures, accepted to MNRAS

Paving the way to simultaneous multi-wavelength astronomy

Authors: M. J. Middleton, P. Casella, P. Gandhi, E. Bozzo, G. Anderson, N. Degenaar, I. Donnarumma, G. Israel, C. Knigge, A. Lohfink, S. Markoff, T. Marsh, N. Rea, S. Tingay, K. Wiersema, D. Altamirano, D. Bhattacharya, W. N. Brandt, S. Carey, P. Charles, M. Diaz Trigo, C. Done, M. Kotze, S. Eikenberry, R. Fender , et al. (27 additional authors not shown)

Abstract: Whilst astronomy as a science is historically founded on observations at optical wavelengths, studying the Universe in other bands has yielded remarkable discoveries, from pulsars in the radio, signatures of the Big Bang at submm wavelengths, through to high energy emission from accreting, gravitationally-compact objects and the discovery of gamma-ray bursts. Unsurprisingly, the result of combinin… ▽ More Whilst astronomy as a science is historically founded on observations at optical wavelengths, studying the Universe in other bands has yielded remarkable discoveries, from pulsars in the radio, signatures of the Big Bang at submm wavelengths, through to high energy emission from accreting, gravitationally-compact objects and the discovery of gamma-ray bursts. Unsurprisingly, the result of combining multiple wavebands leads to an enormous increase in diagnostic power, but powerful insights can be lost when the sources studied vary on timescales shorter than the temporal separation between observations in different bands. In July 2015, the workshop "Paving the way to simultaneous multi-wavelength astronomy" was held as a concerted effort to address this at the Lorentz Center, Leiden. It was attended by 50 astronomers from diverse fields as well as the directors and staff of observatories and spaced-based missions. This community white paper has been written with the goal of disseminating the findings of that workshop by providing a concise review of the field of multi-wavelength astronomy covering a wide range of important source classes, the problems associated with their study and the solutions we believe need to be implemented for the future of observational astronomy. We hope that this paper will both stimulate further discussion and raise overall awareness within the community of the issues faced in a developing, important field. △ Less

Submitted 21 September, 2017; v1 submitted 11 September, 2017; originally announced September 2017.

Comments: 52 pages, 15 figures, accepted, invited review (to appear in New Astronomy Reviews), v3: updated figure and text

Super-Eddington Accretion onto the Neutron Star NGC 7793 P13: Broadband X-ray Spectroscopy and Ultraluminous X-ray Sources

Authors: D. J. Walton, F. Fuerst, F. A. Harrison, D. Stern, M. Bachetti, D. Barret, M. Brightman, A. C. Fabian, M. J. Middleton, A. Ptak, L. Tao

Abstract: We present a detailed, broadband X-ray spectral analysis of the ULX pulsar NGC 7793 P13, a known super-Eddington source, utilizing data from the $XMM$-$Newton$, $NuSTAR$ and $Chandra$ observatories. The broadband $XMM$-$Newton+NuSTAR$ spectrum of P13 is qualitatively similar to the rest of the ULX sample with broadband coverage, suggesting that additional ULXs in the known population may host neut… ▽ More We present a detailed, broadband X-ray spectral analysis of the ULX pulsar NGC 7793 P13, a known super-Eddington source, utilizing data from the $XMM$-$Newton$, $NuSTAR$ and $Chandra$ observatories. The broadband $XMM$-$Newton+NuSTAR$ spectrum of P13 is qualitatively similar to the rest of the ULX sample with broadband coverage, suggesting that additional ULXs in the known population may host neutron star accretors. Through time-averaged, phase-resolved and multi-epoch studies, we find that two non-pulsed thermal blackbody components with temperatures $\sim$0.5 and $\sim$1.5 keV are required to fit the data below 10 keV, in addition to a third continuum component which extends to higher energies and is associated with the pulsed emission from the accretion column. The characteristic radii of the thermal components appear to be similar, and are too large to be associated with the neutron star itself, so the need for two components likely indicates the accretion flow outside the magnetosphere is complex. We suggest a scenario in which the thick inner disc expected for super-Eddington accretion begins to form, but is terminated by the neutron star's magnetic field soon after its onset, implying a limit of $B \lesssim 6 \times 10^{12}$ G for the dipolar component of the central neutron star's magnetic field. Evidence of similar termination of the disc in other sources may offer a further means of identifying additional neutron star ULXs. Finally, we examine the spectrum exhibited by P13 during one of its unusual 'off' states. These data require both a hard powerlaw component, suggesting residual accretion onto the neutron star, and emission from a thermal plasma, which we argue is likely associated with the P13 system. △ Less

Submitted 31 October, 2017; v1 submitted 29 May, 2017; originally announced May 2017.

Comments: 17 pages, 12 figures, accepted for publication in MNRAS

The response of relativistic outflowing gas to the inner accretion disk of a black hole

Authors: Michael L. Parker, Ciro Pinto, Andrew C. Fabian, Anne Lohfink, Douglas J. K. Buisson, William Alston, Erin Kara, Edward M. Cackett, Chia-Ying Chiang, Thomas Dauser, Barbara De Marco, Luigi C. Gallo, Javier Garcia, Fiona A. Harrison, Ashley L. King, Matthew J. Middleton, Jon M. Miller, Giovanni Miniutti, Christopher S. Reynolds, Phil Uttley, Ranjan Vasudevan, Dominic J. Walton, Daniel R. Wilkins, Abderahmen Zoghbi

Abstract: Active galactic nucleus (AGN) feedback is the process by which supermassive black holes in the centres of galaxies may moderate the growth of their hosts. Gas outflows from supermassive black holes release huge quantities of energy into the interstellar medium, clearing the surrounding gas. The most extreme of these, the ultra-fast outflows (UFOs), are the subset of X-ray detected outflows with ve… ▽ More Active galactic nucleus (AGN) feedback is the process by which supermassive black holes in the centres of galaxies may moderate the growth of their hosts. Gas outflows from supermassive black holes release huge quantities of energy into the interstellar medium, clearing the surrounding gas. The most extreme of these, the ultra-fast outflows (UFOs), are the subset of X-ray detected outflows with velocities higher than 10,000 km/s, believed to originate in relativistic disc winds, a few hundred gravitational radii from the black hole. The absorption features produced by these outflows are variable, but no clear link has been found between the behaviour of the X-ray continuum and the energy or equivalent width of the outflow features due to the long time-scales of quasar variability. Here, we present the detection of multiple absorption lines from an extreme ultra-fast gas flow in the X-ray spectrum of the active galactic nucleus IRAS 13224-3809, at 0.236+/-0.006 times the speed of light (71,000 km/s), where the absorption is strongly anti-correlated with the emission from the inner regions of the accretion disk. If the gas flow is identified as a genuine outflow then it is in the fastest 5 per cent of such winds, and its variability is hundreds of times faster than in other variable winds, allowing us to observe in hours what would take months in a quasar. We find signatures of the wind simultaneously in both low and high energy detectors, which are consistent with a single ionized outflow, linking the two phenomena. The detection of the wind responding to the emission from the inner disk demonstrates a connection between accretion processes occurring on very different scales, with the X-rays from within a few gravitational radii of the black hole ionizing the relativistically outflowing gas as the flux rises. △ Less

Submitted 28 February, 2017; originally announced March 2017.

Comments: 13 pages, 7 figures, published in Nature 02/03/17

From ultraluminous X-ray sources to ultraluminous supersoft sources: NGC 55 ULX, the missing link

Authors: C. Pinto, W. Alston, R. Soria, M. J. Middleton, D. J. Walton, A. D. Sutton, A. C. Fabian, H. Earnshaw, R. Urquhart, E. Kara, T. P. Roberts

Abstract: In recent work with high-resolution grating spectrometers (RGS) aboard XMM-Newton Pinto et al. (2016) have discovered that two bright and archetypal ultraluminous X-ray sources (ULXs) have strong relativistic winds in agreement with theoretical predictions of high accretion rates. It has been proposed that such winds can become optically thick enough to block and reprocess the disk X-ray photons a… ▽ More In recent work with high-resolution grating spectrometers (RGS) aboard XMM-Newton Pinto et al. (2016) have discovered that two bright and archetypal ultraluminous X-ray sources (ULXs) have strong relativistic winds in agreement with theoretical predictions of high accretion rates. It has been proposed that such winds can become optically thick enough to block and reprocess the disk X-ray photons almost entirely, making the source appear as a soft thermal emitter or ultraluminous supersoft X-ray source (ULS). To test this hypothesis we have studied a ULX where the wind is strong enough to cause significant absorption of the hard X-ray continuum: NGC 55 ULX. The RGS spectrum of NGC 55 ULX shows a wealth of emission and absorption lines blueshifted by significant fractions of the light speed (0.01 - 0.20)c indicating the presence of a powerful wind. The wind has a complex dynamical structure with the ionization state increasing with the outflow velocity, which may indicate launching from different regions of the accretion disk. The comparison with other ULXs such as NGC 1313 X-1 and NGC 5408 X-1 suggests that NGC 55 ULX is being observed at higher inclination. The wind partly absorbs the source flux above 1 keV, generating a spectral drop similar to that observed in ULSs. The softening of the spectrum at lower (~ Eddington) luminosities and the detection of a soft lag agree with the scenario of wind clumps crossing the line of sight, partly absorbing and reprocessing the hard X-rays from the innermost region. △ Less

Submitted 14 March, 2017; v1 submitted 16 December, 2016; originally announced December 2016.

Comments: 20 pages, 20 figures, 7 tables, accepted for publication in MNRAS after moderate review. Conclusions have not been changed

Crossing the Eddington limit: examining disk spectra at high accretion rates

Authors: Andrew D. Sutton, Douglas A. Swartz, Timothy P. Roberts, Matthew J. Middleton, Roberto Soria, Chris Done

Abstract: The faintest ultraluminous X-ray sources (ULXs), those with 0.3-10 keV luminosities 1 < L_X/10^39 < 3 erg s^-1, tend to have X-ray spectra that are disk-like but broader than expected for thin accretion disks. These `broadened disk' spectra are thought to indicate near- or mildly super-Eddington accretion onto stellar remnant black holes. Here we report that a sample of bright thermal-dominant bla… ▽ More The faintest ultraluminous X-ray sources (ULXs), those with 0.3-10 keV luminosities 1 < L_X/10^39 < 3 erg s^-1, tend to have X-ray spectra that are disk-like but broader than expected for thin accretion disks. These `broadened disk' spectra are thought to indicate near- or mildly super-Eddington accretion onto stellar remnant black holes. Here we report that a sample of bright thermal-dominant black hole binaries, which have Eddington ratios constrained to moderate values, also show broadened disk spectra in the 0.3-10 keV band at an order of magnitude lower luminosities. This broadening would be missed in studies that only look above ~2 keV. While this may suggest that broadened disk ULXs could be powered by accretion onto massive stellar remnant black holes with close to maximal spin, we argue in favor of a scenario where they are at close to the Eddington luminosity, such that radiation pressure would be expected to result in geometrically slim, advective accretion disks. However, this implies that an additional physical mechanism is required to produce the observed broad spectra at low Eddington ratios. △ Less

Submitted 28 November, 2016; originally announced November 2016.

Comments: 20 pages, 7 figures, accepted for publication in ApJ

The Broadband Spectral Variability of Holmberg IX X-1

Authors: D. J. Walton, F. Fuerst, F. A. Harrison, M. J. Middleton, A. C. Fabian, M. Bachetti, D. Barret, J. M. Miller, A. Ptak, V. Rana, D. Stern, L. Tao

Abstract: We present results from four new broadband X-ray observations of the extreme ultraluminous X-ray source Holmberg IX X-1 ($L_{\rm{X}} > 10^{40}$ erg/s), performed by $Suzaku$ and $NuSTAR$ in coordination. Combined with the archival data, we now have broadband observations of this remarkable source from six separate epochs. Two of these new observations probe lower fluxes than seen previously, allow… ▽ More We present results from four new broadband X-ray observations of the extreme ultraluminous X-ray source Holmberg IX X-1 ($L_{\rm{X}} > 10^{40}$ erg/s), performed by $Suzaku$ and $NuSTAR$ in coordination. Combined with the archival data, we now have broadband observations of this remarkable source from six separate epochs. Two of these new observations probe lower fluxes than seen previously, allowing us to extend our knowledge of the broadband spectral variability exhibited. The spectra are well fit by two thermal blackbody components, which dominate the emission below 10 keV, as well as a steep ($Γ\sim 3.5$) powerlaw tail which dominates above $\sim$15 keV. Remarkably, while the 0.3-10.0 keV flux varies by a factor of $\sim$3 between all these epochs, the 15-40 keV flux varies by only $\sim$20%. Although the spectral variability is strongest in the $\sim$1-10 keV band, both of the thermal components are required to vary when all epochs are considered. We also re-visit the search for iron absorption features, leveraging the high-energy $NuSTAR$ data to improve our sensitivity to extreme velocity outflows in light of the ultra-fast outflow recently detected in NGC 1313 X-1. Iron absorption from a similar outflow along our line of sight can be ruled out in this case. We discuss these results in the context of super-Eddington accretion models that invoke a funnel-like geometry for the inner flow, and propose a scenario in which we have an almost face-on view of a funnel that expands to larger radii with increasing flux, resulting in an increasing degree of geometrical collimation for the emission from intermediate temperature regions. △ Less

Submitted 20 March, 2017; v1 submitted 20 October, 2016; originally announced October 2016.

Comments: 11 pages, 6 Figures, accepted for publication in ApJ

Discovery of coherent pulsations from the Ultraluminous X-ray Source NGC 7793 P13

Authors: F. Fuerst, D. J. Walton, F. A. Harrison, D. Stern, D. Barret, M. Brightman, A. C. Fabian, B. Grefenstette, K. K. Madsen, M. J. Middleton, J. M. Miller, K. Pottschmidt, A. Ptak, V. Rana, N. Webb

Abstract: We report the detection of coherent pulsations from the ultraluminous X-ray source NGC 7793 P13. The ~0.42s nearly sinusoidal pulsations were initially discovered in broadband X-ray observations using XMM-Newton and NuSTAR taken in 2016. We subsequently also found pulsations in archival XMM-Newton data taken in 2013 and 2014. The significant (>>5 sigma) detection of coherent pulsations demonstrate… ▽ More We report the detection of coherent pulsations from the ultraluminous X-ray source NGC 7793 P13. The ~0.42s nearly sinusoidal pulsations were initially discovered in broadband X-ray observations using XMM-Newton and NuSTAR taken in 2016. We subsequently also found pulsations in archival XMM-Newton data taken in 2013 and 2014. The significant (>>5 sigma) detection of coherent pulsations demonstrates that the compact object in P13 is a neutron star with an observed peak luminosity of ~1e40 erg/s (assuming isotropy), well above the Eddington limit for a 1.4 M_sun accretor. This makes P13 the second ultraluminous X-ray source known to be powered by an accreting neutron star. The pulse period varies between epochs, with a slow but persistent spin up over the 2013-2016 period. This spin-up indicates a magnetic field of B ~ 1.5e12 G, typical of many accreting pulsars. The most likely explanation for the extreme luminosity is a high degree of beaming, however this is difficult to reconcile with the sinusoidal pulse profile. △ Less

Submitted 17 October, 2016; v1 submitted 22 September, 2016; originally announced September 2016.

Comments: 6 pages, 3 figures, accepted for publication in ApJL