The bright quasar PG1302-102 has been identified as a candidate supermassive black hole binary from its near-sinusoidal optical variability. While the significance of its optical periodicity has been debated due to the stochastic variability of quasars, its multi-wavelength variability in the ultraviolet (UV) and optical bands is consistent with relativistic Doppler boost caused by the orbital motion in a binary. However, this conclusion was based previously on sparse UV data which were not taken simultaneously with the optical data. Here we report simultaneous follow-up observations of PG1302-102 with the Ultraviolet Optical Telescope on the Neil Gehrels Swift Observatory in six optical + UV bands. The additional nine Swift observations produce light curves roughly consistent with the trend under the Doppler boost hypothesis, which predicts that UV variability should track the optical, but with a ~2.2 times higher amplitude. We perform a statistical analysis to quantitatively test this hypothesis. We find that the data are consistent with the Doppler boost hypothesis when we compare the the amplitudes in optical B-band and UV light curves. However, the ratio of UV to V-band variability is larger than expected and is consistent with the Doppler model, only if either the UV/optical spectral slopes vary, the stochastic variability makes a large contribution in the UV, or the sparse new optical data underestimate the true optical variability. We have evidence for the latter from comparison with the optical light curve from ASAS-SN. Additionally, the simultaneous analysis of all four bands strongly disfavors the Doppler boost model whenever Swift V-band is involved. Additional, simultaneous optical + UV observations tracing out another cycle of the 5.2-year proposed periodicity should lead to a definitive conclusion.

We examine the light curves of two quasars, motivated by recent suggestions that a supermassive black hole binary (SMBHB) can exhibit sharp lensing spikes. We model the variability of each light curve as due to a combination of two relativistic effects: the orbital relativistic Doppler boost and gravitational binary self-lensing. In order to model each system we extend previous Doppler plus self-lensing models to include eccentricity. The first quasar is identified in optical data as a binary candidate with a 20-yr period (Ark 120), and shows a prominent spike. For this source, we rule out the lensing hypothesis and disfavor the Doppler-boost hypothesis due to discrepancies in the measured vs. recovered values of the binary mass and optical spectral slope. The second source, which we nickname Spikey, is the rare case of an active galactic nucleus (AGN) identified in Kepler's high-quality, high-cadence photometric data. For this source, we find a model, consisting of a combination of a Doppler modulation and a narrow symmetric lensing spike, consistent with an eccentric SMBHB with a total mass of approximately 30 million solar masses, rest-frame orbital period T=418 days, eccentricity e=0.5, and seen at an inclination of 8 degrees from edge-on. This interpretation can be tested by monitoring Spikey for periodic behavior and recurring flares in the next few years. In preparation for such monitoring we present the first X-ray observations of this object taken by the Neil Gehrels Swift observatory.

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field.

The variability of quasars across multiple wavelengths is a useful probe of physical conditions in active galactic nuclei. In particular, variable accretion rates, instabilities, and reverberation effects in the accretion disk of a supermassive black hole (SMBH) are expected to produce correlated flux variations in UV and optical bands. Recent work has further argued that binary quasars should exhibit strongly correlated UV and optical periodicities. Strong UV-optical correlations have indeed been established in small samples of up to approximately 30 quasars with well-sampled light curves, and have extended the "bluer-when-brighter" trend previously found within the optical bands. Here we further test the nature of quasar variability by examining the observed-frame UV-optical correlations in a large sample of 1,315 bright quasars with overlapping UV and optical light curves for the Galaxy Evolution Explorer (GALEX) and the Catalina Real-time Transient Survey (CRTS), respectively. We find that strong correlations exist in this much larger sample, but we rule out, at approximately 95% confidence, the simple hypothesis that the intrinsic UV and optical variations of all quasars are fully correlated. Our results therefore imply the existence of physical mechanism(s) that can generate uncorrelated optical and UV flux variations.

Pulsar timing array (PTA) collaborations in North America, Australia, and Europe, have been exploiting the exquisite timing precision of millisecond pulsars over decades of observations to search for correlated timing deviations induced by gravitational waves (GWs). PTAs are sensitive to the frequency band ranging just below 1 nanohertz to a few tens of microhertz. The discovery space of this band is potentially rich with populations of inspiraling supermassive black-holes binaries, decaying cosmic string networks, relic post-inflation GWs, and even non-GW imprints of axionic dark matter. This article aims to provide an understanding of the exciting open science questions in cosmology, galaxy evolution, and fundamental physics that will be addressed by the detection and study of GWs through PTAs. The focus of the article is on providing an understanding of the mechanisms by which PTAs can address specific questions in these fields, and to outline some of the subtleties and difficulties in each case. The material included is weighted most heavily towards the questions which we expect will be answered in the near-term with PTAs; however, we have made efforts to include most currently anticipated applications of nanohertz GWs.

