Geometry and kinematics of the broad emission line region in the lensed quasar Q2237+0305

Microlensing-induced distortions of broad emission line profiles observed in the spectra of gravitationally lensed quasars can be used to probe the size, geometry, and kinematics of the broad-line region (BLR). To this end, single-epoch Mg II or Hα line profile distortions observed in five gravitationally lensed quasars, J1131-1231, J1226-0006, J1355-2257, J1339+1310, and HE0435-1223, have been compared with simulated ones. The simulations are based on three BLR models, a Keplerian disk (KD), an equatorial wind (EW), and a polar wind (PW), with different sizes, inclinations, and emissivities. The models that best reproduce the observed line profile distortions were identified using a Bayesian probabilistic approach. We find that the wide variety of observed line profile distortions can be reproduced with microlensing-induced distortions of line profiles generated by our BLR models. For J1131, J1226, and HE0435, the most likely model for the Mg II and Hα BLRs is either KD or EW, depending on the orientation of the magnification map with respect to the BLR axis. This shows that the line profile distortions depend on the position and orientation of the isovelocity parts of the BLR with respect to the caustic network, and not only on their different effective sizes. For the Mg II BLRs in J1355 and J1339, the EW model is preferred. For all objects, the PW model has a lower probability. As for the high-ionization C IV BLR, we conclude that disk geometries with kinematics dominated by either Keplerian rotation or equatorial outflow best reproduce the microlensing effects on the low-ionization Mg II and Hα emission line profiles. The half-light radii of the Mg II and Hα BLRs are measured in the range of 3 to 25 light-days. We also confirm that the size of the region emitting the low-ionization lines is larger than the region emitting the high-ionization lines, with a factor of four measured between the sizes of the Mg II and C IV emitting regions in J1339. Unexpectedly, the microlensing BLR radii of the Mg II and Hα BLRs are found to be systematically below the radius-luminosity (R − L) relations derived from reverberation mapping, confirming that the intrinsic dispersion of the BLR radii with respect to the R − L relations is large, but also revealing a selection bias that affects microlensing-based BLR size measurements. This bias arises from the fact that, if microlensing-induced line profile distortions are observed in a lensed quasar, the BLR radius should be comparable to the microlensing Einstein radius, which varies only weakly with typical lens and source redshifts.

Microlensing of the broad emission line region (BLR) in gravitationally lensed quasars produces line profile distortions that can be used to probe the BLR size, geometry, and kinematics. Based on single-epoch spectroscopic data, we analyzed the C IV line profile distortions due to microlensing in two quasars, SDSS J133907.13+131039.6 (J1339) and SDSS J113803.73+031457.7 (J1138), complementing previous studies of microlensing in the quasars Q2237+0305 and J1004+4112. J1339 shows a strong, asymmetric line profile deformation, while J1138 shows a more modest, symmetric deformation, confirming the rich diversity of microlensing-induced spectral line deformations. To probe the C IV BLR, we compared the observed line profile deformations to simulated ones. The simulations are based on three simple BLR models, a Keplerian disk (KD), an equatorial wind (EW), and a polar wind (PW), of various sizes, inclinations, and emissivities. These models were convolved with microlensing magnification maps specific to the microlensed quasar images, which produced a large number of distorted line profiles. The models that best reproduce the observed line profile deformations were then identified using a Bayesian probabilistic approach. We find that the line profile deformations can be reproduced with the simple BLR models under consideration, with no need for more complex geometries or kinematics. The models with disk geometries (KD and EW) are preferred, while the PW model is definitely less likely. In J1339, the EW model is favored, while the KD model is preferred in Q2237+0305, suggesting that various kinematical models can dominate the C IV BLR. For J1339, we find the C IV BLR half-light radii to be r1/2 = 5.1−2.9+4.6 light-days and r1/2 = 6.7−3.8+6.0 light-days from spectra obtained in 2014 and 2017, respectively. They do agree within uncertainties. For J1138, the amplitude of microlensing is smaller and more dependent on the macro-magnification factor. From spectra obtained in 2005 (single epoch), we find r1/2 = 4.9−2.7+4.9 light-days and r1/2 = 12−8+13 light-days for two extreme values of the macro-magnification factor. Combining these new measurements with those previously obtained for the quasars Q2237+0305 and J1004+4112, we show that the BLR radii estimated from microlensing do follow the C IV radius–luminosity relation obtained from reverberation mapping, although the microlensing radii seem to be systematically smaller, which could indicate either a selection bias or a real offset.

