Abstract We recently discovered a supermassive black hole binary (SMBHB) candidate, PKS J1309+1154, in the combined 46 yr University of Michigan Radio Astronomy Observatory (UMRAO) plus Owens Valley Radio Observatory (OVRO) blazar monitoring programs at 14.5/15 GHz. The light curve of PKS 1309+1154 exhibits a 17.9 yr periodicity. We also reported a hint of a first harmonic with a 9 yr periodicity in this object. Further analysis of the PKS J1309+1154 light curve provides compelling evidence that both the fundamental and the harmonic are real, confirming the existence of real periodicities in blazar light curves. This is the first case, to our knowledge, of watertight evidence for a fundamental periodicity and a harmonic periodicity in a blazar light curve. It makes PKS J1309+1154 astrongSMBHB candidate and thus the third such candidate to be revealed through long-term radio monitoring, the other two being PKS J0805–0111 and PKS 2131–021, both discovered through the OVRO 40 m Telescope monitoring program. It is argued that hundreds of SMBHB candidates will be discovered by the combination of the South Pole Telescope, the Vera Rubin Observatory, and the Simons Observatory. Coherent searches for gravitational waves from a network of SMBHB candidates, starting immediately, are strongly motivated.
more »
« less
This content will become publicly available on March 12, 2027
A Search for Supermassive Black Hole Binary Candidates in 46 yr Radio Light Curves of 83 Blazars
Abstract The combined University of Michigan Radio Astronomy Observatory and Owens Valley Radio Observatory blazar monitoring programs at 14.5/15 GHz provide uninterrupted light curves of ∼ 46–50 yr duration for 83 blazars, selected from among the brightest and most rapidly flaring blazars north of declination −20°. In a search for supermassive black hole binary (SMBHB) candidates, we carried out tests for periodic variability using generalized Lomb–Scargle (GLS), weighted wavelet-Z, and sine-wave fitting analyses of this sample. We used simulations to test the effects of the power-law spectrum of the power spectral density (PSD) on our findings, and show that the irregular sampling in the observed light curves has very little effect on the GLS spectra. Apparent periodicities and putative harmonics appear in all 83 of the GLS spectra of the blazars in our sample. We tested the reality of these apparent periodicities and harmonics with simulations, and found that in the overwhelming majority of cases, they must be due to the steep slope of the PSD and the random nature of blazar flares, implying that, in general, apparent periodicities and harmonics in blazar light curves observed in any energy band should be treated with great caution. We find one new SMBHB candidate: PKS 1309+1154, which exhibits a 17.9 yr periodicity. The fraction of SMBHB candidates in our sample is .
more »
« less
- PAR ID:
- 10693796
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Publisher / Repository:
- Institute of Physics
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 1000
- Issue:
- 1
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 72
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Gamma-ray flares of blazars may be accompanied by high-energy neutrinos due to interactions of high-energy cosmic rays in the jet with photons, as suggested by the detection of the high-energy neutrino IceCube-170922A during a major gamma-ray flare from blazar TXS 0506+056 at the ∼3σsignificance level. In this work, we present a statistical study of gamma-ray emission from blazars to constrain the contribution of gamma-ray flares to their neutrino output. We construct weekly binned light curves for 145 gamma-ray bright blazars in the Fermi Large Area Telescope Monitored Source List adding TXS 0506+056. We derive the fraction of time spent in the flaring state (flare duty cycle) and the fraction of energy released during each flare from the light curves with a Bayesian blocks algorithm. We find that blazars with lower flare duty cycles and energy fractions are more numerous among our sample. We identify a significant difference in flare duty cycles between blazar subclasses at a significance level of 5%. Then using a general scaling relation for the neutrino and gamma-ray luminosities, with a weighting exponent ofγ= 1.0–2.0, normalized to the quiescent gamma-ray or X-ray flux of each blazar, we evaluate the neutrino energy flux of each gamma-ray flare. The gamma-ray flare distribution indicates that blazar neutrino emission may be dominated by flares forγ≳ 1.5. The neutrino energy fluxes for 1 week and 10 yr bins are compared with the decl.-dependent IceCube sensitivity to constrain the standard neutrino emission models for gamma-ray flares. Finally, we present the upper-limit contribution of blazar gamma-ray flares to the isotropic diffuse neutrino flux.more » « less
-
Abstract The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays (PTAs) through excursions from, and breaks in, the expected power law of the GWB strain spectrum. To do this, we create a semianalytic SMBHB population model, fit to North American Nanohertz Observatory for Gravitational Waves (NANOGrav’s) 15 yr GWB amplitude, and with 1000 realizations, we study the populations’ characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power law. The first, at 2 nHz, is below our GWB realizations with ap-value significancep= 0.05–0.06 (≈1.8σ–1.9σ). The second, at 16 nHz, is above our GWB realizations withp= 0.04–0.15 (≈1.4σ–2.1σ). We explore the properties of a loud SMBHB that could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by 3 orders of magnitude, from ∼106to ∼103, between 2 and 20 nHz. This causes a break in the strain spectrum as the stochasticity of the background breaks down at , consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the 26 nHz break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early Universe.more » « less
-
Abstract We study a magnitude-limited sample of 36 broad-lined type Ic supernovae (SNe Ic-BL) from the Zwicky Transient Facility Bright Transient Survey (detected between 2018 March and 2021 August), which is the largest systematic study of SNe Ic-BL done in literature thus far. We present the light curves (LCs) for each of the SNe and analyze the shape of the LCs to derive empirical parameters, along with the explosion epochs for every event. The sample has an average absolute peak magnitude in therband of mag. Using spectra obtained around peak light, we compute expansion velocities from the Feii5169 Å line for each event with high enough signal-to-noise ratio spectra, and find an average value of km s−1. We also compute bolometric LCs, study the blackbody temperature and radii evolution over time, and derive the explosion properties of the SNe. The explosion properties of the sample have average values of , , and erg. Thirteen events have radio observations from the Very Large Array, with eight detections and five non-detections. We find that the populations that have radio detections and radio non-detections are indistinct from one another with respect to their optically inferred explosion properties, and there are no statistically significant correlations present between the events’ radio luminosities and optically inferred explosion properties. This provides evidence that the explosion properties derived from optical data alone cannot give inferences about the radio properties of SNe Ic-BL and likely their relativistic jet formation mechanisms.more » « less
-
Abstract We combine new spectroscopic observations of the ultrafaint dwarf (UFD) galaxy Boötes I (Boo I) from the Southern Stellar Stream Spectroscopic Survey (S5) with ∼15 yr of archival spectroscopic data to create the largest sample of stellar kinematics and metallicities to date for any Milky Way UFD. Our combined sample includes 148 members extending out to ∼7 half-light radii (rh), including 24 newly confirmed members, 18 binary candidates, 15 RR Lyrae stars, and 92 [Fe/H] measurements. Using this larger and more spatially extended sample, we provide updated constraints on Boo I’s systemic properties, including its radial population gradients. Properly accounting for perspective rotation effects in a UFD for the first time, we detect a 4σline-of-sight velocity gradient of 1.2 ± 0.3 km s−1 aligned along Boo I’s orbit and discuss its potential tidal origins. We also infer a metallicity gradient of −0.10 ± 0.02 dex in agreement with previous studies. Using an axisymmetric Jeans model, we provide updated constraints on Boo I’s dark matter density profile, which weakly favors a cusped ( ) dark matter profile. Lastly, we reanalyze Boo I’s metallicity distribution function with a one-zone galactic chemical evolution model and place new constraints on its rapid, inefficient star formation and strong galactic outflows.more » « less
An official website of the United States government
