Pith. sign in

REVIEW 3 major objections 6 minor 74 references

Identifying Compact Chirping SMBHBs in LSST using Bayesian Analysis

T0 review · 3 major / 6 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read LSST quasar light curves alone can reveal the gravitational-wave chirp of a supermassive black-hole binary and measure its properties.

desk verdict Solid injection-recovery study with a genuinely new fast sampler; the headline 'establish on its own' outruns the evidence because the quoted z-scores are within-model posterior signal-to-noise, not calibrated false-alarm rates. read the letter →

arxiv 2506.10846 v1 pith:W2NEJM5R submitted 2025-06-12 astro-ph.HE

Chengcheng Xin , Maximiliano Isi , Will M. Farr , Zoltán Haiman This is my paper · ORCID
classification astro-ph.HE
keywords supermassiveblackholebinariesLSSTchirpBayesianinferencedampedrandomwalkDopplerboostingquasarvariabilitygravitationalwaves
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper aims to show that LSST, the upcoming ten-year survey of the southern sky, can on its own identify the "chirp" -- the accelerating orbital frequency of a supermassive black-hole binary as it is driven together by gravitational waves -- from quasar light curves, without needing a prior gravitational-wave detection by LISA or pulsar timing arrays. The authors generate mock LSST light curves that combine a Doppler-boosted sinusoidal chirp, damped-random-walk quasar noise, and Gaussian photometric errors, and analyze them with a fully Bayesian framework. They find that for chirp amplitudes of $A = 0.5$ mag and times to merger of $t_m = 15$--$10^4$ yr, the amplitude and positive frequency derivative are typically measured with over $5\sigma$ credibility, and that even for $t_m = 50$ yr, chirping is established at $3\sigma$ for $A \gtrsim 0.1$ mag and $5\sigma$ for $A > 0.2$ mag. Because the analysis takes only minutes per light curve, it is scalable to the roughly 100 million quasars LSST is expected to observe, turning the survey into a stand-alone discovery machine for compact binaries that LISA and pulsar timing arrays could later confirm.

What carries the argument

The load-bearing machinery is the "frequency comb" sampler: a hybrid Hamiltonian Monte Carlo and Gibbs scheme in which auxiliary grids of the orbital frequency $f_0$ and its derivative $\dot{f}_0$ are laid down with spacings equal to the known separations between the likelihood's false sideband maxima ($\Delta\omega_0 \approx 7.72525/T$ in frequency and $\Delta\dot{\omega}_0 \approx 15.121/T^2$ in frequency derivative, for an observation span $T$), and the sampler Gibbs-selects a grid point while the grid itself slides smoothly under an HMC proposal. This lets the Markov chain hop between the oscillatory likelihood's many modes instead of getting stuck in a sideband. The linear parameters (chirp amplitude and phase through $a = A\cos\phi$, $b = A\sin\phi$, and the mean magnitude) are analytically marginalized using the Gaussian factorization of the likelihood, which reduces the sampled dimension from seven to four and speeds each light-curve analysis to minutes.

What would settle it

Run the pipeline on simulated pure-noise LSST light curves drawn from a damped harmonic oscillator with no injected chirp and count the fraction that pass the $>5\sigma$ chirp threshold; if that fraction is orders of magnitude above the $\lesssim 10^{-16}$ implied by the paper's $z$-scores under the damped-random-walk model, the central detectability claim fails for realistic noise. Equivalently, apply the pipeline to early real LSST quasar light curves selected to lack known periodicity and check whether $\dot{f}_0 > 0$ detections appear at the predicted rate.

Watch

Extended reading notes

Core claim

The central claim is that a compact supermassive black-hole binary in an LSST quasar can be identified purely from its electromagnetic light curve, by measuring the gravitational-wave-driven frequency evolution of its orbital modulation. Modeling the light curve as a mean magnitude plus a constant-amplitude sinusoidal Doppler-boost chirp (with phase from leading-order post-Newtonian evolution, $\Phi(t) = -\left(\frac{t_m - t}{5 t_M}\right)^{5/8}$), a damped random walk with covariance $\sigma^2 \exp(-|t_1 - t_2|/\tau)$, and Gaussian photometric errors, the authors perform a fully Bayesian joint inference of the seven parameters $\{f_0, \dot{f}_0, A, \phi, \sigma, \tau, m_i\}$. In mock observations with realistic six-day cadence and seasonal gaps, the chirp amplitude $A$ and frequency derivative $\dot{f}_0$ are recovered with almost no correlation with the noise parameters, with $z$-scores typically between 10 and 30 for an injected amplitude of 0.5 mag. The authors report that a non-zero chirp can be measured at $>5\sigma$ for all simulated binaries with $t_m = 15$--$10^4$ yr at $A = 0.5$ mag; the widest system for which the signal itself is detected has $P_0 \approx 1850$ d, and the widest for which the chirp $\dot{f}_0 \neq 0$ is measured has $P_0 \approx 200$ d. For binaries with $t_m = 50$ yr, chirping is established at $3\sigma$ for $A \gtrsim 0.1$ mag and $5\sigma$ for $A > 0.2$ mag, with the redshifted chirp mass and time to merger typically recovered to better than 10% fractional error.

Load-bearing premise

The result stands on the assumption that real quasar stochastic variability is exactly a damped random walk with the empirically calibrated amplitude and timescale scaling, and that real binary light curves are exactly constant-amplitude sinusoidal leading-order post-Newtonian chirps; the paper's own Sections 4.2 and 4.3 state that if either assumption fails, for instance under a damped-harmonic-oscillator noise model or an AGN-disk-precession chirp-like signal, the quoted false-alarm probabilities would not hold.

Editorial extensions

