We present results from the remaining sources in our search for near-infrared (NIR) candidate counterparts to ultraluminous X-ray sources (ULXs) within $\simeq$ 10 Mpc. We observed 23 ULXs in 15 galaxies and detected NIR candidate counterparts to six of them. Two of these have an absolute magnitude consistent with a single red supergiant (RSG). Three counterparts are too bright for a RSG and spatially extended, and thus we classify them as stellar clusters. The other candidate is too faint for a RSG. Additionally, we present the results of our NIR spectroscopic follow-up of five sources: four originally classified as RSG and one as a stellar cluster on the basis of previous photometry. The stellar cluster candidate is actually a nebula. Of the four RSGs candidates, one source has a broad H$\alpha$ emission line redshifted by $\sim z = 1$, making it a background AGN. Two other sources show stellar spectra consistent with them being RSGs. The final RSG candidate is too faint to classify, but does not show strong (nebular) emission lines in its spectrum. After our search for NIR counterparts to 113 ULXs, where we detected a candidate counterpart for 38 ULXs, we have spectroscopically confirmed the nature of 12: five sources are nebulae, one source is not classified, one source is an AGN and five are RSGs. These possible five ULX-RSG binary systems would constitute $\simeq$ $(4 \pm 2)\%$ of the observed ULXs, a fraction almost four times larger than what was predicted by binary evolution simulations.

Several dozen optical echelle spectra demonstrate that HR 6819 is a hierarchical triple. A classical Be star is in a wide orbit with an unconstrained period around an inner 40 d binary consisting of a B3 III star and an unseen companion in a circular orbit. The radial-velocity semi-amplitude of 61.3 km/s of the inner star and its minimum (probable) mass of 5.0 Msun (6.3 +- 0.7 Msun) imply a mass of the unseen object of >= 4.2 Msun (>= 5.0 +- 0.4 Msun), that is, a black hole (BH). The spectroscopic time series is stunningly similar to observations of LB-1. A similar triple-star architecture of LB-1 would reduce the mass of the BH in LB-1 from ~70 Msun to a level more typical of Galactic stellar remnant BHs. The BH in HR 6819 probably is the closest known BH to the Sun, and together with LB-1, suggests a population of quiet BHs. Its embedment in a hierarchical triple structure may be of interest for models of merging double BHs or BH + neutron star binaries. Other triple stars with an outer Be star but without BH are identified; through stripping, such systems may become a source of single Be stars.

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 present intermediate resolution spectroscopy of the optical counterpart to the black hole X-ray transient MAXI J1820+070 (=ASASSN-18ey) obtained with the OSIRIS spectrograph on the 10.4-m Gran Telescopio Canarias. The observations were performed with the source close to the quiescent state and before the onset of renewed activity in August 2019. We make use of these data and K-type dwarf templates taken with the same instrumental configuration to measure the projected rotational velocity of the donor star. We find $v_{rot} \sin i = 84 \pm 5$ km s$^{-1}$ ($1\!-\!\sigma$), which implies a donor to black-hole mass ratio $q = {M_2}/{M_1} = 0.072 \pm 0.012$ for the case of a tidally locked and Roche-lobe filling donor star. The derived dynamical masses for the stellar components are $M_1 = (5.95 \pm 0.22)\sin ^{-3}i$ $M_\odot$ and $M_2 = (0.43 \pm 0.08) \sin^{-3}i$ $M_\odot$. The use of $q$, combined with estimates of the accretion disk size at the time of the optical spectroscopy, allows us to revise our previous orbital inclination constraints to $66^{\circ} < i < 81^{\circ}$. These values lead to 95% confidence level limits on the masses of $5.73 <M_1(M_\odot) < 8.34$ and $0.28 < M_2(M_\odot) < 0.77$. Adopting instead the $63 \pm 3^{\circ}$ orientation angle of the radio jet as the binary inclination leads to $M_1 = 8.48^{+0.79}_{-0.72} M_\odot$ and $M_2 = 0.61^{+0.13}_{-0.12} M_\odot$ ($1\!-\!\sigma$).

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.