Current reverberation mapping (RM) studies primarily focus on single emission lines, particularly the Hβ line, which may not fully reveal the geometry and kinematic properties of the broad-line region (BLR). To overcome this limitation, we have conducted multiline RM observations on two highly variable active galactic nuclei (AGNs), KUG 1141+371 and UGC 3374, using the Lijiang 2.4 m telescope. Our goal was to investigate the detailed structures of different regions within the BLR. We measured the time lags of multiple broad emission lines (Hα, Hβ, Hγ, He i, and He ii) and found clear evidence of radial ionization stratification in the BLRs of both AGNs. Velocity-resolved RM analysis revealed distinct geometry and kinematics between the inner and outer regions of the BLRs. Assuming that the velocity-resolved lags reflect the kinematics of the BLR, our observations indicate that: (1) in KUG 1141+371, the inner BLR exhibits outflow signatures, while the outer region is consistent with virialized motion; and (2) in UGC 3374, the inner region displays virial motion, while the outer region shows inflow. Furthermore, we detected “breathing” behavior in the outer BLR regions of both AGNs, while the inner BLR regions show “anti-breathing,” which may be linked to intrinsic BLR properties. We discuss these findings in the context of various BLR formation models, highlighting the importance of long-term, multiline RM campaigns in understanding of BLR structure and evolution. Additionally, our results suggest that the observed stratification in BLR geometry and kinematics may contribute to the scatter in black hole mass estimates and the rapid changes in velocity-resolved RM signatures reported in recent studies.

Lensed quasars are powerful cosmic laboratories; they are used to simultaneously probe various astrophysical phenomena. Microlensing by stars within distant galaxies acts as strong gravitational lenses of multiply imaged quasars, and provides a unique and direct measurement of the lensed quasar internal structure. Microlensing of the continuum emitting region as well as the broad-line region (BLR) is well characterized by four observable indices, μcont, μBLR, WCI (wing-core), and RBI (red-blue), measured directly from the spectra. During the 2004−2007 monitoring period, image A of the quadruply lensed system Q2237+0305 underwent a strong microlensing amplification, while image D remained unaffected. We used 35 epochs of archival spectrophotometric data of Q2237+0305 obtained with the Very Large Telescope of the European Southern Observatory to develop an independent microlensing method for estimating the geometry and size of the BLR. We measured the index time series for the C IV line and the continuum emission at 1450 Å. We built a library of the simulated microlensing index time series that reproduce the observed times series based on three representative BLR models: Keplerian disk (KD), polar wind (PW), and equatorial wind (EW). After sampling the model parameter space, we find that KD is the predominant model, while PW and EW are less likely. We infer that the system is viewed at an intermediate viewing angle i ∼ 35°, and we estimate the most likely C IV BLR half-light radius r1/2 = 51 ± 23 light days. Our results are in good agreement with previous findings in the literature and extend the validity of the index-based approach to a temporal domain.

Using VLTI/GRAVITY and SINFONI data, we investigate the subparsec gas and dust structure around the nearby type 1 active galactic nucleus (AGN) hosted by NGC 3783. The K-band coverage of GRAVITY uniquely allows simultaneous analysis of the size and kinematics of the broad line region (BLR), the size and structure of the near-infrared(near-IR)-continuum-emitting hot dust, and the size of the coronal line region (CLR). We find the BLR, probed through broad Brγ emission, to be well described by a rotating, thick disc with a radial distribution of clouds peaking in the inner region. In our BLR model, the physical mean radius of 16 light-days is nearly twice the ten-day time-lag that would be measured, which closely matches the ten-day time-lag that has been measured by reverberation mapping. We measure a hot dust full-width at half-maximum (FWHM) size of 0.74 mas (0.14 pc) and further reconstruct an image of the hot dust, which reveals a faint (5% of the total flux) offset cloud that we interpret as an accreting or outflowing cloud heated by the central AGN. Finally, we directly measure the FWHM size of the nuclear CLR as traced by the [Ca VIII] and narrow Brγ line. We find a FWHM size of 2.2 mas (0.4 pc), fully in line with the expectation of the CLR located between the BLR and narrow line region. Combining all of these measurements together with larger scale near-IR integral field unit and mid-IR interferometry data, we are able to comprehensively map the structure and dynamics of gas and dust from 0.01 to 100 pc.

