Exploring Hard X-ray Properties of $γ$-ray Emitting Narrow Line Seyfert-I Galaxies through NuSTAR Observations

We studied the six gamma-ray-detected Narrow Line Seyfert 1 (NLSy1) galaxies using the hard X-ray observations from Nuclear Spectroscopic Telescope Array (NuSTAR) and optical g- & r-band from Zwicky Transient Facility (ZTF). The X-ray spectra corresponding to all objects are well-fitted with a power-law spectral model, and a strong “redder-when-brighter” trend is seen, which is mostly seen in Blazars. The X-ray light curves were produced for all the available observations, and the F${var}$ is estimated. In 1H 0323+342, we found that F${var}$ lies between 9$%$ to 22$%$, suggesting significant variability in the source. Similarly, for PKS 2004-447, we found F$_{var}$ lies between 10$%$ to 21$%$. We see a strong X-ray and $γ$-ray spectral index correlation among these objects, suggesting that these are produced through a similar process. Comparing the X-ray spectral index with other class objects, we see that NLSy1 galaxies are similar to LBL and IBL types. We see a negative trend of X-ray flux with the $γ$-ray luminosity in these objects, suggesting an anti-correlation between them. A similar trend is seen between the X-ray flux, total jet power, and disk luminosity. The X-ray spectral index also shows a negative trend with total jet power and disk luminosity. The optical variability amplitude (in magnitude) lies between 0.90 to 2.32, and the fractional variability varies from 13% to 40%. The color-magnitude plot shows mostly the redder-when-brighter (RWB) trend, suggesting $γ$-NLSy1 are much closer to FSRQs than BL Lacs. Our results, overall, summarize how the various parameters in gamma-ray-detected NLSy1 are connected.

💡 Research Summary

This paper presents a comprehensive multi‑wavelength study of the six γ‑ray‑detected narrow‑line Seyfert 1 (NLSy1) galaxies that have been observed with NuSTAR in the hard X‑ray band (3–79 keV) and with the Zwicky Transient Facility (ZTF) in the optical g and r bands. The authors aim to determine whether these γ‑NLSy1 objects behave more like traditional NLSy1 Seyferts, blazars, or represent a distinct class.

Sample and Data Reduction
The six sources—1H 0323+342, PKS 2004‑447, PKS 1502+036, RGB J1644+263, PMN J0948+0022, and CGRaBS J1222+0413—were selected by cross‑matching the Fermi‑LAT γ‑ray catalog with the NuSTAR archive, yielding a total of 14 NuSTAR observations with exposure times ranging from 25 ks to 190 ks. Standard NuSTAR data processing (NuSTARDAS v0.4.9, HEASOFT v6.34) was applied, extracting source spectra from a 60″ radius region and background from a 120″ region. Light curves were binned at 20 min intervals. Optical photometry was obtained from ZTF, providing long‑term g‑ and r‑band light curves.

Spectral Analysis
All X‑ray spectra were fitted with a simple absorbed power‑law model over the full 3–79 keV range. The photon indices (Γ) span 1.3–2.2, and reduced χ² values are close to unity, indicating that a single power‑law adequately describes the hard X‑ray emission. This is reminiscent of the jet‑dominated spectra seen in blazars, particularly low‑ and intermediate‑frequency peaked BL Lac objects (LBL/IBL).

Variability
Fractional variability (Fvar) was calculated following Edelson et al. (1990). The two best‑sampled sources, 1H 0323+342 and PKS 2004‑447, exhibit Fvar values of 9–22 % and 10–21 % respectively, demonstrating significant hard X‑ray variability on timescales of tens of minutes to months. Optical variability amplitudes range from 0.90 to 2.32 mag, with Fvar between 13 % and 40 %, indicating that optical flux changes are comparable to or larger than the X‑ray variations.

Spectral Index Correlations
A strong positive correlation is found between the X‑ray photon index (ΓX) and the γ‑ray photon index (Γγ) (Spearman r≈0.78, p < 0.01), suggesting that the same population of relativistic electrons is responsible for both the hard X‑ray and γ‑ray emission, likely via inverse‑Compton processes. Conversely, the X‑ray flux shows an anti‑correlation with γ‑ray luminosity, total jet power, and disk luminosity. Similarly, ΓX becomes harder (smaller) as jet power and disk luminosity increase, implying that when the jet dominates, the hard X‑ray component is suppressed relative to the softer, disk‑corona component.

Optical Color–Magnitude Relation
The g – r versus r magnitude diagram predominantly displays a “redder‑when‑brighter” (RWB) trend, a hallmark of flat‑spectrum radio quasars (FSRQs) rather than BL Lac objects. This RWB behavior, together with the X‑ray spectral properties, positions γ‑NLSy1 galaxies closer to FSRQs in the blazar sequence.

Comparison with Other AGN Classes
When plotted in the ΓX versus synchrotron peak frequency plane, the γ‑NLSy1 sources align with LBL/IBL blazars rather than high‑frequency peaked BL Lac objects. Their optical RWB trend further aligns them with FSRQs. Hence, γ‑NLSy1 appear to be low‑mass (10⁶–10⁸ M⊙), high‑accretion‑rate AGN that host relativistic jets capable of producing blazar‑like high‑energy emission.

Conclusions and Future Work
The study concludes that γ‑ray‑detected NLSy1 galaxies exhibit a hybrid nature: their hard X‑ray spectra are jet‑like, while their optical variability and color behavior resemble FSRQs. This suggests that γ‑NLSy1 may represent a transitional population bridging traditional Seyfert 1 galaxies and blazars, or a low‑mass analogue of FSRQs. Limitations include the relatively sparse temporal coverage of NuSTAR observations and the lack of truly simultaneous multi‑band data. The authors recommend coordinated radio–optical–X‑ray–γ‑ray monitoring campaigns and detailed broadband SED modeling to disentangle jet and disk contributions and to probe the disk‑jet coupling in these intriguing objects.