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$).

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.

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 $\Gamma=1.6\pm0.1$, whereas the XMM-Newton detection is slightly softer, with $\Gamma=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.

Ultraluminous X-ray sources (ULXs) are a class of accreting compact objects with X-ray luminosities above 1e39 erg/s. The ULX population counts several hundreds objects but only a minor fraction is well studied. Here we present a detailed analysis of all ULXs hosted in the galaxy NGC 7456. It was observed in X-rays only once in the past (in 2005) by XMM-Newton, but the observation was short and strongly affected by high background. In 2018, we obtained a new, deeper (~90 ks) XMM-Newton observation that allowed us to perform a detailed characterization of the ULXs hosted in the galaxy. ULX-1 and ULX-2, the two brightest objects (Lx~(6-10)e39 erg/s), have spectra that can be described by a two-thermal component model as often found in ULXs. ULX-1 shows also one order of magnitude in flux variability on short-term timescales (hundreds to thousand ks). The other sources (ULX-3 and ULX-4) show flux changes of at least an order of magnitude, and these objects may be candidate transient ULXs although longer X-ray monitoring or further studies are required to ascribe them to the ULX population. In addition, we found a previously undetected source that might be a new candidate ULX (labelled as ULX-5) with a luminosity of ~1e39 erg/s and hard power-law spectral shape, whose nature is still unclear and for which a background Active Galactic Nucleus cannot be excluded. We discuss the properties of all the ULXs in NGC 7456 within the framework of super-Eddington accretion onto stellar mass compact objects. Although no pulsations were detected, we cannot exclude that the sources host neutron stars.

M82 hosts two well-known ultraluminous X-ray sources (ULXs). X-1, an intermediate-mass black hole (IMBH) candidate, and X-2, an ultraluminous X-ray pulsar (ULXP). Here we present a broadband X-ray spectral analysis of both sources based on ten observations made simultaneously with Chandra and NuSTAR. Chandra provides the high spatial resolution to resolve the crowded field in the 0.5--8 keV band, and NuSTAR provides the sensitive hard X-ray spectral data, extending the bandpass of our study above 10 keV. The observations, taken in the period 2015--2016, cover a period of flaring from X-1, allowing us to study the spectral evolution of this source with luminosity. During four of these observations, X-2 was found to be at a low flux level, allowing an unambiguous view of the emission from X-1. We find that the broadband X-ray emission from X-1 is consistent with that seen in other ULXs observed in detail with NuSTAR, with a spectrum that includes a broadened disk-like component and a high-energy tail. We find that the luminosity of the disk scales with inner disk temperature as L~T^-3/2 contrary to expectations of a standard accretion disk and previous results. These findings rule out a thermal state for sub-Eddington accretion and therefore do not support M82 X-1 as an IMBH candidate. We also find evidence that the neutral column density of the material in the line of sight increases with L$_X$, perhaps due to an increased mass outflow with accretion rate. For X-2, we do not find any significant spectral evolution, but we find the spectral parameters of the phase-averaged broadband emission are consistent with the pulsed emission at the highest X-ray luminosities.

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.

We present the results from a monitoring campaign made with the Neil Gehrels Swift Observatory of the M51 galaxies, which contain several variable ultraluminous X-ray sources (ULXs). The ongoing campaign started in May 2018, and we report here on $\sim1.5$ years of observations. The campaign, which consists of 105 observations, has a typical cadence of 3--6 days, and has the goal of determining the long-term X-ray variability of the ULXs. Two of the most variable sources were ULX7 and ULX8, both of which are known to be powered by neutron stars that are exceeding their isotropic Eddington luminosities by factors of up to 100. This is further evidence that neutron star powered ULXs are the most variable. Our two main results are, first, that ULX7 exhibits a periodic flux modulation with a period of 38 days varying over a magnitude and a half in flux from peak to trough. Since the orbital period of the system is known to be 2 days, the modulation is super-orbital, which is a near-ubiquitous property of ULX pulsars. Secondly we identify a new transient ULX, M51 XT-1, the onset of which occurred during our campaign, reaching a peak luminosity of $\sim10^{40}$ erg s$^{-1}$, before gradually fading over the next $\sim200$ days until it slipped below the detection limit of our observations. Combined with the high-quality Swift/XRT lightcurve of the transient, serendipitous observations made with Chandra and XMM-Newton provide insights into the onset and evolution of a likely super-Eddington event.

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.

