Open Access
Issue
A&A
Volume 677, September 2023
Article Number A94
Number of page(s) 11
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202346766
Published online 11 September 2023

© The Authors 2023

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

Accretion onto supermassive black holes (BHs) is tightly linked to their evolution over cosmic time. During the quasar phase, inflowing material forms a radiatively efficient accretion disk (AD), whose nature is still poorly understood. As accretion phenomena in quasars cannot be spatially resolved, indirect methods such as reverberation mapping (RM; Blandford & McKee 1982; Peterson 1993) and gravitational microlensing (Wambsganss 2006) are required to probe them.

Although thin disk models (hereafter SS73 disks) are an extremely useful prescription for modeling ADs and for estimating some of the fundamental properties of growing BHs, it has long been known that they, at least in their simplest form (Shakura & Sunyaev 1973), are inconsistent with quasar observations. Discrepancies between prevalent AD theory and observations are clearly manifested when considering the effective size of the continuum-emitting region in the UV-optical band, which is believed to be dominated by AD emission. While standard theory predicts that the effective wavelength-dependent size of the continuum-emitting region, r(λ), increases with wavelength as λ4/3, microlensing studies, which currently provide the only means to spatially resolve ADs in luminous high-redshift sources, often find a much flatter dependence across the UV-optical range (Blackburne et al. 2011; Jiménez-Vicente et al. 2014). Furthermore, the implied AD sizes are much larger than predicted by theory, especially at rest-UV energies (Fian et al. 2016, 2018, 2021; Motta et al. 2017; Cornachione et al. 2020a,b; Cornachione & Morgan 2020; Rojas et al. 2020). This fact appears to be reinforced by RM studies of low-luminosity low-redshift sources, which employ a different set of assumptions and are prone to different biases (Fausnaugh et al. 2018; Edelson et al. 2019; Chelouche et al. 2019; Fian et al. 2022, 2023). It therefore appears to be a fairly robust conclusion that quasar continuum emission regions are larger than predicted by standard AD theory, and that their implied physics (e.g., the AD temperature profile) is markedly different from that expected from basic energetic arguments. This tension between quasar observations and theory has far-reaching implications for accretion physics in the general astrophysical context (Dexter & Agol 2011; Gardner & Done 2017; Gaskell & Harrington 2018; Hall et al. 2018; Li et al. 2019) and in particular for explaining the census of supermassive BHs; considerable theoretical efforts are being made to remedy this (Dexter & Agol 2011; Hall et al. 2018; Gaskell & Harrington 2018; Li et al. 2019).

Although it is commonly assumed that the UV-optical phenomenology of quasars is a good probe of AD physics, it has been theoretically argued (for more than 20 yr) that the contribution of the broad-line region (BLR) to the continuum signal, albeit poorly constrained by observations, should be non-negligible, especially in low-luminosity sources (Korista & Goad 2001, 2019; Netzer 2022). Evidence has recently emerged for the non-negligible contribution of diffuse non-disk continuum emission to the time-varying signal of some low-luminosity, low-z sources (Cackett et al. 2018; Chelouche et al. 2019), which is consistent with a BLR origin (Korista & Goad 2019; Chelouche et al. 2019; Netzer 2020). Nevertheless, it is as yet unclear whether the findings hold for all low-luminosity sources, or even for the quasar population as a whole, in which case the standard AD paradigm may be correct but our interpretation of the data would require a major revision (Chelouche et al. 2019). Alternatively, the findings may indicate that our understanding of both the AD and the BLR in quasars is highly incomplete, or even fundamentally flawed.

Motivated by this, we aim to explore the effect of diffuse BLR emission on size inferences from microlensing, which, unlike RM, is sensitive to the intrinsically varying and non-varying components of the nuclear signal. Doing so would test current models for the AD and the BLR at a substantially different part of the quasar phase space, for many more objects than are currently available through high-fidelity RM campaigns (Kara et al. 2021, 2023), and using a method whose systematic effects are very different from those inflicting RM. To this end, we employ state-of-the-art emission models for quasars (neglecting host emission) and high-resolution magnification maps, and compare attributes of their magnification-distribution functions to those obtained for pure SS73 disk models.

The paper is organized as follows. In Sect. 2 we describe the microlensing model and discuss its implementation in this work. Section 3 outlines a qualitative approach for modeling nuclear and diffuse quasar emission, and shows the emergence of biases in continuum-region size estimations due to the presence of the BLR. In Sect. 4 a physical model for the BLR is presented, which is calibrated by spectroscopic and RM data, and its effect on the measured microlensing sizes is quantified and compared to published results for the lensed quasar sample of Blackburne et al. (2011). Finally, in Sect. 5 we conclude by summarizing our key findings and discussing their broader implications.

2. The microlensing model

Assessing the effect of diffuse BLR emission on microlensing signals involves several length scales defining the problem: the AD scale and the BLR scale (see Sects. 3 and 4), as well as the relevant scales in the microlensing map associated with the Einstein radius, rE, of the lensing stars and their field density (convergence, or optical depth for lensing), which we now outline.

We simulate the effect of microlenses on the total macro-image flux using magnification maps calculated by the inverse polygon mapping (IPM) method (Mediavilla et al. 2006, 2011). The IPM method is based on the properties of polygon transformation under the inverse mapping defined by the lens equation. The IPM achieves a remarkably high level of accuracy (relative error ∼5 × 10−4), accomplishing this in a significantly shorter computing time (less than 4%) compared to the inverse ray shooting (IRS) method (Wambsganss 1999). The computational efficiency improvement of the IPM over the IRS can exceed two orders of magnitude when both methods are required to deliver the same level of accuracy.

The general characteristics of a magnification map are determined (for each quasar image) by the local convergence, κ, and the local shear, γ, which can be obtained by fitting a model to the coordinates of the macro images. The local convergence is proportional to the total surface mass density κ = κc + κ, where κc is the convergence caused by continuously distributed matter (e.g., dark matter) and κ can be attributed to stellar-mass point lenses. High convergence signifies a large mass concentration and, consequently, significant magnification. Conversely, lower convergence signifies a lower mass concentration, leading to reduced light bending and suppressed magnification. As such, convergence is pivotal in establishing the amplitude of the magnification map. The shear affects the shape of the magnification map, causing a distortion effect that manifests as a stretching along one axis, coupled with a compression along its perpendicular counterpart. Greater shear values result in highly extended features in the magnification map, whereas lower shear values induce smaller distortions, leading to features in the map that maintain a closer resemblance to their original forms.

We used the latest estimate for the fraction of the surface mass density in the form of stars (∼20% near the Einstein radius of the lenses; see Jiménez-Vicente et al. 2015). The maps are 20 000 × 20 000 pixels2 large with sides of length 100 rE (rE ∼ 3 × 1016 cm; cf. Blackburne et al. 2011) and a pixel size of 0.05 light-days, which is smaller than the typical size of the UV AD (see Fig. 1). We then randomly distributed microlenses of mass M = 1 M across the magnification patterns to create the convergence in stars κ. Although the details of the microlensing pattern change with the shape of the initial mass function of the lens galaxy, the magnification probability distributions are insensitive to the choice of the initial mass function (Lewis & Irwin 1995; Wyithe & Turner 2001; Chan et al. 2021). The source sizes can be scaled to a different stellar mass by noting that r E M 1/2 $ r_{\rm E} \propto M_\star^{1/2} $. The value at each point in the map is equal to the magnification of the background source at that point, relative to the average macro-image magnification.

thumbnail Fig. 1.