If this is right

  • LSST can act as a stand-alone electromagnetic discoverer of compact supermassive black-hole binaries: a measured chirp is smoking-gun evidence for gravitational-wave-driven inspiral, independent of LISA or pulsar-timing-array detections.
  • The bulk of the sources that make up the nHz stochastic gravitational-wave background detected by pulsar timing arrays should be identifiable in LSST quasar catalogs on their own.
  • Day-to-week-period binaries detected as chirping quasars can be flagged as LISA targets before LISA launches, enabling targeted searches and multimessenger observations once LISA is operating.
  • The redshifted chirp mass $(1+z)\mathcal{M}$ and the time to merger $t_m$ are typically recovered with fractional errors of $\lesssim 10\%$ for detectable systems, so the survey yields astrophysical population constraints, not just detections.
  • The per-light-curve runtime of typically under 10 minutes makes a search over the full $\sim 10^8$-quasar LSST sample computationally feasible, at an estimated $\sim 16$ million CPU hours.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The quoted false-alarm probabilities (z-scores of 8--30 corresponding to $\lesssim 10^{-16}$) are only valid under the assumed damped-random-walk plus sinusoidal-chirp generative model; the paper itself concedes in Section 4.3 that alternative noise models such as a damped harmonic oscillator would dominate the uncertainty. Applying this pipeline to real LSST data will therefore require model comp
  • If real quasar variability has a high-frequency component that the damped random walk misses, the detectability thresholds for low-amplitude chirps could shift; a natural test is to inject chirps into light curves generated from a damped harmonic oscillator model and re-measure the amplitude thresholds for $3\sigma$ and $5\sigma$ detections.
  • Because the paper assumes pure gravitational-wave-driven inspiral, any circumbinary-disk coupling would change both the true chirp rate and the inferred time to merger; the model would need extra free parameters to absorb that uncertainty.
  • A null result -- applying the full pipeline to the LSST quasar sample and finding no $>5\sigma$ chirping quasars -- would itself be informative: given the expected counts of ${\sim}150$ to $10^5$ binaries with $t_m$ between 15 and $10^4$ yr, it would constrain the fraction of quasars associated with binaries or the amplitude distribution of Doppler-boost variability.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The paper develops a fully Bayesian framework to detect and characterize gravitational-wave-driven 'chirp' signals in simulated LSST quasar lightcurves. The mock data combine a post-Newtonian Doppler-boosting chirp, a damped random walk representing stochastic quasar variability, and Gaussian photometric errors, with realistic cadence and seasonal gaps. The authors demonstrate simultaneous inference of seven parameters (chirp amplitude, phase, frequency, frequency derivative, DRW amplitude and timescale, and mean magnitude) using a custom 'frequency comb' sampling strategy. For injected chirp amplitude A=0.5 mag and merger times t_m = 15-10^4 yr, they report amplitude and frequency-derivative z-scores typically above 5, and they show that for t_m=50 yr, 3-5 sigma chirp detection is possible for A > 0.1-0.2. The paper concludes that LSST could on its own establish the presence of compact supermassive black-hole binaries and identify LISA/PTA-relevant sources.

Significance. If the reported detection statistics are taken at face value, this would be a significant methodological advance: it would demonstrate that a large time-domain survey can identify compact SMBHBs and even measure their orbital frequency evolution, providing an electromagnetic counterpart channel for LISA and PTA sources. The paper's strengths include a careful injection-recovery setup, a novel frequency-comb sampling scheme that addresses multimodal likelihoods, analytic marginalization over linear parameters, and a clear effort to quantify computational cost, with typical runtimes under 10 minutes per lightcurve. The framework is a credible step toward scalable searches in LSST-sized catalogs. However, the headline claim depends crucially on the assumed generative model being exactly right, and the paper's own discussion in Sections 4.2 and 4.3 concedes that model uncertainties would destroy the quoted false-alarm probabilities. The central scientific contribution is therefore better characterized as a proof-of-concept recovery study than as an established discovery claim.

major comments (3)
  1. [Sec. 3 and Sec. 4.3] The reported z-scores are posterior signal-to-noise ratios conditional on the assumed generative model (chirp + DRW + known Gaussian noise), not empirical false-alarm rates. Every simulated lightcurve contains a real injected chirp; no noise-only (A=0) or pure-DRW lightcurves are run through the pipeline. Consequently, the abstract's claim that LSST 'could, on its own, establish the presence of a compact supermassive black-hole binary' and the statement in Sec. 4.3 of 'extremely small false alarm probability' are not supported by the presented statistics. The authors should run a null-injection study (with A=0, and possibly with alternative noise models) and report the distribution of z-scores under the null hypothesis. This is a load-bearing issue because the headline claim is a discovery claim rather than a recovery claim.
  2. [Appendix A.1] The photometric noise variances sigma_{n,k} are fixed to the exact values used to generate the mock data. In a real LSST analysis these variances must be estimated from the data or marginalized over, and fixing them to the truth will generally inflate the reported detection significance. The authors should either marginalize over the noise variances or quantify how much the z-scores degrade when the variances are inferred with realistic priors. This is directly relevant to the 'establish on its own' claim and should be addressed before publication.
  3. [Sec. 4.2 and Sec. 4.3] The paper's own discussion states that alternative noise models (e.g., damped harmonic oscillator) and non-GW chirp-like evolution (e.g., AGN-disk precession in 'tick-tock') could mimic or dominate the signal, and that such uncertainties would preclude the quoted false-alarm probabilities. Because the inference templates include only the GW-driven chirp and DRW, the high z-scores do not demonstrate robustness against these alternatives. The authors should either include alternative templates in the analysis or explicitly frame the results as conditional on the adopted model family, revising the abstract and conclusions accordingly. As written, the central claim overreaches what the simulations establish.
