Using the NOvA neutrino detectors, a broad search has been performed for any signal coincident with 28 gravitational wave events detected by the LIGO/Virgo Collaboration between September 2015 and July 2019. For all of these events, NOvA is sensitive to possible arrival of neutrinos and cosmic rays of GeV and higher energies. For five (seven) events in the NOvA Far (Near) Detector, timely public alerts from the LIGO/Virgo Collaboration allowed recording of MeV-scale events. No signal candidates were found.

Using two years of data from the NOvA Near Detector at Fermilab, we report a seasonal variation of cosmic ray induced multiple-muon event rates which has an opposite phase to the seasonal variation in the atmospheric temperature. The strength of the seasonal multipl$ increase as a function of the muon multiplicity. However, no significant dependence of the strength of the seasonal variation of the multiple-muon variation is seen as a function of the muon zenith angle, or the spatial or angular separation between the correlated muons.

The Laser Ranging Interferometer (LRI) instrument on the Gravity Recovery and Climate Experiment (GRACE) Follow-On mission has provided the first laser interferometric range measurements between remote spacecraft, separated by approximately 220 km. Autonomous controls that lock the laser frequency to a cavity reference and establish the 5 degree of freedom two-way laser link between remote spacecraft succeeded on the first attempt. Active beam pointing based on differential wavefront sensing compensates spacecraft attitude fluctuations. The LRI has operated continuously without breaks in phase tracking for more than 50 days, and has shown biased range measurements similar to the primary ranging instrument based on microwaves, but with much less noise at a level of $1\,{\rm nm}/\sqrt{\rm Hz}$ at Fourier frequencies above 100 mHz.

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

A sub-set of the brightest ultraluminous X-ray sources (ULXs), with X-ray luminosities well above $10^{40}$ erg s$^{-1}$, typically have energy spectra which can be well described as hard power-laws, and short-term variability in excess of $\sim10\%$. This combination of properties suggests that these ULXs may be some of the best candidates to host intermediate mass black holes (IMBHs), which would be accreting at sub-Eddington rates in the hard state seen in Galactic X-ray binaries. In this work, we present a temporal and spectral analysis of all of the available XMM-Newton data from one such ULX, the previously poorly studied 2XMM J143242.1$-$440939, located in NGC 5643. We report that its high quality EPIC spectra can be better described by a broad, thermal component, such as an advection dominated disc or an optically thick Comptonising corona. In addition, we find a hint of a marginal change in the short-term variability which does not appear to be clearly related to the source unabsorbed luminosity. We discuss the implications of these results, excluding the possibility that the source may be host an IMBH in a low state, and favouring an interpretation in terms of super-Eddington accretion on to a black hole of stellar origin. The properties of NGC 5643 ULX1 allow us to associate this source to the population of the hard/ultraluminous ULX class.

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.

Recent evidence - in particular the hard X-ray spectra obtained by NuSTAR, and the large amplitude hard X-ray variability observed when ultraluminous X-ray sources (ULXs) show soft spectra - reveals that common ULX behaviour is inconsistent with known sub-Eddington accretion modes, as would be expected for an intermediate-mass black hole (IMBH). Instead, it appears that the majority of ULXs are powered by super-Eddington accretion onto stellar-mass black holes. Here, we will review work that delves deeper into ULX spectral-timing behaviour, demonstrating it remains consistent with the expectations of super-Eddington accretion. One critical missing piece from this picture is the direct detection of the massive, radiatively-driven winds expected from ULXs as atomic emission/absorption line features in ULX spectra; we will show it is very likely these have already been detected as residuals in the soft X-ray spectra of ULXs. Finally, we will discuss ULXs that do not appear to conform to the emerging ULX behaviour patterns. In particular we discuss the implications of the identification of a good IMBH candidate as a background QSO; and the confirmation of an IMBH/ULX candidate in the galaxy NGC 2276 via the radio/X-ray fundamental plane.

