• We use new Gaia measurements to explore the origin of stars from the Hypervelocity Star Survey as a function of velocity. Stars with >+450 km/s radial velocities, about >100 km/s above Galactic escape velocity at the depth of our survey, originate from the Galactic center. Two bound stars may also originate from the Galactic center. Galactic disk runaways appear at speeds around Galactic escape velocity, but most disk runaways are bound. Halo stars dominate the sample at speeds 100 km/s below Galactic escape velocity. Future Gaia measurements will enable a large, clean sample of Galactic center ejections for measuring the massive black hole ejection rate of HVSs, and for constraining the mass distribution of the Milky Way dark matter halo.
  • Recent Kepler photometry has revealed that about half of white dwarfs (WDs) have periodic, low-level (~ 1e-4 - 1e-3), optical variations. Hubble Space Telescope (HST) ultraviolet spectroscopy has shown that up to about one half of WDs are actively accreting rocky planetary debris, as evidenced by the presence of photospheric metal absorption lines. We have obtained HST ultraviolet spectra of seven WDs that have been monitored for periodic variations, to test the hypothesis that these two phenomena are causally connected, i.e. that the optical periodic modulation is caused by WD rotation coupled with an inhomogeneous surface distribution of accreted metals. We detect photospheric metals in four out of the seven WDs. However, we find no significant correspondence between the existence of optical periodic variability and the detection of photospheric ultraviolet absorption lines. Thus the null hypothesis stands, that the two phenomena are not directly related. Some other source of WD surface inhomogeneity, perhaps related to magnetic field strength, combined with the WD rotation, or alternatively effects due to close binary companions, may be behind the observed optical modulation. We report the marginal detection of molecular hydrogen in WD J1949+4734, only the fourth known WD with detected H_2 lines. We also re-classify J1926+4219 as a carbon-rich He-sdO subdwarf.
  • We report the discovery of two detached double white dwarf (WD) binaries, SDSS J082239.546+304857.19 and SDSS J104336.275+055149.90, with orbital periods of 40 and 46 min, respectively. The 40 min system is eclipsing; it is composed of a 0.30 Msun and a 0.52 Msun WD. The 46 min system is a likely LISA verification binary. The short 20 Myr and ~34 Myr gravitational wave merger times of the two binaries imply that many more such systems have formed and merged over the age of the Milky Way. We update the estimated Milky Way He+CO WD binary merger rate and affirm our previously published result: He+CO WD binaries merge at a rate at least 40 times greater than the formation rate of stable mass-transfer AM~CVn binaries, and so the majority must have unstable mass-transfer. The implication is that spin-orbit coupling in He+CO WD mergers is weak, or perhaps nova-like outbursts drive He+CO WDs into merger as proposed by Shen.
  • We present the results from a Gemini snapshot radial-velocity survey of 44 low-mass white dwarf candidates selected from the Sloan Digital Sky Survey spectroscopy. To find sub-hour orbital period binary systems, our time-series spectroscopy had cadences of 2 to 8 min over a period of 20-30 min. Through follow-up observations at Gemini and the MMT, we identify four double degenerate binary systems with periods ranging from 53 min to 7 h. The shortest period system, SDSS J123549.88+154319.3, was recently identified as a subhour period detached binary by Breedt and collaborators. Here we refine the orbital and physical parameters of this system. High-speed and time domain survey photometry observations do not reveal eclipses or other photometric effects in any of our targets. We compare the period distribution of these four systems with the orbital period distribution of known double white dwarfs; the median period decreases from 0.64 to 0.24 d for $M=0.3-0.5 M_{\odot}$ to $M<0.3 M_{\odot}$ white dwarfs. However, we do not find a statistically significant correlation between the orbital period and white dwarf mass.
  • We address the physical nature of subdwarf A-type (sdA) stars and their possible link to extremely low mass (ELM) white dwarfs (WDs). The two classes of objects are confused in low-resolution spectroscopy. However, colors and proper motions indicate that sdA stars are cooler and more luminous, and thus larger in radius, than published ELM WDs. We demonstrate that surface gravities derived from pure hydrogen models suffer a systematic ~1 dex error for sdA stars, likely explained by metal line blanketing below 9000 K. A detailed study of five eclipsing binaries with radial velocity orbital solutions and infrared excess establishes that these sdA stars are metal-poor ~1.2 Msun main sequence stars with ~0.8 Msun companions. While WDs must exist at sdA temperatures, only ~1% of a magnitude-limited sdA sample should be ELM WDs. We conclude that the majority of sdA stars are metal-poor A-F type stars in the halo, and that recently discovered pulsating ELM WD-like stars with no obvious radial velocity variations may be SX Phe variables, not pulsating WDs.
  • We assess the photometric variability of nine stars with spectroscopic Teff and log(g) values from the ELM Survey that locate them near the empirical extremely low-mass (ELM) white dwarf instability strip. We discover three new pulsating stars: SDSS J135512.34+195645.4, SDSS J173521.69+213440.6 and SDSS J213907.42+222708.9. However, these are among the few ELM Survey objects that do not show radial velocity variations to confirm the binary nature expected of helium-core white dwarfs. The dominant 4.31-hr pulsation in SDSS J135512.34+195645.4 far exceeds the theoretical cutoff for surface reflection in a white dwarf, and this target is likely a high-amplitude delta Scuti pulsator with an overestimated surface gravity. We estimate the probability to be less than 0.0008 that the lack of measured radial velocity variations in four of eight other pulsating candidate ELM white dwarfs could be due to low orbital inclination. Two other targets exhibit variability as photometric binaries. Partial coverage of the 19.342-hr orbit of WD J030818.19+514011.5 reveals deep eclipses that imply a primary radius > 0.4 solar radii--too large to be consistent with an ELM white dwarf. The only object for which our time series photometry adds support to the ELM white dwarf classification is SDSS J105435.78-212155.9, with consistent signatures of Doppler beaming and ellipsoidal variations. We interpret that the ELM Survey contains multiple false positives from another stellar population at Teff < 9000 K, possibly related to the sdA stars recently reported from SDSS spectra.
  • A catalog of 8472 white dwarf (WD) candidates is presented, selected using reduced proper motions from the deep proper motion catalog of Munn et al. 2014. Candidates are selected in the magnitude range 16 < r < 21.5 over 980 square degrees, and 16 < r < 21.3 over an additional 1276 square degrees, within the Sloan Digital Sky Survey (SDSS) imaging footprint. Distances, bolometric luminosities, and atmospheric compositions are derived by fitting SDSS ugriz photometry to pure hydrogen and helium model atmospheres (assuming surface gravities log g = 8). The disk white dwarf luminosity function (WDLF) is constructed using a sample of 2839 stars with 5.5 < M_bol < 17, with statistically significant numbers of stars cooler than the turnover in the luminosity function. The WDLF for the halo is also constructed, using a sample of 135 halo WDs with 5 < M_bol < 16. We find space densities of disk and halo WDs in the solar neighborhood of 5.5 +- 0.1 x 10^-3 pc^-3 and 3.5 +- 0.7 x 10^-5 pc^-3, respectively. We resolve the bump in the disk WDLF due to the onset of fully convective envelopes in WDs, and see indications of it in the halo WDLF as well.
  • We present optical spectroscopy and near-infrared photometry of 57 faint ($g= 19-22$) high proper motion white dwarfs identified through repeat imaging of $\approx3100$ square degrees of the Sloan Digital Sky Survey footprint by Munn et al. (2014). We use $ugriz$ and $JH$ photometry to perform a model atmosphere analysis, and identify ten ultracool white dwarfs with $T_{\rm eff}<4000$ K, including the coolest pure H atmosphere white dwarf currently known, J1657+2638, with $T_{\rm eff}= 3550 \pm 100$K. The majority of the objects with cooling ages larger than 9 Gyr display thick disc kinematics and constrain the age of the thick disc to $\geq11$ Gyr. There are four white dwarfs in our sample with large tangential velocities ($v_{tan} > 120$ km $s^{-1}$) and UVW velocities that are more consistent with the halo than the Galactic disc. For typical $0.6 M_{\odot}$ white dwarfs, the cooling ages for these halo candidates range from 2.3 to 8.5 Gyr. However, the total main-sequence + white dwarf cooling ages of these stars would be consistent with the Galactic halo if they are slightly undermassive. Given the magnitude limits of the current large scale surveys, many of the coolest and oldest white dwarfs remain undiscovered in the solar neighborhood, but upcoming surveys such as GAIA and the Large Synoptic Survey Telescope (LSST) should find many of these elusive thick disc and halo white dwarfs.
  • For sufficiently wide orbital separations {\it a}, the two members of a stellar binary evolve independently. This implies that in a wide double white dwarf (DWD), the more massive WD should always be produced first, when its more massive progenitor ends its main-sequence life, and should therefore be older and cooler than its companion. The bound, wide DWD HS 2220$+$2146 ($a\approx500$ AU) does not conform to this picture: the more massive WD is the younger, hotter of the pair. We show that this discrepancy is unlikely to be due to past mass-transfer phases or to the presence of an unresolved companion. Instead, we propose that HS 2220$+$2146 formed through a new wide DWD evolutionary channel involving the merger of the inner binary in a hierarchical triple system. The resulting blue straggler and its wide companion then evolved independently, forming the WD pair seen today. Although we cannot rule out other scenarios, the most likely formation channel has the inner binary merging while both stars are still on the main sequence. This provides us with the tantalizing possibility that Kozai-Lidov oscillations may have played a role in the inner binary's merger. {\it Gaia} may uncover hundreds more wide DWDs, leading to the identification of other systems like HS 2220$+$2146. There are already indications that other WD systems may have been formed through different, but related, hierarchical triple evolutionary scenarios. Characterizing these populations may allow for thorough testing of the efficiency with which KL oscillations induce stellar mergers.
  • We present Chandra and Swift X-ray observations of four extremely low-mass (ELM) white dwarfs with massive companions. We place stringent limits on X-ray emission from all four systems, indicating that neutron star companions are extremely unlikely and that the companions are almost certainly white dwarfs. Given the observed orbital periods and radial velocity amplitudes, the total masses of these binaries are greater than 1.02 to 1.39 Msun. The extreme mass ratios between the two components make it unlikely that these binary white dwarfs will merge and explode as Type Ia or underluminous supernovae. Instead, they will likely go through stable mass transfer through an accretion disk and turn into interacting AM CVn. Along with three previously known systems, we identify two of our targets, J0811 and J2132, as systems that will definitely undergo stable mass transfer. In addition, we use the binary white dwarf sample from the ELM Survey to constrain the inspiral rate of systems with extreme mass ratios. This rate, 0.00017/year, is consistent with the AM CVn space density estimated from the Sloan Digital Sky Survey. Hence, stable mass transfer double white dwarf progenitors can account for the entire AM CVn population in the Galaxy.
  • We report the discovery of pulsations in three mixed-atmosphere, extremely low-mass white dwarf (ELM WD, M $\leqslant$ 0.3 M$_{\odot}$) precursors. Following the recent discoveries of pulsations in both ELM and pre-ELM WDs, we targeted pre-ELM WDs with mixed H/He atmospheres with high-speed photometry. We find significant optical variability in all three observed targets with periods in the range 320--590 s, consistent in timescale with theoretical predictions of $p$-mode pulsations in mixed-atmosphere $\approx$ 0.18 M$_{\odot}$ He-core pre-ELM WDs. This represents the first empirical evidence that pulsations in pre-ELM WDs can only occur if a significant amount of He is present in the atmosphere. Future, more extensive, time-series photometry of the brightest of the three new pulsators offers an excellent opportunity to constrain the thickness of the surface H layer, which regulates the cooling timescales for ELM WDs.
  • We present the discovery of 15 extremely low mass (5 < log{g} < 7) white dwarf candidates, 9 of which are in ultra-compact double-degenerate binaries. Our targeted ELM Survey sample now includes 76 binaries. The sample has a lognormal distribution of orbital periods with a median period of 5.4 hr. The velocity amplitudes imply that the binary companions have a normal distribution of mass with 0.76 Msun mean and 0.25 Msun dispersion. Thus extremely low mass white dwarfs are found in binaries with a typical mass ratio of 1:4. Statistically speaking, 95% of the white dwarf binaries have a total mass below the Chandrasekhar mass and thus are not Type Ia supernova progenitors. Yet half of the observed binaries will merge in less than 6 Gyr due to gravitational wave radiation; probable outcomes include single massive white dwarfs and stable mass transfer AM CVn binaries.
  • We estimate the merger rate of double degenerate binaries containing extremely low mass (ELM) <0.3 Msun white dwarfs in the Galaxy. Such white dwarfs are detectable for timescales of 0.1 Gyr -- 1 Gyr in the ELM Survey; the binaries they reside in have gravitational wave merger times of 0.001 Gyr -- 100 Gyr. To explain the observed distribution requires that most ELM white dwarf binary progenitors detach from the common envelope phase with <1 hr orbital periods. We calculate the local space density of ELM white dwarf binaries and estimate a merger rate of 3e-3/yr over the entire disk of the Milky Way; the merger rate in the halo is 10 times smaller. The ELM white dwarf binary merger rate exceeds by a factor of 40 the formation rate of stable mass transfer AM CVn binaries, marginally exceeds the rate of underluminous supernovae, and is identical to the formation rate of R CrB stars. On this basis, we conclude that ELM white dwarf binaries can be the progenitors of all observed AM CVn and possibly underluminous supernovae, however the majority of He+CO white dwarf binaries go through unstable mass transfer and merge, e.g. into single massive ~1 Msun white dwarfs.
  • We report evidence for excess blue light from the Type Ia supernova SN 2012cg at fifteen and sixteen days before maximum B-band brightness. The emission is consistent with predictions for the impact of the supernova on a non-degenerate binary companion. This is the first evidence for emission from a companion to a SN Ia. Sixteen days before maximum light, the B-V color of SN 2012cg is 0.2 mag bluer than for other normal SN~Ia. At later times, this supernova has a typical SN Ia light curve, with extinction-corrected M_B = -19.62 +/- 0.02 mag and Delta m_{15}(B) = 0.86 +/- 0.02. Our data set is extensive, with photometry in 7 filters from 5 independent sources. Early spectra also show the effects of blue light, and high-velocity features are observed at early times. Near maximum, the spectra are normal with a silicon velocity v_{Si} = -10,500$ km s^{-1}. Comparing the early data with models by Kasen (2010) favors a main-sequence companion of about 6 solar masses. It is possible that many other SN Ia have main-sequence companions that have eluded detection because the emission from the impact is fleeting and faint.
  • We present the serendipitous discovery of eclipse-like events around the massive white dwarf SDSS J152934.98+292801.9 (hereafter J1529+2928). We selected J1529+2928 for time-series photometry based on its spectroscopic temperature and surface gravity, which place it near the ZZ Ceti instability strip. Instead of pulsations, we detect photometric dips from this white dwarf every 38 minutes. Follow-up optical spectroscopy observations with Gemini reveal no significant radial velocity variations, ruling out stellar and brown dwarf companions. A disintegrating planet around this white dwarf cannot explain the observed light curves in different filters. Given the short period, the source of the photometric dips must be a dark spot that comes into view every 38 min due to the rotation of the white dwarf. Our optical spectroscopy does not show any evidence of Zeeman splitting of the Balmer lines, limiting the magnetic field strength to B<70 kG. Since up to 15% of white dwarfs display kG magnetic fields, such eclipse-like events should be common around white dwarfs. We discuss the potential implications of this discovery on transient surveys targeting white dwarfs, like the K2 mission and the Large Synoptic Survey Telescope.
  • White dwarfs are the end state of most stars, including the Sun, after they exhaust their nuclear fuel. Between 1/4 and 1/2 of white dwarfs have elements heavier than helium in their atmospheres, even though these elements should rapidly settle into the stellar interiors unless they are occasionally replenished. The abundance ratios of heavy elements in white dwarf atmospheres are similar to rocky bodies in the Solar system. This and the existence of warm dusty debris disks around about 4% of white dwarfs suggest that rocky debris from white dwarf progenitors' planetary systems occasionally pollute the stars' atmospheres. The total accreted mass can be comparable to that of large asteroids in the solar system. However, the process of disrupting planetary material has not yet been observed. Here, we report observations of a white dwarf being transited by at least one and likely multiple disintegrating planetesimals with periods ranging from 4.5 hours to 4.9 hours. The strongest transit signals occur every 4.5 hours and exhibit varying depths up to 40% and asymmetric profiles, indicative of a small object with a cometary tail of dusty effluent material. The star hosts a dusty debris disk and the star's spectrum shows prominent lines from heavy elements like magnesium, aluminium, silicon, calcium, iron, and nickel. This system provides evidence that heavy element pollution of white dwarfs can originate from disrupted rocky bodies such as asteroids and minor planets.
  • We present the discovery of 11 new double degenerate systems containing extremely low-mass white dwarfs (ELM WDs). Our radial velocity observations confirm that all of the targets have orbital periods $\leq$ 1 day. We perform spectroscopic fits and provide a complete set of physical and binary parameters. We review and compare recent evolutionary calculations and estimate that the systematic uncertainty in our mass determinations due to differences in the evolutionary models is small ($\approx$ 0.01 M$_{\odot}$). Five of the new systems will merge due to gravitational wave radiation within a Hubble time, bringing the total number of merger systems found in the ELM Survey to 38. We examine the ensemble properties of the current sample of ELM WD binaries, including the period distribution as a function of effective temperature, and the implications for the future evolution of these systems. We also revisit the empirical boundaries of instability strip of ELM WDs and identify new pulsating ELM WD candidates. Finally, we consider the kinematic properties of our sample of ELM WDs and estimate that a significant fraction of the WDs from the ELM Survey are members of the Galactic halo.
  • We measure the three components of velocity dispersion, $\sigma_{R},\sigma_{\theta},\sigma_{\phi}$, for stars within 6 < R < 30 kpc of the Milky Way using a new radial velocity sample from the MMT telescope. We combine our measurements with previously published data so that we can more finely sample the stellar halo. We use a maximum likelihood statistical method for estimating mean velocities, dispersions, and covariances assuming only that velocities are normally distributed. The alignment of the velocity ellipsoid is consistent with a spherically symmetric gravitational potential. From the spherical Jeans equation, the mass of the Milky Way is M(<12 kpc) = $1.3\times10^{11}$ M$_{\odot}$ with an uncertainty of 40%. We also find a region of discontinuity, 15 < R < 25 kpc, where the estimated velocity dispersions and anisotropies diverge from their anticipated values, confirming at high significance the break observed by others. We argue that this break in anisotropy is physically explained by coherent stellar velocity structure in the halo, such as the Sgr stream. To significantly improve our understanding of halo kinematics will require combining radial velocities with future Gaia proper motions.
  • We measure proper motions with the Hubble Space Telescope for 16 extreme radial velocity stars, mostly unbound B stars in the Milky Way halo. Twelve of these stars have proper motions statistically consistent with zero, and thus have radial trajectories statistically consistent with a Galactic center "hypervelocity star" origin. The trajectory of HE 0437-5439 is consistent with both Milky Way and Large Magellanic Cloud origins. A Galactic center origin is excluded at 3-sigma confidence for two of the lowest radial velocity stars in our sample, however. These two stars are probable disk runaways and provide evidence for ~500 km/s ejections from the disk. We also measure a significant proper motion for the unbound sdO star US 708. Its 1,000 km/s motion is in some tension with proposed supernova ejection models, but can be explained if US 708 was ejected from the stellar halo. In the future, we expect Gaia will better constrain the origin of these remarkable unbound stars.
  • We report the discovery of the first millisecond pulsar with a pulsating white dwarf companion. Following the recent discoveries of pulsations in extremely low-mass (ELM, <0.3 Msol) white dwarfs (WDs), we targeted ELM WD companions to two millisecond pulsars with high-speed Gemini photometry. We find significant optical variability in PSR J1738+0333 with periods between roughly 1790-3060 s, consistent in timescale with theoretical and empirical observations of pulsations in 0.17 Msol He-core ELM WDs. We additionally put stringent limits on a lack of variability in PSR J1909-3744, showing this ELM WD is not variable to <0.1 per cent amplitude. Thanks to the accurate distance and radius estimates from radio timing measurements, PSR J1738+0333 becomes a benchmark for low-mass, pulsating WDs. Future, more extensive time-series photometry of this system offers an unprecedented opportunity to constrain the physical parameters (including the cooling age) and interior structure of this ELM WD, and in turn, the mass and spin-down age of its pulsar companion.
  • We predict the distinctive three dimensional space motions of hypervelocity stars (HVSs) and runaway stars moving in a realistic Galactic potential. For nearby stars with distances less than 10~kpc, unbound stars are rare; proper motions alone rarely isolate bound HVSs and runaways from indigenous halo stars. At large distances of 20-100 kpc, unbound HVSs are much more common than runaways; radial velocities easily distinguish both from indigenous halo stars. Comparisons of the predictions with existing observations are encouraging. Although the models fail to match observations of solar-type HVS candidates from SEGUE, they agree well with data for B-type HVS and runaways from other surveys. Complete samples of g <= 20 stars with GAIA should provide clear tests of formation models for HVSs and runaways and will enable accurate probes of the shape of the Galactic potential.
  • We carry out high-speed photometry on 20 of the shortest-period, detached white dwarf binaries known and discover systems with eclipses, ellipsoidal variations (due to tidal deformations of the visible white dwarf), and Doppler beaming. All of the binaries contain low-mass white dwarfs with orbital periods less than 4 hr. Our observations identify the first eight tidally distorted white dwarfs, four of which are reported for the first time here, which we use to put empirical constraints on the mass-radius relationship for extremely low-mass (<0.30 Msun) white dwarfs. We also detect Doppler beaming in several of these binaries, which confirms the high-amplitude radial-velocity variability. All of these systems are strong sources of gravitational radiation, and long-term monitoring of those that display ellipsoidal variations can be used to detect spin-up of the tidal bulge due to orbital decay.
  • We report the discovery of a detached 20 min orbital period binary white dwarf. WD0931+444 (SDSS J093506.93+441106.9) was previously classified as a WD + M dwarf system based on its optical spectrum. Our time-resolved optical spectroscopy observations obtained at the 8m Gemini and 6.5m MMT reveal peak-to-peak radial velocity variations of 400 km/s every 20 min for the WD, but no velocity variations for the M dwarf. In addition, high-speed photometry from the McDonald 2.1m telescope shows no evidence of variability nor evidence of a reflection effect. An M dwarf companion is physically too large to fit into a 20 min orbit. Thus, the orbital motion of the WD is almost certainly due to an invisible WD companion. The M dwarf must be either an unrelated background object or the tertiary component of a hiearchical triple system. WD0931+444 contains a pair of WDs, a 0.32 Msol primary and a >0.14 Msol secondary, at a separation of >0.19 Rsol. After J0651+2844, WD0931+444 becomes the second-shortest period detached binary WD currently known. The two WDs will lose angular momentum through gravitational wave radiation and merge in <9 Myr. The log h ~ -22 gravitational wave strain from WD0931+444 is strong enough to make it a verification source for gravitational wave missions in the milli-Hertz frequency range, e.g. the evolved Laser Interferometer Space Antenna (eLISA), bringing the total number of known eLISA verification sources to nine.
  • A recent cross-correlation between the SDSS DR7 White Dwarf Catalog with the Wide-Field Infrared Survey Explorer ($WISE$) all-sky photometry at 3.4, 4.6, 12, and 22 microns performed by Debes et al. (2011) resulted in the discovery of 52 candidate dusty white dwarfs (WDs). The 6'' $WISE$ beam allows for the possibility that many of the excesses exhibited by these WDs may be due to contamination from a nearby source, however. We present MMT$+$SWIRC $J$- and $H$-band imaging observations (0.5-1.5'' PSF) of 16 of these candidate dusty WDs and confirm that four have spectral energy distributions (SEDs) consistent with a dusty disk and are not accompanied by a nearby source contaminant. The remaining 12 WDs have contaminated $WISE$ photometry and SEDs inconsistent with a dusty disk when the contaminating sources are not included in the photometry measurements. We find the frequency of disks around single WDs in the $WISE$ $\cap$ SDSS sample to be 2.6-4.1%. One of the four new dusty WDs has a mass of $1.04 M_{\odot}$ (progenitor mass $5.4 M_{\odot}$) and its discovery offers the first confirmation that massive WDs (and their massive progenitor stars) host planetary systems.
  • We describe our completed spectroscopic survey for unbound hypervelocity stars (HVSs) ejected from the Milky Way. Three new discoveries bring the total number of unbound HVSs to 21. We place new constraints on the nature of HVSs and on their distances using moderate resolution MMT spectroscopy. Half of the HVSs are fast rotators; they are certain 2.5-4 Msun main sequence stars at 50 - 120 kpc distances. Correcting for stellar lifetime, our survey implies that unbound 2.5-4 Msun stars are ejected from the Milky Way at a rate of 1.5e-6 /yr. The observed HVSs are likely ejected continuously over the past 200 Myr and do not share a common flight time. The anisotropic spatial distribution of HVSs on the sky remains puzzling. Southern hemisphere surveys like SkyMapper will soon allow us to map the all-sky distribution of HVSs. Future proper motion measurements with Hubble Space Telescope and Gaia will provide strong constraints on origin. All existing observations are consistent with HVS ejections from encounters with the massive black hole in the Galactic center.