arxiv: v1 [astro-ph] 1 Oct 2007

Transcription

1 Draft version August 4, 28 Preprint typeset using L A TEX style emulateapj v. 1/9/6 SIGNATURES OF ΛCDM SUBSTRUCTURE IN TIDAL DEBRIS Jennifer M. Siegal-Gaskins 1,2 and Monica Valluri 1,3,4 Draft version August 4, 28 arxiv:71.385v1 [astro-ph] 1 Oct 27 ABSTRACT In the past decade, surveys of the stellar component of the Galaxy such as SDSS and 2MASS have revealed a number of stellar streams, the results of tidal disruptions of dwarf galaxies and globular clusters by the gravitational potential of the Milky Way. Current and future observations of the stellar component of the galaxy are rapidly increasing the precision and quantity of data available, raising the possibility of using tidal streams to constrain the distribution of dark matter in the halo. Simulations of hierarchical structure formation in ΛCDM cosmologies predict that the dark matter halo of a galaxy like the Milky Way contains a smooth component as well as hundreds of subhalos with masses of 1 8 M and greater, and it has been suggested that the existence of coherent tidal streams is incompatible with the expected abundance of substructure. In this study we investigate the properties of tidal streams arising from the disruption of satellites in a variety of dark matter halo models. We select a wide range of orbits in models with both spherical and non-spherical halo shapes, and simulate the disruption of a self-gravitating satellite on these orbits with and without the presence of substructure. In general, we find that the halo shape and the specific orbital path more strongly determine the degree of disruption of the satellite than does the presence or absence of substructure, and that in some scenarios the influence of substructure on the resulting tidal debris is negligible. The presence of substructure typically leads to an increase in the degree of clumpiness of the tidal debris in a sky projection, and differences in the location of the debris compared to the results of the smooth halo model, the latter of which may have important implications for the interpretation and modeling of observed tidal streams. In addition, we identify a unique signature of the presence of substructure in the halo, which should be detectable by current and upcoming surveys. In contrast with the results of previous work, we find that a detection of a coherent tidal stream would not be inconsistent with predicted levels of substructure. Subject headings: cosmology: theory dark matter galaxies: structure methods: N-body simulations 1. INTRODUCTION Recent surveys of the stellar component of the Milky Way have revealed rich structure and a number of new tidal streams (Belokurov et al. 26b; Grillmair 26; Grillmair & Dionatos 26b; Belokurov et al. 27a). In addition to the well-studied Sagittarius (Sgr) stream (Ibata et al. 21a; Majewski et al. 23; Newberg et al. 23; Martínez-Delgado et al. 24) and the Monoceros stream (Newberg et al. 22; Yanny et al. 23; Peñarrubia et al. 25), streams from the tidal disruption of globular clusters such as Palomar 5 (Odenkirchen et al. 21; Rockosi et al. 22; Odenkirchen et al. 23; Grillmair & Dionatos 26a) and NGC 5466 (Grillmair & Johnson 26; Belokurov et al. 26a) have also been identified. With the large increase in data expected from current and upcoming missions such as SDSS-SEGUE (York et al. 2; Newberg et al. 26), RAVE (Steinmetz et al. 26), GAIA (Perryman et al. 21), and SIM-Planetquest (Majewski et al. 26), the detailed properties of the 1 Kavli Institute for Cosmological Physics 2 Department of Physics, University of Chicago, 564 S. Ellis Avenue, Chicago, IL 6637 jsg@uchicago.edu 3 Department of Astronomy and Astrophysics, University of Chicago, 564 S. Ellis Avenue, Chicago, IL 6637 valluri@kicp.uchicago.edu 4 Department of Astronomy, University of Michigan, Church Street, Ann Arbor, MI 4814 stellar component of the Milky Way will soon be measured with unprecedented accuracy. Tidal streams are emerging as a valuable tool for probing the mass distribution of the Galaxy, since the leading and trailing arms of a tidally disrupting satellite trace out its orbital path, which is sensitive to the shape of the gravitational potential. ΛCDM simulations of structure formation find that galaxy-sized dark matter halos are significantly aspherical (see, e.g., Frenk et al. 1988), although the detailed shape of dark matter halos is still an open question. For example, Kazantzidis et al. (24a) argue that simulations incorporating gas cooling lead to more spherical halos than simulations of dark matter alone, and a recent study by Hayashi et al. (27) finds that the triaxiality of the halo is radially dependent, with the inner regions tending to be prolate, and the outer regions more spherical. Since tidal streams probe the gravitational potential over a large range of radii, a number of studies have examined whether observations of tidal streams could be used to constrain the shape of the dark matter halo of the Galaxy, but have in many cases reached different conclusions. For example, Helmi (24) uses radial velocities for a large sample of M giant stars from the Sgr stream to determine that the halo must be quite prolate, while Ibata et al. (21b) conclude from similar data for carbon stars from the Sgr stream that the halo is nearly spherical, and Johnston et al. (25) argue that the Sgr data favors oblate models, ruling out prolate