J1004+4112 is a lensed quasar for which the first broad emission line profile deformations due to microlensing were identified. Detailed interpretations of these features have nevertheless remained controversial. Based on 15 spectra obtained from 2003 to 2018, in this work, we revisit the microlensing effect that distorts the C IV broad emission line profile in J1004+4112. We take advantage of recent measurements of the image macro-magnification ratios, along with the fact that at one epoch, image B was not microlensed, thus constituting a reference spectrum to unambiguously characterize the microlensing effect observed in image A. After disentangling the microlensing in images A and B, we show that the microlensing-induced line profile distortions in image A, although variable, are remarkably similar over a period of 15 years. We find they are characterized by a strong magnification of the blue part of the line profile, a strong demagnification of the red part of the line profile, and a small-to-negligible demagnification of the line core. We used the microlensing effect, characterized by either the full magnification profile of the C IV emission line or a set of four integrated indices, to constrain the broad emission-line region (BLR) size, geometry, and kinematics. For this purpose, we modeled the deformation of the emission lines considering three simple, representative BLR models: a Keplerian disk, an equatorial wind, and a biconical polar wind, with various inclinations with respect to the line of sight. We find that the observed magnification profile of the C IV emission line in J1004+4112 can be reproduced with the simple BLR models we considered, without the need for more complex BLR features. The magnification appears dominated by the position of the BLR with respect to the caustic network – and not by the velocity-dependent size of the BLR. The favored models for the C IV BLR are either the Keplerian disk or the equatorial wind, depending on the orientation of the BLR axis with respect to the caustic network. We also find that the polar wind model can be discarded. We measured the C IV BLR half-light radius as r1/2=2.8−1.7+2.0 light-days. This value is smaller than the BLR radius expected from the radius-luminosity relation derived from reverberation mapping, but it is still in reasonable agreement given the large uncertainties.

We present the results of a four-year velocity-resolved reverberation mapping (RM) campaign of the changing-look active galactic nucleus (CL-AGN) NGC 4151 during its outburst phase. By measuring the time lags of the Hα, Hβ, Hγ, He i, and He ii emission lines, we confirm a stratified broad-line region (BLR) structure that aligns with predictions from photoionization models. Intriguingly, we observed an “anti-breathing” phenomenon, where the lags of broad emission lines decreased with increasing luminosity, contrary to the typical expectation. This anomaly may be attributed to the influence of the ultraviolet-optical lag or nonvirialized motions in the BLR gas. Velocity-resolved RM and ionization mapping analyses revealed rapid and significant changes in the BLR geometry and kinematics on timescales of less than a year, which cannot be interpreted by any single mechanism, such as an inhomogeneous BLR, variations in radiation pressure, or changes in the illuminated ionizing field. Additionally, the Hβ lags of NGC 4151 and other CL-AGNs agree with the radius–luminosity relationship established for AGNs with low accretion rates, implying that the CL phenomenon is more likely driven by intrinsic changes in the accretion rate rather than obscuration. These findings provide new insights into the complex internal processes of CL-AGNs and highlight the importance of long-term, multiline RM for understanding BLR structures, geometry, and kinematics.

This is the first paper in a series devoted to systematic study of the size\nand structure of the broad-line region (BLR) in active galactic nuclei (AGNs)\nusing reverberation mapping (RM) data. We employ a recently developed Bayesian\napproach that statistically describes the variabibility as a damped random walk\nprocess and delineates the BLR structure using a flexible disk geometry that\ncan account for a variety of shapes, including disks, rings, shells, and\nspheres. We allow for the possibility that the line emission may respond\nnon-linearly to the continuum, and we detrend the light curves when there is\nclear evidence for secular variation. We use a Markov Chain Monte Carlo\nimplementation based on Bayesian statistics to recover the parameters and\nuncertainties for the BLR model. The corresponding transfer function is\nobtained self-consistently. We tentatively constrain the virial factor used to\nestimate black hole masses; more accurate determinations will have to await\nvelocity-resolved RM data. Application of our method to RM data with Hbeta\nmonitoring for about 40 objects shows that the assumed BLR geometry can\nreproduce quite well the observed emission-line fluxes from the continuum light\ncurves. We find that the Hbeta BLR sizes obtained from our method are on\naverage ~20% larger than those derived from the traditional cross-correlation\nmethod. Nevertheless, we still find a tight BLR size-luminosity relation with a\nslope of alpha=0.55\\pm0.03 and an intrinsic scatter of ~0.18 dex. In\nparticular, we demonstrate that our approach yields appropriate BLR sizes for\nsome objects (such as Mrk 142 and PG 2130+099) where traditional methods\npreviously encountered difficulties.\n