Precision timing of large arrays (>50) of millisecond pulsars will detect the nanohertz gravitational-wave emission from supermassive binary black holes within the next ~3-7 years. We review the scientific opportunities of these detections, the requirements for success, and the synergies with electromagnetic instruments operating in the 2020s.

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.

LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center.

Supermassive black hole binaries (SMBHBs) should be common in galactic nuclei as a result of frequent galaxy mergers. Recently, a large sample of sub-parsec SMBHB candidates was identified as bright periodically variable quasars in optical surveys. If the observed periodicity corresponds to the redshifted binary orbital period, the inferred orbital velocities are relativistic (v/c~0.1). The optical and UV luminosities are expected to arise from gas bound to the individual BHs, and would be modulated by the relativistic Doppler effect. The optical and UV light curves should vary in tandem with relative amplitudes which depend on the respective spectral slopes. We constructed a control sample of 42 quasars with aperiodic variability, to test whether this Doppler colour signature can be distinguished from intrinsic chromatic variability. We found that the Doppler signature can arise by chance in ~20% (~37%) of quasars in the nUV (fUV) band. These probabilities reflect the limited quality of the control sample and represent upper limits on how frequently quasars mimic the Doppler brightness+colour variations. We performed separate tests on the periodic quasar candidates, and found that for the majority, the Doppler boost hypothesis requires an unusually steep UV spectrum or an unexpectedly large BH mass and orbital velocity. We conclude that at most ~1/3rd of these periodic candidates can harbor Doppler-modulated SMBHBs.

Supermassive black hole binaries (SMBHBs) at sub-parsec separations should be common in galactic nuclei, as a result of frequent galaxy mergers. Hydrodynamical simulations of circumbinary discs predict strong periodic modulation of the mass accretion rate on time-scales comparable to the orbital period of the binary. As a result, SMBHBs may be recognized by the periodic modulation of their brightness. We conducted a statistical search for periodic variability in a sample of 35,383 spectroscopically confirmed quasars in the photometric database of the Palomar Transient Factory (PTF). We analysed Lomb-Scargle periodograms and assessed the significance of our findings by modeling each individual quasar's variability as a damped random walk (DRW). We identified 50 quasars with significant periodicity beyond the DRW model, typically with short periods of a few hundred days. We find 33 of these to remain significant after a re-analysis of their periodograms including additional optical data from the intermediate-PTF and the Catalina Real-Time Transient Survey (CRTS). Assuming that the observed periods correspond to the redshifted orbital periods of SMBHBs, we conclude that our findings are consistent with a population of unequal-mass SMBHBs, with a typical mass ratio as low as q = M2/M1 ~ 0.01.

Graham et al. (2015) discovered a supermassive black hole binary (SMBHB) candidate and identified the detected 5.2 yr period of the optical variability as the orbital period of the binary. Hydrodynamical simulations predict multiple periodic components for the variability of SMBHBs, thus raising the possibility that the true period of the binary is different from 5.2 yr. We analyse the periodogram of PG1302 and find no compelling evidence for additional peaks. We derive upper limits on any additional periodic modulations in the available data, by modeling the light-curve as the sum of red noise and the known 5.2 yr periodic component, and injecting additional sinusoidal signals. We find that, with the current data, we would be able to detect with high significance (false alarm probability <1%) secondary periodic terms with periods in the range predicted by the simulations, if the amplitude of the variability was at least ~0.07 mag (compared to 0.14 mag for the main peak). A three-year follow-up monitoring campaign with weekly observations can increase the sensitivity for detecting secondary peaks ~50%, and would allow a more robust test of predictions from hydrodynamical simulations.

We report a comprehensive catalog of emission episodes within long gamma-ray bursts (GRBs) that are separated by a quiescent period during which gamma-ray emission falls below the background level. We use a fully automated identification method for an unbiased, large scale and expandable search. We examine a comprehensive sample of long GRBs from the BATSE, Swift and Fermi missions, assembling a total searched set of 2710 GRBs, the largest catalog of isolated emission episodes so far. Our search extends out to [-1000s,750s] around the burst trigger, expanding the covered time interval beyond previous studies and far beyond the nominal durations (T90) of most bursts. We compare our results to previous works by identifying pre-peak emission (or precursors), defined as isolated emission periods prior to the episode with the highest peak luminosity of the burst. We also systematically search for similarly defined periods after the burst's peak emission. We find that the pre-peak and post-peak emission periods are statistically similar, possibly indicating a common origin. For the analyzed GRBs, we identify 24% to have more than one isolated emission episode, with 11% having at least one pre-peak event and 15% having at least one post-peak event. We identify GRB activity significantly beyond their T90, which can be important for understanding the central engine activity as well as, e.g., gravitational-wave searches.