🔭 astrophysics

Connecting the long-term variability behaviour of active galactic nuclei to their central engines

This study analyzes long-term radio variability in 54 active galactic nuclei to demonstrate that characteristic timescales derived from power spectral density analysis correspond to the mean duration of individual flares and correlate with central engine parameters, specifically the ratio of black hole mass to normalized mass accretion rate.

Original authors: Sofia Kankkunen, Merja Tornikoski, Talvikki Hovatta

Published 2026-05-07

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Sofia Kankkunen, Merja Tornikoski, Talvikki Hovatta

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, chaotic orchestra. At the center of many galaxies, there are supermassive black holes acting as the conductors. These conductors don't just sit there; they launch massive, high-speed jets of energy (like powerful water hoses) into space. Sometimes, these jets hiccup, creating bright flashes of light called "flares."

This paper is like a detective story where astronomers tried to figure out why these jets hiccup the way they do and what the "rhythm" of these hiccups tells us about the black hole itself.

Here is the breakdown of their investigation in simple terms:

1. The Mystery: The "Beat" of the Black Hole

Astronomers have long known that these black holes vary in brightness over time. If you look at the data mathematically (using something called a Power Spectral Density, or PSD), it looks like a song with a specific "beat" or rhythm. There is a point in the rhythm where the pattern changes (a "bend").

The big question was: What physical thing in the jet creates this specific beat?

Is it how long a single flash (flare) lasts?
Is it the time gap between one flash and the next?

2. The Method: Breaking Down the Song

The researchers took 54 different galaxies and looked at their radio light curves (graphs showing how bright they get over decades). They used a special computer program to break these long, messy graphs down into individual "flares," like separating individual notes from a complex melody.

They measured two things for the brightest flares:

Duration: How long the flash lasted (from the start of the rise to the end of the fade).
Separation: How much time passed between the peak of one flash and the peak of the next.

3. The Discovery: Matching the Rhythm

They compared their measurements to the "beat" (the PSD bend) they found in previous studies.

The Result: They found that the duration of the brightest flares matched the "beat" of the black hole almost perfectly. It's as if the length of the flash is the physical reason for that specific rhythm in the data.
The Gap: The time between flashes also matched somewhat, but the length of the flash was the stronger match.

4. Connecting to the Engine: The Black Hole's Size and Appetite

Once they knew which part of the light curve represented the "beat," they asked: Does this beat change based on the size of the black hole or how much food (gas) it is eating?

The Finding: They found a connection. The "beat" (the timescale of the flares) seems to depend on the mass of the black hole divided by how fast it is eating.
The Analogy: Think of it like a car engine. A massive truck engine (a heavy black hole) that is sipping fuel slowly (low accretion rate) might have a different idle rhythm than a smaller engine revving its fuel intake. The rhythm of the flares tells us about the "size of the engine" and "how hungry it is."

5. The "Core" of the Jet: How Big is the Flash?

The researchers also looked at how fast the flares rose to their peak brightness. They used this speed to estimate the physical size of the region where the light is coming from (the "emission region").

They tried to see if the size of this region was linked to:

How fast the jet is moving (Lorentz factor).
How big the black hole is.
How bright the accretion disk (the food bowl) is.

The Result: They found some weak hints of a connection, but there was a catch. It turned out that distance (redshift) was the real culprit.

The Analogy: Imagine looking at a group of runners. The ones far away (high redshift) seem to run faster and have bigger strides. But when you correct for the distance, the connection disappears. The researchers realized that the apparent links between the jet's size and the black hole's properties were mostly because the distant sources just looked different due to how far away they were, not because of a fundamental physics rule.

Summary of the Takeaways

The Rhythm Match: The length of the brightest flashes in a galaxy's jet matches the "beat" found in the mathematical analysis of the light curve.
The Engine Link: This rhythm is likely tied to the central engine's properties—specifically, the black hole's mass and how fast it is consuming matter.
The Distance Trap: While there seemed to be links between the size of the flash and the speed of the jet, these links were mostly an illusion caused by the fact that the most distant galaxies in the sample behaved differently. When you account for distance, those specific links become much weaker or disappear.

In short: The astronomers successfully identified that the "length of the flash" is the key to understanding the black hole's rhythm, and that rhythm is a signature of the black hole's mass and appetite. However, they also learned to be careful not to confuse "looking different because you are far away" with "actually being different."

Technical Summary: Connecting the long-term variability behaviour of active galactic nuclei to their central engines

