Arxiv:Astro-Ph/0609104V2 31 Jan 2007
Total Page:16
File Type:pdf, Size:1020Kb
Draft version July 23, 2018 A Preprint typeset using LTEX style emulateapj v. 08/22/09 THE VIRGO STELLAR OVER-DENSITY: MAPPING THE INFALL OF THE SAGITTARIUS TIDAL STREAM ONTO THE MILKY WAY DISK David Mart´ınez-Delgado1,2, Jorge Penarrubia˜ 3,4, Mario Juric´5,6, Emilio J. Alfaro2, Zeljko Ivezic´5 Draft version July 23, 2018 ABSTRACT The recently discovered Virgo stellar over-density, which expands over 1000deg2 perpendicularly to the Galactic disk plane (7<Z < 15 kpc, R 7 kpc), is the largest∼ clump of tidal debris ever detected in the outer halo and is likely related with∼ the accretion of a nearby dwarf galaxy by the Milky Way. We carry out N-body simulations of the Sagittarius stream to show that this giant stellar over-density is a confirmation of theoretical model predictions for the leading tail of the Sagittarius stream to cross the Milky Way plane in the Solar neighborhood. Radial velocity measurements are needed to confirm this association and to further constrain the shape of the Milky Way dark matter halo through a new generation of theoretical models. If the identification of Virgo over-density and the Sagittarius leading arm is correct, we predict highly negative radial velocities for the stars of Virgo over-density. The detection of this new portion of the Sagittarius tidal stream would represent an excellent target for the on-going and future kinematic surveys and for dark matter direct detection experiments in the proximity of the Sun. Subject headings: galaxies: individual (Sagittarius) — galaxies:interactions—Galaxy: structure — Galaxy: halo 1. INTRODUCTION dwarf galaxy (Martin et al. 2004; Bellazzini et al. 2004; In the last decade, various large scale CCD surveys Mart´ınez-Delgado et al. 2005), moves on a prograde, (Sloan Digital Sky Survey, 2MASS7, QUEST) have for nearly circular orbit (Pe˜narrubia et al. 2005), according the first time found strong evidence for the presence of to the recent proper motions measured by Dinescu et al. a significant amount of substructure in the halo of the (2005). Additionally, Rocha-Pinto et al.(2004) reported Milky Way in form of long tidal streams or stellar clumps the discovery of a new possible tidal stream in the Trian- (Newberg et al. 2002; Majewski et al. 2003). These have gulum/Andromeda constellations, although subsequent been interpreted as the remnants of latest mergers occur- theoretical simulations identified this tidal debris as a ring as a part of the hierarchical build-up of our Galaxy. part of a more distant, metal poor wrap of the Mono- In addition, N-body simulations of these merger events ceros tidal stream progenitor (Pe˜narrubia et al. 2005; (Helmi 2004; Mart´ınez-Delgado et al. 2004a; Law, John- Mart´ınez-Delgado et al. 2007, in preparation). ston & Majewski 2005; Pe˜narrubia et al. 2006) have In addition to these two large tidal streams, the been extensively used to, firstly, associate the tidal de- analysis of large samples of stellar tracers of halo sub- bris with their progenitor galaxies, secondly, determine structure (RR Lyrae, F-type turnoff stars) have revealed the dynamical history of the progenitor satellites and, the presence of numerous stellar over-densities possibly thirdly, constrain the distribution of dark matter in the associated with unknown merger events in the outer Milky Way. The ultimate goal of these studies is to infer halo. In the Virgo constellation, several groups have arXiv:astro-ph/0609104v2 31 Jan 2007 the formation process of the Milky Way from the prop- reported numerous detections of stellar over-densities erties of the present fossil records. using different tracers: Much of the observational and theoretical effort has (i) The most conspicuous one is a diffuse concentration of metal poor, old stars situated at 20 kpc from the focused specially on the two largest tidal streams dis- ∼ covered so far: the tidal stream of the Sagittarius (Sgr) Sun in the direction of the Virgo constellation, first dwarf galaxy, which wraps around the Galaxy in a highly detected by means of the identification of its associated inclined orbit with respect to the Milky Way disk (Ibata main-sequence turnoff (MSTO) in a color-magnitude et al. 2001; Mart´ınez-Delgado et al. 2001; Majewski diagram (CMD) of a SDSS slice (named S297+63-20.0; et al. 2003) and the Monoceros tidal stream (Yanny et Newberg et al. 2002). Additionally, Vivas & Zinn (2003) reported a clump of 21 RR Lyrae stars located at al. 2003; Ibata et al. 2003), a low-latitude tidal stream ∼ whose progenitor, possibly the controversial Canis Major 20 kpc from the sun in the same region of the sky (see also Ivezic et al. 2005 and Vivas & Zinn 2006). This h 1 Instituto de Astrof´ısica de Canarias, La Laguna, Spain over-density (named the 12 .4 clump) extends over an 2 Instituto de Astrof´ısica de Andaluc´ıa (CSIC), Granada, Spain area of the sky of 15◦ 1◦ (corresponding to to a spatial 3 Max-Planck Institute fuer Astronomie, Heidelberg, Germany × 4 size of 5.3 kpc 0.3 kpc at D 20 kpc) and display a Department of Physics and Astronomy, University of Victoria, line-of-sight depth× of only 1.5kpc≃ (Vivas 2002), which Victoria BC, V8P1A1, Canada ≃ 5 Department of Astrophysical Sciences, Princeton, NJ 08544 is considerably smaller than that of a similar RR Lyrae USA 6 over-density detected at the apocenter of the Sgr stream Institute for Advanced Study, School of Natural Sciences, Ein- ( 4 kpc; Vivas 2002). Owing to its possible tidal origin, stein Drive, Princeton, NJ 08540 ∼ 7 Two Micron All Sky Survey they denoted it as the Virgo stellar stream (VSS). 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 dwarf galaxy. 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 Milky Way 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 metallicity 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.