2 2 models at high confidence. Peñarrubia et al. (25) find that Monoceros stream data is also best fit by an oblate halo. More recently, Fellhauer et al. (26) argue that the bifurcation of the Sgr stream can only be explained by a nearly spherical halo. Law et al. (25) suggest that these conflicting results are evidence of evolution of the Galactic potential. In addition to the uncertainties in the shape of the dark matter halo, the missing satellites problem in ΛCDM cosmologies may present other issues for the interpretation of tidal stream data. On the scale of galaxies like the Milky Way, simulations of structure formation predict as many as several hundred dense clumps with masses 1 8 M in the halo, roughly an order of magnitude more than the observed abundance of dwarf galaxies surrounding the Milky Way (Klypin et al. 1999; Moore et al. 1999). One class of solutions to this problem involves modifications of the properties of the dark matter particle, e.g., warm dark matter (Hogan & Dalcanton 2), or self-interacting dark matter (Spergel & Steinhardt 2), in order to reduce the amount of substructure. Some authors have proposed that suppression in the small-scale end of the primordial power spectrum may reduce the predicted abundance of subhalos (Kamionkowski & Liddle 2; Zentner & Bullock 23). This class of solutions predicts a reduction in the abundance of dark substructure in the halo. An alternative class of solutions to the missing satellites problem suggests mechanisms by which only a fraction of the dark matter subhalos have visible stellar components at the present epoch (Bullock et al. 21b; Kravtsov et al. 24). In these scenarios, a large number of dark substructures are predicted to populate the halo. It is worth noting that in the past few years, the number of known satellites of the Milky Way has increased dramatically, with the discovery of 11 new, faint objects which are likely dwarf galaxies or globular clusters (Willman et al. 25a,b; Zucker et al. 26a,b; Belokurov et al. 26c, 27b; Irwin et al. 27). This rapid discovery of numerous previously unknown satellites of the Milky Way suggests the possibility that a population of very faint, still undiscovered satellites could populate the halo. If such a population of dark or faint subhalos does exist, these structures may be detectable in a number of ways. Gamma rays from the annihilation of dark matter particles in substructures may reveal their presence (e.g., Bergström et al. 1999; Calcáneo-Roldán & Moore 2; Tasitsiomi & Olinto 22; Stoehr et al. 23; Koushiappas et al. 24; Kuhlen et al. 27). Dark substructures may also be detectable via flux anomalies in multi-image gravitational lens systems (Mao & Schneider 1998; Metcalf & Madau 21; Dalal & Kochanek 22; Chiba 22) or via spectroscopic lensing (Moustakas & Metcalf 23), and although such effects have been observed, it has not yet been demonstrated indisputably that substructure is the unique cause. In this study we examine whether substructure may be detectable by its influence on the formation and properties of tidal streams. In particular, we aim to determine the prevalence, characteristics, and detectability of the effects of a population of subhalos similar to that predicted by ΛCDM models on tidal streams of dwarf galaxies. To this end, we consider two additional variables which play a major role in determining the properties of the streams produced: the shape of the dark matter halo and the orbital path of the progenitor. Tidal debris can be dispersed both by precession of the orbit of the progenitor due to a non-spherical potential and by heating from interactions with substructure (Ibata et al. 22; Mayer et al. 22), so it is necessary to evaluate the relative importance of substructure in a non-spherical halo. With regard to the influence of the orbit of the progenitor, since the disruption of a satellite is governed by the tidal forces it experiences along its orbit, the same progenitor satellite can produce very different tidal streams if placed along a different orbital path. Furthermore, the strength and frequency of interactions with a population of substructures will depend on the orbital path of the progenitor in relation to the orbits of the subhalos. To understand the expected range of types of tidal debris arising solely from differences in orbital paths, we select orbits with a variety of characteristics. The goals of this paper are to determine which, if any, observable properties of tidal debris can be used to unambiguously detect or rule out the presence of substructure in the halo when a wide range of allowed orbits and host mass distributions are considered. Several previous studies have investigated the impact a population of substructures would have on tidal streams, using different techniques to study some specific scenarios in detail (Ibata et al. 22; Johnston et al. 22; Mayer et al. 22; Peñarrubia et al. 26), and we discuss our results in relation to the conclusions of those studies in 6. In 2 we discuss the phase space properties of orbits in the context of tidal disruption, describe our host galaxy models, and outline our method for selecting a diverse range of orbits. We give the details of our numerical simulations in 3, and present our results in 4. In 5 we address observational considerations. We evaluate the implications of our work and discuss its relation to previous work on this subject in 6, and finally summarize our conclusions in ORBITS 2.1. Phase space properties of orbits As part of our effort to examine a large number of factors which may influence the formation and characteristics of tidal streams, we consider the role played by the phase space properties of a progenitor s orbit on its tidal disruption. The structure of phase space in Hamiltonian dynamics is a well studied subject (see, e.g., Binney & Tremaine 1987; Lichtenberg & Lieberman 1992). In this study we use some basic principles governing the motion of orbits in phase space to select a range of orbits which may lead to streams with diverse properties. In a spherically symmetric potential, orbits conserve four integrals of motion (the energy, as well as each of the three components of angular momentum), so orbits in such potentials are generically rosettes confined to a single plane in configuration space. In a smooth non-spherical potential, such as that of a galaxy with a disk and possibly a non-spherical halo, a large fraction of orbits are expected to be regular, i.e., they conserve three integrals of motion, such as the energy and two quantities resembling angular momentum,

