A morphological search for stable year-scale modulation in gamma-ray blazars

The information presented in both documents is highly consistent, but by synthesizing them, I can construct a comprehensive and detailed summary suitable for rigorous scientific documentation.

Here is the detailed synthesis:

Comprehensive Summary of Morphological Framework for Year-Scale Modulation in Blazars

This research introduces and calibrates a robust morphological framework designed to characterize the temporal behavior of year-scale quasi-periodic oscillations (QPOs) observed in gamma-ray blazars. The framework was applied to a sample of 120 Fermi-LAT light curves spanning approximately 18.1 years, utilizing Singular Spectrum Analysis (SSA) to select the oscillatory mode and its associated period, alongside Weighted Wavelet Z-transform (WWZ) ridge tracking on the detrended light curve to quantify key morphological parameters: coherence (Q), period drift (epsilon), cycle-to-cycle volatility (sigma P), and persistence.

Methodology and Calibration

The framework's reliability is established through rigorous calibration procedures. The classification system (Rank I, II, and III) was calibrated using 10,000 red-noise simulations and approximately 20,000 injected signals. Crucially, the quality selection for the data ranking process was fixed based on these simulations before applying them to the actual observational data.

Key Methodological Notes:

1. Measurement Technique: The morphological parameters are measured on the detrended light curve. This is a critical methodological choice, as tracing these features directly onto an SSA reconstruction tends to suppress the very drift and jitter that this framework is designed to measure (suppression factors of 3.4 and 44 were noted).

2. Calibration Necessity: Each morphological parameter must be calibrated against injected signals before any physical interpretation is made. Flux-dependent censoring alone can move about one-third of stationary signals out of the Rank I classification, underscoring the need for injection tests to validate the framework's response to temporal evolution.

3. Noise Model Dependence: The final population results—specifically the incidence limits on stable modulation—are inherently dependent on the assumed underlying noise model (e.g., red-noise background). A more complex noise model could alter false-alarm probabilities, leading to tighter or looser upper limits on genuine physical phenomena.

Classification Ranks and Morphological Definitions

The framework classifies sources based purely on their observed morphology, which describes the behavior rather than definitively identifying the underlying physical mechanism. The classification is defined by specific quantitative thresholds for coherence (Q), number of active cycles (N active), period evolution (epsilon), and volatility (sigma P).

Rank I: Stable Clock (Scenario A)

Sources in Rank I exhibit characteristics consistent with a rigid, stable clock mechanism. They are defined by:

• High Coherence: Q > 3.

• Stationarity: The period evolution is negligible (epsilon + epsilon < 0.02 yr-1).

• Low Volatility: sigma P + sigma P 0.90.

Rank II: Coherent, Evolving Modulation (Scenario B)

Sources in Rank II show coherent modulation where the period cannot be established as stationary, failing the rigidity test of Rank I. They are characterized by:

• Failing the stationarity threshold (epsilon + epsilon 0.02 yr-1) OR exhibiting intermediate volatility (0.06 sigma P + sigma P 0.60).

Rank III: Residual Candidates

Sources in Rank III fail to meet the criteria for Rank I or II, typically due to insufficient coherence, a period exceeding the established ceiling, excessive plasma-level volatility, or insufficient persistence.

Key Constraints and Population Results

The analysis imposes strict constraints on these classifications derived from the observational baseline (18.1 years):

• Period Ceiling (P ceil): The requirement for Rank I and II sources necessitates a minimum number of observable cycles (N active > 3). For the 18.1-year baseline, this imposes a hard maximum period ceiling of P ceil 3.1 years (derived from Equation 18).

• Completeness and Limits: Injection–recovery tests demonstrate that Rank I completeness ranges from 63% to 97% at an amplitude Signal-to-Noise Ratio (SNR) of 1, reaching 100% at SNR = 2. Drifting modulations are classified as Rank II in 92% of trials.

• Incidence Limit: Combining the observed counts with the expected false positive rates from red-noise simulations (where controls reach Rank I in 1.1% and Rank II in 2.3%), the framework establishes stringent upper limits on stable, approximately sinusoidal 1–3 yr modulation:

• At **SNR 1 **: At most 7% of the sample hosts stable modulation.

• At **SNR 2 **: At most 5% of the sample hosts stable modulation (95% confidence).

Physical Scenarios and Implications

The morphological framework distinguishes between three primary physical scenarios:

1. Scenario A (Geometric): Characterized by a rigid driver, expected to produce smooth, coherent, and long-lived modulation. This corresponds to Rank I.