In 2010 May, an intermediate luminosity optical transient was discovered in the nearby galaxy NGC 300 by a South African amateur astronomer. In the decade since its discovery, multi-wavelength observations of the misnamed ``SN 2010da'' have continually re-shaped our understanding of this high mass X-ray binary system. In this review, we present an overview of the multi-wavelength observations and attempts to understand the 2010 transient event and, later, the re-classification of this system as NGC~300 ULX-1: a red supergiant + neutron star ultraluminous X-ray source.

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.

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

SN2010da/NGC 300 ULX-1 was first detected as a supernova impostor in May 2010 and was recently discovered to be a pulsating ultraluminous X-ray source. In this letter, we present VLT/X-shooter spectra of this source obtained in October 2018, covering the wavelength range 350-2300 nm. The $J$- and $H$-bands clearly show the presence of a red supergiant donor star that is best matched by a MARCS stellar atmosphere with $T_{\rm eff} = 3650 - 3900$ K and $\log(L_{\rm bol}/L_{\odot}) = 4.25\pm0.10$, which yields a stellar radius $R = 310 \pm 70 R_{\odot}$. To fit the full spectrum, two additional components are required: a blue excess that can be fitted either by a hot blackbody (T $\gtrsim 20,000$ K) or a power law (spectral index $\alpha \approx 4$) and is likely due to X-ray emission reprocessed in the outer accretion disk or the donor star; and a red excess that is well fitted by a blackbody with a temperature of $\sim 1100$ K, and is likely due to warm dust in the vicinity of SN2010da. The presence of a red supergiant in this system implies an orbital period of at least 0.8-2.1 years, assuming Roche lobe overflow. Given the large donor-to-compact object mass ratio, orbital modulations of the radial velocity of the red supergiant are likely undetectable. However, the radial velocity amplitude of the neutron star is large enough (up to 40-60 km s$^{-1}$) to potentially be measured in the future, unless the system is viewed at a very unfavorable inclination.

We present time-resolved 10.4-m GTC and 4.2-m WHT intermediate resolution spectroscopy of the X-ray transient MAXI J1820+070 (=ASASSN-18ey) obtained during its decline to the quiescent state. Cross-correlation of the 21 individual spectra against late-type templates reveals a sinusoidal velocity modulation with a period of 0.68549 +/- 0.00001 d and semi-amplitude of 417.7 +/- 3.9 km/s. We derive a mass function f(M) = 5.18 +/- 0.15 Msun, dynamically confirming the black hole nature of the compact object. Our analysis of the stellar absorption features supports a K3-5 spectral classification for the donor star, which contributes 20% of the total flux at 5200-6800 Angs. The photometric 0.703 +/- 0.003 d periodicity observed during outburst is 2.6% longer than the orbital period supporting the presence of a superhump modulation in the outburst lightcurves. In line with this interpretation, we constrain the binary mass ratio to be q=0.12. In addition, we observe a sharp increase in the Halpha emission line equivalent width during inferior conjunction of the donor star that we interpret as a grazing eclipse of the accretion disc and allows us to constrain the binary inclination to > 69 deg. On the other hand, the absence of X-ray eclipses during outburst imply i < 77 deg. These inclination limits, together with our dynamical solution, lead to a black hole mass in the range 7-8 Msun. We also measure a systemic velocity = -21.6 +/- 2.3 km/s which, combined with the Gaia DR2 proper motion and parallax, implies a large peculiar velocity of 100 km/s.