3 3 or chaotic, in which case they conserve fewer integrals of motion than spatial dimensions. Consequently, orbits are in general no longer confined to a plane, but will densely fill 3-dimensional volumes in configuration space. A resonant orbit is a regular orbit whose fundamental oscillation frequencies in configuration space are related via a linear equation such as lω 1 + mω 2 + nω 3 =, where l, m, n are integers. Such a relation between the fundamental frequencies of motion implies that only two of the three frequencies are independent, with the third frequency a linear combination of the first two. This additional condition on the orbital motion is similar to requiring that a fourth integral of motion be conserved. Consequently, a resonant orbit satisfying one or two resonance equations as above will be confined to a sheet or a closed loop, respectively, in configuration space (Valluri & Merritt 1998; Merritt & Valluri 1999). In most potentials that are near-integrable, the vast majority of orbits are regular orbits which satisfy three integrals of motion. Only a tiny fraction of phase space is expected to be occupied by chaotic orbits and resonant orbits (Lichtenberg & Lieberman 1992). Time-independent perturbations to a spherical potential such as that introduced by a disk or a non-spherical halo have two important effects on the structure of phase space. First, the volume of phase space occupied by resonant islands increases significantly, so that a larger fraction of orbits become associated with resonances. Second, as the strength of the perturbation grows, a larger fraction of the orbits in the potential become chaotic (Lichtenberg & Lieberman 1992). That is to say, the fractions of both chaotic orbits and orbits associated with resonances grow at the expense of the regular orbits when the perturbation strength is small. We note that while most treatments of Hamiltonian dynamical systems assume that the potential is time-independent, the addition of live subhalos, a timedependent perturbation, could be considered pseudorandom noise. Sideris & Kandrup (24) find that the addition of pseudo-random noise has effects on the structure of phase space similar to those induced by timeindependent perturbations as noted above. In particular, they find that pseudo-random noise can trigger parametric resonances that make a significant fraction of stars execute chaotic orbits, but that some fraction of resonant orbits typically survive in the presence of small amounts of noise. In the context of the Milky Way, tidal debris from satellites on resonant orbits would be confined to sheets or loops in configuration space rather than filling 3- dimensional volumes. Furthermore if the satellite is associated with a stable resonant orbit, the tidal debris along the orbit may be less susceptible to perturbations by substructure than streams formed along non-resonant orbits. The properties of tidal debris on chaotic orbits and regular orbits in a non-spherical potential are expected to be similar since both types of orbits fill large volumes of configuration space, so in this study we consider only the distinction between non-resonant (chaotic and regular) orbits and resonant orbits Orbit selection To explore the effect of tidal disruption on a wide variety of orbital paths, we construct libraries of orbits using a method similar to that of Cretton et al. (1999). Each library consists of orbits with a fixed energy E and z- component of angular momentum L z. Fixing these two quantities in a given potential determines a set of coordinates which are turning points (where v R = v z = ) for orbits, termed the zero velocity curve (ZVC). A library is generated by initializing orbits from various positions on the ZVC, and semi-analytically integrating these orbits in the host potential. Our host potential consists of a smooth NFW dark matter halo (Navarro et al. 1997), a Miyamoto- Nagai disk (Miyamoto & Nagai 1975), and a Hernquist bulge (Hernquist 199). We adopt the masses of these components and halo concentration from the favored model of Klypin et al. (22) which gives the halo mass M halo = 1 12 M and concentration c = 12, the disk mass M d = M, and the bulge mass (m 1 + m 2 in Klypin et al. (22)) M b = M. We set the scale lengths of the Miyamoto-Nagai disk to be a = 6.5 kpc and b =.26 kpc, and the Hernquist bulge scale parameter a =.7 kpc (Johnston et al. 1996). Although Klypin et al. (22) use exponential density profiles for the baryonic components of the Galaxy, since these models are considerably more computationally expensive to implement in our simulations, we instead use the Miyamoto-Nagai disk and Hernquist bulge to describe the baryon distribution. These parameters produce a rotation curve which is qualitatively similar, but not identical, to that of Klypin et al. (22). However, we adopt the masses for all of the components from Klypin et al. (22) to ensure consistency in the outer regions of the halo, where tidal debris is often located. We take the NFW halo to be axisymmetric and define the lengths of the principal axes to be (x : y : z) = (a : b : c) with the x and y axes equal. The z-coordinate is assumed to be aligned with the axis of rotation of the Galaxy. We vary the axis ratio q = c/a = c/b of the halo to generate orbit libraries for spherical (q = 1), prolate (q > 1), and oblate (q < 1) halo shapes. We select q values that represent the extremes of the shape distribution: the oblate halo has q =.6 and the prolate halo as q = 1.3. These are not intended to be necessarily representative of real halos but to allow us to explore a wide range of orbit types. The potentials for the oblate and prolate versions of the NFW halo are obtained by defining a variable r which replaces r in the spherically-symmetric NFW potential Φ NFW (r), where r 2 = x 2 + y 2 + (z/q) 2. We set the energy of an orbit by equating it to the energy of a circular orbit in the plane of the disk at a radius of R circ. In this work we focus on orbits with the energy of a circular orbit with R circ = 6 kpc which, for our models, produces debris located at radial distances similar to those at which tidal streams have been observed. We choose L z =.5L z,max, which produces a range of orbits with diverse properties, and integrate orbits in each library. We select orbits close to phase space resonances as well as non-resonant orbits to evaluate the impact that the phase space properties of an orbit have on its disruption. Some properties of the orbits we have selected to simulate are given in Table 1, with orbits closely associated with resonances labeled R in the surface of section (SOS) column, and non-resonant orbits labeled NR. Also given