2. Scenario B (Geometric): Involves an evolving geometric configuration (e.g., Lense–Thirring precession), which is expected to show systematic secular evolution (epsilon > 0) while maintaining coherence. This corresponds to Rank II.

3. Scenario C (Plasma-driven): Manifests as enhanced cycle-to-cycle variability of the characteristic period (sigma P 0.1).

Implications for Gravitational Wave Searches:

• Rank I sources (e.g., PG 1553+113 and B2 1215+30) are consistent with binary orbital motion but do not uniquely identify it; breaking this degeneracy requires longer baselines or independent evidence.

• The gravitational-wave decay timescale (tau GW) for a compact Supermassive Black Hole Binary (SMBBH) at fiducial values is tau GW 7 times 10 cubed years, which is two orders of magnitude shorter than the precision achievable by these measurements.

• Rank I sources are valuable as sky-position and mass priors for continuous-wave searches in Pulsar Timing Arrays (PTAs), but they cannot be directly matched to resolved PTA signals due to the short decay timescale.

Conclusion:

The framework successfully distinguishes between persistent, stationary modulation (Rank I) and coherent variability whose period may evolve or fluctuate (Rank II). While stable year-scale modulation is statistically rare among gamma-ray blazars, the identified Rank I sources remain primary candidates for multiwavelength follow-up and serve as crucial priors for continuous gravitational-wave searches. The framework prioritizes smooth, sinusoidal modulation due to the Morlet basis and 30-day binning, while acknowledging that burst-dominated modulation requires different template approaches.

1. The improved AI system can perform morphological classification of year-scale QPOs in blazar gamma-ray light curves using a framework that quantifies how coherent, persistent and stable a modulation is, distinguishing between Rank I (stable clock), Rank II (evolving period), and Rank III (residual candidates)."

2. The system can identify physical scenarios by mapping the morphology to underlying physics; specifically, it can differentiate between geometric scenarios like Lense–Thirring precession (Scenario B) and plasma-driven scenarios like magnetohydrodynamic instabilities (Scenario C) based on measured parameters such as the fractional drift rate ε and volatility σln P.

3. The AI can constrain the incidence of stable modulation by calculating a robust upper limit, stating that for a stationary period with SNR ≥ 2, the limit is at most 7% of the sample, which is derived from calibrating the framework against both red-noise controls and injected signals.

4. The system can perform source-level physical diagnostics by assessing whether a candidate satisfies multiple criteria simultaneously, such as checking if a Rank I candidate meets conditions: "High coherence, Q − ∆Q > 3, a cycle-count gate, Nactive > 3, dynamical rigidity (ε + ∆ε < 0.02 yr−1) and low volatility (σln P + ∆σln P < 0.06)."

5. The system can provide physically motivated pre-selection for multiwavelength follow-up by calculating a quantified probability, conditional on the noise model, that each is a red-noise interloper, which helps prioritize candidates consistent with binary orbital motion (Rank I).

Year-scale quasi-periodic oscillations (QPOs) in blazar gamma-ray light curves have been proposed as signatures of supermassive binary black holes and precessing jets. Red noise can mimic periodicity, and detection alone does not establish persistence or period stability. We develop and calibrate a morphological framework to characterize year-scale modulation and constrain the incidence of stable modulation in blazars. We analyze 120 Fermi-LAT light curves spanning 18.1 yr. Singular spectrum analysis selects the oscillatory mode and period, while weighted wavelet Z-transform ridge tracking measures persistence, period drift, and cycle-to-cycle volatility. Flux randomization, 10,000 red-noise simulations, and about 20,000 injections calibrate the classification. Of 97 quality-selected sources, two show persistent, coherent modulation consistent with a stationary period (Rank I: PG 1553+113 and B2 1215+30), six show coherent modulation whose stationarity cannot be established (Rank II), and 89 form Rank III. Neither Rank I nor Rank II significantly exceeds red-noise expectations. Stationary 1-3 yr injections reach Rank I in 63-97% of trials at signal-to-noise ratio (SNR) = 1, depending on period, and 100% at SNR = 2; 92% of drifting injections reach Rank II. These results limit stable, approximately sinusoidal 1-3 yr modulation to at most 7% of the sample at SNR >= 1 and 5% at SNR >= 2 (95% confidence). PG 1553+113 has the lowest false-alarm probability, 4.7% after trials. Stable year-scale modulation is uncommon, and no source is established as periodic. The framework distinguishes persistent, stationary modulation from coherent variability whose period may evolve or fluctuate, but not the physical mechanism. PG 1553+113 and B2 1215+30 are candidates for monitoring and continuous gravitational-wave searches, although their stability does not establish a binary origin.