Size and kinematics of the C IV broad emission line region from microlensing-induced line profile distortions in two gravitationally lensed quasars

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.

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.

Line profile distortions are commonly observed in gravitationally lensed quasar spectra. These distortions are caused by microlensing from the stars in the lensing galaxy, which produce differential magnification of spatially and kinematically separated parts of the broad line region (BLR). The quasi-simultaneous visible and near-infrared spectroscopy of the lensed quasar Q2237+0305 reveals strong microlensing-induced line deformations in the high-ionization C IVλ1549 Å and the low-ionization Hα emission lines. We use this effect to constrain the BLR size, geometry, and kinematics in Q2237+0305. For this purpose, we modeled the deformation of the emission lines for three representative BLR models: a Keplerian disk, an equatorial wind, and a biconical polar wind. We considered various inclinations with respect to the line of sight. We find that the observed microlensing effect, characterized by a set of four indices, can only be reproduced by a subsample of the considered BLR models. The microlensing analysis favors a Keplerian disk model for the regions emitting the C IV and the Hα emission lines. A polar wind model remains possible for the C IV BLR, although it is less likely. The equatorial wind model is totally excluded. A preferred inclination of the BLR of 40° is found, in agreement with expectations for a type 1 AGN and past constraints on the accretion disk inclination. The half-light radius of the BLR is r1/2 ≃ 47 ± 19 light-days, with no significant difference between the C IV and Hα BLRs. The size of the C IV BLR agrees with the radius-luminosity relation derived from reverberation mapping, while the size of the Balmer line BLR is one order of magnitude smaller, possibly revealing different quasar properties at high luminosities and high accretion rates.

The quadruply lensed quasar HE0435−1223 shows a clear microlensing effect that affects differently the blue and red wings of the Hα line profile in its image D. To interpret these observations, and constrain the broad emission line region (BLR) properties, the effect of gravitational microlensing on quasar broad emission line profiles and their underlying continuum has been simulated considering representative BLR models and microlensing magnification maps. The amplification and distortion of the Hα line profile, characterized by a set of four indices, can be reproduced by the simulations. Although the constraints on the BLR models set by the observed single-epoch microlensing signal are not very robust, we found that flattened geometries (Keplerian disk and equatorial wind) can more easily reproduce the observed line profile deformations than a biconical polar wind. With an additional independent constraint on the size of the continuum source, the Keplerian disk model of the Hα BLR is slightly favored.

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.

Gravitational microlensing is a powerful tool allowing one to probe the structure of quasars on sub-parsec scale. We report recent results, focusing on the broad absorption and emission line regions. In particular microlensing reveals the intrinsic absorption hidden in the P Cygni-type line profiles observed in the broad absorption line quasar H1413+117, as well as the existence of an extended continuum source. In addition, polarization microlensing provides constraints on the scattering region. In the quasar Q2237+030, microlensing differently distorts the H$\alpha$ and CIV broad emission line profiles, indicating that the low- and high-ionization broad emission lines must originate from regions with distinct kinematical properties. We also present simulations of the effect of microlensing on line profiles considering simple but representative models of the broad emission line region. Comparison of observations to simulations allows us to conclude that the H$\alpha$ emitting region in Q2237+030 is best represented by a Keplerian disk.

Microlensing by stars in the lens galaxy of a gravitationally lensed quasar is a phenomenon that can selectively magnify quasar subregions, producing observable changes in the continuum brightness or distortions in the emission line profiles. Hence, microlensing allows us to probe the inner quasar regions. In this paper, we report measurements of the ratio of the broad emission line region (BLR) radius to the continuum source radius in eight lensed quasars, for the C IV, Mg II, and Hα emission lines and their respective underlying continua at λλ 1550 Å, 2800 Å, and 6563 Å. The microlensing-induced line profile distortions and continuum magnifications were observed in the same single-epoch datasets, and simultaneously compared with microlensing simulations. We found that, on average, the inner radius of the BLR starts at the end of the UV-optical continuum source, independently of the line ionization and the wavelength of the continuum. The half-light radius of the BLR is, on average, a factor of six larger than the half-light radius of the continuum source, independently of the quasar’s bolometric luminosity. We also found a correlation between the BLR radius and the continuum source radius, supporting the idea that the dominant contribution to the UV-optical continuum may come from the BLR itself. Our results independently confirm the results of reverberation mapping studies, and extend them to higher-redshift, higher-luminosity quasars.