Billion-pixel magnification map showing a complex caustic pattern for the positive-parity image, with brighter colors indicating larger magnifications. The insets show successive blow-ups of the maps with relevant quasar emission regions overlaid.

Characterizing the statistics of microlensing magnifications, we examined two typical cases: a positive-parity strongly lensed quasar image (minimum of the time-delay function) with κ = γ = 0.4, and a theoretical average magnification of μ = 5; and a negative-parity quasar image (saddle point) with κ = γ = 0.6, and μ = −5. The positive-parity simulation includes ∼14 300 stellar lenses, and the negative-parity simulation includes ∼21 500 stellar lenses. Microlensing-induced magnifications of our fiducial quasar are obtained by convolving its emission region (see below) with the resulting magnification map (see Fig. 1). In cases where the size of the emission is comparable to the size of the microlensing map and the derived statistic is questionable (this occurs when considering, for example, BLR emission from the most luminous quasars), calculations are reported using a lower resolution map with a pixel size of 0.25 light-days and a size of 400 rE (not shown) while holding the other parameters fixed at their aforementioned values.

Below we quantify the bias in the estimation of the AD size when diffuse (non-AD) emission is present but is not accounted for by the models that are employed to infer the disk size from the data. To this end we convolve models for quasar emission that include the AD and the BLR signals (henceforth SS73+BLR models; see below) with the computed magnification maps, and quantify the variance of the magnification-distribution function (Fig. 2), VarSS73 + BLR. This is, essentially, a measure of the observed microlensing-induced variability for a source that fully samples the microlensing map. While this may not be strictly true for particular sources observed over short times1, it is expected to be a more meaningful measure when compiling the results for a sample of lensed quasars, which better sample the landscape of their microlensing maps. This approach is also motivated by the Blackburne et al. (2011) sample used here, for which only single epoch measurements are available (see Sect. 4). Ignoring prior knowledge of the model used to simulate the data, we next seek a pure SS73-disk model of half-light radius r1/2, whose convolution with the same magnification map yields a variance VarSS73(r1/2) that minimizes ζ(r1/2) = |VarSS73(r1/2)−VarSS73 + BLR|/VarSS73 + BLR, thereby yielding the recovered r1/2. This may then be contrasted with the values for the effective half-light radius and the AD scale used for the AD+BLR model (see below).

thumbnail Fig. 2.

Magnification distributions (relative to the mean magnification of the image) for several assumptions regarding the contribution of a diffuse emission component to the flux in some band, which is parameterized by α (see the main text). Bold distribution curves mark α = 0.0 and α = 1.0 cases for each image parity (see the legend), with the thin curves separated by Δα = 0.2 intervals. The top panels show the emissivity profiles in log-units for three models, with their α-values denoted. Color-coding is arbitrarily but identically normalized in all panels. Note the reduced second moment of the magnification distribution function with increasing α.

We note that in all cases explored here, the approach resulted in a well-defined minimum for ζ, and was able to accurately recover the input value when a pure SS73 disk was used to model the source. We verified that our approach is also consistent with a maximum likelihood scheme that aims to find an optimal match between calculated magnification distribution functions, where a sizeable (±1 mag) uniform prior is used on the mean magnitude, and an extremum is sought (Fian et al. 2021). Our approach is further warranted by the current data quality for lensed quasars; higher moments of the magnification distribution, which are not explored here, may be accessible to future surveys.

3. A qualitative emission model for quasars

To qualitatively assess the effect of diffuse non-AD emission on the recovered AD properties, we first consider a simplified approach wherein the normalized surface brightness (SB) over a narrow range of wavelength, Δλ/λ ≪ 1, satisfies

I ( R ) = 1 α I 1 [ exp ( R 3 / 4 ) 1 ] 1 + α I 2 R β χ [ r in , r out ] , $$ \begin{aligned} I(R)=\frac{1-\alpha }{I_1}\left[{\mathrm{exp}(R^{3/4})-1}\right]^{-1} + \frac{\alpha }{I_2} R^\beta \chi _{[r_{\rm in},r_{\rm out}]}, \end{aligned} $$(1)

where the normalized radial coordinate, R = r/rs, and rs is the wavelength-dependent disk scale. The first/second term in Eq. (1) describes the AD/BLR contribution to the signal, whose weight is parameterized by α ∈ [0, 1], which may be wavelength-dependent (see below). Unless otherwise stated, we neglect the inner and outer disk boundaries. For an SS73 disk, the half-light radius of the disk, r1/2 ≃ 2.44 rs (cf. r1/2 ≃ 1.18 rs for a Gaussian disk), is set by the monochromatic rest optical luminosity, Lopt, of the source up to an unknown inclination angle effect (a face on configuration is assumed throughout this work), and so long as the emission is radiated by the self-similar parts of the disk (see Davis & Laor 2011; Laor & Davis 2011). In this case,

r 1 / 2 1.5 ( L opt 10 45 erg s 1 ) 1 / 2 ( λ 5100 Å ) 4 / 3 light days . $$ \begin{aligned} r_{1/2}\simeq 1.5 \left(\frac{L_{\rm opt}}{10^{45}\,\mathrm{erg\,s^{-1}}}\right)^{1/2} \left(\frac{\lambda }{5100\,\AA }\right)^{4/3}\,\mathrm{light\,days}. \end{aligned} $$(2)

Therefore, for a fiducial quasar with an optical luminosity, Lopt, 45 = Lopt/1045 erg s−1 = 1 (Morgan et al. 2010; Blackburne et al. 2011), rs = rs,3000(λ/3000 Å)/3 with rs, 3000 ≃ 0.3 light-days. The normalization factors in Eq. (1), I 1 0 [exp( R 3/4 )1] 1 RdR $ I_1 \equiv \int_0^\infty [{{\rm exp}(R^{3/4})-1}]^{-1} R\,{\rm d}R $, I 2 0 R β+1 χ [ r in , r out ] (R)dR $ I_2 \equiv \int_0^\infty R^{\beta+1} \chi_{[r_{\rm in},r_{\rm out}]}(R)\,{\rm d}R $, where the step function χ[rin, rout](R) = 1 if r ∈ [rin, rout] and is zero otherwise. The parameter β sets the radial dependence of the diffuse SB component, and is taken here to be in the range [ − 3, 0], but is likely to be closer to β ∼ −2 for most continuum emission components of interest (see Korista & Goad 2019). Specifically, the SB depends on the combined effect of the local emissivity of the gas (as follows from its illumination by the central source), its geometrical configuration (e.g., the covering factor over the quasar sky), and radiative transfer effects along the sightline to the observer (absorption and scattering). In our formalism a three-dimensional emissivity profile, ϵ(r)∝rγ with a radial geometrical covering factor of C(r)∝rδ would translate to a radial dependence of the SB, which is ∝rβ with β = γ + δ + 1. Therefore, a typical emissivity profile with γ = −2 (Korista & Goad 2019), and δ = −1 (Baskin et al. 2014), translates to β = −2.

For the diffuse BLR component, we assumed an inner radius rin = 70 rs, 3000 and an outer radius rout = 2 × 103rs, 3000, which translates to a half-light radius of the BLR emission of ∼20 light-days for β = −2 and a half-light radius of ∼90 light-days for β = 0. The latter values may be more appropriate for low-ionization emission lines (see Korista & Goad 2019) and in particular the Balmer lines, whose emission may be collisionally suppressed at small radii (Baskin et al. 2014). Therefore, the model is qualitatively consistent with the BLR size-luminosity relation (Bentz et al. 2013).