Motivated by the recent discoveries that six Ultraluminous X-ray Sources (ULXs) are powered by highly super-Eddington X-ray pulsars, we searched for additional pulsating ULX (PULX) candidates by identifying sources that exhibit long-term flux variability of at least an order of magnitude (a common feature seen in the 6 known PULXs, which may potentially be related to transitions to the propeller regime). Expanding on previous studies, we used the available fluxes from XMM-Newton, Swift and Chandra, along with carefully computed upper limits in cases of a non-detection, to construct long-term lightcurves for a sample of 296 ULXs selected from the XMM-Newton archive. Among these 296, we find 25 sources showing flux variability larger than a factor of 10, of which 17 show some evidence for (or are at least consistent with) exhibiting bi-modal flux distributions, as would be expected for sources undergoing propeller transitions. These sources are excellent candidates for continued monitoring programs to further test for this behaviour. There are 3 sources in our final sample with fluxes similar to NGC 5907 ULX1, currently the faintest known PULX, which would also be good targets for deeper observations with current facilities to search for pulsations. For the rest of the PULX candidates identified here, the next generation of X-ray telescopes (such as Athena) may be required to determine their nature owing to their lower peak fluxes.

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.

Two recent observations of the nearby galaxy NGC 6946 with NuSTAR, one simultaneous with an XMM-Newton observation, provide an opportunity to examine its population of bright accreting sources from a broadband perspective. We study the three known ultraluminous X-ray sources (ULXs) in the galaxy, and find that ULX-1 and ULX-2 have very steep power-law spectra with $\Gamma=3.6^{+0.4}_{-0.3}$ in both cases. Their properties are consistent with being super-Eddington accreting sources with the majority of their hard emission obscured and down-scattered. ULX-3 (NGC 6946 X-1) is significantly detected by both XMM-Newton and NuSTAR at $L_{\rm X}=(6.5\pm0.1)\times10^{39}$ erg s$^{-1}$, and has a power-law spectrum with $\Gamma=2.51\pm0.05$. We are unable to identify a high-energy break in its spectrum like that found in other ULXs, but the soft spectrum likely hinders our ability to detect one. We also characterise the new source, ULX-4, which is only detected in the joint XMM-Newton and NuSTAR observation, at $L_{\rm X}=(2.27\pm0.07)\times10^{39}$ erg s$^{-1}$, and is absent in a Chandra observation ten days later. It has a very hard cut-off power-law spectrum with $\Gamma=0.7\pm0.1$ and $E_{\rm cut}=11^{+9}_{-4}$ keV. We do not detect pulsations from ULX-4, but its transient nature can be explained either as a neutron star ULX briefly leaving the propeller regime or as a micro-tidal disruption event induced by a stellar-mass compact object.

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.

Intermediate-mass black holes are theoretically predicted but observationally elusive, and evidence for them is often indirect. The nearby face-on spiral galaxy NGC3310 has hosted many supernovae in recent history, and recent Chandra observations have shown a group of strong off-nuclear X-ray sources that are coincident with radio emission seen in archival VLA and MERLIN observations. Their luminosity, spectrum and off-nuclear location make these sources excellent IMBH candidates. To investigate this possibility, we used combined EVN/e-MERLIN observations at both 1.4 and 5 GHz to look for compact radio emission and evidence of jet activity. We detect a compact radio source within one arcsecond of a Chandra source with an estimated mass ${\rm M}_{\rm BH}\sim3\times10^4 {\rm M}_{\odot}$.

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.

We present a broadband X-ray spectral analysis of the M51 system, including the dual active galactic nuclei (AGN) and several off-nuclear point sources. Using a deep observation by NuSTAR, new high-resolution coverage of M51b by Chandra, and the latest X-ray torus models, we measure the intrinsic X-ray luminosities of the AGN in these galaxies. The AGN of M51a is found to be Compton thick, and both AGN have very low accretion rates ($\lambda_{\rm Edd} <10^{-4}$). The latter is surprising considering that the galaxies of M51 are in the process of merging, which is generally predicted to enhance nuclear activity. We find that the covering factor of the obscuring material in M51a is $0.26 \pm 0.03$, consistent with the local AGN obscured fraction at $L_{\rm X}\sim 10^{40}$ erg s$^{-1}$. The substantial obscuring column does not support theories that the torus, presumed responsible for the obscuration, disappears at these low accretion luminosities. However, the obscuration may have resulted from the gas infall driven by the merger rather than the accretion process. We report on several extra-nuclear sources with $L_{\rm X}>10^{39}$ erg s$^{-1}$ and find that a spectral turnover is present below 10 keV in most such sources, in line with recent results on ultraluminous X-ray sources.