In this paper, the sizes of the broad emission line regions (BLRs) and black hole (BH) masses of double-peaked broad low-ionization emission line emitters (DBP emitters) are compared using different methods: virial BH masses versus BH masses from stellar velocity dispersions, the size of BLRs from the continuum luminosity versus the size of BLRs from the accretion disc model. First, the virial BH masses of DBP emitters estimated by the continuum luminosity and linewidth of broad Hβ are about six times (a much larger value, if including another DBP emitters, of which the stellar velocity dispersions are traced by the linewidths of narrow emission lines) larger than the BH masses estimated from the relation MBH‐σ which is a more accurate relation to estimate BH masses. Second, the sizes of the BLRs of DBP emitters estimated by the empirical relation of RBLR‐L 5100 A are about three times (a much larger value, if including another DBP emitters, of which the stellar velocity dispersions are traced by the linewidths of narrow emission lines) larger than the mean flux-weighted sizes of BLRs of DBP emitters estimated by the accretion disc model. The higher electron density of BLRs of DBP emitters would be the main reason which leads to smaller size of BLRs than the predicted value from the continuum luminosity.

In this manuscript, the structure of broad emission line regions (BLRs) of well-mapping double-peaked emitter (AGN with broad double-peaked low-ionization emission lines) 3C390.3 is studied. Besides the best fitted results for double-peaked broad optical balmer lines of 3C390.3 by theoretical disk model, we try to find another way to further confirm the origination of double-peaked line from accretion disk. Based on the long-period observed spectra in optical band around 1995 collected from AGN WATCH project, the theoretical disk parameters of disk-like BLRs supposed by elliptical accretion disk model (Eracleous et al. 1995) have been well determined. Through the theoretical disk-like BLRs, characters of observed light-curves of broad double-peaked H$\alpha$ of 3C390.3 can be well reproduced based on the reverberation mapping technique. Thus the accretion disk model is preferred as one better model for BLRs of 3C390.3. Furthermore, we can find that different disk parameters should lead to some different results about size of BLRs of 3C390.3 from the one measured through observational data, which indicates the measured disk parameters are significantly valid for 3C390.3. After that, the precession of theoretical elliptical disk-like BLRs being considered, we can find that the expected line profile in 2000 by theoretical model is consistent with the observed line profile by HST around 2000. Based on the results, we can further believe that the origination of broad double-peaked balmer emission lines of 3C390.3 are from accretion disk around central black hole.

Here we present observational evidence that the broad emission line region (BELR) of active galactic nuclei (AGN) generally has an outer boundary. This was already clear for sources with an obvious transition between the broad and narrow components of their emission lines. We show that the narrow component of the higher-order Paschen lines is absent in all sources, revealing a broad emission line profile with a broad, flat top. This indicates that the BELR is kinematically separate from the narrow emission line region. We use the virial theorem to estimate the BELR outer radius from the flat top width of the unblended profiles of the strongest Paschen lines, Paα and Paβ, and find that it scales with the ionizing continuum luminosity roughly as expected from photoionization theory. The value of the incident continuum photon flux resulting from this relationship corresponds to that required for dust sublimation. A flat-topped broad emission line profile is produced by both a spherical gas distribution in orbital motion and an accretion disc wind if the ratio between the BELR outer and inner radius is assumed to be less than ∼100–200. On the other hand, a pure Keplerian disc can be largely excluded, since for most orientations and radial extents of the disc the emission line profile is double-horned.