minor comments (6)
  1. [Sec. 3 (Fig. 4, Fig. 10)] The color scale for z-scores is logarithmic in some panels and linear in others (e.g., Fig. 10), and the ranges differ widely across panels; this makes cross-panel comparison difficult and should be clarified in the captions.
  2. [Sec. 2.3, Eq. (6)-(7)] The derivation of the prior widths from Delta log(M) = 1 is clear, but the text should state explicitly that these priors are not meant to represent the full astrophysical uncertainty in mass and Eddington ratio; the choice of a symmetric log-normal spread may affect the tails of the posterior z-scores.
  3. [Sec. 3.3 and Fig. 7] The statement that S[x] 'plunges to ~0' for A_t=0.05 because the posterior is close to the prior is confusing: the prior mean for positive parameters is not zero, and a fractional error near zero for a non-detection is not a sign of accuracy. The text should be reworded to explain that the posterior is dominated by the prior, not that the measurement is accurate.
  4. [Sec. 4.3] The sentence 'generating more DRW realizations will give very similar results' is not demonstrated and is not equivalent to a false-alarm analysis; it would be more appropriate to present a real ensemble of noise-only simulations.
  5. [Sec. 4 (Summary)] The phrase 'z-score = 5 implies a 5 sigma detection' is imprecise because these are posterior z-scores, not frequentist significance levels; the authors should consistently refer to them as posterior signal-to-noise ratios or credible-region statements.
  6. [References] Several references are incomplete or inconsistently formatted (e.g., 'M. C. Davis et al. 2024', 'Kis-Toth & Haiman 2025' has an odd author field, and some arXiv entries lack titles); these should be cleaned up.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: injection-recovery study with independent priors; self-citations are contextual and not load-bearing.

full rationale

The paper is an explicitly labeled injection-recovery study: mock lightcurves are generated from Eq. (2), combining the DB chirp (Eqs. 3-4), a DRW (Eq. 5), and Gaussian noise, and the Bayesian likelihood in Appendix A is the same generative model. This means the reported z-scores quantify parameter-recovery confidence under the assumed model rather than empirical false-alarm rates. That is a statistical limitation, and the paper itself concedes it in Sec. 4.3 ('both the binary model ... and the noise model are inherently uncertain, and these uncertainties will dominate the z-score, precluding such low false alarm probabilities') and Sec. 4.2 (non-GW chirp-like evolution such as AGN-disk precession could mimic the signal). However, no fitted constant is inserted into the claims: the priors on f0 and fdot0 follow from an independent Eddington-ratio scatter estimate (Eqs. 6-7), the DRW parameters come from external empirical relations (MacLeod et al. 2010), and the chirp phase is leading-order PN. The self-citations (XH21, XH24) supply population numbers and implementation details, but the central detectability result does not reduce to them; there is no uniqueness theorem or ansatz smuggled in via self-citation. The discrepancy between 'posterior z-score' and 'false-alarm probability' is a correctness/calibration concern, not a circular derivation. Therefore no circular step is exhibited.

Assumptions & free parameters 6 free parameters · 5 assumptions · 0 invented entities

No new physical entities are postulated. The claims rest on the assumed generative model (post-Newtonian chirp plus DRW plus known noise), on source-count assumptions (f_bin = 1, t_Q = 1e7 yr, f_Edd = 0.3), and on the simplification that photometric errors are known exactly. These are domain assumptions whose uncertainty is acknowledged in Section 4.

free parameters (6)
  • q (binary mass ratio) = 0.1
    Fixed at 0.1 for all simulated binaries in Section 2.1. The chirp rate and amplitude depend on q through the chirp mass M = [q/(1+q)^2]^(3/5) M_bin.
  • f_Edd (Eddington ratio) = 0.3
    Used in Eq. (1) to map i-band magnitude to binary mass and thus to the simulated period and chirp rate distribution.
  • t_Q (quasar lifetime) = 1e7 yr
    Used in the counting argument N proportional to t_m/t_Q in Section 2.1 to estimate the number of compact binaries LSST will see; source counts scale linearly with this choice.
  • f_bin (fraction of quasars hosting SMBHBs) = 1
    Assumed equal to unity in Section 2.1, which maximizes the predicted binary yield and therefore the expected number of observable chirps.
  • Fiducial chirp amplitude A_t = 0.5 mag (scanned down to 0.05 mag)
    The main forecast in Section 3.1 assumes A = 0.5 mag, at the upper end of Doppler boost and hydrodynamical predictions. The amplitude dependence is explicitly explored in Section 3.2, so this is a scanned signal parameter rather than an unexamined fit.
  • Photometric error variances sigma_n,k = Drawn once, then fixed at injected values during inference
    Appendix A.1 treats the per-point error variances as known and fixed at the exact values that generated the mock data. This idealizes the analysis relative to real LSST data, where photometric errors are estimated with uncertainty.
assumptions (5)
  • domain assumption Leading-order post-Newtonian gravitational-wave-driven chirp (f_dot proportional to f^(11/3)) describes the binary orbital evolution over the final 10^4 years.
    Used to derive the chirp phase in Eq. (4), while ignoring gaseous torques and eccentricity; Section 4.1 notes these can accelerate or modify the inspiral.
  • domain assumption Quasar intrinsic variability is a damped random walk with exponential kernel (Eq. 5) and MacLeod et al. (2010) amplitude-timescale scaling.
    Both the simulated data and the likelihood assume DRW; Section 4.3 acknowledges alternative noise models such as the damped harmonic oscillator.
  • domain assumption Binary variability in the optical band is a sinusoid from Doppler boosting of a circular orbit, with constant amplitude A over the 10-year observation.
    Used in Eq. (3); real signals may be saw-toothed (hydrodynamic) or spiky (lensing), and the amplitude may grow toward merger, as discussed in Section 4.1.
  • domain assumption The extrapolated quasar luminosity function of Kulkarni et al. (2019) describes the LSST quasar population down to i = 26.
    Basis for the source sampling in Figure 1 and for the predicted binary counts in Section 2.1.
  • domain assumption LSST i-band cadence is one visit every 6 days with seasonal gaps that skip one third of the data each year.
    Used to simulate observing times in Section 2.2; the real LSST cadence varies by field and band, as noted in Section 4.4.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Identifying Compact Chirping SMBHBs in LSST using Bayesian Analysis." pith.science (2026). https://pith.science/paper/W2NEJM5R

