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. This write-up is an initiative taken within the framework of the European Action on ‘Black holes, Gravitational waves and Fundamental Physics’.

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.39×10 −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.

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 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 (O3). 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 area of 4480 deg2, a median distance of 267 Mpc, and false-alarm rates ranging from 1.5 to 10−25 yr−1. The ZTF coverage in the g and r bands had a median enclosed probability of 39%, median depth of 20.8 mag, and median time lag between merger and the start of observations of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UltraViolet/Optical/InfraRed (UVOIR) photometric points, 64 OIR spectra, and three radio images using 17 different telescopes. We find no promising kilonovae (radioactivity-powered counterparts), and we show how to convert the upper limits to constrain the underlying kilonova luminosity function. Initially, we assume that all GW triggers are bona fide astrophysical events regardless of false-alarm rate and that kilonovae accompanying BNS and NSBH mergers are drawn from a common population; later, we relax these assumptions. Assuming that all kilonovae are at least as luminous as the discovery magnitude of GW170817 (−16.1 mag), we calculate that our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than −16.6 mag (the extrapolated peak magnitude of GW170817) and fade at a rate of 1 mag day−1 (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations based on the online classifications, the joint probability of zero detections, assuming all kilonovae are brighter than −16.6 mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, no more than <57% (<89%) of putative kilonovae could be brighter than −16.6 mag assuming flat evolution (fading by 1 mag day−1) at the 90% confidence level. If we further take into account the online terrestrial probability for each GW trigger, we find that no more than <68% of putative kilonovae could be brighter than −16.6 mag. Comparing to model grids, we find that some kilonovae must have M ej < 0.03 M ⊙, X lan > 10−4, or ϕ > 30° to be consistent with our limits. We look forward to searches in the fourth GW observing run; even 17 neutron star mergers with only 50% coverage to a depth of −16 mag would constrain the maximum fraction of bright kilonovae to <25%.

We analyzed broadband 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 blackbody +power-law (BB+PL) or 2BB model during all three outbursts. Nuclear Spectroscopic Telescope Array observations revealed a hard X-ray tail, with a PL photon index Γ=0.9, extending up to 50 keV, with flux comparable to 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 ∼7 following its strongest 2016 June outburst. A Swift/X-Ray Telescope 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 quasiexponentially with timescales of ∼40 days, while the stronger 2016 May and June outbursts showed a quick short-term decay with timescales of about four days. Our Arecibo radio observations set the deepest limits on the radio emission from a magnetar, with a maximum flux density limit of 14 μJy for the 4.6GHz observations and 7 μJy for the 1.4GHz observations. We discuss these results in the framework of the current magnetar theoretical models.

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 deg 2 field at ∆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, 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.

The accreting millisecond X-ray pulsar SAX J1808.4−3658 shows a peculiar orbital evolution that proceeds at a very fast pace. It is important to identify the underlying mechanism responsible for this behavior because it can help to understand how this system evolves and which physical processes (such asmass loss or spin-orbit coupling) are occurring in the binary. 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 accreting millisecond pulsars and black-widow (BW) radio pulsars. In this work, we report the results of a deep radio pulsation search at 2 GHz using the Green Bank Telescope in 2014 August and an X-ray study of the 2015 outburst with Chandra, SwiftXRT, and INTEGRAL. In quiescence, we detect no radio pulsationsand place the strongest limit to date on the pulsed radio flux density of any accreting millisecond pulsar. We also find that the orbit of SAX J1808.4−3658 continues evolving at a fast pace. We compare the orbital evolution of SAX J1808.4−3658 to that of several other accreting and nonaccreting binaries, including BWs, redbacks, cataclysmic variables, black holes, and neutron stars in lowmass X-ray binaries. We discuss two possible scenarios: either the neutron star has a large moment of inertia and is ablating the donor, generating massloss with an efficiency of 40%, or the donor star has a strong magnetic field of at least 1 kG and is undergoing quasi-cyclic variations due to spin-orbit coupling.

keV) 21,760 cm2 Energy Resolution 85 -175 eV FWHM Time Resolution 100 ns Collimator 4 arcmin FWHM Background Rate 2.2 c/s Count Rate on Crab (0.2-10 keV) 148,000 Large Area Detector (LAD) Energy Range 2-30 keV Effective Area (cm^2 @ 10 keV) 51,000 cm2 Energy Resolution 200 -300 eV FWHM Time Resolution 10 µs Collimator 1° FWHM Count Rate on Crab (2-30 keV) 156,000 Background Rate 822 c/s (5 mcrab)