4 4 TABLE 1 Orbit Parameters Halo Orbit r apo r peri SOS q (kpc) (kpc).6 O NR O R O NR O R 1. S NR S NR S R 1.3 P NR P NR P R P R y (kpc) y (kpc) y (kpc) S1 S3 O2 S1 S3 O2 S r (kpc) S r (kpc) O r (kpc) v R,sos (km/s) v R,sos (km/s) v R,sos (km/s) 3 2 S R sos (kpc) S R sos (kpc) R sos (kpc) O2 y (kpc) O3 O3 O r (kpc) v R,sos (km/s) O R sos (kpc) y (kpc) P1 P1 P r (kpc) v R,sos (km/s) 3 2 P R sos (kpc) Fig. 1. Surface of section for orbits in a model with an oblate dark matter halo (q =.6). All orbits have energy equivalent to R circ = 6 kpc and L z =.5L z,max and were integrated for at least orbital periods. A large family of orbits associated with a resonance is seen in a resonant island, which appears as concentric ellipses centered on (R sos, v R,sos ) (45, ), corresponding to the resonant orbit O2. are the apocenter (r apo ) and pericenter (r peri ) distances of each orbit. We identify orbits in each library associated with phase space resonances for each host potential using surface of section plots (see, e.g., Binney & Tremaine 1987). As an example, a compilation of surfaces of section for all orbits integrated from the oblate halo model library is shown in Figure 1. Since resonant orbits conserve an additional integral of motion, their path collapses to a single point (or a small number of points, for higher-order resonances) in a surface of section. The two concentric ellipses centered at v R,sos = in the surfaces of section for this library are generated by orbits associated with a resonant orbit at the center of these rings (orbit O4). This halo model also admits another resonant orbit family at roughly R sos = 46 kpc, v R,sos = km/s (orbit O2). The non-resonant orbits in Table 1 are selected to represent orbits not strongly associated with any resonant family. Projections of the orbital path and surface of section for some individual orbits are shown in Figure 2. For all orbits identified as resonant we have chosen the orbit in our library corresponding most closely to the resonance based on its surface of section. While our method does not exactly identify any or all resonances in a given y (kpc) P3 P3 P r (kpc) v R,sos (km/s) P R sos (kpc) Fig. 2. Projections in x-y, x-z, r-z, and surface of section R sosv R,sos for selected orbits from each halo model. S1 is a regular rosette orbit and appears as a nicely filled curve which covers a large range of radii and velocities in the SOS plot. Orbit S3 is close to a resonance and occupies a much smaller region in the SOS plot (note that the axes have been rescaled). O2 is extremely close to the strong resonance seen in the upper part of Figure 1, and O3 is an orbit that lies at the edge of the resonant island. P1, a non-resonant orbit, fills a large area of the SOS plot. Orbit P3, on the other hand, is close to a resonance and appears as three points in the SOS plot. potential, we expect that the orbits we select as representative of resonant orbits are sufficiently close to phase space resonances that they will exhibit the resonant behavior discussed in 2.1. As time-dependent perturbations such as subhalos are introduced, the surfaces of section are altered and a legitimate concern is whether any of the resonances we initially identify are still present in the perturbed system. However, as mentioned in 2.1, previous work has shown that perturbations like substructure do not destroy the strongest resonances, therefore we feel confident that even though we choose orbital initial conditions from surfaces of section in static potentials without substructure, when we add substructure to the potentials, the resonant orbits continue to be associated with the same resonances, although their proximity to the resonance may change. 3. NUMERICAL SIMULATIONS

