Pith. sign in

REVIEW 3 major objections 4 minor 2 cited by

On the Variability Features of Active Galactic Nuclei in Little Red Dots

T0 review · 3 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash

Pith's one-line read Little Red Dots show almost no brightness variation over 6–11 observed-frame years, but the paper argues this stillness is consistent with ordinary quasar variability once measurement errors are included.

desk verdict A clean forward-modeling study showing current LRD variability data are noise-limited; the fAGN ≤ 0.3 limit is CHAR-calibrated and the survey requirements (200 sources, 2 yr, 0.07 mag) are the real deliverable. read the letter →

arxiv 2508.16795 v1 pith:MJDXCLBK submitted 2025-08-22 astro-ph.GA

classification astro-ph.GA
keywords LittleRedDotsAGNvariabilityCHARmodelhigh-redshiftgalaxiesaccretiondisksJWSTphotometricmonitoringsupermassiveblackholes
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 asks whether the apparent stillness of Little Red Dots—compact, red, high-redshift sources whose light could come from an active galactic nucleus, a galaxy, or both—is real or simply a measurement limit. Using the Corona-heated Accretion-disk Reprocessing (CHAR) model, a physics-based description of quasar variability, the authors simulate the UV/optical flickering expected from the AGN component in 22 LRDs observed with HST and JWST over rest-frame intervals of about one to three years. They find that the observed magnitude variations are dominated by photometric errors, so the data cannot distinguish a pure-galaxy source from one with a substantial AGN. Within the CHAR model the weak variability is compatible with either an AGN contributing at most about 30% of the observed 1450 Å luminosity (with no extra extinction) or a more luminous AGN hidden behind dust. The paper closes by calculating what observations would settle the question: roughly 200 LRDs, two epochs separated by at least two years, and photometric uncertainty at or below 0.07 mag.

What carries the argument