Using a Gaussian fitting procedure, we have constructed an accurate UV narrow emission line template spectrum for the Seyfert 1 galaxy NGC 5548, from the low-state Hubble Space Telescope/Faint Object Spectrograph UV spectrum taken on 1982 July 5 with the 10 circular aperture. This template spectrum is similar in form to that determined for the prototypical narrow-line Seyfert 2 galaxy NGC 1068. The narrow emission line template spectrum of NGC 5548 has for the first time enabled us to isolate the narrow and broad UV emission lines and thereby to determine firm estimates for the intensities of the broad emission lines for both the low-state spectrum and the HST/FOS archival spectrum obtained from the intensive monitoring campaign undertaken by AGN Watch on 1993 April 19, when the continuum was a factor of ~4.8 brighter. A comparison of the low-state and high-state spectra show that while the narrow UV emission lines are nonvariable over timescales of ~10 months, the broad UV emission lines exhibit large variations in both their strength and shape. By combining a photoionization code with a robust global optimization routine, we have determined global best-fit parameters for the average physical conditions within the broad emission line region gas for both the low- and high-state spectra. By using a best guess estimate for the shape of the ionizing continuum of this source, we find that a single zone photoionization model of the broad emission line region (BLR) cannot simultaneously fit both the emission-line ratios of the strongest UV lines and their variability timescales. However, the line ratios and variability timescales can be reproduced if we assume a stratified BLR, i.e., a BLR that has strong gradients in density (NH ∝ 1/r2, 1011.3-1010.0). These models also suggest that the BLR gas is in a moderately high state of ionization with log10 U ~ -0.6 in the low-state spectrum, rising to ~0.0 in the high state. We find that the observed differences in the broad emission line fluxes and their ratios, between the low- and high-state spectra, are not solely a consequence of changes in the ionizing continuum source luminosity. Rather, they imply in addition, a change in the spectral energy distribution of the ionizing continuum, although changes in either the covering fraction, or composition of the broad emission line region gas, cannot necessarily be ruled out. By constructing a simple two-zone model for the high-ionization lines, we find that in order to reproduce the observed line ratios and line equivalent widths, the gas covering fraction for this source must necessarily be high, ~38% at a radial distance of 2 lt-days and decreasing outward to ~32% at 10 lt-days. This is considerably larger than the typical value of ~10% quoted for active galactic nuclei, derived from the incidence of Lyman edges in high-redshift quasars. Although the statistics for the incidence of Lyman edges in Seyfert 1 galaxies is poorly determined, our derived covering fraction is broadly consistent with the 25% covering fraction estimate obtained from observations with the Hopkins Ultraviolet Telescope.

Fast empirical models of the broad emission line region (BLR) are a powerful tool to interpret velocity-resolved reverberation mapping (RM) data, estimate the mass of the supermassive black holes, and gain insight into its geometry and kinematics. Much of the effort so far has been devoted to describing the emissivity of one emission line at a time. We present here an alternative approach aimed at describing the underlying BLR gas distribution, by exploiting simple numerical recipes to connect it with emissivity. This approach is a step toward describing multiple emission lines originating from the same gas and allows us to clarify some issues related to the interpretation of RM data. We illustrate this approach—implemented in the code CARAMEL-gas—using three data sets covering the Hβ emission line (Mrk 50, Mrk 1511, Arp 151) that have been modeled using the emissivity-based version of the code. As expected, we find differences in the parameters describing the BLR gas and emissivity distribution, but the emissivity-weighted lag measurements and all other model parameters including black hole mass and overall BLR morphology and kinematics are consistent with the previous measurements. We also model the Hα emission line for Arp 151 using both the gas- and emissivity-based BLR models. We find ionization stratification in the BLR with Hα arising at larger radii than Hβ, while all other model parameters are consistent within the uncertainties.

