{"id":"0d91df82-a3c4-4761-a56d-a2c1708346f0","arxiv_id":"2412.19573","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"WISE J0909+0002 shows a 660-689 day rest-frame quasi-periodic optical oscillation, interpreted as Doppler-boosted emission from a supermassive black hole binary.","lead":"Astronomers report a roughly 1.8-year quasi-periodic brightness cycle in the extremely luminous quasar WISE J0909+0002, based on 6.6 years of optical and UV light curves. The pattern points toward a pair of supermassive black holes whose orbital motion boosts the observed light, though red noise cannot be fully excluded.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"DRW false-alarm simulations use amplitude/timescale parameters inconsistent with the measured variability of WISE J0909+0002, so the quoted periodicity significance is not yet trustworthy.","rationale":"The reader identified exactly the same weakest assumption: DRW null parameters (tau 200–600 d, SF_inf 0.1–0.2 mag) do not match the measured variability (JAVELIN tau ~29 d for the combined curve; SF_inf 0.27–1.4 mag). My stress-test agrees and sharpens the issue: the paper explicitly excludes short-tau DRW simulations because they yield extreme amplitudes, but short tau is precisely what the combined light-curve JAVELIN posterior gives. That exclusion implicitly assumes the photometric errors and magnitude range of the actual data would veto such noise, yet those constraints are not folded into the quoted FAP. The BIC values in Table 4 are consistent with this: for DRW (400, 0.1) and (600, 0.1), 0.02–0.08% of mocks beat the sinusoidal model, but no run with a short-tau high-amplitude null is reported. Additional support for periodicity exists (consistent periods across independent bands, ACF, PSD showing SPL+periodic favored), so the claim is not empty; it is simply not yet calibrated under a null that resembles the data. Since the reader's verdict is already CONDITIONAL, my concern does not move the verdict; it reinforces the condition. The paper should either rerun the false-alarm simulation with matched noise parameters or explicitly quantify the effect of the tau > 200 d cut.","tokens_in":22035,"tokens_out":2576,"duration_ms":21678,"concrete_test":"Regenerate the §3.5 null simulations with DRW parameters matched to the data: draw tau and SF_inf from the JAVELIN posterior of the combined light curve (tau ~ 29 d, SF_inf ~ 0.27 mag), and if unconditional draws produce large amplitudes that break the photometric model, use the same conditional/stabilized generation procedure as astroML with clipping at the observed magnitude range or use a CARMA(1,1) or bending-power-law model instead. Count how many of 10,000 mock light curves yield a sinusoidal-fit BIC within 10 of the data (i.e., reproduce the observed quasi-periodic look). If the false-positive fraction rises above ~10%, the paper's quoted 0.02–0.08% FAP cannot be used to support the periodicity claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that the ~660–689 d rest-frame oscillation is real—rests on the false-alarm analysis of §3.5 contrasting the combined light curve against 10,000 SPL and DRW mock light curves. Table 4 restricts DRW parameters to tau = 200–600 d and SF_inf = 0.1 or 0.2 mag. But Table 3 shows the combined light curve has tau = 28.7 d and log SF_inf = -0.57 (SF_inf ~ 0.27 mag), and the three individual bands have log SF_inf between -0.05 and +0.15 (SF_inf ~ 0.9–1.4 mag). The text even states that sims with tau <~100 d were excluded because they produced 'extremely large amplitude (e.g., 10 mag)'. This is a red flag: the simulator evidently cannot reproduce the short-damping, high-amplitude noise that the data themselves imply. As a result the 0.02–0.08% DRW false-positive probabilities are computed under a null that is far smoother and longer-correlated than the observed noise, and the BIC contrast (Fig. 7 and Table 4) likely overstates how rare the observed quasi-sinusoid is. Because the same simulations are the only quantitative significance estimate (the bootstrap FAPs of §3.2 are not weighted by the DRW null), the headline significance can be inflated without any of the reported checks catching it. The PSD analysis (§3.6) partially mitigates this by fitting SPL and DRW+periodic models directly to the combined light curve, but that comparison is on the same data and does not calibrate the false-positive rate under the correct noise parameters.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports optical monitoring of the ELIRG/type-1 quasar WISE J0909+0002 at z=1.87 using archival CRTS, Pan-STARRS, and ZTF data supplemented by new OISTER/MITSuME/MuSaSHI observations. Sinusoidal fits, Lomb-Scargle periodograms, and a ZDCF autocorrelation analysis are used to identify a rest-frame period of roughly 660-689 days in the combined light curve and in the CRTS V, g, and r bands, spanning about 3.6 cycles. The authors then test the periodicity against single power-law and damped random walk red-noise simulations, examine structure functions and DRW parameters, perform PSD model selection, and compare observed amplitude ratios with spectroscopic power-law slopes to test the Doppler-boost scenario for a supermassive black hole binary. They also argue against a circumbinary disk and radio-jet precession, concluding that the periodic variability is likely caused by relativistic Doppler boosting in an SMBHB system.","tokens_in":22474,"tokens_out":13075,"duration_ms":113606,"significance":"If correct, this would be a notable SMBHB candidate: a z=1.87 ELIRG with a ~660-689 day rest-frame optical QPO and a 7.4e9 solar-mass black hole, with multi-band amplitude ratios at least partially consistent with Doppler boosting. The paper's strengths are that it brings new ground-based monitoring data into the period search, applies several independent period-finding diagnostics, and systematically tests competing mechanisms, including a radio-jet precession constraint from the FIRST nondetection. The amplitude-ratio test is not circular because it connects measured photometric amplitudes to independently measured spectral slopes. However, the central significance estimate depends on a red-noise null that is not calibrated to the object's measured variability parameters, and the multiband consistency is weakened by the excluded z and Ic bands and by the small number of cycles in individual bands; the evidence as presented is not yet strong enough to support the 'likely SMBHB' conclusion.","major_comments":[{"comment":"The DRW false-alarm simulations are run with tau = 200-600 d and SF_inf = 0.1-0.2 mag, whereas the JAVELIN fits in Table 3 yield tau = 28.7 d and SF_inf ~ 0.27 mag for the combined light curve and SF_inf ~ 0.9-1.4 mag with tau ~ 350-435 d for the three bands. The statement that simulated DRW curves with tau <~ 100 d produce 'extremely large amplitude (e.g., 10 mag)' indicates that the simulation pipeline is being run outside the parameter regime that describes the data, not that the measured short damping timescale is unphysical. Because the quoted 0.02%-0.08% false-positive probabilities and the Delta-BIC distributions in Figure 7 are conditional on this unrepresentative null, they cannot yet be used as the quantitative significance of the 660-689 d periodicity. I request a reanalysis with CAR(1) parameters drawn from the Table 3 posteriors, including the short-tau combined solution and the per-band SF_inf values, with the resulting false-alarm rates reported; if the simulator cannot realize those parameters, the simulation procedure must be corrected or the significance claim removed.","section":"Section 3.5, Tables 3 and 4"},{"comment":"The presentation of the periodicity as confirmed by 'three-band light curves' is weakened by the small number of cycles and by post-hoc band selection. In the rest frame, the CRTS V-band data span only ~1010 d, i.e., ~1.5 cycles at 654 d; the g-band light curve spans ~1805 d but contains a ~600 d gap and only ~2.7 cycles; and the r-band spans ~1430 d, i.e., ~2.2 cycles. Only the combined light curve reaches the ~3.6 cycles quoted in the text, and that curve is built by rescaling and stitching together different filters. I ask the authors to report per-band cycle counts and gap statistics and to fit a common-period model with per-band amplitudes and phases; as written, the independent multiband confirmation is overstated.","section":"Section 3.2, Table 2, Figures 2-4"},{"comment":"The z- and Ic-band light curves are excluded from the periodicity claim after the fact, yet they are part of the same monitoring program and cover the same source. Their best-fit periods, 1232.69 +/- 74.26 d and 378.56 +/- 40.58 d, are inconsistent with the proposed 660-689 d period, and the text notes that these bands show 'very different trends from the other curves.' Since the abstract and Section 5 describe a multiband quasi-periodic oscillation, the paper should demonstrate quantitatively whether z and Ic are compatible with the common period given their sampling and photometric noise. Without such a test, the possibility that the selected bands are a chance subset is not addressed.","section":"Section 3.1, Table 2"},{"comment":"The Doppler-boost test is only marginally consistent: three of six amplitude-ratio combinations agree with equation (12) within 1 sigma, and the spectral slopes from the two SDSS epochs differ substantially (e.g., the r-band alpha_nu is -0.966 +/- 0.017 at MJD 55532 versus -0.111 +/- 0.021 at MJD 51929), so the agreement depends on which epoch is adopted. The parameter-space plot in Figure 10 also assumes cos(phi) = 1, f2 = 0.8, and the single measured black-hole mass, and it yields i ~ 10 degrees only for a restricted q range, which the paper itself notes is in tension with the usual BAL orientation interpretation. I recommend presenting this as a weak consistency check rather than as 'likely' evidence for Doppler boosting unless additional multi-wavelength amplitude ratios are added.","section":"Section 4.1.1, Figure 9"}],"minor_comments":[{"comment":"The text contains the typo 'BHSH theory' in the discussion of PG 1302-102; this should be 'BSBH.'","section":"Section 1.1"},{"comment":"Section 3.2 says the periodic variability continued for ~6.6 yr in the quasar rest frame, corresponding to ~3.6 cycles, but conclusion (2) of Section 5 says 'at least ~3.6 yr in the rest frame'; the latter should be '6.6 yr' or '3.6 cycles.'","section":"Sections 3.2 and 5"},{"comment":"The phrases 'signal-to-noise ratio of >~ 10' and 'removed those data with a wrong weather condition' need editorial correction for clarity.","section":"Section 2.2"},{"comment":"The caption says all panels include Pan-STARRS, ZTF, and 'our observation data,' but panel (b) shows the CRTS V-band light curve; please clarify which data sets contribute to each panel.","section":"Figure 4 caption"}],"recommendation":"major_revision","confidential_remarks":"The central issue is the DRW simulation calibration, which I believe is fixable by reanalysis with parameters drawn from the measured JAVELIN posteriors and by a common-period test that includes the z and Ic bands. The observational dataset and the scenario tests are valuable, so I would not reject the manuscript; however, the 'likely SMBHB' conclusion should not be advertised until the revised false-alarm rates are available. Please also have the authors correct the Section 5 '3.6 yr' wording."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"WISE J0909+0002 is a new, well-observed periodic quasar candidate, but its significance test is built on a null that does not match the data. The paper does something genuinely new: it adds an ELIRG to the periodic quasar zoo, with ~660–689 day rest-frame oscillations, and it brings fresh monitoring data from OISTER plus archival CRTS/Pan-STARRS/ZTF light curves. The analysis is thorough: Lomb-Scargle, ACF, PSD, structure functions, Doppler boost amplitude ratios, a circumbinary disk SED check, and a radio limit ruling out jet precession. The authors are honest about the limits, explicitly saying pure red noise cannot be fully rejected.\n\nThe soft spot is the false-alarm analysis in §3.5. The DRW simulations use tau=200–600 d and SF_inf=0.1–0.2 mag, while the measured JAVELIN parameters for the combined light curve are tau~29 d and SF_inf~0.27 mag; the individual bands have SF_inf~0.9–1.4 mag. The paper even notes that simulated light curves with tau<~100 d produced amplitudes of 10 mag, which tells you the simulator cannot reproduce the noise properties of the actual data. So the quoted 0.02–0.08% false-positive probabilities are computed under a null that is smoother and longer-correlated than what the data show. The periodicity significance is likely overstated. The detection itself rests on 3.6 cycles, and the band selection (V, g, r) is partly post-hoc; the z and Ic bands show different trends.\n\nThe Doppler boost interpretation is plausible but weak. Only three of six amplitude-ratio combinations are consistent with the model within 1σ, and the inclination angle implied by the model (i~10°) conflicts with the expectation that BAL quasars are viewed closer to edge-on. The authors discuss this, but it does not strengthen the case.\n\nOverall, this is a legitimate candidate report with solid data and an honest discussion. The central claim, however, is not yet trustworthy because the null simulations are mis-specified. It deserves a serious referee: the fix is clear (run the simulations with the measured tau and SF_inf, or use a more flexible noise model), and the object is worth monitoring. I would not cite it as evidence for a binary SMBH until that is done. For a reading group, it is a useful example of how a significance claim can depend on the assumed noise.","headline":"A genuinely new periodic quasar candidate in an ELIRG, but the claimed significance rests on DRW simulations that do not match the measured noise.","tokens_in":23027,"tokens_out":3080,"would_cite":false,"duration_ms":45933,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A distant quasar repeats its brightness every ~660 to 689 days, hinting at a binary supermassive black hole.","keywords":["quasi-periodic oscillations","supermassive black hole binary","Doppler boost","WISE J090924.01+000211.1","extremely luminous infrared galaxy","damped random walk","AGN variability","Lomb-Scargle periodogram"],"falsifier":"Regenerate the 10,000 mock damped-random-walk light curves with the measured noise parameters ($SF_\\infty \\approx 0.27$ to $1.40$ mag, $\\tau\\simeq29$ day) and count how often they match the best-fit sinusoid; if the false-positive fraction exceeds a few percent, the quasi-periodic signal is not distinguishable from red noise. Continued monitoring that shows the ~660–689 day phase drifting or the amplitude fading would also falsify the binary and Doppler-boost interpretation.","tokens_in":2287,"feed_emoji":"🔭","tokens_out":3819,"duration_ms":113288,"temperature":0.7,"pith_summary":"This paper claims that WISE J0909+0002, an extremely luminous infrared quasar at redshift 1.87, has been oscillating in brightness with a quasi-sinusoidal period of roughly 660 to 689 days in the quasar's own frame, sustained for about 3.6 cycles. The periodicity appears in archival CRTS, Pan-STARRS, and ZTF light curves and in new ground-based monitoring, and signal-to-noise, Lomb-Scargle, autocorrelation, and power-spectrum analyses all point to the same period. The authors argue this is not ordinary quasar flickering: the structure functions do not follow the damped random walk that describes most AGN variability, and mock damped-random-walk light curves produce such a clean sinusoid only rarely. They conclude that the most likely cause is relativistic Doppler boosting of emission from a supermassive black hole binary, and they rule out circumbinary disk features and radio jet precession as alternatives. If right, this object becomes one of the clearest periodic-quasar candidates found inside an extremely luminous galaxy merger, a rare place to look for binary black hole signatures.","feed_headline":"Quasar pulses every ~660 days for 6.6 years","feed_subtitle":"The repeating brightness pattern points to a supermassive black hole binary, if it is not red noise.","key_machinery":"The argument rides on the Doppler boost amplitude-ratio identity and on a carefully merged multi-band light curve. For a power-law spectrum $F_\\nu\\propto\\nu^{\\alpha_\\nu}$, Doppler boosting gives $\\Delta F_\\nu/F_\\nu = (3-\\alpha_\\nu)\\beta\\cos\\varphi\\sin i$, so the ratio of variability amplitudes between two bands equals $(3-\\alpha_s)/(3-\\alpha_l)$; this identity is what separates the boost scenario from intrinsic disk variability. Around it the paper wraps a sinusoidal model $F(t)=A\\sin(2\\pi t/P+\\varphi)+b$ with a signal-to-noise criterion $\\xi=A^2/(2\\sigma_{\\rm res}^2)$, Lomb-Scargle periodograms with bootstrap false-alarm probabilities, a z-transformed discrete correlation function for the autocorrelation, JAVELIN CAR(1) fits and structure functions to test the damped random walk null, and 10,000 simulated SPL and DRW light curves compared by BIC. The same data and SED are then used to test and reject the circumbinary disk and jet precession alternatives.","core_discovery":"The central claim is that the combined optical/UV light curve of WISE J0909+0002, assembled from CRTS V-band, Pan-STARRS, ZTF, and new ground-based observations, is quasi-periodic with a rest-frame period of 660–689 days, significant at a bootstrap false-alarm probability below 0.001, with a z-transformed autocorrelation period of 659.6 days and a power spectrum that favors a power-law-plus-periodic model. The paper identifies the cause as a relativistic Doppler boost: in the model $\\Delta F_\\nu/F_\\nu = (3-\\alpha_\\nu)\\beta\\cos\\varphi\\sin i$, the measured ratios of variability amplitudes between bands agree with the ratios of $(3-\\alpha_\\nu)$ predicted from SDSS spectral slopes in three of six band-pair cases, and a viable parameter space allows inclination near $10^\\circ$ for mass ratios $q=0.11$ to $0.43$. The circumbinary disk model is disfavored because the spectral energy distribution shows no predicted UV cut-off or notch, and radio jet precession is excluded by a FIRST upper limit implying radio loudness $R\\leq0.4$. The paper therefore proposes that WISE J0909+0002 is a supermassive black hole binary candidate, while acknowledging that the damped random walk null is not completely rejected because rare mock DRW curves can mimic the signal.","pith_inferences":["A direct robustness test the paper does not carry out: rerun the 10,000 mock DRW simulations with the measured noise parameters ($SF_\\infty \\approx 0.27$ to $1.40$ mag, $\\tau\\simeq29$ day for the combined curve) instead of the adopted $SF_\\infty=0.1$ to $0.2$ mag and $\\tau=200$ to $600$ day; the reported 0.02% to 0.08% false-positive rates could rise substantially.","A Doppler-boosted binary predicts that all bands oscillate in phase with amplitude ratios fixed by the spectral slope; simultaneous multi-band photometry across one full observed cycle (~1900 days) could discriminate this from a single-disk instability producing a similar period.","If confirmed, the ~660 day rest-frame period places the pair at sub-parsec separation, offering a rare environmental test of merger-driven fueling in hyperluminous infrared systems, which the paper leaves implicit.","The z- and Ic-band light curves deviate from the common sinusoid; more frequent, well-sampled infrared photometry could determine whether these deviations are real band-dependent behavior or artifacts of sparse epochs, which matters for the Doppler boost interpretation."],"forward_implications":["If the periodicity is real, WISE J0909+0002 is a supermassive black hole binary candidate with an orbital period near 660–689 days, one of the few found in an extremely luminous infrared galaxy.","The Doppler boost explanation predicts that variability amplitude should scale with wavelength through the $(3-\\alpha_\\nu)$ factor; future simultaneous UV, optical, and infrared monitoring can confirm or reject this within a few cycles.","The absence of damped-random-walk structure in such a luminous quasar would mark it as a variability outlier, strengthening the view that ELIRGs are in a distinct, merger-driven accretion phase.","Continued monitoring should show the same phase and amplitude persisting; if they do, the QPO becomes a robust binary candidate suitable for gravitational-wave follow-up.","If a binary is confirmed, the system joins the small set of periodic quasars that can be used to probe sub-parsec supermassive black hole pairs and their merger-driven fueling."],"supporting_citations":[{"why":"Establishes the archetypal periodic quasar PG 1302-102 and motivates searches for quasi-periodic oscillations in AGN.","marker":"Graham et al. 2015b"},{"why":"Supplies the relativistic Doppler boost model and the variability-amplitude equation that the paper tests against multiband light curves.","marker":"D'Orazio et al. 2015b"},{"why":"Provides the search methodology, detection rates, and false-positive benchmarks for periodic quasars.","marker":"Charisi et al. 2016"},{"why":"Gives the damped random walk variability baseline and structure-function comparison values used to argue the light curves are not DRW.","marker":"MacLeod et al. 2010"},{"why":"Defines the CAR(1)/damped random walk stochastic model that JAVELIN fits to the light curves.","marker":"Kelly, Bechtold, and Siemiginowska 2009"},{"why":"Supplies the SNR criterion, BIC threshold, and Doppler boost amplitude-ratio test applied here.","marker":"Chen et al. 2020"},{"why":"Provides a comparison periodic quasar whose Doppler boost and circumbinary disk diagnostics are directly compared with WISE J0909+0002.","marker":"Liao et al. 2021"},{"why":"Identifies WISE J0909+0002 as an ELIRG and provides the infrared luminosity, black hole mass, and SED used in all model tests.","marker":"Toba et al. 2021"},{"why":"Provides the Lomb-Scargle implementation and false-alarm treatment used to establish the ~660–689 day periodicity.","marker":"VanderPlas 2018"},{"why":"Supplies the JAVELIN code used to estimate the damped random walk parameters and their uncertainties.","marker":"Zu, Kochanek, and Peterson 2011"}],"fun_headline_variants":["Quasar's 660-day cycle hints at black hole binary","Distant quasar's periodic flicker points to binary black hole","Quasar brightness rhythm may reveal supermassive pair","660-day quasar pulses: Doppler boost or just noise?","ELIRG quasar shows quasi-periodic 660-day variability"],"cache_read_input_tokens":24960,"weakest_assumption_plain":"The load-bearing premise is that the red-noise baseline used to judge the periodicity—a damped random walk with $SF_\\infty = 0.1$ to $0.2$ mag and $\\tau = 200$ to $600$ day—is representative of this object's noise; the measured $SF_\\infty$ (up to about 1.4 mag) and combined-light-curve $\\tau\\simeq29$ day suggest the real noise is stronger and faster, so the reported false-positive probability of 0.02% to 0.08% could be optimistic.","fun_headline_variants_meta":{"raw":{"variants":["Quasar's 660-day cycle hints at black hole binary","Distant quasar's periodic flicker points to binary black hole","Quasar brightness rhythm may reveal supermassive pair","660-day quasar pulses: Doppler boost or just noise?","ELIRG quasar shows quasi-periodic 660-day variability"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000248,"raw_usage":{"total_tokens":1704,"prompt_tokens":1258,"completion_tokens":446,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":874,"completion_tokens_details":{"reasoning_tokens":360}},"tokens_in":874,"tokens_out":446,"duration_ms":5357,"temperature":1.0,"reasoning_tokens":360,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T00:11:58.702917+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Regenerate the 10,000 mock damped-random-walk light curves with the measured noise parameters ($SF_\\infty \\approx 0.27$ to $1.40$ mag, $\\tau\\simeq29$ day) and count how often they match the best-fit sinusoid; if the false-positive fraction exceeds a few percent, the quasi-periodic signal is not distinguishable from red noise. Continued monitoring that shows the ~660–689 day phase drifting or the amplitude fading would also falsify the binary and Doppler-boost interpretation.","supporting_citations":[{"cited_title":"L., et al","cited_arxiv_id":null,"evidence_quote":"Gives the damped random walk variability baseline and structure-function comparison values used to argue the light curves are not DRW."},{"cited_title":"2020, MNRAS, 499, 2245","cited_arxiv_id":null,"evidence_quote":"Supplies the SNR criterion, BIC threshold, and Doppler boost amplitude-ratio test applied here."},{"cited_title":"2021, MNRAS, 500, 4025","cited_arxiv_id":null,"evidence_quote":"Provides a comparison periodic quasar whose Doppler boost and circumbinary disk diagnostics are directly compared with WISE J0909+0002."},{"cited_title":"2021, A&A, 649, L11","cited_arxiv_id":null,"evidence_quote":"Identifies WISE J0909+0002 as an ELIRG and provides the infrared luminosity, black hole mass, and SED used in all model tests."}],"review_version":1}