Search

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., Comment: SETI meets time domain astronomy. Accepted into the Astronomical Journal

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., Comment: 11 pages, 10 figures, accepted to MNRAS

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., Comment: 13 pages, 5 figures, 2 tables. Accepted for publication in ApJ

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 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., Comment: 28 pages, 16 figures, 5 tables, Accepted to ApJ

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., Comment: 8 preprint pages, 4 figures, accepted by ApJ Letters

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.