The discovery of multiply-imaged gravitationally lensed QSOs is fundamental to many astronomical and cosmological studies. However, these objects are rare and challenging to discover due to requirements of high-angular resolution astrometric, multiwavelength photometric and spectroscopic data. This has limited the number of known systems to a few hundred objects. We aim to reduce the constraints on angular resolution and discover multiply-imaged QSO candidates by using new candidate selection principles based on unresolved photometric time-series and ground-based images from public surveys. We selected candidates for multiply-imaged QSOs based on low levels of entropy computed from Catalina unresolved photometric time-series or Euclidean similarity to known lenses in a space defined by the wavelet power spectra of Pan-STARSS DR2 or DECaLS DR7 images, combined with multiple {\it Gaia} DR2 sources or large astrometric errors and supervised and unsupervised learning methods. We then confirmed spectroscopically some candidates with the Palomar Hale, Keck-I, and ESO/NTT telescopes. Here we report the discovery and confirmation of seven doubly-imaged QSOs and one likely double quasar. This demonstrates the potential of combining space-astrometry, even if unresolved, with low spatial-resolution photometric time-series and/or low-spatial resolution multi-band imaging to discover multiply-imaged lensed QSOs.

Aims: In this work, we aim to provide a reliable list of gravitational lens (GL) candidates based on a search performed over the entire Gaia Data Release 2 (Gaia DR2). We also show that the sole astrometric and photometric informations coming from the Gaia satellite yield sufficient insights for supervised learning methods to automatically identify GL candidates with an efficiency that is comparable to methods based on image processing. Methods: We simulated 106,623,188 lens systems composed of more than two images, based on a regular grid of parameters characterizing a non-singular isothermal ellipsoid lens model in the presence of an external shear. These simulations are used as an input for training and testing our supervised learning models consisting of Extremely Randomized Trees. The latter are finally used to assign to each of the 2,129,659 clusters of celestial objects a discriminant value that reflects the ability of our simulations to match the observed relative positions and fluxes from each cluster. Once complemented with additional constraints, these discriminant values allowed us to identify GL candidates out of the list of clusters. Results: We report the discovery of 15 new quadruply-imaged lens candidates with angular separations less than 6" and assess the performance of our approach by recovering 12 out of the 13 known quadruply-imaged systems with all their components detected in Gaia DR2 with a misclassification rate of fortuitous clusters of stars as lens systems that is below one percent. Similarly, the identification capability of our method regarding quadruply-imaged systems where three images are detected in Gaia DR2 is assessed by recovering 10 out of the 13 known quadruply-imaged systems having one of their constituting images discarded. The associated misclassification rate varying then between 5.8% and 20%, depending on the image we decided to remove.

We report the spectroscopic confirmation and modeling of the quadruply imaged quasar GRAL113100-441959, the first gravitational lens (GL) to be discovered mainly from astrometric considerations. Follow-up spectra obtained with Keck/LRIS reveal the lensing nature of this quadruply-imaged quasar with redshift $z_s = 1.090 \pm 0.002$, but show no evidence of the central lens galaxy. Using the image positions and $G$-band flux ratios provided by Gaia Data Release 2 as constraints, we model the system with a singular power-law elliptical mass distribution (SPEMD) plus external shear, to different levels of complexity. We show that relaxing the isothermal constraint of the SPEMD is not statistically significant, and thus we simplify the SPEMD to a singular isothermal ellipsoid to estimate the Einstein radius of the main lens galaxy $\theta_{\text{E}} = 0.\!\!^{\prime\prime}851$, the intensity and position angle of the external shear $(\gamma,\theta_{\gamma}) = (0.044,11.\!\!^{\circ}5)$, and we predict the lensing galaxy position to be $(x_{\text{gal}},y_{\text{gal}}) = (-0.\!\!^{\prime\prime}424,-0.\!\!^{\prime\prime}744)$ with respect to image A. We provide time delay predictions for pairs of images, assuming a plausible range of lens redshift values $z_l$ between $0.5$ and $0.9$. We finally examine the impact on time delays of the so-called Source Position Transformation, a family of degeneracies existing between different lens density profiles that reproduce most of the lensing observables equally well. We show that this effect contributes significantly to the time delay error budget and cannot be ignored during the modeling. This has implications for robust cosmography applications of lensed systems. GRAL113100-441959 is the first in a series of seven new spectroscopically confirmed GLs discovered from Gaia Data Release 2.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

