Minijets and Broken Stationarity in a Blazar : Novel Insights into the Origin of $γ$-ray Variability in CTA 102

This paper analyzes 18 years of Fermi-LAT data for the blazar CTA 102 to demonstrate that its 2017 giant flare marked a transition from a log-normal, frequently flaring state to a more stable plateau, a phenomenon explained by magnetic relaxation and successfully reproduced by a modified minijets-in-a-jet Monte Carlo simulation.

The Gist

Imagine a distant galaxy, 10 billion light-years away, acting like a cosmic lighthouse. This galaxy, called CTA 102, shoots a powerful beam of energy straight at Earth. Inside that beam, things are usually chaotic and unpredictable, flashing brightly and then dimming.

For a long time, scientists thought these flashes happened in a very specific, predictable mathematical pattern (called "log-normal"), similar to how a crowd of people might vary slightly in height around an average. But this new paper, looking at 18 years of data, says: "Actually, it's not that simple."

Here is the story of what the researchers found, explained without the heavy math.

1. The Great "Super-Flash" of 2017

For 18 years, the team watched this galaxy. Most of the time, it had small, frequent flickers. But in 2017, something massive happened. The galaxy didn't just flicker; it went into a "super-bright" mode, becoming 100 times brighter in high-energy gamma rays than usual. It was like a candle suddenly turning into a searchlight.

The researchers split their 18-year data into two groups:

• Before the Flash: The chaotic, flickering era.
• After the Flash: The calmer era that followed.

2. The "Skewness" Mystery

The scientists looked at the shape of the data. Imagine a hill of sand.

• Before 2017: The hill had a very long, spiky tail on one side. This meant there were many "outlier" events—sudden, massive bursts of energy that were rare but extreme.
• After 2017: That long, spiky tail got chopped off. The hill became more rounded and stable. The galaxy wasn't having those wild, extreme bursts as often anymore.

In simple terms: The galaxy went from being a wild, unpredictable rollercoaster to a steady, predictable train. The researchers call this a change in "skewness."

3. The "Mini-Jets" Analogy

How does a galaxy do this? The paper suggests the main jet isn't just one solid stream. Instead, think of the main jet as a giant highway. Inside this highway, there are thousands of tiny, fast cars called "minijets."

• The Rules of the Road: These minijets are zipping around in random directions.
• The Flash: A massive gamma-ray flash only happens when a huge number of these tiny cars all accidentally line up perfectly. They must be pointing toward a specific target (a cloud of gas near the galaxy) and pointing straight at Earth at the exact same time.
• The Result: When they line up, their combined speed and direction create a massive boost in brightness (like a lighthouse beam focusing perfectly). This is rare, which is why big flares are rare.

4. The "Magnetic Hair" Theory

So, why did the behavior change after 2017? The authors propose a theory involving magnetic fields.

Imagine the magnetic fields inside the jet are like a tangled ball of yarn.

• Before 2017: The yarn was a mess. The tangles were constantly snapping and reconnecting (like static electricity). Every time a piece of yarn snapped, it created a tiny "minijet" that zipped off. Because the yarn was so tangled, these snaps happened often and chaotically, creating those wild, long tails in the data.
• The 2017 Event: The "Super-Flash" was caused by a massive, violent untangling event. It was like someone violently shaking the ball of yarn, causing a huge burst of energy.
• After 2017: After that big shake, the yarn settled down. It became neat and ordered. Because the magnetic field was now smooth and organized, there were fewer "snaps" and fewer chaotic minijets. The galaxy became calmer, and the wild, long tails in the data disappeared.

5. The Computer Simulation

To prove this idea, the scientists built a computer model. They programmed thousands of virtual "minijets" to move randomly inside a jet.

• When they let the model run, it naturally produced flares that looked just like the real data.
• The model showed that the "messy" state (tangled magnetic fields) creates the wild, long-tailed distribution.
• The model showed that when the system "relaxes" (becomes ordered), the distribution smooths out, just like the real galaxy did after 2017.

The Bottom Line

The paper concludes that CTA 102 isn't just a random flickering light. It is a system where magnetic fields get tangled and then untangle.

• When the fields are tangled, the galaxy is wild, with frequent, extreme flares.
• When a massive event untangles the fields (like the 2017 flare), the galaxy settles into a calmer, more stable state.

The researchers used a special mathematical formula (a "Modified Log-Normal Power-Law") to describe this transition, proving that their "minijet" theory fits the real-world data perfectly. It's a story of cosmic chaos turning into cosmic order.

Technical Summary

Technical Summary: Minijets and Broken Stationarity in Blazar CTA 102

Problem Statement
High-energy blazar light curves have historically been modeled as following a log-normal flux distribution, a signature attributed to multiplicative processes within the jet or connections to accretion. However, the physical origin of variability in specific sources, particularly regarding the transition between different activity states, remains an open question. This study investigates the flat spectrum radio quasar (FSRQ) CTA 102 ($z=1.037$), which exhibited a massive $\gamma$-ray flare in 2017 (a $\sim$100-fold flux jump) alongside X-ray flaring. The central problem addressed is whether the pre-flare and post-flare states of CTA 102 adhere to the standard log-normal distribution and, if not, what physical mechanisms explain the observed statistical deviations and the transition between these states.