If ultraluminous X-ray sources (ULXs) are powered by accretion onto stellar remnant black holes, then many must be accreting at super-Eddington rates. It is predicted that such high accretion rates should give rise to massive, radiatively-driven winds. However, observational evidence of a wind, in the form of absorption or emission features, has remained elusive. As such, the reported detection of X-ray spectral residuals in XMM-Newton spectra of NGC 5408 X-1, which could be related to absorption in a wind is potentially very exciting. However, it has previously been assumed by several authors that these features simply originate from background diffuse plasma emission related to star-formation in the ULX's host galaxy. In this work we utilise the spatial resolving power of Chandra to test whether we can rule out this latter interpretation. We demonstrate that the majority of the luminosity in these spectral features is emitted from a highly localised region close to the ULX, and appears point-like even with Chandra. It is therefore highly likely that the spectral features are associated with the ULX itself, and little of the flux in this spectral component originates from spatially extended emission in the host galaxy. This may be consistent with the suggestion of absorption in an optically thin phase of a super-Eddington wind. Alternatively, we could be seeing emission from collisionally ionised material close to the black hole, but critically this would be difficult to reconcile with models where the source inclination largely determines the observed X-ray spectral and timing properties.

The lack of unambiguous detections of atomic features in the X-ray spectra of ultraluminous X-ray sources (ULXs) has proven a hindrance in diagnosing the nature of the accretion flow. The possible association of spectral residuals at soft energies with atomic features seen in absorption and/or emission and potentially broadened by velocity dispersion could therefore hold the key to understanding much about these enigmatic sources. Here we show for the first time that such residuals are seen in several sources and appear extremely similar in shape, implying a common origin. Via simple arguments we assert that emission from extreme colliding winds, absorption in a shell of material associated with the ULX nebula and thermal plasma emission associated with star formation are all highly unlikely to provide an origin. Whilst CCD spectra lack the energy resolution necessary to directly determine the nature of the features (i.e. formed of a complex of narrow lines or intrinsically broad), studying the evolution of the residuals with underlying spectral shape allows for an important, indirect test for their origin. The ULX NGC 1313 X-1 provides the best opportunity to perform such a test due to the dynamic range in spectral hardness provided by archival observations. We show through highly simplified spectral modelling that the strength of the features (in either absorption or emission) appears to anti-correlate with spectral hardness, which would rule out an origin via reflection of a primary continuum and instead supports a picture of atomic transitions in a wind or nearby material associated with such an outflow.

The known members of the class of hyperluminous X-ray sources (HLXs) are few in number, yet they are of great interest as they are regarded as the likeliest intermediate-mass black hole (IMBH) candidates amongst the wider population of ultraluminous X-ray sources (ULXs). Here we report optical photometry and spectroscopy of a HLX candidate associated with the galaxy IC 4320, that reveal it is a background AGN. We discuss the implications of the exclusion of this object from the small number of well-studied HLXs, that appears to accentuate the difference in characteristics between the good IMBH candidate ESO 243-49 HLX-1 and the small handful of other HLXs.

Jet ejection by accreting black holes is a mass invariant mechanism unifying stellar and supermassive black holes (SMBHs) that should also apply for intermediate-mass black holes (IMBHs), which are thought to be the seeds from which SMBHs form. We present the detection of an off-nuclear IMBH of $\sim$5 $\times$ 10$^{4}$ M$_\odot$ located in an unusual spiral arm of the galaxy NGC 2276 based on quasi-simultaneous \textit{Chandra} X-ray observations and European VLBI Network (EVN) radio observations. The IMBH, NGC2276-3c, possesses a 1.8 pc radio jet that is oriented in the same direction as large-scale ($\sim$650 pc) radio lobes and whose emission is consistent with flat to optically thin synchrotron emission between 1.6 GHz and 5 GHz. Its jet kinetic power ($4 \times 10^{40}$ erg s$^{-1}$) is comparable to its radiative output and its jet efficiency ($\geq$ 46\%) is as large as that of SMBHs. A region of $\sim$300 pc along the jet devoid of young stars could provide observational evidence of jet feedback from an IMBH. The discovery confirms that the accretion physics is mass invariant and that seed IMBHs in the early Universe possibly had powerful jets that were an important source of feedback.

