arXiv:astro-ph/0609104v2 31 Jan 2007 SonDgtlSySre,2MASS Survey, Sky Digital (Sloan ti rv,Pictn J08540 NJ Princeton, Drive, stein l 03 bt ta.20) o-aiuetdlstream tidal low-latitude Major a et Canis (Yanny 2003), controversial the stream al. Majewski possibly tidal progenitor, et 2001; Monoceros whose Ibata the 2003; and al. al. 2003) Mart´ınez-Delgado et (Ibata al. disk 2001; Way et Milky the al. to (Sgr) highly et respect a Sagittarius in with dis- the orbit of streams the inclined around stream tidal wraps tidal which largest galaxy, the dwarf two far: the so covered on specially prop- focused the records. from fossil Way present Milky infer the to the of is of erties studies process the these of formation in and, goal the matter ultimate satellites dark The of Way. progenitor Milky distribution the the of determine constrain de- thirdly, secondly, history have tidal , dynamical the 2006) progenitor the associate their firstly, al. with to, bris Pe˜narrubia et used 2005; extensively been Majewski events John- & Law, merger 2004a; ston Galaxy. these al. our Mart´ınez-Delgado of et of 2004; simulations build-up (Helmi N-body hierarchical the addition, of In occur- part have mergers a These latest as of 2003). ring remnants al. the the et as clumps interpreted Majewski of stellar been 2002; or halo streams al. tidal the of et long (Newberg presence in of form the substructure in Way for of Milky evidence amount strong significant found a time first the USA itraB,VPA,Canada V8P1A1, BC, Victoria H IG TLA VRDNIY APN H NALO H SA THE OF INFALL THE MAPPING OVER-DENSITY: STELLAR THE 7 2 6 5 4 3 rpittpstuigL using 2018 typeset 23, Preprint July version Draft ntels eae aiu ag cl C surveys CCD scale large various decade, last the In 1 has effort theoretical and observational the of Much nttt eAto´sc eAdlcı CI) rnd,S Granada, Andaluc´ıa (CSIC), Astrof´ısica de de Instituto nttt o dacdSuy colo aua cecs E Sciences, Natural of School Study, Advanced for Institute eateto srpyia cecs rneo,N 08544 NJ Princeton, Sciences, Astrophysical of Department eateto hsc n srnm,Uiest fVictori of University Astronomy, and Physics of Department a-lnkIsiuefe srnme edleg German Heidelberg, Astronomie, fuer Institute Max-Planck nttt eAto´sc eCnra,L aua Spain Laguna, La Canarias, Astrof´ısica de de Instituto w irnAlSySurvey Sky All Micron Two ai Mart David oteGlci ikpae(7 plane disk Galactic the to xeln agtfrteo-on n uueknmtcsresan surveys . Sagittarius headings: kinematic the the of future Subject proximity of and the portion on-going in the new experiments identifi for this rad the target negative of If highly excellent detection predict we The models. sh correct, Radial theoretical is over-density. the of arm constrain neighborhood. leading generation further Sagittarius Solar new to the the a and in through n association plane halo this a Way fo confirm of Milky predictions to accretion the model needed the stream cross theoretical Sagittarius with to of the related stream confirmation of likely a simulations is N-body is out over-density and carry halo We Way. outer Milky the in detected h eetydsoee ig tla vrdniy hc expands which over-density, stellar Virgo discovered recently The ´ 1. ınez-Delgado A INTRODUCTION T E tl mltajv 08/22/09 v. emulateapj style X aais niiul(aitru)—glxe:neatosGlx:st galaxies:interactions—Galaxy: — (Sagittarius) individual galaxies: aay halo Galaxy: 1,2 og Pe Jorge , 7 US)hv for have QUEST) , < Z < NOTEMLYWYDISK WAY MILKY THE ONTO narrubia ˜ 5kc R kpc, 15 rf eso uy2,2018 23, July version Draft ABSTRACT pain in- y a, 3,4 ∼ ai Juri Mario , p) stelretcupo ia ersever debris tidal of clump largest the is kpc), 7 fmtlpo,odsasstae at situated concentration stars diffuse old a is poor, one metal conspicuous of have most over-densities groups The stellar (i) several outer of tracers: , different the detections using Virgo in numerous the events reported In merger possibly unknown halo. over-densities sub- with stellar halo numerous associated revealed of of have tracers stars) presence turnoff stellar the F-type of Lyrae, (RR samples structure large of Mono- analysis 2005; the of preparation). al. in wrap (Pe˜narrubia et 2007, poor progenitor al. a metal Mart´ınez-Delgado stream et as distant, tidal debris more ceros a tidal of this part identified Trian- simulations the subsequent in theoretical although stream tidal , possible new gulum/Andromeda reported a al. al.(2004) of et discovery et Dinescu the Rocha-Pinto by Additionally, measured motions prograde, (2005). according proper a recent 2005), the on al. 2004; to (Pe˜narrubia moves et al. orbit et circular 2005), Bellazzini nearly 2004; al. al. Mart´ınez-Delgado et et (Martin galaxy dwarf loIei ta.20 n ia in20) This 2006). Zinn (see & sky Vivas the and of at region 2005 the located same (named stars al. the over-density Lyrae in et sun Ivezic RR the also 21 from of kpc (2003) Zinn 20 clump & S297+63-20.0; Vivas a Additionally, (named reported 2002). color-magnitude slice al. SDSS first a et a Newberg constellation, in of associated (CMD) its Virgo (MSTO) of diagram the identification turnoff the of of main-sequence means direction by the detected in Sun rao h k f15 of sky the of area ieo 5 of size ieo-ih et fonly of depth line-of-sight vrdniydtce tteaoetro h g stream Lyrae Sgr RR the of similar apocenter a ( the of at that detected than over-density smaller considerably is hydntdi steVroselrsra (VSS). stream stellar Virgo the as it denoted they ∼ nadto oteetolretdlsras the streams, tidal large two these to addition In p;Vvs20) wn oispsil ia origin, tidal possible its to Owing 2002). Vivas kpc; 4 . c ´ kpc 3 5,6 mloJ Alfaro J. Emilio , h edn alo h Sagittarius the of tail leading the r p fteMlyWydr matter dark Way Milky the of ape o akmte ietdetection direct matter dark for d a eoiisfrtesaso Virgo of stars the for velocities ial × aino ig vrdniyand over-density Virgo of cation ia temwudrpeetan represent would stream tidal over 0 . p at kpc 3 oso htti in stellar giant this that show to ◦ ∼ ab wr aayb the by galaxy dwarf earby eoiymaueet are measurements velocity × 1000deg 1 ITRU IA STREAM TIDAL GITTARIUS ◦ 12 ≃ crepnigt oaspatial a to to (corresponding D h 1 4clump .4 . 2 kc(ia 02,which 2002), (Vivas 5kpc ≃ 2 ejoIvezi Zeljko , perpendicularly 0kc n ipa a display and kpc) 20 ∼ utr — ructure xed vran over extends ) 0kcfo the from kpc 20 c ´ 5 ∼ 2 Mart´ınez-Delgado et al.

Mart´ınez-Delgado et al. (2004a) suggested that the where several theoretical models predict the presence of VSS might be related to the complex debris structure Sgr galaxy debris. In order to explore a possible associa- of the Sgr stream predicted by the theoretical models in tion of the VSS and VOD with the Sgr tidal stream, we this region of the sky. However, spectroscopic follow-up shall compare their derived spatial and kinematic prop- of a sample of these RR Lyrae stars by Duffau et al. erties with those obtained from up-to-date theoretical (2006) showed a fairly narrow velocity distribution that models. pointed to a progenitor galaxy less massive than the Our Galaxy model consists of a Miyamoto-Nagai Sgr dwarf and, therefore, to a hitherto unknown tidally (1975) disk, a Hernquist (1990) bulge and a Navarro, disrupted . Using SDSS photometric data, Frenk & White (1996) dark matter halo (hereafter, NFW these authors estimated that the VSS covers roughly an halo). The gravitational potential of each of those com- area of 106 deg2 in the sky. ponents in cylindrical coordinates is: GM (ii) More recently, Juri´c et al. (2006) detected Φ (r)= d , (1) a remarkable excess of main-sequence stars towards d − R2 + (a + √z2 + b2)2 the direction of Virgo constellation, named the Virgo q over-density (VOD). These authors construct a three- GMb Φb(r)= , (2) dimensional density distribution of stars −r + c based on the position of 48 million stars obtained by GM q the SDSS in order to identify∼ stellar over-densities. The r h h Φh( )= rvir rvir (3) ln(1 + ) 2rs × VOD density map reveals a giant, diffuse lump of stars rs − rs+rvir without any evidence of a distinct core that extends over ∞ m(u) du more than 1000 deg2 (see Fig. 1) towards the direc- 2 . ◦ ◦  1+ m(u) 2 −  tion (l,b) =∼ (300 , 65 ). This stellar over-density shows Z0 (1 + u) qh + u an remarkably extended vertical structure perpendicular p There, rvir, rs are the virial and scale radii, respectively, to the Galactic plane (∆Z 8 kpc, with heliocentric dis- 2 2 2 tances between 5-17 kpc) that≃ possibly crosses the Galac- Mh = Mh(rvir) and r = R + z . The potential of tic plane, extending into the southern Galactic hemi- an axi-symmetric NFW halo was calculated from Chan- sphere (see also Fuchs, Phleps & Meisenheimer 2006). drasekhar (1960) using elliptic coordinates The density of VOD main-sequence stars peaks at 16 ∼ R2 z2 kpc from the Sun (or RR Lyrae magnitude of 16.7), in m2(u)= + , (4) ∼ 2 2 2 good agreement with the RR Lyrae distance determina- rs (1 + u) rs (qh + u) tion. These authors suggest that the giant stellar over- density might be related with a hitherto unknown merger where qh is the axis-ratio of iso-density surfaces. In this event of a nearby, low dwarf galaxy .Alterna- contribution, we shall denote oblate and prolate halos as tively, Xu, Deng & Hu (2006) suggest that this stellar those with qh < 1 and qh > 1, respectively. For spherical lump could be related to a non-axisymmetric Galactic halos Eq. (3) reduces to stellar spheroid (see Sec. 3.1). Although VOD and VSS are coincident in position and GMh ln(1 + r/rs) Φh = rvir rvir (5) display a huge sky-projected size, the different vertical −ln(1 + rs ) rs+rvir r ≡ structure suggests that they are independent systems. In − 2 ln(1 + r/rs) absence of radial velocity measurements for the VOD, it Vc . is not possible to address whether they share a common − r/rs origin or are completely unrelated. For this reason, we will refer to them with different names in this paper, Following Johnston et al. (1999) we fix the disk and 11 following the nomenclature given to each over-density in bulge parameters as Md = 1.0 10 M⊙, Mb = 3.4 10 × × the original papers. 10 M⊙, a =6.5kpc, b =0.26 kpc and c = 0.7kpc. The The main aim of this paper is to investigate the nature Milky Way halo parameters at z = 0 were taken from of the VOD and to analyze whether this stellar lump Klypin, Zhao & Somerville (2002) being Mh = 1.0 12 × could be related to the Sgr stream through a comparison 10 M⊙, rvir = 258kpc, rs = 21.5kpc, which leads to a between up-to-date theoretical models and the available concentration at the present epoch of c = rvir/rs = 12. observational data. The structure, kinematics and pos- According to the results of Pe˜narrubia et al. (2006), the sible origin of the VSS will be analyzed in a companion evolution of the host Galaxy potential is not reflected in paper. the present properties of tidal streams, so that we use a In Sec. 2 the theoretical models of Sgr tidal stream static Milky Way potential for simplicity. are described. In Sec.3 we discuss the possible origins The actual disruption of satellites is modelled by ‘live’ of the VOD. In particular, we explore the presence of (i.e. self-consistent, self-gravitating) N-body realiza- tidal debris from the Sgr tidal stream in Virgo by means tions of a King (1960) model, with a dimensionless of update theoretical models and compare its expected central potential W0 = 4, or concentration parameter structural parameter with those derived for Virgo. Fi- c log10(rt/rk) 0.84, where rk and rt are the King and nally, conclusions are given in Sec. 4. tidal≡ radii, respectively.≃ Our satellite models have N = 105 particles. The initial and final (i.e present) masses 2. THEORETICAL MODELS 9 8 are Ms(t ) = 10 M⊙ and Ms(tf )=5 10 M⊙, respec- 0 × The giant stellar over-density detected in the Virgo tively. The initial King and tidal radius are rk(t0)=0.58 constellation is located in a region of the Milky Way kpc and rt(t0)=4.01 kpc. The origin of remnants in Virgo 3

The equation of motion for each satellite particle is: With those inputs, we have evolved the N-body satel- 2 lite model from 4 Gyr in the past to the present. d ri = (Φs +Φd +Φb +Φh)i (6) dt2 −∇ 3. THE ORIGIN OF THE VIRGO OVER-DENSITY where Φs is the self-gravitational potential of the satellite In this Section we analyse the likelihood of various sce- galaxy and i =1, ..., N. We use superbox (Fellhauer et narios that may shed light on the nature of the VOD. al. 2000) to calculate Φs at each time-step and solve Eq. 6 through a leap-frog scheme with a constant time-step of 3.1. The signature of a tri-axial Galactic stellar halo ∆t = 0.65 Myr, which is about 1/100th the satellite’s An alternative explanation for the presence of the un- dynamical time. expected MS feature in the CMD reported in several Note that the Galactic potential is static and, therefore, studies (and so-called the VOD: Newberg et al. 2002; Eq. 6 does not implement the response of the Milky Way Duffau et al. 2006; Juric et al. 2006; see Sec.1) is to to the presence of Sgr. We also neglect the effects of postulate the existence of a non-axisymmetry Galactic dynamical friction on the Sgr’s orbit (as suggested by model component, such as a triaxial halo (Newberg & Law et al. 2005). Yanny 2005; Xu et al. 2006). We have explored this We have reproduced the N-body models presented in possibility by inspecting a sample of deep CMDs (with Law et al. (2005). These models reproduce the following 8 limiting magnitude Rlim 24.5) of Galactic fields taken observational constraints : with the Wide Field Camera∼ (WFC) at the prime focus of 1. The heliocentric position and radial velocity of the Isaac Newton 2.5-meter telescope (La Palma,Spain) dur- Sgr dwarf are (D,l,b) = (24kpc, 5.6◦, 14.2◦) and ing different observing runs devoted to the search of tidal − vr = 171 km/s, respectively. debris associated to known Galactic satellites or globular clusters (Mart´ınez-Delgado et al. 2004b). Table 1 lists 2. The orbital plane of Sgr has an inclination with the position of our fields with positive and negative de- respect to the Milky Way disk of i 76◦. tections of the MS feature in the CMD. These negative ≃ detections can be interpreted as the lack of tidal debris 3. The averaged heliocentric distance for Sgr leading or that its is extremely low to be de- debris is D 50 kpc. tected in the small field of view of the WFC (35′ 35′)9. ∼ × Items 1. and 2. above show that we have observational Fig. 1 shows the position of these negative detection measurements for five of the six coordinates that deter- (open circles) over-plotted in the density map of VOD derived by Juri´cet al (2006) from SDSS observations. mine the Sgr orbit. Law et al. (2005) surveyed the un- ◦ known coordinate (namely, the tangential velocity com- The most interesting result is a negative detection 20 ∼ ponent, vtan) so that the resulting N-body models repro- North of the center of our target field (see Table 1) that duced all the observational constraints listed. clearly shows a decrease of stream star counts in a small Among other free parameters, Law et al. (2005) also area of the sky. Although our sparse pencil beam survey explored which halo axis-ratio (qh) would produce the cannot constrain the spatial extent of this stellar popu- best-fitting model to the available observational data. lation in the sky, the distribution of negative detections They found that, whereas the geometry and kinematics strongly indicates that this MS population is not likely to of the trailing arm are scarcely sensitive to the adopted be a manifestation of a non-axisymmetric Galactic model qh, radial velocity measurements of the leading arm were component, but is probably associated with the remnants best matched by prolate (qh > 1) halo models, in agree- of a dwarf galaxy under tidal disruption, in agreement ment with Helmi (2004). However, Johnston, Law and with previous studies. Juri´cet al. (2006) have reached Majewski (2005) showed that Sgr models in prolate ha- the same conclusion by analyzing the vertical (Z) distri- los cannot reproduce the observed precession rate in bution of VOD debris from their sample of SDSS G stars the youngest pieces of the Sgr stream. Furthermore, (0.10 < r i< 0.15). − they were able to constrain the halo axis-ratio to be qh = 0.83 0.92, excluding models with qh > 1 at 3.2. A tidal debris of the Sgr dwarf? − − a 3σ level. In view of these results, we have performed 3.2.1. Sky projected density maps N-body simulations where the halo axis-ratio was fixed It is well understood both from observations and the- either to qh =0.8 orto qh =1.4, excluding for the sake of brevity spherical halos, which neither reproduce the ra- oretical modeling that the Sgr leading arm turns around dial velocity variation along the stream nor its precession and falls towards the Galactic plane in the proximity of rate. the Solar neighborhood, yielding a vertical debris struc- In a coordinate system where the Sun is located at ture that might cross the Galactic plane and possibly (X,Y,Z) = (8, 0, 0) kpc, with a velocity of (U,V,W ) = extend into the southern Galactic hemisphere (Majew- ( 10, 220, 7) km/s (Binney & Merrifield 1998), the Sgr ski et al. 2003; Mart´ınez-Delgado et al. 2004a; Law et velocities− − that best fit the kinematic and spatial distri- al. 2005). This downward flow of the Sgr stars onto the butions of stream debris are (237, 35, 220) km/s for Galactic disk may share similar characteristic to those − derived for the VOD in the stripe covered by the SDSS. qh = 0.8 (oblate halo) and (244, 39, 249) km/s for q =1.4 (prolate halo). − As a guide for the reader, we illustrate this putative sce- h nario in Fig. 2 together with the identification of the 8 It is also important to remark that we have not input any different parts of the Sgr tidal stream obtained from our observational data of the VOD in the fitting of these simulations, that are only based in the fitting of an independent set of data of 9 Mart´ınez-Delgado et al. 2004a estimated the limit SB detection the Sgr tidal stream. of our CMD technique is Σ ∼ 32 mag arcsec−2 4 Mart´ınez-Delgado et al. theoretical models (see Sec. 5 for a three-dimensional rection in a localized area of the sky.10 comparison). The prolate (qh=1.4) halo model (bottom panel) pro- To examine whether the VOD is a part of the Sgr tidal duces a slightly worse estimate of the location of VOD structure, we compare the predictions from the theoreti- center and its general shape. It does predicts a substan- cal models outlined in Sec. 2 to the SDSS data presented tial amount of faint, diffuse debris in the 180

At present there is no kinematic information available In this contribution, we have investigated the origin on the stars of Virgo over-density. Radial velocities are of the Virgo over-density discovered from the analysis of necessary in order to establish unambiguously whether the SDSS data. The VOD is the largest clump of tidal the Virgo over-density is associated to the Sgr stream. debris ever detected in the outer halo of the Milky Way. Theoretical models predict that, independently of the We find that the position, apparent angular size, stellar halo shape, the leading arm of the Sgr stream is present density and three dimensional geometry match the the- in this area of the sky and moves on average with nega- oretical predictions from the present Sgr stream models. tive radial velocities as it falls onto the Milky Way disk. According to these models, the Virgo over-density cor- However, the exact magnitude of the mean radial velocity responds to Sgr leading arm falling onto the Milky Way is a vital discriminator for the halo shape, as prolate and disk. In agreement with other portions of the Sgr stream, oblate halos lead to fairly distinct trends of Vr along the the spatial distribution of the VOD is best reproduced by leading arm (Mart´ınez-Delgado et al. 2004a; Helmi 2004, models with oblate halos qh < 1. However, we must re- Law et al. 2005). We must remark that at present there mark that a unambiguous association between the VOD is no information on the kinematics of the Sgr leading and the Sgr stream cannot be established in absence of arm in this area of the sky, since this part of the stream . was not covered by the spectroscopic survey of M-giants If the Virgo over-density proves to belong to the Sgr carried out by Majewski et al. (2004). stream it might represent the key piece to solve the para- There have been radial velocity measurements of RR dox described above, as we have for the first time the Lyrae stars in the Virgo constellation that have shown opportunity to measure the radial velocity of millions of the presence of a spatial over-density (called the 12h.4 stream stars just in the region where observations and clump, see Sec. 1) located at D 20 kpc. Duffau et theoretical models do not match. al. (2006) find a mean radial velocity≃ in the galactic rest Additionally, the detection of a new portion of the frame of (Vr)gsr 99.8 km/s (see Table 1) with a fairly stream’s leading arm provides strong constraints on fu- low velocity dispersion≃ ( σ=17.3 km s−1), unexpected for ture theoretical models and will help to address whether the smooth population of RR Lyrae in the Milky Way, the Sgr stream passes through the solar neighbourhood that suggests a possible tidal stream detection. Accord- as it crosses the Milky Way plane (as predicted by Kundu ing to the theoretical models of the Sgr stream, the posi- et al. 2002, Majewski et al. 2003; Mart´ınez-Delgado et tive value of (Vr)gsr excludes a membership of these stars al. 2004). In that case, the Sgr stream would repre- to the leading arm. However, it is unclear whether the sent an important target for the on-going RAVE sur- 12h.4 clump is not associated to the Sgr stream at all, vey (Steinmetz et al. 2006), whose goal is to measure 6 since the oblate halo models (qh =0.8 in this work) also the radial velocity of around 10 stars in both hemi- predict the presence of the trailing arm stars in that vol- spheres. Thanks to the large sample of radial velocities, if ume of the Milky Way (see Fig. 5, where the blue box the Milky Way crossing occurs at a heliocentric distance represents the 12h.4 clump position) with an averaged D < 3 kpc, its presence will be detected in form of an radial velocity of (Vr)gsr= 90 km/s, fairly similar to that excess of stars with vr < 0 in the Northern hemisphere reported by Duffau et al. (2006) ((V ) =+83 or +100 (b 0) and with vr > 0 for b < 0 with respect to the r gsr ≥ km/s, depending of their sample selection). Therefore, if background population of Milky Way stars (Pe˜narrubia, the 12h.4 clump proves to belong to the Sgr trailing arm Navarro & Helmi, in preparation). The presence of a tidal stream in the solar neighbourhood also represents it would favor the models with qh < 1. an excellent target for Weakly Interacting Massive Parti- cles (WIMPs) direct detection experiments, since a sub- 4. DISCUSSION structure of dark matter in the proximity of the Sun Theoretical models of the Sgr stream are not able to would yield a “cold” flow of WIMPs thought the detector reproduce all the observational constraints. All models of either of these dark matter experiments (see Savage, agree that the geometry of the tidal stream (i.e its 3D Freese and Gondolo 2006 and Gondolo et al. 2005 for spatial distribution) clearly rules out halo models with a details). This would make the solar vicinity an excellent prolate shape (Mart´ınez-Delgado et al. 2004a, Johnston lab for probing the physical properties and nature of the et al. 2005, Law et al. 2005) and that Milky Way halos dark matter. with qh < 1 are preferred at a 3σ level. Paradoxically, the models that best explain the stream structure fail to reproduce the radial velocity trend along the leading arm unless prolate halo models (qh > 1) are invoked (Helmi We thanks to Brian Yanny and Kathy Vivas for making 2004). Apparently, the present status of the Sgr stream available to us some resultson the Virgo tidal stream dur- modeling is that of a fish biting its tail. ing the preparation of this paper. We also thanks to the It is not clear how to reconcile observations and the- anonymous referee for useful comments that helped to ory. Possible explanations for this mismatch are: (i) the improve this manuscript. DMD acknowledges funding by Milky Way models are over-simplistic (although similar the Spanish Ministry of Education and Science (Ramon models successfully reproduce other tidal streams such y Cajal contract and research project AYA 2001-3939- as Monoceros, Pe˜narrubia et al 2005, and the Magellanic C03-01). JP thanks Julio Navarro for financial support. stream, Ruzicka, Palous & Theis 2006, both streams be- Funding for the SDSS and SDSS-II has been pro- ing best reproduced by models with oblate halos) or (ii) vided by the Alfred P. Sloan Foundation, the Partic- the presence of possible systematic errors on the data of ipating Institutions, the National Science Foundation, the Sgr tidal stream used in the computation of these the U.S. Department of Energy, the National Aeronau- simulations (as suggested by Chou et al. 2006). tics and Space Administration, the Japanese Monbuka- The origin of remnants in Virgo 7 gakusho, the Max Planck Society, and the Higher Educa- the Joint Institute for Nuclear Astrophysics, the Kavli tion Funding Council for England. The SDSS Web Site Institute for Particle Astrophysics and Cosmology, the is http://www.sdss.org/. Korean Scientist Group, the Chinese Academy of Sci- The SDSS is managed by the Astrophysical Research ences (LAMOST), Los Alamos National Laboratory, the Consortium for the Participating Institutions. The Par- Max-Planck-Institute for Astronomy (MPIA), the Max- ticipating Institutions are the American Museum of Nat- Planck-Institute for Astrophysics (MPA), New Mexico ural History, Astrophysical Institute Potsdam, Univer- State University, Ohio State University, University of sity of Basel, Cambridge University, Case Western Re- Pittsburgh, University of Portsmouth, Princeton Uni- serve University, University of Chicago, Drexel Univer- versity, the United States Naval Observatory, and the sity, Fermilab, the Institute for Advanced Study, the University of Washington. Japan Participation Group, Johns Hopkins University, REFERENCES Bellazzini, M., Ibata, R., Monaco, L., Martin, N., Irwin, M. J., Majewski, S. R., Kunkel, W. E., Law, D. R., Patterson, R. J., Lewis, G. F. 2004, MNRAS, 354, 1263 Polak, A. A.; Rocha-Pinto, H. J.; Crane, J. D., Frinchaboy, P. Binney J., Merrifield M., 1998, Galactic Astronomy. Princeton M.; Hummels, C. B., Johnston, K. V.; Rhee, J.; Skrutskie, M. University Press, Princeton, New Jersey F.; Weinberg, M. 2004, AJ, 128, 245 Chandrasekhar, S. 1960, Principles of Stellar Dynamics, Dover, Martin, N., Ibata, R. A., Bellazzini, M., Irwin, M. J., Lewis, G. F., New York Dehnen, W. 2004,MNRAS, 348, 12 Chou. M.-Y., Majewski, S. R., Cunha, K, Smith, V. V., Patterson, Mart´ınez-Delgado, D.,Aparicio, A., G´omez-Flechoso, M. A., R. J., Mart´ınez-Delgado,D., Law, D. R., Crane, D. J., Mu˜noz, Carrera, R. 2001, ApJ, 549, L199 R. R., Garc´ıa L´opez, R., Geisler, D., Skrutskie, M. F. 2006, ApJ, submitted (astro-ph/0605101) Mart´ınez-Delgado, D., G´omez-Flechoso, M. A., Aparicio, A., Dinescu, D. I., Mart´ınez-Delgado, D., Girard, T. M., Pe˜narrubia, Carrera, R. 2004a, ApJ, 601, 242 J., Rix, H.-W., Butler, D., & van Altena, W. F. 2005, ApJ, 631, Mart´ınez-Delgado, D., Dinescu, D. I, Zinn, R., Tutsoff, A.; Cot´e, P., L49 Boyarchuck, A. 2004b, in Satellites and Tidal Streams, Prada, Duffau, S., Zinn, R, Vivas, A. K., Carraro, G., M´endez, R. A., Mart´ınez-Delgado and Mahoney eds. PASP vol. 327,255 Winnick, R., Gallart, C. 2006, ApJ, 636, 97 Mart´ınez-Delgado, D., Butler, D. J., Rix, H.-W., Franco, I., Fellhauer M., Kroupa P., Baumgardt H., Bien R., Boily C. M., Pe˜narrubia, J., Alfaro, E. J., Dinescu, D. I. 2005, ApJ, 633,205 Spurzem R., Wassmer N., 2000, NewA, 5, 305 Miyamoto M., Nagai R., 1975, PASJ, 27, 533 Fuchs, B., Phleps, S., Meisenheimer, K. 2006, A&A, 457, 541 Nakatura et al. 2000, ApJ, 540, 1172 Gondolo, P., Freese, K., Newberg, H. J., Lewis, M. 2005, New Navarro J., Frenk C.S., White S.D.M., 1996, ApJ, 462, 563 Astronomy Review, 49, 193 Newberg, H.J. et al. 2002, ApJ, 569,245 Helmi, A. 2004, ApJ, 610, L97 Newberg, H.J., Yanny, B. 2005, in Physics at the end of the Hernsquit, L. 1990, ApJ, 356, 359 Galactic Cosmic Ray Spectrum Juri´cet al. 2006, ApJ, submitted (astro-ph/0510520) , Journal of Physics Conference Ibata, R. A , Irwin, M.J., Lewis, G. F., Stolte, A. 2001, ApJ, 547, Series, Thomsonan & Sokolsky eds. in press (astro-ph/0507671) L133 Pe˜narrubia, J. et al. 2005, ApJ, 626, 128 Ibata, R. A., Irwin, M. J., Lewis, G. F., Ferguson, A., Tanvir, N. Pe˜narrubia, J., Benson, A. J., Mart´ınez-Delgado, D., & Rix, H. W. 2003, MNRAS, 340, L21 2006, ApJ, 645, 240 Iv´ezic, Z., Lupton, R. H., Schlegel, D., Smolcic, V., Johnston, D., Rocha-Pinto, H. J.; Majewski, S. R.; Skrutskie, M. F.; Crane, J. Gunn, J. E., Knapp, G. R., Strauss, M. A., Rockosi, C. M. 2004 D.; Patterson, R. J. , 2004, ApJ, 615, 732 in Satellites and Tidal Streams, Prada, Mart´ınez-Delgado and Ruzicka A., Palous J., Theis Ch., 2006, A&A, in press Mahoney eds. PASP vol. 327,104 (astro-ph/0608175) Johnston K.V., Majewski S.R., Siegel M.H., Reid I.N., Kunkel Savage, C., Freese, K., Gondolo, P. 2006, astro-ph/0607121 W.E., 1999, AJ, 118, 1719 Steinmezt et al., 2006,AJ, in press (astro-ph/0606211) Johnston K.V., Law D.R., Majewski S.R., 2005, ApJ, 619, 800 King I.R., 1966, AJ, 71, 65 Vivas, K. V. 2002, Ph. D. thesis Klypin A., Zhao H., Somerville R., 2002, ApJ, 573, 597 Vivas, K. V., Zinn, R. 2003, Memorie della Societa Astronomica Kundu, A. et al. 2002, ApJ, 576, L125 Italiana, vol. 74, p.928 Law, D. R., Johnston, K. V., Majewski, S. R. 2005, ApJ, 619, 807 Vivas, K. V., Zinn, R. 2006, AJ, 132, 714 reference Xu, Y., Deng, Majewski, S. R., Skrutskie, M., Weinberg, M., Ostheimer, J. 2003, L. C., Hu, J. Y. 2006, MNRAS, submitted (astro-ph/0602565) ApJ, 599, 1082 Yanny, B. et al. 2003, ApJ, 588, 824 8 Mart´ınez-Delgado et al.

Fig. 1.— Density maps from SDSS of the North Galactic hemisphere. Position of our positive (solid triangles) and negative (open circles) detection of the Virgo stellar over-density given in Table 1, overplotted on a Lambert projection of the density map of Virgo given in Juri´cet al. (2005). This map was generated by selecting MS stars with b < 0 deg, 0.2 < g − r < 0.3 and 20 30 is the Milky Way disk protruding into the survey region.

TABLE 1 Positive and negative detections of Virgo from colour-magnitude diagrams

RA(J2000) Dec(J2000) l b Detection 12:24:00.00 -1:00:00.0 288.62 61.11 Positive 12:24:00.00 -0:30:00.0 288.39 61.59 Positive 10:20:00.00 40:00:00.0 180.71 56.24 Negative 13:48:16.00 52:08:32.0 103.69 62.81 Negative 15:37:49.00 69:17:31.0 104.68 41.44 Negative 12:40:00.00 18:20:00.0 285.67 80.82 Negative 13:56:21.00 -27:10:04.0 320.28 33.51 Negative 13:07:53.00 29:21:00.0 65.42 85.75 Negative 10:49:22.00 51:03:04.0 158.57 56.78 Negative 14:07:37.00 11:50:28.0 356.10 66.48 Negative 16:11:05.00 14:57:29.0 28.76 42.18 Negative 14:44:27.00 -1:00:00.0 351.54 50.88 Negative The origin of remnants in Virgo 9

Fig. 2.— Galactic X-Z projection of our N-body simulation of the Sgr tidal stream in a Milky Way dark matter halo with qh = 0.8. For simplicity, only particles unbound less than 2.5 Gyr ago have been plotted. It displays the Sgr leading arm falling onto the Galactic disk, crossing the proximity of the Sun (yellow circle). The particles with a position consistent with that of the VOD are marked with red colour. 10 Mart´ınez-Delgado et al.

Fig. 3.— Density maps of the Sgr leading tail stars with heliocentric distance d<15kpc from our N-body simulation of the Sgr tidal stream in a Milky Way dark matter halo with qh = 0.8 (top panel) andqh =1.4 (bottom panel). The origin of remnants in Virgo 11

Fig. 4.— Top row: comparison of Galactic X-Y projections of the Virgo over-density obtained from SDSS data (left column), with two Sgr stream models, qh = 0.8 (center) and qh = 1.4 (right column). In the left column, we plot the residual stellar density after subtraction of the Galactic stellar density model in a thin ∆Z = 0.5kpc slice at Z = 10 kpc above the Galactic − plane. The color coding represents the stellar overdensity of 0.1

Fig. 5.— Two different 3-dimensional projections (top and bottom rows)of the Sagittarius tidal stream models computed for an oblate halo qh = 0.8 (left column) and a prolate qh = 1.4 halo (right column). The column density of simulation particles is represented by rainbow color-map, with the red being the areas of highest, and blue the lowest density. Over-plotted on the column density map as black dots are the simulation particles themselves. The core of Sgr dwarf is clearly visible as a dense, dark region, due to the high density of simulation particles which are still gravitationally bound to the dwarf’s core. The origin of the coordinate system axes is the Galactic center, with the red axis (X) pointing towards the Sun (yellow circle), the blue axis (Z) pointing towards the north Galactic pole. The region in which Juri´cet al. (2006) detect the Virgo overdensity is marked by the white wire-frame box. The blue wire-frame marks the location of Virgo Stellar Stream as reported by Duffau et al. (2006). .