We present Very Large Telescope/X-shooter and Chandra X-ray observatory/ACIS observations of the ULX [SST2011] J110545.62+000016.2 in the galaxy NGC 3521. The source identified as a candidate near-infrared counterpart to the ULX in our previous study shows an emission line spectrum of numerous recombination and forbidden lines in the visible and near-infrared spectral regime. The emission from the candidate counterpart is spatially extended ($\sim$ 34 pc) and appears to be connected with an adjacent H II region, located $\sim$ 138 pc to the NE. The measured velocities of the emission lines confirm that both the candidate counterpart and H II region reside in NGC 3521. The intensity ratios of the emission lines from the ULX counterpart show that the line emission originates from the combined effect of shock and photoionisation of low metallicity (12 + log (O/H) = 8.19 $\pm$ 0.11) gas. Unfortunately, there is no identifiable spectral signature directly related to the photosphere of the mass-donor star in our spectrum. From the archival Chandra data, we derive the X-ray luminosity of the source in the 0.3-7 keV range to be (1.9 $\pm$ 0.8) $\times$ 10$^{40}$ erg cm$^{-2}$ s$^{-1}$, almost a factor of four higher than what is previously reported.

We report optical and infrared observations of the X-ray source CXOGBS J174623.5-310550. This Galactic object was identified as a potential quiescent low-mass X-ray binary accreting from an M-type donor on the basis of optical spectroscopy and the broad Halpha emission line. The analysis of X-shooter spectroscopy covering 3 consecutive nights supports an M2/3-type spectral classification. Neither radial velocity variations nor rotational broadening is detected in the photospheric lines. No periodic variability is found in I- and r'-band light curves. We derive r' = 20.8, I = 19.2 and Ks = 16.6 for the optical and infrared counterparts with the M-type star contributing 90% to the I-band light. We estimate its distance to be 1.3-1.8 kpc. The lack of radial velocity variations implies that the M-type star is not the donor star in the X-ray binary. This could be an interloper or the outer body in a hierarchical triple. We constrain the accreting binary to be a < 2.2 hr orbital period eclipsing cataclysmic variable or a low-mass X-ray binary lying in the foreground of the Galactic Bulge.

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.

We present a mid-infrared (IR) sample study of nearby ultraluminous X-ray sources (ULXs) using multi-epoch observations with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. Spitzer/IRAC observations taken after 2014 were obtained as part of the Spitzer Infrared Intensive Transients Survey (SPIRITS). Our sample includes 96 ULXs located within 10 Mpc. Of the 96~ULXs, 12 have candidate counterparts consistent with absolute mid-IR magnitudes of supergiants, and 16 counterparts exceeded the mid-IR brightness of single supergiants and are thus more consistent with star clusters or non-ULX background active galactic nuclei (AGN). The supergiant candidate counterparts exhibit a bi-modal color distribution in a Spitzer/IRAC color-magnitude diagram, where "red" and "blue" ULXs fall in IRAC colors $[3.6] - [4.5]\sim0.7$ and $[3.6] - [4.5]\sim0.0$, respectively. The mid-IR colors and absolute magnitudes of 4 "red" and 5 "blue" ULXs are consistent with that of supergiant B[e] (sgB[e]) and red supergiant (RSG) stars, respectively. While "blue", RSG-like mid-IR ULX counterparts likely host RSG mass donors, we propose the "red" counterparts are ULXs exhibiting the "B[e] phenomenon'' rather than hosts of sgB[e] mass donors. We show that the mid-IR excess from the "red" ULXs is likely due to thermal emission from circumstellar or circumbinary dust. Using dust as a probe for total mass, we estimate mass-loss rates of $\dot{M}\sim1\times10^{-4}$ M$_\odot$ yr$^{-1}$ in dust-forming outflows of red ULXs. Based on the transient mid-IR behavior and its relatively flat spectral index, $\alpha=-0.19\pm0.1$, we suggest that the mid-IR emission from Holmberg IX X-1 originates from a variable jet.

What are the most important conditions and processes governing the growth of stellar-origin compact objects? The identification of compact object type as either black hole (BH) or neutron star (NS) is fundamental to understanding their formation and evolution. To date, time-domain determination of compact object type remains a relatively untapped tool. Measurement of orbital periods, pulsations, and bursts will lead to a revolution in the study of the demographics of NS and BH populations, linking source phenomena to accretion and galaxy parameters (e.g., star formation, metallicity). To perform these measurements over sufficient parameter space, a combination of a wide-field (>5000 deg^2) transient X-ray monitor over a dynamic energy range (~1-100 keV) and an X-ray telescope for deep surveys with <5 arcsec PSF half-energy width (HEW) angular resolution are required. Synergy with multiwavelength data for characterizing the underlying stellar population will transform our understanding of the time domain properties of transient sources, helping to explain details of supernova explosions and gravitational wave event rates.

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.

