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.

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.

Optical and radio observations of shock-ionised bubble nebulae surrounding ultraluminous X-ray sources (ULXs) suggest that they are powered by jets or super-critical outflows presumably launched from the ULX accretion disc. Recent simulations of these systems have shown that the shocked wind can emit thermal X-rays with estimated luminosities $\sim 10^{36}$ erg s$^{-1}$. In this work, we investigated whether it is possible to detect and spatially resolve the X-ray emission from these systems using archival Chandra observations of the ULX Holmberg IX X-1. This source is an ideal target to study for two reasons: it is surrounded by an optical bubble nebula with a large spatial extent ($\sim$ 400 pc) that can easily be resolved with Chandra. Further, it has a hard X-ray continuum that is easily distinguishable from the expected soft thermal emission from the nebula. However, a spectral and photometric analysis on stacked Chandra observations of the source reveals that there is no strong evidence for an X-ray bubble associated with it, to a limiting luminosity of $\sim 2 \times 10^{36}$ erg s$^{-1}$. The detection of such X-ray nebulae may be possible with future X-ray missions such as Athena, which would provide useful constraints on the kinematics of the outflow. Finally, our observations also emphasise that the nebular emission does not contribute significantly to the residuals in the X-ray spectrum of the source, which are more likely to be linked to processes localised to the ULX.

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.

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.