Precision astrometry from the second Gaia data release has allowed astronomers to predict 5,787 microlensing events, with 528 of these having maximums within the extended Gaia mission (J2014.5 - J2026.5). Future analysis of the Gaia time-series astrometry of these events will, in some cases, lead to precise gravitational mass measurements of the lens. We find that 61% of events predicted during the extended Gaia mission with sources brighter than G = 18 are likely to be spurious, with the background source in these cases commonly being either a duplicate detection or a binary companion of the lens. We present quality cuts to identify these spurious events and a revised list of microlensing event candidates. Our findings imply that half of the predictable astrometric microlensing events during the Gaia mission have yet to be identified.

The second data release of the Gaia mission contained astrometry and photometry for an incredible 1,692,919,135 sources, but how many sources did Gaia miss and where do they lie on the sky? The answer to this question will be crucial for any astronomer attempting to map the Milky Way with Gaia DR2. We infer the completeness of Gaia DR2 by exploiting the fact that it only contains sources with at least five astrometric detections. The odds that a source achieves those five detections depends on both the number of observations and the probability that an observation of that source results in a detection. We predict the number of times that each source was observed by Gaia and assume that the probability of detection is either a function of magnitude or a distribution as a function of magnitude. We fit both these models to the 1.7 billion stars of Gaia DR2, and thus are able to robustly predict the completeness of Gaia across the sky as a function of magnitude. We extend our selection function to account for crowding in dense regions of the sky, and show that this is vitally important, particularly in the Galactic bulge and the Large and Small Magellanic Clouds. We find that the magnitude limit at which Gaia is still 99% complete varies over the sky from $G=18.9$ to $21.3$. We have created a new Python package selectionfunctions (https://github.com/gaiaverse/selectionfunctions) which provides easy access to our selection functions.

The Gaia space mission is crafting revolutionary astrometric, photometric and spectroscopic catalogues that will allow us to map our Galaxy, but only if we know the completeness of this Gaia-verse of catalogues: what stars does it contain and what stars is it missing? We argue that the completeness is driven by Gaia's spinning-and-precessing scanning law and will apply this principle to the Gaia-verse over this series. We take a first step by identifying the periods in time that did not contribute any measurements to Gaia DR2; these gaps create ribbons of incompleteness across the sky that will bias any study that ignores them, although some of these gaps may be filled in future data releases. Our first approach was to use the variable star photometry to identify the 94 gaps longer than 1% of a day. Our second approach was to predict the number of observations of every point on the sky, which in comparison to the reported number of detections revealed additional gaps in the astrometry and spectroscopy. Making these predictions required us to make the most precise, publicly-available determination of the Gaia scanning law. Using this scanning law, we further identified that most stars fainter than $G=22$ in DR2 have spurious magnitudes due to a miscalibration resulting from a thunderstorm over Madrid. Our list of gaps and precision scanning law will allow astronomers to know when Gaia's eye was truly on their binary star, exoplanet or microlensing event during the time period of the second data release.

Most binaries are undetected. Astrometric reductions of a system using the assumption that the object moves like a single point mass can be biased by unresolved binary stars. The discrepancy between the centre of mass of the system (which moves like a point mass) and the centre of light (which is what we observe) introduces additional motion. We explore the extent to which binary systems affect single object models fit to astrometric data. This tells us how observations are diluted by binaries and which systems cause the largest discrepancies - but also allows us to make inferences about the binarity of populations based on observed astrometric error. By examining a sample of mock observations, we show that binaries with periods close to one year can mimic parallax and thus bias distance measurements, whilst long period binaries can introduce significant apparent proper motion. Whilst these changes can soak up some of the error introduced by the binary, the total deviation from the best fitting model can be translated into a lower limit on the on-sky separation of the pair. Throughout we link these predictions to data from the Gaia satellite, whilst leaving the conclusions generalizable to other surveys.

For stars with unresolved companions, motions of the centre of light and that of mass decouple, causing a single-source astrometric model to perform poorly. We show that such stars can be easily detected with the reduced chi2 statistic, or RUWE, provided as part of Gaia DR2. We convert RUWE into the amplitude of the image centroid wobble, which, if scaled by the source distance, is proportional to the physical separation between companions (for periods up to several years). We test this idea on a sample of known spectroscopic binaries and demonstrate that the amplitude of the centroid perturbation scales with the binary period and the mass ratio as expected. We apply this technique to the Gaia DR2 data and show how the binary fraction evolves across the Hertzsprung--Russell diagram. The observed incidence of unresolved companions is high for massive young stars and drops steadily with stellar mass, reaching its lowest levels for white dwarfs. We highlight the elevated binary fraction for the nearby Blue Stragglers and Blue Horizontal Branch stars. We also illustrate how unresolved hierarchical triples inflate the relative velocity signal in wide binaries. Finally, we point out a hint of evidence for the existence of additional companions to the hosts of extrasolar hot jupiters.

Selection functions are vital for understanding the observational biases of spectroscopic surveys. With the wide variety of multi-object spectrographs currently in operation and becoming available soon, we require easily generalisable methods for determining the selection functions of these surveys. Previous work, however, has largely been focused on generating individual, tailored selection functions for every data release of each survey. Moreover, no methods for combining these selection functions to be used for joint catalogues have been developed. We have developed a Poisson likelihood estimation method for calculating selection functions in a Bayesian framework, which can be generalised to any multi-object spectrograph. We include a robust treatment of overlapping fields within a survey as well as selection functions for combined samples with overlapping footprints. We also provide a method for transforming the selection function that depends on the sky positions, colour, and apparent magnitude of a star to one that depends on the galactic location, metallicity, mass, and age of a star. This `intrinsic' selection function is invaluable for chemodynamical models of the Milky Way. We demonstrate that our method is successful at recreating synthetic spectroscopic samples selected from a mock galaxy catalogue.

The Gaia Radial Velocity Spectrometer (RVS) provides a sample of 7,224,631 stars with full six-dimensional phase space information. Bayesian distances of these stars are available from the catalogue of Sch\"onrich et al. (2019). We exploit this to map out the behaviour of the velocity ellipsoid within 5 kpc of the Sun. We find that the tilt of the disc-dominated RVS sample is accurately described by the relation $\alpha = (0.952 \pm 0.007)\arctan (|z|/R)$, where ($R,z$) are cylindrical polar coordinates. This corresponds to velocity ellipsoids close to spherical alignment (for which the normalising constant would be unity) and pointing towards the Galactic centre. Flattening of the tilt of the velocity ellipsoids is enhanced close to the plane and Galactic centre, whilst at high elevations far from the Galactic center the population is consistent with exact spherical alignment. Using the LAMOST catalogue cross-matched with Gaia DR2, we construct thin disc and halo samples of reasonable purity based on metallicity. We find that the tilt of thin disc stars straddles $\alpha = (0.909-1.038)\arctan (|z|/R)$, and of halo stars straddles $\alpha = (0.927-1.063)\arctan (|z|/R)$. We caution against the use of reciprocal parallax for distances in studies of the tilt, as this can lead to serious artefacts.