The quasar 3C 273 has been observed with infrared spectroastrometry (SA) of the broad Paα line and optical reverberation mapping (RM) of the broad Hβ line. SA delivers information about the angular size and structure of the Paα broad-line region (BLR), while RM delivers information about the physical size and structure of the Hβ BLR. Based on the fact that the two BLRs share the mass of the supermassive black hole (SMBH) and viewing inclination, a combination of SA and velocity-resolved RM (SARM) thereby allows us to simultaneously determine the SMBH mass and geometric distance through dynamically modeling the two BLRs. We construct a suite of dynamical models with different geometric configurations and apply a Bayesian approach to obtain the parameter inferences. Overall the obtained masses and distances are insensitive to specific BLR configurations but more or less depend on parameterizations of the vertical distributions. The most probable model, chosen in light of the Bayes factor, yields an angular-size distance log ( D A / Mpc ) = 2.83 − 0.28 + 0.32 and an SMBH mass log ( M • / M ⊙ ) = 9.06 − 0.27 + 0.21 , which agree with the relationships between SMBH masses and bulge properties. The BLRs have an inclination of 5 − 1 + 1 degrees, consistent with that of the large-scale jet in 3C 273. Our approach reinforces the capability of SARM analysis to measure SMBH masses and distances of active galactic nuclei and quasars even though SA and RM observations are undertaken with different emission lines and/or in different periods.

Unveiling the structure of the Broad Line Region (BLR) of AGN is critical to\nunderstand the quasar phenomenon. Detail study of the geometry and kinematic of\nthese objects can answer the basic questions about the central BH mass,\naccretion mechanism and rate, growth and evolution history. Observing the\nresponse of the BLR clouds to continuum variations, Reverberation Mapping (RM)\nprovides size vs luminosity and mass vs luminosity relations for QSOs and Sy1\nAGNs with the goal to use these objects as standard candles and mass tags.\nHowever, the RM size can receive different interpretations depending on the\nassumed geometry and the corresponding mass depends on an unknown geometrical\nfactor as well on the possible confusion between local and global velocity\ndispersion. From RM alone, the scatter around the mean mass is as large as a\nfactor 3. Though BLRs are expected to be much smaller than the current spatial\nresolution of large optical interferometers (OI), we show that differential\ninterferometry with AMBER, GRAVITY and successors can measure the size and\nconstrain the geometry and kinematics on a large sample of QSOs and Sy1 AGNs.\nAMBER and GRAVITY (K around 10.5) could be easily extended up to K equal to 13\nby an external coherencer or by advanced incoherent data processing. Future\nVLTI instrument could reach K around 15. This opens a large AGN BLR program\nintended to obtain a very accurate calibration of mass, luminosity and distance\nmeasurements from RM data which will allow using many QSOs as standard candles\nand mass tags to study the general evolution of mass accretion in the Universe.\nThis program is analyzed with our BLR model allowing predicting and\ninterpreting RM and OI measures together and illustrated with the results of\nour observations of 3C273 with the VLTI.\n

ABSTRACTThe structure of the broad-line region (BLR) is an essential ingredient in the determination of active galactic nucleus (AGN) virial black hole masses, which in turn are important to study the role of black holes in galaxy evolution. Constraints on the BLR geometry and dynamics can be obtained from velocity-resolved studies using reverberation mapping data (i.e. monitoring data). However, monitoring data are observationally expensive and only available for a limited sample of AGNs, mostly confined to the local Universe. Here, we explore a new version of a Bayesian inference, physical model of the BLR that uses an individual spectrum and prior information on the BLR size from the radius–luminosity relation, to model the AGN BLR geometry and dynamics. We apply our model to a sample of 11 AGNs, which have been previously modelled using monitoring data. Our single-epoch BLR model is able to constrain some of the BLR parameters with inferred parameter values that agree within the uncertainties with those determined from the modelling of monitoring data. We find that our model is able to derive stronger constraints on the BLR for AGNs with broad emission lines that qualitatively have more substructure and more asymmetry, presumably as they contain more information to constrain the physical model. The performance of this model makes it a practical and cost-effective tool to determine some of the BLR properties of a large sample of low- and high-redshift AGNs, for which monitoring data are not available.

Recent studies have shown that line profile distortions are commonly observed in gravitationally lensed quasar spectra. We investigate the effect of gravitational microlensing on quasar broad emission line profiles and their underlying continuum, combining the emission from simple representative BLR models with generic microlensing magnification maps. Specifically, we considered Keplerian disk, polar, and equatorial wind BLR models of various sizes. The effect of microlensing has been quantified with four observables: $\mu^{BLR}$, the total magnification of the broad emission line; $\mu^{cont}$, the magnification of the underlying continuum; as well as red/blue, RBI and wings/core, WCI, indices that characterize the line profile distortions. The simulations showed that distortions of line profiles, such as those recently observed in lensed quasars, can indeed be reproduced and attributed to the differential effect of microlensing on spatially separated regions of the BLR. While the magnification of the emission line $\mu^{BLR}$ sets an upper limit on the BLR size and, similarly, the magnification of the continuum $\mu^{cont}$ sets an upper limit on the size of the continuum source, the line profile distortions mainly depend on the BLR geometry and kinematics. We thus built (WCI,RBI) diagrams that can serve as diagnostic diagrams to discriminate between the various BLR models on the basis of quantitative measurements. It appears that a strong microlensing effect puts important constraints on the size of the BLR and on its distance to the high-magnification caustic. In that case, BLR models with different geometries and kinematics are more prone to produce distinctive line profile distortions for a limited number of caustic configurations, which facilitates their discrimination.

