Constraining black hole spin in PG 1535+547 amidst complex multi-layered absorption

We present a spectroscopic analysis of XMM-Newton and NuSTAR observations of the ‘complex’ NLS1 PG 1535+547 at redshift $z=0.038$. These observations span three epochs: 2002 and 2006 with XMM-Newton alone, covering the $0.3-10$ keV energy range, and a coordinated XMM-Newton and NuSTAR observation in 2016, covering the $0.3-60$ keV energy range. The X-ray spectra across all epochs exhibit both neutral and ionized absorption, along with reflection features from the accretion disc, including a prominent Compton hump in the broadband data. Notably, the spectral shape varies across epochs. Our analysis suggests this variability is attributed to changes in both line-of-sight absorption and the intrinsic emission from PG 1535+547. The source is obscured by multiple layers of partially and/or fully covering neutral and ionized absorbers, with neutral column densities ranging from undetectable levels in the least obscured phase to $\sim0.3-5\times10^{23}\mathrm{cm^{-2}}$ in the most obscured phase. A clear warm absorber is revealed during the least obscured phase. The continuum remains fairly consistent ($Γ\approx 2.2\pm0.1$) during the first two observations, followed by a substantial flux decrease (by a factor of $\sim7$ in the $2-10$ keV band) in 2016 compared to 2006. The 2016 data indicates the source is in a reflection-dominated state during this epoch, with a reflection fraction of $R>7$ and an X-ray source located at a height $\leq 1.72r_g$. Simultaneous fitting of the multi-epoch data suggests a rapidly rotating black hole with a spin parameter, $a>0.99$. These findings imply that strong light-bending effects may account for the observed continuum flux reduction.

💡 Research Summary

This paper presents a comprehensive broadband X‑ray spectral analysis of the narrow‑line Seyfert 1 galaxy PG 1535+547 (z = 0.038) using XMM‑Newton data from 2002 and 2006 and a coordinated XMM‑Newton + NuSTAR observation from 2016. The three epochs together provide coverage from 0.3 keV up to 60 keV, allowing the authors to disentangle the contributions of neutral and ionised absorption, distant reflection, and relativistic disc reflection.

The authors reduced the data with standard SAS and nustardas pipelines, extracted source and background spectra, and combined the MOS1 and MOS2 data for each epoch. Spectra were binned to a minimum signal‑to‑noise of 5 per bin (XMM) and 3 per bin (NuSTAR). Visual inspection shows a deep soft‑X‑ray dip (0.7–2 keV) in the 2006 epoch, indicative of a warm absorber, while the 2002 and 2016 spectra are dominated by stronger neutral absorption that masks the warm component. All epochs display a clear Fe Kα emission complex around 6–7 keV and a pronounced Compton hump peaking at ~20–30 keV in the NuSTAR data, confirming the presence of relativistic reflection.

To model the spectra the authors employed two variants of the RELXILL family: relxillCp (broken‑power‑law emissivity) and relxilllpCp (lamppost geometry). Both models include a thermal Comptonisation continuum (cut‑off power law) and a relativistic reflection component with free iron abundance, ionisation parameter, and reflection fraction R. The lamppost version directly yields the height h of the X‑ray source above the black hole and the self‑consistent reflection fraction R_frac. Neutral absorption was modeled with tbnew, allowing for partial covering via the partcov convolution. Ionised absorption was represented by XSTAR tables (Γ = 2, free column density N_H and ionisation ξ, solar abundances except Fe and O which were free). Distant reflection was added with the borus02 model, and a soft thermal plasma (mekal) accounted for host‑galaxy emission.

Spectral fitting shows that the primary photon index remains roughly constant (Γ ≈ 2.2 ± 0.1) between 2002 and 2006, but the 2–10 keV flux drops by a factor of ~7 in 2016. This drop is driven by an increase in neutral column density (N_H ≈ (0.3–5) × 10^23 cm⁻²) and a transition to a reflection‑dominated state with R > 7. The lamppost fit yields a source height h ≤ 1.72 r_g and a reflection fraction R_frac ≈ 8–10, indicating that the X‑ray corona is extremely compact and located very close to the event horizon. In such a geometry, strong gravitational light‑bending focuses the primary emission onto the disc, suppressing the observed continuum while enhancing the reflected component.

Crucially, simultaneous fitting of all three epochs constrains the black‑hole spin to a > 0.99, implying an innermost stable circular orbit (ISCO) that lies essentially at the horizon. The high spin is robust against the complex, variable absorption because the broadband data (especially the NuSTAR coverage of the Compton hump) breaks the degeneracy between absorption curvature and relativistic broadening.

The study demonstrates that even in moderately obscured NLS1s—where partial covering neutral clouds, warm absorbers, and outflows are common—reliable spin measurements are achievable with modern relativistic reflection models and high‑energy coverage. The low coronal height and extreme light‑bending provide a natural explanation for the dramatic flux decline observed in 2016. The authors conclude that PG 1535+547 hosts a rapidly rotating supermassive black hole, and that variability in this source is dominated by changes in line‑of‑sight absorption combined with relativistic light‑bending effects. Future observations with higher spectral resolution (e.g., XRISM, Athena) will further clarify the geometry of the absorbers and the dynamics of the compact corona.