Assuming a gravitational origin for the Fe III$\lambda\lambda$2039-2113 redshift and using microlensing based estimates of the size of the region emitting this feature, we obtain individual measurements of the virial factor, $f$, in 10 quasars. The average values for the Balmer lines, $\langle f_{H\beta}\rangle={\bf 0.43\pm 0.20}$ and $\langle f_{H\alpha}\rangle={\bf 0.50\pm 0.24}$, are in good agreement with the results of previous studies for objects with lines of comparable widths. In the case of Mg II, consistent results, $f_{Mg II} \sim {\bf 0.44}$, can be also obtained accepting a reasonable scaling for the size of the emitting region. The modeling of the cumulative histograms of individual measurements, $CDF(f)$, indicates a {relatively} high value for the ratio between isotropic and cylindrical motions, $a\sim {\bf 0.4}-0.7$. On the contrary, we find very large values of the virial factor associated to the Fe III$\lambda\lambda$2039-2113 blend, $f_{FeIII}=14.3\pm2.4$, which can be explained if this feature arises from a flattened nearly face-on structure, similar to the accretion disk.

We use spectroscopic observations of the gravitationally lensed systems SDSS0924+0219(BC), Q1355-2257(AB), and SDSS1029+2623(BC) to analyze microlensing and dust extinction in the observed components. We detect chromatic microlensing effects in the continuum and microlensing in the broad emission line profiles of the systems SDSS0924+0219(BC), and Q1355-2257(AB). Using magnification maps to simulate microlensing and modeling the emitting region as a Gaussian intensity profile with size $r_s \propto \lambda ^p$, we obtain the probability density functions for a logarithmic size prior at $\lambda_{rest-frame}=3533$ {\AA}. In the case of SDSS0924+0219, we obtain: $r_s = 4^{+3}_{-2}$ $\sqrt{M/M_{\odot}}$ light-days (at $1 \sigma$), which is larger than the range of other estimates, and $p = 0.8 \pm 0.2$ (at $1 \sigma$), which is smaller than predicted by the thin disk theory, but still in agreement with previous results. In the case of Q1355-2257 we obtain (at $1 \sigma$): $r_s = 3.6^{+3.0}_{-1.6}$ $\sqrt{M/M_{\odot}}$ light-days, which is also larger than the theoretical prediction, and $p = 2.0 \pm 0.7$ that is in agreement with the theory within errors. SDSS1029+2326 spectra show evidence of extinction, probably produced by a galaxy in the vicinity of image C. Fitting an extinction curve to the data we estimate $\Delta E \sim 0.2$ in agreement with previous results. We found no evidence of microlensing for this system.

We present new spectroscopic and polarimetric observations of the gravitational lens SDSS J1004+4112 taken with the 6m telescope of the Special Astrophysical Observatory (SAO, Russia). In order to explain the variability that is observed only in the blue wing of the C IV emission line, corresponding to image A, we analyze the spectroscopy and polarimetry of the four images of the lensed system. Spectra of the four images were taken in 2007, 2008, and 2018, and polarization was measured in the period 2014-2017. Additionally, we modeled the microlensing effect in the polarized light, assuming that the source of polarization is the equatorial scattering in the inner part of the torus. We find that a blue enhancement in the CIV line wings affects component A in all three epochs. We also find that the UV continuum of component D was amplified in the period 2007-2008, and that the red wings of CIII] and CIV appear brighter in D than in the other three components. We report significant changes in the polarization parameters of image D, which can be explained by microlensing.Our simulations of microlensing of an equatorial scattering region in the dusty torus can qualitatively explain the observed changes in the polarization degree and angle of image D. We do not detect significant variability in the polarization parameters of the other images (A, B, and C), although the averaged values of the polarization degree and angle are different for the different images. Microlensing of a broad line region model including a compact outflowing component can qualitatively explain the CIV blue wing enhancement (and variation) in component A. However, to confirmed this hypothesis, we need additional spectroscopic observation in future.

