Search

Soft Ultra-Luminous X-ray (ULXs) sources are a subclass of the ULXs that can switch from a supersoft spectral state, where most of the luminosity is emitted below 1 keV, to a soft spectral state with significant emission above 1 keV. In a few systems, dips have been observed. The mechanism behind this state transition and the dips nature are still debated. To investigate these issues, we obtained a long XMM-Newton monitoring campaign of a member of this class, NGC 247 ULX-1. We computed the hardness-intensity diagram for the whole dataset and identified two different branches: the normal branch and the dipping branch, which we study with four and three hardness-intensity resolved spectra, respectively. All seven spectra are well described by two thermal components: a colder ($kT_{\rm bb}$ $\sim$ 0.1-0.2 keV) black-body, interpreted as emission from the photo-sphere of a radiatively-driven wind, and a hotter ($kT_{\rm disk}$ $\sim$ 0.6 keV) multicolour disk black-body, likely due to reprocessing of radiation emitted from the innermost regions. In addition, a complex pattern of emission and absorption lines has been taken into account based on previous high-resolution spectroscopic results. We studied the evolution of spectral parameters and the flux of the two thermal components along the two branches and discuss two scenarios possibly connecting the state transition and the dipping phenomenon. One is based on geometrical occultation of the emitting regions, the other invokes the onset of a propeller effect., 13 pages, accepted by MNRAS on 2021/08/16

Swift J1858.6-0814 is a recently discovered X-ray binary notable for extremely strong variability (by factors $>100$ in soft X-rays) in its discovery state. We present the detection of five thermonuclear (Type I) X-ray bursts from Swift J1858.6-0814, implying that the compact object in the system is a neutron star. Some of the bursts show photospheric radius expansion, so their peak flux can be used to estimate the distance to the system. The peak luminosity, and hence distance, can depend on several system parameters; for the most likely values, a high inclination and a helium atmosphere, $D=12.8_{-0.6}^{+0.8}$ kpc, although systematic effects allow a conservative range of $9-18$ kpc. Before one burst, we detect a QPO at $9.6\pm0.5$ mHz with a fractional rms amplitude of $2.2\pm0.2$% ($0.5-10$ keV), likely due to marginally stable burning of helium; similar oscillations may be present before the other bursts but the light curves are not long enough to allow their detection. We also search for burst oscillations but do not detect any, with an upper limit in the best case of 15% fractional amplitude (over $1-8$ keV). Finally, we discuss the implications of the neutron star accretor and this distance on other inferences which have been made about the system. In particular, we find that Swift J1858.6-0814 was observed at super-Eddington luminosities at least during bright flares during the variable stage of its outburst., 14 pages, 8 figures, 2 tables, MNRAS accepted