We report on our search for the optical counterparts of two ultraluminous X-ray pulsars with known orbital periods, M82 X-2 and NGC 5907 X-1, in new and archival HST observations, in an effort to characterize the donor stars in these systems. We detect five near-infrared sources consistent with the position of M82 X-2 that are too bright to be single stars. We also detect seven sources in the WFC3/UVIS F336W image whose photometry matches that of 10-15 M$_\odot$ stars turning off the main sequence. Such stars have densities consistent with the properties of the donor star of M82 X-2 as inferred from X-ray timing analysis, although it is also possible that the donor is a lower mass star below our detection limit or that there is a significant contribution from the accretion disc to the optical emission. We detect three candidate counterparts to NGC 5907 X-1 in the near-infrared. All of these are too bright to be the donor star of the ULX, which based on its orbital period is a red giant. The high background at the location of NGC 5907 X-1 precludes us from detecting this expected donor star. The recently discovered NGC 5907 ULX-2 also falls within the field of view of the near-infrared imaging; we detect four sources in the error circle, with photometry that matches AGB stars. The star suggested to be the counterpart of NGC 5907 ULX-2 by Pintore et al. (2018) falls outside our 2-$\sigma$ error circle.

Accreting supermassive black holes (SMBHs) can exhibit variable emission across the electromagnetic spectrum and over a broad range of timescales. The variability of active galactic nuclei (AGNs) in the ultraviolet and optical is usually at the few tens of per cent level over timescales of hours to weeks. Recently, rare, more dramatic changes to the emission from accreting SMBHs have been observed, including tidal disruption events, 'changing look' AGNs and other extreme variability objects. The physics behind the 're-ignition', enhancement and 'shut-down' of accretion onto SMBHs is not entirely understood. Here we present a rapid increase in ultraviolet-optical emission in the centre of a nearby galaxy, marking the onset of sudden increased accretion onto a SMBH. The optical spectrum of this flare, dubbed AT 2017bgt, exhibits a mix of emission features. Some are typical of luminous, unobscured AGNs, but others are likely driven by Bowen fluorescence - robustly linked here with high-velocity gas in the vicinity of the accreting SMBH. The spectral features and increased ultraviolet flux show little evolution over a period of at least 14 months. This disfavours the tidal disruption of a star as their origin, and instead suggests a longer-term event of intensified accretion. Together with two other recently reported events with similar properties, we define a new class of SMBH-related flares. This has important implications for the classification of different types of enhanced accretion onto SMBHs.

Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwide campaign to observe AT2018cow (ATLAS18qqn), the first fast-luminous optical transient to be found in real time at low redshift. Our first spectra (<2 days after discovery) are entirely featureless. A very broad absorption feature suggestive of near-relativistic velocities develops between 3-8 days, then disappears. Broad emission features of H and He develop after >10 days. The spectrum remains extremely hot throughout its evolution, and the photospheric radius contracts with time (receding below R<10^14 cm after 1 month). This behaviour does not match that of any known supernova, although a relativistic jet within a fallback supernova could explain some of the observed features. Alternatively, the transient could originate from the disruption of a star by an intermediate-mass black hole, although this would require long-lasting emission of highly super-Eddington thermal radiation. In either case, AT2018cow suggests that the population of fast luminous transients represents a new class of astrophysical event. Intensive follow-up of this event in its late phases, and of any future events found at comparable distance, will be essential to better constrain their origins.