We present results obtained from two broadband X-ray observations of the extreme ultraluminous X-ray source (ULX) NGC5907 ULX1, known to have a peak X-ray luminosity of ~5e40 erg/s. These XMM-Newton and NuSTAR observations, separated by only ~4 days, revealed an extreme level of short-term flux variability. In the first epoch, NGC5907 ULX1 was undetected by NuSTAR, and only weakly detected (if at all) with XMM-Newton, while in the second NGC5907 ULX1 was clearly detected at high luminosity by both missions. This implies an increase in flux of ~2 orders of magnitude or more during this ~4 day window. We argue that this is likely due to a rapid rise in the mass accretion rate, rather than to a transition from an extremely obscured to an unobscured state. During the second epoch we observed the broadband 0.3-20.0 keV X-ray luminosity to be (1.55+/-0.06)e40 erg/s, similar to the majority of the archival X-ray observations. The broadband X-ray spectrum obtained from the second epoch is inconsistent with the low/hard accretion state observed in Galactic black hole binaries, but is well modeled with a simple accretion disk model incorporating the effects of photon advection. This strongly suggests that, when bright, NGC5907 ULX1 is a high-Eddington accretor.

Although attempts have been made to constrain the stellar types of optical counterparts to ULXs, the detection of optical variability instead suggests that they may be dominated by reprocessed emission from X-rays which irradiate the outer accretion disc. Here, we report results from a combined X-ray and optical spectral study of a sample of ULXs, which were selected for having broadened disc-like X-ray spectra, and known optical counterparts. We simultaneously fit optical and X-ray data from ULXs with a new spectral model of emission from an irradiated, colour-temperature-corrected accretion disc around a black hole, with a central Comptonising corona. We find that the ULXs require reprocessing fractions of $\sim 10^{-3}$, which is similar to sub-Eddington thermal dominant state BHBs, but less than has been reported for ULXs with soft ultraluminous X-ray spectra. We suggest that the reprocessing fraction may be due to the opposing effects of self-shielding in a geometrically thick super-critical accretion disc, and reflection from far above the central black hole by optically thin material ejected in a natal super-Eddington wind. Then, the higher reprocessing fractions reported for ULXs with wind-dominated X-ray spectra may be due to enhanced scattering onto the outer disc via the stronger wind in these objects. Alternatively, the accretion discs in these ULXs may not be particularly geometrically thick, rather they may be similar in this regard to the thermal dominant state BHBs.

The most extreme ultraluminous X-ray sources (ULXs), with $L_\mathrm{X} > 5 \times 10^{40}$ erg s$^{-1}$, are amongst the best candidates for hosting intermediate-mass black holes (IMBHs) in the haloes of galaxies. Jet radio emission is expected from a sub-Eddington accreting IMBH in the low/hard (radio bright) state. In a search for such IMBH jet radio emission, we have observed with the Very Large Array (VLA) at 5 GHz a sample of seven extreme ULXs whose X-ray properties indicate they are in the hard state. Assuming they remain in this state, the non-detection of radio emission for six of the target sources allows us to constrain their black hole mass to the IMBH regime, thus ruling out a supermassive black hole nature. For the extreme ULX in the galaxy NGC\,2276, we detect extended radio emission formed by two lobes of total flux density 1.43 $\pm$ 0.22 mJy and size $\sim$650 pc. The X-ray counterpart is located between the two lobes, suggesting the presence of a black hole with jet radio emission. The radio luminosity allows us to constrain the black hole mass of this source to the IMBH regime; hence, the extreme ULX in NGC\,2276 could be the first detection of extended jet radio emission from an IMBH. The radio emission could also possibly come from a radio nebula powered by the ULX with a minimum total energy of $5.9 \times 10^{52}$ erg, thus constituting the most powerful and largest ULX radio nebula ever observed.

