Hypervelocity stars (HVSs) travel from the Galactic Centre across the dark matter halo of the Milky Way, where they are observed with velocities in excess of the Galactic escape speed. Because of their quasi-radial trajectories, they represent a unique probe of the still poorly constrained dark matter component of the Galactic potential. In this paper, we present a new method to produce such constraints. Our likelihood is based on the local HVS density obtained by back-propagating the observed phase space position and quantifies the ejection probability along the orbit. To showcase our method, we apply it to simulated Gaia samples of $\sim200$ stars in three realistic Galactic potentials with dark matter components parametrized by spheroidal NFW profiles. We find that individual HVSs exhibit a degeneracy in the scale mass-scale radius plane ($M_s-r_s$) and are able to measure only the combination $\alpha = M_s/r_s^2$. Likewise, a degeneracy is also present between $\alpha$ and the spheroidal axis-ratio $q$. In the absence of observational errors, we show the whole sample can nail down both parameters with {\it sub-per cent} precision (about $1\%$ and $0.1\%$ for $\alpha$ and $q$ respectively) with no systematic bias. This remarkable power to constrain deviations from a symmetric halo is a consequence of the Galactocentric origin of HVSs. To compare our results with other probes, we break the degeneracy in the scale parameters and impose a mass-concentration relation. The result is a competitive precision on the virial mass $M_{200}$ of about $10\%$.

We search for the fastest stars in the subset of stars with radial velocity measurements of the second data release (DR2) of the European Space Agency mission Gaia. Starting from the observed positions, parallaxes, proper motions, and radial velocities, we construct the distance and total velocity distribution of more than $7$ million stars in our Milky Way, deriving the full 6D phase space information in Galactocentric coordinates. These information are shared in a catalogue, publicly available at http://home.strw.leidenuniv.nl/~marchetti/research.html. To search for unbound stars, we then focus on stars with a probability greater than $50 \%$ of being unbound from the Milky Way. This cut results in a clean sample of $125$ sources with reliable astrometric parameters and radial velocities. Of these, $20$ stars have probabilities greater than 80 $\%$ of being unbound from the Galaxy. On this latter sub-sample, we perform orbit integration to characterize the stars' orbital parameter distributions. As expected given the relatively small sample size of bright stars, we find no hypervelocity star candidates, stars that are moving on orbits consistent with coming from the Galactic Centre. Instead, we find $7$ hyper-runaway star candidates, coming from the Galactic disk. Surprisingly, the remaining $13$ unbound stars cannot be traced back to the Galaxy, including two of the fastest stars (around $700$ km/s). If conformed, these may constitute the tip of the iceberg of a large extragalactic population or the extreme velocity tail of stellar streams.

Hypervelocity stars (HVSs) are amongst the fastest objects in our Milky Way. These stars are predicted to come from the Galactic center (GC) and travel along unbound orbits across the Galaxy. In the coming years, the ESA satellite Gaia will provide the most complete and accurate catalogue of the Milky Way, with full astrometric parameters for more than $1$ billion stars. In this paper, we present the expected sample size and properties (mass, magnitude, spatial, velocity distributions) of HVSs in the Gaia stellar catalogue. We build three Gaia mock catalogues of HVSs anchored to current observations, exploring different ejection mechanisms and GC stellar population properties. In all cases, we predict hundreds to thousands of HVSs with precise proper motion measurements within a few tens of kpc from us. For stars with a relative error in total proper motion below $10 \%$, the mass range extends to ~$10 M_{\odot}$ but peaks at ~$1$ $M_\odot$. The majority of Gaia HVSs will therefore probe a different mass and distance range compared to the current non-Gaia sample. In addition, a subset of a few hundreds to a few thousands of HVSs with $M$ ~ $3$ $M_\odot$ will be bright enough to have a precise measurement of the three-dimensional velocity from Gaia alone. Finally, we show that Gaia will provide more precise proper motion measurements for the current sample of HVS candidates. This will help identifying their birthplace narrowing down their ejection location, and confirming or rejecting their nature as HVSs. Overall, our forecasts are extremely encouraging in terms of quantity and quality of HVS data that can be exploited to constrain both the Milky Way potential and the GC properties.

The paucity of hypervelocity stars (HVSs) known to date has severely hampered their potential to investigate the stellar population of the Galactic Centre and the Galactic Potential. The first Gaia data release gives an opportunity to increase the current sample. The challenge is the disparity between the expected number of hypervelocity stars and that of bound background stars. We have applied a novel data mining algorithm based on machine learning techniques, an artificial neural network, to the Tycho-Gaia astrometric solution (TGAS) catalogue. With no pre-selection of data, we could exclude immediately $\sim 99 \%$ of the stars in the catalogue and find 80 candidates with more than $90\%$ predicted probability to be HVSs, based only on their position, proper motions, and parallax. We have cross-checked our findings with other spectroscopic surveys, determining radial velocities for 30 and spectroscopic distances for 5 candidates. In addition, follow-up observations have been carried out at the Isaac Newton Telescope for 22 stars, for which we obtained radial velocities and distance estimates. We discover 14 stars with a total velocity in the Galactic rest frame > 400 km/s, and 5 of these have a probability $>50\%$ of being unbound from the Milky Way. Tracing back their orbits in different Galactic potential models we find one possible unbound HVS with velocity $\sim$ 520 km/s, 5 bound HVSs, and, notably, 5 runaway stars with median velocity between 400 and 780 km/s. At the moment, uncertainties in the distance estimates and ages are too large to confirm the nature of our candidates by narrowing down their ejection location, and we wait for future Gaia releases to validate the quality of our sample. This test successfully demonstrates the feasibility of our new data mining routine.

The mass assembly history of the Milky Way can inform both theory of galaxy formation and the underlying cosmological model. Thus, observational constraints on the properties of both its baryonic and dark matter contents are sought. Here we show that hypervelocity stars (HVSs) can in principle provide such constraints. We model the observed velocity distribution of HVSs, produced by tidal break-up of stellar binaries caused by Sgr A*. Considering a Galactic Centre (GC) binary population consistent with that inferred in more observationally accessible regions, a fit to current HVS data with significance level > 5\% can only be obtained if the escape velocity from the GC to 50 kpc is $V_G < 850$ km/s, regardless of the enclosed mass distribution. When a NFW matter density profile for the dark matter halo is assumed, haloes with $V_G < 850$ km/s are in agreement with predictions in the $\Lambda$CDM model and that a subset of models around $M_{200} \sim 0.5-1.5 \times 10^{12}$ solar masses and $r_s < 35$ kpc can also reproduce Galactic circular velocity data. HVS data alone cannot currently exclude potentials with $V_G > 850$ km/s. Finally, specific constraints on the halo mass from HVS data are highly dependent on the assumed baryonic mass potentials. This first attempt to simultaneously constrain GC and dark halo properties is primarily hampered by the paucity and quality of data. It nevertheless demonstrates the potential of our method, that may be fully realised with the ESA Gaia mission.