As expected, we find that the magnification distributions narrow down with rising α (Fig. 2). The details of the magnification distributions are different for the positive and negative parity images of the quasar, and yet they lead to similar conclusions. Hence, only the negative parity results will be discussed below. Not unexpectedly, larger emitting regions average over more extended parts of the magnification map2, which suppress significant (de-)magnification events. Qualitatively, a best-match criterion between the magnification-distribution function that characterizes AD-BLR models and those of pure-disk models that are (erroneously) used to recover a disk size, would lead to overestimations. More quantitatively, the implied disk scale, rs, in the presence of diffuse BLR emission can be several times larger than the true AD scale of the system for our fiducial quasar with an optical luminosity, Lopt = 1045 erg s−1, even for modest values of α ∼ 0.2 (see Fig. 3), which are motivated by photoionization calculations (see Korista & Goad 2019 and the inset of our Fig. 4). Interestingly, non-negligible biases also emerge when comparing the calculated half-light radius for the AD-BLR model to that implied by the best-matched pure-disk one, with overestimations factors of about two for typical values of α. This result contrasts previous works, which often conjecture that the emissivity profile has little effect on half-light size estimations (see Mortonson et al. 2005).

thumbnail Fig. 3.

Simulations quantifying the overestimation factors of the disk scale and the half-light radius as a function of the relative contribution of a diffuse BLR emission component to the total continuum signal, α. The disk-scale overestimation factor is defined as the ratio of the value recovered via the magnification-distribution function that matches the input value, Rs, in our AD-BLR model. Results for the half-light radius are defined by the ratio of the value recovered for the best-matched pure AD model to that which characterizes the AD-BLR model.

The bias level in source-size estimations depends on the interplay between the size of quasar constituents and the Einstein radius for a fixed α. To quantify this effect, we consider high-luminosity (Lopt = 1046 erg s−1) and low-luminosity (Lopt = 1044 erg s−1) quasars. We assume that the AD and BLR sizes scale as L opt 1 / 2 $ L_{\mathrm{opt}}^{1/2} $ (see Bentz et al. 2013), and show the results in Fig. 3. An overestimation of the disk size by factors of ≲10 is obtained for low-luminosity sources. Size inferences for high-luminosity sources are less biased and are overestimated by ∼50% for otherwise similar model parameterizations. Similar behavior with source luminosity is seen for the half-light size estimations’ bias.

More generally, the source size and the relative contribution of the diffuse BLR component to the signal are wavelength-dependent for a given quasar, as is the bias. We demonstrate this for our fiducial quasar by considering two BLR models that are implemented using Eq. (1) while holding all other BLR properties fixed at their aforementioned values: a phenomenological model based on the quasar composite spectrum of Vanden Berk et al. (2001) and a physical model for the BLR emission that is based on photoionization calculations (Korista & Goad 2001, 2019). In the former model, a broken power-law is fit to the data, and the residuals, normalized to the local power-law fit, are attributed to the BLR emission (Vanden Berk et al. 2001) and hence set α(λ). A smoothed version of the model is shown in the inset of Fig. 4 to allow an easier comparison with photometric data. In the second model, an independent estimate for α(λ) is provided, which relies on the physical locally optimally emitting cloud model for the BLR (Baldwin et al. 1995), whose α(λ) predictions factor in many relevant processes, but do not include emission lines and line blends (Korista & Goad 2019). In both of the models, α(λ)≳0.1 over much of the UV-optical range and may reach α ≃ 0.4 near the radiative recombination-edge thresholds of the Balmer and Paschen series. These models are used in conjunction with Eq. (1) and the simulated magnification map to predict the biases shown in Fig. 4.

thumbnail Fig. 4.

Wavelength-dependent overestimation factors of disk sizes in the presence of a BLR component. The effect of the power-law index of the emissivity profile (Eq. (1)) is small within the range [ − 3, 0] (shaded regions). The inset shows the relative contribution of the BLR to the flux in the band, α, according to two models: the observationally motivated phenomenological line emission model of Vanden Berk et al. (2001) is shown in blue, and the more physical BLR model of Korista & Goad (2019) in red.

The two models for α(λ), despite being very different in construction, lead to qualitatively similar conclusions. The wavelength-dependent overestimation factor of the inferred disk size for our fiducial quasar is very substantial at rest-UV energies (Fig. 4), where the disk emission is more compact. The bias is clearly affected by major spectral BLR features such as radiative-recombination edges and strong lines and line blends. Specifically, overestimation factors of ≳4 are expected at ≲4000 Å where the small blue bump dominates (Wills et al. 1985). The contribution from the Rayleigh-scattering wings of Lyα (see the wavelength range below ∼1500 Å), may lead to significant biases according to these models. The Paschen edge at ∼8000 Å leads to a smaller bias than the Balmer edge because the underlying disk scale at long wavelengths is larger. As demonstrated above, we expect sources whose luminosity is greater (smaller) than the fiducial value to be prone to smaller (larger) biases.

Varying the emissivity profile of the diffuse emission component (i.e., β) has a modest effect on the bias level for α < 0.5. For example, a β = 0 model shifts the emission to larger radii, thereby increasing the bias by ∼10% relative to our fiducial case. Conversely, a steeper emissivity profile results in a somewhat reduced bias level. For narrowband data, where α > 0.5 values are encountered at the spectral locations of prominent emission lines, the implied half-light radius may be underestimated by a factor of a few for β > −2 (not shown).

4. A quantitative assessment of reported AD sizes

Here we aim to check whether reported size estimations for the UV-optical disk in quasars, which largely disagree with AD models for these objects, could be reconciled with SS73 theory when considering the biases introduced by diffuse BLR emission. A refined BLR model is presented below, followed by a comparison between AD size estimations for the Blackburne et al. (2011) sample and the simulations that include BLR emission. We emphasize that it is not our intention to solve the full observational problem for particular lenses, which is exacerbated by (unknown) intrinsic quasar variations with time-delays between the images, and is subject to uncertainties in the macro- and micro-lens models.

4.1. A realistic model for the BLR

We aim to incorporate a realistic model for the diffuse emission from the BLR, which includes all major emission lines, line-blends, and continuum emission, and is consistent with BLR phenomenology of major emission lines across a wide energy range. The BLR model considered here assumes radiation-pressure–confined (RPC) non-accelerating clouds (Dopita et al. 2002), which are ionization-bound, and for which emission and scattering from deep within the cold neutral regions of the clouds is neglected (Baskin et al. 2014). Specifically, a one-dimensional slab is illuminated from one side by the ionizing disk and corona radiation, thereby exerting radiation pressure force on the illuminated layers. Under hydrostatic equilibrium, pressure increases within the deeper cloud layers to counter-balance the radiation pressure force, thereby leading to high compression rates, and to the gas pressure in the optically thick cloud layers being comparable to the radiation pressure, as observed (Baskin et al. 2014). A salient feature of such models is a wide range of ionization levels characterizing individual clouds: highly ionized gas near the illuminated surface of the cloud to neutral gas deep in its dense and optically thick layers.