Thanks to its spatial resolution the ESA/Gaia space mission offers a unique opportunity to discover new multiply-imaged quasars and to study the already known lensed systems at sub-milliarcsecond astrometric precisions. In this paper, we address the detection of the known multiply-imaged quasars from the Gaia Data Release 2 and determine the astrometric and photometric properties of the individually detected images found in the Gaia DR2 catalogue. We have compiled an exhaustive list of quasar gravitational lenses from the literature to search for counterparts in the Gaia Data Release 2. We then analyze the astrometric and photometric properties of these Gaia's detections. To highlight the tremendous potential of Gaia at the sub-milliarcsecond level we finally perform a simple Bayesian modeling of the well-known gravitational lens system HE0435-1223, using Gaia Data Release 2 and HST astrometry. From 478 known multiply imaged quasars, 200 have at least one image found in the Gaia Data Release 2. Among the 41 known quadruply-imaged quasars of the list, 26 have at least one image in the Gaia Data Release 2, 12 of which are fully detected, 6 have three counterparts, 7 have two and 1 has only one. As expected, the modeling of HE0435-1223 shows that the model parameters are significantly better constrained when using Gaia astrometry compared to HST astrometry, in particular the relative positions of the background quasar source and the centroid of the deflector. The Gaia sub-milliarcsecond astrometry also significantly reduces the parameter correlations. Besides providing an up-to-date list of multiply imaged quasars and their detection in the Gaia DR2, this paper shows that more complex modeling scenarios will certainly benefit from Gaia sub-milliarcsecond astrometry.

Context. Strong gravitationally lensed quasars are among the most interesting and useful observable extragalactic phenomena. Because their study constitutes a unique tool in various fields of astronomy, they are highly sought, not without difficulty. Indeed, even in this era of all-sky surveys, their recognition remains a great challenge, with barely a few hundred currently known systems. Aims. In this work we aim to detect new strongly lensed quasar candidates in the recently published Gaia Data Release 2 (DR2), which is the highest spatial resolution astrometric and photometric all-sky survey, attaining effective resolutions from 0.4" to 2.2". Methods. We cross-matched a merged list of quasars and candidates with the Gaia DR2 and found 1,839,143 counterparts within 0.5". We then searched matches with more than two Gaia DR2 counterparts within 6". We further narrowed the resulting list using astrometry and photometry compatibility criteria between the Gaia DR2 counterparts. A supervised machine learning method, Extremely Randomized Trees, is finally adopted to assign to each remaining system a probability of being lensed. Results. We report the discovery of three quadruply-imaged quasar candidates that are fully detected in Gaia DR2. These are the most promising new quasar lens candidates from Gaia DR2 and a simple singular isothermal ellipsoid lens model is able to reproduce their image positions to within $\sim$1 mas. This letter demonstrates the gravitational lens discovery potential of Gaia.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

Before the publication of the Gaia Catalogue, the contents of the first data release have undergone multiple dedicated validation tests. These tests aim at analysing in-depth the Catalogue content to detect anomalies, individual problems in specific objects or in overall statistical properties, either to filter them before the public release, or to describe the different caveats of the release for an optimal exploitation of the data. Dedicated methods using either Gaia internal data, external catalogues or models have been developed for the validation processes. They are testing normal stars as well as various populations like open or globular clusters, double stars, variable stars, quasars. Properties of coverage, accuracy and precision of the data are provided by the numerous tests presented here and jointly analysed to assess the data release content. This independent validation confirms the quality of the published data, Gaia DR1 being the most precise all-sky astrometric and photometric catalogue to-date. However, several limitations in terms of completeness, astrometric and photometric quality are identified and described. Figures describing the relevant properties of the release are shown and the testing activities carried out validating the user interfaces are also described. A particular emphasis is made on the statistical use of the data in scientific exploitation.

The Gaia mission started its regular observing program in the summer of 2014, and since then it is regularly obtaining observations of asteroids. This paper draws the outline of the data processing for Solar System objects, and in particular on the daily "short-term" processing, from the on-board data acquisition to the ground-based processing. We illustrate the tools developed to compute predictions of asteroid observations, we discuss the procedures implemented by the daily processing, and we illustrate some tests and validations of the processing of the asteroid observations. Our findings are overall consistent with the expectations concerning the performances of Gaia and the effectiveness of the developed software for data reduction.

The ESA Gaia mission, to be launched during 2013, will observe billions of objects, among which many galaxies, during its scanning of the sky. This will provide a large space-based dataset with unprecedented spatial resolution. Because of its natural Galactic and Astrometric priority, Gaia's observational strategy was optimized for point sources. Nonetheless, it is expected that 10^6 sources will be extragalactic, and a large portion of them will be angularly small galaxies. Although the mission was designed for point sources, a dedicated analysis of the raw data will allow the recovery of morphology of those objects at a 0.2" level. This may constitute a unique all-sky survey of such galaxies. We describe the conceptual design of the method we created for performing the morphological analysis of these objects as well as first results obtained from data simulations of low-resolution, highly binned, satellite data. Based on the obtained results we conclude that it is possible to push the limits of the Gaia space mission by analyzing galaxy morphology. (Abridged)