Connecting the long-term variability behaviour of active galactic nuclei to their central engines
Pith reviewed 2026-05-08 17:03 UTC · model grok-4.3
The pith
The inverse of the PSD bend frequency in AGN radio light curves matches the mean duration of the brightest flares and correlates with black hole mass over accretion rate.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors find that the inverse of the PSD bend frequency of radio light curves best corresponds to the mean duration of the brightest flares from the decomposition. For some sources the mean flare separation has a comparable timescale. Treating flare durations and separations as proxies for the PSD timescale yields a positive correlation with black hole mass divided by the normalised mass accretion rate, which indicates that the variability timescales are associated with the central engine.
What carries the argument
Decomposition of radio light curves into individual flares, which isolates rise times, durations, and separations to compare directly against PSD bend frequencies and central engine parameters such as black hole mass and accretion rate.
Load-bearing premise
The flare decomposition accurately isolates the intrinsic variability components without bias from the fitting choices or data sampling, and the reported correlations are not driven primarily by redshift or selection effects.
What would settle it
A sample of AGNs where independently measured PSD bend frequencies fail to match the inverse of mean brightest flare durations from light curve decomposition, or where the correlation with black hole mass over accretion rate vanishes after controlling for redshift.
Figures
read the original abstract
Analysing the long-term radio variability of active galactic nuclei (AGNs) is essential to understanding the physics of relativistic jets launched by supermassive black holes. We aim to connect the characteristic timescales obtained from a prior power spectral density (PSD) analysis to the decomposed timescales of the light curves. In addition, we probe for potential associations between the timescales and the physical characteristics of the relativistic jet as well as the central engine. We decomposed the long-term radio light curves of 54 sources observed at the Aalto University Mets\"ahovi Radio Observatory into individual flares to understand which timescale of variability is related to the low-frequency bend in the PSD. In addition, we used the obtained rise times of the brightest flares to look for associations between the emission-region size in the jet and different central engine parameters. We found that the inverse of the PSD bend frequency of radio light curves best corresponds to the mean duration of the brightest flares. For some sources, the mean flare separation had a similar timescale. Using the flare durations and separations as proxies for the PSD timescale, we found a positive correlation with black hole mass divided by the normalised mass accretion rate. This suggests that the variability timescales obtained from the PSDs of radio light curves are associated with the central engine. Furthermore, when comparing the obtained rise times of the brightest flares to the jet and central engine parameters, we found weak tentative correlations, but they may be driven by a common dependency on redshift.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript analyzes long-term radio light curves of 54 AGNs observed at Metsähovi, decomposes them into individual flares, and connects the resulting timescales to a prior PSD analysis. It claims that the inverse of the PSD bend frequency corresponds to the mean duration of the brightest flares (with mean flare separation similar in some sources), reports a positive correlation between these timescales and black hole mass divided by normalized mass accretion rate, and finds only weak tentative correlations between brightest-flare rise times and jet/central-engine parameters that may be redshift-driven.
Significance. If the reported correspondence and correlation survive quantitative statistical tests and explicit controls for decomposition choices and redshift, the work would supply an empirical bridge between radio variability timescales and central-engine parameters, which is of interest for jet physics in AGNs. The approach of using flare decomposition as a proxy for PSD bend frequency is conceptually straightforward, but the absence of error bars, robustness checks, and bias quantification leaves the result preliminary rather than definitive.
major comments (3)
- [Abstract and results] Abstract and results sections: the claimed best correspondence between 1/f_bend and mean brightest-flare duration is stated without any statistical measure (correlation coefficient, p-value, or uncertainty on the mean durations), error bars on the flare parameters, or robustness tests against variations in the decomposition algorithm. This makes the central claim difficult to evaluate quantitatively.
- [Correlation analysis] Correlation analysis: the positive correlation between the PSD-derived timescales (via flare durations/separations) and M_BH/ṁ is presented without a control for redshift or luminosity selection effects, even though the abstract itself flags a possible redshift dependence for the rise-time correlations. A partial-correlation test or redshift-binned analysis is required to determine whether the reported trend is physical.
- [Methods (flare decomposition)] Flare decomposition method: no quantitative description is given of the flare-identification thresholds, baseline-subtraction procedure, or amplitude/duration selection criteria. Because the mean duration of the brightest flares is used as the direct proxy for 1/f_bend, any systematic bias introduced by these choices would propagate directly into the claimed correspondence and the subsequent correlation.
minor comments (2)
- [Introduction or methods] The phrase 'normalised mass accretion rate' is used without an explicit definition or reference to the normalization convention adopted; this should be clarified in the text or a dedicated equation.
- [Figures and results] Figure captions and text should state the exact number of sources for which both PSD bend frequencies and flare decompositions are available, and whether any sources were excluded from the correlation analysis.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed report. We agree that the original manuscript would benefit from additional quantitative statistical support, explicit controls for selection effects, and more precise methodological descriptions. We have revised the manuscript to address each point and provide the following point-by-point responses.
read point-by-point responses
-
Referee: [Abstract and results] Abstract and results sections: the claimed best correspondence between 1/f_bend and mean brightest-flare duration is stated without any statistical measure (correlation coefficient, p-value, or uncertainty on the mean durations), error bars on the flare parameters, or robustness tests against variations in the decomposition algorithm. This makes the central claim difficult to evaluate quantitatively.
Authors: We acknowledge that the submitted version lacked explicit statistical quantification and robustness checks for the claimed correspondence. In the revised manuscript we have added the Spearman rank correlation coefficient and p-value between 1/f_bend and the mean duration of the brightest flares, together with uncertainties on the per-source mean durations. Error bars are now shown on the relevant figure. We have also performed robustness tests by varying the flare-identification threshold and baseline window by ±20 % and confirmed that the correspondence remains stable. These additions appear in the revised results section and a new supplementary figure. revision: yes
-
Referee: [Correlation analysis] Correlation analysis: the positive correlation between the PSD-derived timescales (via flare durations/separations) and M_BH/ṁ is presented without a control for redshift or luminosity selection effects, even though the abstract itself flags a possible redshift dependence for the rise-time correlations. A partial-correlation test or redshift-binned analysis is required to determine whether the reported trend is physical.
Authors: We agree that redshift and luminosity selection effects require explicit control. The revised manuscript now includes a partial-correlation analysis that controls for redshift and shows the correlation between the flare-duration proxy and M_BH/ṁ remains significant. We have additionally performed a redshift-binned analysis (low-z and high-z subsamples) in which the positive trend is recovered in both bins. These controls are described in the updated correlation subsection. revision: yes
-
Referee: [Methods (flare decomposition)] Flare decomposition method: no quantitative description is given of the flare-identification thresholds, baseline-subtraction procedure, or amplitude/duration selection criteria. Because the mean duration of the brightest flares is used as the direct proxy for 1/f_bend, any systematic bias introduced by these choices would propagate directly into the claimed correspondence and the subsequent correlation.
Authors: We have expanded the methods section with quantitative specifications: the baseline is removed via a 3-year running median, flares are identified when the flux exceeds 2.5 times the local rms for at least three epochs, and the brightest flares are defined as those with peak flux >0.3 Jy and fitted duration between 0.2 and 10 yr. The mean duration is computed from the top 20 % brightest flares per source. These details are now provided to allow reproducibility and to quantify possible systematic effects. revision: yes
Circularity Check
No significant circularity in empirical timescale matching
full rationale
The paper reports an observational correspondence between the inverse PSD bend frequency (from a prior analysis) and the mean duration of brightest flares after decomposing 54 radio light curves, plus a secondary correlation with M_BH / normalized accretion rate. No mathematical derivation, first-principles prediction, or fitted parameter is presented that reduces by construction to the inputs; the flare durations are treated as independent proxies whose match to 1/f_bend is measured directly. The prior PSD work is referenced but does not bear the load of the new connection claim, which remains falsifiable against the decomposition results. Redshift caveats are explicitly flagged, confirming the analysis is self-contained against external benchmarks rather than self-referential.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The PSD bend frequency represents a physically meaningful characteristic variability timescale comparable to flare durations.
Reference graph
Works this paper leans on
-
[1]
Dissecting magnetar variability with Bayesian hierarchical models , author=. ApJ , volume=. 2015 , publisher=
work page 2015
-
[2]
The relation between accretion rate and jet power in X-ray luminous elliptical galaxies , author=. MNRAS , volume=. 2006 , publisher=
work page 2006
-
[3]
Planck 2018 results. VI. Cosmological parameters , author=. A&A , volume=
work page 2018
-
[4]
The structure of the M87 jet: a transition from parabolic to conical streamlines , author=. ApJL , volume=. 2012 , publisher=
work page 2012
-
[5]
Evolution of the 0.01--25 Hz power spectral components in Cygnus X-1 , author=. A&A , volume=. 2005 , publisher=
work page 2005
-
[6]
Limits of X-ray variability in active galactic nuclei , author=. Nature , volume=. 1986 , publisher=
work page 1986
- [7]
-
[8]
Astrophysical Journal, Part 1, vol
Relativistic jets as compact radio sources , author=. Astrophysical Journal, Part 1, vol. 232, Aug. 15, 1979, p. 34-48. , volume=
work page 1979
-
[9]
DNest4: Diffusive nested sampling in C++ and Python , author=. J. Stat. Softw , volume=
-
[10]
Diffusive nested sampling , author=. Stat. Comput. , volume=. 2011 , publisher=
work page 2011
-
[11]
Jet collimation in NGC 315 and other nearby AGN , author=. A&A , volume=. 2021 , publisher=
work page 2021
-
[12]
The effective acceleration of plasma outflow in the paraboloidal magnetic field , author=. MNRAS , volume=. 2006 , publisher=
work page 2006
-
[13]
Relativistic jets from active galactic nuclei , author=. 2012 , publisher=
work page 2012
- [14]
-
[15]
DISK--JET CONNECTION IN THE RADIO GALAXY 3C 120 , author=. ApJ , volume=. 2009 , publisher=
work page 2009
-
[16]
Connection between the accretion disk and jet in the radio galaxy 3C 111 , author=. ApJ , volume=. 2011 , publisher=
work page 2011
- [17]
-
[18]
A Cutoff in the x-ray fluctuation power density spectrum of the Seyfert 1 galaxy NGC 3516 , author=. ApJ , volume=. 1999 , publisher=
work page 1999
-
[19]
A scheme to unify low-power accreting black holes-Jet-dominated accretion flows and the radio/X-ray correlation , author=. A&A , volume=. 2004 , publisher=
work page 2004
-
[20]
Jet Collimation Profile of the Low-luminosity Active Galactic Nucleus M84: Insight into the Jet Formation in the Low-accretion Regime , author=. ApJ , volume=. 2025 , publisher=
work page 2025
-
[21]
Towards a unified model for black hole X-ray binary jets , author=. MNRAS , volume=. 2004 , publisher=
work page 2004
-
[22]
Jets from black hole x--ray binaries: Testing, refining and extending empirical models for the coupling to x--rays , author=. MNRAS , volume=. 2009 , publisher=
work page 2009
-
[23]
A high signal-to-noise ratio composite quasar spectrum , author=. 1991 , publisher=
work page 1991
-
[24]
On the location of the supermassive black hole in CTA 102 , author=. A&A , volume=. 2015 , publisher=
work page 2015
-
[25]
A universal radio--X-ray correlation in low/hard state black hole binaries , author=. MNRAS , volume=. 2003 , publisher=
work page 2003
-
[26]
A wide and collimated radio jet in 3C84 on the scale of a few hundred gravitational radii , author=. Nature Astronomy , volume=. 2018 , publisher=
work page 2018
-
[27]
On the nature of rapid X-ray variability in active galactic nuclei , author=. MNRAS , volume=. 1993 , publisher=
work page 1993
-
[28]
X-ray variability of 104 active galactic nuclei-XMM-Newton power-spectrum density profiles , author=. A&A , volume=. 2012 , publisher=
work page 2012
-
[29]
A general mechanical model for| f| spectral density random noise with special reference to flicker noise 1/| f| , author=. Proc. IEEE , volume=. 1968 , publisher=
work page 1968
-
[30]
International Journal of Modern Physics: Conference Series , volume=
Physical properties of jets in AGN , author=. International Journal of Modern Physics: Conference Series , volume=. 2012 , organization=
work page 2012
-
[31]
The University of Michigan radio astronomy data base. I-Structure function analysis and the relation between BL Lacertae objects and quasi-stellar objects , author=. ApJ , volume=
-
[32]
Relation between X-ray emission and optical spectra , author=
A combined optical and X-ray study of unobscured type 1 active galactic nuclei--II. Relation between X-ray emission and optical spectra , author=. MNRAS , volume=. 2012 , publisher=
work page 2012
-
[33]
Long-term radio variability of active galactic nuclei at 37 GHz , author=. A&A , volume=. 2025 , publisher=
work page 2025
-
[34]
Active galactic nucleus time-variability analysis and its caveats , author=. A&A , volume=. 2025 , publisher=
work page 2025
-
[35]
Are the variations in quasar optical flux driven by thermal fluctuations? , author=. ApJ , volume=. 2009 , publisher=
work page 2009
-
[36]
ppcor: an R package for a fast calculation to semi-partial correlation coefficients , author=. Commun. Stat. Appl. Methods , volume=
-
[37]
Galaxy group-associated distances to very high energy gamma-ray emitting BL Lacs KUV 00311- 1938 and S2 0109+ 22 , author=. MNRAS , volume=. 2024 , publisher=
work page 1938
-
[38]
A transition from parabolic to conical shape as a common effect in nearby AGN jets , author=. MNRAS , volume=. 2020 , publisher=
work page 2020
-
[39]
The core shift effect in the blazar 3C 454.3 , author=. MNRAS , volume=. 2014 , publisher=
work page 2014
-
[40]
X-ray variability of active galactic nuclei-A universal power spectrum with luminosity-dependent amplitude , author=. ApJ , volume=
-
[41]
Constraining the limiting brightness temperature and Doppler factors for the largest sample of radio-bright blazars , author=. ApJ , volume=. 2018 , publisher=
work page 2018
-
[42]
Synchrotron flaring behaviour of Cygnus X-3 during the February-March 1994 and September 2001 outbursts , author=. A&A , volume=. 2007 , publisher=
work page 1994
-
[43]
MOJAVE: monitoring of jets in active galactic nuclei with VLBA experiments. VI. Kinematics analysis of a complete sample of blazar jets , author=. The Astronomical Journal , volume=. 2009 , publisher=
work page 2009
-
[44]
Relativistic beaming and flux variability in active galactic nuclei , author=. ApJ , volume=. 2001 , publisher=
work page 2001
-
[45]
MOJAVE. XV. VLBA 15 GHz total intensity and polarization maps of 437 parsec-scale AGN jets from 1996 to 2017 , author=. ApJS , volume=. 2018 , publisher=
work page 1996
-
[46]
MOJAVE. XVII. Jet kinematics and parent population properties of relativistically beamed radio-loud blazars , author=. ApJ , volume=. 2019 , publisher=
work page 2019
-
[47]
Total flux density variations in extragalactic radio sources. III. Doppler boosting factors, lorentz factors, and viewing angles for active galactic nuclei , author=. ApJ , volume=. 1999 , publisher=
work page 1999
-
[48]
Asymptotic structure of Poynting-dominated jets , author=. ApJ , volume=. 2009 , publisher=
work page 2009
-
[49]
A jet model for the fast IR variability of the black hole X-ray binary GX 339-4 , author=. MNRAS , volume=. 2018 , publisher=
work page 2018
-
[50]
Observational evidence for the accretion-disk origin for a radio jet in an active galaxy , author=. Nature , volume=. 2002 , publisher=
work page 2002
-
[51]
Active galactic nuclei as scaled-up Galactic black holes , author=. Nature , volume=. 2006 , publisher=
work page 2006
-
[52]
The Jet Paradigm: From Microquasars to Quasars , pages=
X-Ray variability of AGN and relationship to Galactic black hole binary systems , author=. The Jet Paradigm: From Microquasars to Quasars , pages=. 2009 , publisher=
work page 2009
-
[53]
Time correlation between the radio and gamma-ray activity in blazars and the production site of the gamma-ray emission , author=. MNRAS , volume=. 2014 , publisher=
work page 2014
-
[54]
A jet model for the broadband spectrum of XTE J1118+ 480-Synchrotron emission from radio to X-rays in the Low/Hard spectral state , author=. A&A , volume=. 2001 , publisher=
work page 2001
-
[55]
As above, so below: exploiting mass scaling in black hole accretion to break degeneracies in spectral interpretation , author=. ApJL , volume=. 2015 , publisher=
work page 2015
-
[56]
A Fundamental Plane of black hole activity , author=. MNRAS , volume=. 2003 , publisher=
work page 2003
-
[57]
The accretion disc-jet connection in blazars , author=. MNRAS , volume=. 2019 , publisher=
work page 2019
-
[58]
Physical parameters of active galactic nuclei derived from properties of the jet geometry transition region , author=. MNRAS , volume=. 2020 , publisher=
work page 2020
-
[59]
Periodogram analysis and continuous spectra , author=. Biometrika , volume=. 1950 , publisher=
work page 1950
-
[60]
On the fundamental dichotomy in the local radio-AGN population: accretion, evolution and host galaxy properties , author=. MNRAS , volume=. 2012 , publisher=
work page 2012
-
[61]
The transition between BL Lac objects and flat spectrum radio quasars , author=. MNRAS , volume=. 2011 , publisher=
work page 2011
-
[62]
The power of relativistic jets is larger than the luminosity of their accretion disks , author=. Nature , volume=. 2014 , publisher=
work page 2014
-
[63]
Collimation, acceleration, and recollimation shock in the jet of gamma-ray emitting radio-loud narrow-line Seyfert 1 Galaxy 1H0323+ 342 , author=. ApJ , volume=. 2018 , publisher=
work page 2018
-
[64]
Spectrum and spectral density estimation by the Discrete Fourier transform (DFT), including a comprehensive list of window functions and some new at-top windows , author=
-
[65]
Doppler factors, Lorentz factors and viewing angles for quasars, BL Lacertae objects and radio galaxies , author=. A&A , volume=. 2009 , publisher=
work page 2009
-
[66]
The jet--disc connection in AGN , author=. MNRAS , volume=. 2014 , publisher=
work page 2014
-
[67]
On characterizing the variability properties of X-ray light curves from active galaxies , author=. MNRAS , volume=. 2003 , publisher=
work page 2003
-
[68]
Astronomical Journal , volume=
VLA observations of objects in the Palomar Bright Quasar Survey , author=. Astronomical Journal , volume=
-
[69]
Astrophysics and Space Science , volume=
Doppler boosting, superluminal motion, and the kinematics of AGN jets , author=. Astrophysics and Space Science , volume=. 2007 , publisher=
work page 2007
-
[70]
Five-year wilkinson microwave anisotropy probe* observations: cosmological interpretation , author=. ApJS , volume=. 2009 , publisher=
work page 2009
-
[71]
Opacity, variability, and kinematics of AGN jets , author=. MNRAS , volume=. 2019 , publisher=
work page 2019
-
[72]
Synchrotron flaring in the jet of 3C 279 , author=. A&A , volume=. 2006 , publisher=
work page 2006
-
[73]
Reverberation Measurements for 17 Quasars and theSize-Mass-Luminosity Relations in Active GalacticNuclei , author=. ApJ , volume=. 2000 , publisher=
work page 2000
-
[74]
Optical and radio variability of the northern VHE gamma-ray emitting BL Lacertae objects , author=. A&A , volume=. 2016 , publisher=
work page 2016
-
[75]
The jet power, radio loudness, and black hole mass in radio-loud active galactic nuclei , author=. ApJ , volume=. 2006 , publisher=
work page 2006
-
[76]
Identifying changing jets through their radio variability , author=. A&A , volume=. 2021 , publisher=
work page 2021
-
[77]
Models for high-frequency radio outbursts in extragalactic sources, with application to the early 1983 millimeter-to-infrared flare of 3C 273 , author=. ApJ , volume=
work page 1983
-
[78]
The jet-disk connection and blazar unification , author=. ApJ , volume=. 2003 , publisher=
work page 2003
-
[79]
arXiv preprint physics/0204033 , year=
1/f noise: a pedagogical review , author=. arXiv preprint physics/0204033 , year=
-
[80]
Collimation of the Relativistic Jet in the Quasar 3C 273 , author=. ApJ , volume=. 2022 , publisher=
work page 2022
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.