The equilibrium structure of the cloud is self-consistently calculated here using the Cloudy v17.02 photoionization code (Ferland et al. 2017) including all relevant bound-bound, bound-free, and free-free opacities, as well as electron scattering. A dust-free medium is assumed, as appropriate for much of the line-emitting BLR (Netzer & Laor 1993). The cloud’s structure depends little on the boundary conditions at the location of the illuminated surface so long as the gas pressure at the boundary satisfies, Pgas(r)≪F(r)/c, where F(r) is the total flux and c is the speed of light. In our calculations, we consider geometrically thin cloud configurations so that the thickness of the emitting layer, δr, satisfies δr/r < 0.1 (Baskin et al. 2014). In this limit geometrical attenuation inside the cloud is negligible. In all models considered here, the illuminated surface of the cloud is at or near the Compton temperature, which is ≳106 K for a standard spectral energy distribution (SED) of a type-I quasar (Chelouche et al. 2019; see our Fig. 5). The gas composition is characterized by a metal-to-hydrogen ratio of 2.5 times the solar ratio (Warner et al. 2003). The potential dependence of the ionizing SED on BH mass and accretion rate is ignored. Upon the calculation of the cloud’s radial structure, its emissivity through the illuminated surface (i.e., the inwardly emitted spectrum) is obtained from infrared to X-ray energies (Fig. 5). We assume that sightlines toward the observer do not pass through BLR clouds, and hence the adopted BLR configuration is akin to the bowl geometry (Goad et al. 2012), which is viewed face-on and is consistent with recent findings for the BLR geometry in low-luminosity sources (Williams et al. 2018).

thumbnail Fig. 5.

Inwardly emitted spectra through the illuminated surface of a BLR cloud. Different curves correspond to the emission of clouds at different distances, r, from the ionizing source. The numbers correspond to log(r/rin). Dotted lines show the ionizing SED used for the calculations, relative to which full coverage BLR emission is assumed.

4.2. Model calibration

The combined spectrum of many individual clouds depends on their spatial distribution. Here we assume that the gas covers a fraction of the quasar sky, which is quantified by a radially dependent covering fraction, C(r), and that the quasar emits isotropically. Specifically, the cumulative covering fraction of the gas up to distance rin ≤ r ≤ rout is given by

C ( r ) = C 1 + C 2 1 r / r in ξ δ d ξ , $$ \begin{aligned} C(r) = C_1 + C_2\int _{1}^{r/r_{\rm in}}\xi ^\delta \mathrm{d}\xi , \end{aligned} $$(3)

where C1 and C2 are the parameters of the model. The power-law index, δ is rather loosely constrained by observations (certainly if C1 > 0), and we take δ = −1 (Baskin et al. 2014), so that C(r) = C1 + C2ln(r/rin). Cases for which C1 > 0 allow for the inclusion of an inner BLR “funnel”, which is motivated by recent RM results (Chelouche et al. 2019), while C1 = 0 is the more commonly employed model (Baskin et al. 2014; Netzer 2020). Specifically, we define the RPC model with C1 = 0 and the modified RPC (mRPC) model with C1 = 0.35, which results in an inner funnel at the radius where a disk wind is launched (Czerny & Hryniewicz 2011; Chelouche et al. 2019). This model is reminiscent of RPC models that are characterized by δ < −1 values (Netzer 2020).

The outer BLR boundary, rout, is set by the sublimation radius for dust, such that r out 10 3 L opt , 45 1 / 2 $ r_{\mathrm{out}}\simeq 10^3\,L_{\mathrm{opt,45}}^{1/2} $ light-days (see Fig. 6)3. Realistically, a range of sublimation radii exists for different grain sizes and compositions, and yet a more accurate treatment of the problem is beyond the scope of the present work. Our neglect of BLR emission beyond rout is justified by the large dust opacity at extreme-UV energies, which effectively competes with that of the gas over ionizing photons (Netzer & Laor 1993). The suppression of gas emission is further enhanced by virtue of RPC models experiencing a much faster compression and hence being characterized by smaller emitting columns in the presence of dust, as these are inversely proportional to the effective absorption cross-section (Baskin et al. 2014). We neglect dust emission in our model, which is expected to contribute at wavelengths ≳1 μm. Further, with rout being typically significantly larger than the Einstein radius for our targets, the results are insensitive to the particular value assumed. Here we take rout/rin ∼ 30, which is based on arguments related to the effective launching of disk material to great heights above the AD surface beyond the dust sublimation radius in the disk (Czerny & Hryniewicz 2011).

thumbnail Fig. 6.

Quasar SB as a function of the radial coordinate and rest wavelength for an Lopt = 1043 erg s−1 source. Color-coding is proportional to the logarithm of the SB (see the colorbar using arbitrarily normalized units). Note the presence of the inner stable co-rotating orbit, the innermost stable circular orbit (ISCO), which is included here for a ∼107 solar-mass BH.

As our focus in this paper is on broadband emission from BLR gas, the particular kinematic model is immaterial, and simplified gas kinematics is assumed, whereby the velocity dispersion of the medium is set by the virial theorem, and Gaussian line profiles are assumed, except for Lyα where an approximate treatment of the extended wings is included (Ferland et al. 2013). The backward emitted spectrum emerging from the illuminated face of a cloud placed at several distances from the ionizing source shows emission lines, radiative recombination edges, and free-free emission covering a broad energy range (Fig. 5). Notably, free-free emission combined with bound-free emission is substantial in the near-UV and optical bands, especially from gas that is located at r ≳ rin, where the substantial compression occurs, and gas densities exceeding 1012 cm−3 are reached deep within the cloud layers. From such regions, Balmer line emission is collisionally suppressed (Baskin et al. 2014; Chelouche et al. 2019). Iron emission is substantial across the full spectral range (including the Kα line at X-ray wavelengths) but is difficult to predict in the optical and UV bands, and should be cautiously interpreted.

The free parameters of the model are therefore C1, C2, which are calibrated using ample reverberation and spectroscopic data pertaining to active galactic nuclei (AGN) with Lopt, 45 ≳ 0.01, and are given in Table 1. Specifically, the equivalent width (EW) of rest-UV emission lines was taken from the Dietrich et al. (2002) sample for intermediate- to high-redshift sources of comparable luminosity. For the shortest wavelength emission lines (e.g., O VIλ1035) an extrapolation to low-luminosity sources was used, which is consistent with the data for local AGN (Scott et al. 2005). We note the lack of discernible redshift dependence in the EWs of prominent emission lines (Dietrich et al. 2002) and take the scatter in the EW-redshift data as a measure of the EW uncertainty. The EW of prominent optical lines is affected by the host-galaxy contribution (Stern & Laor 2012; Shen 2016), and we adopt EW(Hβ) = 80 ± 20 Å as appropriate over a broad range of quasar luminosities (Dietrich et al. 2002; Marziani et al. 2009; Shen 2016), and with an uncertainty that is consistent with the scatter seen in luminous sources for which the host galaxy contribution is small (Marziani et al. 2009). We take EW(Hα) = 390 Å, which is based on an IR sample of luminous quasars (Espey et al. 1989) and consistent with recent composites for objects with a small host contribution to the flux (Glikman et al. 2006; Stern & Laor 2012; Selsing et al. 2016). We scale the relative uncertainty in Hα to that of Hβ. The EW of Paβ is based on recent composites (Glikman et al. 2006) after correcting for a ∼50% host (and dust) contribution to the flux (Selsing et al. 2016) and conservatively estimating the uncertainty at the 50% level. EW values for the X-ray lines were taken from X-ray grating observations of low-luminosity AGN (Kaspi et al. 2002; Scott et al. 2005). For more luminous sources, the Baldwin effect is observed, whereby the EWs of most emission lines decrease such that EW L UV ϵ $ {\rm EW} \propto L_{\rm UV}^\epsilon $, where ϵ ∼ −0.2 (Baldwin 1977; Dietrich et al. 2002; Page et al. 2004).

Table 1.

Emission line constraints for an Lopt = 1043 erg s−1 source.