@misc{pith2026250610846,
  author       = {Pith},
  title        = {Pith review of: Identifying Compact Chirping SMBHBs in LSST using Bayesian Analysis},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/W2NEJM5R}},
  note         = {Machine review of arXiv:2506.10846}
}
abstract

The Legacy Survey of Space and Time (LSST) is expected to observe up to ${\sim}100$ million quasars in the next decade. In this work, we show that it is possible to use such data to measure the characteristic frequency evolution of a "chirp" induced by gravitational waves, which can serve as robust evidence for the presence of a compact supermassive black-hole binary. Following the LSST specifications, we generate mock lightcurves consisting of (i) a post-Newtonian chirp produced by orbital motion through, e.g., relativistic Doppler boosting, (ii) a damped random walk representing intrinsic quasar variability, and (iii) Gaussian photometric errors, while assuming non-uniform observations with extended gaps over a period of 10 yr. Through a fully-Bayesian analysis, we show that we can simultaneously measure the chirp and noise parameters with little degeneracy between the two. For chirp signals with an amplitude of $A = 0.5$ mag and a range of times to merger ($t_m = 15{-}10^4$ yr), we can typically measure a non-zero amplitude and positive frequency derivative with over $5\sigma$ credibility. For binaries with $t_m = 50$ yr, we achieve $3\sigma$ ($5\sigma$) confidence that the signal is chirping for $A \gtrsim 0.1$ ($A > 0.2$). Our analysis can take as little as 35 s (and typically $<$ 10 min) to run, making it scalable to a large number of lightcurves. This implies that LSST could, on its own, establish the presence of a compact supermassive black-hole binary, and thus discover gravitational wave sources detectable by LISA and by Pulsar Timing Arrays.

Figures

Figures reproduced from arXiv: 2506.10846 by the authors.