We present near-infrared and optical imaging of fifteen candidate black hole X-ray binaries. In addition to quiescent observations for all sources, we also observed two of these sources (IGR J17451-3022 and XTE J1818-245) in outburst. We detect the quiescent counterpart for twelve out of fifteen sources, and for the remaining three we report limiting magnitudes. The magnitudes of the detected counterparts range between $K_s$ = 17.59 and $K_s$ = 22.29 mag. We provide (limits on) the absolute magnitudes and finding charts of all sources. Of these twelve detections in quiescence, seven represent the first quiescent reported values (for MAXI J1543-564, XTE J1726-476, IGR J17451-3022, XTE J1818-245, MAXI J1828-249, MAXI J1836-194, Swift J1910.2-0546) and two detections show fainter counterparts to XTE J1752-223 and XTE J2012+381 than previously reported. We used theoretical arguments and observed trends, for instance between the outburst and quiescent X-ray luminosity and orbital period $P_{orb}$ to derive an expected trend between $\Delta K_s$ and $P_{orb}$ of $\Delta K_s \propto \log P_{orb}^{0.565}$. Comparing this to observations we find a different behaviour. We discuss possible explanations for this result.

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

We present Chandra and optical observations of a candidate dual AGN discovered serendipitously while searching for recoiling black holes via a cross-correlation between the serendipitous XMM source catalog (2XMMi) and SDSS-DR7 galaxies with a separation no larger than ten times the sum of their Petrosian radii. The system has a stellar mass ratio M$_{1}$/M$_{2}\approx 0.7$. One of the galaxies (Source 1) shows clear evidence for AGN activity in the form of hard X-ray emission and optical emission-line diagnostics typical of AGN ionisation. The nucleus of the other galaxy (Source 2) has a soft X-ray spectrum, bluer colours, and optical emission line ratios dominated by stellar photoionisation with a "composite" signature, which might indicate the presence of a weak AGN. When plotted on a diagram with X-ray luminosity vs [OIII] luminosity both nuclei fall within the locus defined by local Seyfert galaxies. From the optical spectrum we estimate the electron densities finding n$_{1} < 27$ e$^{-}$ cm$^{-3}$ and n$_{2} \approx 200$ e$^{-}$ cm$^{-3}$. From a 2D decomposition of the surface brightness distribution we infer that both galaxies host rotationally supported bulges (Sersic index $< 1$). While the active nature of Source 1 can be established with confidence, whether the nucleus of Source 2 is active remains a matter of debate. Evidence that a faint AGN might reside in its nucleus is, however, tantalising.

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

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.

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.

With the first direct detection of merging black holes in 2015, the era of gravitational wave (GW) astrophysics began. A complete picture of compact object mergers, however, requires the detection of an electromagnetic (EM) counterpart. We report ultraviolet (UV) and X-ray observations by Swift and the Nuclear Spectroscopic Telescope ARray (NuSTAR) of the EM counterpart of the binary neutron star merger GW170817. The bright, rapidly fading ultraviolet emission indicates a high mass ($\approx0.03$ solar masses) wind-driven outflow with moderate electron fraction ($Y_{e}\approx0.27$). Combined with the X-ray limits, we favor an observer viewing angle of $\approx 30^{\circ}$ away from the orbital rotation axis, which avoids both obscuration from the heaviest elements in the orbital plane and a direct view of any ultra-relativistic, highly collimated ejecta (a gamma-ray burst afterglow).

We obtained 16 VLT/X-shooter observations of GX 339-4 in quiescence in the period May - September 2016 and detected absorption lines from the donor star in its NIR spectrum. This allows us to measure the radial velocity curve and projected rotational velocity of the donor for the first time. We confirm the 1.76 day orbital period and we find that $K_2$ = $219 \pm 3$ km s$^{-1}$, $\gamma = 26 \pm 2$ km s$^{-1}$ and $v \sin i = 64 \pm 8$ km s$^{-1}$. From these values we compute a mass function $f(M) =1.91 \pm 0.08~M_{\odot}$, a factor $\sim 3$ lower than previously reported, and a mass ratio $q = 0.18 \pm 0.05$. We confirm the donor is a K-type star and estimate that it contributes $\sim 45-50\%$ of the light in the $J$- and H-band. We constrain the binary inclination to $37^\circ < i < 78^\circ$ and the black hole mass to $2.3~M_{\odot} < M_\mathrm{BH} < 9.5~M_{\odot}$. GX 339-4 may therefore be the first black hole to fall in the 'mass-gap' of $2-5~M_{\odot}$.