The load-bearing machinery is the CHAR model, an accretion-disk reprocessing model in which a magnetically coupled corona perturbs the disk heating rate and produces temperature fluctuations on the thermal timescale. Its two governing scalings are a variability amplitude sigma proportional to Lbol,AGN,int^-0.25 and a characteristic timescale tau proportional to alpha^-1 lambda^1.19 Mdot^0.65, where alpha is viscosity, lambda is rest-frame wavelength, and Mdot is accretion rate. The paper maps each LRD's observed luminosity through two parameters—fAGN, the AGN's fraction of observed 1450 Å light, and beta, the ratio of intrinsic to observed AGN luminosity (a dust-extinction or scattering corr

What would settle it

Monitor about 200 spectroscopically confirmed LRDs with two epochs at least two years apart in the observed frame and photometric errors at or below 0.07 mag; if the observed |Δm| distribution is statistically inconsistent with the CHAR-predicted distribution for every fAGN and beta, the null hypothesis fails. A complementary test: measure the rest-frame structure function of a single well-observed LRD and check whether its amplitude and timescale match the CHAR prediction at its estimated AGN luminosity.

Watch

Extended reading notes

Core claim

Under the null hypothesis that LRDs vary like ordinary quasars, the CHAR model predicts that a 22-LRD sample with the luminosities and time baselines of Tee et al. (2025) should show a measurable spread in absolute magnitude changes if the AGN contributes a large fraction of the observed light. Comparing 2500 simulated light curves per parameter choice to the HST/JWST measurements with an Anderson–Darling test, the paper finds that the observed distribution is indistinguishable from models spanning AGN fractions fAGN = 0 to fAGN = 1 once measurement errors are included; only an upper limit fAGN ≲ 0.3 at zero extinction (log beta = 0) survives. The paper's central claim is a dichotomy: the mi

Load-bearing premise

The argument rests on the assumption that Little Red Dots flicker by the same accretion-disk mechanism as ordinary low-redshift quasars; if they instead accrete in a different mode, such as a puffed-up super-Eddington disk, the derived limits on how much light the black hole contributes do not apply.

Editorial extensions

If this is right

  • With the current 22-LRD data, variability alone cannot reject either a pure-galaxy LRD or one powered mostly by an AGN; the only zero-extinction constraint is an AGN fraction upper limit near 30%.
  • A sample of about 200 LRDs each observed twice at least two years apart in the observed frame can distinguish galaxy-dominated from AGN-dominated models, provided photometric errors are tightened to ≤0.07 mag (roughly 0.03 mag at JWST/F444W).
  • If such a sample still shows no variability beyond noise, and galaxy light is not the cause, the time domain would supply independent evidence that LRD accretion differs from conventional low-redshift quasar accretion—for example, super-Eddington disks.
  • The existing ~300-LRD sample already has the needed size and ~0.07 mag precision, but its ~100-day observed-frame separations are too short to constrain fAGN and beta; extending its baselines to two years would make it decisive.
  • The fAGN–beta degeneracy means some models, such as fAGN = 1 with log beta = 2 and fAGN = 0.3 with log beta = 0, produce nearly identical variability, so photometric monitoring alone may never uniquely fix both parameters.

Reading between the lines

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

  • If future two-epoch monitoring of ~200 LRDs detects variability consistent with CHAR but at amplitudes above the fAGN ≤ 0.3 limit, the AGN-fraction interpretation would need revision upward, and the host-dilution picture would weaken.
  • Super-Eddington puffed-up disks have much shorter viscous timescales than standard thin disks, so a dedicated search for fast rest-frame variability in bright LRDs could directly test that scenario—a test the paper discusses only qualitatively.
  • Combining the proposed variability campaign with spectral classification to remove the ~20–30% of photometrically selected LRDs that are galaxies, brown dwarfs, or line-driven fake V-shapes would sharpen the constraints beyond what the paper's pure sample achieves.
  • Because CHAR predicts wavelength-dependent inter-band time delays that scale differently with extinction than with AGN fraction, adding time-delay or polarization measurements could break the fAGN–beta degeneracy that photometric variability alone cannot resolve.
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 / 4 minor

Summary. The paper applies the Corona-heated Accretion-disk Reprocessing (CHAR) model, a physically motivated model of AGN UV/optical variability, to 22 Little Red Dots (LRDs) from Tee et al. (2025) with HST/JWST observations spanning observed-frame intervals of 6–11 years. The authors simulate magnitude variations under different assumed AGN fractions fAGN and dust-correction factors β, compare the simulated and observed |Δm| distributions with Anderson–Darling tests, and find that the observed variability is dominated by measurement uncertainties. Within the CHAR model, the weak observed variability can be explained either by a small AGN contribution (fAGN ≲ 0.3 for log β = 0) or by intrinsically luminous AGNs with high dust extinction. The paper also uses bootstrapped 200-object mock samples to recommend observational requirements—about 200 LRDs, ≥2 yr observed-frame separations, and photometric uncertainty ≤0.07 mag—to distinguish these scenarios in the future. The analysis is explicitly framed under the hypothesis that LRDs share the same variability mechanism as normal quasars.

Significance. The paper is a thoughtful and methodologically careful contribution to the LRD debate. Its strengths are the use of a physically motivated variability model rather than purely empirical scaling relations, the explicit statement of the testable null hypothesis, the forward-modeling approach that avoids fitting the conclusion, and the concrete, testable recommendations for future observing programs. The CHAR model has external validation against SDSS S82 and Kepler data, which reduces circularity concerns. If the CHAR-model calibration holds at LRD luminosities, the derived constraints and the proposed observational requirements are valuable for planning JWST time-domain programs. However, the quantitative results are conditional on an extrapolation of the model to lower luminosities and higher redshift than its validation sample, and on a statistical threshold that is not fully characterized.

major comments (3)
  1. [Sections 2.2 and 3, Eq. (1), Fig. 3] The headline constraint fAGN ≤ 0.3 for log β = 0 is model-dependent in a way that is acknowledged but not quantified. The CHAR model was validated on SDSS S82 quasars with Lbol ~10^45.5–10^47 erg/s (Fig. 1), while the LRDs here have Lbol one to two orders of magnitude lower; the scalings σ ∝ Lbol^-0.25 and τ ∝ α^-1 λ^1.19 Mdot^0.65 are extrapolated without low-luminosity validation, and α is fixed to 0.4. A sensitivity test varying α over a plausible range (e.g., 0.1–1), or a comparison using an empirical scaling relation (e.g., MacLeod et al. 2010), would show how much the fAGN upper limit moves. Without such a test, the numerical constraints in Figures 3–4 should be labeled as illustrative rather than robust.
  2. [Section 2.3 and Figure 3] The statistical decision rule uses an arbitrary threshold P=0.5, where P is the percentage of 2500 simulations with A-D p-value > 0.01. The reported 'upper limit' fAGN=0.3 is the boundary where P crosses this threshold, not a standard confidence limit. The paper does not report how the boundary changes for other threshold choices (e.g., P=0.1 or 0.9), nor is the power of the binned A-D test on |Δm| histograms with N=22 objects characterized beyond the P metric. Because the central conclusion—that observed variability is dominated by measurement uncertainties—rests on the inability to distinguish models, a calibration of the test under the null hypothesis and a threshold sensitivity analysis would materially strengthen the claim.
  3. [Abstract and Section 3, Figure 5] The abstract's 'two scenarios' framing is an oversimplification. Figure 5 shows that fAGN=1/logβ=2, fAGN=0.3/logβ=0, and fAGN=0.5/logβ=1 produce similar structure functions; the allowed region in Figure 3 is a continuous degenerate band, not two discrete alternatives. The paper should state explicitly that variability alone cannot currently separate a low-fAGN/low-β solution from a high-fAGN/high-β solution, and should avoid implying that the two scenarios are mutually exclusive. This is important because the degeneracy is a central message of the paper.
minor comments (4)
  1. [Section 2.3] The choice of P=0.5 as the threshold is introduced without justification. Please provide a sentence explaining why this value is appropriate or report the sensitivity to it.
  2. [Section 3, observational requirements] The sentence about photometric uncertainty '0.1/√2 = 0.07 mag' is unclear. If each epoch has uncertainty σ, the magnitude-difference uncertainty is √2 σ; please spell this out explicitly to avoid confusion.
  3. [Figure 3] The colorbar label 'Percentage [p-value>0.01]' could be confused with a p-value itself. Consider labeling it 'P (%)' with a clear definition in the caption.
  4. [Title and general text] The title displays 'V ariability' with a spurious space (likely a LaTeX issue). There are also minor grammar issues such as 'the model and statistical approach' in Section 3; a light copyedit is recommended.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the CHAR model is self-cited but externally validated; the fAGN constraints are forward-modeled, and the fAGN–β degeneracy is explicitly acknowledged.

full rationale

The paper's central quantitative results—the fAGN ≤ 0.3 upper limit for log β = 0 and the degeneracy with an intrinsically luminous AGN—are obtained by forward simulation, not by fitting the model to LRDs. For each grid point in (fAGN, log β), the authors simulate CHAR light curves, inject the Tee et al. (2025) photometric uncertainties, and compare the resulting |Δm| distributions with observations using Anderson–Darling tests. The inputs are the CHAR model's physical scalings (e.g., σ ∝ Lbol^-0.25, τ ∝ α^-1 λ^1.19 Mdot^0.65), which originate in self-cited work (Sun et al. 2020a; Zhou et al. 2024), but the manuscript also cites external validations (SDSS S82, Kepler, AGN STORM II, Tang et al. 2023, Arévalo et al. 2024). No parameter is adjusted to the LRD data to force the conclusion. The fAGN–log β degeneracy is a genuine consequence of Eq. (1) and the CHAR scaling, and the paper explicitly displays it in Figure 5 and discusses it in Section 3. The acknowledged extrapolation from luminous SDSS S82 quasars to lower-luminosity, z>4 LRDs is a calibration/validity risk, not a circularity; the authors repeatedly flag this assumption ('Our analysis is based on the testable hypothesis…', 'We emphasize that the aforementioned observational requirements are based on the CHAR model…'). Therefore no circular step reduces the result to its inputs by construction; score 1 reflects only the presence of central self-citation with adequate external support.

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

The central constraints depend on the CHAR model (external, self-cited) and on the scanned parameters fAGN and beta, plus fixed alpha. No new physical entities are introduced.

free parameters (4)
  • fAGN = 0 to 1 (grid step 0.1)
    Fraction of observed 1450 Angstrom luminosity contributed by the AGN; scanned, not fitted, and is the primary parameter constrained.
  • log beta = 0 to 2 (grid step 0.2)
    Logarithmic dust correction factor relating observed to intrinsic AGN UV luminosity; scanned.
  • alpha (viscosity parameter) = 0.4
    Fixed to 0.4 following King et al. (2007) and Sun (2023); affects CHAR model timescales and amplitudes.
  • A-D test threshold P = 0.5
    Critical percentage of simulations with p>0.01 used to define distinguishable; chosen by hand and arbitrary.
assumptions (5)
  • domain assumption The CHAR model accurately describes AGN UV/optical variability from low to high luminosity, including the sigma ~ Lbol^-0.25 amplitude scaling.
    The paper uses CHAR predictions as the null hypothesis; if the model fails for LRDs, the inferred fAGN constraints are invalid. The model is taken from Sun et al. (2020a).
  • domain assumption The host galaxy light in LRDs is non-variable during the observed intervals.
    Any variability from stellar processes, supernovae, or other galaxy variability is ignored; the total light is modeled as AGN plus constant galaxy.
  • ad hoc to paper The observed L1450 is the sum of AGN and host, with AGN extincted by dust up to log beta = 2 (AV up to 1 mag) following an SMC extinction curve.
    The wide beta range is motivated by submillimeter/IR observations (Casey et al. 2025; Chen et al. 2025b) but is a free modeling assumption that drives the degeneracy.
  • domain assumption Measurement uncertainties on magnitude variations are Gaussian and accurately estimated by Tee et al. (2025).
    The entire conclusion that variability is noise-dominated relies on the quoted delta merr being correct; underestimated errors would change the result.
  • domain assumption The 22 LRD sample is representative, and bootstrap resampling preserves the relevant distributions.
    Mock samples of 200 are drawn from the 22 with replacement to estimate survey requirements.

how reviews work

0 comments
Cite this review

Pith. "Pith review of On the Variability Features of Active Galactic Nuclei in Little Red Dots." pith.science (2026). https://pith.science/paper/MJDXCLBK

@misc{pith2026250816795,
  author       = {Pith},
  title        = {Pith review of: On the Variability Features of Active Galactic Nuclei in Little Red Dots},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/MJDXCLBK}},
  note         = {Machine review of arXiv:2508.16795}
}
abstract

