OBSERVING REQUEST University of Arizona Observatories

Year: 2006 Term: May - Aug Proposal Number:

Internal Kinematics and Stellar Populations of Compact Elliptical Discovered in the

PI: David Burstein (ASU)

CoI(s): Ronald Marzke (SFSU), Luiz da Costa, Paulo Pellegrini (Observatorio Nacional, Brazil)

Proposal Type: short-term Abstract of Scientific Justification We request one night with the MMT Blue Channel in the cross-dispersed echellette mode to obtain spectra of eleven compact elliptical galaxies discovered in the Sloan Digital Sky Survey. The high spectral resolution offered by this mode is required to obtain precise velocity dispersions of these low-mass dwarfs, which resemble the rare compact elliptical M32 at the bright end and the brightest ultracompact dwarf galaxies at the faint end. In addition, the broad wavelength coverage provides stellar absorption line strengths from NH at 3360A˚ to the Na doublet at 5890A˚ as well as nuclear emission lines from OII to SII . Together, these observations will allow the first systematic study of the formation history of compact elliptical galaxies and of the relationship between compact ellipticals and other spheroidal stellar systems.

Summary of observing runs requested for this project Run Telescope Instrument No. Nights Moon Cage Optimal months Accept. months 1 MMT Blue/Echellette 1 gray f9 May May-June 2 3 4

Scheduling constraints and non-usable dates (up to two lines). May 24-30 conflicts with KPNO 2.1m allocation for the -survey component of the program proposed here. Is sharing nights possible and/or advisable? If scheduled after 6/3/2006, two first halves would be appropriate Approval for Instrument Use from PI? For MMT requests only: Percentage of time to be assigned to UAO and CfA 1.0, 0.0 Page 2

Observing List (attach list if longer than 10 objects) Object RA DEC Magniude SDSSJ110404.4+451619 11:04:04.39 +45:16:19.1 rAB = 16.35 NGC 3979A 11:56:00.56 –02:43:12.3 rAB = 17.38 NGC 4073A 12:04:17.32 +10:15:02.8 rAB = 16.69 NGC 4073B 12:04:36.73 +01:53:33.5 rAB = 17.04 SDSSJ121456.31+473106.1 12:14:56.31 +47:31:06.1 rAB = 16.58 NGC 4621A 12:42:11.04 +11:38:41.2 rAB = 16.08 NGC 4621B 12:41:55.33 +11:40:03.7 rAB = 17.21 SDSSJ144049.71+032803.2 14:40:49.71 +03:28:03.2 rAB = 17.83 SDSSJ163214.81+480825.2 16:32:14.81 +48:08:25.2 rAB = 16.92 NGC 5846B 15:06:34.27 +01:33:31.6 rAB = 15.01

CONTACT INFORMATION: Observer(s) Institution: David Burstein Email(s): [email protected] Phone Num- ber(s) : 480-965-4336

GRADUATE STUDENTS: Whose second year project is this? Whose thesis is this?

For each graduate student named on the front page, provide the following information: Student’s Name Advisor’s Name Advisor’s Signature Page 3

