[論文レビュー] ZTF Monitoring of $γ-$ray emitting Narrow Line Seyfert 1 Galaxies

Aman Kumar, Suvas Chandra Chaudhary+3 more|arXiv (Cornell University)|Feb 25, 2026

Astrophysics and Cosmic Phenomena被引用数 0

ひとこと要約

The paper analyzes multi-year ZTF optical light curves of γ-ray emitting NLSy1 galaxies using DRW and CARMA models to map optical variability timescales to disk and jet processes, probing disc–jet coupling.

ABSTRACT

The $γ$-ray-emitting narrow-line Seyfert-I ($γ$-NLSy1) are among the most interesting systems for studying disk-jet coupling. The soft X-ray properties of these systems suggest the presence of a disc component, which peaks in the optical/UV regime, in addition to the active jet. In this work, we investigate the optical emission from $γ$-NLSy1 using long-term Zwicky Transient Facility (ZTF) observations and discussed in the context of blazars. We have reported the long-term flux and color variability in the g- and r-bands. The fractional variability ($F_{ m var}$) goes as high as 72\%, with a mean value of 23\%, while the amplitude of variability ($ψ$) values range from 0.24 to 3.20, which is consistent with the long-term Swift-UVOT variability studies. The color-magnitude diagrams exhibit an RWB or BWB trend similar to that of blazars. The $t_{ m var}$ suggests an emitting region size of $10^{15-17}$ cm, aligned with emissions coming from the inner accretion disk or base of the jet. The PSD analysis using both DRW and CARMA modeling exhibits a characteristic break timescale of a few days to hundreds of days, which is likely linked to fundamental physical timescales in the system, such as thermal or viscous timescales in the accretion disk or timescales for acceleration and energy dissipation in the jet. The existence of these timescales acts as another signature of the disc-jet connection. These time scales are correlated with black hole mass, and the relation is consistent with previous studies.

研究の動機と目的

• Investigate how optical variability in γ-ray emitting NLSy1 galaxies arises from disk versus jet processes.
• Quantify long-term and intraday variability using ZTF g- and r-band data.
• Characterize variability timescales and structure via PSD analysis and time-domain modeling.
• Explore color variability and interband correlations to infer emission components (disk vs jet).

提案手法

• Compute fractional variability (F_var) and variability amplitude (ψ) from long-term ZTF light curves.
• Model light curves with DRW (CAR(1)) and higher-order CARMA(p,q) processes to derive PSDs.
• Use AIC to select optimal CARMA model and obtain analytical and MCMC-based PSD uncertainties.
• Estimate break timescales from CARMA PSD fits and compare across g- and r-bands.
• Apply Discrete Correlation Function (DCF/ZDCF) to assess interband lags between g- and r-band light curves.
• Analyze flux distributions with normal, lognormal, or bimodal Gaussian mixtures and assess bimodality via BIC and BC.
• Examine flux–rms relations to infer multiplicative versus additive variability processes.

Figure 1: ZTF lightcurve for $\gamma$ -NLSY1 J142105+385522 in $g$ and $r$ band. Other lightcurves are shown in Appendix A .

実験結果

リサーチクエスチョン

• RQ1 What are the characteristic optical variability timescales in γ-ray emitting NLSy1 galaxies as traced by ZTF data?
• RQ2 Do DRW or higher-order CARMA models better describe the optical variability in these sources?
• RQ3 What do interband (g vs r) correlations and lags reveal about the origin of variability (disk vs jet)?
• RQ4 How do color–magnitude trends (RWB vs BWB) relate to the relative disk and jet contributions?
• RQ5 Are the observed PSD breaks and timescales correlated with black hole mass or jet/disk properties?

主な発見

• F_var in the r-band ranges from 1.75% to 53.49% with a mean around 23%, while in the g-band it ranges from 1.03% to 72.59%, indicating significant optical variability.
• ψ ranges from 0.21 to 2.86 in r-band and 0.25 to 3.20 in g-band, showing strong peak-to-peak variability similar to jet-dominated sources.
• g- and r-band damping timescales (τ_damping) span from minutes to 1100 days (r-band mean ~131 days) and from an hour to 838 days (g-band mean ~83 days), implying wavelength-dependent variability.
• CARMA modeling favors higher-order processes (p ≥ 2, q up to 3; often CARMA(4,2) or CARMA(4,3)) over simple CAR(1), revealing multiple characteristic timescales and complex PSD shapes.
• PSD analysis shows breaks corresponding to physical timescales; double power-law fits yield a broad range of break timescales from ~0.03 days to >250 days across sources and bands.
• Discrete Correlation Function analyses reveal interband lags (negative mean lag for g behind r; some sources show near-zero or measurable lags) indicating mixed jet-disk variability and possible single-zone jet origin in some cases.]
• table_headers: []
• table_rows: []

Figure 2: Colour Magnitude Plot of J142105+385522. Other colour-magnitude plots are shown in Appendix A .