KAKaran Akbari

Play from the cross-band paper

Spot the binary

If two supermassive black holes orbit each other, the quasar around them should brighten and fade once per orbit. The catch is that every quasar also flickers randomly, and random flicker is good at faking a period. I searched 1,369 real ones for this. Here you get 8 at a time.

Round 1 of 3

Dealing the quasars.

Found so far, all rounds
0 of 0
Fakes so far
0

This round's quasars

tap one to look closer
Quasar A Drag to turn it
Light curve date 2012.2
Optical
+0.00 mag
X-ray
not measured

Light curve

dots: what a telescope caught · thin line: what the quasar did

Periodogram: how well a repeating wave fits at each period

dashed: textbook false-alarm line
Setting the false-alarm lines

Setting up.

Free play run hundreds of quasars at once and count how often the search is fooled

0of 0 quasars

counted as having a period

  • Optical went over the line0
  • Counted as a period0
  • nothing over the line
  • counted, and it's fake

500 d
0.13 mag
0.08 mag
Binary signal
0.20 mag
400 d

Changing a knob empties the jar, since the odds change with it. Real Swift-BAT points are mostly noise, so push the X-ray noise above σ to get closer to them.

What you're seeing

Every quasar in the game is simulated. The random flicker is a damped random walk, which is the usual model for how quasars vary (Kelly et al. 2009; MacLeod et al. 2010): the brightness wanders and slowly forgets where it was, over about 500 days here. Vaughan et al. (2016) argue that most periods claimed from single light curves are this kind of noise. A real binary adds a sine wave with a period somewhere between 150 and 750 days, 0.10 to 0.20 mag high, which is about the size of the flicker itself. You never see the answer until you press Check.

The 3D quasar is drawn from the same numbers as the plots. The disc's glow follows the optical light curve and the blue-violet hot spot above the black hole follows the X-ray one, so when the disc brightens, the moving dot in the light curve goes up. After Check, a real binary shows its second black hole going round at the injected period. The optical dots are spaced like ASAS-SN, about 390 points from 2012 to 2023 with a gap every year when the Sun is in the way. The X-ray dots are spaced like Swift-BAT, 157 monthly points from late 2004 to 2017.

The periodogram (Lomb 1976; Scargle 1982) fits a sine wave at every period from 100 days up to a third of the light curve and shows how much of the wiggle each one explains. The dashed line is a false-alarm line, set so only 1 in 100 noise-only curves goes above it anywhere. The textbook line assumes white noise, where every point is independent of the last. Quasar flicker wanders slowly instead, and slow wandering piles up power at long periods, which is where a binary would show up too. The red-noise line in round 3 is set by simulating that wandering, 1,000 curves per band. The grey stripes near 180 and 365 days are skipped, because the yearly gaps in the data make fake peaks there.

The real version ran on 1,369 active galaxies: 1,194 hard X-ray AGN from Swift-BAT and 175 from 4XMM-DR14, against optical light curves from ZTF and ASAS-SN. Both bands had to pass a red-noise line at 1 in 1,000 and agree on the period within 5%. None did. Noise alone predicted about 0.015 fake candidates across the BAT sample, and I got 0. Putting fake signals into the real light curves and running everything again says at most about 3% of the BAT galaxies can have a shared period between 100 and 900 days with a hard X-ray modulation of 30% or more (95% confidence). Below about 15% the limit says nothing, because monthly X-ray points can't pull out a signal that small.

What this toy leaves out

The false-alarm line here is 1 in 100, from 1,000 simulated curves and the tallest peak anywhere. The paper used 1 in 1,000, a statistic over 25 period bins and 2,000 simulations per light curve, with 100,000 for anything flagged in both bands. So the toy's line is 10 times looser and more fakes get over it.

The toy knows the true flicker timescale and size, because it made them. The paper had to fit them to every light curve, and for BAT that was messy: 89.9% of the raw fits fell back to a default value until the count rates were rescaled. The toy's X-ray points are also cleaner than real BAT data (0.08 mag of noise each by default). In the paper BAT got back at most about 18% of injected signals even at the largest amplitude, so if the toy's X-ray band finds a binary easily, the real one probably wouldn't have.

There's one optical band here. The paper used up to four (ZTF g and r, ASAS-SN V and g), kept the best and corrected for picking the best, and then ran a model-independent test (Robnik et al. 2024) on anything that got through. Nothing did. Both bands flicker independently here by construction, and for real AGN that's an assumption. Signal sizes are in magnitudes and the same in both bands, where the paper uses fractional flux amplitude (A = 0.3 is about 0.33 mag). The 3D quasar is an illustration: nothing is to scale and the orbit is sped up to match the playback.