Although we are nearing a consensus that most ULXs are stellar-mass black holes in a super-Eddington state, little is yet established of the physics of this accretion mode. Here, we use a combined X-ray spectral and timing analysis of a sample of ULXs to investigate this new accretion regime. We suggest a spectral classification scheme that separates ULXs into three classes: a broadened disc class, and two-component hard and soft ultraluminous regimes. At the lowest luminosities the ULX population is dominated by sources with broadened disc spectra, whilst two component spectra are seen at higher luminosities, suggestive of a distinction between ~ Eddington and super-Eddington accretion modes. We find high levels of variability are limited to ULXs with soft ultraluminous spectra, and a few broadened disc sources. Furthermore, the variability is strongest at high energies, suggesting it originates in the harder spectral component. These properties are consistent with current models of super-Eddington emission, where a wind forms a funnel around the central regions of the accretion flow. As the wind provides the soft component this suggests that inclination is the key determinant in the observed X-ray spectra, which is very strongly supported by the variability results if this originates due to clumpy material at the edge of the wind intermittently obscuring our line-of-sight to the central regions of the ULX. The pattern of spectral variability with luminosity in two ULXs that straddle the hard/soft ultraluminous regime boundary is consistent with the wind increasing at higher accretion rates, and thus narrowing the opening angle of the funnel. Hence, this work suggests that most ULXs can be explained as stellar-mass black holes accreting at and above the Eddington limit, with their observed characteristics dominated by two variables: accretion rate and inclination. (abridged)

We present a multi-mission X-ray analysis of a bright (peak observed 0.3-10 keV luminosity of ~ 6x10^{40} erg s^{-1}), but relatively highly absorbed ULX in the edge-on spiral galaxy NGC 5907. The ULX is spectrally hard in X-rays (Gamma ~ 1.2-1.7, when fitted with an absorbed power-law), and has a previously-reported hard spectral break consistent with it being in the ultraluminous accretion state. It is also relatively highly absorbed for a ULX, with a column of ~ 0.4-0.9x10^{22} atom cm^{-2} in addition to the line-of-sight column in our Galaxy. Although its X-ray spectra are well represented by accretion disc models, its variability characteristics argue against this interpretation. The ULX spectra instead appear dominated by a cool, optically-thick Comptonising corona. We discuss how the measured 9 per cent rms variability and a hardening of the spectrum as its flux diminishes might be reconciled with the effects of a very massive, radiatively-driven wind, and subtle changes in the corona respectively. We speculate that the cool disc-like spectral component thought to be produced by the wind in other ULXs may be missing from the observed spectrum due to a combination of a low temperature (~ 0.1 keV), and the high column to the ULX. We find no evidence, other than its extreme X-ray luminosity, for the presence of an intermediate mass black hole (~ 10^2 - 10^4 Msun) in this object. Rather, the observations can be consistently explained by a massive (greater than ~ 20 Msun) stellar remnant black hole in a super-Eddington accretion state.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

We present the results from an X-ray and optical study of a new sample of eight extreme luminosity ultraluminous X-ray source (ULX) candidates, which were selected as the brightest ULXs (with L_X > 5x10^40 erg/s) located within 100 Mpc identified in a cross correlation of the 2XMM-DR1 and RC3 catalogues. These objects are so luminous that they are difficult to describe with current models of super-Eddington accretion onto all but the most massive stellar remnants; hence they are amongst the most plausible candidates to host larger, intermediate-mass black holes (IMBHs). Two objects are luminous enough in at least one observation to be classed as hyperluminous X-ray source (HLX) candidates, including one persistent HLX in an S0 galaxy that (at 3x10^41 erg/s) is the second most luminous HLX yet detected. The remaining seven sources are located in spiral galaxies, and several appear to be closely associated with regions of star formation as is common for many less luminous ULXs. However, the X-ray characteristics of these extreme ULXs appear to diverge from the less luminous objects. They are typically harder, possessing absorbed power-law continuum spectra with photon indexes ~ 1.7, and are potentially more variable on short timescales, with data consistent with ~ 10-20 per cent rms variability on timescales of 0.2-2 ks. These properties appear consistent with the sub-Eddington hard state, which given the observed luminosities of these objects suggests the presence of IMBHs with masses in the range 10^3-10^4 M_Sun. As such, this strengthens the case for these brightest ULXs as good candidates for the eventual conclusive detection of the highly elusive IMBHs. However, we caution that a combination of the highest plausible super-Eddington accretion rates and the largest permitted stellar black hole remnants cannot be ruled out without future, improved observations.