A further set of constraints on our models involves the characteristic size scale, RBLR, of the region that gives rise to the emission of a particular line (see Table 1). Thus, different emission lines are characterized by different RBLR values with that of Hβ often being taken as a representative size for the entire BLR (Bentz et al. 2013). Here, reverberation sizes were collected from the literature (Bentz et al. 2013) and rescaled by L opt 1 / 2 $ L_{\mathrm{opt}}^{1/2} $ to correspond to a Lopt, 45 = 0.01 source. Specifically, RBLR(Hβ) = 10 light-days, and RBLR(Hα) = 15 light-days were used (Bentz et al. 2010, 2013). The size of the CIV λ1548 emission region was taken to be identical to that of Lyα and SiIV λ1403, and four times smaller than the Hβ’s (De Rosa et al. 2015; Lira et al. 2018). The CIII] λ1900 emission line region was taken to be twice that of the CIV line (Lira et al. 2018). Sizes for HeII lines were rescaled from recent reverberation studies (Bentz et al. 2010; De Rosa et al. 2015). Constraints for the MgII line are unreliable on accounts of blending with the iron blends, and hence not included in the study. We refrain from using the recent GRAVITY data for the size of the Paschen-series emitting region since a proper comparison between reverberation data and interferometric data is yet to be performed. Model predictions for the emissivity-weighted radius are then compared to the measured RBLR values for the emission lines listed in Table 1, as are the EW predictions. An optimal solution, in the χ2 sense, is sought in the phase space, whose predictions are given in Table 1.

Final model parameterization yields C(rout) = 0.35 for the RPC model, and C(rout) = 0.7 for the mRPC model. The main difference between the models is the larger contribution to the UV-optical continuum from the BLR in the mRPC case, and therefore to slightly lower EW values of prominent lines. The SB of the best-fit mRPC model is shown in Fig. 6 as a function of the radial coordinate and rest wavelength for a Lopt, 45 = 1 source.

Scaling the model to luminous quasars is done here by assuming that the system size scales with L1/2 (Bentz et al. 2013; Fian et al. 2023, and below). Unless otherwise stated, we assume that the covering fractions, C1 and C2, are proportional to Lϵ with ϵ = −0.2, which is motivated by the observed Baldwin effect (Baldwin et al. 1978; see our Table 1). Thus, in our model, the relative contribution of both line emission and non-disk continuum emission to the total flux is lower in more luminous quasars since the total covering fraction of the BLR over the quasar sky is smaller4. This is qualitatively consistent with the decreasing fraction of obscured to non-obscured sources with source luminosity (Merloni et al. 2014) that may be akin to the receding torus model (Lawrence 1991). Specifically, the mRPC model with its inner funnel predicts a covering fraction of ∼0.3 for our fiducial quasar, which is comparable to the observed fraction of X-ray and optically obscured sources among the quasar population (Merloni et al. 2014). We note, however, that we do not attempt to provide a physical explanation for the Baldwin effect, which is likely driven by additional factors relating, for example, to the shape of the ionizing spectrum and the gas metallicity (Korista et al. 1998).

Assuming a face-on azimuthally symmetric source configuration, we convolve the wavelength-dependent SB for our fiducial, Lopt, 45 = 1 quasar with a magnification map, and estimate the source size as per the approach defined in Sect. 2. The recovered wavelength-dependent r1/2 is shown in Fig. 7 after degrading the spectral resolution to qualitatively match broadband data. Predictions were calculated for several microlensing maps whose Einstein radii were set to an observationally motivated range of values (log(rE)∈[15.5, 17.0]). Significant deviations from the input AD’s half-light radius are apparent in all cases. Due to the different scales characterizing the quasar’s SB profile, which involve the AD and the BLR, estimations of the wavelength-dependent half-light radius that are based on SS73 models do not simply scale with the Einstein radius. Specifically, the larger the Einstein radius is compared to the relevant physical scales in the SB map, the larger is the overestimation factor for the disk size. This is particularly noticeable at short wavelengths, which tends to flatten the size-wavelength relation, as is observed in quasars with comparable luminosity (see Fig. 7; Blackburne et al. 2011). As a corollary, good understanding of quasar structure may be useful for assessing the mean mass of stars in the lens.

thumbnail Fig. 7.

Dependence of the deduced wavelength-dependent half-light radius, r1/2(λ), on the Einstein radius. The inferred r1/2(λ) from microlensing simulations are shown as black curves for different values of the Einstein radius, rE (see the legend). Overlaid are the microlensed-based size estimates for two Lopt, 45 ≃ 1 quasars as well as the predicted SS73 wavelength-dependent disk size (dashed line).

4.3. Comparison to the Blackburne et al. (2011) sample

Microlensing-based size estimates for the continuum-emitting region in quasars were collected from the literature, where ∼90 wavelength-dependent size measurements for 12 strongly lensed quasars were augmented by ∼20 additional measurements for the same sources from various sources (Table 2). All measurements were scaled to an Einstein radius of a solar-mass star to allow the data sets to be compared. Specifically, the main sample analyzed here is based on the data provided in Blackburne et al. (2011). Table 2 presents additional size measurements for eight systems in that sample, along with their corresponding rest wavelength. Conversions have been made when necessary to translate all reported scales to SS73-disk half-light radii (as opposed to, e.g., Gaussian or uniform disk), and these are shown in Fig. 8. We find that the size estimates drawn from different studies using different methods are consistent in most cases (compare the red and black points in Fig. 8).

thumbnail Fig. 8.

Microlensing-based size estimates versus theoretical predictions for 12 strongly lensed quasars. The expected half-light radius for an SS73 disk for each source is shown as a blue line. The red curves show model predictions that include all relevant bound-bound, bound-free, and free-free emission processes and are convolved with a typical broadband throughput kernel to mimic broadband data. Three models are shown: a RPC model (Baskin et al. 2014) that accounts for the Baldwin relation (solid line), a similar model but with the addition of an inner BLR funnel (Chelouche et al. 2019; dotted-dashed curve), and a similar model to the latter but ignoring the effect of the Baldwin relation and setting the BLR properties to those given in Table 1, which characterize low-luminosity sources (dashed curve). Tiles are sorted by increasing object luminosity from left to right and from top to bottom. Microlensing data are shown with respective error bars (68% confidence intervals) as red squares (Blackburne et al. 2011) and as black points (see Table 2). Upper and lower limits are denoted by downward and upward pointing triangles, respectively. The Einstein radius for each lens is denoted by the horizontal black dotted line. Shaded regions mark the scale range subtended by the BLR (rout exceeds the axis limits in all cases).

Table 2.

Half-light radius, r1/2, estimates (at the 68% confidence level) from the literature (scaled to microlenses with a mean mass of 1 M).

We next compare detailed model predictions to published microlensing estimations for the half-light radius of the continuum-emitting region for the above sources. Models are rescaled according to the magnification corrected luminosity for individual lensed quasars (Eq. (2)), and are then convolved with a magnification map whose Einstein radius matches the one provided for each source in the sample, which are all in agreement within 0.2 dex. Model predictions for the recovered wavelength-dependent sizes for each source are shown in Fig. 8 (red curves), where the salient features are qualitatively reproduces, namely, a flatter size-wavelength dependence than assumed by pure AD theory (blue line), and implied AD sizes that are considerably larger than predicted by theory, especially at short wavelengths due to the BLR contribution to the bands. Predictions for three model variants for the BLR are provided, which differ by their treatment of the Baldwin effect (ϵ = −0.2 vs. ϵ = 0) and/or the presence of an inner BLR funnel as motivated by recent reverberation results (RPC model with C1 = 0 vs. an mRPC model with C1 = 0.35).

