No Arabic abstract
We present and analyze the positions, distances, and radial velocities for over 4000 blue horizontal-branch (BHB) stars in the Milky Ways halo, drawn from SDSS DR8. We search for position-velocity substructure in these data, a signature of the hierarchical assembly of the stellar halo. Using a cumulative close pair distribution (CPD) as a statistic in the 4-dimensional space of sky position, distance, and velocity, we quantify the presence of position-velocity substructure at high statistical significance among the BHB stars: pairs of BHB stars that are close in position on the sky tend to have more similar distances and radial velocities compared to a random sampling of these overall distributions. We make analogous mock-observations of 11 numerical halo formation simulations, in which the stellar halo is entirely composed of disrupted satellite debris, and find a level of substructure comparable to that seen in the actually observed BHB star sample. This result quantitatively confirms the hierarchical build-up of the stellar halo through a signature in phase (position-velocity) space. In detail, the structure present in the BHB stars is somewhat less prominent than that seen in most simulated halos, quite possibly because BHB stars represent an older sub-population. BHB stars located beyond 20 kpc from the Galactic center exhibit stronger substructure than at $rm r_{gc} < 20$ kpc.
We use $N$-body simulations to investigate the excitation of bending waves in a Milky Way-like disc-bulge-halo system. The dark matter halo consists of a smooth component and a population of subhaloes while the disc is composed of thin and thick components. Also considered is a control simulation where all of the halo mass is smoothly distributed. We find that bending waves are more vigorously excited in the thin disc than the thick one and that they are strongest in the outer regions of the disc, especially at late times. By way of a Fourier decomposition, we find that the complicated pattern of bending across the disc can be described as a superposition of waves, which concentrate along two branches in the radius-rotational frequency plane. These branches correspond to vertical resonance curves as predicted by a WKB analysis. Bending waves in the simulation with substructure have a higher amplitude than those in the smooth-halo simulation, though the frequency-radius characteristics of the waves in the two simulations are very similar. A cross correlation analysis of vertical displacement and bulk vertical velocity suggests that the waves oscillate largely as simple plane waves. We suggest that the wave-like features in astrometric surveys such as the Second Data Release from textit{Gaia} may be due to long-lived waves of a dynamically active disc rather than, or in addition to, perturbations from a recent satellite-disc encounter.
Tidal debris from infalling satellites can leave observable structure in the phase-space distribution of the Galactic halo. Such substructure can be manifest in the spatial and/or velocity distributions of the stars in the halo. This paper focuses on a class of substructure that is purely kinematic in nature, with no accompanying spatial features. To study its properties, we use a simulated stellar halo created by dynamically populating the Via Lactea II high-resolution N-body simulation with stars. A significant fraction of the stars in the inner halo of Via Lactea share a common speed and metallicity, despite the fact that they are spatially diffuse. We argue that this kinematic substructure is a generic feature of tidal debris from older mergers and may explain the detection of radial-velocity substructure in the inner halo made by the Sloan Extension for Galactic Understanding and Exploration. The GAIA satellite, which will provide the proper motions of an unprecedented number of stars, should further characterize the kinematic substructure in the inner halo. Our study of the Via Lactea simulation suggests that the stellar halo can be used to map the speed distribution of the local dark-matter halo, which has important consequences for dark-matter direct-detection experiments.
In the $Gaia$ era stellar kinematics are extensively used to study Galactic halo stellar populations, to search for halo structures, and to characterize the interface between the halo and hot disc populations. We use distribution function-based models of modern datasets with 6D phase space data to qualitatively describe a variety of kinematic spaces commonly used in the study of the Galactic halo. Furthermore, we quantitatively assess how well each kinematic space can separate radially anisotropic from isotropic halo populations. We find that scaled action space (the ``action diamond) is superior to other commonly used kinematic spaces at this task. We present a new, easy to implement selection criterion for members of the radially-anisotropic $Gaia$-Enceladus merger remnant, which we find achieves a sample purity of 82 per cent in our models with respect to contamination from the more isotropic halo. We compare this criterion to literature criteria, finding that it produces the highest purity in the resulting samples, at the expense of a modest reduction in completeness. We also show that selection biases that underlie nearly all contemporary spectroscopic datasets can noticeably impact the $E-L_{z}$ distribution of samples in a manner that may be confused for real substructure. We conclude by providing recommendations for how authors should use stellar kinematics in the future to study the Galactic stellar halo.
We have measured the amount of kinematic substructure in the Galactic halo using the final data set from the Spaghetti project, a pencil-beam high latitude sky survey. Our sample contains 101 photometrically selected and spectroscopically confirmed giants with accurate distance, radial velocity and metallicity information. We have developed a new clustering estimator: the 4distance measure, which when applied to our data set leads to the identification of 1 group and 7 pairs of clumped stars. The group, with 6 members, can confidently be matched to tidal debris of the Sagittarius dwarf galaxy. Two pairs match the properties of known Virgo structures. Using models of the disruption of Sagittarius in Galactic potentials with different degrees of dark halo flattening, we show that this favors a spherical or prolate halo shape, as demonstrated by Newberg et al. (2007) using SDSS data. One additional pair can be linked to older Sagittarius debris. We find that 20% of the stars in the Spaghetti data set are in substructures. From comparison with random data sets we derive a very conservative lower limit of 10% to the amount of substructure in the halo. However, comparison to numerical simulations shows that our results are also consistent with a halo entirely built up from disrupted satellites, provided the dominating features are relatively broad due to early merging or relatively heavy progenitor satellites.
We present a new, high-resolution chronographic (age) map of the Milky Ways halo, based on the inferred ages of ~130,000 field blue horizontal-branch (BHB) stars with photometry from the Sloan Digital Sky Survey. Our map exhibits a strong central concentration of BHB stars with ages greater than 12 Gyr, extending up to ~15 kpc from the Galactic center (reaching close to the solar vicinity), and a decrease in the mean ages of field stars with distance by 1-1.5 Gyr out to ~45-50 kpc, along with an apparent increase of the dispersion of stellar ages, and numerous known (and previously unknown) resolved over-densities and debris streams, including the Sagittarius Stream. These results agree with expectations from modern LambdaCDM cosmological simulations, and support the existence of a dual (inner/outer) halo system, punctuated by the presence of over-densities and debris streams that have not yet completely phase-space mixed.