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.

In this paper we report the current status of a new research program. The primary goal of the "Vanishing & Appearing Sources during a Century of Observations" (VASCO) project is to search for vanishing and appearing sources using existing survey data to find examples of exceptional astrophysical transients. The implications of finding such objects extend from traditional astrophysics fields to the more exotic searches for evidence of technologically advanced civilizations. In this first paper we present new, deeper observations of the tentative candidate discovered by Villarroel et al. (2016). We then perform the first searches for vanishing objects throughout the sky by comparing 600 million objects from the US Naval Observatory Catalogue (USNO) B1.0 down to a limiting magnitude of $\sim 20 - 21$ with the recent Pan-STARRS Data Release-1 (DR1) with a limiting magnitude of $\sim$ 23.4. We find about 150,000 preliminary candidates that do not have any Pan-STARRS counterpart within a 30 arcsec radius. We show that these objects are redder and have larger proper motions than typical USNO objects. We visually examine the images for a subset of about 24,000 candidates, superseding the 2016 study with a sample ten times larger. We find about $\sim$ 100 point sources visible in only one epoch in the red band of the USNO which may be of interest in searches for strong M dwarf flares, high-redshift supernovae or other catagories of unidentified red transients.

We construct a sample of 10,680 wall galaxies and 3,064 void galaxies with $M_R \lesssim -20$ by cross referencing a void catalog from literature with Baryon Oscillation Spectroscopic Survey (BOSS) CMASS and WiggleZ galaxies, where the CMASS survey targets redder galaxies and the WiggleZ survey targets bluer galaxies. Comparing the density profiles of the red and blue galaxies as a function of the void radius, we find that the number ratio of red-to-blue galaxies increases with distances from the void centers, suggesting a deficit of luminous and normal red galaxies in voids. We find a mean (g -- r) magnitude color of 1.298 and 1.210 for the wall and void galaxies, respectively, when considering the combined red and blue sample, which is found to be a significant difference. However, when considering the blue and red samples separately, we find no significant color difference. We conclude that the constituents galaxies of each population, rather than intrinsic color difference, is the main driver in the apparent average color difference of galaxies in voids and walls, indicating a deficit of luminous and normal red galaxies in voids. Our analysis suggests that the primary environmental-dependence effect on galaxy evolution for normal and luminous galaxies between void and wall regions is manifested in the number of red galaxies, which depends on the environmental-dependent merger history. Using a semi-analytic simulation model, we can successfully reproduce the apparent color difference between the void and wall galaxies.

Quasar microlensing serves as a unique probe of discrete objects within galaxies and galaxy clusters. Recent advancement of the technique shows that it can constrain planet-scale objects beyond our native galaxy by studying their induced microlensing signatures, the energy shift of emission lines originated in the vicinity of the black hole of high redshift background quasars. We employ this technique to exert effective constraints on the planet-mass object distribution within two additional lens systems, Q J0158$-$4325 ($z_l = 0.317$) and SDSS J1004+4112 ($z_l = 0.68$) using Chandra observations of the two gravitationally-lensed quasars. The observed variations of the emission line peak energy can be explained as microlensing of the FeK$\alpha$ emission region induced by planet-mass microlenses. To corroborate this, we perform microlensing simulations to determine the probability of a caustic transiting the source region and compare this with the observed line shift rates. Our analysis yields constraints on the sub-stellar population, with masses ranging from Moon ($10^{-8} M_{\odot}$) to Jupiter ($10^{-3} M_{\odot}$) sized bodies, within these galaxy or cluster scale structures, with total mass fractions of $\sim 3\times10^{-4}$ and $\sim 1\times10^{-4}$ with respect to halo mass for Q J0158$-$4325 and SDSS J1004+4112, respectively. Our analysis suggests that unbound planet-mass objects are universal in galaxies, and we surmise the objects to be either free-floating planets or primordial black holes. We present the first-ever constraints on the sub-stellar mass distribution in the intra-cluster light of a galaxy cluster. Our analysis yields the most stringent limit for primordial black holes at the mass range.