5 5 We simulate the tidal disruption of our satellites on the orbits selected in 2 using the N-body simulation code GADGET-2 (Springel 25), which calculates gravitational forces using a hierarchical multipole expansion algorithm. The host galaxy model described in 2.2 is incorporated with a static potential using a modification to GADGET-2 kindly provided by B. Robertson. The progenitor satellite is modeled with particles following a NFW density profile with virial mass M vir = 1 1 M and concentration parameter c = 12. We generate the satellite initial conditions using the exact distribution function as detailed in Kazantzidis et al. (24b), and acknowledge S. Kazantzidis for providing us with code which implements this method. The gravitational softening length for the satellite particles is set to ǫ sat =.7 kpc. This initialization scheme sets up a satellite which is in equilibrium when isolated, but will not be in equilibrium when placed in the host potential. As a result, the satellite will initially lose mass at an artificially high rate due to the sudden introduction of an external tidal field, until it has relaxed into the host potential. Since we expect that satellites of the Milky Way fall into the potential of the Galaxy gradually from large distances or are formed in quasi-equilibrium with the halo, the mass loss due to instantaneously adding an external potential is unphysical. To minimize this effect, we attempt to create a satellite in quasi-equilibrium with the host potential at a large distance. We place the progenitor satellite on circular orbits in the plane of the disk at 3, 4, and kpc, and track the mass loss rate over time. For all three cases there is an initial period of rapid mass loss, but the mass loss rate stabilizes after less than.2 orbital periods. We generate a remnant satellite to use as our progenitor by selecting all particles still bound to the satellite after.2 orbital periods for the satellite on the kpc orbit. The closest distance at which we initialize any of our orbits is 54 kpc, and satellites initialized at a larger radius experience a less intense shock than those initialized closer to the Galactic center, so we choose the kpc orbit for generating the remnant to temper the most intense shock our satellites would experience. The remant satellite contains 14k particles and has a mass of M. For our simulations in halos with substructure, the subhalos are added as softened point masses. We use 253 subhalos from halo G 2 in the cosmological N-body simulation of Kravtsov et al. (24) which have masses of M, and thank A. Kravtsov for making this data available to us. We initialize the subhalos in our simulations at their final spatial coordinates in the cosmological simulation, and transform their velocities by requiring that their current position is a turning point of the orbit (v R = v z = ) and that the z-component of their angular momentum is preserved. Navarro et al. (1997) showed that density profiles of CDM halos are well-fit by a NFW profile over a wide range of masses. Furthermore, Kazantzidis et al. (24b) found that dark matter subhalos retain their steep inner density profile even when tidally stripped. Ideally, we should model the subhalos in our simulation using a NFW profile, however, for numerical simplicity we model the subhalos as spline-softened point masses. The splinesoftened acceleration is maximized at r.8ǫ, while TABLE 2 Subhalo Group Properties Mass Range Number of ǫ sub (1 1 M ) Subhalos (kpc) the NFW acceleration continues to increase at arbitrarily small radii, so softening the subhalos in this way necessarily places a limit on the maximum acceleration possible that is not present in the NFW model. This softening method also makes the subhalos fluffy, since the splinesoftened acceleration reaches a maximum at r >, and the radius of the maximum acceleration increases with increasing ǫ. We also note that the softened acceleration can exceed the NFW acceleration for some range of radii. To approximate the effect of adding substructures which internally are well-described by a NFW density profile, we soften the subhalos in our simulation by choosing the softening length ǫ sub such that the maximum acceleration due to the spline-softened subhalo is equal to the acceleration due to the NFW profile at r =.1r s (a radius well-within that shown to be stable against tidal disruption by Kazantzidis et al. 24b). While this choice increases the effective size of our subhalos, thus increasing the cross-section for interaction, it also artificially limits the acceleration due to the steep inner density profile. Ibata et al. (22) chose to soften subhalos in their work by requiring the softened acceleration to be smaller than the NFW acceleration at all radii, thereby underestimating the effect of the substructures. We choose to allow the softened acceleration to exceed the NFW acceleration for a small range of radii to better model the effect that the predicted inner slope of the NFW profile may have on interactions with tidal streams. We sort our subhalo particles by mass to form four groups with different softenings. We determine the concentration parameter c for the most massive subhalo in each group using the c(m vir ) relation determined by Bullock et al. (21a), which gives 18 c 4 for our subhalo masses, and from c and M vir determine r s, which sets the softening length as described above. Table 2 summarizes the characteristics of each group. We comment on the implications the choice of softening has on our results in 6. We do not explicitly include a stellar component in our satellite model, however, we attempt to follow the distribution of a stellar component by identifying a subset of dark matter particles as star particles. In the initial equilibrium halo, we mark star particles according to a Hernquist distribution. After we generate the remnant satellite there are 12, star particles remaining. For each selected orbit, we simulate the disruption of the remnant satellite in the host model both with and without dark matter subhalos, evolving our simulations for 4.9 Gyr. 4. RESULTS 4.1. Sky distribution of the debris

