skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Strongly Lensed Supermassive Black Hole Binaries as Nanohertz Gravitational-wave Sources
Abstract Supermassive black hole binary systems (SMBHBs) should be the most powerful sources of gravitational waves (GWs) in the universe. Once pulsar timing arrays (PTAs) detect the stochastic GW background from their cosmic merger history, searching for individually resolvable binaries will take on new importance. Since these individual SMBHBs are expected to be rare, here we explore how strong gravitational lensing can act as a tool for increasing their detection prospects by magnifying fainter sources and bringing them into view. Unlike for electromagnetic waves, when the geometric optics limit is nearly always valid, for GWs the wave-diffraction-interference effects can become important when the wavelength of the GWs is larger than the Schwarzchild radius of the lens, i.e., M lens 10 8 f mHz 1 M . For the GW frequency range explored in this work, the geometric optics limit holds. We investigate GW signals from SMBHBs that might be detectable with current and future PTAs under the assumption that quasars serve as bright beacons that signal a recent merger. Using the black hole mass function derived from quasars and a physically motivated magnification distribution, we expect to detect a few strongly lensed binary systems out toz≈ 2. Additionally, for a range of fixed magnifications 2 ≤μ≤ 100, strong lensing adds up to ∼30 more detectable binaries for PTAs. Finally, we investigate the possibility of observing both time-delayed electromagnetic signals and GW signals from these strongly lensed binary systems—that will provide us with unprecedented multi-messenger insights into their orbital evolution.  more » « less
Award ID(s):
2020265
PAR ID:
10462151
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
955
Issue:
1
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 25
Size(s):
Article No. 25
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The merger of two galaxies, each hosting a supermassive black hole (SMBH) of mass 106Mor more, could yield a bound SMBH binary. For the early-type galaxy NGC 4472, we study how astrometry with a next-generation Very Large Array could be used to monitor the reflex motion of the primary SMBH of massMpri, as it is tugged on by the secondary SMBH of mass M sec . Casting the orbit of the putative SMBH binary in terms of its periodP, semimajor axisabin, and mass ratio q = M sec / M pri 1 , we find the following: (1) Orbits with fiducial periods ofP= 4 yr and 40 yr could be spatially resolved and monitored. (2) For a 95% accuracy of 2μas per monitoring epoch, subparsec values ofabincould be accessed over a range of mass ratios notionally encompassing major q > 1 4 and minor q < 1 4 galaxy mergers. (3) If no reflex motion is detected forMpriafter 1 (10) yr of monitoring, an SMBH binary with periodP= 4 (40) yr and mass ratioq> 0.01 (0.003) could be excluded. This would suggest no present-day evidence for a past major merger like that recently simulated, where scouring by aq∼ 1 SMBH binary formed a stellar core with kinematic traits like those of NGC 4472. (4) Astrometric monitoring could independently check the upper limits onqfrom searches for continuous gravitational waves from NGC 4472. 
    more » « less
  2. Abstract We present the results of an all-sky search for continuous gravitational waves in the public LIGO O3 data. The search covers signal frequencies 20.0 Hz ≤f≤ 800.0 Hz and a spin-down range down to −2.6 × 10−9Hz s−1, motivated by detectability studies on synthetic populations of Galactic neutron stars. This search is the most sensitive all-sky search to date in this frequency/spin-down region. The initial search was performed using the first half of the public LIGO O3 data (O3a), utilizing graphical processing units provided in equal parts by the volunteers of the Einstein@Home computing project and by the ATLAS cluster. After a hierarchical follow-up in seven stages, 12 candidates remain. Six are discarded at the eighth stage, by using the remaining O3 LIGO data (O3b). The surviving six can be ascribed to continuous-wave fake signals present in the LIGO data for validation purposes. We recover these fake signals with very high accuracy with our last stage search, which coherently combines all O3 data. Based on our results, we set upper limits on the gravitational-wave amplitudeh0and translate these into upper limits on the neutron star ellipticity and on ther-mode amplitude. The most stringent upper limits are at 203 Hz, withh0= 8.1 × 10−26at the 90% confidence level. Our results exclude isolated neutron stars rotating faster than 5 ms with ellipticities greater than 5 × 10 8 d 100 pc within a distancedfrom Earth andr-mode amplitudes α 10 5 d 100 pc for neutron stars spinning faster than 150 Hz. 
    more » « less
  3. Abstract LetXbe ann-element point set in thek-dimensional unit cube$$[0,1]^k$$ [ 0 , 1 ] k where$$k \ge 2$$ k 2 . According to an old result of Bollobás and Meir (Oper Res Lett 11:19–21, 1992) , there exists a cycle (tour)$$x_1, x_2, \ldots , x_n$$ x 1 , x 2 , , x n through thenpoints, such that$$\left( \sum _{i=1}^n |x_i - x_{i+1}|^k \right) ^{1/k} \le c_k$$ i = 1 n | x i - x i + 1 | k 1 / k c k , where$$|x-y|$$ | x - y | is the Euclidean distance betweenxandy, and$$c_k$$ c k is an absolute constant that depends only onk, where$$x_{n+1} \equiv x_1$$ x n + 1 x 1 . From the other direction, for every$$k \ge 2$$ k 2 and$$n \ge 2$$ n 2 , there existnpoints in$$[0,1]^k$$ [ 0 , 1 ] k , such that their shortest tour satisfies$$\left( \sum _{i=1}^n |x_i - x_{i+1}|^k \right) ^{1/k} = 2^{1/k} \cdot \sqrt{k}$$ i = 1 n | x i - x i + 1 | k 1 / k = 2 1 / k · k . For the plane, the best constant is$$c_2=2$$ c 2 = 2 and this is the only exact value known. Bollobás and Meir showed that one can take$$c_k = 9 \left( \frac{2}{3} \right) ^{1/k} \cdot \sqrt{k}$$ c k = 9 2 3 1 / k · k for every$$k \ge 3$$ k 3 and conjectured that the best constant is$$c_k = 2^{1/k} \cdot \sqrt{k}$$ c k = 2 1 / k · k , for every$$k \ge 2$$ k 2 . Here we significantly improve the upper bound and show that one can take$$c_k = 3 \sqrt{5} \left( \frac{2}{3} \right) ^{1/k} \cdot \sqrt{k}$$ c k = 3 5 2 3 1 / k · k or$$c_k = 2.91 \sqrt{k} \ (1+o_k(1))$$ c k = 2.91 k ( 1 + o k ( 1 ) ) . Our bounds are constructive. We also show that$$c_3 \ge 2^{7/6}$$ c 3 2 7 / 6 , which disproves the conjecture for$$k=3$$ k = 3 . Connections to matching problems, power assignment problems, related problems, including algorithms, are discussed in this context. A slightly revised version of the Bollobás–Meir conjecture is proposed. 
    more » « less
  4. 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 f GW 2 / 3 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 26 19 + 28 nHz , 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
  5. A<sc>bstract</sc> We present a quantum M2 brane computation of the instanton prefactor in the leading non-perturbative contribution to the ABJM 3-sphere free energy at largeNand fixed levelk. Using supersymmetric localization, such instanton contribution was found earlier to take the form$$ {F}^{inst}\left(N,k\right)=-{\left({\sin}^2\frac{2\pi }{k}\right)}^{-1}\exp \left(-2\pi \sqrt{\frac{2N}{k}}\right)+.\dots $$ F inst N k = sin 2 2 π k 1 exp 2 π 2 N k + . The exponent comes from the action of an M2 brane instanton wrapped onS3/ℤk, which represents the M-theory uplift of the ℂP1instanton in type IIA string theory on AdS4× ℂP3. The IIA string computation of the leading largekterm in the instanton prefactor was recently performed in arXiv:2304.12340. Here we find that the exact value of the prefactor$$ {\left({\sin}^2\frac{2\pi }{k}\right)}^{-1} $$ sin 2 2 π k 1 is reproduced by the 1-loop term in the M2 brane partition function expanded near theS3/ℤkinstanton configuration. As in the Wilson loop example in arXiv:2303.15207, the quantum M2 brane computation is well defined and produces a finite result in exact agreement with localization. 
    more » « less