REVIEW 3 major objections 2 cited by
Optical Variability Structure Function of Low-Luminosity AGN using ATLAS Lightcurves
T0 review · 3 major / 0 minor · reviewed 2026-07-01 · grok-4.3
Pith's one-line read The slope of the optical variability structure function in low-luminosity AGN rises with black hole mass from ~0.1 to ~0.3.
desk verdict The mass-dependent SF slope in low-luminosity AGN is the main new claim, but it rests on host subtraction that the abstract leaves untested. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The ensemble variability structure function, which tracks the average squared flux difference as a function of time separation after host subtraction.
What would settle it
An independent sample of low-luminosity AGN with accurate host fractions showing either a flat slope-mass relation or clear breaks in the structure function at long lags would falsify the central claims.
Extended reading notes
Core claim
In a sample of 246 low-luminosity AGN at z < 0.1, the ensemble structure function slope increases with black hole mass from approximately 0.1 at log M_BH/M_⊙ ~ 6.5 to 0.3 at log M_BH/M_⊙ ~ 8. No breaks appear in the structure function, and spectra taken 20 years apart indicate that the function keeps rising into decadal timescales. The amplitude shows an anticorrelation with luminosity and a positive correlation with black hole mass. Variability patterns suggest extinction does not drive the range of Seyfert subtypes.
Load-bearing premise
Host-AGN spectral decomposition on recent spectra yields an accurate host fraction that isolates intrinsic AGN variability in the light curves without residual contamination or bias.
Editorial extensions
If this is right
- Variability grows continuously with time lag without a preferred break timescale.
- More massive black holes produce steeper growth of variability amplitude with lag.
- Amplitude decreases at higher luminosities while increasing with black hole mass.
- Seyfert subtype differences arise from factors other than dust extinction.
Reading between the lines
- Models of AGN variability must treat black hole mass as a controlling parameter rather than a secondary effect.
- Multi-decade monitoring campaigns could directly test whether the structure function continues its rise or eventually saturates.
- The mass dependence offers a new observable to distinguish between competing physical drivers such as disk instabilities versus reprocessing.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript analyzes the ensemble optical variability structure function (SF) for a sample of 246 low-luminosity AGN at z < 0.1 drawn from the 6dFGS, using ATLAS lightcurves with ~2-day cadence over eight years. After host-AGN decomposition on recent spectra to isolate the AGN flux, the authors report that the SF power-law slope increases with black hole mass (from ~0.1 at log M_BH/M_⊙ ~6.5 to ~0.3 at log M_BH/M_⊙ ~8), find no evidence for breaks in the SF, infer continued rise on decadal timescales from two-epoch spectra separated by 20 years, measure an anticorrelation of SF amplitude with luminosity and positive correlation with black hole mass, and argue that extinction is not the primary driver of Seyfert subtype differences.
Significance. If the host subtraction and SF measurements are robust, the reported mass dependence of the SF slope in the low-luminosity regime would extend prior high-luminosity studies and provide a new observational constraint on variability mechanisms. The absence of breaks and evidence for continued rise on long timescales would challenge some earlier results. The large ensemble from ATLAS and 6dFGS enables statistical power not available in smaller samples.
major comments (3)
- [Abstract] Abstract: The exact mathematical definition of the structure function is not provided (e.g., whether it follows the standard form SF(τ) = sqrt(<[m(t+τ)−m(t)]²>) or a normalized variant), nor are uncertainties or error bars reported on the binned SF values or the fitted slopes in mass bins; without these, the statistical significance of the claimed increase from ~0.1 to ~0.3 cannot be evaluated.
- [Abstract] Abstract: Sample selection criteria are not specified, including the precise definition of 'low-luminosity', any cuts on spectral quality, variability amplitude, or host dominance, and how the final 246 objects were chosen from the parent 6dFGS catalog; such details are required to assess potential selection biases that could correlate with black hole mass and affect the reported SF-mass trend.
- [Abstract] Abstract: The host-AGN decomposition procedure used to derive the host fraction for subtraction from the ATLAS lightcurves is described only at the level of 'perform host-AGN decomposition on recent spectra'; no information is given on the spectral templates, fitting method, aperture corrections, or any validation (e.g., via simulations or comparison to independent decompositions), which is load-bearing for the central SF slope vs. mass result because mass-correlated residuals would directly bias the normalized variability amplitudes.
Simulated Author's Rebuttal
We thank the referee for their careful and constructive review of our manuscript. Their comments highlight areas where additional clarity is needed, particularly in the abstract and supporting details. We address each major comment below and will revise the manuscript to incorporate the requested information, improving transparency without altering the core results.
read point-by-point responses
-
Referee: [Abstract] Abstract: The exact mathematical definition of the structure function is not provided (e.g., whether it follows the standard form SF(τ) = sqrt(<[m(t+τ)−m(t)]²>) or a normalized variant), nor are uncertainties or error bars reported on the binned SF values or the fitted slopes in mass bins; without these, the statistical significance of the claimed increase from ~0.1 to ~0.3 cannot be evaluated.
Authors: We agree that the abstract should explicitly state the structure function definition and include uncertainties to allow readers to assess the significance of the slope trend. The analysis uses the standard definition SF(τ) = sqrt(<[m(t+τ)−m(t)]²>), consistent with prior AGN variability literature. In the revised manuscript, we will add this definition to the abstract, report error bars on binned SF values (derived from bootstrap resampling of the ensemble), and include uncertainties on the fitted slopes per mass bin. The reported increase from ~0.1 to ~0.3 remains statistically significant given the sample size and binning, but the added details will make this evaluation straightforward. revision: yes
-
Referee: [Abstract] Abstract: Sample selection criteria are not specified, including the precise definition of 'low-luminosity', any cuts on spectral quality, variability amplitude, or host dominance, and how the final 246 objects were chosen from the parent 6dFGS catalog; such details are required to assess potential selection biases that could correlate with black hole mass and affect the reported SF-mass trend.
Authors: The abstract is necessarily concise, but we acknowledge that key selection details should be summarized for context. The full manuscript details the parent 6dFGS sample, with low-luminosity defined via bolometric luminosity cuts (log L_bol < 44 erg s^{-1}) and black hole mass estimates from spectral fitting, plus quality cuts on spectra (S/N > 10) and reliable variability measurements. No explicit cuts on variability amplitude or host dominance were applied beyond the decomposition step. In revision, we will add a concise summary of these criteria to the abstract or a new methods subsection, including the step-by-step reduction from the parent catalog to the final 246 objects, and explicitly discuss checks for mass-correlated selection biases. revision: yes
-
Referee: [Abstract] Abstract: The host-AGN decomposition procedure used to derive the host fraction for subtraction from the ATLAS lightcurves is described only at the level of 'perform host-AGN decomposition on recent spectra'; no information is given on the spectral templates, fitting method, aperture corrections, or any validation (e.g., via simulations or comparison to independent decompositions), which is load-bearing for the central SF slope vs. mass result because mass-correlated residuals would directly bias the normalized variability amplitudes.
Authors: We recognize that the host-AGN decomposition is central to the results and requires fuller documentation. The procedure uses recent 6dFGS spectra with stellar population templates from Bruzual & Charlot (2003), fitted via pPXF to separate AGN and host components, with aperture corrections based on the fiber size and seeing. Validation includes comparisons to independent decompositions for a subset and mock spectrum tests to quantify residuals. In the revised manuscript, we will expand the description in the methods section (and reference it in the abstract) to include these specifics, along with quantitative validation results, to address concerns about potential mass-dependent biases. revision: yes
Circularity Check
No circularity: purely observational measurement of SF from external lightcurves after host subtraction
full rationale
The paper performs ensemble structure-function analysis on ATLAS lightcurves of 246 low-luminosity AGN after subtracting a host fraction obtained from independent recent spectra. The reported trends (SF slope vs. black-hole mass, amplitude vs. luminosity, absence of breaks) are direct empirical correlations extracted from these data products. No step fits a parameter to a subset and then presents a closely related quantity as a prediction; no result is defined in terms of itself; no uniqueness theorem or ansatz is imported via self-citation; and no known empirical pattern is merely renamed. The derivation chain therefore remains self-contained against external benchmarks and does not reduce to its inputs by construction.
Assumptions & free parameters
free parameters (1)
- SF slope in mass bins =
0.1 to 0.3
assumptions (1)
- domain assumption Host-galaxy light subtraction via spectral decomposition isolates the variable AGN component without bias.
Cite this review
Pith. "Pith review of Optical Variability Structure Function of Low-Luminosity AGN using ATLAS Lightcurves." pith.science (2026). https://pith.science/paper/7CLZHVWR
@misc{pith2026260503577,
author = {Pith},
title = {Pith review of: Optical Variability Structure Function of Low-Luminosity AGN using ATLAS Lightcurves},
year = {2026},
howpublished = {\url{https://pith.science/paper/7CLZHVWR}},
note = {Machine review of arXiv:2605.03577}
}
abstract
The origin of the optical flux variability in active galactic nuclei (AGN) is largely unknown. Previous studies have correlated features of the variability structure function (SF) with AGN properties, though they mostly involved high-luminosity AGN to avoid biases from host galaxy flux. In this work, we characterise optical variability in a sample of 246 low-luminosity AGN at $z < 0.1$ from the Six-degree Field Galaxy Survey (6dFGS) through the ensemble variability SF. We use lightcurves from the Asteroid Terrestrial-impact Last Alert System (ATLAS) with a cadence of $\sim$2 days over eight years, and perform host-AGN decomposition on recent spectra to obtain the host fraction. We find that the slope of the SF depends on black hole mass, increasing from $\sim 0.1$ at $\log M_{\mathrm{BH}}/M_\odot \sim 6.5$ to $\sim 0.3$ at $\log M_{\mathrm{BH}}/M_\odot \sim 8$. Contrary to some earlier work, we do not find breaks in the SF, and two-epoch spectra taken 20 years apart suggest that the SF keeps rising into decadal timescales. In addition, we measure an anticorrelation of the amplitude with the luminosity and a positive correlation with the black hole mass. The variability behaviour also suggests that extinction is not the main driver of the variety in Seyfert subtypes.
Figures
Figures from the paper (7 more)
Forward citations
Cited by 2 Pith papers
-
Automatically distinguishing Rubin transients from AGN using variability metrics
A two-dimensional cut on detection flux divided by pre-detection standard deviation and mean flux separates non-AGN transients from AGN with F1 scores of 0.76 to 0.81 on ZTF and simulated LSST data.
-
OzSSy1: The Australian Southern Seyfert-1 Spectroscopic Atlas and Catalogue at z < 0.1
Presents a public spectroscopic atlas of 887 Southern broad-line AGNs at low redshift, including continuum decompositions and emission-line fits.
Reference graph
Works this paper leans on
-
[1]
A., et al., 2009, Lsst science book, version 2.0 Amrutha N., Wolf C., Onken C
Abell P. A., et al., 2009, Lsst science book, version 2.0 Amrutha N., Wolf C., Onken C. A., Hon W. J., Lai S., Tonry J. L., Webster R., 2024, MNRAS, 535, 2322 AmruthaN.,WolfC.,OnkenC.,HonW.,LaiS.,RaithelD.,TanA.,Webster R., 2026, Nature Communications Arévalo P., Uttley P., Lira P., Breedt E., McHardy I. M., Churazov E., 2009, MNRAS, 397, 2004 Arévalo P.,...
work page 2009
Reviewed July 1, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.