---
abstract: '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.'
author:
- Mariangela Lisanti
- 'David N. Spergel'
- Piero Madau
bibliography:
- 'starsBibFile.bib'
nocite: '[@*]'
title: Signatures of Kinematic Substructure in the Galactic Stellar Halo
---

Introduction
============

The process of galaxy formation alters the phase-space distribution of the dark matter (DM) and stellar components of the Milky Way (MW) halo. The nature of the residual phase-space structure in the halo depends on the details of its formation history, and is sensitive to whether the galaxy formed from smooth collapse [@Eggen:1962dj] or from the merger of many protogalactic fragments [@Searle:1978gc]. The $\Lambda$CDM paradigm currently provides the most well-motivated picture of MW formation, including both the dark and baryonic matter as basic ingredients. In the $\Lambda$CDM framework, the MW halo forms from the hierarchical merging of smaller subhalos [@White:1977jf]. The subhalos are tidally disrupted as they fall into the host, and DM is torn off, along with stars that formed in the dense cores of the subhalos.

Tidal remnants from a completely disrupted subhalo eventually come into equilibrium with the host halo. An incomplete merging event, however, leaves tidal debris with phase-space structure distinguishable from that of the smooth equilibrated halo. Dwarf galaxies are examples of infalling satellites that have not been completely disrupted. These dwarfs orbit about the MW, leaving tidal debris in their wake, especially near the turning points of their orbits where the tidal forces are strongest. After tidal stripping, the debris exhibits distinctive structure in both position and velocity. With time, the debris comes into equilibrium with the host, and any distinctive phase-space features are washed out.

The accretion events that build the MW stellar halo leave their imprint on the phase-space distribution of its constituent stars [@Johnston:1996sb; @Johnston:1997fv; @Venn:2004hk]. This structure persists for some time because stars exchange energy and momenta on timescales that are much longer than the age of the Galaxy [@Helmi:2008eq]. As a result, structure in the stellar halo serves as a fossil record of the MW’s formation history and kinematic or spatial features may be indicative of one or more merger events.

The chemical composition of stars provides additional information about their origin [@McWilliam:1997ua; @Robertson:2005gv; @Font:2005qs]. The metal content is particularly indicative because iron is introduced into the interstellar medium from supernova explosions, and is thus related to the total integrated star formation. The chemical properties of stars brought into the MW halo depend on the mass of their subhalo host because a more massive subhalo has had more time to form stars, and also retains more metal. The stars that are deposited in the MW by such mergers are typically more metal-rich than those deposited earlier by smaller subhalos [@Unavane:1995an].

Evidence for stellar substructure has been accumulating with the advent of large-scale surveys, such as the Sloan Digital Sky Survey (SDSS) [@Fukugita:1996qt; @Gunn:1998vh; @York:2000gk; @Smith:2002pca; @Pier:2002iq; @Ivezic:2004bf; @Tucker:2006dv; @Gunn:2006tw], the Sloan Extension for Galactic Understanding and Exploration (SEGUE) [@AllendePrieto:2007ed; @Lee:2007mf; @Lee:2007ec; @Yanny:2009kg], the Spaghetti survey [@Morrison:2000gp], the Two Micron All Sky Survey [@Skrutskie:2006wh], the Quasar Equatorial Survey Team [@1999AAS...195.1502C; @2012MNRAS.427.3374M], and the Radial Velocity Experiment [@Steinmetz:2006qt; @Zwitter:2008vc]. The Sagittarius dwarf [@Ibata:1994fv] is one of the most stunning examples of stellar substructure from an on-going accretion event; the dwarf is in the midst of a merger with the MW and the tidal stream it has left in its orbital wake has been mapped to amazing precision [@Johnston:1995vd; @Ivezic:2000ua; @Yanny:2000ty; @Ibata:2000pu; @Vivas:2001dn; @Majewski:2003ux]. Many other examples of stellar substructure have been documented [@Majewski:1996zz; @Chiba:1997ta; @Totten:1998; @Helmi:1999uj; @Helmi:1999ks; @Totten:2000ab; @Chiba:2000vu; @Newberg:2001sx; @Gilmore:2002jv; @Helmi:2002iu; @RochaPinto:2004ru; @Grillmair:2006bd; @Vivas:2006nh; @Kepley:2007vx; @Juric:2005zr; @Klement:2008ws; @Seabroke:2008; @Watkins:2009im; @Grillmair:2008fv; @Starkenburg:2009nd; @An:2009hj; @Smith:2009kr; @Klement:2009km; @Schlaufman:2009jv; @Harrigan:2010pd; @Schlaufman:2011kf; @Johnston:2012yh; @Sheffield:2012yg]—*e.g.*, the Monoceros “Ring" [@Yanny:2003zu; @Ibata:2003di; @Ivezic:2008wk], the Orphan Stream [@Belokurov:2006ms; @Grillmair:2006nx; @Belokurov:2006kc], the Virgo Stellar Stream [@Duffau:2005ta; @Zinn:2013saa], and tidal tails near the  [@Odenkirchen:2000zx; @Rockosi:2002wu] and  [@Belokurov:2005ad; @Grillmair:2006se] globular clusters.