The high-redshift ($z>4$) compact sources with ``V-shaped" spectral energy distributions (SEDs), known as Little Red Dots (LRDs), are discovered by the James Webb Space Telescope and provide valuable clues to the physics of active galactic nuclei (AGNs) in the early universe. The nature of LRDs is controversial. Recently, several studies have investigated LRDs through variability, a characteristic feature of AGNs. These studies explore LRD variability by extrapolating empirical relationships from local quasars. Here, we adopt the Corona-heated Accretion-disk Reprocessing (CHAR) model, which is motivated by accretion physics and applicable to reproduce AGN conventional variability, to study the variability of $22$ LRDs in \citet{Tee2025}. Our results indicate that the observed variability in LRDs is dominated by measurement uncertainties. Within the CHAR model, the lack of variability in LRDs can be explained by two scenarios: either AGNs contribute $\lesssim30\%$ of the observed luminosities, or they are intrinsically luminous AGNs. We use simulations to demonstrate the observational requirements to effectively investigate LRDs via variability: first, a sample of about $200$ LRDs; second, each LRD has two observations separated by at least two years in the observed frame; third, the photometric uncertainty is $\leq 0.07$ mag. If the LRDs still lack variability under these conditions, the time-domain study would provide independent evidence that the accretion mode of LRDs differs significantly from low-redshift quasars.

