Arxiv:1612.06398V1 [Astro-Ph.GA] 19 Dec 2016 Thousands of Objects Per System (E.G., Battaglia Et Al
Total Page:16
File Type:pdf, Size:1020Kb
Draft version October 1, 2018 Preprint typeset using LATEX style emulateapj v. 5/2/11 CRATER 2: AN EXTREMELY COLD DARK MATTER HALO Nelson Caldwell1, Matthew G. Walker2, Mario Mateo3, Edward W. Olszewski4, Sergey Koposov5, Vasily Belokurov5, Gabriel Torrealba5, Alex Geringer-Sameth2 and Christian I. Johnson1 Draft version October 1, 2018 ABSTRACT We present results from MMT/Hectochelle spectroscopy of 390 red giant candidate stars along the line of sight to the recently-discovered Galactic satellite Crater 2. Modelling the joint distribution of stellar positions, velocities and metallicities as a mixture of Crater 2 and Galactic foreground populations, we identify 62 members of Crater 2, for which we resolve line-of-sight velocity dispersion +0:3 1 ∼ +0:4 1 σvlos =2:7 0:3 km s− about mean velocity of vlos =87:5 0:4 km s− (solar rest frame). We also − +0:04h i − +0:1 resolve a metallicity dispersion σ[Fe=H]=0:22 0:03 dex about a mean of [Fe/H] = 1:98 0:1 dex that is 0:28 0:14 dex poorer than is estimated− from photometry. Despiteh Crateri − 2's relatively− large ± size (projected halflight radius Rh 1 kpc) and intermediate luminosity (MV 8), its velocity dispersion is the coldest that has been∼ resolved for any dwarf galaxy. These properties∼ − make Crater 2 the most extreme low-density outlier in dynamical as well as structural scaling relations among the Milky Way's dwarf spheroidals. Even so, under assumptions of dynamical equilibrium and negligible contamination by unresolved binary stars, the observed velocity distribution implies a gravitationally +1:2 6 +15 dominant dark matter halo, with dynamical mass 4:4 0:9 10 M and mass-to-light ratio 53 11 − × +0:5 − M =LV; enclosed within a radius of 1 kpc, where the equivalent circular velocity is 4:3 0:5 km 1 ∼ − s− . Subject headings: galaxies: dwarf | galaxies: individual (Crater 2) | (galaxies:) Local Group | galaxies: kinematics and dynamics | methods: data analysis | techniques: spec- troscopic 1. INTRODUCTION parameters from photometric and spectroscopic surveys, The Milky Way's dwarf-galactic satellites represent an respectively, reveals empirical scaling relations that serve extremum of galaxy formation. They include the small- as targets for models of dwarf galaxy formation and evo- est, least luminous, least chemically enriched and most lution. For example, the size/luminosity relation traced dark-matter dominated galaxies known. Their abun- by dwarf spheroidals extends the one observed for disk dances as functions of size, luminosity, metallicity and galaxies but is offset from relations traced by central dark matter density are sensitive to the nature of dark bulges and larger elliptical galaxies (Kormendy & Bender matter (e.g., the ability of dark matter particles to cluster 2012). Moreover, the least luminous dwarf galaxies also on scales of tens of pc) as well as to how galaxy for- exhibit the lowest metallicities, extending the galactic ∼ luminosity/metallicity relation over 5 orders of mag- mation proceeds within low-mass dark matter halos (e.g., ∼ Pontzen & Governato 2014; Di Cintio et al. 2014; Read nitude in luminosity (Kirby et al. 2013). Finally, corre- et al. 2016). Empirical information has increased dra- lations among internal line-of-sight velocity dispersion, matically during the last decade, as deep imaging surveys size, and luminosity have been interpreted as evidence have grown the number of known Galactic satellites from for similarity among the dark matter halos that dwarf 10 to several tens (e.g., Belokurov et al. 2007; Koposov galaxies inhabit (Mateo et al. 1993; Strigari et al. 2008; et∼ al. 2015; Bechtol et al. 2015; The DES Collaboration Walker et al. 2009b). et al. 2015) and spectroscopic campaigns have measured While it is clear that dwarf galaxies' structural and stellar kinematics and chemical abundances for tens to chemodynamical properties are interconnected, physical interpretations are complicated by our ignorance regard- arXiv:1612.06398v1 [astro-ph.GA] 19 Dec 2016 thousands of objects per system (e.g., Battaglia et al. 2006; Martin et al. 2007; Simon & Geha 2007; Walker, ing the extent to which these empirical correlations are Mateo & Olszewski 2009; Kirby et al. 2010). driven by selection effects. At a given distance and lumi- The combination of structural and chemodynamical nosity, smaller galaxies have higher surface densities and hence are detected more easily than are larger galaxies. [email protected] Careful quantification of sensitivities for various imag- 1 Harvard-Smithsonian Center for Astrophysics, 60 Garden ing surveys (e.g., Koposov et al. 2008) leaves open the Street, Cambridge, MA 02138, USA possibility that large numbers of large, low-luminosity 2 McWilliams Center for Cosmology, Department of Physics, galaxies exist undetected, perhaps modifying what are Carnegie Mellon University, 5000 Forbes Ave., Pittsburgh, PA 15213, United States apparent size/luminosity correlations and/or other rela- 3 Department of Astronomy, University of Michigan, 311 West tions involving quantities that correlate with size and Hall, 1085 S. University Ave., Ann Arbor, MI 48109 4 luminosity. Steward Observatory, The University of Arizona, 933 N. In fact dwarf galaxies continue to be discovered down Cherry Ave., Tucson, AZ 85721 5 Institute of Astronomy, University of Cambridge, Madingley to the faintest surface brightnesses to which current Road, Cambridge, CB3 0HA, United Kingdom imaging surveys are sensitive. Most recently, (Torrealba 2 Caldwell et al et al. 2016, T16 hereafter) used data from the ATLAS east (see right panel of Figure 1). Exposure times ranged survey (Shanks et al. 2015) to discover Crater 2, a dwarf between 2 3 hours per field. Table 1 provides details galaxy at distance D 117 kpc from the Sun. With ab- of the exposures,− which provided spectra for 610 unique ∼ solute magnitude MV 8:2 and projected halflight ra- science targets. In each field, 40 fibers observe regions dius R 1:1 kpc, Crater∼ − 2 exhibits the faintest surface of blank sky in order to facilitate∼ sky subtraction. h ∼ brightness of any known dwarf galaxy, with µV 31 mag In addition to science exposures, we also acquired ex- 2 ∼ arcsec− . Relative to other Galactic satellites of similar posures of ThAr arc lamps, taken immediately before and luminosity (or size), Crater 2 is exceptionally large (or after science exposures, for the purpose of wavelength underluminous). calibration. At the beginning of each night, we also ob- In order to gauge how this outlier status translates into tained twilight sky exposures in order to correct for fiber chemodynamical scaling relations, and in order to study throughput variations and calibrate our measurements the dark matter content of Crater 2 in general, here we of stellar-atmospheric parameters. For all exposures, the present results from initial MMT/Hectochelle spectro- Hectochelle CCDs were binned by 3 pixels in the spectral scopic observations of 390 stars along the line of sight direction and by 2 in the spatial direction. The effective to the galaxy. From the spectroscopic data, we estimate resolution is R 32; 000. All data frames were processed line-of-sight velocity as well as effective temperature, sur- using the TDC/Hectochelle≈ pipeline, which provides 1D, face gravity and metallicity for each individual star. We throughput-corrected, wavelength-calibrated spectra, as then use these measurements to characterise Crater 2's well as their corresponding variance spectra, as described internal chemodynamics, to study its dark matter con- in detail by Caldwell et al. (2009). tent, and to place Crater 2 in the context of the popula- Given the processed spectra for a given exposure, we tion of known Milky Way satellites. follow the procedure of Koposov et al. (2011) to estimate the mean sky spectrum and to subtract it from each sci- 2. OBSERVATIONS AND DATA REDUCTION ence spectrum. All processed and sky-subtracted spectra We observed Crater 2 using the Hectochelle multi- for this program are publicly available and included in object fiber spectrograph (Szentgyorgyi et al. 2011) on the online database associated with this article. the 6.5m MMT on Mt. Hopkins, Arizona. Hectochelle's We analyze each sky-subtracted spectrum following field of view subtends 1◦ in diameter and thus is well- the procedure of Walker et al. (2015). Briefly, we fit a matched to Crater 2's projected halflight radius of R smoothed library of synthetic template spectra, obtain- h ∼ 310 (T16). Hectochelle has been widely used for the study ing Bayesian inferences for the line-of-sight velocity (vlos) of MW dwarf galaxies (e.g., Mateo et al. 2008; Belokurov as well as stellar-atmospheric parameters including effec- et al. 2009; Walker et al. 2009a), and achieves velocity tive temperature (Teff ), surface gravity (log g) and metal- 1 precision better than 1 km s− (Walker et al. 2015). licity ([Fe/H]). As in previous work, we use the synthetic library originally generated for the SEGUE Stellar Pa- 2.1. Target Selection rameter Pipeline (Lee et al. 2008a,b, `SSPP' hereafter), which contains continuum-normalized spectra computed In order to identify stars for spectroscopic targeting, over a regular grid in Teff , log g and [Fe/H] and assumes we selected red giant branch (RGB) candidates from the a piecewise-linear relation between [α/Fe] and [Fe/H]. photometric catalog provided by the VLT/ATLAS sur- We perform the fits using the nested-sampling algo- vey (Shanks et al. 2015). The left panel of Figure 1 shows rithm MultiNest6 (Feroz & Hobson 2008; Feroz et al. the color-magnitude diagram for point sources near the 2009), which returns random samples from the posterior center of Crater 2. The overplotted isochrone (Dotter probability distribution function. From these samples, et al. 2008) corresponds to age=12 Gyr, [Fe/H]= 1:7 we record the mean, variance, skewness and kurtosis that and distance modulus m M = 20:3, closely following− − summarize the 1D posterior PDF for each free parame- the ridge line identified by T16.