It is well-established that stellar substructure in position-space is washed out before substructure in velocity-space [@Johnston:1997fv; @Helmi:1999ks; @Freeman:2002wq]. A stream is an example of substructure that is coherent in both position- and velocity-space. Compared to streams, purely kinematic substructure is older and thus provides a way to probe earlier mergers in the halo. In this paper, we focus on a particular class of purely kinematic substructure—termed “debris flow"—whose constituents are not coherent in velocity, but still exhibit a characteristic *speed*. Debris flow is intermediate between an equilibrated halo and dynamically cold tidal stream.

Debris flow has already been studied in the context of DM [@Lisanti:2011as; @Kuhlen:2012fz] and this work is the first to explore its signatures in the Galactic stellar halo. In particular, we study a simulated stellar halo created by dynamically populating a DM-only $N$-body simulation with stars (Sec. \[sec: VL2tagging\]). The simulated halo exhibits structure in velocity and metallicity that distinguishes it from the smooth-halo expectation (Sec. \[sec: StellarFlows\]). In particular, the stars share a common speed and metallicity, determined by the mass and orbital properties of its host satellite, despite being spatially diffuse.

The observational signatures of debris flow, as well as its implications for current and future surveys, are described in Sec. \[sec: ECHOS\]. A recent study of the SEGUE data [@Schlaufman:2009jv; @Schlaufman:2011kf] found evidence for radial-velocity substructure in the inner halo. The SEGUE study uses a very large sample of metal-poor stars and finds several high-confidence detections that are not associated with any known streams. As we show, the SEGUE findings could potentially be explained by debris flow, although confirmation will require more complete proper-motion measurements from the GAIA satellite [@Perryman:2001sp; @2008IAUS..248..217L].

The mounting evidence for phase-space structure in the stellar halo supports the picture that the MW formed through hierarchical mergers. In addition, it strongly suggests that the Galaxy’s DM is not smoothly distributed in phase space. Ideally, the identification of stellar substructure in the Solar neighborhood can be used to infer information about the local DM speed distribution, which is relevant for direct-detection experiments. We conclude by discussing the implications of stellar debris flow for DM (Sec. \[sec:DMImplications\]).

