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]
ApJ , volume=
Dissecting magnetar variability with Bayesian hierarchical models , author=. ApJ , volume=. 2015 , publisher=
2015
-
[2]
MNRAS , volume=
The relation between accretion rate and jet power in X-ray luminous elliptical galaxies , author=. MNRAS , volume=. 2006 , publisher=
2006
-
[3]
Planck 2018 results. VI. Cosmological parameters , author=. A&A , volume=
2018
-
[4]
ApJL , volume=
The structure of the M87 jet: a transition from parabolic to conical streamlines , author=. ApJL , volume=. 2012 , publisher=
2012
-
[5]
A&A , volume=
Evolution of the 0.01--25 Hz power spectral components in Cygnus X-1 , author=. A&A , volume=. 2005 , publisher=
2005
-
[6]
Nature , volume=
Limits of X-ray variability in active galactic nuclei , author=. Nature , volume=. 1986 , publisher=
1986
-
[7]
A&A , volume=
Variability in the noise properties of Cygnus X-1 , author=. A&A , volume=
-
[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=
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=
2011
-
[11]
A&A , volume=
Jet collimation in NGC 315 and other nearby AGN , author=. A&A , volume=. 2021 , publisher=
2021
-
[12]
MNRAS , volume=
The effective acceleration of plasma outflow in the paraboloidal magnetic field , author=. MNRAS , volume=. 2006 , publisher=
2006
-
[13]
2012 , publisher=
Relativistic jets from active galactic nuclei , author=. 2012 , publisher=
2012
-
[14]
ApJ , volume=
The nature of the BL Lacertae object AO 0235+ 164 , author=. ApJ , volume=
-
[15]
ApJ , volume=
DISK--JET CONNECTION IN THE RADIO GALAXY 3C 120 , author=. ApJ , volume=. 2009 , publisher=
2009
-
[16]
ApJ , volume=
Connection between the accretion disk and jet in the radio galaxy 3C 111 , author=. ApJ , volume=. 2011 , publisher=
2011
-
[17]
ApJ , volume=
The gasdynamics of compact relativistic jets , author=. ApJ , volume=
-
[18]
ApJ , volume=
A Cutoff in the x-ray fluctuation power density spectrum of the Seyfert 1 galaxy NGC 3516 , author=. ApJ , volume=. 1999 , publisher=
1999
-
[19]
A&A , volume=
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=
2004
-
[20]
ApJ , volume=
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=
2025
-
[21]
MNRAS , volume=
Towards a unified model for black hole X-ray binary jets , author=. MNRAS , volume=. 2004 , publisher=
2004
-
[22]
MNRAS , volume=
Jets from black hole x--ray binaries: Testing, refining and extending empirical models for the coupling to x--rays , author=. MNRAS , volume=. 2009 , publisher=
2009
-
[23]
1991 , publisher=
A high signal-to-noise ratio composite quasar spectrum , author=. 1991 , publisher=
1991
-
[24]
A&A , volume=
On the location of the supermassive black hole in CTA 102 , author=. A&A , volume=. 2015 , publisher=
2015
-
[25]
MNRAS , volume=
A universal radio--X-ray correlation in low/hard state black hole binaries , author=. MNRAS , volume=. 2003 , publisher=
2003
-
[26]
Nature Astronomy , volume=
A wide and collimated radio jet in 3C84 on the scale of a few hundred gravitational radii , author=. Nature Astronomy , volume=. 2018 , publisher=
2018
-
[27]
MNRAS , volume=
On the nature of rapid X-ray variability in active galactic nuclei , author=. MNRAS , volume=. 1993 , publisher=
1993
-
[28]
A&A , volume=
X-ray variability of 104 active galactic nuclei-XMM-Newton power-spectrum density profiles , author=. A&A , volume=. 2012 , publisher=
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=
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=
2012
-
[31]
I-Structure function analysis and the relation between BL Lacertae objects and quasi-stellar objects , author=
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=
2012
-
[33]
A&A , volume=
Long-term radio variability of active galactic nuclei at 37 GHz , author=. A&A , volume=. 2025 , publisher=
2025
-
[34]
A&A , volume=
Active galactic nucleus time-variability analysis and its caveats , author=. A&A , volume=. 2025 , publisher=
2025
-
[35]
ApJ , volume=
Are the variations in quasar optical flux driven by thermal fluctuations? , author=. ApJ , volume=. 2009 , publisher=
2009
-
[36]
ppcor: an R package for a fast calculation to semi-partial correlation coefficients , author=. Commun. Stat. Appl. Methods , volume=
-
[37]
MNRAS , volume=
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=
1938
-
[38]
MNRAS , volume=
A transition from parabolic to conical shape as a common effect in nearby AGN jets , author=. MNRAS , volume=. 2020 , publisher=
2020
-
[39]
MNRAS , volume=
The core shift effect in the blazar 3C 454.3 , author=. MNRAS , volume=. 2014 , publisher=
2014
-
[40]
ApJ , volume=
X-ray variability of active galactic nuclei-A universal power spectrum with luminosity-dependent amplitude , author=. ApJ , volume=
-
[41]
ApJ , volume=
Constraining the limiting brightness temperature and Doppler factors for the largest sample of radio-bright blazars , author=. ApJ , volume=. 2018 , publisher=
2018
-
[42]
A&A , volume=
Synchrotron flaring behaviour of Cygnus X-3 during the February-March 1994 and September 2001 outbursts , author=. A&A , volume=. 2007 , publisher=
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=
2009
-
[44]
ApJ , volume=
Relativistic beaming and flux variability in active galactic nuclei , author=. ApJ , volume=. 2001 , publisher=
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=
1996
-
[46]
MOJAVE. XVII. Jet kinematics and parent population properties of relativistically beamed radio-loud blazars , author=. ApJ , volume=. 2019 , publisher=
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=
1999
-
[48]
ApJ , volume=
Asymptotic structure of Poynting-dominated jets , author=. ApJ , volume=. 2009 , publisher=
2009
-
[49]
MNRAS , volume=
A jet model for the fast IR variability of the black hole X-ray binary GX 339-4 , author=. MNRAS , volume=. 2018 , publisher=
2018
-
[50]
Nature , volume=
Observational evidence for the accretion-disk origin for a radio jet in an active galaxy , author=. Nature , volume=. 2002 , publisher=
2002
-
[51]
Nature , volume=
Active galactic nuclei as scaled-up Galactic black holes , author=. Nature , volume=. 2006 , publisher=
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=
2009
-
[53]
MNRAS , volume=
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=
2014
-
[54]
A&A , volume=
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=
2001
-
[55]
ApJL , volume=
As above, so below: exploiting mass scaling in black hole accretion to break degeneracies in spectral interpretation , author=. ApJL , volume=. 2015 , publisher=
2015
-
[56]
MNRAS , volume=
A Fundamental Plane of black hole activity , author=. MNRAS , volume=. 2003 , publisher=
2003
-
[57]
MNRAS , volume=
The accretion disc-jet connection in blazars , author=. MNRAS , volume=. 2019 , publisher=
2019
-
[58]
MNRAS , volume=
Physical parameters of active galactic nuclei derived from properties of the jet geometry transition region , author=. MNRAS , volume=. 2020 , publisher=
2020
-
[59]
Biometrika , volume=
Periodogram analysis and continuous spectra , author=. Biometrika , volume=. 1950 , publisher=
1950
-
[60]
MNRAS , volume=
On the fundamental dichotomy in the local radio-AGN population: accretion, evolution and host galaxy properties , author=. MNRAS , volume=. 2012 , publisher=
2012
-
[61]
MNRAS , volume=
The transition between BL Lac objects and flat spectrum radio quasars , author=. MNRAS , volume=. 2011 , publisher=
2011
-
[62]
Nature , volume=
The power of relativistic jets is larger than the luminosity of their accretion disks , author=. Nature , volume=. 2014 , publisher=
2014
-
[63]
ApJ , volume=
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=
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]
A&A , volume=
Doppler factors, Lorentz factors and viewing angles for quasars, BL Lacertae objects and radio galaxies , author=. A&A , volume=. 2009 , publisher=
2009
-
[66]
MNRAS , volume=
The jet--disc connection in AGN , author=. MNRAS , volume=. 2014 , publisher=
2014
-
[67]
MNRAS , volume=
On characterizing the variability properties of X-ray light curves from active galaxies , author=. MNRAS , volume=. 2003 , publisher=
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=
2007
-
[70]
ApJS , volume=
Five-year wilkinson microwave anisotropy probe* observations: cosmological interpretation , author=. ApJS , volume=. 2009 , publisher=
2009
-
[71]
MNRAS , volume=
Opacity, variability, and kinematics of AGN jets , author=. MNRAS , volume=. 2019 , publisher=
2019
-
[72]
A&A , volume=
Synchrotron flaring in the jet of 3C 279 , author=. A&A , volume=. 2006 , publisher=
2006
-
[73]
ApJ , volume=
Reverberation Measurements for 17 Quasars and theSize-Mass-Luminosity Relations in Active GalacticNuclei , author=. ApJ , volume=. 2000 , publisher=
2000
-
[74]
A&A , volume=
Optical and radio variability of the northern VHE gamma-ray emitting BL Lacertae objects , author=. A&A , volume=. 2016 , publisher=
2016
-
[75]
ApJ , volume=
The jet power, radio loudness, and black hole mass in radio-loud active galactic nuclei , author=. ApJ , volume=. 2006 , publisher=
2006
-
[76]
A&A , volume=
Identifying changing jets through their radio variability , author=. A&A , volume=. 2021 , publisher=
2021
-
[77]
ApJ , volume=
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=
1983
-
[78]
ApJ , volume=
The jet-disk connection and blazar unification , author=. ApJ , volume=. 2003 , publisher=
2003
-
[79]
arXiv preprint physics/0204033 , year=
1/f noise: a pedagogical review , author=. arXiv preprint physics/0204033 , year=
-
[80]
ApJ , volume=
Collimation of the Relativistic Jet in the Quasar 3C 273 , author=. ApJ , volume=. 2022 , publisher=
2022
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.