All models are qualitatively consistent with the half-light source sizes as inferred from microlensing data (red and black data points), although the model that does not account for the Baldwin effect and assumes as strong emission lines as in low-luminosity sources appears to lie above the short-wavelength data for the two brightest objects in the sample. At long wavelengths, microlensing size estimates are less biased and tend to converge to the predicted SS73 disk sizes based on the sources’ luminosity in the more luminous objects. We note that in the case of SDSS J1330, the theoretical predictions lie significantly above the upper limits on the disk size, perhaps due to an overestimation of the quasar luminosity, or due to a less probable microlensing configuration, or a combination thereof (this source is unique in the sample of Blackburne et al. 2011 as it does not have corresponding X-ray data to further constrain the microlensing model).

5. Conclusions

The comparison of our simulations to available data in the literature points to a possible solution to the oversized AD conundrum, which is related to the interpretation of the observed signals rather than to disk physics. Current AD and BLR models, limited and uncertain as they may be, are in qualitative agreement with microlensing-deduced continuum emission-region sizes. While additional causes for size overestimation exist, such as the effect of the BLR reverberation signal (see Paic et al. 2022) on forward modeling schemes (Kochanek 2004), we show that the mere contribution of the BLR to the continuum signal is able to largely account for the implied disk overestimation factors as well as to give rise to the flat size-wavelength relations observed in a sample of 12 lenses. From this, combined with recent RM studies of low-luminosity sources (Chelouche et al. 2019), a coherent picture emerges wherein the BLR contribution to the signal must be accounted for before any constraints can be placed on AD physics. This conclusion seems to apply whether the total nuclear emission from the quasar is being analyzed, as is the case with microlensing, or only its variable component, as is the case for RM studies.

Our findings appear to significantly weaken the tension between AD theory and observations. They suggest that the steep disk temperature profiles found by microlensing studies are erroneous as the data are largely affected by the BLR, which does not emit as an aggregate of black bodies and hence does not obey a temperature-wavelength relation. The BLR models considered here further suggest that microlensing is a promising tool for constraining the BLR physics, and in particular for assessing the elusive contribution of the BLR to continuum emission from quasars. Microlensing may therefore shed light on the long-sought origin of the Baldwin effect (Baldwin 1977; Baldwin et al. 1978) by probing the dependence of the diffuse BLR continuum emission on the source luminosity and emission-line phenomenology, and revealing the mechanisms underlying the quasar eigenvector phenomenology (Sulentic et al. 2000). Shedding light on the above is crucial for revealing the inner workings of quasars, with implications for their usage as proxies of the BH census and, perhaps, as standard cosmological rulers. Our results indicate that microlensing may already provide useful constraints on AD physics in sources whose BLR emission is weak and its extent much larger than the typical Einstein radius, such as in the most luminous sources (see, e.g., MG J0414 in Fig. 8). In particular, microlensing could shed light on the role of metallicity in setting the inner disk properties (Jiang et al. 2016) as well as on the effects of both outflows from the disk atmosphere and the radiative transfer in the disk atmosphere on the disk phenomenology (Laor & Davis 2014; Hall et al. 2018). Microlensing may also help elucidate the recently proposed connections between the AD, the BLR, and the torus in quasars (Czerny & Hryniewicz 2011; Czerny et al. 2016; Baskin & Laor 2018; Dorodnitsyn & Kallman 2021). Conversely, with quasar physics better established, it may be possible to independently estimate the mean mass of microlenses in the lensing galaxy from the shape and amplitude of the size-wavelength relation for quasars. This has implications for probing stellar populations in the outskirts of the lensing objects. With upcoming all-sky surveys, mapping quasar interiors using microlensing is expected to become widespread, and data may be able to tap into higher moments of the magnification distribution, potentially alleviating some of the biases and uncertainties discussed here.


1

A compact source could, for example, inhabit a particularly calm region of the microlensing map, leading to little microlensing induced variability, thereby behaving as a much larger source occupying a more typical region of the map.

2

A finite dispersion of the magnification distributions is also apparent in the limiting case of α = 1, thereby distinguishing the BLR model from highly extended sources, such as the host galaxy, whose contribution has not been included in the model.

3

Here we used the fact that the bolometric luminosity is ≃8 times the optical luminosity for the chosen SED.

4

The EW of the Balmer lines is less affected since the dominant free-free emission and bound-free emission at optical wavelengths is similarly suppressed.

Acknowledgments

We thank E. Mediavilla for generously providing the IPM code that enabled us to compute the microlensing magnification maps. We thank an anonymous referee for helpful suggestions that improved the conclusions section. This research was supported by the grant PID2020-118687GB-C32, financed by the Spanish Ministerio de Ciencia e Innovación. Research by D.C. and S.K. is financially supported by the German Science Foundation (DFG) grants HA3555-14/1, CH71-34-3, and by the Israeli Science Foundation grant no. 2398/19. Computations made use of high-performance computing facilities at the University of Haifa, which are partly supported by a grant from the Israeli Science Foundation (grant 2155/15).

