Draft version October 14, 2014 Preprint typeset using LATEX style emulateapj v. 5/2/11 THE MASSIVE SURVEY { I. A VOLUME-LIMITED INTEGRAL-FIELD SPECTROSCOPIC STUDY OF THE MOST MASSIVE EARLY-TYPE GALAXIES WITHIN 108 MPC Chung-Pei Ma1, Jenny E. Greene2, Nicholas McConnell3, Ryan Janish4, John P. Blakeslee5, Jens Thomas6 and Jeremy D. Murphy2 1 Department of Astronomy, University of California, Berkeley, CA 94720, USA; [email protected] 2 Department of Astrophysical Sciences,Princeton University, Princeton, NJ 08544, USA 3 Institute for Astronomy, University of Hawaii at Manoa, Honolulu, HI 96822, USA 4 Department of Physics, University of California, Berkeley, CA 94720, USA 5 Dominion Astrophysical Observatory, NRC Herzberg Institute of Astrophysics, Victoria, BC V9E 2E7, Canada 6 Max Planck-Institute for Extraterrestrial Physics, Giessenbachstr. 1, D-85741 Garching, Germany Draft version October 14, 2014 ABSTRACT Massive early-type galaxies represent the modern-day remnants of the earliest major star formation episodes in the history of the universe. These galaxies are central to our understanding of the evolution of cosmic structure, stellar populations, and supermassive black holes, but the details of their complex formation histories remain uncertain. To address this situation, we have initiated the MASSIVE Survey, a volume-limited, multi-wavelength, integral-field spectroscopic (IFS) and photometric survey of the structure and dynamics of the ∼ 100 most massive early-type galaxies within a distance of ∗ 11:5 108 Mpc. This survey probes a stellar mass range M & 10 M and diverse galaxy environments that have not been systematically studied to date. Our wide-field IFS data cover about two effective radii of individual galaxies, and for a subset of them, we are acquiring additional IFS observations on sub-arcsecond scales with adaptive optics. We are also acquiring deep K-band imaging to trace the extended halos of the galaxies and measure accurate total magnitudes. Dynamical orbit modeling of the combined data will allow us to simultaneously determine the stellar, black hole, and dark matter halo masses. The primary goals of the project are to constrain the black hole scaling relations at high masses, investigate systematically the stellar initial mass function and dark matter distribution in massive galaxies, and probe the late-time assembly of ellipticals through stellar population and kinematical gradients. In this paper, we describe the MASSIVE sample selection, discuss the distinct demographics and structural and environmental properties of the selected galaxies, and provide an overview of our basic observational program, science goals and early survey results. Subject headings: galaxies: elliptical and lenticular, cD | galaxies: evolution | galaxies: kinematics and dynamics | galaxies: stellar content | galaxies: structure | dark matter 1. INTRODUCTION range provides an effective tool to study the spatial The most massive early-type galaxies in the local uni- and kinematic structure, star formation histories, and verse are powerful probes of galaxy evolution. They stellar and dark matter halo masses of local galax- formed most of their stars rapidly at redshifts z > 2 ies. While many IFS surveys are ongoing, such as (e.g., Blakeslee et al. 2003; Thomas et al. 2005) but have VENGA/VIXENS (Blanc et al. 2013), CALIFA (S´anchez grown in number and size by a factor of two or more since et al. 2012), SLUGGS (Brodie et al. 2014), and eventu- z ≈ 1 (e.g., Daddi et al. 2005; Trujillo et al. 2006; Faber ally MaNGA (Bundy et al. in preparation) and SAMI et al. 2007; van der Wel et al. 2008; Damjanov et al. 2009; (Croom et al. 2012), none of them probes the volume, Cappellari et al. 2009; van Dokkum et al. 2010; van de mass range, or spatial scales required to systematically Sande et al. 2011), probably in large part through dis- study the most massive elliptical galaxies, a regime that arXiv:1407.1054v2 [astro-ph.GA] 13 Oct 2014 is critical for understanding the assembly of galaxies and sipationless merging and accretion (e.g., De Lucia et al. 3D 2006; Boylan-Kolchin et al. 2006; Naab et al. 2009; Ko- supermassive black holes. The ATLAS project (Cap- rmendy & Bender 2009; Oser et al. 2010; Thomas et al. pellari et al. 2011) surveyed 260 galaxies within 42 Mpc. Because of their relatively small survey volume, only a 2014). They contain nuclear black holes whose masses ∗ 11:5 handful of galaxies had stellar masses M ∼> 10 M . MBH are correlated with properties of the stellar bulge 00 00 (e.g., Magorrian et al. 1998; Ferrarese & Merritt 2000; Their field-of-view of 33 × 41 provided coverage within Gebhardt et al. 2000; Tremaine et al. 2002; Marconi & one half-light radius of most of their sample galaxies. Hunt 2003; H¨aring& Rix 2004; G¨ultekin et al. 2009; In this paper we describe MASSIVE, a volume-limited, Beifiori et al. 2012; McConnell & Ma 2013; Kormendy multi-wavelength, spectroscopic and photometric survey & Ho 2013). These scaling relations between black holes of the most massive galaxies in the local universe. The and their host galaxies imply co-evolution between the sample includes 116 candidate galaxies in the northern two components over the lifetime of a galaxy, but the sky with distance D < 108 Mpc and absolute K-band detailed mechanisms remain uncertain. magnitude MK < −25:3, corresponding to stellar masses ∗ > 11:5 Integral field spectroscopy (IFS) over a wide radial M ∼ 10 M . MASSIVE is designed to address a wide range of outstanding problems in elliptical galaxy 2 Ma et al. formation, including the variation in dark matter frac- TABLE 1 tion and stellar initial mass function (IMF) within and Selection criteria for MASSIVE galaxies among early-type galaxies, the connection between black hole accretion and galaxy growth, and the late-time as- Distance D < 108 Mpc Absolute K magnitude MK < −25:3 sembly of galaxy outskirts. We combine comprehensive Declination δ > −6◦ ground-based NIR imaging with IFS data to measure Galactic extinction AV < 0:6 stellar populations and kinematics out to ∼ 2 effective Morphology E and S0 radii. Using the Mitchell Spectrograph (formerly called VIRUS-P; Hill et al. 2008) at McDonald Observatory, we cover a 10700 × 10700 field of view with 400 fibers. Thus, we are sensitive to low surface-brightness emission in the outer parts of the galaxies (e.g., Murphy et al. 2011; 2. SAMPLE SELECTION Adams et al. 2012; Greene et al. 2012). 2.1. Overview For a subset of galaxies, we are obtaining additional adaptive-optics assisted IFS data to map the stellar kine- The main selection criteria of our survey are summa- matics on ∼ 100 pc scales, within the sphere of influ- rized in Table 1. The survey volume of radius D < 108 ence of nuclear black holes. The high-resolution data Mpc is chosen to be large enough to encompass the Coma are required to detect the gravitational effect of the su- cluster. This volume is more than an order of magnitude permassive black holes on the stellar orbits. These data larger than that probed by ATLAS3D , enabling us to alone, however, are insufficient for removing the degen- obtain a statistical sample of early-type galaxies at the eracy among the dark matter halo, the stellar mass-to- highest end of the galaxy mass function. The correspond- light ratio, and the central black hole mass (Gebhardt & ing redshift limit is cz < 7560 km s−1 or z < 0:025 (for −1 −1 Thomas 2009). We therefore combine the high-resolution H0 = 70 km s Mpc ). and wide-field IFS data for simultaneous modelling of Within the survey volume, our goal is to select galax- the three mass components (e.g., Schulze & Gebhardt ies above a given total stellar mass. Since selection based 2011; McConnell et al. 2011a,b, 2012; Rusli et al. 2013b; on K-band luminosities is close to a stellar mass selec- Thomas et al. 2014). tion, particularly for these red galaxies, we use the near- We are also acquiring deep K-band data at UKIRT infrared K-band magnitude from the Extended Source and CFHT to measure the extended stellar halos of these Catalog (XSC; Jarrett et al. 2000) of the Two Micron All luminous galaxies and determine more accurately their Sky Survey (2MASS; Skrutskie et al. 2006). This catalog K-band magnitudes. contains photometric measurements in the J, H, and K The wide-field IFU portion of the survey is currently bands of ∼ 1:6 × 106 objects with K ≤ 13:5 mag. The complete at MK < −25:7 and 75% complete at MK < data have a mean photometric accuracy better than 0.1 −25:5. Seven galaxies in the survey have published val- mag and are mostly unaffected by interstellar extinction ues of black hole mass (Sec 3.6); 15 additional galax- and stellar confusion, although the 2MASS luminosities ies have existing or incoming high-resolution kinematic may be systematically underestimated for very extended data. Deep K-band photometry has been obtained for objects (see Sec 2.3). 45 galaxies thus far. Peculiar velocities add uncertainties to the determina- The selection of the galaxy sample for the MASSIVE tion of distances and absolute magnitudes, and conse- survey is described in Section 2. Since this survey is quently the selection of our sample. We use the 2MASS volume-limited and defined by the stellar masses of the Galaxy Redshift Survey (2MRS; Huchra et al. 2012) and galaxies via their K-band luminosities, we discuss in de- the group catalog based on 2MRS (Crook et al. 2007) to tail the determinations of distance and absolute K−band correct the radial velocity-derived distances. We begin magnitude. Basic properties of the survey galaxies such with an initial velocity cutoff corresponding to a redshift- as stellar mass, size, velocity dispersion, shape, color, distance of 140 Mpc and MK < −22 and correct for and central black holes are presented in Section 3.