{"id":"f28d2985-2701-4865-9d63-6e324b66bcab","arxiv_id":"1908.04280","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"WISE J030654.88+010833.6 is a rare mid-infrared flaring AGN in a ULIRG candidate at redshift about 0.3, with variability likely from accretion-rate changes reprocessed by dust.","lead":"We found a distant dusty galaxy whose infrared glow roughly doubled and then faded over eight years, caused by a supermassive black hole at its center that is feeding on gas. This rare flaring galactic nucleus inside an ultra-luminous infrared galaxy candidate gives astronomers a new testbed for how black holes and their host galaxies grow together.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"ULIRG classification and BH mass rest on an unconstrained Mrk231 template extrapolation; the existing IRAS non-detection is used only as a rough SFR consistency check, not as a quantitative template test.","rationale":"The paper's core discovery—a WISE-selected, MIR-variable, red extended galaxy with radio and X-ray detections—is credible. The W1 and W2 light curves show correlated factor-of-two variability over 8 years with quoted per-epoch errors of about 0.03 mag, and the source has independent AGN indicators: W1-W2 greater than 0.8, NVSS radio luminosity above 10^24 W/Hz, and Chandra X-ray luminosity around 2e43 erg/s. I do not see a reason to doubt the flaring-AGN identification. The load-bearing weakness is the translation of that detection into physical classification. The ULIRG label, the SFR_IR of about 70 M_sun/yr, and the BH mass upper limit all come from a single scaling of the Mrk231 template to MIR data, with no FIR anchor. The paper is transparent about this, explicitly calling the SFR_IR and BH mass estimates highly uncertain and noting the absence of FIR photometry or spectra. It also carefully labels the source a ULIRG candidate rather than a confirmed ULIRG. Still, the IRAS non-detection is already in hand and could either support or falsify the template extrapolation if the predicted 60/100 micron fluxes were compared quantitatively with the IRAS noise. The paper's consistency check using a 1-sigma SFR equivalent is too coarse to settle that. Therefore the reader's conditional verdict is appropriate, and the concern is best resolved by a quantitative IRAS-based test or by future FIR observations.","tokens_in":14282,"tokens_out":11607,"duration_ms":126643,"concrete_test":"Use SCANPI or IRIS/IRAS 60 and 100 micron maps to measure or upper-limit the source flux at the position of WISE-J0301. Predict the 60/100 micron flux densities from the best-fit Mrk231 template (scaled 0.45, z=0.28) by convolving the template SED with the IRAS bandpasses. If the predicted combined 60+100 micron flux exceeds roughly 3-sigma of the measured noise, the template extrapolation is internally inconsistent and the LIR range in Table 3 cannot be maintained. If the prediction is below 1-sigma, the template is consistent with all existing data but still unconstrained, and the conditional verdict stands pending actual FIR observations.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 5.3 derives L_IR = 9.0e11–1.5e12 L_sun and SFR_IR ~70 M_sun/yr, and Section 5.4 derives M_BH <= 1e8 M_sun, all by scaling the Mrk231 template by 0.45 to the MIR photometry. The fit uses only WISE W1–W4 and Spitzer IRAC data: the rest-frame 8–1000 micron SED, particularly the FIR peak, is assumed from Mrk231 with no FIR photometry or spectrum. Because the quoted L_IR range straddles the 10^12 L_sun ULIRG boundary, a modest change in FIR shape or normalization removes the ULIRG classification. The IRAS non-detection is not a validation: the paper quotes a 1-sigma FIR limit equivalent to SFR 67 M_sun/yr and calls it broadly consistent with SFR_IR, but this compares a star-formation-only conversion of an upper limit to a template-derived total-infrared SFR, and the predicted 60/100 micron fluxes from the scaled template are never compared to the actual IRAS noise. The photo-z tension (SDSS KD-tree z=0.215 vs adopted 0.28–0.31) compounds the issue, since a lower z also lowers L_IR. The variability and AGN identification themselves do not depend on the template and are well supported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports the discovery and multiwavelength characterization of WISE J030654.88+010833.6, a WISE-selected, red extended galaxy in Stripe 82. From WISE and NEOWISE single-exposure photometry, the authors construct 8-year W1/W2 light curves that show correlated factor-of-~2 variability at 3.4 and 4.6 micron, with per-epoch photometric uncertainties at or below 0.03 mag. They assemble GALEX, SDSS, PanSTARRS, 2MASS, WISE, Spitzer IRAC, NVSS, and archival Chandra data. FSPS/EAZY SED fitting gives a photometric redshift z ~ 0.28-0.31 and a stellar mass near 1.9e11 Msun; scaling the Mrk231 template to the MIR photometry yields L_IR ~ 9e11-1.5e12 Lsun, SFR_IR ~ 70 Msun/yr, and a BH mass range ~1e7-1e8 Msun. The paper interprets the MIR variability as accretion-rate variations whose energy is reprocessed by dust, with a possible lag to be tested by future Spitzer and PanSTARRS data.","tokens_in":14562,"tokens_out":9758,"duration_ms":100454,"significance":"If the variability detection and AGN identification hold, WISE-J0301 is a scientifically valuable object: a mid-IR variable, radio- and X-ray-detected AGN at z ~ 0.3 whose host may be a ULIRG. The variability claim is well supported by the light-curve quality cuts, the small photometric errors, and the correlated behavior between W1 and W2, and the radio and X-ray detections independently support AGN activity. The analysis uses publicly available archival data and standard external templates, so the variability claim is not circular. The main weakness is that the ULIRG classification, SFR, and BH mass depend on a scaled Mrk231 template without FIR photometry; those parts of the paper need to be either quantitatively supported or explicitly downgraded.","major_comments":[{"comment":"The ULIRG classification and SFR_IR are derived by scaling the Mrk231 template by 0.45 to the MIR photometry, with no FIR detection (IRAS non-detection only; no Akari or Herschel data). The quoted L_IR = 9.0e11-1.5e12 Lsun straddles the 10^12 Lsun ULIRG boundary, so a modest change in FIR shape or normalization changes the classification. The IRAS consistency check does not validate the template: the paper quotes a 1-sigma 60+100 micron limit converted to SFR ~67 Msun/yr through a star-formation conversion, but it never compares the scaled template's predicted 60 and 100 micron flux densities with the per-band SCANPI noise. I request that quantitative comparison, or a clear statement in the title and conclusions that the ULIRG nature is unconfirmed. In addition, the template scaling should be applied to AGN-only W1/W2 fluxes after subtracting the stellar contribution estimated in Section 3.1, rather than to the total WISE photometry, and the sensitivity to this choice should be reported.","section":"Section 5.3, Figure 3"},{"comment":"The adopted redshift z = 0.28-0.31 is load-bearing for L_IR, because the input optical photometry from SDSS and PanSTARRS was taken at different epochs and is variable (Table 1), and because the SDSS KD-tree photo-z is z = 0.215. Since L_IR scales as the square of the luminosity distance, adopting z = 0.215 would lower L_IR by roughly a factor of two, which would move the source below the ULIRG threshold. Please provide the SED fit at z = 0.215 and the resulting L_IR, SFR_IR, and BH mass, and discuss how the epoch-dependent photometry affects the photo-z preference.","section":"Section 4, Table 2; Section 5.3"},{"comment":"The BH mass upper limit of ~1e8 Msun is obtained by multiplying the Mrk231 BH mass by the same 0.45 luminosity scaling factor used in the MIR template fit, which assumes that the Eddington ratio and the MIR-to-bolometric correction of WISE-J0301 are identical to those of Mrk231; the X-ray lower limit assumes Eddington-limited accretion. These are assumptions rather than measurements, so the statement that the BH-to-stellar mass ratio is consistent with broad-line AGN is a consistency check of the template, not an independent constraint. I recommend presenting the BH mass as an illustrative range under explicit assumptions, as Section 6 does, and qualifying the corresponding sentence in the abstract.","section":"Section 5.4, Eq. (6)"}],"minor_comments":[{"comment":"The wavelength entry '2.14 x 10^5' for NVSS should be labeled as 21 cm (2.14e5 micron) so that the column units are unambiguous.","section":"Table 1, NVSS row"},{"comment":"The text says '9 bands of GALEX+SDSS+2MASS, and 10 bands of GALEX+SDSS u+Pan-STARRS+2MASS,' but the listed surveys contain 10 and 11 bands, respectively; please clarify which bands were excluded.","section":"Section 4, Table 2"},{"comment":"The Spitzer aperture correction '1.047' should state whether it is a multiplicative flux factor or an additive magnitude correction.","section":"Section 3.1"},{"comment":"The statement that PanSTARRS DR2 shows 'a factor of 3 variability' in the y-band is difficult to reconcile with the 0.25-0.39 mag variability quoted for i/r/z; please correct or clarify this value.","section":"Section 3.2"},{"comment":"The Introduction says the Mrk231 fit 'identify the host galaxy as a ULIRG,' while Sections 5.3 and 6 use 'possibly' and 'highly uncertain'; the wording should be harmonized with the actual uncertainty so the body and abstract do not overstate the classification.","section":"Introduction / Section 5.3"}],"recommendation":"major_revision","confidential_remarks":"The core variability detection and the AGN identification are solid and should survive revision. The main risk is that the ULIRG and BH-mass claims are presented more strongly than the available constraints support; the authors themselves acknowledge this in Section 6. I would be satisfied after either a quantitative IRAS test of the Mrk231 template or an explicit downgrade of these claims. The stress-test concern about the IRAS non-detection is valid: it is used as a consistency check, not as a template validation, and the per-band predicted fluxes are never compared with the IRAS noise."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi [Name], here's my read on Prakash et al. 1908.04280.\n\nThe paper's core claim is credible: WISE-J0301 is a strongly MIR-variable AGN (factor ~2 in W1/W2 over 8 years) in a red, extended galaxy at z~0.3, with supporting radio and X-ray detections. That is the real payoff, and it is well supported. The WISE single-exposure photometry is careful—quality cuts, aperture corrections, stellar contribution subtracted—and the correlated variability and steep spectral index point to dust-reprocessed accretion variability. The comparison to Mrk231 is honest, including the radio excess.\n\nWhat is genuinely new is the object itself: a rare mid-IR variable AGN hosted in a likely ULIRG. It is a case study, not a new mechanism, and the paper doesn't overclaim that.\n\nThe soft spots are where the stress-test lands. The ULIRG classification, SFR~70 M_sun/yr, and BH mass all come from scaling the Mrk231 template by 0.45 to the MIR photometry. There is no FIR detection and no spectrum. The IRAS non-detection is used as a consistency check, but that check is not quantitative: the predicted 60/100 um fluxes from the scaled template are never compared to the IRAS noise, and the quoted SFR limit assumes a star-formation-only conversion of an upper limit. A modest change in the FIR shape or normalization moves L_IR across the 10^12 L_sun ULIRG boundary. The adopted photo-z (0.28–0.31) is also not bulletproof—SDSS KD-tree says 0.215—and a lower z would lower L_IR. The paper acknowledges most of this, which I credit, but the uncertainties are not propagated into the tabulated physical parameters.\n\nStill, the central variability result does not depend on any of that. The AGN identification is solid from radio luminosity and X-ray detection, and the MIR variability is direct. So the discovery stands as a useful addition to the small population of such objects. The paper is well-cited and the authors are appropriately careful.\n\nWho should read this: people working on AGN variability, dusty torus reprocessing, and the AGN-starburst connection. It is a good referee candidate: the object is worth publishing, but the referee should push for a more rigorous treatment of the FIR template uncertainty, perhaps alternative SED fits and a direct comparison of the template predictions to the IRAS limits.\n\nMy recommendation: send to peer review. The central claim is strong enough to warrant referee time, and the specific concerns are addressable in revision.","headline":"A credible discovery of a rare MIR-variable AGN in a ULIRG candidate; the variability is solid, but the ULIRG classification and BH mass rest on a template extrapolation that needs better scrutiny.","tokens_in":15168,"tokens_out":2460,"would_cite":true,"duration_ms":23480,"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":"Mid-infrared flux doubles in a Stripe 82 galaxy, marking a flaring AGN.","keywords":["mid-infrared variability","active galactic nuclei","ultraluminous infrared galaxies","WISE","accretion rate","dust reprocessing","Stripe 82","SED fitting"],"falsifier":"A single far-infrared measurement of WISE-J0301 at, say, 160–500 microns would settle the ULIRG question: if the observed flux is inconsistent with the scaled Mrk231 template and the integrated $L_{\\rm IR}$ falls below $10^{12}\\,L_\\odot$, the ULIRG classification fails, while a measured optical-to-mid-infrared lag would test the accretion-rate reprocessing scenario.","tokens_in":14070,"feed_emoji":"🔭","tokens_out":18220,"duration_ms":143515,"temperature":0.7,"pith_summary":"This paper identifies WISE-J0301, a red extended galaxy in the Stripe 82 field, as one of the few known mid-infrared variable active galactic nuclei hosted in an ultraluminous infrared galaxy (ULIRG) candidate. Using WISE single-exposure images, the authors show that its 3.4 and 4.6 micron flux varied by roughly a factor of two over eight years, with correlated changes in both bands and no strong optical variability in the bluer bands. A stellar-population fit to the ultraviolet, optical, and near-infrared photometry gives a photometric redshift of $0.28\\le z\\le0.31$ and a stellar mass near $2\\times10^{11}\\,M_\\odot$. Scaling the spectral template of Mrk231, a nearby prototypical obscured AGN, to the mid-infrared photometry yields a total infrared luminosity of about $9\\times10^{11}$ to $1.5\\times10^{12}\\,L_\\odot$, placing the galaxy at the ULIRG boundary, with a star-formation rate near $70\\,M_\\odot\\,\\mathrm{yr}^{-1}$ and a black hole mass upper limit near $10^8\\,M_\\odot$. The paper's central physical claim is that the long-term mid-infrared variability is caused by changes in the black hole's accretion rate, with the energy reprocessed by dust near the AGN.","feed_headline":"Stripe 82 galaxy's mid-infrared flux doubles, marking a flaring AGN","feed_subtitle":"Eight years of WISE data suggest accretion-rate changes, with the emission reprocessed by dust near the black hole.","key_machinery":"The load-bearing mechanism is the spectral template of Mrk231, a nearby prototypical obscured AGN, scaled by a factor of 0.45 to the observed WISE mid-infrared photometry. This template converts the measured 3.4 and 4.6 micron fluxes into an estimate of the total 8–1000 micron luminosity, from which the ULIRG classification, the star-formation rate of roughly $70\\,M_\\odot\\,\\mathrm{yr}^{-1}$, and the upper limit on the black hole mass all follow. The companion mechanism is the construction of the eight-year WISE light curves from quality-filtered single-exposure images, with the stellar contribution subtracted using the SED model; that procedure isolates a factor-of-two variability that is correlated between W1 and W2 and absent in the bluest optical bands.","core_discovery":"The paper's discovery is that WISE J030654.88+010833.6 (WISE-J0301), a red extended galaxy at photometric redshift $0.28\\le z\\le0.31$, is a mid-infrared variable AGN hosted in an ultraluminous infrared galaxy candidate. In the WISE W1 and W2 light curves spanning eight years, the source brightened and faded by roughly a factor of two, with the changes correlated between the two bands, while the bluest optical light showed no comparable variability. The ultraviolet-to-near-infrared SED is well fit by a stellar population, giving $M_*\\sim1.9\\times10^{11}\\,M_\\odot$ and an unobscured star-formation rate of about $0.7\\,M_\\odot\\,\\mathrm{yr}^{-1}$; the mid-infrared photometry is fit by a scaled Mrk231 template (factor 0.45), implying $L_{\\rm IR}$ between $9.0\\times10^{11}$ and $1.5\\times10^{12}\\,L_\\odot$, a reprocessed star-formation rate near $70\\,M_\\odot\\,\\mathrm{yr}^{-1}$, and a black hole mass upper limit near $10^8\\,M_\\odot$. The paper interprets the correlated long-term mid-infrared variability as variations in the accretion rate, with the luminosity thermally reprocessed by dust near the AGN, and notes that the host may be at a very late stage of a merger.","pith_inferences":["The ULIRG classification depends entirely on the Mrk231 scaling, so a single far-infrared measurement that falsifies the template would leave the variability result intact but move the source below $10^{12}\\,L_\\odot$; the paper's physical interpretation would then need revision.","The observed range of W1–W2 spectral index, from about $-2.7$ to $-1.2$, may reflect a varying mix of hot dust and 3.3 $\\mu$m PAH emission; decomposing the color variability by epoch could separate dust reprocessing from intrinsic changes in the accretion disk.","Because the y-band minimum appears to fall in the gap of the WISE light curve, the currently claimed lack of optical-mid-infrared correlation is provisional; denser sampling across a future outburst could either confirm the expected lag or rule out the thermal-reprocessing scenario.","If the host is confirmed as a ULIRG at $z\\sim0.3$, it becomes a strong target for molecular-gas mapping to test whether AGN feedback is clearing the nucleus and whether the AGN dominates the total infrared output."],"forward_implications":["The variability detection is independent of the SED template: the correlated factor-of-two change in W1 and W2 over eight years stands on the photometry alone, so the object is a robust mid-infrared variable AGN even if the ULIRG classification is later revised.","If the accretion-rate interpretation is correct, simultaneous optical and mid-infrared monitoring should reveal a time delay set by the light-travel time from the accretion disk to the hot dust, a test that planned infrared and optical monitoring can carry out.","With $L_{\\rm IR}$ near $10^{12}\\,L_\\odot$, the source joins the rare low-redshift class of AGN-dominated ULIRGs, and its black-hole-to-stellar-mass ratio becomes consistent with local broad-line AGN, supporting the Mrk231-like SED.","A confirmed ULIRG at $z\\approx0.3$ with no obvious merger morphology would indicate either a very late-stage merger or a non-merger trigger for the simultaneous starburst and AGN activity."],"supporting_citations":[{"why":"Supplies the Mrk231 template that is scaled to the mid-infrared photometry; this template drives the L_IR, SFR_IR, and black-hole mass estimates.","marker":"Chary & Elbaz (2001)"},{"why":"Provides the original Mrk231 spectral energy distribution used as the prototype obscured-AGN template.","marker":"Sanders et al. (1988)"},{"why":"Provides the template-fitting method used to derive the photometric redshift and stellar population parameters from the UV-to-NIR photometry.","marker":"Brammer et al. (2008)"},{"why":"Supplies the stellar population synthesis models underlying the galaxy templates fit to the UV-NIR SED.","marker":"Conroy et al. (2009)"},{"why":"Establishes the rarity of mid-infrared variable AGN, providing the comparison population of 162 objects with radio detections.","marker":"Assef et al. (2018)"},{"why":"Provides the W1-W2 color selection criterion used to choose AGN candidates in the WISE catalog.","marker":"Stern et al. (2012)"},{"why":"Gives the mid-infrared to X-ray luminosity correlation used to correct the X-ray luminosity and set a lower limit on black hole mass.","marker":"Lutz et al. (2004)"},{"why":"Supplies the Mrk231 black hole mass used as the scaling reference for the WISE-J0301 mass upper limit.","marker":"Leighly et al. (2014)"},{"why":"Provides the broad-line AGN black-hole-to-stellar-mass relation used to argue that the $10^8\\,M_\\odot$ estimate is plausible.","marker":"Reines & Volonteri (2015)"}],"fun_headline_variants":["Variable AGN found hiding in ultraluminous galaxy","WISE spots flaring AGN in dusty ULIRG candidate","Mid-infrared clues reveal flaring black hole in Stripe 82","Galaxy's 8-year flicker points to active black hole","Flaring AGN detected in distant ULIRG candidate"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's physical parameters—ULIRG classification, star-formation rate, and black hole mass upper limit—all rest on the assumption that the Mrk231 spectral template, scaled by 0.45, accurately represents the unseen rest-frame 8–1000 micron SED of WISE-J0301, since no far-infrared detection or spectrum exists.","fun_headline_variants_meta":{"raw":{"variants":["Variable AGN found hiding in ultraluminous galaxy","WISE spots flaring AGN in dusty ULIRG candidate","Mid-infrared clues reveal flaring black hole in Stripe 82","Galaxy's 8-year flicker points to active black hole","Flaring AGN detected in distant ULIRG candidate"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000832,"raw_usage":{"total_tokens":3727,"prompt_tokens":1134,"completion_tokens":2593,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":750,"completion_tokens_details":{"reasoning_tokens":2507}},"tokens_in":750,"tokens_out":2593,"duration_ms":18826,"temperature":1.0,"reasoning_tokens":2507,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:47:00.551381+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A single far-infrared measurement of WISE-J0301 at, say, 160–500 microns would settle the ULIRG question: if the observed flux is inconsistent with the scaled Mrk231 template and the integrated $L_{\\rm IR}$ falls below $10^{12}\\,L_\\odot$, the ULIRG classification fails, while a measured optical-to-mid-infrared lag would test the accretion-rate reprocessing scenario.","supporting_citations":[{"cited_title":"2001, ApJ, 556, 562","cited_arxiv_id":null,"evidence_quote":"Supplies the Mrk231 template that is scaled to the mid-infrared photometry; this template drives the L_IR, SFR_IR, and black-hole mass estimates."},{"cited_title":"J., Stern, D., Noirot, G., et al","cited_arxiv_id":null,"evidence_quote":"Establishes the rarity of mid-infrared variable AGN, providing the comparison population of 162 objects with radio detections."},{"cited_title":"M., Terndrup, D","cited_arxiv_id":null,"evidence_quote":"Supplies the Mrk231 black hole mass used as the scaling reference for the WISE-J0301 mass upper limit."}],"review_version":1}