We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

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.

IGR J17591$-$2342 is a new accreting millisecond X-ray pulsar (AMXP) that was recently discovered in outburst in 2018. Early observations revealed that the source's radio emission is brighter than that of any other known neutron star low-mass X-ray binary (NS-LMXB) at comparable X-ray luminosity, and assuming its likely $\gtrsim 6$ kpc distance. It is comparably radio bright to black hole LMXBs at similar X-ray luminosities. In this work, we present the results of our extensive radio and X-ray monitoring campaign of the 2018 outburst of IGR J17591$-$2342. In total we collected 10 quasi-simultaneous radio (VLA, ATCA) and X-ray (Swift-XRT) observations, which make IGR J17591$-$2342 one of the best-sampled NS-LMXBs. We use these to fit a power-law correlation index $\beta = 0.37^{+0.42}_{-0.40}$ between observed radio and X-ray luminosities ( $L_\mathrm{R}\propto L_\mathrm{X}^{\beta}$). However, our monitoring revealed a large scatter in IGR J17591$-$2342's radio luminosity (at a similar X-ray luminosity, $L_\mathrm{X} \sim 10^{36}$ erg s$^{-1}$, and spectral state), with $L_\mathrm{R} \sim 4 \times 10^{29}$ erg s$^{-1}$ during the first three reported observations, and up to a factor of 4 lower $L_\mathrm{R}$ during later radio observations. Nonetheless, the average radio luminosity of IGR J17591$-$2342 is still one of the highest among NS-LMXBs, and we discuss possible reasons for the wide range of radio luminosities observed in such systems during outburst. We found no evidence for radio pulsations from IGR J17591$-$2342 in our Green Bank Telescope observations performed shortly after the source returned to quiescence. Nonetheless, we cannot rule out that IGR J17591$-$2342 becomes a radio millisecond pulsar during quiescence.

We present results from a search for a radio transient associated with the LIGO/Virgo source S190814bv, a likely neutron star-black hole (NSBH) merger, with the Australian Square Kilometre Array Pathfinder. We imaged a $30\,{\rm deg}^2$ field at $\Delta T$=2, 9 and 33 days post-merger at a frequency of 944\,MHz, comparing them to reference images from the Rapid ASKAP Continuum Survey observed 110 days prior to the event. Each epoch of our observations covers $89\%$ of the LIGO/Virgo localisation region. We conducted an untargeted search for radio transients in this field, resulting in 21 candidates. For one of these, \object[AT2019osy]{AT2019osy}, we performed multi-wavelength follow-up and ultimately ruled out the association with S190814bv. All other candidates are likely unrelated variables, but we cannot conclusively rule them out. We discuss our results in the context of model predictions for radio emission from neutron star-black hole mergers and place constrains on the circum-merger density and inclination angle of the merger. This survey is simultaneously the first large-scale radio follow-up of an NSBH merger, and the most sensitive widefield radio transients search to-date.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

Transitional millisecond pulsars (tMSPs), which are systems that harbor a pulsar in the throes of the recycling process, have emerged as a new source class since the discovery of the first such system a decade ago. These systems switch between accretion-powered low-mass X-ray binary (LMXB) and rotation-powered radio millisecond pulsar (RMSP) states, and provide exciting avenues to understand the physical processes that spin-up neutron stars to millisecond periods. During the last decade, three tMSPs, as well as a candidate source, have been extensively probed using systematic, multi-wavelength campaigns. Here we review the observational highlights from these campaigns and our general understanding of tMSPs.

We present coordinated Chandra X-ray Observatory and Karl G. Jansky Very Large Array observations of the transitional millisecond pulsar PSR J1023+0038 in its low-luminosity accreting state. The unprecedented five hours of strictly simultaneous X-ray and radio continuum coverage for the first time unambiguously show a highly reproducible, anti-correlated variability pattern. The characteristic switches from the X-ray high mode into a low mode are always accompanied by a radio brightening with duration that closely matches the X-ray low mode interval. This behavior cannot be explained by a canonical inflow/outflow accretion model where the radiated emission and the jet luminosity are powered by, and positively correlated with, the available accretion energy. We interpret this phenomenology as alternating episodes of low-level accretion onto the neutron star during the X-ray high mode that are interrupted by rapid ejections of plasma by the active rotation-powered pulsar, possibly initiated by a reconfiguration of the pulsar magnetosphere, that cause a transition to a less luminous X-ray mode. The observed anti-correlation between radio and X-ray luminosity has an additional consequence: transitional MSPs can make excursions into a region of the radio/X-ray luminosity plane previously thought to be occupied solely by black hole X-ray binary sources. This complicates the use of this luminosity relation to identify candidate black holes, suggesting the need for additional discriminants when attempting to establish the true nature of the accretor.