We present the results of our continued systematic search for near-infrared (NIR) candidate counterparts to ultraluminous X-ray sources (ULXs) within 10 Mpc. We observed 42 ULXs in 24 nearby galaxies and detected NIR candidate counterparts to 15 ULXs. Fourteen of these ULXs appear to have a single candidate counterpart in our images and the remaining ULX has 2 candidate counterparts. Seven ULXs have candidate counterparts with absolute magnitudes in the range between -9.26 and -11.18 mag, consistent with them being red supergiants (RSGs). The other eight ULXs have candidate counterparts with absolute magnitudes too bright to be a single stellar source. Some of these NIR sources show extended morphology or colours expected for Active Galactic Nuclei (AGN), strongly suggesting that they are likely stellar clusters or background galaxies. The red supergiant candidate counterparts form a valuable sample for follow-up spectroscopic observations to confirm their nature, with the ultimate goal of directly measuring the mass of the compact accretor that powers the ULX using binary Doppler shifts.

We present mid-infrared (IR) light curves of the Ultraluminous X-ray Source (ULX) Holmberg II X-1 from observations taken between 2014 January 13 and 2017 January 5 with the \textit{Spitzer Space Telescope} at 3.6 and 4.5 $\mu$m in the \textit{Spitzer} Infrared Intensive Transients Survey (SPIRITS). The mid-IR light curves, which reveal the first detection of mid-IR variability from a ULX, is determined to arise primarily from dust emission rather than from a jet or an accretion disk outflow. We derived the evolution of the dust temperature ($T_\mathrm{d}\sim600 - 800$ K), IR luminosity ($L_\mathrm{IR}\sim3\times10^4$ $\mathrm{L}_\odot$), mass ($M_\mathrm{d}\sim1-3\times10^{-6}$ $\mathrm{M}_\odot$), and equilibrium temperature radius ($R_\mathrm{eq}\sim10-20$ AU). A comparison of X-1 with a sample spectroscopically identified massive stars in the Large Magellanic Cloud on a mid-IR color-magnitude diagram suggests that the mass donor in X-1 is a supergiant (sg) B[e]-star. The sgB[e]-interpretation is consistent with the derived dust properties and the presence of the [Fe II] ($\lambda=1.644$ $\mu$m) emission line revealed from previous near-IR studies of X-1. We attribute the mid-IR variability of X-1 to increased heating of dust located in a circumbinary torus. It is unclear what physical processes are responsible for the increased dust heating; however, it does not appear to be associated with the X-ray flux from the ULX given the constant X-ray luminosities provided by serendipitous, near-contemporaneous X-ray observations around the first mid-IR variability event in 2014. Our results highlight the importance of mid-IR observations of luminous X-ray sources traditionally studied at X-ray and radio wavelengths.

We report the detection of a $78.1\pm0.5$ day period in the X-ray lightcurve of the extreme ultraluminous X-ray source NGC 5907 ULX1 ($L_{\rm{X,peak}}\sim5\times10^{40}$ erg s$^{-1}$), discovered during an extensive monitoring program with Swift. These periodic variations are strong, with the observed flux changing by a factor of $\sim$3-4 between the peaks and the troughs of the cycle; our simulations suggest that the observed periodicity is detected comfortably in excess of 3$\sigma$ significance. We discuss possible origins for this X-ray period, but conclude that at the current time we cannot robustly distinguish between orbital and super-orbital variations.

