Size and kinematics of the low-ionization broad emission line region from microlensing-induced line profile distortions in gravitationally lensed quasars

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.

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.

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.

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.

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.

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.

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.

We analyze the properties of the broad line region (BLR) in low luminosity\nAGN by using HST/STIS spectra. We consider a sample of 24 nearby galaxies in\nwhich the presence of a BLR has been reported from their Palomar ground-based\nspectra. Following a widely used strategy, we used the [SII] doublet to\nsubtract the contribution of the narrow emission lines to the H-alpha+[NII]\ncomplex and to isolate the BLR emission. Significant residuals that suggest a\nBLR, are present. However, the results change substantially when the [OI]\ndoublet is used. Furthermore, the spectra are also reproduced well by just\nincluding a wing in the narrow H-alpha and [NII] lines, thus not requiring the\npresence of a BLR. We conclude that complex structure of the narrow line region\n(NLR) is not captured with this approach and that it does not lead to general\nrobust constraints on the properties of the BLR in these low luminosity AGN.\nNonetheless, the existence of a BLR is firmly established in 5 Seyferts, and 5\nLINERs. However, the measured BLR fluxes and widths in the 5 LINERs differ\nsubstantially with respect to the ground-based data. The BLR sizes in LINERs,\nwhich are estimated by using the virial formula from the line widths and the\nblack hole mass, are about 1 order of magnitude greater than the extrapolation\nto low luminosities of the relation between the BLR radius and AGN luminosity\nobserved in more powerful active nuclei. We ascribe the larger BLR radius to\nthe lower accretion rate in LINERs when compared to the Seyfert, which causes\nthe formation of an inner region dominated by an advection-dominated accretion\nflow (ADAF). The estimated BLR sizes in LINERs are comparable to the radius\nwhere the transition between the ADAF and the standard thin disk occurs due to\ndisk evaporation.\n

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.

Accretion onto a blackhole can indue extreme radiation stimulating emission lines in nearby gases. In the deep gravitational well, these clouds move at speeds up to 10,000 km/s, and so broadens the emission lines via the relativistic Doppler Effect. These broad emission lines on spectra are the most prominent feature of AGNs. The region around central black hole, where such broad emission lines are produced, is known as the broad-line region (BLR). We developed a self-consistent model to describe both the geometry of the BLR and the dynamics of the gas clouds in the BLR. The motion of clouds is described by a series of Keplerian orbits. The spectral shift of emission line can be derived by the radial velocity at each point in Keplerian orbits. Assuming that the BLR is stable, that is, the structure and distribution of gas in BLR doesn’t change over a considerable period of time, the accumulated emission line profile can be constructed by simple addition of all Keplerian orbits in the BLR. By comparing the calculated profiles to spectra observed in SDSS (Sloan Digital Sky Survey), we found that our model can satisfactorily match most observed profiles. We attempted to review and summarize potential binary AGNs. We applied the model to the study of binary AGNs and their profiles of emission lines. We considered the velocity offset due to mutual Keplerian motion of black holes in binary and added their respective profiles to produce the final, resulting emission line profile from binary AGNs. Three kinds of emission line profiles are featured in our simulations: twin-peak structures where two strong broad-lines are present; sub-peak structures where a smaller but still visible emission line accompanies the dominant line; single broad emission peak but with significant velocity offset from the systematic redshift that is derived from the narrow emission-lines on the same spectrum. We then analyzed 1,348 AGN spectra with high S/N ratio in SDSS database, and selected 26 as candidate binary supermassive blackhole systems. These candidates are valuable for following-up observation, their binary nature could be confirmed by detecting the profile variation due to orbital motions.

Spectroscopic observations of 13 Seyfert 1 galaxies made from 1979 to 1984 at Palomar and Steward Observatories were analyzed for Hβ line profile variations. Significant profile changes were detected in five galaxies. These variations are often associated with changes in continuum and/or Hβ strength, suggesting a causal connection between variations in the continuum source and changes in emission-line profiles. Moreover, the overall strength of Hβ line flux variations is strongly correlated with line width in the far wings, whereas no such correlation was found in the line core. This finding is consistent with previous observations (Shuder 1982) which suggest that higher velocity gas lies closer to the central continuum source, since this inner gas would be expected to track variations more closely. In each of the five galaxies in which Hβ line profile variations were detected, particularly NGC 5548 and Mrk 6, only two distinct line profile shapes were observed. Since intermediate profile types were not detected, the transition from one type to the other must occur rapidly (within a few months). These rapid line profile variations are thus consistent with numerous other studies which strongly point to a small broad-line region (BLR) size in Seyfert galaxies. Observational models of BLR kinematics and structure for each of the galaxies in the present sample are presented. Line asymmetries suggest that radial motions exist in the BLRs of most of these objects. Chaotic motions also appear to be quite common. However, evidence for rotational motion was observed in only one case. A model that is consistent with the majority of this sample incorporates both radial and chaotic motions combined with some form of external obscuration such as an optically and geometrically thick torus.

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.