In this work, we calculated the sizes of unresolved X-ray emission regions in three gravitationally lensed radio-loud quasars, B\,1422+231, MG\,J0414+0534 and Q\,0957+561, using a combination of imaging and spectral analysis on the X-ray data taken from the \textit{Chandra X-Ray Observatory}. We tentatively detected FeK$\alpha$ emission lines in MG\,J0414+0534 and Q\,0957+561 with over 95\% significance, whereas, we did not significantly detect FeK$\alpha$ emission in B\,1422+231. We constructed differential microlensing light curves from absorption corrected count rates. We subsequently performed a microlensing analysis on the X-ray microlensing light curves to measure the X-ray source sizes in soft (0.83--3.6 keV), hard (3.6--21.8 keV), and full (0.83--21.8 keV) bands, based on either Bayesian or maximum likelihood probabilities. For B\,1422+231, sizes from the two methods are consistent with each other, e.g. $R_X^{hard}/R_G = 6.17\pm5.48 \text{ (Bayesian), } 11.81\pm3.75 \text{ (maximum likelihood)}$, where $R_G=GM_{BH}/c^2)$. However, for MG\,J0414+0534 and Q\,0957+561, the two methods yield completely different results suggesting that more frequently sampled data with better signal-to-noise ratio are needed to measure the source size for these two objects. Comparing the acquired size values with the radio-quiet sample in the literature we found that our results are consistent with X-ray source size scaling approximately as $R_X \propto M_{BH}$ with the mass of the central supermassive black hole. Our results also indicate that radio-loud quasars tend to have larger unresolved X-ray emission sizes compared to the radio-quiet ones.

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.

We present an analysis of Chandra spectra of five gravitationally lensed active galactic nuclei. We confirm the previous detections of FeK$\alpha$ emission lines in most images of these objects with high significance. The line energies range from 5.8 to 6.8 keV with widths from unresolved to 0.6 keV, consistent with emission close to spinning black holes viewed at different inclination angles. We also confirm the positive offset from the Iwasawa-Taniguchi effect, the inverse correlation between the FeK$\alpha$ equivalent width and the X-ray luminosity in AGN, where our measured equivalent widths are larger in lensed quasars. We attribute this effect to microlensing, and perform a microlensing likelihood analysis to constrain the emission size of the relativistic reflection region and the spin of supermassive black holes, assuming that the X-ray corona and the reflection region, responsible for the iron emission line, both follow power-law emissivity profiles. The microlensing analysis yields strong constraints on the spin and emissivity index of the reflection component for Q2237+0305, with $a > 0.92$ and $n > 5.4$. For the remaining four targets, we jointly constrain the two parameters, yielding $a=0.8\pm0.16$ and an emissivity index of $n=4.0\pm 0.8$, suggesting that the relativistic X-ray reflection region is ultra-compact and very close to the innermost stable circular orbits of black holes, which are spinning at close to the maximal value. We successfully constrain the half light radius of the emission region to $< 2.4$ $r_g$ ($r_g = GM/c^2$) for Q2237+0305 and in the range 5.9--7.4 $r_g$ for the joint sample.

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 define a second-order moment of the observational differential microlensing curves that can be used to impose constraints on physical properties of lensed quasars. We show that this quantity is sensitive both to variations in the source size and the deflector mass. We formulize a methodology to recover the source size from the observational measurements when the mass spectrum is fixed. As a case study, we test it with a sample of four quadruple lenses, both in simulated scenarios and with real data from the $\textit{Chandra X-ray Observatory}$. In our simulations with a uniform stellar population the method works best to detect sources around $0.1$ Einstein radii, giving correct upper/lower limits for much smaller/bigger sizes without requiring a big leap in additional computational effort as compared to a single-epoch approach, yet taking advantage of multi-epoch information. We apply the method to a small sample of X-ray data from four objects assuming a range of star masses, and obtain a degeneracy relation between the source size and deflector mass. Combined with previous estimates for the size of the X-ray corona, the degeneracy relation suggests that X-ray microlensing is mainly induced by planetary mass objects.

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.