[Context] It is thought that ultraluminous X-ray sources (ULXs) are mainly powered by super-Eddington accreting neutron stars or black holes as shown by the recent discovery of X-ray pulsations and relativistic winds. [Aims] This work presents a follow-up study of the spectral evolution over two decades of the pulsing ULX NGC 1313 X-2 in order to understand the structure of the accretion disc. The primary objective is to determine the shape and nature of the dominant spectral components by investigating their variability with the changes in the source luminosity. [Methods[ We performed a spectral analysis over the canonical 0.3-10.0 keV energy band of all the high signal-to-noise XMM-Newton observations (96% of the available data), and we tested a number of different spectral models, which should approximate super-Eddington accretion discs. The baseline model consists of two thermal blackbody components with different temperatures plus an exponential cutoff powerlaw. [Results] The baseline model provides a good description of the X-ray spectra. In particular, the hotter and brighter (LXa., a., 6-9a ×a 1039 erg s-1) thermal component describes the emission from the super-Eddington inner disc and the cutoff powerlaw describes the contribution from the accretion column of the neutron star. Instead, the cooler component describes the emission from the outer region of the disc close to the spherisation radius and the wind. The luminosity-temperature relation for the cool component follows a negative trend, which is not consistent with La, T4, as is expected from a sub-Eddington thin disc of Shakura-Sunayev. This is not consistent with La, T2 either, as is expected for an advection-dominated disc. However, this would rather agree with a wind-dominated X-ray emitting region. Instead, the (Lx, Tdisk) relation for the hotter component is somewhere in between the first two theoretical scenarios. [Conclusions] Our findings agree with the super-Eddington scenario and provide further detail on the disc structure. The source spectral evolution is qualitatively similar to that seen in NGC 1313 X-1 and Holmberg IX X-1, indicating a common structure and evolution among archetypal ULXs.

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

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., Comment: 17 pages, 9 figures, 5 tables, accepted for publication in MNRAS

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 $\alpha < 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., Comment: Accepted to MNRAS. 19 pages, 17 figures

Most Ultraluminous X-ray sources (ULXs) are thought to be powered by super-Eddington accretion onto stellar-mass compact objects. Accretors in this extreme regime are naturally expected to ionise copious amounts of plasma in their vicinity and launch powerful radiation-driven outflows from their discs. High spectral resolution X-ray observations (with RGS gratings onboard XMM-Newton) of a few ULXs with the best datasets indeed found complex line spectra and confirmed such extreme (0.1-0.3c) winds. However, a search for plasma signatures in a large ULX sample with a rigorous technique has never been performed, thereby preventing us from understanding their statistical properties such as the rate of occurrence, to constrain the outflow geometry and its duty cycle. We developed a fast method for automated line detection in X-ray spectra and applied it to the full RGS ULX archive, rigorously quantifying the statistical significance of any candidate lines. Collecting the 135 most significant features detected in 89 observations of 19 objects, we created the first catalogue of spectral lines detected in soft X-ray ULX spectra. We found that the detected emission lines are concentrated around known rest-frame elemental transitions and thus originate from low-velocity material. The absorption lines instead avoid these transitions, suggesting they were imprinted by blueshifted outflows. Such winds therefore appear common among the ULX population. Additionally, we found that spectrally hard ULXs show fewer line detections than soft ULXs, indicating some difference in their accretion geometry and orientation, possibly causing over-ionisation of plasma by the harder spectral energy distributions of harder ULXs., Accepted for publication in MNRAS. 22 pages, 6 figures and 3 tables

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., Comment: Accepted MNRAS version following minor comments. 8 pages, 6 figures

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., 11 pages, 7 figures; accepted for publication in MNRAS

We present the results of the latest NuSTAR monitoring campaign of the Compton-thick Seyfert 2 galaxy NGC 1068, composed of four ∼50 ks observations performed between 2017 July and 2018 February to search for flux and spectral variability on time-scales from 1 to 6 months. We detect one unveiling and one eclipsing event with time-scales less than 27 and 91 d, respectively, ascribed to Compton-thick material with NH = (1.8 ± 0.8) × 1024 and ≥ (2.4 ± 0.5) × 1024 cm−2 moving across our line of sight. This gas is likely located in the innermost part of the torus or even further inward, thus providing further evidence of the clumpy structure of the circumnuclear matter in this source. Taking advantage of simultaneous Swift–XRT observations, we also detected a new flaring ULX, at a distance d ∼ 30 arcsec (i.e. ∼2 kpc) from the nuclear region of NGC 1068, with a peak X-ray intrinsic luminosity of (3.0 ± 0.4) × 1040 erg s−1 in the 2–10 keV band., AZ, SB, and GM acknowledge financial support from the Italian Space Agency under grant no. 2017-14-H.O. SB acknowledges financial support from the Italian Space Agency under grant ASI-INAF I/037/12/0. AM and GM acknowledge financial support from the Italian Space Agency under grant ASI-INAF I/037/12/0-011/13. FEB acknowledges support from CONICYT-Chile (Basal AFB-170002) and the Ministry of Economy, Development, and Tourism’s Millennium Science Initiative through grant IC120009, awarded to the Millennium Institute of Astrophysics, MAS. CR acknowledges support from the CONICYT + PAI Convocatoria Nacional subvencion a instalacion en la academia convocatoria año 2017 PAI77170080.

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., Re-submitted to MNRAS after minor revisions, 7 pages, 5 figures

Relativistically blueshifted absorption features of highly ionised ions, the so-called ultra-fast outflows (UFOs), have been detected in the X-ray spectra of a number of accreting supermassive black holes. If these features truly originate from accretion disc winds accelerated to more than 10 per cent of the speed of light, their energy budget is very significant and they can contribute to or even drive galaxy-scale feedback from active galactic nuclei (AGN). However, the UFO spectral features are often weak due to high ionisation of the outflowing material, and the inference of the wind physical properties can be complicated by other spectral features in AGN such as relativistic reflection. Here we study a highly accreting Narrow Line Seyfert 1 galaxy PG 1448+273. We apply an automated, systematic routine for detecting outflows in accreting systems and achieve an unambiguous detection of a UFO in this AGN. The UFO absorption is observed in both soft and hard X-ray bands with the XMM-Newton observatory. The velocity of the outflow is (26900 +- 600) km/s (~0.09c), with an ionisation parameter of log ({\xi} / erg cm s^-1)=4.03_{-0.08}^{+0.10} and a column density above 10^23 cm^-2. At the same time, we detect weak warm absorption features in the spectrum of the object. Our systematic outflow search suggests the presence of further multi-phase wind structure, but we cannot claim a significant detection considering the present data quality. The UFO is not detected in a second, shorter observation with XMM-Newton, indicating variability in time, observed also in other similar AGN., Comment: Accepted for publication in MNRAS. 14 pages, 11 figures, 1 table

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

We present two new NuSTAR observations of the narrow line Seyfert 1 (NLS1) galaxy Mrk 766 and give constraints on the two scenarios previously proposed to explain its spectrum and that of other NLS1s: relativistic reflection and partial covering. The NuSTAR spectra show a strong hard (>15 keV) X-ray excess, while simultaneous soft X-ray coverage of one of the observations provided by XMM-Newton constrains the ionised absorption in the source. The pure reflection model requires a black hole of high spin ($a>0.92$) viewed at a moderate inclination ($i=46^{+1}_{-4}$ degrees). The pure partial covering model requires extreme parameters: the cut-off of the primary continuum is very low ($22^{+7}_{-5}$ keV) in one observation and the intrinsic X-ray emission must provide a large fraction (75%) of the bolometric luminosity. Allowing a hybrid model with both partial covering and reflection provides more reasonable absorption parameters and relaxes the constraints on reflection parameters. The fractional variability reduces around the iron K band and at high energies including the Compton hump, suggesting that the reflected emission is less variable than the continuum., 14 pages, 10 figures, MNRAS accepted

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., 14 pages, 9 figures, accepted to MNRAS

We present results from a coordinated $XMM$-$Newton$+$NuSTAR$ observation of the type 1.8 Seyfert galaxy IRAS 13197-1627. This is a highly complex source, with strong contributions from relativistic reflection from the inner accretion disk, neutral absorption and further reprocessing by more distant material, and ionised absorption from an outflow. We undertake a detailed spectral analysis combining the broadband coverage provided by $XMM$-$Newton$+$NuSTAR$ with a multi-epoch approach incorporating archival observations performed by $XMM$-$Newton$ and $Suzaku$. Our focus is on characterising the reflection from the inner accretion disk, which previous works have suggested may dominate the AGN emission, and constraining the black hole spin. Using lamppost disk reflection models, we find that the results for the inner disk are largely insensitive to assumptions regarding the geometry of the distant reprocessor and the precise form of the illuminating X-ray continuum. However, these results do depend on the treatment of the iron abundance of the distant absorber/reprocessor. The multi-epoch data favour a scenario in which the AGN is chemically homogeneous, and we find that a rapidly rotating black hole is preferred, with $a^* \geq 0.7$, but a slowly-rotating black hole is not strongly excluded. In addition to the results for the inner disk, we also find that both the neutral and ionised absorbers vary from epoch to epoch, implying that both have some degree of inhomogeneity in their structure., 15 pages, 7 figures, accepted for publication in MNRAS

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state., Comment: 50 pages, 30 figures, uses aastex62.cls. Accepted for publication in ApJ, 07/06/2019. High-level products page in MAST will go live after 7/15/2019. Replaced Figure 4 on 7/12/2019 to be more red/green color-blind friendly

We present results from a deep, coordinated $XMM$-$Newton$+$NuSTAR$ observation of the Seyfert 2 galaxy IRAS 00521-7054. The $NuSTAR$ data provide the first detection of this source in high-energy X-rays ($E > 10$ keV), and the broadband data show this to be a highly complex source which exhibits relativistic reflection from the inner accretion disc, further reprocessing by more distant material, neutral absorption, and evidence for ionised absorption in an extreme, ultrafast outflow ($v_{\rm{out}} \sim 0.4c$). Based on lamppost disc reflection models, we find evidence that the central supermassive black hole is rapidly rotating ($a > 0.77$), consistent with previous estimates from the profile of the relativistic iron line, and that the accretion disc is viewed at a fairly high inclination ($i \sim 59^{\circ}$). Based on extensive simulations, we find the ultrafast outflow is detected at $\sim$4$��$ significance (or greater). We also estimate that the extreme outflow should be sufficient to power galaxy-scale feedback, and may even dominate the energetics of the total output from the system., 12 pages, 8 figures; accepted for publication in MNRAS

The reported discovery of a cyclotron resonance scattering feature (CRSF) in the spectrum of M51 ULX-8 may provide an important clue as to the nature of the magnetic field in those ultraluminous X-ray sources hosting neutron stars. In this paper we present the covariance (linearly correlated variability) spectrum of M51 ULX-8 on long (> 2000s) timescales. This allows us to unambiguously decompose the spectrum which requires multiple components in order to explain the broad-band emission. Having a clearer picture of the spectral decomposition leads to various tests for the dipole field strength of the neutron star which can be extended to other ULXs when certain criteria are met. In the case of M51 ULX-8, we rule out a very strong (~10^15 G) dipole solution with either a sub- or super-critical disc. Instead, our tests indicate an upper limit on the dipole field of ~10^12 G, and a classical super-critical inflow, similar to that inferred in other ULXs found to harbour neutron stars, although we do not rule out the presence of an additional, strong (~10^15 G) multipole field falling off steeply with distance from the neutron star., 9 pages, 4 figures, accepted to MNRAS, updated following corrections

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., 6 pages, 4 figures, accepted for publication by MNRAS. Comments are welcomed

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., 20 pages, 14 figures, 7 tables. Accepted for publication by MNRAS

We discovered 2.8 s pulsations in the X-ray emission of the ultraluminous X-ray source (ULX) M51 ULX-7 within the UNSEeN project, which was designed to hunt for new pulsating ULXs (PULXs) with XMM-Newton. The pulse shape is sinusoidal and large variations of its amplitude were observed even within single exposures (pulsed fraction from less than 5% to 20%). M51 ULX-7 is a variable source, generally observed at an X-ray luminosity between $10^{39}$ and $10^{40}$ erg s$^{-1}$, located in the outskirts of the spiral galaxy M51a at a distance of 8.6 Mpc. According to our analysis, the X-ray pulsar orbits in a 2-d binary with a projected semi-major axis $a_\mathrm{X} \sin i \simeq$ 28 lt-s. For a neutron star (NS) of 1.4 $M_{\odot}$, this implies a lower limit on the companion mass of 8 $M_{\odot}$, placing the system hosting M51 ULX-7 in the high-mass X-ray binary class. The barycentric pulse period decreased by $\simeq$0.4 ms in the 31 d spanned by our May -- June 2018 observations, corresponding to a spin-up rate $\dot{P} \simeq -1.5\times10^{-10}\text{s s}^{-1}$. In an archival 2005 XMM-Newton exposure, we measured a spin period of $\sim$3.3 s, indicating a secular spin-up of $\dot{P}_{\mathrm{sec}}\simeq -10^{-9}\text{ s s}^{-1}$, a value in the range of other known PULXs. Our findings suggest that the system consists of an OB giant and a moderately magnetic (dipole field component in the range $10^{12}$ G $\lesssim B_{\mathrm{dip}}\lesssim 10^{13}$G) accreting NS with weakly beamed emission ($1/12\lesssim b\lesssim1/4$)., Comment: Accepted for publication in The Astrophysical Journal

We collate active galactic nuclei (AGN) with reported detections of both relativistic reflection and ultra-fast outflows. By comparing the inclination of the inner disc from reflection with the line-of-sight velocity of the outflow, we show that it is possible to meaningfully constrain the geometry of the absorbing material. We find a clear relation between the velocity and inclination, and demonstrate that it can potentially be explained either by simple wind geometries or by absorption from the disc surface. Due to systematic errors and a shortage of high- quality simultaneous measurements our conclusions are tentative, but this study represents a proof-of-concept that has great potential., 5 pages, 3 figures, accepted to MNRAS letters

We report on the NuSTAR observations of two bright Seyfert 1 galaxies, namely MCG +8-11-11 (100 ks) and NGC 6814 (150 ks). The main goal of these observations was to investigate the Comptonization mechanisms acting in the innermost regions of AGN which are believed to be responsible for the UV/X-ray emission. The spectroscopic analysis of the NuSTAR spectra of these two sources revealed that although they had different properties overall (black hole masses, luminosity and Eddington ratios) they had very similar coronal properties. Both presented a power law spectrum with a high-energy cutoff at $\sim 150-200$ keV, a relativistically broadened Fe K$\alpha$ line and the associated disk reflection component, plus a narrow iron line likely emitted in Compton thin and distant matter. The intrinsic continuum was well described by Comptonization models that show for MCG +8-11-11 a temperature of the coronal plasma of kT$_{\rm e} \sim$ 60 keV and an extrapolated optical depth $\tau$=1.8; for NGC 6814 the coronal temperature was kT$_{\rm e} \sim$ 45 keV with an extrapolated optical depth of $\tau$=2.5. We compare and discuss these values to some most common Comptonization models which aim at explaining the energy production and stability of coronae in active galactic nuclei.

International audience; Based on phase-resolved broadband spectroscopy using XMM-Newton and NuSTAR, we report on a potential cyclotron resonant scattering feature (CRSF) at E ∼ 13 keV in the pulsed spectrum of the recently discovered ultraluminous X-ray source (ULX) pulsar NGC 300 ULX1. If this interpretation is correct, the implied magnetic field of the central neutron star is B ∼ 1012 G (assuming scattering by 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.

We report on the serendipitous discovery of a new transient in NGC 5907, at a peak luminosity of 6.4x10^{39} erg/s. The source was undetected in previous 2012 Chandra observations with a 3 sigma upper limit on the luminosity of 1.5x10^{38} erg/s, implying a flux increase of a factor of >35. We analyzed three recent 60ks/50ks Chandra and 50ks XMM-Newton observations, as well as all the available Swift observations performed between August 2017/March 2018. Until the first half of October 2017, Swift observations do not show any emission from the source. The transient entered the ULX regime in less than two weeks and its outburst was still on-going at the end of February 2018. The 0.3-10 keV spectrum is consistent with a single multicolour blackbody disc (kT~1.5 keV). The source might be a ~30 solar mass black hole accreting at the Eddington limit. However, although we did not find evidence of pulsations, we cannot rule-out the possibility that this ULX hosts an accreting neutron star., Accepted on MNRAS, 5 pages, 2 figure, 1 table

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, 10 pages, 6 figures, accepted for publication in A&A

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-{\sigma} 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., Comment: Author version of the paper published in Nature Astronomy on 26 February 2018. See free-to-view published version here: http://rdcu.be/HQpR

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., 13 pages, 6 figures; accepted for publication in ApJ

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 ($\Gamma \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., Comment: 11 pages, 6 Figures, accepted for publication in ApJ

PG1247+267 is one of the most luminous known quasars at z similar to 2 and is a strongly super-Eddington accreting supermassive black hole (SMBH) candidate. We obtained NuSTAR data of this intriguing source in December 2014 with the aim of studying its high-energy emission, leveraging the broad band covered by the new NuSTAR and the archival XMM-Newton data. Several measurements are in agreement with the super-Eddington scenario for PG1247+267: the soft power law (Gamma = 2.3 +/- 0.1); the weak ionized Fe emission line; and a hint of the presence of outflowing ionized gas surrounding the SMBH. The presence of an extreme reflection component is instead at odds with the high accretion rate proposed for this quasar. This can be explained with three different scenarios; all of them are in good agreement with the existing data, but imply very different conclusions: i) a variable primary power law observed in a low state, superimposed on a reflection component echoing a past, higher flux state; ii) a power law continuum obscured by an ionized, Compton thick, partial covering absorber; and iii) a relativistic disk reflector in a lamp-post geometry, with low coronal height and high BH spin. The first model is able to explain the high reflection component in terms of variability. The second does not require any reflection to reproduce the hard emission, while a rather low high-energy cutoff of similar to 100 keV is detected for the first time in such a high redshift source. The third model require a face-on geometry, which may affect the SMBH mass and Eddington ratio measurements. Deeper X-ray broad-band data are required in order to distinguish between these possibilities.

We present a broad-band (~0.3-70 keV) spectral and temporal analysis of NuSTAR observations of the luminous infrared galaxy NGC 6240, combined with archival Chandra, XMM-Newton and BeppoSAX data. NGC 6240 is a galaxy in a relatively early merger state with two distinct nuclei separated by ~1."5. Previous Chandra observations have resolved the two nuclei, showing that they are both active and obscured by Compton-thick material. Although they cannot be resolved by NuSTAR, thanks to the unprecedented quality of the NuSTAR data at energies >10 keV, we clearly detect, for the first time, both the primary and the reflection continuum components. The NuSTAR hard X-ray spectrum is dominated by the primary continuum piercing through an absorbing column density which is mildly optically thick to Compton scattering (tau ~ 1.2, N_H ~ 1.5 x 10^(24) cm^-2). We detect moderate hard X-ray (> 10 keV) flux variability up to 20% on short (15-20 ksec) timescales. The amplitude of the variability is maximum at ~30 keV and is likely to originate from the primary continuum of the southern nucleus. Nevertheless, the mean hard X-ray flux on longer timescales (years) is relatively constant. Moreover, the two nuclei remain Compton-thick, although we find evidence of variability of the material along the line of sight with column densities N_H, Comment: 14 pages, 10 figures, accepted for publication in A&A

(Abridged) Our deep multiwavelength campaign on NGC 5548 revealed an unusually strong X-ray obscuration. The resulting dramatic decrease in incident ionizing flux allowed us to construct a comprehensive physical, spatial and temporal picture for the long-studied AGN wind in this object. Here we analyze the UV spectra of the outflow acquired during the campaign as well as from four previous epochs. We find that a simple model based on a fixed total column-density absorber, reacting to changes in ionizing illumination, matches the very different ionization states seen in five spectroscopic epochs spanning 16 years. Absorption troughs from C III* appeared for the first time during our campaign. From these troughs, we infer that the main outflow component is situated at 3.5+-1 pc from the central source. Three other components are situated between 5-70 pc and two are further than 100 pc. The wealth of observational constraints and the disparate relationship of the observed X-ray and UV flux between different epochs make our physical model a leading contender for interpreting trough variability data of quasar outflows., Comment: Submitted to A&A, 13 pages and 7 figures, plus 27 pages of 3 appendix figures

MCG-6-30-15, at a distance of 37 Mpc (z=0.008), is the archetypical Seyfert 1 galaxy showing very broad Fe K$\alpha$ emission. We present results from a joint NuSTAR and XMM-Newton observational campaign that, for the first time, allows a sensitive, time-resolved spectral analysis from 0.35 keV up to 80 keV. The strong variability of the source is best explained in terms of intrinsic X-ray flux variations and in the context of the light bending model: the primary, variable emission is reprocessed by the accretion disk, which produces secondary, less variable, reflected emission. The broad Fe K$\alpha$ profile is, as usual for this source, well explained by relativistic effects occurring in the innermost regions of the accretion disk around a rapidly rotating black hole. We also discuss the alternative model in which the broadening of the Fe K$\alpha$ is due to the complex nature of the circumnuclear absorbing structure. Even if this model cannot be ruled out, it is disfavored on statistical grounds. We also detected an occultation event likely caused by BLR clouds crossing the line of sight., Comment: 19 pages, 16 figures, accepted for publication on ApJ

Results are presented from a 500ks long XMM-Newton observation of the Narrow-Line Seyfert 1 galaxy IRAS13224-3809. The source is rapidly variable on timescales down to a few 100s. The spectrum shows strong broad Fe-K and L emission features which are interpreted as arising from reflection from the inner parts of an accretion disc around a rapidly spinning black hole. Assuming a power-law emissivity for the reflected flux and that the innermost radius corresponds to the innermost stable circular orbit, the black hole spin is measured to be 0.988 with a statistical precision better than one per cent. Systematic uncertainties are discussed. A soft X-ray lag of 100s confirms this scenario. The bulk of the power-law continuum source is located at a radius of 2-3 gravitational radii., 7 pages, 14 figures, submitted to MNRAS

As part of a long term monitoring campaign of Mrk 335, deep XMM-Newton observations catch the narrow-line Seyfert 1 galaxy (NLS1) in a complex, intermediate flux interval as the active galaxy is transiting from low- to high-flux. Other works on these same data examined the general behaviour of the NLS1 (Grupe et al.) and the conditions of its warm absorber (Longinotti et al.). The analysis presented here demonstrates the X-ray continuum and timing properties can be described in a self-consistent manner adopting a blurred reflection model with no need to invoke partial covering. The rapid spectral variability appears to be driven by changes in the shape of the primary emitter that is illuminating the inner accretion disc around a rapidly spinning black hole (a > 0.7). While light bending is certainly prominent, the rather constant emissivity profile and break radius obtained in our spectral fitting suggest that the blurring parameters are not changing as would be expected if the primary source is varying its distance from the disc. Instead changes could be intrinsic to the power law component. One possibility is that material in an unresolved jet above the disc falls to combine with material at the base of the jet producing the changes in the primary emitter (spectral slope and flux) without changing its distance from the disc., 10 pages. Accepted for publication in MNRAS

We present a X-ray spectral analysis of a large sample of 25 'bare' active galactic nuclei, sources with little or no complicating intrinsic absorption, observed with Suzaku. Our work focuses on studying the potential contribution from relativistic disc reflection, and examining the implications of this interpretation for the intrinsic spectral complexities frequently displayed by AGN in the X-ray bandpass. During the analysis, we take the unique approach of attempting to simultaneously undertake a systematic analysis of the whole sample, as well as a detailed treatment of each individual source, and find that disc reflection has the required flexibility to successfully reproduce the broadband spectrum observed for all of the sources considered. Where possible, we use the reflected emission to place constraints on the black hole spin for this sample of sources. Our analysis suggests a general preference for rapidly rotating black holes, which if taken at face value is most consistent with the scenario in which SMBH growth is dominated by prolonged, ordered accretion. However, there may be observational biases towards AGN with high spin in the compiled sample, limiting our ability to draw strong conclusions for the general population at this stage. Finally, contrary to popular belief, our analysis also implies that the dichotomy between radio loud/radio quiet AGN is not solely related to black hole spin., Comment: 21 pages, 10 figures. Accepted for publication in MNRAS

We report on Chandra observations of the bright ultra-luminous X-ray (ULX) source in NGC3921. Previous XMM observations reported in the literature show the presence of a bright ULX at a 0.5-10 keV luminosity of 2x10^40 erg/s. Our Chandra observation finds the source at a lower luminosity of ~8x10^39 erg/s, furthermore, we provide a Chandra position of the ULX accurate to 0.7" at 90% confidence. The X-ray variability makes it unlikely that the high luminosity is caused by several separate X-ray sources. In 3 epochs of archival Hubble Space Telescope (HST) observations we find a candidate counterpart to the ULX. There is direct evidence for variability between the two epochs of WFPC2 F814W observations with the observation obtained in 2000 showing a brighter source. Furthermore, converting the 1994 F336W and 2000 F300W WFPC2 and the 2010 F336W WFC3 observations to the Johnson U-band filter assuming a spectral type of O7I we find evidence for a brightening of the U-band light in 2000. Using the higher resolution WFC3 observations the candidate counterpart is resolved into two sources of similar color. We discuss the nature of the ULX and the probable association with the optical counterpart(s). Finally, we investigate a potential new explanation for some (bright) ULXs as the decaying stages of flares caused by the tidal disruption of a star by a recoiled supermassive black hole. However, we find that there should be at most only 1 of such systems within z=0.08., 9 pages, 2 figures, accepted for publication in ApJ

The Circinus galaxy is one of the nearest obscured AGN, making it an ideal target for detailed study. Combining archival Chandra and XMM-Newton data with new NuSTAR observations, we model the 2-79 keV spectrum to constrain the primary AGN continuum and to derive physical parameters for the obscuring material. Chandra's high angular resolution allows a separation of nuclear and off-nuclear galactic emission. In the off-nuclear diffuse emission we find signatures of strong cold reflection, including high equivalent-width neutral Fe lines. This Compton-scattered off-nuclear emission amounts to 18% of the nuclear flux in the Fe line region, but becomes comparable to the nuclear emission above 30 keV. The new analysis no longer supports a prominent transmitted AGN component in the observed band. We find that the nuclear spectrum is consistent with Compton-scattering by an optically-thick torus, where the intrinsic spectrum is a powerlaw of photon index $\Gamma = 2.2-2.4$, the torus has an equatorial column density of $N_{\rm H} = (6-10)\times10^{24}$cm$^{-2}$ and the intrinsic AGN $2-10$ keV luminosity is $(2.3-5.1)\times 10^{42}$ erg/s. These values place Circinus along the same relations as unobscured AGN in accretion rate-vs-$\Gamma$ and $L_X$-vs-$L_{IR}$ phase space. NuSTAR's high sensitivity and low background allow us to study the short time-scale variability of Circinus at X-ray energies above 10 keV for the first time. The lack of detected variability favors a Compton-thick absorber, in line with the the spectral fitting results.

We present X-ray spectral analyses for three Seyfert 2 active galactic nuclei (AGNs), NGC 424, NGC 1320, and IC 2560, observed by NuSTAR in the 3-79 keV band. The high quality hard X-ray spectra allow detailed modeling of the Compton reflection component for the first time in these sources. Using quasi-simultaneous NuSTAR and Swift/XRT data, as well as archival XMM-Newton data, we find that all three nuclei are obscured by Compton-thick material with column densities in excess of similar to 5 x 10(24) cm(-2), and that their X-ray spectra above 3 keV are dominated by reflection of the intrinsic continuum on Compton-thick material. Due to the very high obscuration, absorbed intrinsic continuum components are not formally required by the data in any of the sources. We constrain the intrinsic photon indices and the column density of the reflecting medium through the shape of the reflection spectra. Using archival multi-wavelength data we recover the intrinsic X-ray luminosities consistent with the broadband spectral energy distributions. Our results are consistent with the reflecting medium being an edge-on clumpy torus with a relatively large global covering factor and overall reflection efficiency of the order of 1%. Given the unambiguous confirmation of the Compton-thick nature of the sources, we investigate whether similar sources are likely to be missed by commonly used selection criteria for Compton-thick AGNs, and explore the possibility of finding their high-redshift counterparts.