Scientific Justification The observed structure of elliptical galaxies motivated early models of galaxy formation and remains one of the definitive tests of modern theory. Relationships among luminosity, surface brightness, effective radius, velocity dispersion and metallicity have played a particularly important role in the the development of the cold dark matter model (Faber, Blumenthal & Primack 1984, Dekel & Silk 1986). These scaling relations divide elliptical galaxies into two families with markedly different structural properties (Kormendy 1985, Binggeli, Sandage & Tarenghi 1984, Bender, Burstein & Faber 1992). Among giant ellipticals and spiral bulges, the faintest are the most compact, while the vast majority of dwarf ellipticals (dEs) follow the opposite trend and are, at their faintest, among the most diffuse galaxies known. The rare and extremely compact dwarf ellipticals M32, NGC 4486B, and NGC 5846A are puzzling exceptions to this rule and appear more aligned with the giant-elliptical and bulge scaling relations than with those of typical dEs (Wirth & Gallagher 1984). Because the faint bulges of late-type spirals are often obscured by bright, dusty disks, the less encumbered compact ellipticals offer unique constraints on the formation of low-mass spheroids. Although the structural similarity between compact and giant ellipticals hints at a common origin, the proximity of the few known compact ellipticals to massive companions suggests that tidal interactions play a role (King 1962), either during their formation (Burkert 1994) or over the course of their evolution (Faber 1973a). In one view (Nieto 1990, Bekki et al. 2001), M32 is interpreted as the remnant bulge of a spiral galaxy stripped of its disk by M31. With so few compact ellipticals available, the detailed measurements of their structure (Kent 1987, Choi et al. 2001, Graham 2002), stellar populations (Faber et al. 1973b, Davidge 1991, Bender et al. 1993, Worthey 2004, Rose et al. 2005 ) and kinematics (Bender et al. 1992) are difficult to place into the broader context of formation. At distances beyond the Virgo cluster, compact ellipticals are difficult to distinguish from in typical seeing conditions, and their measured abundance is thus a weak lower limit (Disney 1976). Extrapolation of the surface-brightness distribution of early-type galaxies from the Sloan Digital Sky Survey (SDSS, Gunn et al. 2006) suggests that the population of undiscovered compact ellipticals could be substantial (Shen et al. 2003). Indeed, the recent discovery of two M32 counterparts in HST/ACS imaging of Abell 1689 (Mieske et al. 2005) and the discovery of much fainter, ultra- compact dwarfs in the Fornax cluster (Phillips et al. 2001) suggest that compact stellar systems are, at the very least, under-represented in traditional galaxy surveys. Recently, Hasegan et al. (2005) and Jones et al. (2006) discovered faint, ultracompact dwarfs in the Virgo cluster, but the field abundance of compact and ultracompact dwarfs remains completely unknown. In order to compile the first statistically meaningful sample of compact elliptical galaxies, we are conducting a redshift survey of compact candidates drawn from the SDSS/DR4 (Adelman- McCarthy et al., 2006). We have selected bright sources with r-band Petrosian half-light radii 00 r50 < 1.6 and ugriz colors consistent with low-redshift galaxies (Fig 1). Most of our candidates are rejected as possible stellar contaminants during spectroscopic targeting of the primary SDSS galaxy samples (Strauss et al. 2002, Eisenstein et al. 2001), but we have identified a subsample of nine galaxies with SDSS that are confirmed compact ellipticals. Remarkably, one of these lies in the halo of the field elliptical NGC 5846, which already hosts the prototypical compact elliptical NGC 5846A. Two of the other confirmed candidates have no sizable companions and thus provide unique constraints on the tidal origin of compact elliptical galaxies. Among the sources rejected by the SDSS targeting algorithm, we have identified two high-probability candidates in the halo of the Virgo cluster elliptical NGC 4621. One of these candidates appears on an archival HST/ACS image, and its compact surface brightness profile places it between the brighter compact elliptical prototypes and the Fornax ultracompact dwarfs observed with HST/STIS by de Propris et Page 4 al. (2005). Already, our sample nearly doubles the number of previously known compact ellipticals, and another 124 candidates are targeted for spectroscopy at the KPNO 2.1m in May 2006. Here, we request one night with the MMT Blue Channel in its echellette mode to obtain high- resolution spectra of eleven compact elliptical galaxies with absolute magnitudes ranging from Mr = −13.9 to Mr = −18.3. Nine of these galaxies are spectroscopically confirmed in the SDSS, and the other two are the newly-discovered companions to NGC 4621. Our goals are to measure precise velocity dispersions, absortion-line strengths, and (in some cases) emission-line ratios of these galaxies in order to compare their kinematics, stellar populations and nuclear activity to those of other spheroidal stellar systems. We have submitted a complementary proposal to obtain surface brightness profiles and color gradients of these same galaxies using the ACS High Resolution Camera on HST. We outline our strategy below. Kinematics: Correlations among luminosity, surface brightness and velocity dispersion constrain the dark-matter content of hot stellar systems (Burstein et al. 1997). At the high-mass end, giant 0.25 ellipticals and bulges follow the Faber-Jackson relation, σ0 ∝ L (Faber & Jackson 1976). At the low-mass end, globular clusters appear to contain little dark matter and follow a steeper luminosity- 0.5 velocity dispersion relation: σ0 ∝ L (McLaughlin & van der Marel 2005). Little is known about the mass range between these extremes. The internal kinematics of seven faint, ultracompact dwarfs (UCDs) in the Virgo cluster align more with massive globular clusters (Hasegan et al. 2005), while a handful of brighter UCDs in Fornax appear to follow an extrapolation of the Faber- Jackson relation (Drinkwater et al. 2003). Although brighter UCDs have been identified as possible remnants of stripped, nucleated dwarf ellipticals, their kinematics may indicate a closer kinship with compact elliptical galaxies like M32. Precise velocity dispersions of our sample of compact ellipticals will begin to fill the essentially unexplored region of the fundamental plane between M32 and the ultracompact dwarfs and will thus constrain the link between globular clusters and galaxies. Figure 2 shows simulated spectra of our targets, which require both high resolution and high SNR for precise velocity dispersion measurements (see Experimental Design). Line Strengths: Faber (1973a) notes that the compact ellipticals M32, NGC 4486B and NGC 5846A have heavy-element abundances characteristic of more luminous galaxies and suggests that compact ellipticals are the stripped cores of brighter ellipiticals. On the other hand, Nieto (1990) and Bekki et al. (2001) suggest that M32 was once the bulge of a spiral galaxy. The nuclear profile of M32 is a classic power-law (Lauer et al. 1998), while NGC 4486B has a core that is larger than expected given the low luminosity of NGC 4486B (Kormendy et al. 1997). Together, the core properties and stellar populations of known compact ellipticals suggest that their progenitors may be diverse. In spiral bulges, Thompson & Davies (2006) show that luminosity-weighted age, total metallicity and α/Fe ratio correlate strongly with central velocity dispersion and that these quantities coincide with those of elliptical galaxies at fixed mass. Echellete spectra, which would cover the range from NH at 3360Ato˚ the Na doublet at 5890A,˚ would allow us to place our compact elliptical candidates on these relations and thus to constrain the mass of their progenitors. Nuclear Activity: The spectrum of one of our targets is dominated by Balmer emission from an HII nucleus (Fig. 3). The proposed echellette spectra would include the full Balmer series as well as lines of OII, OIII, NII and SII. The high SNR would allow detection and classification of weak nuclear emission in the other 10 sources along the lines of Ho, Fillipenko & Sargent (1997).