Previously, planets have been detected only in the Milky Way galaxy. Here, we show that quasar microlensing provides a means to probe extragalactic planets in the lens galaxy, by studying the microlensing properties of emission close to the event horizon of the supermassive black hole of the background quasar, using the current generation telescopes. We show that a population of unbound planets between stars with masses ranging from Moon to Jupiter masses is needed to explain the frequent Fek line energy shifts observed in the gravitationally lensed quasar RXJ1131-1231 at a lens redshift of $z=0.295$ or 3.8 billion light-years away. We constrain the planet mass fraction to be larger than 0.0001 of the halo mass, which is equivalent to 2,000 objects ranging from Moon to Jupiter mass per main sequence star.

We present an X-ray photometric analysis of six gravitationally lensed quasars, with observation campaigns spanning from 5 to 14 years, measuring the total (0.83 - 21.8 keV restframe), soft (0.83 - 3.6 keV), and hard (3.6 - 21.8 keV) band image flux ratios or each epoch. Using the ratios of the model-predicted macro-magnifications as baselines, we build differential microlensing light curves and obtain joint likelihood functions for the average X-ray emission region sizes. Our analysis yields a Probability Distribution Function for the average half-light radius of the X-Ray emission region in the sample that peaks slightly above 1 gravitational radius and with nearly indistinguishable 68% confidence (one-sided) upper limits of 17.8 and 18.9 gravitational radii for the soft and hard X-ray emitting regions, assuming a mean stellar mass of 0.3 solar masses. We see hints of energy dependent microlensing between the soft and hard bands in two of the objects. In a separate analysis on the root-mean-square of the microlensing variability, we find significant differences between the soft and hard bands but the sign of the difference is not consistent across the sample. This suggests the existence of some kind of spatial structure to the X-ray emission in an otherwise extremely compact source. We also discover a correlation between the RMS microlensing variability and the average microlensing amplitude.

We measure the differential microlensing of the UV Fe II and Fe III emission line blends between 14 quasar image pairs in 13 gravitational lenses. We find that the UV iron emission is strongly microlensed in 4 cases with amplitudes comparable to that of the continuum. Statistically modeling the magnifications we infer a typical size of r ~ 4*sqrt(M/Msun) light-days for the Fe line emitting regions which is comparable to the size of the region generating the UV continuum (3 ~ 7 light-days). This may indicate that a significant part of the UV Fe II and Fe III emission originates in the quasar accretion disk.

We measure the differential microlensing of the broad emission lines between 18 quasar image pairs in 16 gravitational lenses. We find that high ionization lines such as CIV are more strongly microlensed than low ionization lines, indicating that the high ionization line emission regions are more compact. If we statistically model the distribution of microlensing magnifications, we obtain estimates for the broad line region radius of 24 (-15/+22) and 55 (-35/+150) light-days (90% confidence) for the high and low ionization lines, respectively. When the sample is divided attending to quasar luminosity, we find that the line emission regions of more luminous quasars are larger, with a slope consistent with the expected scaling from photoionization models. Our estimates also agree well with the results from local reveberation mapping studies.

We estimate the fraction of mass that is composed of compact objects in gravitational lens galaxies. This study is based on microlensing measurements (obtained from the literature) of a sample of 29 quasar image pairs seen through 20 lens galaxies. We determine the baseline for no microlensing magnification between two images from the ratios of emission line fluxes. Relative to this baseline, the ratio between the continua of the two images gives the difference in microlensing magnification. The histogram of observed microlensing events peaks close to no magnification and is concentrated below 0.6 magnitudes, although two events of high magnification, $\Delta m \sim 1.5$, are also present. We study the likelihood of the microlensing measurements using frequency distributions obtained from simulated microlensing magnification maps for different values of the fraction of mass in compact objects, $\alpha$. The concentration of microlensing measurements close to $\Delta m \sim 0$ can be explained only by simulations corresponding to very low values of $\alpha$ (10% or less). A maximum likelihood test yields $\alpha=0.05_{-0.03}^{+0.09}$ (90% confidence interval) for a quasar continuum source of intrinsic size $r_{s_0}\sim 2.6 \cdot 10^{15} \rm cm$. Regarding the current controversy about Milky Way/LMC and M31 microlensing studies, our work supports the hypothesis of a very low content in MACHOS (Massive Compact Halo Objects).