Figures

Figures reproduced from arXiv: 2508.16795 by the authors.

Figure 1
Figure 1. Rest-frame wavelengths and bolometric lumi￾nosities for LRDs and SDSS S82 quasars. The red dots and purple squares are the results for the JWST/F115W and JWST/F150W filters, respectively. The bolometric lumi￾nosities for LRDs are estimated under the assumption that the observed L1450,obs are entirely due to AGN activities and without considering dust extinction for AGNs. The blue con￾tours and dots are SDSS S82 quas… view at source ↗
Figure 2
Figure 2. Distributions of the absolute magnitude variation (|∆mF115W| or |∆mF150W|) for models with different AGN contributions (fAGN) and log β = 0, and real observations of Tee et al. (2025). The first and second rows show the results for the JWST/F115W and JWST/F150W filters, respectively. In all panels, the black square, red dots, green diamonds, and blue pentagons represent the AGN contributions (i.e., fAGN) of 0, 0.1, … view at source ↗
Figure 3
Figure 3. The percentage (P) of simulations with p-values from the A-D test between the observed and simulated magnitude variations that are greater than 0.01. The right y-axis shows the V-band extinction (AV) corresponding to log β, which is based on the SMC extinction curve. The left and right panels represent the results for the JWST/F115W and JWST/F150W filters, respectively. A lower percentage (redder in color) suggests … view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: The values of P for models with different AGN contributions and luminosity correction factors compared to a reference model. The right y-axis shows the V-band extinction (AV) corresponding to log β, which is based on the SMC extinction curve. The reference models (blac…
Figure 5
Figure 5. Figure 5: Examples of structure functions for an LRD at a redshift of 7.04 with different AGN contributions and lumi￾nosity correction factors. The vertical dashed line indicates the rest-frame time corresponding to an observed time inter￾val of two years. The horizontal dashed …

Discussion (0). Sign in to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. NEXUS: A Search for Nuclear Variability with the First Two JWST NIRCam Epochs

    astro-ph.GA 2025-09 conditional novelty 6.0 of 10

    Using two JWST NIRCam epochs, difference imaging finds 465 nuclear variable sources and sets tight F444W variability upper limits of 3 to 10 percent for ten Little Red Dots.

  2. Do Little Red Dots Vary?

    astro-ph.GA 2025-09 conditional novelty 6.0 of 10

    Super-Eddington accretion models can explain why little red dots show almost no variability, whereas standard sub-Eddington AGN variability models predict changes that should already have been seen.

Reference graph

Works this paper leans on

2 extracted references · 1 canonical work pages · cited by 2 Pith papers

  1. [1]

    A., Czerny, B., Lasota, J

    Abramowicz, M. A., Czerny, B., Lasota, J. P., & Szuszkiewicz, E. 1988, ApJ, 332, 646, doi: 10.1086/166683 Akins, H. B., Casey, C. M., Lambrides, E., et al. 2024, arXiv e-prints, arXiv:2406.10341, doi: 10.48550/arXiv.2406.10341 Ananna, T. T., Bogd´ an,´A., Kov´ acs, O. E., Natarajan, P., & Hickox, R. C. 2024, ApJL, 969, L18, doi: 10.3847/2041-8213/ad5669 A...

  2. [918]

    I., & Sunyaev, R

    http://www.jstor.org/stable/2288805 Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 Shen, Y., Richards, G. T., Strauss, M. A., et al. 2011, ApJS, 194, 45, doi: 10.1088/0067-0049/194/2/45 Simmonds, C., Tacchella, S., Hainline, K., et al. 2024, MNRAS, 535, 2998, doi: 10.1093/mnras/stae2537 Son, S., Kim, M., & Ho, L. C. 2025, A&A, 695, A268, doi: 10.1051...

Pith tools

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