Constraining the Physical Properties of Quadruply Lensed Quasars using Optical-to-X-Ray Data

Technical Summary: Constraining the Physical Properties of Quadruply Lensed Quasars using Optical-to-X-Ray Data

Problem Statement
Determining the bolometric luminosity ($L_{\text{bol}}$) of Active Galactic Nuclei (AGN) is a fundamental challenge in astrophysics, as it is required to calculate accretion rates and Eddington ratios ($\lambda_{\text{Edd}}$), which are crucial for understanding black hole growth. While $L_{\text{bol}}$ is best derived from broadband Spectral Energy Distributions (SEDs) spanning optical, ultraviolet (UV), and X-ray bands, such multiwavelength data are often unavailable. This limitation is particularly acute for gravitationally lensed AGN, where existing literature consists largely of source-specific studies or small samples with limited data. Furthermore, while X-ray observations offer advantages for lensed systems (being less contaminated by host galaxy light and foreground lenses), they are susceptible to microlensing due to the small physical scale of the X-ray emitting region. Conversely, optical/UV data are prone to host galaxy contamination and dust extinction but offer better spatial resolution. There is a lack of a homogeneous, systematic study constraining the physical properties of a large sample of lensed AGN using consistent multiwavelength data.

Methodology
The authors compiled a sample of 27 quadruply lensed AGN systems (redshift $0.6 < z < 3.1$) from the JWST GO-2046 proposal and archival data. The study employed a multi-step analysis:

1. Data Compilation: The team gathered optical/UV photometry and spectroscopy, as well as Chandra X-ray observations. For 17 systems, they constructed broadband SEDs by combining optical/UV photometry (from literature or new HST/WFC3 observations) with X-ray spectra. For 9 systems lacking sufficient optical/UV bands, they relied on X-ray data. Two systems were excluded due to insufficient data.
2. SED Fitting: For the 17 systems with multi-band data, the authors fitted the optical/UV SEDs using the DISKPN model (a thermal accretion disk model) in XSPEC, accounting for dust extinction. X-ray spectra (0.3–10 keV) were fitted with a power-law model (POWERLAW) including Galactic absorption (TBABS). The total $L_{\text{bol}}$ was calculated by summing the intrinsic optical/UV disk luminosity ($L_{\text{UV}}$) and the intrinsic X-ray luminosity ($L_X$).
3. Bolometric Corrections: For the 9 systems without full SEDs, $L_{\text{bol}}$ was estimated using luminosity-dependent bolometric corrections. Specifically, the 2–10 keV X-ray bolometric correction ($\kappa_{2-10}$) was applied, utilizing a relation derived from non-lensed AGN. One system was estimated using an optical bolometric correction at 3000 Å.
4. Black Hole Mass Estimation: For 19 systems with available rest-frame optical/UV spectra, black hole masses ($M_{\text{BH}}$) were estimated using the virial method. This involved fitting broad emission lines (Mg II, C IV, or H$\alpha$) to determine the Full Width at Half Maximum (FWHM) and applying empirical relations between line width, continuum luminosity, and BLR radius.
5. Microlensing Diagnostics: The authors utilized the well-known non-linear correlation between monochromatic UV ($L_{2500\text{\AA}}$) and X-ray ($L_{2\text{keV}}$) luminosities to test for systematic uncertainties in lens magnification corrections, specifically looking for deviations caused by macro-model biases or micro/milli-lensing.

Key Contributions and Results

• Sample Size: This work presents one of the largest homogeneous multiwavelength studies of quadruply lensed AGN to date, analyzing 27 systems. For most of these systems, measurements of $L_{\text{bol}}$, $M_{\text{BH}}$, and $\lambda_{\text{Edd}}$ are reported for the first time.
• Bolometric Luminosity Estimation: The study validates the use of luminosity-dependent 2–10 keV X-ray bolometric corrections ($\kappa_{2-10}$) for lensed AGN. The authors find that these corrections provide unbiased and reliable predictions of $L_{\text{bol}}$ with a maximum scatter of $\sim 0.5$ dex. In contrast, optical bolometric corrections tend to overestimate $L_{\text{bol}}$ but show lower scatter once microlensing effects are accounted for.
• Microlensing Detection: By plotting the lensed images on the UV-X-ray luminosity relation, the authors identified outliers. They suggest that deviations in systems like DESJ0405 and J1042 are likely caused by differential microlensing (magnification in X-rays vs. demagnification in UV, or vice versa), rather than macro-model errors. This demonstrates the utility of the UV-X-ray relation as a diagnostic tool for identifying micro/milli-lensing events in lensed quasars.
• Physical Properties: The study provides a distribution of physical properties for the sample, showing a range of bolometric luminosities ($10^{44}$ to $10^{48}$ erg s$^{-1}$), black hole masses ($10^7$ to $10^{10} M_\odot$), and Eddington ratios.

Significance
The paper claims to address a critical gap in the literature by providing a systematic, multiwavelength characterization of lensed AGN. By demonstrating that luminosity-dependent X-ray bolometric corrections are reliable estimators for $L_{\text{bol}}$ in the absence of full SEDs, the authors provide a robust method for studying lensed AGN where optical data is incomplete. Furthermore, the work establishes the UV-X-ray luminosity relation as a novel and effective tool for quantifying uncertainties in magnification corrections induced by microlensing, thereby improving the accuracy of physical property measurements for lensed quasars. The results support the use of lensed AGN as "natural telescopes" to probe inner AGN structures, provided that microlensing effects are properly constrained.