References

  1. Baldwin, J. A. 1977, ApJ, 214, 679 [NASA ADS] [CrossRef] [Google Scholar]
  2. Baldwin, J. A., Burke, W. L., Gaskell, C. M., & Wampler, E. J. 1978, Nature, 273, 431 [NASA ADS] [CrossRef] [Google Scholar]
  3. Baldwin, J., Ferland, G., Korista, K., & Verner, D. 1995, ApJ, 455, L119 [NASA ADS] [Google Scholar]
  4. Baskin, A., & Laor, A. 2018, MNRAS, 474, 1970 [Google Scholar]
  5. Baskin, A., Laor, A., & Stern, J. 2014, MNRAS, 438, 604 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bate, N. F., Webster, R. L., & Wyithe, J. S. B. 2007, MNRAS, 381, 1591 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bate, N. F., Floyd, D. J. E., Webster, R. L., & Wyithe, J. S. B. 2008, MNRAS, 391, 1955 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bate, N. F., Vernardos, G., O’Dowd, M. J., et al. 2018, MNRAS, 479, 4796 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bentz, M. C., Walsh, J. L., Barth, A. J., et al. 2010, ApJ, 716, 993 [Google Scholar]
  10. Bentz, M. C., Denney, K. D., Grier, C. J., et al. 2013, ApJ, 767, 149 [Google Scholar]
  11. Blackburne, J. A., Pooley, D., Rappaport, S., & Schechter, P. L. 2011, ApJ, 729, 34 [Google Scholar]
  12. Blandford, R. D., & McKee, C. F. 1982, ApJ, 255, 419 [Google Scholar]
  13. Cackett, E. M., Chiang, C.-Y., McHardy, I., et al. 2018, ApJ, 857, 53 [NASA ADS] [CrossRef] [Google Scholar]
  14. Chan, J. H. H., Rojas, K., Millon, M., et al. 2021, A&A, 647, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Chelouche, D., Pozo Nuñez, F., & Kaspi, S. 2019, Nat. Astron., 3, 251 [NASA ADS] [CrossRef] [Google Scholar]
  16. Cornachione, M. A., & Morgan, C. W. 2020, ApJ, 895, 93 [Google Scholar]
  17. Cornachione, M. A., Morgan, C. W., Burger, H. R., et al. 2020a, ApJ, 905, 7 [NASA ADS] [CrossRef] [Google Scholar]
  18. Cornachione, M. A., Morgan, C. W., Millon, M., et al. 2020b, ApJ, 895, 125 [Google Scholar]
  19. Czerny, B., & Hryniewicz, K. 2011, A&A, 525, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Czerny, B., Du, P., Wang, J.-M., & Karas, V. 2016, ApJ, 832, 15 [NASA ADS] [CrossRef] [Google Scholar]
  21. Dai, X., Kochanek, C. S., Chartas, G., et al. 2010, ApJ, 709, 278 [Google Scholar]
  22. Davis, S. W., & Laor, A. 2011, ApJ, 728, 98 [NASA ADS] [CrossRef] [Google Scholar]
  23. De Rosa, G., Peterson, B. M., Ely, J., et al. 2015, ApJ, 806, 128 [NASA ADS] [CrossRef] [Google Scholar]
  24. Dexter, J., & Agol, E. 2011, ApJ, 727, L24 [Google Scholar]
  25. Dietrich, M., Hamann, F., Shields, J. C., et al. 2002, ApJ, 581, 912 [NASA ADS] [CrossRef] [Google Scholar]
  26. Dopita, M. A., Groves, B. A., Sutherland, R. S., Binette, L., & Cecil, G. 2002, ApJ, 572, 753 [NASA ADS] [CrossRef] [Google Scholar]
  27. Dorodnitsyn, A., & Kallman, T. 2021, ApJ, 910, 67 [NASA ADS] [CrossRef] [Google Scholar]
  28. Edelson, R., Gelbord, J., Cackett, E., et al. 2019, ApJ, 870, 123 [NASA ADS] [CrossRef] [Google Scholar]
  29. Espey, B. R., Carswell, R. F., Bailey, J. A., Smith, M. G., & Ward, M. J. 1989, ApJ, 342, 666 [NASA ADS] [CrossRef] [Google Scholar]
  30. Fausnaugh, M. M., Starkey, D. A., Horne, K., et al. 2018, ApJ, 854, 107 [NASA ADS] [CrossRef] [Google Scholar]
  31. Ferland, G. J., Porter, R. L., van Hoof, P. A. M., et al. 2013, Rev. Mex. Astron. Astrofis., 49, 137 [Google Scholar]
  32. Ferland, G. J., Chatzikos, M., Guzmán, F., et al. 2017, Rev. Mex. Astron. Astrofis., 53, 385 [NASA ADS] [Google Scholar]
  33. Fian, C., Mediavilla, E., Hanslmeier, A., et al. 2016, ApJ, 830, 149 [Google Scholar]
  34. Fian, C., Mediavilla, E., Jiménez-Vicente, J., Muñoz, J. A., & Hanslmeier, A. 2018, ApJ, 869, 132 [NASA ADS] [CrossRef] [Google Scholar]
  35. Fian, C., Mediavilla, E., Jiménez-Vicente, J., et al. 2021, A&A, 654, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Fian, C., Chelouche, D., Kaspi, S., et al. 2022, A&A, 659, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Fian, C., Chelouche, D., Kaspi, S., et al. 2023, A&A, 672, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Floyd, D. J. E., Bate, N. F., & Webster, R. L. 2009, MNRAS, 398, 233 [NASA ADS] [CrossRef] [Google Scholar]
  39. Gardner, E., & Done, C. 2017, MNRAS, 470, 3591 [NASA ADS] [CrossRef] [Google Scholar]
  40. Gaskell, C. M., & Harrington, P. Z. 2018, MNRAS, 478, 1660 [NASA ADS] [CrossRef] [Google Scholar]
  41. Glikman, E., Helfand, D. J., & White, R. L. 2006, ApJ, 640, 579 [NASA ADS] [CrossRef] [Google Scholar]
  42. Goad, M. R., Korista, K. T., & Ruff, A. J. 2012, MNRAS, 426, 3086 [Google Scholar]
  43. Hall, P. B., Sarrouh, G. T., & Horne, K. 2018, ApJ, 854, 93 [Google Scholar]
  44. Jiang, Y.-F., Davis, S. W., & Stone, J. M. 2016, ApJ, 827, 10 [NASA ADS] [CrossRef] [Google Scholar]
  45. Jiménez-Vicente, J., Mediavilla, E., Muñoz, J. A., & Kochanek, C. S. 2012, ApJ, 751, 106 [Google Scholar]
  46. Jiménez-Vicente, J., Mediavilla, E., Kochanek, C. S., et al. 2014, ApJ, 783, 47 [Google Scholar]
  47. Jiménez-Vicente, J., Mediavilla, E., Kochanek, C. S., & Muñoz, J. A. 2015, ApJ, 799, 149 [Google Scholar]
  48. Kara, E., Mehdipour, M., Kriss, G. A., et al. 2021, ApJ, 922, 151 [NASA ADS] [CrossRef] [Google Scholar]
  49. Kara, E., Barth, A. J., Cackett, E. M., et al. 2023, ApJ, 947, 62 [NASA ADS] [CrossRef] [Google Scholar]
  50. Kaspi, S., Brandt, W. N., George, I. M., et al. 2002, ApJ, 574, 643 [NASA ADS] [CrossRef] [Google Scholar]
  51. Kochanek, C. S. 2004, ApJ, 605, 58 [Google Scholar]
  52. Korista, K. T., & Goad, M. R. 2001, ApJ, 553, 695 [Google Scholar]
  53. Korista, K. T., & Goad, M. R. 2019, MNRAS, 489, 5284 [CrossRef] [Google Scholar]
  54. Korista, K., Baldwin, J., & Ferland, G. 1998, ApJ, 507, 24 [NASA ADS] [CrossRef] [Google Scholar]
  55. Laor, A., & Davis, S. W. 2011, MNRAS, 417, 681 [NASA ADS] [CrossRef] [Google Scholar]
  56. Laor, A., & Davis, S. W. 2014, MNRAS, 438, 3024 [NASA ADS] [CrossRef] [Google Scholar]
  57. Lawrence, A. 1991, MNRAS, 252, 586 [NASA ADS] [Google Scholar]
  58. Lewis, G. F., & Irwin, M. J. 1995, MNRAS, 276, 103 [NASA ADS] [Google Scholar]
  59. Li, Y.-P., Yuan, F., & Dai, X. 2019, MNRAS, 483, 2275 [NASA ADS] [CrossRef] [Google Scholar]
  60. Lira, P., Kaspi, S., Netzer, H., et al. 2018, ApJ, 865, 56 [NASA ADS] [CrossRef] [Google Scholar]
  61. MacLeod, C. L., Morgan, C. W., Mosquera, A., et al. 2015, ApJ, 806, 258 [NASA ADS] [CrossRef] [Google Scholar]
  62. Marziani, P., Sulentic, J. W., Stirpe, G. M., Zamfir, S., & Calvani, M. 2009, A&A, 495, 83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Mediavilla, E., Muñoz, J. A., Lopez, P., et al. 2006, ApJ, 653, 942 [NASA ADS] [CrossRef] [Google Scholar]
  64. Mediavilla, E., Mediavilla, T., Muñoz, J. A., et al. 2011, ApJ, 741, 42 [NASA ADS] [CrossRef] [Google Scholar]
  65. Merloni, A., Bongiorno, A., Brusa, M., et al. 2014, MNRAS, 437, 3550 [Google Scholar]
  66. Morgan, C. W., Kochanek, C. S., Morgan, N. D., & Falco, E. E. 2010, ApJ, 712, 1129 [Google Scholar]
  67. Mortonson, M. J., Schechter, P. L., & Wambsganss, J. 2005, ApJ, 628, 594 [NASA ADS] [CrossRef] [Google Scholar]
  68. Mosquera, A. M., & Kochanek, C. S. 2011, ApJ, 738, 96 [Google Scholar]
  69. Motta, V., Mediavilla, E., Rojas, K., et al. 2017, ApJ, 835, 132 [NASA ADS] [CrossRef] [Google Scholar]
  70. Netzer, H. 2020, MNRAS, 494, 1611 [NASA ADS] [CrossRef] [Google Scholar]
  71. Netzer, H. 2022, MNRAS, 509, 2637 [NASA ADS] [Google Scholar]
  72. Netzer, H., & Laor, A. 1993, ApJ, 404, L51 [NASA ADS] [CrossRef] [Google Scholar]
  73. Page, K. L., O’Brien, P. T., Reeves, J. N., & Turner, M. J. L. 2004, MNRAS, 347, 316 [NASA ADS] [CrossRef] [Google Scholar]
  74. Paic, E., Vernardos, G., Sluse, D., et al. 2022, A&A, 659, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Peterson, B. M. 1993, PASP, 105, 247 [NASA ADS] [CrossRef] [Google Scholar]
  76. Rojas, K., Motta, V., Mediavilla, E., et al. 2020, ApJ, 890, 3 [Google Scholar]
  77. Scott, J. E., Kriss, G. A., Lee, J. C., et al. 2005, ApJ, 634, 193 [NASA ADS] [CrossRef] [Google Scholar]
  78. Selsing, J., Fynbo, J. P. U., Christensen, L., & Krogager, J. K. 2016, A&A, 585, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  79. Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 500, 33 [NASA ADS] [Google Scholar]
  80. Shen, Y. 2016, ApJ, 817, 55 [NASA ADS] [CrossRef] [Google Scholar]
  81. Stern, J., & Laor, A. 2012, MNRAS, 423, 600 [NASA ADS] [CrossRef] [Google Scholar]
  82. Sulentic, J. W., Marziani, P., & Dultzin-Hacyan, D. 2000, ARA&A, 38, 521 [Google Scholar]
  83. Vanden Berk, D. E., Richards, G. T., Bauer, A., et al. 2001, AJ, 122, 549 [Google Scholar]
  84. Vernardos, G. 2018, MNRAS, 480, 4675 [NASA ADS] [CrossRef] [Google Scholar]
  85. Wambsganss, J. 1999, J. Comput. Appl. Math., 109, 353 [NASA ADS] [CrossRef] [Google Scholar]
  86. Wambsganss, J. 2006, in Saas-Fee Advanced Course 33: Gravitational Lensing: Strong, Weak and Micro, eds. G. Meylan, P. Jetzer, P. North, et al. (Berlin: Springer), 453 [Google Scholar]
  87. Warner, C., Hamann, F., & Dietrich, M. 2003, ApJ, 596, 72 [NASA ADS] [CrossRef] [Google Scholar]
  88. Williams, P. R., Pancoast, A., Treu, T., et al. 2018, ApJ, 866, 75 [NASA ADS] [CrossRef] [Google Scholar]
  89. Wills, B. J., Netzer, H., & Wills, D. 1985, ApJ, 288, 94 [NASA ADS] [CrossRef] [Google Scholar]
  90. Wyithe, J. S. B., & Turner, E. L. 2001, MNRAS, 320, 21 [Google Scholar]

