A Multimessenger Search for the Supermassive Black Hole Binary in 3C 66B with the Parkes Pulsar Timing Array

Using the Parkes Pulsar Timing Array's third data release, this study fails to confirm the supermassive black hole binary candidate in 3C 66B, instead establishing upper limits that partially rule out electromagnetic parameter estimates while proposing a new joint likelihood framework for multimessenger analysis and standard siren applications.

The Gist

Imagine the universe is a giant, silent ocean. For a long time, we thought this ocean was perfectly still. But recently, scientists realized the whole ocean is actually rippling with a low, constant hum caused by massive objects crashing into each other. This is the Gravitational Wave Background.

Now, imagine you are trying to hear a single, specific voice singing a note in the middle of that humming crowd. That is what this paper is about.

Here is the story of the search for that specific voice, told simply:

1. The Target: A Cosmic Dance Floor

Deep in space, in a galaxy called 3C 66B, astronomers think there are two supermassive black holes (the heaviest, most gravity-filled objects in the universe) locked in a slow, tight dance. They are so close they are less than a light-year apart (a "subparsec" distance).

Because they are so heavy and moving in a circle, they should be creating a very specific, loud "squeak" in the fabric of space-time. Based on observations of light and radio waves (electromagnetic data), scientists predicted exactly what this "squeak" should sound like: its pitch (frequency) and how loud it should be.

2. The Microphone: The Parkes Pulsar Timing Array

To hear this cosmic squeak, we can't use a normal microphone. We need a telescope that listens to pulsars.

Think of pulsars as the universe's most perfect metronomes. They are dead stars that spin hundreds of times a second, sending out radio beams like lighthouse beacons. If space-time is perfectly smooth, these beacons arrive at Earth with perfect, clockwork precision.

But if a gravitational wave passes through, it stretches and squeezes space. This makes the "ticks" of the metronome arrive a tiny fraction of a second early or late. The Parkes Pulsar Timing Array (PPTA) is a team of astronomers using the giant "Murriyang" radio telescope in Australia to listen to 32 of these cosmic metronomes simultaneously, looking for that tiny wobble.

3. The Search: Tuning the Radio

The scientists in this paper took the latest data from their 18-year listening session (called "DR3"). They didn't just listen to the whole ocean; they tuned their radio specifically to the frequency predicted by the 3C 66B black hole dance.

They used two methods to check their results:

• The Bayesian Method: This is like asking, "How much more likely is it that we are hearing the black hole voice compared to just random static?"
• The Frequentist Method: This is like asking, "If there were no black hole, how often would we see a signal this strong just by pure luck?"

4. The Result: A "Maybe" Answer

Here is the twist: They didn't find the voice.

• The Verdict: The data didn't show a clear signal. The "Bayes Factor" (a score of confidence) was essentially zero. It's like listening to a crowded room and not hearing the specific singer you were looking for.
• The Good News: They also couldn't say the singer isn't there. The signal might just be too quiet for their current microphones.
• The Bad News: They did rule out the "loudest" versions of the theory. If the black holes were as massive or as close as some earlier predictions suggested, they should have heard them. Since they didn't, those specific predictions are likely wrong.

The Analogy: Imagine you are looking for a specific person in a crowd. You know they are wearing a red hat. You scan the crowd and don't see anyone in a red hat. You can't say for sure the person isn't there (maybe they are hiding behind a pillar), but you can say for sure they aren't standing right in front of you in plain sight.

5. The New Idea: Using "Silence" as a Tool

Even though they didn't find the black hole, the paper proposes a brilliant new way to use this "silence."

If we do eventually find a confirmed black hole pair like this, we can use it as a "Standard Siren."

• How it works: The gravitational waves tell us exactly how far away the black hole is (like knowing how loud a siren is tells you how far away a fire truck is). The light from the galaxy tells us how fast it's moving away from us.
• The Goal: By combining these two, we can measure how fast the universe is expanding (the Hubble Constant) with incredible precision.

This paper suggests that instead of waiting for random black holes to be found, we should specifically hunt for the ones we already know about from light observations. If we find them, they become the perfect rulers for measuring the universe.

Summary

• What they did: Listened for a specific pair of dancing black holes in galaxy 3C 66B using 18 years of data from Australian radio telescopes.
• What they found: No clear signal. The black holes are either quieter than expected, or they don't exist in the way we thought.
• Why it matters: They narrowed down the possibilities, ruling out the "loudest" theories. They also showed a new path forward: using these specific cosmic dances to measure the expansion of the universe, turning a failed search into a future cosmological breakthrough.

Technical Summary

Here is a detailed technical summary of the paper "A Multimessenger Search for the Supermassive Black Hole Binary in 3C 66B with the Parkes Pulsar Timing Array."

1. Problem Statement

The detection of continuous gravitational waves (CGWs) from individual supermassive black hole binaries (SMBHBs) is a primary goal for Pulsar Timing Arrays (PTAs). The galaxy 3C 66B is a leading candidate for such a source, based on electromagnetic (EM) observations suggesting a subparsec binary with a specific orbital frequency and chirp mass.

• The Conflict: Previous searches by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) have placed upper limits on the chirp mass that are compatible with, but do not confirm, the EM-derived parameters.
• The Goal: This study aims to perform the first targeted search for a CGW from 3C 66B using the Parkes Pulsar Timing Array (PPTA) Data Release 3 (DR3). The objective is to either confirm the EM model, rule it out, or place tighter constraints on the binary's parameters, while exploring the feasibility of using such targeted sources for cosmology (Standard Sirens).