We present H-band spectra of the candidate counterparts of five ULXs (two in NGC 925, two in NGC 4136, and Holmberg II X-1) obtained with Keck/MOSFIRE. The candidate counterparts of two ULXs (J022721+333500 in NGC 925 and J120922+295559 in NGC 4136) have spectra consistent with (M-type) red supergiants (RSGs). We obtained two epochs of spectroscopy of the candidate counterpart to J022721+333500, separated by 10 months, but discovered no radial velocity variations with a 2-$\sigma$ upper limit of 40 km/s. If the RSG is the donor star of the ULX, the most likely options are that either the system is seen at low inclination (< 40$^\circ$), or the black hole mass is less than 100 M$_\odot$, unless the orbital period is longer than 6 years, in which case the obtained limit is not constraining. The spectrum of the counterpart to J120922+295559 shows emission lines on top of its stellar spectrum, and the remaining three counterparts do not show absorption lines that can be associated with the atmosphere of a star; their spectra are instead dominated by emission lines. Those counterparts with RSG spectra may be used in the future to search for radial velocity variations, and, if those are present, determine dynamical constraints on the mass of the accretor.

We report on six Chandra and one HST/WFC3 observation of CXO J122518.6+144545, discovered by Jonker et al. (2010) as a candidate hyperluminous X-ray source (HLX), X-ray bright supernova or recoiling supermassive black hole at $L_X = 2.2 \times 10^{41}$ erg/s (if associated with the galaxy at 182 Mpc). We detect a new outburst of the source in a Chandra image obtained on Nov 20, 2014 and show that the X-ray count rate varies by a factor $> 60$. New HST/WFC3 observations obtained in 2014 show that the optical counterpart is still visible at $g' = 27.1 \pm 0.1$, $1 \pm 0.1$ magnitude fainter than in the discovery HST/ACS observation from 2003. This optical variability strongly suggests that the optical and X-ray source are related. Furthermore, these properties strongly favour an HLX nature of the source over the alternative scenarios. We therefore conclude that CXO J122518.6+144545 is most likely an outbursting HLX. It is only the second such object to be discovered, after HLX-1 in ESO 243-49. Its high X-ray luminosity makes it a strong candidate to host an intermediate mass black hole.

We present two epochs of near-infrared spectroscopy of the candidate red supergiant counterpart to RX~J004722.4-252051, a ULX in NGC 253. We measure radial velocities of the object and its approximate spectral type by cross-correlating our spectra with those of known red supergiants. Our VLT/X-shooter spectrum is best matched by that of early M-type supergiants, confirming the red supergiant nature of the candidate counterpart. The radial velocity of the spectrum, taken on 2014, August 23, is $417 \pm 4$ km/s. This is consistent with the radial velocity measured in our spectrum taken with Magellan/MMIRS on 2013, June 28, of $410 \pm 70$ km/s, although the large error on the latter implies that a radial velocity shift expected for a black hole of tens of $M_\odot$ can easily be hidden. Using nebular emission lines we find that the radial velocity due to the rotation of NGC 253 is 351 $\pm$ 4 km/s at the position of the ULX. Thus the radial velocity of the counterpart confirms that the source is located in NGC 253, but also shows an offset with respect to the local bulk motion of the galaxy of 66 $\pm$ 6 km/s. We argue that the most likely origin for this displacement lies either in a SN kick, requiring a system containing a $\gtrsim 50$ $M_\odot$ black hole, and/or in orbital radial velocity variations in the ULX binary system, requiring a $\gtrsim 100$ $M_\odot$ black hole. We therefore conclude that RX~J004722.4-252051 is a strong candidate for a ULX containing a massive stellar black hole.

In this paper we present the results of the first systematic search for counterparts to nearby ultraluminous X-ray sources (ULXs) in the near-infrared (NIR). We observed 62 ULXs in 37 galaxies within 10 Mpc and discovered 17 candidate NIR counterparts. The detection of 17 out of 62 ULX candidates points to intrinsic differences between systems that show and those that do not show infrared emission. For six counterparts we conclude from the absolute magnitudes and - - in some cases - additional information such as morphology and previously reported photometric or spectroscopic observations, that they are likely background active galactic nuclei or ULXs residing in star clusters. Eleven counterparts have absolute magnitudes consistent with them being single red supergiant stars. Alternatively, these systems may have larger accretion discs that emit more NIR light than the systems that we do not detect. Other scenarios such as emission from a surrounding nebula or from a compact radio jet are also possible, although for Holmberg II X-1 the NIR luminosity far exceeds the expected jet contribution. The eleven possible red supergiant counterparts are excellent candidates for spectroscopic follow-up observations. This may enable us to measure the mass function in these systems if they are indeed red supergiant donor stars where we can observe absorption lines.