This paper proposes a simple, empirically derived, unifying structure for the inner regions of quasars. This structure is constructed to explain the broad absorption line regions (BALRs) and the narrow associated ultraviolet and X-ray (NALs) and is also found to explain the broad emission line regions (BELRs) and several scattering features, including a substantial fraction of the broad X-ray Fe-K emission line and the biconical extended narrow emission line region (ENLR) structures seen on large kiloparsec scales in Seyfert images. The model proposes that a funnel-shaped thin shell outflow creates all of these features. The wind arises vertically from a narrow range of radii on a disk at BELR velocities. Radiation force then accelerates the flow radially, so that it bends outward to a cone angle of ~60° and has a divergence angle of ~6°, to give a covering factor of ~10%. When the central continuum is viewed from the side, through this wind, narrow high-ionization associated ultraviolet absorption lines and the X-ray absorbers are seen, as in many low-luminosity active galactic nuclei (AGNs). When viewed end-on, the full range of velocities is seen in absorption with a large total column density, giving rise to the broad absorption line systems seen in a minority of quasars, the BAL QSOs. The wind is both warm (~106 K) and highly ionized. This warm highly ionized medium (WHIM) has a density of ~109 cm-3, putting it in pressure equilibrium with the BELR clouds; the BELR is then a cool phase embedded in the overall outflow, avoiding cloud destruction through shear. The wind has the correct ionization parameter and filling factor for this. The high- and low-ionization zones of the BELR correspond to the cylindrical and conical regions of the wind, since the former is exposed to the full continuum while the latter receives only the continuum filtered by the former. The warm wind is significantly Thomson thick along the radial flow direction, producing the polarized optical continuum found in BALs, but is only partially ionized, creating a broad fluorescent 6.4 keV Fe-K emission line and greater than 10 keV Compton hump. The conical shell outflow can produce a biconical matter-bounded NELR. Luminosity-dependent changes in the structure, reducing the cylindrical part of the flow or increasing the mean angle to the disk axis and decreasing the wind opening angle, may explain the UV and X-ray Baldwin effects and the greater prevalence of obscuration in low-luminosity AGNs.

We have obtained new spectropolarimetric observations at visible wavelengths of the changing-look active galactic nucleus (AGN) Mrk 1018. The AGN direct spectrum shows an extremely weak continuum with faint broad Hβ and Hα emission lines. Both lines can be fit with a single very broad emission line component of full width at half maximum FWHM ≃ 7200 km s−1, with no evidence of the additional 3000 km s−1-wide component that was previously detected. While this is in agreement with line formation in a Keplerian disk, the line profile variability suggests that the broad emission line region is likely more complex. The continuum polarization of Mrk 1018 is low; it is not higher in the current faint state compared to the past bright state, confirming that dust obscuration is not the mechanism at the origin of the change of look. The polarization profile of the Hα line is asymmetric with no rotation of the polarization angle, which possibly reveals line formation in a polar outflow. Alternatively, the polarization profile may be the consequence of a time delay between the direct and the polarized light. Interestingly, the polarization signatures predicted for broad lines emitted around supermassive binary black holes are not observed.

We derive a theoretical relation between RBLR, the size of the broad emission line region (BLR) of active galactic nuclei (AGNs), and the observed soft X-ray luminosity and spectrum. We show that in addition to the well-known RBLR~L1/2 scaling, RBLR should depend also on the soft X-ray spectral slope, and we derive the expected relation between RBLR and the X-ray luminosity and spectral index. Applying this relation to calculate a predicted BLR radius for 10 AGNs with reverberation data, we show that including the dependence on the spectrum improves the agreement between the calculated BLR radius and the radius independently determined from reverberation mapping. Similarly, we evaluate an expression for the line width and show that including the dependence on the spectrum significantly improves the agreement between the calculated BLR velocity dispersion and the observed FWHM of the Hβ line. The theoretical expression for the line width also provides a physical explanation to the anticorrelation between the soft X-ray slope and the emission-line width observed in narrow-line Seyfert galaxies.