It has long remained an open question as to the geometry of the broad line region (BLR) in active galactic nuclei (AGNs). The reverberation mapping technique which measures the response of the broad emission lines to the ionizing continuum, when combined with multiwavelength continuum fitted by sophisticated accretion disks, provides a way of probing the BLR geometry. We analyze a sample of 35 AGNs, which have been monitored by the reverberation mapping campaign. In view of energy budget, the reverberation-based BH masses are found to be in agreement with those obtained by accretion disk models in two thirds of the present sample while the reverberation mapping methods underestimate the BH masses in about one third of objects, as also suggested by Collin et al. in a recent work. We point out that there are obviously two kinds of BLR geometry, which are strongly dependent on the Eddington ratio, and separated by the value LBol/LEdd∼0.1. These results prefer a scenario of the disk and wind configuration of the BLR and identify the Eddington ratio as the physical driver regulating the wind in the BLR.

Understanding the geometry and kinematics of the broad line region (BLR) of active galactic nuclei (AGN) is important to estimate black hole masses in AGN and study the accretion process. The technique of reverberation mapping (RM) has provided estimates of BLR size for more than 100 AGN now; however, the structure of the BLR has been studied for only a handful number of objects. Towards this, we investigated the geometry of the BLR for a large sample of 57 AGN using archival RM data. We performed systematic modelling of the continuum and emission line light curves using a Markov chain Monte Carlo method based on Bayesian statistics implemented in PBMAP (Parallel Bayesian code for reverberation−MAPping data) code to constrain BLR geometrical parameters and recover velocity integrated transfer function. We found that the recovered transfer functions have various shapes such as single-peaked, double-peaked, and top-hat suggesting that AGN have very different BLR geometries. Our model lags are in general consistent with that estimated using the conventional cross-correlation methods. The BLR sizes obtained from our modelling approach is related to the luminosity with a slope of 0.583 ± 0.026 and 0.471 ± 0.084 based on H β and H α lines, respectively. We found a non-linear response of emission line fluxes to the ionizing optical continuum for 93 ${{\ \rm per\ cent}}$ objects. The estimated virial factors for the AGN studied in this work range from 0.79 to 4.94 having a mean at 1.78 ± 1.77 consistent with the values found in the literature.

Photometric reverberation mapping employs a wide band pass to measure the AGN continuum variations and a suitable narrow band to trace the echo of an emission line in the broad line region (BLR). The narrow band catches both the emission line and the underlying continuum, and one needs to extract the pure emission line light curve. We performed a test on two local AGNs, PG0003+199 and Ark120, by observing well-sampled broad- (B, V) and narrow-band light curves with the robotic 15 cm telescope VYSOS-6 on Cerro Armazones, Chile. We find that, as long as the emission line contributes 50% to the band pass, the pure emission line light curve can be reconstructed from photometric monitoring data so that the time lag τ can be measured. For both objects the lags are consistent with spectroscopic reverberation results. We calculated virial black hole masses in agreement with literature values, by combining the BLR size RBLR (τ) from photometric monitoring with the velocity dispersion of a single contemporaneous spectrum. Applying the flux variation gradient method, we estimate the host galaxy contribution in the apertures used and the host-subtracted restframe 5100 Å luminosity LAGN. Our LAGN differs significantly from previous estimates, placing both sources ~50% closer to the RBLR − LAGN relation. This suggests that the scatter in the current RBLR − LAGN relation is largely caused by uncertainties in RBLR due to undersampled light curves and by uncertainties in the host-subtracted AGN luminosities inferred so far. If the scatter can be reduced, then two quasar samples matching in RBLR should also match in intrinsic LAGN, independent of redshift, thus offering the prospect of probing cosmological models. Photometric reverberation mapping opens the door to efficiently measuring hundreds of BLR sizes and host-subtracted AGN luminosities even with small telescopes, but also routinely with upcoming large survey telescopes like the LSST.