We present the discovery of a new type of explosive X-ray flash in Chandra images of the old elliptical galaxy M86. This unique event is characterised by the peak luminosity of 6x10^42 erg/s for the distance of M86, the presence of precursor events, the timescale between the precursors and the main event (~4,000 s), the absence of detectable hard X-ray and gamma-ray emission, the total duration of the event and the detection of a faint associated optical signal. The transient is located close to M86 in the Virgo cluster at the location where gas and stars are seen protruding from the galaxy probably due to an ongoing wet minor merger. We discuss the possible mechanisms for the transient and we conclude that the X-ray flash could have been caused by the disruption of a compact white dwarf star by a ~10^4 Msun black hole. Alternative scenarios such that of a foreground neutron star accreting an asteroid or the detection of an off-axis (short) gamma-ray burst cannot be excluded at present.

We obtained VLT/FORS2 spectra of the optical counterparts of four high-luminosity (L_X >= 10^40 erg/s) ULX candidates from the catalog of Walton (2011). We first determined accurate positions for the X-ray sources from archival Chandra observations and identified counterparts in archival optical observations that are sufficiently bright for spectroscopy with an 8 meter telescope. From the spectra we determine the redshifts to the optical counterparts and emission line ratios. One of the candidate ULXs, in the spiral galaxy ESO 306-003, appears to be a bona fide ULX in an HII region. The other three sources, near the elliptical galaxies NGC 533 and NGC 741 and in the ring galaxy AM 0644-741, turn out to be background AGN with redshifts of 1.85, 0.88 or 1.75 and 1.40 respectively. Our findings confirm the trend of a high probability of finding background AGN for systems with a ratio of log(F_X/F_opt) in the range of -1 to 1.

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.

In this paper we report accurate Chandra positions for two ultraluminous X-ray sources: NGC 7319-X4 at Right Ascension (RA) = 339.02917(2) deg, Declination (Dec) = 33.97476(2) deg and NGC 5474-X1 at RA = 211.24859(3) deg, Dec = 53.63584(3) deg. We perform bore-sight corrections on the Chandra X-ray Satellite observations of these sources to get to these accurate positions of the X-ray sources and match these positions with archival optical data from the Wide Field and Planetary Camera 2 on board the Hubble Space Telescope. We do not find the optical counterparts: the limiting absolute magnitudes of the observations in the WFPC2 standard magnitude system are B = -7.9, V = -8.7 and I = -9.3 for NGC 7319-X4 and U = -6.4 for NGC 5474-X1. We report on the X-ray spectral properties and we find evidence for X-ray variability in NGC 5474-X1. Finally, we briefly discuss several options for the nature of these ULXs.

In this Paper we report the discovery of CXO J122518.6+144545; a peculiar X-ray source with a position 3.6+-0.2",off-nuclear from an SDSS DR7 z=0.0447 galaxy. The 3.6" offset corresponds to 3.2 kpc at the distance of the galaxy. The 0.3-8 keV X-ray flux of this source is 5x10^-14 erg cm^-2 s^-1 and its 0.3-8 keV luminosity is 2.2x10^41 erg/s (2.7x10^41 erg/s; 0.5-10 keV) assuming the source belongs to the associated galaxy. We find a candidate optical counterpart in archival HST/ACS g'-band observations of the field containing the galaxy obtained on June 16, 2003. The observed magnitude of g'=26.4+-0.1 corresponds to an absolute magnitude of -10.1. We discuss the possible nature of the X-ray source and its associated candidate optical counterpart and conclude that the source is either a very blue type IIn supernova, a ULX with a very bright optical counterpart or a recoiling super-massive black hole.