2. Methodology

Data and Signal Model

• Dataset: PPTA DR3, comprising 18 years of timing data from 31 millisecond pulsars (excluding PSR J1741+1351 due to low priority/short span).
• Signal Model: The authors model a non-spinning SMBHB in a circular orbit at zeroth post-Newtonian order. The signal includes both the Earth term (metric perturbation at the Solar System Barycenter) and the Pulsar term (perturbation at the pulsar location).
• Approximation: The "stationary phase approximation" is used, neglecting frequency evolution over the 18-year span because the expected frequency drift ($\sim 0.04$ nHz/yr) is smaller than the PPTA frequency resolution ($\sim 1.75$ nHz).

Bayesian and Frequentist Analysis

• Bayesian Inference: The authors compute the posterior probability distributions for parameters (chirp mass $\mathcal{M}$, GW frequency $f_{GW}$, strain amplitude $h_0$, etc.) using parallel-tempered Markov-Chain Monte-Carlo (MCMC).
  • Priors: Two approaches were used:
    1. EM-Matched Priors: Priors centered on values derived from EM observations (Iguchi et al. 2010) to test the specific hypothesis.
    2. Broad/Conservative Priors: Uniform priors over wide ranges to derive model-independent upper limits.
• Model Selection: The Bayes Factor ($B$) is calculated to compare the signal model against the null hypothesis (noise only).
• Frequentist Search: An $F$-statistic ($F_e$) is calculated to provide a frequentist measure of evidence, yielding a false alarm probability (FAP) and $p$-value.
• Noise Modeling: The analysis incorporates white noise, red noise (power-law), and a common-spectrum process associated with the Gravitational Wave Background (GWB).
• Multimessenger Reanalysis: The authors independently re-analyzed the 3 mm flux monitor data from Iguchi et al. (2010) using a Bayesian framework to generate an EM posterior for the chirp mass, allowing for a direct comparison with GW posteriors.

Joint Likelihood Framework

The paper outlines a methodology for constructing a joint likelihood of EM and GW data. Since the chirp mass ($\mathcal{M}$) is a common parameter, the total likelihood is the product of the GW likelihood and the EM likelihood, enabling a unified constraint on the system.

3. Key Results

Detection Status

• No Detection: The search found no evidence for a CGW signal from 3C 66B.
• Bayes Factor: $\ln B = -0.0027(7)$ (or $B \approx 0.997$). This value indicates that the data is consistent with the noise-only hypothesis; the EM model is neither confirmed nor ruled out, but the specific parameter space is constrained.
• Frequentist Statistic: $F_e = 4.90$ (or $2F_e \approx 9.8$), corresponding to a$p$-value of$0.04$($\sim 2.1\sigma$). This is consistent with noise fluctuations.

Upper Limits (95% Credibility)

Using the PPTA DR3 data, the authors placed the following upper limits:

• Chirp Mass ($\mathcal{M}$):
  • With Earth term only: $\mathcal{M} < 8.25 \times 10^8 M_\odot$.
  • With Earth + Pulsar terms: $\mathcal{M} < \mathbf{6.90 \times 10^8 M_\odot}$.
  • Comparison: This limit rules out 58% of the uncertainty region of the EM model (Iguchi et al. 2010), which predicted $\mathcal{M} = 7.9^{+3.8}_{-4.5} \times 10^8 M_\odot$.
• Characteristic Strain Amplitude ($h_0$):
  • With Earth + Pulsar terms: $\log_{10}(h_0) < \mathbf{-14.44}$.
  • This partially rules out the strain amplitude predicted by the EM model ($7.3^{+6.8}_{-5.8} \times 10^{-15}$).

Frequency Dependence

The analysis showed that the PPTA is most sensitive to the frequency range suggested by the EM model ($f_{GW} \approx 60$ nHz). The upper limits on mass are tighter in this targeted search compared to all-sky searches.

4. Significance and Contributions

1. Tightest Constraints to Date: The PPTA DR3 analysis provides the most stringent upper limits on the chirp mass and strain amplitude for 3C 66B to date, surpassing previous NANOGrav results in sensitivity for this specific target.
2. Multimessenger Framework: The paper successfully demonstrates a new methodology for joint-likelihood analysis, combining EM flux data and GW timing data. This allows for a more rigorous test of SMBHB models than either dataset could provide alone.
3. Standard Siren Cosmology: The authors propose that targeted searches for known EM candidates (like 3C 66B) could bypass the localization difficulties of all-sky searches.
   • If a signal were detected, the known sky position and redshift would allow the GW signal to act as a "Standard Siren" to measure the Hubble-Lemaître constant ($H_0$).
   • Although no detection was made, the study showed that the absence of a signal can constrain the luminosity distance, thereby placing limits on $H_0$ (assuming the EM model is correct).
4. Validation of EM Models: The results highlight a tension between the EM-predicted parameters and the GW upper limits, suggesting that either the EM model overestimates the mass/strain, or the binary is not in the configuration assumed (e.g., eccentricity or inclination effects).

Conclusion

The study concludes that while the specific EM model for the 3C 66B binary is not confirmed, the PPTA DR3 data has significantly narrowed the viable parameter space. The work establishes a robust framework for future multimessenger searches and demonstrates the potential of targeted PTA searches to contribute to precision cosmology once a confident detection is achieved.