We test the recently proposed (Mediavilla et al. 2018) black hole mass scaling relationship based on the redshift {with respect to the quasar's rest frame} of the Fe III$\lambda\lambda$2039-2113 line blend. To this end, we fit this feature in the spectra of a well suited sample of quasars, observed with X-shooter at the Very Large Telescope (VLT), whose masses have been independently estimated using the virial theorem. For the quasars of this sample we consistently confirm the redshift of the Fe III$\lambda\lambda$2039-2113 blend and find that it correlates with the squared widths of H$\beta$, H$\alpha$ and Mg II, which are commonly used as a measure of $M_{BH}/R$ to determine masses from the virial theorem. The average differences between virial and Fe III$\lambda\lambda$2039-2113 redshift based masses are 0.18$\pm 0.21$ dex, 0.18$\pm 0.22$ dex and 0.14$\pm 0.21$ dex, when the full widths at half maximum (FWHM) of the H$\beta$, H$\alpha$ and MgII lines are, respectively, used. The difference is reduced to 0.10$\pm 0.16$ dex when the standard deviation, $\sigma$, of {the} MgII line is used, instead. We also study the high S/N composite quasar spectra of the Baryon Oscillation Spectroscopic Survey (BOSS), finding that the Fe III$\lambda\lambda$2039-2113 redshifts and Mg II squared widths, $FWHM_{MgII}^2$, match very well the correlation found for the individual quasar spectra observed with X-shooter. This correlation is expected if the redshift is gravitational.

Microlensing by stars within distant galaxies acting as strong gravitational lenses of multiply-imaged quasars, provides a unique and direct measurement of the internal structure of the lensed quasar on nano-arcsecond scales. The measurement relies on the temporal variation of high-magnification caustic crossings which vary on timescales of days to years. Multiwavelength observations provide information from distinct emission regions in the quasar. Through monitoring of these strong gravitational lenses, a full tomographic view can emerge with Astronomical-Unit scale resolution. Work to date has demonstrated the potential of this technique in about a dozen systems. In the 2020s there will be orders of magnitude more systems to work with. Monitoring of lens systems for caustic-crossing events to enable triggering of multi-platform, multi-wavelength observations in the 2020s will fulfill the potential of quasar microlensing as a unique and comprehensive probe of active black hole structure and dynamics.

Optical and X-ray observations of strongly gravitationally lensed quasars (especially when four separate images of the quasar are produced) determine not only the amount of matter in the lensing galaxy but also how much is in a smooth component and how much is composed of compact masses (e.g., stars, stellar remnants, primordial black holes, CDM sub-halos, and planets). Future optical surveys will discover hundreds to thousands of quadruply lensed quasars, and sensitive X-ray observations will unambiguously determine the ratio of smooth to clumpy matter at specific locations in the lensing galaxies and calibrate the stellar mass fundamental plane, providing a determination of the stellar $M/L$. A modest observing program with a sensitive, sub-arcsecond X-ray imager, combined with the planned optical observations, can make those determinations for a large number (hundreds) of the lensing galaxies, which will span a redshift range of $\sim$$0.25<z<1.5$

We present a measurement of the accretion disk size of the quadruple lensed quasar HE 0435-1223 from well-sampled 13-yr COSMOGRAIL optical light curves. Using accurate time delays for the images A, B, C, and D, we modeled and removed the intrinsic quasar variability, and found microlensing events of amplitude up to 0.6, 0.4, and 0.5 mag in the images A, C and D respectively. From the statistics of microlensing magnifications in these images we use Bayesian methods to estimate the size of the quasar accretion disk. We have inferred the half-light radius for the accretion disk using two different methods, $R_{1/2} = 7.6_{-1.1}^{+12.0} \sqrt{M/0.3M_{\odot}}$ light-days (histogram product) and $R_{1/2} = 7.7_{-1.2}^{+7.0} \sqrt{M/0.3M_{\odot}}$ light-days ($\chi^2$ criterion). The results are self-consistent and in good agreement with the continuum size predicted by single-epoch spectroscopy and previous studies making use of narrow-band photometry of HE 0435-1223.

We find that the Fe III$\lambda\lambda$2039-2113 spectral feature in quasars appears systematically redshifted by amounts accountable under the hypothesis of gravitational redshift induced by the central supermassive black hole. Our analysis of 27 composite spectra from the BOSS survey indicates that the redshift and the broadening of the lines in the Fe III$\lambda\lambda$2039-2113 blend roughly follow the expected correlation in the weak limit of Schwarzschild geometry for virialized kinematics. Assuming that the Fe III UV redshift provides a measure of $M_{BH}\over R$ (${\Delta \lambda\over \lambda}\simeq{3\over2}{G\over c^2} {M_{BH}\over R}$) and using different estimates of the emitting region size, $R$ (either from gravitational microlensing, reverberation mapping or from the scaling of size with intrinsic quasar luminosity), we obtain masses for 10 objects which are in agreement within uncertainties with previous mass estimates based on the virial theorem. Reverberation mapping estimates of the size of the Fe III$\lambda\lambda$2039-2113 emitting region in a sample of objects would be needed to confirm the gravitational origin of the measured redshifts. Meanwhile, we present a tentative black hole mass scaling relationship based on the Fe III$\lambda\lambda$2039-2113 redshift useful to measure the black hole mass of one individual object from a single spectrum.

We study the BELs in a sample of 11 gravitationally lensed quasars with at least two epochs of observation to identify intrinsic variability and to disentangle it from microlensing. To improve our statistical significance we also include 15 systems with single-epoch spectra. MgII and CIII] emission lines are only weakly affected by microlensing, but CIV shows strong microlensing in some cases, even for regions of the line core, associated with small projected velocities. However, excluding the strongly microlensed cases, there is a strikingly good match between the red wings of the CIV and CIII] profiles. Analysis of these results supports the existence of two regions in the BLR, one that is insensitive to microlensing (of size $\gtrsim$ 50 light-days and kinematics not confined to a plane) and another that shows up only when it is magnified by microlensing (of size of a few light-days, comparable to the accretion disk). Both regions can contribute in different proportions to the emission lines of different species and, within each line profile, to different velocity bins, all of which complicates detailed studies of the BLR based on microlensing size estimates. The strength of the microlensing indicates that some spectral features that make up the pseudo-continuum, such as the shelf-like feature at {\lambda}1610 or several FeIII blends, may in part arise from an inner region of the accretion disk. In the case of FeII, microlensing is strong in some blends but not in others. This opens up interesting possibilities to study quasar accretion disk kinematics.

We present eight monitoring seasons of the four brightest images of the gravitational lens SDSS J1004+4112 observed between December 2003 and October 2010. Using measured time delays for the images A, B and C and the model predicted time delay for image D we have removed the intrinsic quasar variability, finding microlensing events of about 0.5 and 0.7 mag of amplitude in the images C and D. From the statistics of microlensing amplitudes in images A, C, and D, we have inferred the half-light radius (at {\lambda} rest = 2407 {\AA}) for the accretion disk using two different methods, $R_{1/2}=8.7^{+18.5}_{-5.5} \sqrt{M/0.3 M_\odot}$ (histograms product) and $R_{1/2} = 4.2^{+3.2}_{-2.2} \sqrt{M/0.3 M_\odot}$ light-days ($\chi^2$). The results are in agreement within uncertainties with the size predicted from the black hole mass in SDSS J1004+4112 using the thin disk theory.