REFERENCES: Bender,R., Burstein,D., & Faber,S. 1992, ApJ, 399, 462 • Bekki,K., Couch,W., Drinkwater,M., & Gregg, M. 2001, ApJ, 557, L39 • Binggeli,B., Sandage, A., & Tammann,G., 1985, AJ, 90, 1681 • Burkert, A., 1994, MNRAS, 266, 877 • Burstein, D., Bender, R., Faber, S., & Nothenius, R., 1997, AJ, 114, 1365 • Choi,P., Guhathakurta, R., & Johnston, K. 2002, AJ, 124, 310 • Cˆot´e,P. et al. 2004, ApJS, 153, 223 • Davidge, T. 1991, AJ, 102, 896 • Davies, R.L. & Thomas, D., 2006, MNRAS, 366, 510 • de Propris, R., et al., 2005, ApJ, 623, L105 • Disney, M., 1976, Nature, 263, 573 • Drinkwater, M.J., et Page 5

al., 2003, Nature, 432, 519 • Eisenstein, D., et al., 2001, AJ, 122, 2267 • Faber, S.M. 1973a, ApJ, 179, 423 • Faber, S.M. 1973b, ApJ, 179, 731 • Graham, A., 2002, ApJ, 568, 13 • Gunn, J., et al., 2006, AJ, in press • Hasegan, M. et al. 2005, ApJ, 627, 203 • Kent, S.M., 1987, AJ, 94, 306 • King,I. 1962, AJ, 67, 471 • Kormendy,J. 1985, ApJ, 295, 73 • Lauer, T. et al. 1998, 116, 2263 • Mieske, S. et al., 2005, A&A, 430, L25 • 258 • Phillips,S., et al., 2001, ApJ, 560, 201 • Shen, S. et al., 2003, MNRAS, 343, 978 • Strauss, M.A. et al., 2002, AJ, 124, 1810 • Wirth,A. & Gallagher, J. 1984, ApJ, 282, 85