A new, transient ultraluminous X-ray source (ULX) was recently discovered by Chandra in M31 with a luminosity at ~ 5 x 10^39 erg/s. Here we analyse a series of five subsequent XMM-Newton observations. The steady decline in X-ray luminosity over 1.5 months gives an observed e-fold timescale of ~40 days; similar to the decay timescales seen in soft X-ray transients in our own Galaxy. This supports the interpretation of this ULX as a stellar mass black hole in a low-mass X-ray binary (LMXB), accreting at super Eddington rates. This is further supported by the lack of detection of an O/B star in quiescence and the XMM-Newton spectral data being dominated by a disc-like component rather than the power-law expected from a sub-Eddington intermediate-mass black hole. These data give the best sequence of high Eddington fraction spectra ever assembled due to the combination of low absorption column to M31 and well calibrated bandpass down to 0.3 keV of XMM-Newton. The spectra can be roughly described by our best current disc model, BHSPEC, assuming a 10 M_solar black hole with best fit spin of ~0.4. However, the data are better described by a two component model, where the disc emission is significantly affected by advection, and with an additional low temperature Comptonisation component at high energies. This could simply indicate the limitations of current disc models, though the energy-dependent variability also weakly supports a two component interpretation of the data. Irrespective of the detailed interpretation of the spectral properties, these data support the presence of accretion onto a stellar mass black hole in a LMXB accreting in the Eddington-regime. This allows an unambiguous connection of this object, and, by extension, similar low luminosity ULXs, to `standard' X-ray binaries.

The bright ultraluminous X-ray source (ULX), M33 X-8, has been observed several times by XMM-Newton, providing us with a rare opportunity to `flux bin' the spectral data and search for changes in the average X-ray spectrum with flux level. The aggregated X-ray spectra appear unlike standard sub-Eddington accretion state spectra which, alongside the lack of discernible variability at any energy, argues strongly against conventional two-component, sub-Eddington models. Although the lack of variability could be consistent with disc-dominated spectra, sub-Eddington disc models are not sufficiently broad to explain the observed spectra. Fits with a ~ Eddington accretion rate slim disc model are acceptable, but the fits show that the temperature decreases with flux, contrary to expectations, and this is accompanied by the appearance of a harder tail to the spectrum. Applying a suitable two-component model reveals that the disc becomes cooler and less advection dominated as the X-ray flux increases, and this is allied to the emergence of an optically-thick Comptonisation medium. We present a scenario in which this is explained by the onset of a radiatively-driven wind from the innermost regions of the accretion disc, as M33 X-8 exceeds the Eddington limit. Furthermore, we argue that the direct evolution of this spectrum with increasing luminosity (and hence radiation pressure) leads naturally to the two-component spectra seen in more luminous ULXs.

We present highlights from a study of a sample of 10 extreme-luminosity candidate ultraluminous X-ray sources (L_X > 5 x 10^(40) erg sec^(-1)), all at distances < 100 Mpc, identified from a cross-correlation of the RC3 catalogue of galaxies with the 2XMM catalogue. Five of the sample have also been observed by Chandra. Of the 10 sources, seven reside in the disc or arms of spiral galaxies, and the remaining three are close to large elliptical galaxies. Unlike many less luminous ultraluminous X-ray sources, temporal variability is observed on short (ks) and long (year) timescales for most sources in our sample. Long term spectral variability is also evident in some sources. In one case, we use archival Chandra data to demonstrate that a hyperluminous X-ray source candidate identified by XMM-Newton is actually resolved into multiple point sources at high spatial resolution, but note that the other candidates remain unresolved under Chandra's intense scrutiny.