Search

We present new VLTI/GRAVITY near-infrared interferometric measurements of the angular size of the innermost hot dust continuum for 14 type 1 AGNs. The angular sizes are resolved on scales of ~0.7 mas and the inferred ring radii range from 0.028 to 1.33 pc, comparable to those reported previously and a factor 10-20 smaller than the mid-infrared sizes in the literature. Combining our new data with previously published values, we compile a sample of 25 AGN with bolometric luminosity ranging from $10^{42}$ to $10^{47} \rm erg~s^{-1}$, with which we study the radius-luminosity (R-L) relation for the hot dust structure. Our interferometric measurements of radius are offset by a factor 2 from the equivalent relation derived through reverberation mapping. Using a simple model to explore the dust structure's geometry, we conclude that this offset can be explained if the 2 um emitting surface has a concave shape. Our data show that the slope of the relation is in line with the canonical $R \propto L^{0.5}$ when using an appropriately non-linear correction for bolometric luminosity. In contrast, using optical luminosity or applying a constant bolometric correction to it results in a significant deviation in the slope, suggesting a potential luminosity dependence on the spectral energy distribution. Over four orders of magnitude in luminosity, the intrinsic scatter around the R-L relation is 0.2 dex, suggesting a tight correlation between innermost hot dust structure size and the AGN luminosity., Comment: 14 pages, 5 figures, 4 tables, accepted for publication in A&A