3.0

2.5 150

2.0

u - g 100 1.5 Photons per pixel 1.0 50 0.5 -0.2 0.0 0.2 0.4 0.6 5160 5170 5180 5190 i - z Wavelength (Angstroms) Figure 1: SDSS ugiz colors of previously known Figure 2: Simulated 1200s MMT/Blue Echellette spectra compact elliptical galaxies (the three red squares), of an rAB = 17.8 stellar template (black) broadened with bright elliptical galaxies from the RC3 (crosses), dif- 30 km s−1 (purple) and 75 km s−1 (red) Gaussian velocity fuse Virgo dEs (open circles), and sources identified distributions. The small region highlighted here includes as stars by the SDSS photometric pipeline (colored the Mg b triplet and is drawn from the tenth echelle order. density contours). The lines indicate color cuts im- In order to highlight the line widths, we do not show the posed on our sample of compact candidates. noise in this figure (but include it in the simulations).

Figure 3: SDSS gri images and spectra (where available) of the 11 compact elliptical candidates targeted for MMT spectroscopy. For comparison, we include in the lower right corner the SDSS image of one of the three prototypical compact elliptical galaxies, NGC 5846A. Page 6

Experimental Design & Technical Description Describe your overall observational program. How will these observations contribute toward the accomplishment of the goals outlined in the science justification? If you’ve requested long-term status, justify why this is necessary for successful completion of the science. (limit text to one page) Our choice of instrument is driven by two conflicting needs: high spectral resolution (for precise velocity dispersions of low-mass galaxies) and broad wavelength coverage (for the analysis of ab- sorption features). The velocity dispersion of M32 (which is comparable to the typical galaxy in our sample) is 75 km s−1. Our faintest candidate is similar to the brightest ultracompact dwarf in Fornax, which has σ = 37 km s−1. Measurement of velocity dispersions below approximately 90 km s−1 requires higher-resolution spectra than are available from the SDSS. The Blue Channel in its cross-dispersed echellette mode yields ∼30 km s−1 resolution (100 slit) over narrow wavelength ranges while retaining nearly full wavelength coverage between approximately 3300A˚ and 8100A.˚ The velocity dispersion measurements set the resolution and signal-to-noise requirements of the program. We expect the velocity dispersions of our lowest-mass candidates to be greater than 30 km s−1. We have simulated spectra of these low-mass galaxies using a high-resolution solar spectrum broadened with a Gaussian velocity distribution and convolved with the instrumental profile of the Blue Channel echellette (assuming a 100 slit). In most cases, approximately half the total light from these compact sources is lost at the slit. The throughput in the cross-dispersed mode is poorly known , but the relative blaze functions can be measured from a publicly-available spectrum of a flux standard (taken through clouds). During engineering runs, the throughput of the echelle grating was estimated to be more than a factor of two less sensitive than other gratings at the peak of their blaze. We have therefore used the flux standard to establish the relative transmission and then normalized it using the published throughput of the 1200l grating at 5000A˚ in first order (dividing by a factor of 2.5 to mimic the lower throughput of the echelle). Figure 2 shows the mean number of photons expected in a 20-minute exposure for a galaxy with rAB = 17.8 and two different velocity dispersions: σ = 30 km s−1 and σ = 75 km s−1. In order to estimate the expected uncertainties in our velocity dispersions, we simulate 100 spectra per exposure time including noise from the source, the sky and the readout (to facilitate cosmic-ray rejection, we take three exposures per target, and the effective read nosie is 4.8e−). At each exposure time, we extract the 10th echelle order (which includes Mg b) and then calculate velocity dispersions using the IRAF package RVSAO and a simulated stellar template. The standard deviation of the 100 measured velocity dispersions provides a first estimate of the uncertainties at each exposure time. In 20-minute exposures, the SNR per 0.19A˚ pixel exceeds 30 at the Mg b triplet, and the uncertainty in velocity dispersion is 2.8 km s−1 at σ = 30 km s−1 and 8.4 km s−1 at σ = 75 km s−1. In practice, we can improve the SNR by combining multiple echelle orders, but given that real measurements will also be affected by focus variations and template mismatch, we choose to set the exposure time using a single order and then use other orders to check for systematic errors. The broad wavelength coverage allows measurement of absorption-line strengths from NH at 3360A˚ to the Na doublet at 5890A.˚ These measurements can be made at much lower resolution and are therefore less compromised by the unavoidable focus variations across the CCD. For this analysis,the data can also be rebinned in the dispersion direction to achieve very high SNR even in the near-UV. Table 1 lists the eleven targets along with their physical properties. Distances and luminosities −1 −1 assume H0 = 70 km s Mpc . The column labeled Offset gives the distance to the nearest bright neighbor, and the last column identifies the neighbor. We can obtain 20 minutes of total exposure time for each target in a total of six hours (including readout time and telescope overhead). For calibration, we will also observe several standards as well as three bright galaxies with well-measured line indices. In total, these observations will require one night at the MMT. Page 7