We search for transient sources in a sample of ULXs from the 3XMM-DR4 release of the XMM-Newton Serendipitous Source Catalogue in order to find candidate neutron star ULXs alternating between an accreting state and the propeller regime, in which the luminosity drops dramatically. By examining their fluxes and flux upper limits, we identify five ULXs that demonstrate long-term variability of over an order of magnitude. Using Chandra and Swift data to further characterise their light curves, we find that two of these sources are detected only once and could be X-ray binaries in outburst that only briefly reach ULX luminosities. Two others are consistent with being super-Eddington accreting sources with high levels of inter-observation variability. One source, M51 ULX-4, demonstrates apparent bimodal flux behaviour that could indicate the propeller regime. It has a hard X-ray spectrum, but no significant pulsations in its timing data, although with an upper limit of 10% of the signal pulsed at ~1.5 Hz a pulsating ULX cannot be excluded, particularly if the pulsations are transient. By simulating XMM-Newton observations of a population of pulsating ULXs, we predict that there could be approximately 200 other bimodal ULXs that have not been observed sufficiently well by XMM-Newton to be identified as transient.

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.

The luminosity range at and just below the 10^39 erg/s cut-off for defining ultraluminous X-ray sources (ULXs) is a little-explored regime. It none-the-less hosts a large number of X-ray sources, and has great potential for improving our understanding of sources with ~Eddington accretion rates. We select a sample of four sources in this Eddington Threshold regime with good data for further study; these objects possess a variety of soft spectral shapes. We perform X-ray spectral and timing analysis on the XMM-Newton and Chandra data for these objects to gain insight into their accretion mechanisms, and also examine their optical counterparts using HST images. NGC 300 X-1 is a highly luminous and well-known example of the canonical steep power-law accretion state. M51 ULS exhibits a cool blackbody-like spectrum and is consistent with being an ultraluminous supersoft source (ULS), possibly a super-Eddington accreting object viewed at a high inclination through an optically thick outflowing wind. NGC 4395 ULX-1 and NGC 6946 ULX-1 have unusually steep power-law tails, for which we discuss a variety of possible physical mechanisms and links to similar features in Galactic microquasars, and we conclude that these sources are likely intermediate objects between the soft ultraluminous regime of ULXs and classic ULSs.

ULX-7, in the northern spiral arm of M51, demonstrates unusual behaviour for an ultraluminous X-ray source, with a hard X-ray spectrum but very high short-term variability. This suggests that it is not in a typical ultraluminous state. We analyse the source using archival data from XMM-Newton, Chandra and NuSTAR, and by examining optical and radio data from HST and VLA. Our X-ray spectral analysis shows that the source has a hard power-law spectral shape with a photon index Gamma~1.5, which persists despite the source's X-ray luminosity varying by over an order of magnitude. The power spectrum of the source features a break at 6.5^{+0.5}_{-1.1}x10^-3 Hz, from a low-frequency spectral index of alpha_1=-0.1^{+0.5}_{-0.2} to a high-frequency spectral index of alpha_2=0.65^{+0.05}_{-0.14}, making it analogous to the low-frequency break found in the power spectra of low/hard state black holes (BHs). We can take a lower frequency limit for a corresponding high-frequency break to calculate a BH mass upper limit of 1.6x10^3 solar masses. Using the X-ray/radio fundamental plane we calculate another upper limit to the BH mass of 3.5x10^4 solar masses for a BH in the low/hard state. The hard spectrum, high rms variability and mass limits are consistent with ULX-7 being an intermediate-mass BH; however we cannot exclude other interpretations of this source's interesting behaviour, most notably a neutron star with an extreme accretion rate.

We present the current results of an investigation into M51 ULX-7, using archival data from XMM-Newton, Chandra and NuSTAR, and optical and radio data from HST and VLA. The source has a consistently hard power-law X-ray spectrum and high short-term variability. This is unusual variability behaviour for a ULX, as we would expect highly variable ULXs to have soft energy spectra. The power spectrum features a break at ~1e-3 Hz, from low frequency spectral index alpha=0.1 to high frequency spectral index alpha=0.8, analogous to the low frequency break found in power spectra of black holes accreting in the low/hard state. We do not observe a corresponding high frequency break, however taking the white noise level as a frequency lower limit of the break, we can calculate a black hole mass upper limit of 9.12e4 solar masses, assuming that the ULX is in the low/hard state. While there is no radio detection, we find a flux density upper limit of 87 micro-Jy/beam. Using the X-ray/radio fundamental plane, we calculate a black hole mass upper limit of 1.95e5 solar masses. Therefore, this ULX is consistent with being an IMBH accreting in the low/hard state.