All Tables

Table 1.

Emission line constraints for an Lopt = 1043 erg s−1 source.

Table 2.

Half-light radius, r1/2, estimates (at the 68% confidence level) from the literature (scaled to microlenses with a mean mass of 1 M).

All Figures

thumbnail Fig. 1.

Billion-pixel magnification map showing a complex caustic pattern for the positive-parity image, with brighter colors indicating larger magnifications. The insets show successive blow-ups of the maps with relevant quasar emission regions overlaid.

In the text
thumbnail Fig. 2.

Magnification distributions (relative to the mean magnification of the image) for several assumptions regarding the contribution of a diffuse emission component to the flux in some band, which is parameterized by α (see the main text). Bold distribution curves mark α = 0.0 and α = 1.0 cases for each image parity (see the legend), with the thin curves separated by Δα = 0.2 intervals. The top panels show the emissivity profiles in log-units for three models, with their α-values denoted. Color-coding is arbitrarily but identically normalized in all panels. Note the reduced second moment of the magnification distribution function with increasing α.

In the text
thumbnail Fig. 3.

Simulations quantifying the overestimation factors of the disk scale and the half-light radius as a function of the relative contribution of a diffuse BLR emission component to the total continuum signal, α. The disk-scale overestimation factor is defined as the ratio of the value recovered via the magnification-distribution function that matches the input value, Rs, in our AD-BLR model. Results for the half-light radius are defined by the ratio of the value recovered for the best-matched pure AD model to that which characterizes the AD-BLR model.

In the text
thumbnail Fig. 4.

Wavelength-dependent overestimation factors of disk sizes in the presence of a BLR component. The effect of the power-law index of the emissivity profile (Eq. (1)) is small within the range [ − 3, 0] (shaded regions). The inset shows the relative contribution of the BLR to the flux in the band, α, according to two models: the observationally motivated phenomenological line emission model of Vanden Berk et al. (2001) is shown in blue, and the more physical BLR model of Korista & Goad (2019) in red.

In the text
thumbnail Fig. 5.

Inwardly emitted spectra through the illuminated surface of a BLR cloud. Different curves correspond to the emission of clouds at different distances, r, from the ionizing source. The numbers correspond to log(r/rin). Dotted lines show the ionizing SED used for the calculations, relative to which full coverage BLR emission is assumed.

In the text
thumbnail Fig. 6.

Quasar SB as a function of the radial coordinate and rest wavelength for an Lopt = 1043 erg s−1 source. Color-coding is proportional to the logarithm of the SB (see the colorbar using arbitrarily normalized units). Note the presence of the inner stable co-rotating orbit, the innermost stable circular orbit (ISCO), which is included here for a ∼107 solar-mass BH.

In the text
thumbnail Fig. 7.

Dependence of the deduced wavelength-dependent half-light radius, r1/2(λ), on the Einstein radius. The inferred r1/2(λ) from microlensing simulations are shown as black curves for different values of the Einstein radius, rE (see the legend). Overlaid are the microlensed-based size estimates for two Lopt, 45 ≃ 1 quasars as well as the predicted SS73 wavelength-dependent disk size (dashed line).

In the text
thumbnail Fig. 8.

Microlensing-based size estimates versus theoretical predictions for 12 strongly lensed quasars. The expected half-light radius for an SS73 disk for each source is shown as a blue line. The red curves show model predictions that include all relevant bound-bound, bound-free, and free-free emission processes and are convolved with a typical broadband throughput kernel to mimic broadband data. Three models are shown: a RPC model (Baskin et al. 2014) that accounts for the Baldwin relation (solid line), a similar model but with the addition of an inner BLR funnel (Chelouche et al. 2019; dotted-dashed curve), and a similar model to the latter but ignoring the effect of the Baldwin relation and setting the BLR properties to those given in Table 1, which characterize low-luminosity sources (dashed curve). Tiles are sorted by increasing object luminosity from left to right and from top to bottom. Microlensing data are shown with respective error bars (68% confidence intervals) as red squares (Blackburne et al. 2011) and as black points (see Table 2). Upper and lower limits are denoted by downward and upward pointing triangles, respectively. The Einstein radius for each lens is denoted by the horizontal black dotted line. Shaded regions mark the scale range subtended by the BLR (rout exceeds the axis limits in all cases).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.