Table 1. Compact Elliptical Galaxy Candidates 00 ID Distance (Mpc) Mr g − r r50 ( ) r50 (kpc) Offset (kpc) Neighbor 1 66 -17.5 0.85 1.4 0.43 27 NGC 677 2 82 -17.2 0.61 0.8 0.32 2 NGC 3979 3 85 -18.0 0.77 1.3 0.53 67 NGC 4073 4 82 -17.5 0.77 1.4 0.55 55 NGC 4073 5 27 -17.3 0.77 1.5 0.20 24 NGC 5846 6 113 -17.5 0.75 1.0 0.52 70 NGC 5718 7 97 -18.3 0.69 1.4 0.64 490 NGC 4231 8 87 -18.3 0.81 1.3 0.52 47 MCG +08-20-077 9 98 -18.0 0.70 1.3 0.59 220 SDSS J163136.1+480320 10 17 -15.1 0.78 0.8 0.07 10 NGC 4621 11 17 -13.9 0.78 0.9 0.08 13 NGC 4621

Previous Use of Steward Facilities List allocations of telescope time on facilities available through Steward to the investigators during the past 2 years, together with the current status of the data (cite publi- cations where appropriate). Mark with an asterisk those allocations of time related to the current proposal. Burstein et al. 2004, ApJ, 614, 158, “Globular Cluster and Galaxy Formation: M31, The Milky Way, and Implications for Globular Cluster Systems of Spiral Galaxies”, based on MMT observa- tions from 9/27-28/2003. 11/24-25/2005 MMT/Blue Channel: Burstein, Li & Dolan, “From Whence do Old Stellar Popu- lations Come?” We took fainter M31 globular clusters, comparable in luminosity to Galactic GCs, and find the fainter M31 GCs have comparable NH absorption feature as do the Galactic GCs.

Revised version of UAO Observing Proposal LATEX macros v2.5.