We discuss the effects of microlensing on the broad emission lines (BELs) of QSOs in the light of recent determinations of the size of the broad-line region (BLR) and its scaling with QSO luminosity. Microlensing by star-sized objects can produce significant amplifications in the BEL of some multiple-imaged QSOs, and could be very relevant for high-ionization lines. We have identified a group of 10 gravitational lens systems (~30% of the selected sample) in which microlensing could be observed. Using standard kinematic models for active galactic nuclei, we have studied the changes induced in the line profile by a microlens located at different positions with respect to the center of the BLR. We found that microlensing could produce important effects such as the relative enhancement of different parts of the line profile or the displacement of the peak of the line. The study of BEL profiles of different ionization in a microlensed QSO image could be an alternative method for probing the BLR structure and size.

Context.We analyze the broad Hβline profile variability of a “changing look” active galactic nucleus (CL-AGN) NGC 3516 over an extensive period of 25 years (from 1996 to 2021). The observed change in the broad line profile may indicate a change in the geometry of the broad line region (BLR). The main objective is to follow and understand the change in the BLR over a long period as well as its connection to the CL mechanism.Aims.Using spectral line profiles, we aim to explore changes in the kinematics and dimensions of the BLR in NGC 3516. We consider two possible scenarios: the changes in the broad-line emission are either caused by a decrease of ionisation continuum emission or by the BLR obscuration by outer dusty regions. With this investigation, we aim to clarify the CL mechanism of this AGN.Methods.We analyzed the spectral band around the Hβline as well as the broad Hβline parameters and how they change over time. We modelled the broad-line profiles, assuming that there is an emission from the accretion disc superposed with emission from a surrounding region that is outside the disc.Results.We find that in the type 1 activity phase occurring when the strong broad emission lines are observed, the BLR is very complex. There is a clear disc-like BLR that contributes to the broad line wings and an additional intermediate line region (ILR) that contributes to the line core. In the high-activity phase, the ILR emission is close to the center of the line, although in some cases, it is slightly shifted to the red. In the low-activity phase (i.e. type 2 phase), the ILR component has a significant shift to the blue, indicating an outflow.Conclusions.We propose that the changing-look mechanism in NGC 3516 is rather connected with the intrinsic effects than with an outer obscuring region. It may still be possible that the dust plays an important role in the low-activity phase when it is coming from within the BLR, leading to a dusty BLR. In this way, it would cause a decrease in the ionisation and recombination rates.

I study the location of the γ-ray emission in blazar jets by creating a Compton-scattering approximation that is valid for all anisotropic radiation fields in the Thomson through Klein–Nishina regimes, is highly accurate, and can speed up numerical calculations by up to a factor of ∼10. I apply this approximation to synchrotron self-Compton, external Compton scattering of photons from the accretion disk, broad line region (BLR), and dust torus. I use a stratified BLR model and include detailed Compton-scattering calculations of a spherical and flattened BLR. I create two dust torus models, one where the torus is an annulus and one where it is an extended disk. I present detailed calculations of the photoabsorption optical depth using my detailed BLR and dust torus models, including the full angle dependence. I apply these calculations to the emission from a relativistically moving blob traveling through these radiation fields. The ratio of γ-ray to optical flux produces a predictable pattern that could help locate the γ-ray emission region. I show that the bright flare from 3C 454.3 in 2010 November detected by the Fermi Large Area Telescope is unlikely to originate from a single blob inside the BLR. This is because it moves outside the BLR in a time shorter than the flare duration, although emission by multiple blobs inside the BLR is possible. Also, γ-rays are unlikely to originate from outside of the BLR, due to the scattering of photons from an extended dust torus, since the cooling timescale would be too long to explain the observed short variability.