Problem Statement
Active Galactic Nuclei (AGNs) exhibit multi-wavelength variability driven by relativistic jets launched by supermassive black holes. While Power Spectral Density (PSD) analysis is a standard method for identifying characteristic variability timescales (specifically the low-frequency bend in the red noise spectrum), the physical interpretation of these timescales remains debated. Previous studies on X-ray binaries and AGNs suggest PSD bends scale with black hole mass ( $M_{BH}$ ) and accretion rate ( $\dot{M}$ ), but the specific light-curve feature corresponding to the radio PSD bend is unclear. It is uncertain whether the bend frequency inverse corresponds to the duration of individual flares or the temporal separation between them. Furthermore, the relationship between these variability timescales and the physical size of the emission region (often associated with the jet transition region) and central engine parameters requires further investigation using a larger sample than previously available.

Methodology
The authors analyzed long-term radio light curves of 54 AGNs monitored at the Aalto University Metsähovi Radio Observatory (MRO) at 37 GHz. The sample included 11 sources with previously constrained PSD timescales (from Kankkunen et al. 2025, "Paper I") and 43 additional sources (including 3 borderline cases) selected for high data density and prominent flares.

The core methodology involved:

Light-Curve Decomposition: The Magnetron code was employed to decompose light curves into individual flares using Bayesian hierarchical modelling and diffusive nested sampling. An Ornstein-Uhlenbeck (OU) process was used to model underlying background variability, preventing the misidentification of minor turbulent variations as flares.
Timescale Extraction: For each source, the authors calculated the mean duration and mean peak separation of the brightest flares (defined as those exceeding the mean amplitude in a posterior draw). Rise times were also extracted for the brightest flares (amplitude > 50% of the maximum) to serve as proxies for the emission region size.
Parameter Estimation:
- Doppler Factors ( $\delta_{var}$ ): Estimated using variability brightness temperatures derived from flare rise times and amplitudes, compared against a Gaussian distribution of intrinsic brightness temperatures.
- Lorentz Factors ( $\Gamma$ ): Calculated using $\delta_{var}$ and apparent jet speeds ( $\beta_{app}$ ) from the MOJAVE survey.
- Central Engine Parameters: Black hole masses ( $M_{BH}$ ) and accretion disk luminosities ( $L_{acc}$ ) were estimated using broad-line luminosities (H $\beta$ , MgII, or [OIII]) and full width at half maximum (FWHM) values. Normalized mass accretion rates ( $\dot{m}$ ) were derived assuming a radiative efficiency of $\eta = 0.4$ .
Statistical Analysis: Kendall rank correlation tests were performed to assess associations between the derived timescales (flare durations, separations, and rise times) and physical parameters ( $M_{BH}$ , $\dot{m}$ , $L_{acc}$ , $\Gamma$ ). Partial correlation analyses were conducted to account for potential dependencies on redshift ( $z$ ).

Key Results

PSD Timescale vs. Light Curve Features: The inverse of the PSD bend frequency corresponds best to the mean duration of the brightest flares. While mean flare separations also correlate with the PSD timescale, the flare durations show a stronger one-to-one agreement, particularly for the 11 constrained sources.
Connection to Central Engine: Using flare durations and separations as proxies for the PSD timescale, the authors found a moderate positive correlation with the ratio of black hole mass to normalized mass accretion rate ( $M_{BH}/\dot{m}$ ). This correlation is significant in the observer frame but weakens or becomes non-significant when corrected to the emission frame. The correlation appears driven primarily by $M_{BH}$ , though the role of $\dot{m}$ cannot be definitively excluded due to sample size limitations.
Emission Region Size: The rise times of the brightest flares, used as proxies for the jet transition region width, showed weak tentative correlations with $M_{BH}$ , $L_{acc}$ , and $\Gamma$ . However, these associations were found to be driven by a common dependency on redshift. When redshift was accounted for via partial correlation, the significance of the relationships between rise times and central engine parameters disappeared.
Redshift Dependence: A significant positive correlation was found between redshift and emission-frame flare durations, while a negative correlation existed between redshift and emission-frame flare separations (implying a higher ejection rate for higher redshift sources). These opposing trends may explain the lack of correlation between observed flare durations/separations and Doppler factors.

Significance and Claims
The paper claims that the characteristic timescales derived from radio PSD bends are physically associated with the central engine of AGNs, specifically scaling with $M_{BH}/\dot{m}$ . This supports the hypothesis that the central engine modulates the jet variability, similar to findings in the X-ray domain. The study suggests that the mean duration of the brightest flares is the most appropriate light-curve proxy for the PSD bend timescale.

The authors modestly note that while correlations with central engine parameters exist, they are often weak and potentially confounded by redshift effects and sample biases (e.g., the selection of sources with discernible flares). The results regarding the emission region size and Lorentz factor are described as "weak tentative correlations" that may not hold once redshift dependencies are removed. Consequently, the paper concludes that while a connection to the central engine is suggested, larger samples and more accurate parameter estimates are required to definitively characterize these relationships. The work does not propose new experimental setups but rather validates the use of light-curve decomposition as a method to extend PSD-based analyses to a broader sample of AGNs.