Methodology
The authors analyzed an 18-year archival Fermi-LAT light curve (0.1–100 GeV) of CTA 102, spanning from August 2008 to November 2025. Data points with a test statistic (TS) < 25 were rejected. The light curve was divided into "pre-flare" and "post-flare" epochs based on the 2017 event (July 2016–August 2017).

The analysis employed several statistical techniques:

1. Distribution Analysis: Kolmogorov-Smirnov (KS) tests were used to compare pre- and post-flare distributions.
2. Skewness Quantification: The authors computed skewness for 1-year bins (100 points) and calculated cumulative and "rolling" skewness to track temporal evolution. Uncertainties were derived via bootstrapping (1,000 samples).
3. Statistical Testing: A Mann-Whitney U test compared pre- and post-flare skewness distributions, and a t-test evaluated the hypothesis of a return to Gaussianity.
4. Modeling: A modified "minijets-in-a-jet" Monte Carlo simulation was developed. This model assumes 50 randomly oriented minijets within the jet, interacting with Broad Line Region (BLR) photons via External Compton (EC) scattering. The simulation varied minijet Lorentz factors (drawn from a power law) and viewing angles to reproduce flux distributions.
5. Distribution Fitting: Both observed and simulated data were fitted using a Modified Log-Normal Power-Law (MLP) distribution (Basu et al. 2015), characterized by a parameter $\alpha = \Delta/\eta$ (decay rate over growth rate).

Key Results

• Non-Log-Normality: Contrary to the standard assumption, neither the pre-flare nor the post-flare GeV light curves follow a strictly log-normal distribution. Both exhibit extended tails.
• Statistically Significant State Transition: A Mann-Whitney test confirmed a highly significant difference ($p \sim 10^{-9}$) between pre- and post-flare skewness. The pre-flare state exhibited higher skewness (mean $\sim$0.7) with a long tail indicative of frequent flaring. The post-flare state showed reduced skewness (mean $\sim$0.5) and a "transfer" of flux from the tail toward the median, indicating a transition to a more plateaued state with occasional flaring.
• Photon Index Stability: Despite the drastic change in flux distribution and skewness, the $\gamma$-ray photon indices remained virtually unchanged between the two states (KS statistic $\sim$0.04), suggesting that the radiative conditions (e.g., electron energy distribution) did not fundamentally alter, but rather the geometry or magnetic topology did.
• Simulation Success: The modified minijets model successfully reproduced the GeV flares and the non-log-normal flux distributions. The simulation demonstrated that massive flares occur only when a maximum number of minijets align toward both the BLR and the line of sight, a low-probability event.
• MLP Fitting: The MLP distribution provided a robust fit for both observed and simulated data. The parameter $\alpha$ was found to be higher in the post-flare state ($1.21 \pm 0.09$) compared to the pre-flare state ($1.12 \pm 0.10$), differing at a $\gtrsim 1\sigma$ level.

Physical Interpretation and Mechanism
The authors propose that the observed transition is driven by magnetic relaxation.

• Pre-Flare State: Characterized by a tangled, disordered magnetic field. Frequent small-scale magnetic reconnection events injected energy, creating minijets that occasionally aligned to produce flares. This resulted in a high injection rate ($\eta$) and a long stopping timescale, manifesting as a high-skewness, long-tailed distribution.
• The 2017 Flare: Attributed to a massive reconnection event (or a cluster of events) that released a burst of energy.
• Post-Flare State: The reconnection event "relaxed" the magnetic field into a more ordered, lower-energy configuration. This ordering reduced the rate of reconnection events (lowering $\eta$) and the number of active minijets. Consequently, the flux distribution stabilized, the tail diminished, and the skewness decreased. The constancy of the photon index supports the view that the change was geometric/magnetohydrodynamic rather than a change in the fundamental particle acceleration mechanism.

Significance and Claims
The paper claims to provide novel insights into the origin of $\gamma$-ray variability in blazars by demonstrating that:

1. Broken Stationarity: Blazar light curves can undergo statistically significant transitions in their flux distribution properties (specifically skewness) without changing their spectral indices.
2. Minijets as the Unifying Mechanism: A modified minijets-in-a-jet model, incorporating random orientations and alignment with the BLR, can reproduce both the rare, massive flares and the specific non-log-normal flux distributions observed in CTA 102.
3. Magnetic Relaxation: The transition from a high-variance, flaring state to a steady state can be physically explained by magnetic relaxation following a major reconnection event, offering a mechanism for the "broken stationarity" observed in the data.
4. Model Validation: The successful fit of the MLP distribution to both real data and the minijet simulation cements the link between minijet dynamics and the statistical nature of blazar variability, suggesting that the standard log-normal model is insufficient for sources like CTA 102.