6 6 Fig. 3. Sky projection in Galactic coordinates of star particles after 4.9 Gyr for selected orbits in spherical halo models. Orbits in models with subhalos produce clumpier debris than in models with a smooth halo component. In the absence of subhalos the debris from the resonant orbit (S3) appears more coherent than the debris from the non-resonant orbit (S1). However, the additional disruption caused by the subhalos is of a similar degree for the resonant (S3) and non-resonant (S1) orbits. The S1 orbit produces a broad band of debris that appears somewhat clumpy even without the presence of subhalos. The S1 orbit with subhalos has clearly been perturbed, displacing a large amount of the debris relative to the smooth halo simulation. The projection of the debris in Galactic coordinates after 4.9 Gyr of evolution for a subset of orbits is shown in Figure 3 (spherical halo model), Figure 4 (oblate halo model), and Figure 5 (prolate halo model). The particles plotted are only those originally marked as star particles as described in 3. Comparing the simulations in all three halo models, it is immediately notable that even those without subhalos do not consistently produce thin, coherent streams. While simulations with subhalos sometimes lead to slightly more scattered debris or clumpier structures in these coordinates, the degree of additional scattering or clumping for the same orbit simulated with and without subhalos is small compared to the variation seen between different orbits and different halo models without subhalos. In Galactic coordinates, the S3 orbit without subhalos produces the most coherent stream, which is not surprising given that this orbit has r apo r peri 6 kpc. The S1 orbit, on the other hand, is a rosette which is largely confined to the plane of the disk of the Galaxy, and has a large apocenter to pericenter ratio (r apo /r peri 6). In this projection the S1 orbit in the smooth halo produces clumpy debris clustered along the Galactic plane. For both orbits in the spherical potential the presence of substructure produces some additional dispersion. Substructure also leads to noticeably clumpier debris for the S3 orbit, and results in debris displaced from the Galactic plane for orbit S1. In the oblate and prolate halo models, the resulting streams are quite dispersed compared to the S3 stream and somewhat clumpy both in the simulations with and without subhalos. The high degree of dispersion is expected from the precession of the orbit due to the extreme flattening (q =.6) or elongation (q = 1.3) of the halo. For the orbits shown in the oblate halo, O2 and O3, as well as the P1 orbit in the prolate halo, the debris actually appears less dispersed and more clustered in the simulations with subhalos. Orbit P3 produces somewhat smooth streams which are qualitatively very similar both with and without subhalos, suggesting it has not had any significant interactions with the subhalo population. In a sky projection, the simulations with substructure often produce clumpier debris than those without substructure, regardless of the shape of the dark matter halo. While this increase in clustering appears to be strongly correlated with the presence of substructure, further study is necessary to identify the mechanism responsible for the clumping, and to determine whether particles which appear clustered in a sky projection are truly clustered in configuration space or in phase space. Resonant orbits S3 and P3 produce slightly more coherent streams than the non-resonant orbits in the same halos. Both the highly resonant orbit O2 and the nonresonant orbit O3, however, produce qualitatively similar streams. We attribute this to the fact that while O3 is strictly non-resonant it is associated with the same resonance island as O2 and therefore has a similar orbital shape (as seen in Figure 2), although it lies further away from the center of the resonant island seen in the SOS plot in Figure 1. The higher degree of coherence can be expected for orbit S3, since it has r apo /r peri 1, and does not probe the inner region of the Galactic potential. While in this sample of orbits there is some evidence that resonant orbits may produce more coherent streams than non-resonant orbits, from the similar increase in clumping and dispersion of the debris in the simulations with subhalos, it does not appear that these resonant orbits are less susceptible to the effects of substructure than the non-resonant orbits. One important result of comparing the simulations

7 7 Fig. 4. Sky projection in Galactic coordinates of star particles after 4.9 Gyr for selected orbits in oblate halo models. The debris for orbits in models with subhalos appears slightly more clustered in this projection than for the models without subhalos, with most of the debris following a single path. Both the resonant orbit (O2) and the non-resonant orbit (O3) produce debris that is fairly dispersed, regardless of the presence of subhalos. The O3 simulation with subhalos produces an interesting feature at l 45, b 3 resembling a bifurcation in the stream. Fig. 5. Sky projection in Galactic coordinates of star particles after 4.9 Gyr for selected orbits in prolate halo models. The debris from orbit P1 (a non-resonant orbit) is more highly dispersed than the debris from P3 (a resonant orbit) regardless of the presence of subhalos. The addition of subhalos results in slightly more clustered debris in the case of P1, while the debris from orbit P3 qualitatively shows little difference between models with and without subhalos, clearly producing streams in both cases. with and without subhalos is that, in some cases, orbits in models with subhalos result in streams and debris in very different locations than in the smooth halo models. The S1 orbit with subhalos produces a stream that is clearly offset from the Galactic plane for l, while the stream from the same orbit in the model without substructure wraps around the Galaxy at b. Similar formations can be seen in most of the simulations with substructure, suggesting that disturbances in tidal streams due to interactions with substructures (both dark and luminous) in the halo may be an important consideration when constructing models to fit observations of tidal streams. Of particular interest is a small stream of stars slightly offset from the main stream at l 45, b 3 in the O3 simulation with subhalos which in these coordinates qualitatively resembles the bifurcation seen in the Sgr stream, which Fellhauer et al. (26) conclude could only be explained if the Galaxy s

