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.

We study the stellar kinematic properties and spatial distribution of the association Cyg OB3 using precise astrometric data from Gaia DR2. All known O- and B-type stars in Cyg OB3 region with positions, parallaxes and proper motions available are included, comprising a total of 41 stars. The majority of stars are found to be concentrated at a heliocentric distance of 2.0 +/- 0.3 kpc. The mean peculiar velocity of the sample after removing Galactic rotation and solar motion is ~22 km/s, dominated by the velocity component towards the Galactic center. The relative position and velocity of the black hole X-ray binary Cyg X-1 with respect to the association suggest that Cyg OB3 is most likely its parent association. The peculiar kinematic properties of some of the stars are revealed and are suggestive of past stellar encounters. The sample includes a previously known runaway star HD 227018, and its high peculiar velocity of ~50 km/s is confirmed with Gaia. We estimated the velocities of stars relative to the association and the star HD 225577 exhibits peculiar velocity smaller than its velocity relative to the association. The star has lower value of proper motion than the rest of the sample. The results suggest a slowly expanding nature of the association, which is supported by the small relative speeds <20 km/s with respect to the association for a majority of the sample stars.

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.

Stars slingshotted by the supermassive black hole at the Galactic centre will escape the Milky Way so quickly that their trajectories will be almost straight lines. Previous works have shown how these `hypervelocity stars' are subsequently deflected by the gravitational field of the Milky Way and the Large Magellanic Cloud (LMC), but have neglected to account for the reflex motion of the Milky Way in response to the fly by of the LMC. A consequence of this motion is that the hypervelocity stars we see on the outskirts of the Milky Way today were ejected from where the Milky Way centre was hundreds of millions of years ago. This change in perspective causes large apparent deflections in the trajectories of the hypervelocity stars, which are of the same order as the deflections caused by the gravitational force of the Milky Way and LMC. We quantify these deflections by simulating the production of hypervelocity stars in an isolated Milky Way (with a spherical or flattened dark matter halo), in a fixed-in-place Milky Way with a passing LMC, and in a Milky Way which responds to the passage of the LMC. The proper motion precision necessary to measure these deflections will be possible with the combination of Gaia with the proposed GaiaNIR successor mission, and these measurements will unlock the hypervelocity stars as probes of the shape of the Milky Way, the mass of the LMC, and of the dance of these two galaxies.

We present the serendipitous discovery of the fastest Main Sequence hyper-velocity star (HVS) by the Southern Stellar Stream Spectroscopic Survey (S5). The star S5-HVS1 is a $\sim 2.35$ M$_\odot$ A-type star located at a distance of $\sim 9$ kpc from the Sun and has a heliocentric radial velocity of $1017\pm 2.7$ km/s without any signature of velocity variability. The current 3-D velocity of the star in the Galactic frame is $1755\pm50$ km/s. When integrated backwards in time, the orbit of the star points unambiguously to the Galactic Centre, implying that S5-HVS1 was kicked away from Sgr A* with a velocity of $\sim 1800$ km/s and travelled for $4.8$ Myr to its current location. This is so far the only HVS confidently associated with the Galactic Centre. S5-HVS1 is also the first hyper-velocity star to provide constraints on the geometry and kinematics of the Galaxy, such as the Solar motion $V_{y,\odot}= 246.1\pm 5.3$ km/s or position $R_0=8.12\pm 0.23$ kpc. The ejection trajectory and transit time of S5-HVS1 coincide with the orbital plane and age of the annular disk of young stars at the Galactic centre, and thus may be linked to its formation. With the S5-HVS1 ejection velocity being almost twice the velocity of other hyper-velocity stars previously associated with the Galactic Centre, we question whether they have been generated by the same mechanism or whether the ejection velocity distribution has been constant over time.

We report spectroscopic measurements of stars in the recently discovered young stellar association Price-Whelan 1 (PW 1), which was found in the vicinity of the Leading Arm (LA) of the Magellanic Stream. We obtained Magellan+MIKE high-resolution spectra of the 28 brightest stars in PW 1 and used The Cannon to determine their stellar parameters. We find that the mean metallicity of PW 1 is [Fe/H]=-1.23 with a small scatter of 0.06 dex and the mean radial velocity is Vhelio=276.7 km/s with a dispersion of 11.0 km/s. Our results are consistent in Teff, logg, and [Fe/H] with the young and metal-poor characteristics (116 Myr and [Fe/H]=-1.1) determined for PW 1 from our discovery paper. We find a strong correlation between the spatial pattern of the PW 1 stars and the LA II gas with an offset of -10.15 deg in L_MS and +1.55 deg in B_MS. The similarity in metallicity, velocity, and spatial patterns indicates that PW 1 likely originated in LA II. We find that the spatial and kinematic separation between LA II and PW 1 can be explained by ram pressure from Milky Way gas. Using orbit integrations that account for the LMC and MW halo and outer disk gas, we constrain the halo gas density at the orbital pericenter of PW 1 to be n_halo (17 kpc) = 2.7 (3) x 10^-3 atoms/cm^3 and the disk gas density at the midplane at 20 kpc to be n_disk (20 kpc,0) = 6.0 (1.8) x 10^-2 atoms/cm^3. We, therefore, conclude that PW 1 formed from the LA II of the Magellanic Stream, making it a powerful constraint on the Milky Way-Magellanic interaction.

Gaia DR2 5932173855446728064 was recently proposed to be unbound from the Milky Way based on the $-614.3\pm2.5\;\mathrm{km}\;\mathrm{s}^{-1}$ median radial velocity given in Gaia DR2. We obtained eight epochs of spectroscopic follow-up and find a very different median radial velocity of $-56.5 \pm 5.3\;\mathrm{km}\;\mathrm{s}^{-1}$. If this difference were to be explained by binarity, then the unseen companion would be an intermediate-mass black hole; we therefore argue that the Gaia DR2 radial velocity must be in error. We find it likely that the spectra obtained by Gaia were dominated by the light from a star $4.3\;\mathrm{arcsec}$ away, and that, due to the slitless, time delay integration nature of Gaia spectroscopy, this angular offset corresponded to a spurious $620\;\mathrm{km}\;\mathrm{s}^{-1}$ shift in the calcium triplet of the second star. We argue that such unanticipated alignments between stars may account for 105 of the 202 stars with radial velocities faster than $500\;\mathrm{km}\;\mathrm{s}^{-1}$ in Gaia DR2 and propose a quality cut to exclude stars that are susceptible. We propose further cuts to remove stars where the colour photometry is suspect and stars where the radial velocity measurement is based on fewer than four transits, and thus produce an unprecedentedly clean selection of Gaia RVS stars for use in studies of Galactic dynamics.

We use Gaia Data Release 2 to search for possible surviving binary companions to three of the best studied historical Milky Way core-collapse supernovae. Consistent with previous work, we find there to be no plausible binary companion to either the Crab or Cas A supernovae. For the first time, we present a systematic search for a former companion to the Vela supernova, and rule out essentially any surviving luminous ($>L_\odot$) companion. Based on parallax and proper motion, we identify a faint source (Star A; Gaia Source ID 5521955992667891584) which is kinematically consistent with being a former binary companion to the Vela SN progenitor. However, the inferred absolute magnitude of this source is extremely faint, raising the possibility that it may in fact be a background interloper. In addition, we derive a new distance ($3.37^{+4.04}_{-0.97}$ kpc) to the Crab SN based on the Gaia parallax measurements, which is significantly further than the 2 kpc distance typically adopted. Finally, we demonstrate that Gaia can be used to measure the secular decline in the luminosity of the Crab pulsar, and provide a new test of pulsar models.

Iorio et al. (2018) mapped out the Milky Way halo using a sample of RR Lyrae stars drawn from a cross-match of Gaia with 2MASS. We investigate the significant residual in their model which we constrain to lie at Galactocentric radii $12<R<27\;\mathrm{kpc}$ and extend over $2600\;\mathrm{deg}^2$ of the sky. A counterpart of this structure exists in both the Catalina Real Time Survey and the sample of RR Lyrae variables identified in Pan-STARRS by Hernitschek et al. (2016), demonstrating that this structure is not caused by the spatial inhomogeneity of Gaia. The structure is likely the Virgo Stellar Stream and/or Virgo Over-Density. We show the structure is aligned with the Magellanic Stream and suggest that it is either debris from a disrupted dwarf galaxy that was a member of the Vast Polar Structure or that it is SMC debris from a tidal interaction of the SMC and LMC $3\;\mathrm{Gyr}$ ago. If the latter then the sub-structure in Virgo may have a Magellanic origin.

Double detonations in double white dwarf (WD) binaries undergoing unstable mass transfer have emerged in recent years as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. One potential outcome of this "dynamically driven double-degenerate double-detonation" (D^6) scenario is that the companion WD survives the explosion and is flung away with a velocity equal to its > 1000 km/s pre-SN orbital velocity. We perform a search for these hypervelocity runaway WDs using Gaia's second data release. In this paper, we discuss seven candidates followed up with ground-based instruments. Three sources are likely to be some of the fastest known stars in the Milky Way, with total Galactocentric velocities between 1000 and 3000 km/s, and are consistent with having previously been companion WDs in pre-SN Ia systems. However, although the radial velocity of one of the stars is > 1000 km/s, the radial velocities of the other two stars are puzzlingly consistent with 0. The combined five-parameter astrometric solutions from Gaia and radial velocities from follow-up spectra yield tentative 6D confirmation of the D^6 scenario. The past position of one of these stars places it within a faint, old SN remnant, further strengthening the interpretation of these candidates as hypervelocity runaways from binary systems that underwent SNe Ia.

Hypervelocity stars are intriguing rare objects traveling at speeds large enough to be unbound from the Milky Way. Several mechanisms have been proposed for producing them, including the interaction of the Galaxy's super-massive black hole (SMBH) with a binary; rapid mass-loss from a companion to a star in a short-period binary; the tidal disruption of an infalling galaxy and finally ejection from the Large Magellanic Cloud. While previously discovered high-velocity early-type stars are thought to be the result of an interaction with the SMBH, the origin of high-velocity late type stars is ambiguous. The second data release of Gaia (DR2) enables a unique opportunity to resolve this ambiguity and determine whether any late-type candidates are truly unbound from the Milky Way. In this paper, we utilize the new proper motion and velocity information available from DR2 to re-evaluate a collection of historical data compiled on the newly-created Open Fast Stars Catalog. We find that almost all previously-known high-velocity late-type stars are most likely bound to the Milky Way. Only one late-type object (LAMOST J115209.12+120258.0) is unbound from the Galaxy. Performing integrations of orbital histories, we find that this object cannot have been ejected from the Galactic centre and thus may be either debris from the disruption of a satellite galaxy or a disc runaway.

Using proper motion measurements from Gaia DR2, we probe the origin of 26 previously known hypervelocity stars (HVSs) around the Milky Way. We find that a significant fraction of these stars have a high probability of originating close to the Milky Way centre, but there is one obvious outlier. HVS3 is highly likely to be coming almost from the centre of the Large Magellanic Cloud (LMC). During its closest approach, $21.1^{+6.1}_{-4.6}$ Myr ago, it had a relative velocity of $870^{+69}_{-66}$ kms$^{-1}$ with respect to the LMC. This large kick velocity is only consistent with the Hills mechanism, requiring a massive black hole at the centre of the LMC. This provides strong direct evidence that the LMC itself harbours a massive black hole of at least $4\times 10^3 -10^4 M_\odot$.

We explore the hypothesis that B type emission-line stars (Be stars) have their origin in mass-transfer binaries by measuring the fraction of runaway Be stars. We assemble the largest-to-date catalogue of 632 Be stars with 6D kinematics, exploiting the precise astrometry of the Tycho-Gaia Astrometric Solution (TGAS) from the first Gaia Data Release. Using binary stellar evolution simulations, we make predictions for the runaway and equatorial rotation velocities of a runaway Be star population. Accounting for observational biases, we calculate that if all classical Be stars originated through mass transfer in binaries, then $17.5\%$ of the Be stars in our catalogue should be runaways. The remaining $82.5\%$ should be in binaries with subdwarfs, white dwarfs or neutron stars, because those systems either remained bound post-supernova or avoided the supernova entirely. Using a Bayesian methodology, we compare the hypothesis that each Be star in our catalogue is a runaway to the null hypothesis that it is a member of the Milky Way disc. We find that $13.1^{+2.6}_{-2.4}\%$ of the Be stars in our catalogue are runaways, and identify a subset of 40 high-probability runaways. We argue that deficiencies in our understanding of binary stellar evolution, as well as the degeneracy between velocity dispersion and number of runaway stars, can explain the slightly lower runaway fraction. We thus conclude that all Be stars could be explained by an origin in mass-transfer binaries. This conclusion is testable with the second Gaia data release (DR2).

We search for runaway former companions of the progenitors of nearby Galactic core-collapse supernova remnants (SNRs) in the Tycho-Gaia astrometric solution (TGAS). We look for candidates for a sample of ten SNRs with distances less than $2\;\mathrm{kpc}$, taking astrometry and $G$ magnitude from TGAS and $B,V$ magnitudes from the AAVSO Photometric All-Sky Survey (APASS). A simple method of tracking back stars and finding the closest point to the SNR centre is shown to have several failings when ranking candidates. In particular, it neglects our expectation that massive stars preferentially have massive companions. We evolve a grid of binary stars to exploit these covariances in the distribution of runaway star properties in colour - magnitude - ejection velocity space. We construct an analytic model which predicts the properties of a runaway star, in which the model parameters are the properties of the progenitor binary and the properties of the SNR. Using nested sampling we calculate the Bayesian evidence for each candidate to be the runaway and simultaneously constrain the properties of that runaway and of the SNR itself. We identify four likely runaway companions of the Cygnus Loop, HB 21, S147 and the Monoceros Loop. HD 37424 has previously been suggested as the companion of S147, however the other three stars are new candidates. The favoured companion of HB 21 is the Be star BD+50 3188 whose emission-line features could be explained by pre-supernova mass transfer from the primary. There is a small probability that the $2\;\mathrm{M}_{\odot}$ candidate runaway TYC 2688-1556-1 associated with the Cygnus Loop is a hypervelocity star. If the Monoceros Loop is related to the on-going star formation in the Mon OB2 association, the progenitor of the Monoceros Loop is required to be more massive than $40\;\mathrm{M}_{\odot}$ which is in tension with the posterior for HD 261393.

We explore the possibility that the observed population of Galactic hypervelocity stars (HVSs) originate as runaway stars from the Large Magellanic Cloud (LMC). Pairing a binary evolution code with an N-body simulation of the interaction of the LMC with the Milky Way, we predict the spatial distribution and kinematics of an LMC runaway population. We find that runaway stars from the LMC can contribute Galactic HVSs at a rate of $3 \times 10^{-6}\;\mathrm{yr}^{-1}$. This is composed of stars at different points of stellar evolution, ranging from the main-sequence to those at the tip of the asymptotic giant branch. We find that the known B-type HVSs have kinematics which are consistent with an LMC origin. There is an additional population of hypervelocity white dwarfs whose progenitors were massive runaway stars. Runaways which are even more massive will themselves go supernova, producing a remnant whose velocity will be modulated by a supernova kick. This latter scenario has some exotic consequences, such as pulsars and supernovae far from star-forming regions, and a small rate of microlensing from compact sources around the halo of the LMC.

There is growing theoretical and observational evidence that protoplanetary disc evolution may be significantly affected by the canonical levels of far ultraviolet (FUV) radiation found in a star forming environment, leading to substantial stripping of material from the disc outer edge even in the absence of nearby massive stars. In this paper we perform the first full radiation hydrodynamic simulations of the flow from the outer rim of protoplanetary discs externally irradiated by such intermediate strength FUV fields, including direct modelling of the photon dominated region (PDR) which is required to accurately compute the thermal properties. We find excellent agreement between our models and the semi-analytic models of Facchini et al. (2016) for the profile of the flow itself, as well as the mass loss rate and location of their "critical radius". This both validates their results (which differed significantly from prior semi-analytic estimates) and our new numerical method, the latter of which can now be applied to elements of the problem that the semi--analytic approaches are incapable of modelling. We also obtain the composition of the flow, but given the simple geometry of our models we can only hint at some diagnostics for future observations of externally irradiated discs at this stage. We also discuss the potential for these models as benchmarks for future photochemical-dynamical codes.

We predict the distribution of hypervelocity stars (HVSs) ejected from the Large Magellanic Cloud (LMC), under the assumption that the dwarf galaxy hosts a central massive black hole (MBH). For the majority of stars ejected from the LMC the orbital velocity of the LMC has contributed a significant fraction of their galactic rest frame velocity, leading to a dipole density distribution on the sky. We quantify the dipole using spherical harmonic analysis and contrast with the monopole expected for HVSs ejected from the Galactic Center. There is a tendril in the density distribution that leads the LMC which is coincident with the well-known and unexplained clustering of HVSs in the constellations of Leo and Sextans. Our model is falsifiable, since it predicts that Gaia will reveal a large density of HVSs in the southern hemisphere.