We present a multifrequency analysis of the radio galaxy 3CR 196.1 ($z = 0.198$), associated with the brightest galaxy of the cool core cluster CIZAJ0815.4-0303. This nearby radio galaxy shows a hybrid radio morphology and an X-ray cavity, all signatures of a turbulent past activity, potentially due to merger events and AGN outbursts. We present results of the comparison between $Chandra$ and VLT/MUSE data for the inner region of the galaxy cluster, on a scale of tens of kpc. We discovered H$\alpha$ + [N II]$\lambda6584$ emission spatially associated with the X-ray cavity (at $\sim$10 kpc from the galaxy nucleus) instead of with its rim. This result differs from previous discoveries of ionized gas surrounding X-ray cavities in other radio galaxies harbored in galaxy clusters and could represent the first reported case of ionized gas filling an X-ray cavity, either due to different AGN outbursts or to the cooling of warm ($10^4

We present EMIR, a powerful near-infrared (NIR) camera and multi-object spectrograph (MOS) installed at the Nasmyth focus of the 10.4 m GTC. EMIR was commissioned in mid-2016 and is offered as a common-user instrument. It provides spectral coverage of 0.9 to 2.5 $\mu m$ over a field of view (FOV) of 6.67x6.67 squared arcmin in imaging mode, and 6.67x4 squared arcmin in spectroscopy. EMIR delivers up to 53 spectra of different objects thanks to a robotic configurable cold slit mask system that is located inside the cryogenic chamber, allowing rapid reconfiguration of the observing mask. The imaging mode is attained by moving all bars outside the FOV and then leaving an empty space in the GTC focal surface. The dispersing suite holds three large pseudo-grisms, formed by the combination of high-efficiency FuSi ion-etched ruled transmission grating sandwiched between two identical ZnSe prisms, plus one standard replicated grism. These dispersing units offer the spectral recording of an atmospheric window $J,H,K$ in a single shot with resolving powers of 5000, 4250, 4000, respectively for a nominal slit width of 0.6\arcsec, plus the combined bands $YJ$ or $HK$, also in a single shot, with resolution of $\sim$ 1000. The original Hawaii2 FPA detector, which is prone to instabilities that add noise to the signal, is being replaced by a new Hawaii2RG detector array, and is currently being tested at the IAC. This paper presents the most salient features of the instrument, with emphasis on its observing capabilities and the functionality of the configurable slit unit. Sample early science data is also shown., Comment: 10 pages, 11 figures

This work focuses on active galactic nuclei (AGNs), and the relation between the sizes of the hot dust continuum and the broad-line region (BLR). We find that the continuum size measured using optical/near-infrared interferometry (OI) is roughly twice that measured by reverberation mapping (RM). Both OI and RM continuum sizes show a tight relation with the H$\beta$ BLR size with only an intrinsic scatter of 0.25 dex. The masses of supermassive black holes (BHs) can hence be simply derived from a dust size in combination with a broad line width and virial factor. Since the primary uncertainty of these BH masses comes from the virial factor, the accuracy of the continuum-based BH masses is close to those based on the RM measurement of the broad emission line. Moreover, the necessary continuum measurements can be obtained on a much shorter timescale than those required monitoring for RM, and are also more time efficient than those needed to resolve the BLR with OI. The primary goal of this work is to demonstrate measuring the BH mass based on the dust continuum size with our first calibration of the $R_\mathrm{BLR}$-$R_\mathrm{d}$ relation. The current limitation and caveats are discussed in detail. Future GRAVITY observations are expected to improve the continuum-based method and have the potential to measure BH masses for a large sample of AGNs in the low-redshift Universe., Comment: Accepted for publication in A&A; 11 pages, 5 figures, 3 tables

Using bona-fide black hole (BH) mass estimates from reverberation mapping and the line ratio [SiVI] 1.963$\mu$m/Br$\gamma_{\rm broad}$ as tracer of the AGN ionising continuum, a novel BH-mass scaling relation of the form log($M_{\rm BH}) = (6.40\pm 0.17) - (1.99\pm 0.37) \times$ log ([SiVI]/Br$\gamma_{\rm broad})$, dispersion 0.47 dex, over the BH mass interval, $10^6 - 10^8$ M$_{\odot}$ is found. Following on the geometrically thin accretion disc approximation and after surveying a basic parameter space for coronal lines production, we believe one of main drivers of the relation is the effective temperature of the disc, which is effectively sampled by the [SiVI] 1.963$\mu$m coronal line for the range of BH masses considered. By means of CLOUDY photoionisation models, the observed anti-correlation appears to be formally in line with the thin disc prediction T_disc $\propto {M_{\rm BH}}^{-1/4}$., Comment: 22 pages, 16 figures. Accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:2010.00075

We report observations of a relatively long period of 3He-rich solar energetic particles (SEPs) measured by Solar Orbiter. The period consists of several well-resolved ion injections. The high-resolution STEREO-A imaging observations reveal that the injections coincide with EUV jets/brightenings near the east limb, not far from the nominal magnetic connection of Solar Orbiter. The jets originated in two adjacent, large, and complex active regions as observed by the Solar Dynamics Observatory when the regions rotated to the Earth's view. It appears that the sustained ion injections were related to the complex configuration of the sunspot group and the long period of 3He-rich SEPs to the longitudinal extent covered by the group during the analyzed time period., Comment: accepted for publication in A&A, Letters to the Editor

The angular size of the broad line region (BLR) of the nearby active galactic nucleus (AGN) NGC 3783 has been spatially resolved by recent observations with VLTI/GRAVITY. A reverberation mapping (RM) campaign has also recently obtained high quality light curves and measured the linear size of the BLR in a way that is complementary to the GRAVITY measurement. The size and kinematics of the BLR can be better constrained by a joint analysis that combines both GRAVITY and RM data. This, in turn, allows us to obtain the mass of the supermassive black hole in NGC3783 with an accuracy that is about a factor of two better than that inferred from GRAVITY data alone. We derive $M_\mathrm{BH}=2.54_{-0.72}^{+0.90}\times 10^7\,M_\odot$. Finally, and perhaps most notably, we are able to measure a geometric distance to NGC 3783 of $39.9^{+14.5}_{-11.9}$ Mpc. We are able to test the robustness of the BLR-based geometric distance with measurements based on the Tully-Fisher relation and other indirect methods. We find the geometric distance is consistent with other methods within their scatter. We explore the potential of BLR-based geometric distances to directly constrain the Hubble constant, $H_0$, and identify differential phase uncertainties as the current dominant limitation to the $H_0$ measurement precision for individual sources., Comment: 9 pages and 5 figures in main text, Accepted for publication in A&A

The nuclear stellar disc (NSD) is, together with the nuclear star cluster (NSC) and the central massive black hole, one of the main components in the central parts of our Milky Way. However, until recently, only few studies of the stellar content of the NSD have been obtained due to extreme extinction and stellar crowding. With a dedicated KMOS (VLT, ESO) spectroscopic survey, we study the kinematics and global metallicities of the NSD based on the observations of K/M giant stars. We trace radial velocities and metallicities which were derived based on spectral indices (Na I and CO) along the NSD and compare those with a Galactic Bulge sample of APOGEE (DR16) and data from the NSC. We find that the metallicity distribution function and the fraction of metal-rich and metal-poor stars in the NSD are different from the corresponding distributions and ratios of the NSC and the Galactic Bulge. By tracing the velocity dispersion as a function of metallicity, we clearly see that the NSD is kinematically cool and that the velocity dispersion decreases with increasing metallicity contrary to the inner Bulge sample of APOGEE ($\rm |b| < 4^{o}$). Using molecular gas tracers ($\rm H_{2}CO$, CO(4-3)) of the Central Molecular Zone (CMZ) we find an astonishing agreement between the gas rotation and the rotation of the metal-rich population indicating that the metal-rich stars could have formed from gas in the CMZ. On the other hand, the metal-poor stars show a much slower rotation profile with signs of counter-rotation indicating a different origin of these stars. Coupling kinematics with global metallicities, our results demonstrate that the NSD is chemically and kinematically distinct with respect to the inner Bulge indicating a different formation scenario., Comment: accepted for publication in A\&A

We present the analysis of X-ray and optical observations of gas filaments observed in the radio source 3CR 318.1, associated with NGC 5920, the Brightest Cluster Galaxy (BCG) of MKW 3s, a nearby cool core galaxy cluster. This work is one of the first X-ray and optical analyses of filaments in cool core clusters carried out using MUSE observations. We aim at identifying the main excitation processes responsible for the emission arising from these filaments. We complemented the optical VLT/MUSE observations, tracing the colder gas phase, with X-ray $\textit{Chandra}$ observations of the hotter highly ionized gas phase. Using the MUSE observations, we studied the emission line intensity ratios along the filaments to constrain the physical processes driving the excitation, and, using the $\textit{Chandra}$ observations, we carried out a spectral analysis of the gas along these filaments. We found a spatial association between the X-ray and optical morphology of these filaments, which are colder and have lower metal abundance than the surrounding intra-cluster medium (ICM), as already seen in other BCGs. Comparing with previous results from the literature for other BCGs, we propose that the excitation process that is most likely responsible for these filaments emission is a combination of star formation and shocks, with a likely contribution from self-ionizing, cooling ICM. Additionally, we conclude that the filaments most likely originated from AGN-driven outflows in the direction of the radio jet., Comment: 9 pages, 5 figures, ApJL accepted, pre-proof version

Using VLTI/GRAVITY and SINFONI data, we investigate the sub-pc gas and dust structure around the nearby type 1 AGN hosted by NGC 3783. The K-band coverage of GRAVITY uniquely allows a simultaneous analysis of the size and kinematics of the broad line region (BLR), the size and structure of the near-IR continuum emitting hot dust, and the size of the coronal line region (CLR). We find the BLR probed through broad Br$\gamma$ emission is well described by a rotating, thick disk with a radial distribution of clouds peaking in the inner region. In our BLR model the physical mean radius of 16 light days is nearly twice the 10 day time lag that would be measured, which matches very well the 10 day time lag that has been measured by reverberation mapping. We measure a hot dust FWHM size of 0.74 mas (0.14 pc) and further reconstruct an image of the hot dust which reveals a faint (5% of the total flux) offset cloud which we interpret as an accreting cloud heated by the central AGN. Finally, we directly measure the FWHM size of the nuclear CLR as traced by the [CaVIII] and narrow Br$\gamma$ line. We find a FWHM size of 2.2 mas (0.4 pc), fully in line with the expectation of the CLR located between the BLR and narrow line region. Combining all of these measurements together with larger scale near-IR integral field unit and mid-IR interferometry data, we are able to comprehensively map the structure and dynamics of gas and dust from 0.01--100 pc., Comment: 20 pages and 18 figures in main text, Accepted for publication in A&A

In the central few degrees of the bulge of the Milky Way there is a flattened structure of gas, dust and stars (the central molecular zone) similar to nuclear disks in other galaxies. Due to extreme foreground extinction we possess only sparse information about the (mostly old) stellar population of the nuclear disc. Here we present our KMOS spectroscopic survey of the stars in the nuclear disk reaching the old populations. In order to obtain an unbiased data set, we sample stars in the full extinction range along each line-of-sight. We also observe reference fields in neighboring regions of the Galactic bulge. Here we describe the design and execution of the survey and present first results. We obtain spectra and five spectral indices of 3113 stars with a median S/N of 67 and measure radial velocities for 3051 stars. Of those, 2735 sources have sufficient S/N to estimate temperatures and metallicities from indices. We derive metallicities using the CO 2-0 and Na I K-band spectral features, where we derive our own empirical calibration using metallicities obtained with higher resolution observations. We use 183 giant stars for calibration spanning in metallicity from -2.5 to 0.6 dex and covering temperatures of up to 5500 K. The derived index based metallicities deviate from the calibration values with a scatter of 0.32 dex. The internal uncertainty of our metallicities is likely smaller. We use these metallicity measurements together with the CO index to derive effective temperatures using literature relations. We publish the catalog here. Our data set complements Galactic surveys such as Gaia and APOGEE for the inner 200 pc radius of the Milky Way which is not readily accessible by those surveys due to extinction. We will use the derived properties in future papers for further analysis of the nuclear disk., Comment: 19 pages, 14 figures, 3 tables, accepted for publication in A&A

We present a systematic analysis of the extended X-ray emission discovered around 35 FR II radio galaxies from the revised Third Cambridge catalog (3CR) Chandra Snapshot Survey with redshifts between 0.05 to 0.9. We aimed to (i) test for the presence of extended X-ray emission around FR II radio galaxies, (ii) investigate if the extended emission origin is due to Inverse Compton scattering of seed photons arising from the Cosmic Microwave Background (IC/CMB) or to thermal emission from an intracluster medium (ICM) and (iii) test the impact of this extended emission on hotspot detection. We investigated the nature of the extended X-ray emission by studying its morphology and compared our results with low-frequency radio observations (i.e., $\sim$150 MHz), in the TGSS and LOFAR archives, as well as with optical images from Pan-STARRS. In addition, we optimized a search for X-ray counterparts of hotspots in 3CR FR II radio galaxies. We found statistically significant extended emission ($>$3$\sigma$ confidence level) along the radio axis for $\sim$90%, and in the perpendicular direction for $\sim$60% of our sample. We confirmed the detection of 7 hotspots in the 0.5 - 3 keV. In the cases where the emission in the direction perpendicular to the radio axis is comparable to that along the radio axis, we suggest that the underlying radiative process is thermal emission from ICM. Otherwise, the dominant radiative process is likely non-thermal IC/CMB emission from lobes. We found that non-thermal IC/CMB is the dominant process in $\sim$70% of the sources in our sample, while thermal emission from the ICM dominates in $\sim$15% of them., Comment: 29 pages, 11 figures, ApJS accepted, pre-proof version

The Vanishing & Appearing Sources during a Century of Observations (VASCO) project investigates astronomical surveys spanning a time interval of 70 years, searching for unusual and exotic transients. We present herein the VASCO Citizen Science Project, which can identify unusual candidates driven by three different approaches: hypothesis, exploratory, and machine learning, which is particularly useful for SETI searches. To address the big data challenge, VASCO combines three methods: the Virtual Observatory, user-aided machine learning, and visual inspection through citizen science. Here we demonstrate the citizen science project and its improved candidate selection process, and we give a progress report. We also present the VASCO citizen science network led by amateur astronomy associations mainly located in Algeria, Cameroon, and Nigeria. At the moment of writing, the citizen science project has carefully examined 15,593 candidate image pairs in the data (ca. 10% of the candidates), and has so far identified 798 objects classified as "vanished". The most interesting candidates will be followed up with optical and infrared imaging, together with the observations by the most potent radio telescopes., Comment: 16 pages, 7 figures, contribution to IAC 2020 conference

We present new near-infrared VLTI/GRAVITY interferometric spectra that spatially resolve the broad Br$\gamma$ emission line in the nucleus of the active galaxy IRAS 09149-6206. We use these data to measure the size of the broad line region (BLR) and estimate the mass of the central black hole. Using an improved phase calibration method that reduces the differential phase uncertainty to 0.05 degree per baseline across the spectrum, we detect a differential phase signal that reaches a maximum of ~0.5 degree between the line and continuum. This represents an offset of ~120 $\mu$as (0.14 pc) between the BLR and the centroid of the hot dust distribution traced by the 2.3 $\mu$m continuum. The offset is well within the dust sublimation region, which matches the measured ~0.6 mas (0.7 pc) diameter of the continuum. A clear velocity gradient, almost perpendicular to the offset, is traced by the reconstructed photocentres of the spectral channels of the Br$\gamma$ line. We infer the radius of the BLR to be ~65 $\mu$as (0.075 pc), which is consistent with the radius-luminosity relation of nearby active galactic nuclei derived based on the time lag of the H$\beta$ line from reverberation mapping campaigns. Our dynamical modelling indicates the black hole mass is $\sim 1\times10^8\,M_\odot$, which is a little below, but consistent with, the standard $M_{\rm BH}$-$\sigma_*$ relation., Comment: Accepted for publication in A&A. 14 pages, 11 figures in main text

We present the analysis of nine radio sources belonging to the Third Cambridge Revised catalog (3CR) observed with $Chandra$ during Cycle 20 in the redshift range between 1.5 and 2.5. This study completes the 3CR $Chandra$ Snapshot Survey thus guaranteeing the X-ray coverage of all 3CR sources identified to date. This sample lists two compact steep spectrum sources, four radio galaxies and three quasars. We detected X-ray emission from all nuclei, with the only exception of 3C 326.1 and 3C 454.1 and from radio lobes in 6 out of 9 sources at level of confidence larger than $\sim$5$\sigma$. We measured X-ray fluxes and luminosities for all nuclei and lobes in the soft (0.5 - 1 keV), medium (1 - 2 keV) and hard (2 - 7 keV) X-ray bands. Since the discovered X-ray extended emission is spatially coincident with the radio structure in all cases, its origin could be due to Inverse Compton scattering of the Cosmic Microwave Background (IC/CMB) occurring in radio lobes., Comment: 15 pages, 11 figures, 3 tables, ApJS accepted, pre-proof version

We study the optical linear and circular polarisation in the optically-thin regime of the core and jet of M87. Observations were acquired two days before the Event Horizon Telescope (EHT) campaign in early April 2017. A high degree ($\sim 20$ per cent) of linear polarisation (P$_{\rm lin}$) is detected in the bright jet knots resolved at $\sim 10\, \rm{arcsec}$ to $23\, \rm{arcsec}$ ($0.8$-$1.8\, \rm{kpc}$) from the centre, whereas the nucleus and inner jet show P$_{\rm lin} \lesssim 5$ per cent. The position angle of the linear polarisation shifts by $\sim 90$ degrees from each knot to the adjacent ones, with the core angle perpendicular to the first knot. The nucleus was in a low level of activity (P$_{\rm lin} \sim 2$-$3$ per cent), and no emission was detected from HST-1. No circular polarisation was detected either in the nucleus or the jet above a $3\sigma$ level of P$_{\rm circ} \leq 1.5$ per cent, discarding the conversion of P$_{\rm lin}$ into P$_{\rm circ}$. A disordered magnetic field configuration or a mix of unresolved knots polarised along axes with different orientations could explain the low P$_{\rm lin}$. The latter implies a smaller size of the core knots, in line with current interferometric observations. Polarimetry with EHT can probe this scenario in the future. A steep increase of both P$_{\rm lin}$ and P$_{\rm circ}$ with increasing frequency is expected for the optically-thin domain, above the turnover point. This work describes the methodology to recover the four Stokes parameters using a $\lambda/4$ wave-plate polarimeter., Comment: Accepted for publication in MNRAS. 10 pages, 8 figures

A statistical study of intermediate Palomar Transient Factory supernovae (SNe) in Type 1 AGN has shown a major deficit of supernovae around Type 1 AGN host galaxies, with respect to Type 2 AGN hosts. The aim of this work is to test whether there is any preference for Type 1 AGN to host SN of a specific kind. Through the analysis of SN occurrence and their type (thermonuclear vs core-collapse), we can directly link the type of stars producing the SN events, thus this is an indirect way to study host galaxies in Type 1 AGN. We examine the detection fractions of SNe, the host galaxies and compare the sample properties to typical host galaxies in the Open Supernova Catalog (OSC; Guillochon et al. 2017). The majority of the host galaxies in the AGN sample are late-type, similar to typical galaxies hosting SN within the OSC. The findings are supportive of a deficiency of SNe near Type 1 AGN, although we cannot with certainty assess the overall detection fractions of SNe in Type 1 AGN relative to other SN host galaxies. We can state that Type 1 AGN has equal detection fractions of thermonuclear vs core-collapse SNe. However, we note the possibility of a higher detection rate of core-collapse supernovae in Type-1 AGN with insecure AGN classifications., Comment: 12 pages, 5 figures, 3 tables. Accepted for publication in MNRAS

We report results on deep Jansky Very Large Array A-configuration observations at 22 GHz of the hotspots of the radio galaxies 3C227 and 3C445. Synchrotron emission in the optical on scales up to a few kpc was reported for the four hotspots. Our VLA observations point out the presence of unresolved regions with upper limit to their linear size of about 100 pc. This is the first time that such compact components in hotspots have been detected in a mini-sample, indicating that they are not a peculiar characteristic of a few individual hotspots. The polarization may reach values up to 70 per cent in compact (about 0.1 kpc scale) regions within the hotspot, indicating a highly ordered magnetic field with size up to a hundred parsecs. On larger scales, the average polarization of the hotspot component is about 30-45 per cent, suggesting the presence of a significant random field component, rather than an ordered magnetic field. This is further supported by the displacement between the peaks in polarized intensity and in total intensity images that is observed in all the four hotspots. The electric vector position angle is not constant, but changes arbitrarily in the central part of the hotspot regions, whereas it is usually perpendicular to the total intensity contours of the outermost edge of the hotspot structure, likely marking the large-scale shock front. The misalignment between X-ray and radio-to-optical emission suggests that the former is tracing the current particle acceleration, whereas the latter marks older shock fronts., Comment: 12 pages, 13 figures, published in MNRAS

The acceleration and radiative processes active in low-power radio hotspots are investigated by means of new deep near-infrared (NIR) and optical VLT observations, complemented with archival, high-sensitivity VLT, radio VLA and X-ray Chandra data. For the three studied radio galaxies (3C 105, 3C 195 and 3C 227), we confirm the detection of NIR/optical counterparts of the observed radio hotspots. We resolve multiple components in 3C 227 West and in 3C 105 South and characterize the diffuse NIR/optical emission of the latter. We show that the linear size of this component ($\gtrsim$4 kpc) makes 3C 105 South a compelling case for particles' re-acceleration in the post-shock region. Modeling of the radio-to-X-ray spectral energy distribution (SED) of 3C 195 South and 3C 227 W1 gives clues on the origin of the detected X-ray emission. In the context of inverse Compton models, the peculiarly steep synchrotron curve of 3C 195 South sets constraints on the shape of the radiating particles' spectrum that are testable with better knowledge of the SED shape at low ($\lesssim$GHz) radio frequencies and in X-rays. The X-ray emission of 3C 227 W1 can be explained with an additional synchrotron component originating in compact ($<$100 pc) regions, such those revealed by radio observations at 22 GHz, provided that efficient particle acceleration ($\gamma\gtrsim$10$^7$) is ongoing. The emerging picture is that of systems in which different acceleration and radiative processes coexist., Comment: 16 pages, 10 figures, accepted for publication in MNRAS

We present near-infrared interferometric data on the Seyfert 2 galaxy NGC 1068, obtained with the GRAVITY instrument on the European Southern Observatory Very Large Telescope Interferometer. The extensive baseline coverage from 5 to 60 M\lambda allowed us to reconstruct a continuum image of the nucleus with an unrivaled 0.2 pc resolution in the K-band. We find a thin ring-like structure of emission with a radius r = 0.24+/-0.03 pc, inclination i = 70+/-5 deg, position angle PA = -50+/-4 deg, and h/r < 0.14, which we associate with the dust sublimation region. The observed morphology is inconsistent with the expected signatures of a geometrically and optically thick torus. Instead, the infrared emission shows a striking resemblance to the 22 GHz maser disc, which suggests they share a common region of origin. The near-infrared spectral energy distribution indicates a bolometric luminosity of (0.4-4.7) x 10^45 erg/s, behind a large A_K ~ 5.5 (A_V ~ 90) screen of extinction that also appears to contribute significantly to obscuring the broad line region., Comment: Accepted for publication in Astronomy and Astrophysics. 13 pages with 7 figures

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