8 8 halo is nearly spherical, but which appears in our oblate model with substructure. Similar features can be found in many of our simulations. We do not in any way suggest that our orbit O3 is an appropriate description of the orbit of Sagittarius, but mention this feature as an example of the type of effect substructure can have on tidal streams. We note that since the dispersion of the streams in the non-spherical halo models is considerable, and for a given orbit is seen both with and without the presence of subhalos, heating due to subhalos is not the dominant cause of the dispersion in our prolate and oblate halo models Structure in phase space coordinates The distribution of the star particles in line-of-sight velocity in the rest frame of the Galaxy and heliocentric distance is shown in Figures 6, 7, and 8. We project the simulated data into these coordinates to show that coherent structure is evident in the debris in observable coordinates. Ideally, we would view the debris in terms of radial velocity and distance with respect to the Galactic center, since a Galactocentric coordinate system reflects the symmetries of the potential, so a wealth of structure is typically visible (see, e.g., Bullock & Johnston 25). However, even upcoming experiments are unlikely to be able to measure proper motions of stars at very large distances, making it impossible to determine Galactocentric velocities for a subset of stars in the survey. For stars at large distances, the comparatively small offset of the sun from the Galactic center leads to small differences between the Galactocentric and heliocentric distances, however, the relative motion of the sun leads to a wide spread in the observed line-of-sight velocities for debris with similar Galactocentric radial velocities but a large range in angular position. By subtracting the line-ofsight component of the sun s velocity v sun,los, a known quantity, from each star s line-of-sight velocity v los, we recover the star s velocity along the line-of-sight relative to the Galactic rest frame v los,rest, which for stars at large distances is a good approximation to the radial velocity with respect to the Galactic center. As a result, structure is visible in these observable coordinates. It is notable that simulations with subhalos tend to produce debris which is less coherent in these coordinates than debris in their corresponding simulations in the smooth halo models. However, the variation in the patterns formed in the absence of substructure is quite large compared to the changes induced by the addition of substructure. For example, orbits S1 and S3 in the smooth halo models produce drastically different patterns: the S1 simulation produces clear bands of debris in these coordinates, while the S3 simulation appears as two compact clumps with little internal structure. Many simulations with subhalos are distinguishable in these coordinates by a conspicuous trail of debris at large radii strongly correlated in v los,rest. Furthermore, v los,rest for the stars in this feature is positive and increases with increasing radius, unlike patterns formed by the wrapping of debris in smooth halo models (e.g., in the S1 and P3 orbits without subhalos). This feature is present to some degree in most of our simulations with subhalos in all halo shapes, and is consistently absent when the satellite is evolved in a smooth halo. Specifying that the stars in this feature must have increasing positive velocity at larger radii eliminates the degeneracy between patterns formed by subhalo interactions and those resulting from the development of streams in a smooth halo as in the S1 and P3 orbits. These arms of stars typically extend far beyond the debris seen in the smooth halo simulations. Many orbits in simulations with substructure produce debris extending to 2 kpc, while debris from these orbits in smooth halo simulations is typically confined to kpc. The stars in these arms are generally not very spatially clustered, so their strong correlation in velocity may prove essential for identifying stars as originating from the same progenitor, as well as for cleanly identifying this feature as the result of substructure. We emphasize that this observable signature appears to be unique to simulations with subhalos, making a detection of such correlated stars a potential smoking gun indicating the presence of substructure in the halo, regardless of the shape of the dark matter halo. To confirm that this striking feature does not result from lack of precision of the integration scheme, we repeated the S3 simulation with smaller timesteps to prevent unphysical two-body interactions, and found the same feature present. Additionally, we repeated a number of simulations using the subhalos taken from host G3 in Kravtsov et al. (24), and again found that many of the simulations developed similar trails of stars at large radii, suggesting that the feature does not arise from a rare, serendipitous arrangement of the subhalos. To identify the mechanism leading to tidal debris that is highly coherent in velocity at large radii, we examine the acceleration produced by the subhalos on selected star particles. We estimate the maximum acceleration produced by a subhalo on a given star particle over the course of the simulation by determining the closest approach of the particle to the subhalo in any snapshot. We emphasize that due to the limited time resolution of our snapshots, the distance of closest approach determined in this way is likely to be significantly overestimated, and hence the maximum acceleration may be significantly underestimated. During periods of strong acceleration the position of the particle changes quickly, and so it is unlikely that the true closest approach of the particle to each subhalo was recorded in a snapshot. Figure 9 shows the maximum acceleration due to each subhalo as a function of subhalo mass for four selected star particles in the O3 simulation with substructure. Each of the four particles was selected from a different region of the r sun -v los,rest plot in Figure 7: a particle from the shorter, upper arm (black crosses), the tip of the longer arm (red diamonds), the middle of the longer arm (purple squares), and the lower part of the main region (blue triangles). A clear trend is evident, with the most massive subhalos typically producing the largest accelerations. A large acceleration does not guarantee a substantial energy exchange, but indicates the strength of an individual encounter. The vertical axis is scaled logarithmically to show the maximum accelerations of the lower mass subhalos. All three particles from the arms had a strong encounter with a very massive subhalo ( 1 1 M ), while the particle not located in an arm (marked with a blue triangle) did not have any encounters producing accelerations greater than 3 km/s