The Via Lactea II Stellar Halo {#sec: VL2tagging}
==============================

Via Lactea-II (VL2) is a DM-only $N$-body simulation that employs about a billion $4.1\times10^3 \text{ M}_\odot$ particles to model the formation of a $M_{200}=1.9\times10^{12}$ M$_{\odot}$ MW-size halo and its substructure [@Diemand:2008in; @Zemp:2008gw]. It is initialized at $z=104.3$ to a WMAP3 $\Lambda$CDM cosmology [@Spergel:2006hy] and evolved to the present day. Twenty-seven “snapshots" of the simulation from to , spaced roughly 680 Myr apart, were analyzed in detail. The subhalos in each snapshot were identified using the 6DFOF group finder [@Diemand:2006ey]. The progression of all 3200 subhalos with infall mass greater than $10^7$ M$_\odot$ was tracked from snapshot to snapshot to build the evolutionary history.

To study the stellar counterpart of the VL2 DM halo, the simulation output was dynamically populated with stars [@Rashkov:2011cq]. The tagging prescription labeled the most tightly bound DM particles in each subhalo as stars, assigning a mass, ${\text{m}_{\text{sp}}}$,[^1] and metallicity, ${\left[\text{Fe/H}\right]}$, to each. The tagging prescription was tuned to reproduce the luminosity function, the half-light radii, and the metallicity-luminosity relation of the MW’s observed satellites. We briefly review the tagging procedure here; for a more detailed discussion, see [@Rashkov:2011cq] and references therein.

The total stellar mass of the subhalo at infall, $\text{M}_*$, was assumed to follow the power-law $$\frac{\text{M}_*}{{\text{M}_{\text{h}}}} = 1.6\times10^{-5} \left( \frac{{\text{M}_{\text{h}}}}{10^9 {\text{ M}_{\odot}}}\right)^{1.8} \, ,
\label{eq:totalmass}$$ where $M_h$ is the maximum mass of the subhalo’s host. $\text{M}_*/{\text{M}_{\text{h}}}$ is the stellar formation efficiency and determines the satellite luminosity today. The total stellar mass was distributed amongst the 1% most-bound particles in the subhalo to get ${\text{m}_{\text{sp}}}$. This number determines the concentration of the stellar system at infall, and governs the amount of stellar material stripped at later times as well as the present-day structural properties of the surviving satellites. Tagging the 1% highest total binding energy particles with stellar population provides a good fit to the distribution of half-light radii in MW dSphs [@Rashkov:2011cq].

The metallicity of each star was taken to be $${\left[\text{Fe/H}\right]}= -7.87 + 0.9\times \log \left( \frac{{\text{M}_{\text{h}}}}{10^3 {\text{ M}_{\odot}}} \, \right).
\label{eq:RashkovMet}$$ Stars in more massive subhalos have larger metallicities because their hosts retain more enriched material. All stars in a subhalo are assigned the same metallicity. Two processes—stellar mass loss from tidal stripping and the dimming of the stellar population with age—turn the above assumed scaling into the observed present-day luminosity-metallicity relation of [@Kirby:2008ab].

![The mass density distribution of all unbound stars in VL2 as a function of galactocentric distance (solid black). The individual contributions from stars with (dotted blue), (dashed green), and (dot-dashed red) are also shown. The shaded gray region indicates density fall-offs $\rho \propto r^{-\alpha}$ with $2<\alpha<4$, consistent with [@Bell:2007ts].](density_distribution.pdf){width="3.25in"}

\[fig:halodensity\]

![image](echosSP.pdf){width="6.5in"}

\[fig:scatter\]

The stellar tagging procedure results in a list of 1.6 million stars in the VL2 host halo, along with the position, velocity, and metallicity of each. About $31\%$ of the stars in the VL2 stellar halo are gravitationally bound to satellites at $z=0$ within , where $R_\text{gc}$ is galactocentric distance. For , the bound fraction is only about 4$\%$. The bound stars in the VL2 halo reproduce the sizes and brightness profiles of the observed dSph population, as well as their metallicities, velocity dispersions, and spatial distributions, as studied in [@Rashkov:2011cq].

The primary focus of this paper is the unbound population of stars in VL2. The total mass and density distribution of the VL2 stellar halo reproduce the measured quantities from SDSS. For example, as noted by [@Rashkov:2011cq], the $5\times10^8{\text{ M}_{\odot}}$ mass of the VL2 stellar halo compares well with observations in the same volume [@Bell:2007ts]. Figure \[fig:halodensity\] shows the mass density of all VL2 stars as a function of galactocentric distance (solid black). The distribution exhibits a power-law behavior consistent with that observed by SDSS (indicated by the shaded gray band) [@Bell:2007ts]. Below , the VL2 density distribution is dominated by stars with (dashed green). From , the radial fall-off is determined by stars with (dot-dashed red). Note that the maximum metallicity found in the VL2 stellar halo is .

Figure \[fig:scatter\] shows the galactocentric speeds of a random sampling of 3% of the stars in VL2 from , as a function of infall redshift (left panel) and metallicity (right panel). Each point denotes a single star “particle"; its color and relative size indicates the fraction of the total stellar mass in the volume from ($f_\text{GC}$) that it constitutes. The VL2 stars can be separated into approximately two populations based on their metallicities. The first population arises from small-mass subhalos that fell into the MW at $z_\text{infall} \gtrsim 4.5$ and is more metal-poor than the population arising from more massive progenitors that merged at .[^2]

   $z_\text{infall}$     $f_\text{GC}$        $f_\text{IH}$
  ------------------- -------------------- --------------------
         2.00          $8.5\times10^{-1}$   $8.1\times10^{-1}$
         2.50          $1.8\times10^{-2}$   $2.4\times10^{-2}$
         3.24          $6.4\times10^{-3}$   $3.0\times10^{-3}$
         4.56          $9.7\times10^{-2}$   $1.4\times10^{-1}$
         7.77          $2.9\times10^{-2}$   $1.4\times10^{-2}$
         9.14          $1.1\times10^{-3}$   $5.8\times10^{-4}$
         11.2          $6.2\times10^{-5}$   $4.8\times10^{-5}$
         12.7          $3.8\times10^{-9}$   $7.2\times10^{-9}$

  : The fraction of stellar mass that fell into the VL2 Milky Way host at redshift $z_\text{infall}$. $f_\text{GC}$ is the fraction within of the Galactic center; $f_\text{IH}$ is the fraction of stellar mass within the inner halo, defined in . Over $\sim$80$\%$ of the stellar mass was accreted at $z_\text{infall} = 2$ in either volume.[]{data-label="tab:snapshots"}

The late-accreting and more metal-rich population in the VL2 stellar halo dominates its mass density. Table \[tab:snapshots\] summarizes the fraction of stellar mass between ($f_\text{GC}$) that fell in at redshift $z_\text{infall}$. For comparison, Table \[tab:snapshots\] also provides the fraction ($f_\text{IH}$) of stellar mass within the inner VL2 halo, defined, as in [@Schlaufman:2009jv], to be $$\text{Inner Halo} = \left\{
     \begin{array}{l} 
        \vspace{5pt}
       |z| >  4 \text{ kpc}\\ \vspace{5pt}
       R_\text{gc} > 10 \text{ kpc} \\ 
       d_\text{hel} < 17.5 \text{ kpc} \, ,
     \end{array}
          \right. 
          \label{eq:innerhalo}$$ where $z$ is the vertical distance from the Galactic plane and $d_\text{hel}$ is the distance from the Sun. We will come back to this volume in Sec. \[sub:RVsubstructure\], when we compare the VL2 stellar halo with observations of the inner MW halo.

Nearly 97% of the stellar mass from was accreted after , with $\sim$85% accreted around . Indeed, the vast majority of the stellar halo in VL2 originates from a subhalo of stellar mass $\text{M}_* = 8.8\times10^8 \text{ M}_\odot$ that fell in at $z=2$ [@Rashkov:2011cq]. Using (\[eq:totalmass\]) and (\[eq:RashkovMet\]), we estimate that this subhalo has a total mass of ${\text{M}_{\text{h}}}= 5\times10^{10} \text{ M}_\odot$ and its stellar constituents have ${\left[\text{Fe/H}\right]}= -0.94$. This is in agreement with the large stellar mass density observed for ${\left[\text{Fe/H}\right]}\sim-1$ in the right panel of Fig. \[fig:scatter\].

From this point forward, we will divide the VL2 stars into two separate populations with $${\left[\text{Fe/H}\right]}< -1.8 \quad \text{and} \quad {\left[\text{Fe/H}\right]}> -1.8  \, .$$ A metallicity of ${\left[\text{Fe/H}\right]}\sim-1.8$ is roughly where the average galactocentric speed of the stars changes from to . The metallicity distribution for each population of VL2 stars within is shown in Figure \[fig:velmetdist\]. The more metal-rich population is narrowly peaked, with mean . The more metal-poor population has a broader distribution with mean .

There is a clear correlation between the kinematic and chemical properties of the stars in the VL2 halo. Whether the actual stellar halo in the MW can be divided into two distinct populations—and what the relative fraction of each population is—depends entirely on its accretion history. The stellar halo in VL2 is the manifestation of one possible accretion history and the observed features do not need to be generic. Indeed, as was pointed out by [@Schlaufman:2012ki], the observed metallicity distributions in the Inner Halo are not consistent with VL2. That being said, VL2 still serves as a valuable guide to better understand the correlations between the kinematic and chemical properties of stars in the halo.

Stellar Debris Flow {#sec: StellarFlows}
===================

Next, we analyze the metallicity and galactocentric velocity distributions of the stars in the VL2 halo. The top row of Fig. \[fig:velocitypanel\] shows the galactocentric speed distributions for the VL2 halo stars with ${\left[\text{Fe/H}\right]}<-1.8$ (left panel) and ${\left[\text{Fe/H}\right]}>-1.8$ (middle panel). In the radial bin, the ${\left[\text{Fe/H}\right]}>-1.8$ population exhibits a narrow peak in galactocentric speed about the mean . The ${\left[\text{Fe/H}\right]}<-1.8$ population has a much broader distribution in speed, but also exhibits a peak about . We focus primarily on the kinematic substructure observed in the more metal-rich population because the total mass density of these stars is much greater than that of the metal-poor population.

To compare the speed distributions obtained from VL2 to that expected for a smooth stellar halo, we simulate a mock catalog of halo stars whose distributions are consistent with observations of the inner MW halo (see [@Ivezic:2013aja] for a review, and references therein). The galactocentric position coordinates of the stars are assumed to follow a spherically-symmetric power-law with [@Yanny:2000ty; @Morrison:2000gp; @Bell:2007ts] $$\rho \propto r^{-3.5} \, .
\label{eq:rhoexpected}$$ The galactocentric spherical velocity distribution is modeled as a multivariate normal distribution with mean and covariance matrix $$\mu_{r,\theta,\phi} = \left( 
\begin{array}{c}
0 \\
0 \\
0 
\end{array} \right) \quad \quad
\Sigma_{r,\theta,\phi} = \left( 
\begin{array}{c c c}
120^2   & 0             & 0\\
0           & 100^2 & 0\\
0           & 0         & 100^2
\end{array} \right)  \, 
\label{eq:velexpected}$$ (with velocity measured in units of km/s) [@SommerLarsen:1996bn; @Sirko:2003uq; @Sirko:2003up; @Xue:2008se]. The mock catalog is constructed by selecting a star’s galactocentric position and velocity coordinates from these distributions. The density and velocity distributions in (\[eq:rhoexpected\]) and (\[eq:velexpected\]) match those used in [@Schlaufman:2009jv].

![The metallicity distribution for VL2 stars with . The inset shows the separately normalized distributions for the (dashed blue) and (dotted red) populations.](metdistributionHist_new.pdf){width="3.25in"}

\[fig:velmetdist\]

The top right panel of Fig. \[fig:velocitypanel\] shows the galactocentric speed distribution for the mock star catalog. Neither the ${\left[\text{Fe/H}\right]}>-1.8$ nor the ${\left[\text{Fe/H}\right]}<-1.8$ distributions resemble the smooth-halo expectation. The more metal-rich population is significantly different, as it exhibits much higher speeds than what is expected for the smooth halo. The metal-poor population bears more of a resemblance to the smooth-halo expectation, with lower speeds and a larger dispersion than its metal-rich VL2 counterpart, however it too exhibits a peaked feature indicative of kinematic substructure. Note that the differences between the smooth halo and VL2 distributions cannot be accounted for by small variations in the parameters of (\[eq:rhoexpected\]) and (\[eq:velexpected\]).

![image](DMstarcomp.pdf){width="7in"}

\[fig:velocitypanel\]

The VL2 stellar halo has a high-metallicity, high-speed population of stars that is not consistent with observations of the smooth inner halo of the MW. While the kinematic substructure is clear, there is no associated structure in position-space. To illustrate this, we randomly select 100 spheres of radius and centered at . Figure \[fig:vsampling\] shows the mean galactocentric speed distributions (dotted line), as well as the $\pm1\sigma$ spread (shaded band) of these samples. There is some variation in the speed distributions, however the peak at (for the ${\left[\text{Fe/H}\right]}<-1.8$ population) and $\sim$330 km/s (for the ${\left[\text{Fe/H}\right]}>-1.8$ population) are present over all the sampled regions. This shows that this kinematic substructure is not associated with any spatial substructure, and is thus indicative of debris flow rather than a stream.

As a point of comparison, notice that the metal-poor population has another, smaller, peak at . This peak is not present in all the sampled regions; at these speeds, a vanishing distribution is consistent to within one standard deviation. As it turns out, the stars within this peak are associated with the stream that is visible in the inset of Fig. \[fig:vsampling\] (left panel). When we look at the positions associated with the other kinematic structure in the VL2 stellar halo, there is no clear position-dependence as in this case.

![image](densityplots.pdf){width="6.5in"}

The speed distribution of stars with ${\left[\text{Fe/H}\right]}>-1.8$ shifts to lower values as one moves to larger galactocentric distances, as illustrated in Fig. \[fig:velocitypanel\]. Stars with have a mean speed of , while those with have a mean speed of . The decrease in speed is set by energy conservation, as we now explain.

When a subhalo falls into the MW, it makes numerous orbits, losing mass from tidal stripping. A rough estimate of the speed of the debris flow can be obtained from the energy conservation requirement: $$v^2\left(8.5\text{ kpc}\right) - v^2\left(D_\text{apo}^f\right) = 2 \left[\Phi\left(D_\text{apo}^f\right) - \Phi\left(8.5 \text{ kpc}\right) \right] \, ,
\label{eq:energyconservation}$$ where $\Phi$ is the gravitational potential of the host halo and $D_\text{apo}^f$ is the distance of the final apocenter. In [@Kuhlen:2012fz], it was found that the five most representative subhalos in the VL2 DM debris flow had a mean final apocenter distance of and speed of . We will assume these values here for the stellar debris flow.

![image](vlosvTHist.pdf){width="6in"}

\[fig:vlosvT\]

Taking the NFW-like profile that best-fits the VL2 host halo [@Diemand:2009bm], we find that , , and . Therefore, if the tidal debris is dominated by a single satellite (or several satellites with similar properties at final apocenter), then its speed is primarily determined by the transfer of gravitational potential energy into kinetic energy as the satellite falls towards the Galactic center (GC). This simple and intuitive picture reproduces the behavior of the high-metallicity population in Fig. \[fig:velocitypanel\] very well.

![image](EchosNew.pdf){width="7in"}

\[fig:echos1\]

The breakdown of the galactocentric velocity into its radial and tangential components depends on the orbital path of the infalling subhalo(s). The second and third rows in Fig. \[fig:velocitypanel\] show the radial and tangential-velocity distributions, respectively, for the ${\left[\text{Fe/H}\right]}<-1.8$ population (left panel), ${\left[\text{Fe/H}\right]}>-1.8$ population (middle panel), and the smooth halo (right panel). The kinematic substructure associated with the metal-rich population is primarily radial for and becomes primarily tangential for .

Debris flow may be thought of as a collection of overlapping streams [@Helmi:1999ks]. In an infinitesimal phase-space volume, one will observe $f_i(\textbf{x}, \textbf{v}) \, d^3\textbf{x} \, d^3\textbf{v}$ stars associated with stream $i$. Because the stars in each stream are collisionless, their phase-space density is conserved as a function of time. Therefore, as the stars in a given stream spread-out in position, they become more coherent in velocity. However, the constituent stars share a common speed due to energy conservation. Therefore, as each individual stream becomes more coherent in velocity, the dispersion in the average speed of its stars decreases. The net result is that the dispersion in the stars’ speeds (over all streams in the infinitesimal volume) also decreases. The velocity direction of the stars varies from stream to stream, so the total debris (from all overlapping streams) is not coherent in total velocity, even though it is coherent in speed.

This simple picture leads to a concrete prediction for the observation of debris flow. If the mass density of the debris from one or more subhalos of a given mass and infall redshift dominates the MW’s inner halo, then one expects to observe stars that are spread out over large regions of sky (*i.e.,* no spatial coherence), but share a common speed and metallicity. Whether the kinematic substructure is evident in the stars’ radial or tangential-velocity components depends on the properties of the satellites’ orbits.

Observational Evidence {#sec: ECHOS}
======================

Thus far, we have studied the properties of the VL2 debris flow in the Galactic frame. To compare with observational studies, however, we need the distributions in the heliocentric frame. To boost into this frame, we apply the following transformation: $$\left(\tilde{v}_x, \tilde{v}_y, \tilde{v}_z\right) = (v_x - v_x^\odot, \, v_y - v_y^\odot, \, v_z - v_z^\odot) \, ,$$ where $\left( v_x^\odot, v_y^\odot, v_z^\odot \right) = (-11.1, -232, 7.25) \text{ km/s}$ [@Schoenrich:2009bx]. (We denote all heliocentric velocities with a tilde.) It follows that the radial, $\tilde{v}_\text{r}$, and tangential, $\tilde{v}_\text{T}$, velocity components in the heliocentric frame are $$\begin{aligned}
\tilde{v}_\text{r}&=& \tilde{v}_z \sin b - ( \tilde{v}_x \cos l + \tilde{v}_y \sin l) \cos b \\ \nonumber
\tilde{v}_\text{T} &=& \sqrt{ \tilde{v}_\text{b}^2 + \tilde{v}_\text{l}^2} \,,\end{aligned}$$ where $$\begin{aligned}
\tilde{v}_\text{b} &=& \tilde{v}_z \cos b + (\tilde{v}_x \cos l + \tilde{v}_y \sin l) \sin b\\ \nonumber
\tilde{v}_\text{l} &=& \tilde{v}_x \sin l - \tilde{v}_y \cos l \, ,\end{aligned}$$ and $l$ and $b$ are the Galactic longitude and latitude, respectively.

Figure \[fig:vlosvT\] shows the radial and tangential velocities of the (solid green) and (dotted orange) populations in the VL2 inner halo, as defined by (\[eq:innerhalo\]). Note that the orientation of the Galactic plane is arbitrary in VL2 because there is no baryonic disk, so the results for two disk orientations are presented.

The velocity substructure persists in the heliocentric frame. For Disk Orientation A, which is perpendicular to the z-axis, the radial velocity is peaked at . This is quite different from the smooth-halo expectation (solid black line), which is distributed about zero. In addition, the more metal-rich sample is skewed towards higher tangential velocities than its smooth-halo counterpart. The kinematic substructure is primarily radial for Disk Orientation A. This changes, however, with the orientation of the disk plane. For example, when the disk is rotated by 90$^\circ$ (Disk Orientation B), the radial-velocity structure is still present, but suppressed, and the high-$\tilde{v}_\text{T}$ tail of the metal-rich distribution is more pronounced. In this case, more stars have smaller radial velocities (explaining the peak at ), but larger tangential velocities.

The peaks in the radial-velocity distributions of Fig. \[fig:vlosvT\] are symmetric about zero, indicating that an equivalent amount of debris is being stripped from the satellites as they move towards and away from the GC. One should keep in mind that this symmetry may be affected by the presence of the Galactic plane (not accounted for in VL2), which will affect the orbit of the satellites. A complete characterization of this behavior requires studying debris flow in an $N$-body simulation that contains baryons.

Radial-Velocity Substructure in SEGUE {#sub:RVsubstructure}
-------------------------------------

[@Schlaufman:2009jv] claim to detect radial-velocity substructure in the inner MW halo using a SEGUE catalog of 43,000 metal-poor main-sequence turnoff (MPMSTO) stars, 10,739 of which are located in the inner halo. In general, MPMSTO stars are good tracers of the properties of the inner halo because they provide sufficient number density and luminosity. SEGUE has 137 lines-of-sight, each designated by a unique latitude and longitude. Each line-of-sight corresponds to a $\sim$$7^{\circ}$-square patch of sky and is separated from its nearest neighbor by $\sim$10$^{\circ}$–20$^{\circ}$. As a result, one can assume that the lines-of-sight are independent of one another.

[@Schlaufman:2009jv] obtained the radial-velocity distribution along each line-of-sight, and then compared it to the expected distribution for a smooth halo along that same line-of-sight. The smooth halo was modeled using a mock star catalog that satisfied empirical measurements of the inner halo—the density and velocity distributions of the mock stars were sampled from distributions like (\[eq:rhoexpected\]) and (\[eq:velexpected\]). Significant deviations in the measured line-of-sight velocities from those in the mock catalog point to the presence of kinematic substructure. For a complete discussion of the statistical approaches used to identify deviations from the smooth-halo expectations, we refer the reader to [@Schlaufman:2009jv].

[@Schlaufman:2009jv] report 10 high-confidence detections, in which no false-positives are expected. Relaxing their statistical requirements, the authors also compile a list of 21 lower-confidence detections, of which less than three may be false positives. The algorithm used to aggregate this collection of lower-confidence detections is more adept at picking out diffuse substructure along a line-of-sight, but at the expense of a greater false-positive rate. Three each of the high-confidence and lower-confidence detections are associated with known streams [@Schlaufman:2011kf]. However, the remaining detections appear to be distributed roughly isotropically along all the lines-of-sight and do not exhibit a stream-like morphology.

In a follow-up study, [@Schlaufman:2011kf] found that these detections are chemically distinct from the smooth halo, strongly suggesting that they are the tidal remnants of a merger event.[^3] Figure \[fig:echos1\] shows the radial velocity of the high- and lower-confidence detections, versus their metallicity. The detections with the largest number densities are clustered around ${\left[\text{Fe/H}\right]}\sim-1$. In general, the radial-velocity substructure detections in the SEGUE data are more metal-rich than typical halo stars, but more metal-poor than typical thick and thin-disk stars. The kinematic and chemical properties of these stars are indicative of the tidal disruption of a merging satellite. It is not clear whether the observed substructure is due to one or more merging satellites, however. If it is only due to one merging event, then [@Schlaufman:2011kf] estimate that the progenitor mass is at least and that the merger took place around $z\lesssim0.5$, based on the metallicity and velocity dispersion of the detections.

![image](EchosHist.pdf){width="7in"}

\[fig:echoshist\]

To be consistent with debris flow, the substructure must be characterized by a coherent speed. Only the heliocentric radial velocities of these detections was observed, so we start there to look for any patterns. Figure \[fig:echoshist\] shows the distribution of the heliocentric radial velocities for the high- and lower-confidence detections (left and middle panels, respectively). The expectation for the smooth halo is shown in the right panel (solid black), along with the distributions for stars with ${\left[\text{Fe/H}\right]}>-1$ in the VL2 inner halo with Disk Orientation A (dashed green) and B (dotted orange). Due to the selection function, the SEGUE detections are clustered around $130^\circ \leq l \leq 225^\circ$ and $|b| \leq 66^\circ$, so we add this additional requirement to the stars in the smooth halo and VL2 sample. There are hints that the distribution of the lower-confidence detections is bimodal. If this is the case, then the SEGUE detections would be more consistent with the VL2 distributions, than with the smooth-halo expectation.

As discussed earlier, the constituent stars of stellar debris flow share a common speed, set by the energetics of the satellite orbit. The corresponding radial and tangential components of the stellar constituents depend on the orbital path of the satellite, but are always consistent with this speed. The current data does suggest that more substructure might indeed be present along the lines-of-sight. Figure \[fig:echosmet1\] shows the metallicity of the stars labeled as being part of radial-velocity substructure and the metallicity of those that are not (along the same lines-of-sight), reproduced from [@Schlaufman:2011kf]. One cannot distinguish the latter set of stars from the smooth halo using the SEGUE data alone. However, notice that the lines-of-sight that have a significant fraction of stars in radial-velocity substructure also have higher metallicities for all the other stars. This correlation suggests that these other stars may not actually be part of the smooth halo, but are themselves in tangential-velocity substructure. Our study of stellar debris flow in the VL2 simulated halo suggests that this should indeed be the case.

![The metallicity of the high-confidence detections in [@Schlaufman:2009jv; @Schlaufman:2011kf], compared to the metallicity of the stars along the same line-of-sight that are not identified as being part of radial-velocity substructure (filled circles). The color of the points indicates the $n_s/N_s$ along that line-of-sight, where $n_s$ is the number of stars in identified radial-velocity substructure and $N_s$ is the total number of other stars. The black diamonds correspond to substructure detections that are associated with known streams. The data is taken from [@Schlaufman:2011kf]. ](EchosSmoothMet.pdf){width="3.25in"}

\[fig:echosmet1\]

While current observational data from SEGUE suggests the presence of local kinematic substructure in the stellar halo, a complete characterization requires the full three-dimensional velocities of the stars. The GAIA satellite, which launched in December 2013, will provide high-precision radial and proper motions for an unprecedented number of stars in the MW [@Perryman:2001sp; @2008IAUS..248..217L]. This data will be critical for mapping the local stellar phase-space distribution.

If debris flow is present in the Solar neighborhood, one would expect to observe a group of stars that share a common metallicity and speed, but are not localized spatially. The kinematic and chemical properties of the debris flow will depend on the subhalo(s) from which the stars were tidally stripped. The fraction of the inner halo in kinematic substructure can be substantial—as evidenced by our study of VL2. The GAIA satellite will provide the first opportunity to fully characterize this structure in detail, helping to decode information about older mergers in the MW halo and, potentially, the local DM distribution, as we now discuss.

Dark Matter Implications & Conclusions {#sec:DMImplications}
======================================

The presence of substructure in the phase-space distribution of the Galactic halo provides a fossil record of its evolutionary history. In particular, the kinematic and chemical properties of the substructure hints at the time of the merger, as well as its size. For instance, a recently accreted satellite may leave tidal debris in the form of streams, whose constituent stars exhibit structure in both position and velocity-space. In contrast, the debris from an older merger event loses coherence in position, even as it retains structure in velocity. The debris from very old mergers is fully isotropized and indistinguishable in position and velocity from the host. The metal content of the stellar debris is correlated with the mass of the merging satellite; more massive satellites lead to more metal-rich debris.

In this paper, we focus on older mergers whose tidal debris is only characterized by kinematic substructure—specifically, a coherent speed. We refer to this class of substructure as “debris flow" to distinguish it from tidal streams, which are coherent in both position and velocity. We study the properties of stellar debris flow in the VL2 halo dynamically populated with stars. The VL2 halo contains a subset of high-speed stars with metallicities ${\left[\text{Fe/H}\right]}\sim -1$ that are distinct from the smooth halo. Because the mass density distribution for this stellar subset exhibits no spatial features, it is an example of debris flow rather than a stream.

Having studied the stellar debris flow in VL2, we can now ask whether the DM in the simulated halo exhibits similar phase-space features. A study of the DM debris flow in VL2 was completed in [@Lisanti:2011as] and [@Kuhlen:2012fz]. There, evidence was presented for kinematic substructure in the VL2 DM halo with the following properties. First, it comprises $\sim$22% of all the DM within of the GC, but almost $\sim$80% of the DM with speeds greater than $\sim$450 km/s. Second, the density profile of the debris is consistent with that of the background halo within $\sim$50 kpc of the GC. Third, the distribution of speeds is peaked $\sim$340 km/s, and the velocities are not coherent.

The speed distribution for the DM debris flow in VL2 is remarkably similar to that of its stellar counterpart. (Compare, for example, Fig. \[fig:velocitypanel\] of this work with Fig. 2 of [@Lisanti:2011as].) In addition, both the DM and stellar debris appear to have similar origin. The DM debris flow in VL2 is dominated by tidal debris from $\sim$5 satellites with infall time ${z_{\text{infall}}}\sim2$ and mass $\sim10^9$–$10^{10}$ M$_\odot$. The stellar halo in VL2 is dominated by a subhalo of similar mass that fell in at $z=2$. This suggests that the stars are indeed a good tracer for the underlying DM distribution in the dynamically populated stellar halo of VL2.

To see whether the features of the VL2 debris flow hold more generically requires repeating this analysis with other $N$-body simulations. The speed and overall density of the debris flow depends on the merger history, which varies between simulations. One would want to understand, for instance, the range of possibilities allowed for a set of realistic accretion histories. In addition, a study of DM and stellar debris flow in a simulation that contains baryons is crucial for understanding how feedback mechanisms that inject energy redistribute the DM. We expect that the tagging procedure used for VL2 should fail in regions dominated by baryons, such as near the GC or in the Galactic plane. It has been shown, for example, that such tagging prescriptions can induce systematic uncertainties in the concentration and density structure of halos [@Bailin:2014qfa]. Regarding the properties of debris flow specifically, the presence of baryons can affect the orbits of subhalos as they pass through baryon-rich areas, which could change the properties of the flow. To better understand the effect of baryons, we plan to study debris flow in Eris, a full hydrodynamic simulation that includes the effects of star formation and feedback [@Guedes:2011ux; @Kuhlen:2013tra].

Identifying the presence of stellar debris flow in the Solar neighborhood through its kinematic behavior and chemical abundance would have profound implications for DM searches. The presence of local kinematic substructure would strongly suggest that a fraction of the local halo is comprised of debris from a merging satellite. This would indicate that there are more high-speed DM particles in the solar neighborhood than would be expected for a fully equilibrated halo. It can explain, for instance, deviations from the Maxwell-Boltzmann distribution along the high-velocity tail, as observed in numerical simulations [@Fairbairn:2008gz; @Vogelsberger:2008qb; @Kuhlen:2009vh; @Lisanti:2010qx]. These deviations can affect the annual modulation spectrum at direct-detection experiments [@Freese:2012xd], as discussed in detail in [@Kuhlen:2012fz]. In particular, debris flow may result in a larger modulation amplitude at high recoil energies (than would be expected for a Maxwell-Boltzmann distribution), especially for DM with large minimum scattering thresholds.

In Sec. \[sec: ECHOS\], we discussed evidence from SEGUE for radial-velocity substructure in the stellar halo. The stars that are identified as being part of this kinematic substructure have distinctive metallicities (${\left[\text{Fe/H}\right]}\sim-1$) from the smooth halo and are isotropically distributed. The morphology of the detections is consistent with debris flow. A more complete understanding would require knowing the full proper motions of the stars, which will be provided by the upcoming GAIA satellite. If the SEGUE results are confirmed, then they would provide evidence that a significant fraction of the local halo is not in equilibrium. This would strongly suggest that, like the stellar debris in the halo, the local DM is in kinematic substructure as well.

Acknowledgments {#acknowledgments .unnumbered}
===============

We thank M. Kuhlen and V. Rashkov for collaboration at the start of this work and for providing the stellar catalog. We also thank J. Bovy, S. Gardner, A. Kepley, K. Johnston, H. Morrison, C. Rockosi, and M. Valluri for helpful discussions. This research was supported in part by the NSF under grants PHY11-25915 and PHY-1066293. Support for this work was also provided by the NSF through grants OIA-1124453 and AST-1229745 and by NASA through grant NNX12A587G(P.M.).\

[^1]: Just as one DM “particle" in the simulation does not actually correspond to an elementary particle, one star “particle" corresponds to a conglomerate of stars with total mass ${\text{m}_{\text{sp}}}$.

[^2]: In the right panel of Fig. \[fig:scatter\], there is a tail of high-metallicity stars that extends to low speeds. This contribution to the high-metallicity sample is subdominant and was accreted around $z_\text{infall} \sim 8$.

[^3]: The ECHOS are also less $\alpha$-enhanced than typical smooth-halo stars. We do not focus on this here because the measured errors on \[$\alpha$/Fe\] are typically larger than those for ${\left[\text{Fe/H}\right]}$. In addition, the VL2 stellar tagging prescription does not include information on the \[$\alpha$/Fe\] of the stars.