We analyzed broad-band X-ray and radio data of the magnetar SGR J1935+2154 taken in the aftermath of its 2014, 2015, and 2016 outbursts. The source soft X-ray spectrum <10 keV is well described with a BB+PL or 2BB model during all three outbursts. NuSTAR observations revealed a hard X-ray tail, $\Gamma=0.9$, extending up to 79 keV, with flux larger than the one detected <10 keV. Imaging analysis of Chandra data did not reveal small-scale extended emission around the source. Following the outbursts, the total 0.5-10 keV flux from SGR J1935+2154 increased in concordance to its bursting activity, with the flux at activation onset increasing by a factor of $\sim7$ following its strongest June 2016 outburst. A Swift/XRT observation taken 1.5 days prior to the onset of this outburst showed a flux level consistent with quiescence. We show that the flux increase is due to the PL or hot BB component, which increased by a factor of $25$ compared to quiescence, while the cold BB component $kT=0.47$ keV remained more or less constant. The 2014 and 2015 outbursts decayed quasi-exponentially with time-scales of $\sim40$ days, while the stronger May and June 2016 outbursts showed a quick short-term decay with time-scales of $\sim4$ days. Our Arecibo radio observations set the deepest limits on the radio emission from a magnetar, with a maximum flux density limit of 14 $\mu$Jy for the 4.6 GHz observations and 7 $\mu$Jy for the 1.4 GHz observations. We discuss these results in the framework of the current magnetar theoretical models.

The accreting millisecond X-ray pulsar (AMXP) SAX J1808.4-3658, shows a peculiar orbital evolution that proceeds at a much faster pace than predicted by conservative binary evolution models. It is important to identify the underlying mechanism responsible for this behavior because it can help to understand how this system evolves. It has also been suggested that, when in quiescence, SAX J1808.4-3658 turns on as a radio pulsar, a circumstance that might provide a link between AMXPs and black-widow radio pulsars. In this work we report the results of a deep radio pulsation search at 2 GHz using the Green Bank Telescope in August 2014 and an X-ray monitoring of the 2015 outburst with Chandra, Swift, and INTEGRAL. In particular, we present the X-ray timing analysis of a 30-ks Chandra observation executed during the 2015 outburst. We detect no radio pulsations, and place the strongest limit to date on the pulsed radio flux density of any AMXP. We also find that the orbit of SAX J1808.4-3658 continues evolving at a fast pace and we compare it to the bhevior of other accreting and non-accreting binaries. We discuss two scenarios: either the neutron star has a large moment of inertia (I>1.7x10^45 g cm^2) and is ablating the donor (by using its spin-down power) thus generating mass-loss with an efficiency of 40% or the donor star is undergoing quasi-cyclic variations due to a varying mass-quadrupole induced by either a strong (1 kG) field or by some unidentified mechanism probably linked to irradiation.

Transitional millisecond pulsars (tMSPs) switch, on roughly multi-year timescales, between rotation-powered radio millisecond pulsar (RMSP) and accretion-powered low-mass X-ray binary (LMXB) states. The tMSPs have raised several questions related to the nature of accretion flow in their LMXB state and the mechanism that causes the state switch. The discovery of coherent X-ray pulsations from PSR J1023+0038 (while in the LMXB state) provides us with the first opportunity to perform timing observations and to compare the neutron star's spin variation during this state to the measured spin-down in the RMSP state. Whereas the X-ray pulsations in the LMXB state likely indicate that some material is accreting onto the neutron star's magnetic polar caps, radio continuum observations indicate the presence of an outflow. The fraction of the inflowing material being ejected is not clear, but it may be much larger than that reaching the neutron star's surface. Timing observations can measure the total torque on the neutron star. We have phase-connected nine XMM-Newton observations of PSR J1023+0038 over the last 2.5 years of the LMXB state to establish a precise measurement of spin evolution. We find that the average spin-down rate as an LMXB is 26.8+/-0.4% faster than the rate (-2.39x10^-15 Hz s-1) determined during the RMSP state. This shows that negative angular momentum contributions (dipolar magnetic braking and outflow) exceed positive ones (accreted material), and suggests that the pulsar wind continues to operate at a largely unmodified level. We discuss implications of this tight observational constraint in the context of possible accretion models.