9 9 Fig. 6. Heliocentric line-of-sight velocity in the rest frame of the Galaxy versus heliocentric distance of star particles in spherical halo models after 4.9 Gyr for the same orbits as in Figure 3. Orbits in simulations with subhalos produce arms of debris in these coordinates at relatively large radii (r sun 1 kpc) strongly correlated in velocity. This feature is also seen for many orbits in oblate and prolate models with subhalos. The S3 simulation with subhalos produces multiple arms of debris, the longest of which extends to 2 kpc, however, the debris from the S3 simulation without subhalos scarcely extends beyond kpc, and shows little structure in these coordinates. Orbit S1 in the smooth halo shows two clear bands of stars, whereas only one fuzzy band is seen for the same orbit simulated with substructure. /Gyr recorded in a snapshot. The time evolution of the total energy of selected particles is shown in Figure 1, along with the proximity of each particle to the subhalo which caused it the greatest acceleration. The particles selected correspond to those in Figure 9. The solid line shows total energy of the particle, and the dotted line shows the subhalo distance. In the top left panel, the distance is shown for the second most massive subhalo; in the other panels, the distance to the most massive subhalo is plotted. For the three particles in the arms at the end of the simulation (the two top panels and the bottom left panel), a large increase in energy is evident, clearly corresponding to the time of closest approach of the subhalo. The energy of the particle in the bottom right panel, which had a lesser maximum acceleration in Figure 9, shows smaller deviations, but which still appear to be correlated with the proximity of the subhalo. The net change in total energy for this particle over the course of the simulation is small compared to the increases in the other particles energies. For the three particles which did gain a substantial amount of energy, the majority of the increase occurred during a single interaction with the subhalo producing the largest recorded acceleration. While the formation of an arm of stars at large radii from the disruption of a satellite would thus seem to be dependent upon the specific orbital paths of the most massive substructures and the satellite, this feature is frequently present in our simulations. We expect that the magnitude of the energy transfer is somewhat sensitive to our choice of softening, since it is likely that during these strong encounters the particle probes the softened potential in the innermost regions of the subhalo. In addition, we note that the fluffiness of the subhalos may to some degree artificially increase the frequency of interactions capable of producing such large accelerations. Are these arms of particles at large distances the result of special conditions, or do these features arise naturally from scattering interactions with subhalos? Figure 11 shows the distribution of the fractional change in total energy over the course of the simulation for all satellite particles (both dark matter and star particles) in the O3 simulation with subhalos. The solid line shows the change in energy for all particles, while the red dashed line represents only the particles in the longer arm at the end of the simulation (r sun > 17 kpc). The small positive offset from zero seen for the distribution for all particles shows that satellite particles on average gained a small amount of energy ( 8%). This suggests that interactions with the subhalo population have a slight tendency to increase the energy of the stream particles. The lower panel shows the distribution of only the particles from the arm in relation to the distribution of all the satellite particles. The arm particles appear to naturally fill the positive end of the E/ E distribution, accounting for almost all of the particles which experienced an energy increase of 45% or more. In addition, all of the 55 star particles in this feature experienced an

10 1 Fig. 7. Heliocentric line-of-sight velocity in the rest frame of the Galaxy versus heliocentric distance of star particles in oblate halo models after 4.9 Gyr for the same orbits as in Figure 4. As in the other halo models, arms of star particles at large radii are present for orbits in simulations with subhalos, with the O2 arm stretching past 2 kpc. Neither of these orbits shows significant structure in these coordinates when simulated in a smooth halo. Although O2 is a resonant orbit, its debris does not appear significantly more coherent in these coordinates than the debris from orbit O3. acceleration from a single subhalo of at least 6 km/s /Gyr, and for all but one of these particles the most massive subhalo was responsible for the strongest recorded encounter. Thus it seems that this feature generally results from interactions of satellite particles with the most massive substructures, although the details of these interactions may in fact be more complex. A number of mechanisms for transferring large amounts of energy to stream particles are possible, including two-body encounters (e.g., scattering of a particle by a subhalo), three-body encounters (e.g., simultaneous interaction of two subhalos and a stream particle), and less easily categorized interactions (e.g., a particlesubhalo encounter during which the acceleration of the subhalo by the host potential is a non-negligible contribution to the particle s energy gain). While we defer a thorough investigation of the energy exchange mechanism to a future study, we examine the processes involved for a few selected particles which experienced large gains in energy in a relatively short time. Due to the large number of influences on stream particles in the simulation at any given time, the simple interactions we consider should be taken as only approximate descriptions of the underlying dynamics. For clarity in the context of understanding the energy transfer to stream particles, we define a two-body interaction to be a simple scattering event during which the acceleration and energy gain of the particle is due primarily to the influence of a single subhalo. In contrast, we take a three-body interaction to be one in which two subhalos induce comparable accelerations to the particle during the time the majority of the energy transfer occurs. We find examples of an isolated two-body interaction between a particle and a subhalo, and of a two-body interaction during which the acceleration of the host potential appears to play an important role. While we do not find clear evidence of three-body interactions involving two subhalos, we emphasize that our inquiry has been very restricted, and thus in no way rules out the possibility of such interactions. Our simulations follow the disruption of a single satellite, and in most cases we find that some fraction of stars are found at large radii due to interactions with the more massive substructures. In general, the presence and strength of this feature is dependent on the volume filled by debris stripped from the satellite, which varies based on its orbital path. Galaxies like the Milky Way contain debris from a large number of accreted satellites, and so this feature can be expected to be quite prevalent if the halo is populated with massive substructures, since the large volume occupied by stellar debris would all but ensure some number of reasonably strong interactions. As a result, current and future sky surveys should be able to discover this signature of the presence of substructures in the halo by identifying stars at large radii with strongly correlated velocities. Previous work (Ibata et al. 22; Johnston et al. 22; Mayer et al. 22; Peñarrubia et al. 26) has suggested that heating of the streams would result from interactions with substructure, increasing the velocity disper-