Figure 1
Figure 1. Expected number density of binary quasars in LSST (color) as a function of total binary mass (abscissa) and redshift (ordinate), using the extrapolated quasar luminosity function from G. Kulkarni et al. (2019) and assuming that the LSST co-added i-band magnitude threshold is mi = 26 (left edge of the contour). The total number of expected LSST quasars is ∼120 million. the quasar, and (iii) Gaussian photometric error… view at source ↗
Figure 2
Figure 2. Example mock LSST lightcurve for a SMBHB with Mbin = 107.8M⊙ and z = 2.9. The lightcurve (black points) has a mean magnitude mi = 24.2 per Eq. (1) (line), modulated by the binary through relativistic DB with a chirping sinusoid, Eq. (3) with A = 0.5 mag (blue points); additionally, the quasar varies stochastically follow￾ing a DRW, Eq. (5) with σ = 0.5 mag and τ = 136 d (orange points); finally, the photometric erro… view at source ↗
Figure 3
Figure 3. Posterior for all seven parameters for the binary lightcurve in [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (8 more)
Figure 4
Figure 4. Figure 4: Summary of the inference on the 100 lightcurves simulated for tm =15, 50, 100, 103 and 104 yr (panels). We show posterior z-scores, mean(A)/std(A), for the chirp amplitude A (right colorbar), as a function of binary mass (abscissa) and redshift (ordinate), as well as t…
Figure 5
Figure 5. Figure 5: The detectability of binary chirp amplitude (A) as a function of the true injected amplitude At, for nine binaries in our sample. The z-score increases as At increases. It has a weak dependence on orbital periods (quoted in the legend). The chirp rate is fixed at 7.5% …
Figure 6
Figure 6. Figure 6: The significance at which the frequency chirp ˙f0 can be measured as a function of variability amplitude At (abscissa) and period P0 (color). The y-axis shows the posterior z-score for ˙f0, i.e., mean( ˙f0)/std( ˙f0), calculated for the same set of At and P0 as in [PI…
Figure 7
Figure 7. Figure 7: Precision and accuracy with which we can measure the observed chirp mass (1 + z)M (top left panel) and time– to-merger tm (bottom left panel), as a function of the true binary amplitude At (abscissa) and period P0 (color). The left column shows precision through log-st…
Figure 8
Figure 8. Figure 8: Illustration of the frequency comb construction in frequencies and frequency derivative space (ω0, ω˙ 0). The HMC sampler first picks a lower-left corner (w init 0 , w˙ init 0 ) near (ωmin, ω˙ min), following Eq. (A8), which is then used to construct an Nω × Nω˙ grid w…
Figure 9
Figure 9. Figure 9: Effective priors for the frequency ω0 (left) and chirp amplitude A (right). (a) The effective prior on ω0 results from the sampling strategy associated with Eq. (A10) which defines a frequency grid of Nω0 points separated by ∆ω0 and starting at values not smaller than …
Figure 10
Figure 10. Figure 10: We show the chirp ( ˙f0) posterior z-score’s for the same binaries presented in [PITH_FULL_IMAGE:figures/full_fig_p016_10.png]
Figure 11
Figure 11. Figure 11: The runtime of our Bayesian inference model for each lightcurve in the fiducial binary population (top right panel in [PITH_FULL_IMAGE:figures/full_fig_p017_11.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

74 extracted references · 9 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

  3. [3]

    - [1] #1 = = ^ ^ ^ .\!\!^ d .\!\!^ h .\!\!^ m .\!\!^ s .\!\!^ @mss

    thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...

    work page 2017

  4. [4]

    M., et al

    Agazie, G., Anumarlapudi, A., Archibald, A. M., et al. 2023, title The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Background , The Astrophysical Journal Letters, 951, L8, 10.3847/2041-8213/acdac6

    doi:10.3847/2041-8213/acdac6 2023

  5. [5]

    M., & Hickox , R

    Alexander , D. M., & Hickox , R. C. 2012, title What drives the growth of black holes? , , 56, 93, 10.1016/j.newar.2011.11.003

    doi:10.1016/j.newar.2011.11.003 2012

  6. [6]

    C., Blandford , R

    Begelman , M. C., Blandford , R. D., & Rees , M. J. 1980, title Massive black hole binaries in active galactic nuclei , , 287, 307, 10.1038/287307a0

    doi:10.1038/287307a0 1980

  7. [7]

    P., Jankowiak, M., et al

    Bingham, E., Chen, J. P., Jankowiak, M., et al. 2019, title Pyro: Deep Universal Probabilistic Programming, J. Mach. Learn. Res., 20, 28:1. http://jmlr.org/papers/v20/18-403.html

    work page 2019

  8. [8]

    Charisi, M., Bartos, I., Haiman, Z., et al. 2016, title A population of short-period variable quasars from PTF as supermassive black hole binary candidates , Monthly Notices of the Royal Astronomical Society, 463, 2145, 10.1093/mnras/stw1838

    doi:10.1093/mnras/stw1838 2016

Show all 74 references
  1. [9]

    R., Runnoe, J., Bogdanovic, T., & Trump, J

    Charisi, M., Taylor, S. R., Runnoe, J., Bogdanovic, T., & Trump, J. R. 2022, title Multimessenger time-domain signatures of supermassive black hole binaries , Monthly Notices of the Royal Astronomical Society, 510, 5929. https://doi.org/10.1093/mnras/stab3713

    2022 doi

  2. [10]

    Chen, Y.-C., Liu, X., Liao, W.-T., et al. 2020, title Candidate Periodically Variable Quasars from the Dark Energy Survey and the Sloan Digital Sky Survey , Monthly Notices of the Royal Astronomical Society, 21, 1, 10.1093/mnras/staa2957

    2020 doi

  3. [11]

    2024, title Electromagnetic signatures from accreting massive black hole binaries in time domain photometric surveys , A&A, 691, 1, 10.1051/0004-6361/202449598

    Cocchiararo, F., Franchini, A., Lupi, A., & Sesana, A. 2024, title Electromagnetic signatures from accreting massive black hole binaries in time domain photometric surveys , A&A, 691, 1, 10.1051/0004-6361/202449598

    2024 doi

  4. [12]

    2024, title LISA Definition Study Report , ESA LISA Definition Study Report; e-print arXiv:2402.07571

    Colpi , M., Danzmann , K., Hewitson , M., et al. 2024, title LISA Definition Study Report , ESA LISA Definition Study Report; e-print arXiv:2402.07571

    2024 arXiv

  5. [13]

    Cutler, C., & Flanagan, E. E. 1994, title Gravitational waves from merging compact binaries: How accurately can one extract the binary's parameters from the inspiral wave form? , Phys. Rev. D, 49, 2658, 10.1103/PhysRevD.49.2658

    1994 doi

  6. [14]

    Davelaar, J., & Haiman, Z. 2022 a , title Self-Lensing Flares from Black Hole Binaries: Observing Black Hole Shadows via Light Curve Tomography , Physical Review Letters, 128, 191101, 10.1103/PhysRevLett.128.191101

    2022 doi

  7. [15]

    2022 b , title Self-lensing flares from black hole binaries: General-relativistic ray tracing of black hole binaries , Physical Review D, 105, 103010, 10.1103/PhysRevD.105.103010

    Davelaar, J., & Haiman, Z. 2022 b , title Self-lensing flares from black hole binaries: General-relativistic ray tracing of black hole binaries , Physical Review D, 105, 103010, 10.1103/PhysRevD.105.103010

    2022 doi

  8. [16]

    C., Grace, K

    Davis, M. C., Grace, K. E., Trump, J. R., et al. 2024, title Reliable Identification of Binary Supermassive Black Holes from Rubin Observatory Time-domain Monitoring , The Astrophysical Journal, 965, 34, 10.3847/1538-4357/ad276e

    2024 doi

  9. [17]

    J., Ryan, G., & Miller, M

    Dittmann, A. J., Ryan, G., & Miller, M. C. 2023, title The Decoupling of Binaries from Their Circumbinary Disks , The Astrophysical Journal Letters, 949, L30, 10.3847/2041-8213/acd183

    2023 doi

  10. [18]

    J., & Di Stefano , R

    D'Orazio, D. J., & Di Stefano , R. 2018, title Periodic self-lensing from accreting massive black hole binaries , Monthly Notices of the Royal Astronomical Society, 474, 2975, 10.1093/mnras/stx2936

    2018 doi

  11. [19]

    J., Duffell, P

    D'Orazio, D. J., Duffell, P. C., & Tiede, C. 2024, title Fast Methods for Computing Photometric Variability of Eccentric Binaries: Boosting, Lensing, and Variable Accretion , The Astronomical Journal, 977, 244, 10.3847/1538-4357/ad938b

    2024 doi

  12. [20]

    J., Haiman, Z., & Schiminovich, D

    D'Orazio, D. J., Haiman, Z., & Schiminovich, D. 2015, title Relativistic boost as the cause of periodicity in a massive black-hole binary candidate , Nature, 525, 351, 10.1038/nature15262

    2015 doi

  13. [21]

    2023, title Optical follow-up of the tick-tock massive black hole binary candidate , Mnras, 518, 4172, 10.1093/mnras/stac3344

    Dotti, M., Bonetti, M., Rigamonti, F., et al. 2023, title Optical follow-up of the tick-tock massive black hole binary candidate , Mnras, 518, 4172, 10.1093/mnras/stac3344

    2023 doi

  14. [22]

    D., Duffell, P., MacFadyen, A

    Farris, B. D., Duffell, P., MacFadyen, A. I., & Haiman, Z. 2015, title Binary black hole accretion during inspiral and merger , Monthly Notices of the Royal Astronomical Society: Letters, 447, L80, 10.1093/mnrasl/slu184

    2015 doi

  15. [23]

    2018, title Scalable Backpropagation for Gaussian Processes using Celerite , Research Notes of the American Astronomical Society, 2, 31, 10.3847/2515-5172/aaaf6c

    Foreman-Mackey , D. 2018, title Scalable Backpropagation for Gaussian Processes using Celerite , Research Notes of the American Astronomical Society, 2, 31, 10.3847/2515-5172/aaaf6c

    2018 doi

  16. [24]

    2017, title Fast and Scalable Gaussian Process Modeling with Applications to Astronomical Time Series , , 154, 220, 10.3847/1538-3881/aa9332

    Foreman-Mackey , D., Agol , E., Ambikasaran , S., & Angus , R. 2017, title Fast and Scalable Gaussian Process Modeling with Applications to Astronomical Time Series , , 154, 220, 10.3847/1538-3881/aa9332

    2017 doi

  17. [25]

    A., Mignon-Risse, R., et al

    Foustoul, V., Webb, N. A., Mignon-Risse, R., et al. 2025, title A catalogue of candidate milli-parsec separation massive black hole binaries from long term optical photometric monitoring , 2505.06656

    2025 arXiv

  18. [26]

    2024, title Gas-induced perturbations on the gravitational wave in-spiral of live post-Newtonian LISA massive black hole binaries , arXiv:2410.17305v1

    Garg, M., Franchini, A., Lupi, A., Bonetti, M., & Mayer, L. 2024, title Gas-induced perturbations on the gravitational wave in-spiral of live post-Newtonian LISA massive black hole binaries , arXiv:2410.17305v1

    2024

  19. [27]

    D., Greene , J

    Goulding , A. D., Greene , J. E., Bezanson , R., et al. 2018, title Galaxy interactions trigger rapid black hole growth: An unprecedented view from the Hyper Suprime-Cam survey , , 70, S37, 10.1093/pasj/psx135

    2018 doi

  20. [28]

    J., Djorgovski, S

    Graham, M. J., Djorgovski, S. G., Stern, D., et al. 2015 a , title A systematic search for close supermassive black hole binaries in the Catalina real-time transient survey , Monthly Notices of the Royal Astronomical Society, 453, 1562, 10.1093/mnras/stv1726

    2015 doi

  21. [29]

    J., Djorgovski, S

    Graham, M. J., Djorgovski, S. G., Stern, D., et al. 2015 b , title A possible close supermassive black-hole binary in a quasar with optical periodicity , Nature, 518, 74, 10.1038/nature14143

    2015 doi

  22. [30]

    2017, title Electromagnetic chirp of a compact binary black hole: A phase template for the gravitational wave inspiral , Physical Review D, 2, 023004, 10.1103/PhysRevD.96.023004

    Haiman, Z. 2017, title Electromagnetic chirp of a compact binary black hole: A phase template for the gravitational wave inspiral , Physical Review D, 2, 023004, 10.1103/PhysRevD.96.023004

    2017 doi

  23. [31]

    Haiman, Z., Kocsis, B., & Menou, K. 2009, title The population of viscosity- and gravitational wave-driven supermassive black hole binaries among luminous active galactic nuclei , Astrophysical Journal, 700, 1952, 10.1088/0004-637X/700/2/1952

    2009 doi

  24. [32]

    W., Price-Whelan, A

    Hogg, D. W., Price-Whelan, A. M., & Leistedt, B. 2020, title Data Analysis Recipes: Products of multivariate Gaussians in Bayesian inferences , 2005.14199

    2020 arXiv

  25. [33]

    F., & Hernquist, L

    Hopkins, P. F., & Hernquist, L. 2009, title Quasars are not light bulbs: Testing models of quasar lifetimes with the observed Eddington ratio distribution , The Astrophysical Journal, 698, 1550, 10.1088/0004-637X/698/2/1550

    2009 doi

  26. [34]

    2025, title Periodic Variability in Space Photometry of 181 New Supermassive Black Hole Binary Candidates , , 1 2505.16884

    Huijse, P., Davelaar, J., De Ridder , J., Jannsen, N., & Aerts, C. 2025, title Periodic Variability in Space Photometry of 181 New Supermassive Black Hole Binary Candidates , , 1 2505.16884

    2025 arXiv

  27. [35]

    E., Aigrain, S., & Karastergiou, A

    Ingram, A., Motta, S. E., Aigrain, S., & Karastergiou, A. 2021, title A self-lensing binary massive black hole interpretation of quasi-periodic eruptions , Monthly Notices of the Royal Astronomical Society, 503, 1703, 10.1093/mnras/stab609

    2021 doi

  28. [36]

    2023, title Parametrizing gravitational-wave polarizations , Class

    Isi, M. 2023, title Parametrizing gravitational-wave polarizations , Class. Quant. Grav., 40, 203001, 10.1088/1361-6382/acf28c

    2023 doi

  29. [37]

    Ivezic, Z., & MacLeod, C. L. 2013, title Optical variability of quasars: a damped random walk , Proceedings of the International Astronomical Union, 9, 395 , 10.1017/S1743921314004396

    2013 doi

  30. [38]

    2022, title Tick-Tock: The Imminent Merger of a Supermassive Black Hole Binary , 2201.11633

    Jiang, N., Yang, H., Wang, T., et al. 2022, title Tick-Tock: The Imminent Merger of a Supermassive Black Hole Binary , 2201.11633

    2022 arXiv

  31. [39]

    2000, title A unified model for the evolution of galaxies and quasars , , 311, 576, 10.1046/j.1365-8711.2000.03077.x

    Kauffmann , G., & Haehnelt , M. 2000, title A unified model for the evolution of galaxies and quasars , , 311, 576, 10.1046/j.1365-8711.2000.03077.x

    2000

  32. [40]

    C., Bechtold, J., & Siemiginowska, A

    Kelly, B. C., Bechtold, J., & Siemiginowska, A. 2009, title Are the variations in quasar optical flux driven by thermal fluctuations? , , 698, 895, 10.1088/0004-637X/698/1/895

    2009 doi

  33. [41]

    2025, title Can quasars, triggered by mergers, account for NANOGrav's stochastic gravitational wave background? , Classical and Quantum Gravity, 42, 19pp

    Kis-T \' o th, \' A ., & Haiman, Zoltan Frei, Z. 2025, title Can quasars, triggered by mergers, account for NANOGrav's stochastic gravitational wave background? , Classical and Quantum Gravity, 42, 19pp

    2025

  34. [42]

    S., Eisenstein, D

    Kochanek, C. S., Eisenstein, D. J., Cool, R. J., et al. 2012, title AGES: The AGN and galaxy evolution survey , The Astrophysical Journal Supplement Series, 200, 10.1088/0067-0049/200/1/8

    2012 doi

  35. [43]

    A., Onken , C

    Kollmeier , J. A., Onken , C. A., Kochanek , C. S., et al. 2006, title Black Hole Masses and Eddington Ratios at 0.3 < z < 4 , The Astrophysical Journal, 648, 128, 10.1086/505646

    2006 doi

  36. [44]

    M., Davelaar , J., Haiman , Z., et al

    Krauth , L. M., Davelaar , J., Haiman , Z., et al. 2024, title Self-lensing flares from black hole binaries: General-relativistic ray tracing of circumbinary accretion simulations , Physical Review D, 109, 103014, 10.1103/PhysRevD.109.103014

    2024 doi

  37. [45]

    Kulkarni, G., Worseck, G., & Hennawi, J. F. 2019, title Evolution of the AGN UV luminosity function from redshift 7.5 , Monthly Notices of the Royal Astronomical Society, 31, 1035, 10.1093/mnras/stz1493

    2019 doi

  38. [46]

    Li, K., Bogdanovi \' c , T., & Ballantyne, D. R. 2020, title Pairing of Massive Black Holes in Merger Galaxies Driven by Dynamical Friction , The Astrophysical Journal, 896, 113, 10.3847/1538-4357/ab93c6

    2020 doi

  39. [47]

    R., & Bonetti, M

    Li, K., Bogdanovi \' c , T., Ballantyne, D. R., & Bonetti, M. 2022, title Massive Black Hole Binaries from the TNG50-3 Simulation. I. Coalescence and LISA Detection Rates , The Astrophysical Journal, 933, 104, 10.3847/1538-4357/ac74b5

    2022 doi

  40. [48]

    Lomb, N. R. 1976, title Least-Squares Frequency Analysis of Unequally Spaced Data , Astrophysics and Space Science, 39, 447, 10.1007/BF00648343

    1976 doi

  41. [49]

    A., Allison, J., et al

    LSST Science Collaboration , Abell, P. A., Allison, J., et al. 2009, title LSST Science Book, Version 2.0, e-print arXiv:0912.0201. 0912.0201

    2009 arXiv

  42. [50]

    L., Ivezi \' c , Z ., Kochanek, C

    MacLeod, C. L., Ivezi \' c , Z ., Kochanek, C. S., et al. 2010, title Modeling the time variability of SDSS Stripe 82 quasars as a damped random walk , Astrophysical Journal, 721, 1014, 10.1088/0004-637X/721/2/1014

    2010 doi

  43. [51]

    2004, in Coevolution of Black Holes and Galaxies, ed

    Martini , P. 2004, in Coevolution of Black Holes and Galaxies, ed. L. C. Ho , 169. astro-ph/0304009

    2004 arXiv

  44. [52]

    F., Edelson, R., Baumgartner, W., & Gandhi, P

    Mushotzky, R. F., Edelson, R., Baumgartner, W., & Gandhi, P. 2011, title Kepler observations of rapid optical variability in active galactic nuclei , Astrophysical Journal Letters, 743, 6, 10.1088/2041-8205/743/1/L12

    2011 doi

  45. [53]

    J., Haiman, Z., & Macfadyen, A

    O'Neill, D., Tiede, C., D'Orazio, D. J., Haiman, Z., & Macfadyen, A. 2025, title Gravitational Wave Decoupling in Retrograde Circumbinary Disks , arXiv:2501.11679v1

    2025 arXiv

  46. [54]

    P. C. Peters , & Mathews, J. 1963, title Gravitational Radiation from Point Masses in a Keplerian Orbit , Physical Review, 131, 435. papers2://publication/uuid/746C5D1C-C087-4723-AA8D-3C20B5B506E7

    1963

  47. [55]

    2024, title Self-lensing flares from black hole binaries IV: the number of detectable shadows , , 1 2409.04583

    Park, K., Xin, C., Davelaar, J., & Haiman, Z. 2024, title Self-lensing flares from black hole binaries IV: the number of detectable shadows , , 1 2409.04583

    2024 arXiv

  48. [56]

    2019, title Composable Effects for Flexible and Accelerated Probabilistic Programming in NumPyro , arXiv e-prints, arXiv:1912.11554, 10.48550/arXiv.1912.11554

    Phan , D., Pradhan , N., & Jankowiak , M. 2019, title Composable Effects for Flexible and Accelerated Probabilistic Programming in NumPyro , arXiv e-prints, arXiv:1912.11554, 10.48550/arXiv.1912.11554

    2019 doi arXiv

  49. [57]

    E., & Williams , C

    Rasmussen , C. E., & Williams , C. K. I. 2006, Gaussian Processes for Machine Learning

    2006

  50. [58]

    E., Charisi, M., et al

    Robnik, J., Bayer, A. E., Charisi, M., et al. 2024, title Periodicity significance testing with null-signal templates: reassessment of PTF's SMBH binary candidates , Mnras, 534, 1609, 10.1093/mnras/stae2220

    2024 doi

  51. [59]

    Scargle, J. D. 1982, title STUDIES IN ASTRONOMICAL TIME SERIES ANALYSIS. II. STATISTICAL ASPECTS OF SPECTRAL ANALYSIS OF UNEVENLY SPACED DATA , Astrophysical Journal, 263, 835, 10.1086/160554

    1982 doi

  52. [60]

    2023, title Orbital evolution of binaries in circumbinary discs , Monthly Notices of the Royal Astronomical Society, 522, 2707, 10.1093/mnras/stad1131

    Siwek, M., Weinberger, R., & Hernquist, L. 2023, title Orbital evolution of binaries in circumbinary discs , Monthly Notices of the Royal Astronomical Society, 522, 2707, 10.1093/mnras/stad1131

    2023 doi

  53. [61]

    L., Mushotzky, R

    Smith, K. L., Mushotzky, R. F., Boyd, P. T., et al. 2018, title The Kepler Light Curves of AGN: A Detailed Analysis , The Astrophysical Journal, 857, 141, 10.3847/1538-4357/aab88d

    2018 doi

  54. [62]

    Tang, Y., Haiman, Z., & MacFadyen, A. 2018, title The late inspiral of supermassive black hole binaries with circumbinary gas discs in the LISA band , Monthly Notices of the Royal Astronomical Society, 476, 2249, 10.1093/mnras/sty423

    2018 doi

  55. [63]

    Tiede, C., & Orazio, D. J. D. 2024, title Eccentric binaries in retrograde discs , Monthly Notices of the Royal Astronomical Society, 527, 6021, 10.1093/mnras/stad3551

    2024 doi

  56. [64]

    Van Rossum, G., & Drake, F. L. 2009, Python 3 Reference Manual (Scotts Valley, CA: CreateSpace)

    2009

  57. [65]

    G., et al

    Vaughan, S., Uttley, P., Markowitz, A. G., et al. 2016, title False periodicities in quasar time-domain surveys , Monthly Notices of the Royal Astronomical Society, 461, 1, 10.1093/mnras/stw1412

    2016 doi

  58. [66]

    2024, title TomWagg/software-citation-station: v1.2, , v1.2 Zenodo, 10.5281/zenodo.13225824

    Wagg, T., Broekgaarden, F., & Gültekin, K. 2024, title TomWagg/software-citation-station: v1.2, , v1.2 Zenodo, 10.5281/zenodo.13225824

    2024 doi

  59. [67]

    Wagg , T., & Broekgaarden , F. S. 2024, title Streamlining and standardizing software citations with The Software Citation Station , arXiv e-prints, arXiv:2406.04405. 2406.04405

    2024

  60. [68]

    R., Zrake, J., MacFadyen, A., & Haiman, Z

    Westernacher-Schneider, J. R., Zrake, J., MacFadyen, A., & Haiman, Z. 2022, title Multiband light curves from eccentric accreting supermassive black hole binaries , Physical Review D, 106, 10.1103/physrevd.106.103010

    2022 doi

  61. [69]

    A., Charisi, M., Taylor, S

    Witt, C. A., Charisi, M., Taylor, S. R., & Burke-Spolaor, S. 2022, title Quasars with Periodic Variability: Capabilities and Limitations of Bayesian Searches for Supermassive Black Hole Binaries in Time-domain Surveys , The Astrophysical Journal, 936, 89, 10.3847/1538-4357/ac8356

    2022 doi

  62. [70]

    Xin, C., Charisi, M., Haiman, Z., et al. 2020, title Testing the relativistic Doppler boost hypothesis for the binary candidate quasar PG1302-102 with multiband Swift data , Monthly Notices of the Royal Astronomical Society, 496, 1683, 10.1093/mnras/staa1643

    2020 doi

  63. [71]

    Xin, C., & Haiman, Z. 2021, title Ultra-short-period massive black hole binary candidates in LSST as LISA 'verification binaries' , Monthly Notices of the Royal Astronomical Society, 506, 2408, 10.1093/mnras/stab1856

    2021 doi

  64. [72]

    Xin, C., & Haiman, Z. 2024, title Identifying the electromagnetic counterparts of LISA massive black hole binaries in archi v al LSST data , Monthly Notices of the Royal Astronomical Society, 533, 3164, 10.1093/mnras/stae2009

    2024 doi

  65. [73]

    T., Vogeley, M

    Yu, W., Richards, G. T., Vogeley, M. S., Moreno, J., & Graham, M. J. 2022, title Examining AGN UV / Optical Variability beyond the Simple Damped Random Walk , The Astrophysical Journal, 936, 132, 10.3847/1538-4357/ac8351

    2022 doi

  66. [74]

    2021, title Equilibrium Eccentricity of Accreting Binaries , The Astrophysical Journal Letters, 909, L13, 10.3847/2041-8213/abdd1c

    Zrake, J., Tiede, C., MacFadyen, A., & Haiman, Z. 2021, title Equilibrium Eccentricity of Accreting Binaries , The Astrophysical Journal Letters, 909, L13, 10.3847/2041-8213/abdd1c

    2021 doi

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.