Arxiv:1805.08326V2 [Astro-Ph.GA] 7 Dec 2018 Dwarf Galaxies (Belokurov Et Al
Total Page:16
File Type:pdf, Size:1020Kb
Draft version December 10, 2018 Typeset using LATEX twocolumn style in AASTeX62 Rotating halo traced by the K giant stars from LAMOST and Gaia Hao Tian1 (LAMOST Fellow) Chao Liu,1 Yan Xu,1 and Xiangxiang Xue1 1Key Laboratory of Optical Astronomy, National Astronomical Observatories, Chinese Academy of Sciences, Datun Road 20A, Beijing 100012, PR China; (Received December 10, 2018; Revised December 10, 2018; Accepted |{) Submitted to |- ABSTRACT With the help of Gaia DR2, we are able to obtain the full 6-D phase space information for stars from LAMOST DR5. With high precision of position, velocity, and metallicity, the rotation of the local stellar halo is presented using the K giant stars with [Fe/H]< −1 dex within 4 kpc from the Sun. By fitting the rotational velocity distribution with three Gaussian components, stellar halo, disk, and a +4 −1 counter-rotating hot population, we find that the local halo progradely rotates with VT = +27−5 km s −1 providing the local standard of rest velocity of VLSR = 232 km s . Meanwhile, we obtain the dispersion +4 −1 of rotational velocity is σT = 72−4 km s . Although the rotational velocity strongly depends on the choice of VLSR, the trend of prograde rotation is substantial even when VLSR is set at as low as 220 km s−1. Moreover, we derive the rotation for subsamples with different metallicities and find that the rotational velocity is essentially not correlated with [Fe/H]. This may hint a secular evolution origin of the prograde rotation. It shows that the metallicity of the progradely rotating halo is peaked within -1.9<[Fe/H]<-1.6 without considering the selection effect. We also find a small fraction of counter- rotating stars with larger dispersion and lower metallicity. Finally, the disk component rotates with +6 −1 +3 −1 VT = +182−6 km s and σT = 45−3 km s , which is quite consistent with the metal-weak thick disk population. Keywords: galaxies: individual (Milky Way) | Galaxy: kinematics and dynamics | Galaxy: halo | Stars: individual (K-giants) 1. INTRODUCTION stream, discovered by Ibata et al.(1994), many streams It has been widely accepted that the Milky Way, at have been discovered with matched filter method (Rock- least the halo, was formed through hierarchically merg- osi et al. 2002; Grillmair & Carlin 2016; Grillmair & ing smaller stellar systems, such as globular clusters and Dionatos 2006a,b,c; Grillmair 2009) and other methods (Mateu et al. 2018). But the streams are not the final arXiv:1805.08326v2 [astro-ph.GA] 7 Dec 2018 dwarf galaxies (Belokurov et al. 2006). With a long time evolving, the stellar members of those merged sys- fate of those merged clusters or dwarf galaxies, a longer tems will be continuously stripped and then form stellar evolving time will make the streams more and more re- streams (Johnston et al. 1999). Such phenomena have laxed and finally become part of the smoothing halo, been widely searched with photometric sky surveys, e.g. i.e. it is not visible in the geometric space (Helmi et SDSS (York et al. 2000) and Pan-STARRS (Bernard al. 1999), even blurred in the phase space (Helmi & de et al. 2016). Since the first stellar stream, Sagittarius Zeeuw 2000; Sanderson et al. 2015). Assuming those merged stellar systems fell in with random orbits (Sanderson et al. 2015), the contribu- Corresponding author: Chao Liu tion of z−axis angular momentum Lz overall should [email protected] be around 0. The net angular mometum is directly related to the merging history, not only including the 2 Tian et al. minor mergers of the globular cluster and dwarf galax- mates led to a debate (Sch¨onrich et al. 2011; Beers et al. ies, such as the Sagittarius Stream, but also includ- 2012). This will be highly improved by the Gaia mission ing the major merger events, such as the recently dis- (Perryman et al. 2001). covered Gaia-Enceladus(Helmi et al. 2018), also named GDR2 (Gaia Collaboration et al. 2018a,b,c,d,e,f) con- Gaia-Sausage(Belokurov et al. 2018), which contibutes tains ∼1.3 billion stars high accuracy proper motions a larger fraction of retrogradely rotating stars. Another and parallaxes, which have brought a series of amazing possibility is that the merged halo was driven to slowly results about the halo, such as the velocity determina- rotate by the non-axisymmetric structure in the disk tion of globular clusters and rotation curve of the Large (e.g. bar) in a long time scale (Athanassoula et al. 2013). Magellanic Cloud (Gaia Collaboration et al. 2018e). Be- The studies on dynamics of the halo have been done sides, Gaia Collaboration et al.(2018e) also derived a since 1980s (Frenk & White 1980) who used 66 clus- lower limit total mass for the Milky Way. ters to constrained the rotation velocity 60±26 km s−1. Guoshoujing Telescope (the Large Sky Area Multi- As the velocity ellipsoid is directly related to the mass Object Fiber Spectroscopic Telescope, hereafter LAM- distribution of the halo (Binney & Tremaine 2008), OST) is a 4-meter, quasi-meridian, reflecting Schmidt more work focused on the distribution of the dispersions telescope which makes it an efficient spectroscopic tele- (Sommer-Larsen et al. 1997; Sirko et al. 2004; Smith et scope (Cui et al. 2012; Deng et al. 2012; Zhao et al. al. 2009; Hattori et al. 2013; Kafle et al. 2013; Loebman 2012). Liu et al.(2017) and Xu et al.(2018) used metal et al. 2018). poor K giants from LAMOST to trace the density pro- Because of the difficulty of observation of proper mo- file of the halo and found an inner oblate, outer nearly tions, especially for distant stars. Morrison et al.(1990) spherical profile. Recently, LAMOST has released its studied the kinematics of the halo and disk using the 5th data (DR5). As low resolution of R = 1800, the G and K giant stars within 4 kpc from the Sun, and uncertainties of Teff and [Fe/H] are typical ∼ 100 K and −1 claimed a prograde rotation of VT = 25 ± 15 km s . ∼ 0:1 dex, respectively, for the stellar spectra with signal They also showed that the metallicity of the thick disk to noise ratio larger than 30. stars can be as low as -1.6. Smith et al.(2009) studied The paper studies the rotation of the local stellar halo the kinematics of the local halo using ∼1700 solar neigh- using the conbination of the GDR2 and LAMOST DR5 bourhood subdwarfs and found that the stellar halo ex- K-giant stars and is constructed as follow. In Section2, hibits no rotation and the rotational velocity Vφ disper- the sample stars are selected. The model of the distibu- sion 82 ± 2 km s−1 , what should be noticed is that the tion of rotational velociy is developed in Section3. The typical velocity errors are large, about 30 to 50 km s−1. results are discussed in Section4. The conclusions are More recently, the progress of the observations and listed in Section5. data analysis techniques, high precision proper motions have been obtained through combining catalogues ob- 2. DATA SAMPLE tained at different epochs cross a long time baseline In this work, we use K giant stars selected from LAM- (Tian et al. 2017; Belokurov et al. 2018). Deason et al. OST DR5 based on the criteria suggested by Liu et al. (2017) used the 3D velocity calculated from the proper (2014) to study the rotation of the nearby halo. After motion of different tracers, including RR Lyrae, Blue cross-matching with GDR2, we obtained the parallax, Horizontal Branch (BHB) and K giant stars, in a distant proper motions, radial velocities and metallicities for the volume from Gaia Data Release (hereafter GDR) 1 and tracer stars. A robust metallicity cut with [Fe/H]< −1 SDSS dataset to study the halo rotation. A prograde is used to remove the majority of disk stars, which leaves rotation was found with V = 14±2±10 km s−1. What T 6660 K giant stars. In this case most of the thin disk is more, the metal richer stars are showing a stronger stars and many thick disk stars will be removed. As we prograde rotation which is explained as a small fraction will discuss in the following sections, the left disk com- of halo stars formed in situ. Kafle et al.(2017) also ponent will be considered during the analysis. With the found a similar trend with samples of Main Sequence distances from Bailer-Jones et al.(2018), the samples Turnoff and K giant stars only with radial velocity aval- within 4 kpc and relative distance error D =D < 0:2 ible. What is different is that the metal poorer K giant error are further selected so that the distance estimates are stars show a retrograde rotation. reliable, and 3827 stars are left in the sample. Com- Since Carollo et al.(2007) discovered the inner and paring the radial velocities between GDR2 and LAM- outer halo and claimed that the inner halo is rotating OST, we find an offset RV − RV of 5.38 km progradely, in contrast the outer components is rotating Gaia LAMOST s−1, which was also mentioned by Tian et al.(2015) and retrogradely. But the systematic bias in distance esti- Sch¨onrich & Aumer(2017). To avoid this systematics, Halo Rotation 3 0.8 103 pmra error 400 pmdec error 0.6 ) 300 * ( t M 2 n D 0.4 10 200 u o c count(*) 0.2 100 0 1 0.0 10 0 1 2 3 4 0.00 0.04 0.08 0.12 0.16 0.20 0 2 4 DM proper motion error RV Figure 1.