Fossil Evidence for the Two-phase Formation of Elliptical Galaxies
Massive early-type galaxies have undergone dramatic structural evolution over the last 10 Gyr. A companion paper shows that nearby elliptical galaxies with generically contain three photometric subcomponents: a compact inner component with effective radius 1 kpc, an intermediate-scale middle component with kpc, and an extended outer envelope with kpc. Here we attempt to relate these substructures with the properties of early-type galaxies observed at higher redshifts. We find that a hypothetical structure formed from combining the inner plus the middle components of local ellipticals follows a strikingly tight stellar mass-size relation, one that resembles the distribution of early-type galaxies at . Outside of the central kpc, the median stellar mass surface density profiles of this composite structure agree closest with those of massive galaxies that have similar cumulative number density at within the uncertainty. We propose that the central substructures in nearby ellipticals are the evolutionary descendants of the “red nuggets” formed under highly dissipative (“wet”) conditions at high redshifts, as envisioned in the initial stages of the two-phase formation scenario recently advocated for massive galaxies. Subsequent accretion, plausibly through dissipationless (“dry”) minor mergers, builds the outer regions of the galaxy identified as the outer envelope in our decomposition. The large scatter exhibited by this component on the stellar mass-size plane testifies to the stochastic nature of the accretion events.
Recent observations have established that high-redshift early-type galaxies (ETGs) are more compact (Daddi et al. 2005; Trujillo et al. 2006; Damjanov et al. 2011) and have higher stellar velocity dispersions (Cappellari et al. 2009; Onodera et al. 2009) than nearby ETGs of the same stellar mass. Since , “red nuggets” on average have doubled in stellar mass and increased their size by a factor of (Buitrago et al. 2008; van Dokkum et al. 2010; Papovich et al. 2012; Szomoru et al. 2012). Some appear to have a disk-like morphology (van der Wel et al. 2011), and they are always bluer in the outskirts (Gargiulo et al. 2012). The accumulated evidence suggests that massive ETGs build up inside-out through non-dissipational processes (Bezanson et al. 2009; van Dokkum et al. 2010).
These developments severely challenge classical models of elliptical (E) galaxy formation such as monolithic collapse (Larson 1975) and binary major mergers (Toomre & Toomre 1972; Negroponte & White 1983). Instead, the current evidence favors a “two-phase” scenario (Oser et al. 2010; Johansson et al. 2012). Intense dissipational processes such as cold accretion (Dekel et al. 2009) or gas-rich mergers rapidly build up an initially compact progenitor. After star formation is quenched, a second phase of slower, more protracted evolution is dominated by non-dissipational processes such as dry, minor mergers.
For such a scenario to work, the rate of minor mergers (Bluck et al. 2012; Newman et al. 2012) has to be consistent with the prevalence of faint companion galaxies at high redshifts. Perhaps other mechanisms may also be important, such as major merger (Bernardi et al. 2011; Prieto et al. 2013) or AGN-induced expansion (Fan et al. 2008). Given that large uncertainties remain in identifying the relevant physics and quantifying their detailed balance, a key sanity check is to see whether one can separate out the fossil bodies of compact galaxies at the core of nearby galaxies. This separation is successful only when several clear and solid predictions, especially the similarity between the core components and the high- compact progenitors in mass versus size relation and detailed mass profile, are met. These stringent predictions are the basis of this study to verify or falsify. We assume a CDM cosmology with , , and kms Mpc.
2 Observational Material
This work uses the three-component models of nearby Es from the Carnegie-Irvine Galaxy Survey (CGS; Ho et al. 2011), a systematic study of 605 bright ( mag) southern () galaxies. The currently completed photometric part of CGS includes BVRI images obtained using the 100-inch du Pont telescope at Las Campanas Observatory, using a 8989 detector with a pixel scale of , under 1 seeing conditions. More details on the observations and data reductions can be found in Ho et al. (2011) and Li et al. (2011).
Huang et al. (2013) apply GALFIT (Peng et al. 2010) to perform a detailed two-dimensional analysis of the -band images of 94 CGS Es. Their decomposition reveals that, contrary to popular perception, the global light distribution of most nearby Es have more complicated structures than can be described by a single Sérsic (1968) profile. Instead, 70 out of the 94 Es can be fit with three distinct subcomponents: (1) an inner, compact center with effective radius 1 kpc comprising a small fraction of the light (); (2) an intermediate-scale, middle structure with kpc and ; and (3) an outer, dominant, extended, moderately flattened envelope with kpc and .
The image decomposition was performed without ascribing any physical significance to the substructure. Nevertheless, the fact that the individual subcomponents define distinct sequences on the size-luminosity plane, and that three-component fits often seemed naturally preferred, hint at the possibility that the subcomponents reflect underlying physical reality. Here we compare some key properties of the three-component models of local Es with observations of high- ETGs to shed new light on the structural evolution of massive galaxies.
For this purpose, we compiled stellar masses and effective radii for 1323 high- ETGs with stellar masses , as follows: 352 () from Damjanov et al. (2011); 910 () from COSMOS (Scoville et al. 2007; using morphological classification from Scarlata et al. 2007 and Tasca et al. 2009); 32 (, choosing Sérsic indices ) from the GOODS-NICMOS survey (Conselice et al. 2011); 8 (, ) from Papovich et al. (2012); and 21 () from CANDELS (Grogin et al. 2011, Koekemoer et al. 2011), as analyzed by Szomoru et al. (2012). Although the sample is in many respects heterogeneous (in terms of observed bandpass, selection criteria, modeling method), the general trends we explore in this paper should not be strongly affected, as demonstrated by Damjanov et al. (2011).
3 Linking Nearby and High- Early-type Galaxies
Figure 1(a) illustrates the redshift evolution of galaxy size for massive ETGs. Although the scatter is substantial, the size increase with decreasing is confirmed (e.g., Damjanov et al. 2011; Cimatti et al. 2012). The median circularized effective radius increases from kpc at , to kpc at , to kpc at . By , Es have reached a median kpc according to the single-component fits applied to the subsample of 35 massive CGS Es that satisfy (solid black points). The CGS mass cut roughly corresponds to twice the mass limit of the high- sample.
How are the subcomponents of nearby Es related to the high- systems? Judging by the distribution of sizes for the subcomponents shown on the left side of Figure 1(a), the typical effective radius of the inner component (median kpc) does overlap substantially with that of high- red nuggets; however, their stellar masses, only 10%–15% of the total, are about a factor of 3–5 lower than the high- systems (Figure 1b). The middle component, with median kpc, matches well the sizes and masses of ETGs at , but they are still on average more extended than those at . On the other hand, with median kpc, the outer component of local Es is clearly larger than the globally averaged size of most ETGs at any redshift.
The mass-size relation offers another vantage point. As described in Huang et al. (2013), single-component fits to the CGS Es produce sizes that follow well the mass-size relation of local ETGs selected from the Sloan Digital Sky Survey (SDSS; Guo et al. 2009, dashed line in Figure 1b). As expected, high- ETGs, especially those with , fall systematically below the local relation. At fixed stellar mass, high- systems are a factor more compact than their counterparts. While the inner component of the CGS galaxies is generally not massive enough to directly map onto high- objects, the middle component can be viewed as a slightly modified version of the compact quiescent galaxies. The mass-size relation for the middle component runs roughly parallel to, but is offset slightly below, that of local ETGs. It lies approximately in between the locus of high- and low- points. In dramatic contrast with the inner and middle components, the outer envelope displays significantly larger scatter on the plane. While some of the scatter undoubtedly arises from measurement uncertainty, which, as indicated at the bottom-right corner of Figure 1(b), is considerable for extended, low-surface brightness features, we believe that most of the scatter is intrinsic.
As discussed in Huang et al. (2013), the central core may have been reshaped more recently (Hopkins et al. 2009a, b, c) than the overall compact stellar body of interest, thus requiring a separate component to model well in order to sensitively extract the lower lying outer component. Regardless of the exact process that shaped the small central structure, it is a small perturbation to the overall structure of interest, hence we group the innermost two components (designated henceforth as innermiddle) into a single entity. To measure the overall size, mass, and profile of the composite, we generated model images from best fit parameters that summed the innermost subcomponents, accounting for the PSF. Then, we fit a single-component model to extract the effective radius of this composite structure. As the inner component has higher density than the middle one, we intuitively expect the effective radius of the combined structure to be smaller than that of the middle component alone. This is indeed the case, as Figure 1(b) shows. Surprisingly, the composite innermiddle structure (orange diamonds) not only forms a very tight mass-size relation with a slope similar to that of Guo et al. (2009) for SDSS ETGs, but it also lies on top of the distribution of points for most high- ETGs, especially those with .
As advocated by van Dokkum et al. (2010), galaxies selected at constant cumulative number density provide a more physical connection between progenitors and descendants than those selected above a constant stellar mass. We therefore compare our observations with the sample of Patel et al. (2013), which is selected at a cumulative number density ( ) that corresponds to massive ETGs () similar to the CGS sources at . Figure 1(b) shows that the median and of our innermiddle components agree quite well with the properties of galaxies at from Patel et al. (2013). Apart from global quantities such as mass and effective radius, the surface density profiles themselves provide a more direct comparison of the structural connection between high- ETGs and local Es. For this purpose, we convert the surface brightness profiles of the innermiddle components, after removing the point-spread function smearing, into stellar mass density profiles using the mass-to-light ratio adopted in Huang et al. (2013). Their median profile is shown in Figure 2(a), and in Figure 2(b) we compare it to the median profiles of the cumulative number density-selected galaxies in different redshift bins from Patel et al. (2013). The median profile of the innermiddle components of local Es matches quite closely the profiles of massive galaxies at over most of their radial extent. Unlike the median and , the median profile of massive galaxies only marginally close to the profile of innermiddle component. This seemingly difference relates to the higher density of high- objects at kpc. It is confirmed in the median density profile of the massive ETGs at from Szomoru et al. (2012), which show even higher central densities than the sample from Patel et al. at comparable redshifts. This apparent discrepancy may be due to different selection criteria employed in the two studies.
The two-phase formation scenario of ETGs (Oser et al. 2010; Johansson et al. 2012) offers an appealing framework to interpret the above trends and to connect the distant and local populations of massive galaxies. Dissipative processes (“in situ” star formation, in the language of Oser et al. 2010) connected with early gas-rich events naturally lead to high central densities. We suggest that the central, more compact structures of nearby Es—namely the “inner” or “middle” components, possibly both—are the natural outcomes of this initial, dissipative phase. We propose that they are the direct remnants of the high- massive quiescent galaxies. The combined luminosity (mass) fraction of the innermiddle components, are roughly comparable to the factor of increase necessary to grow the high- objects, if they are the progenitors of local Es. Their overall physical scales are also consistent with their originally compact nature. In detail, the current evidence suggests that the high- objects tend to have somewhat higher central densities than their local remnants. This is an important constraint that future theoretical models should try to reproduce.
We equate the extended outer envelope with material accreted through dissipationless minor mergers (“ex situ” stars; Oser et al. 2010). The large intrinsic scatter on the plane is likely related to the stochastic nature by which the outer envelope is expected to have accumulated. The minor merger hypothesis has previously been questioned on the basis of the apparently small observed scatter of the mass-size relation of nearby Es (Nipoti et al. 2012). For those two observations to be consistent, dry mergers would have to preferentially add material to large radii, realized through many stochastic events spread over a protracted period. We propose that this history is represented by the extended, outer component identified through our photometric decomposition.
We end by briefly comparing our decomposition results with predictions from recent numerical simulations for the formation of massive ETGs. After separating the final stellar system into its in situ (dissipative) and ex situ (accreted) constituents in their simulations, Oser et al. (2010) show that, in comparison to the in situ material, the typical size of the accreted component should correlate with the total stellar mass of the system. Furthermore, the stellar mass fraction of the accreted component may also correlate with the total stellar mass and effective radius of the whole system (Oser et al. 2010; Johansson et al. 2012; Lackner et al. 2012). However, their numerical simulations produce a wide range of possible outcomes: the predictions for the variation of accreted fraction with total stellar mass in Oser et al. (2010) disagree strongly with those from the work of Lackner et al. (2012).
We compare our analysis with the simulations to look for correlations in the accreted components. We assume that the outer and inner+middle components of our multi-component fits approximate the “accreted” and “in situ” material in the theoretical models. Figure 3(a) shows that the effective radius of the outer component of nearby Es does indeed increase strongly with the total stellar mass of the system, much more so than either of the inner or middle components. While the physical scales and the slope of the correlation are not the same as predicted in the simulations of Oser et al. , there is qualitative agreement. However, the comparison between the accreted fraction, approximated as /, and effective radius (Figures 3b) and total stellar mass (Figures 3c) do not show much, if any, correlation. Apart from the uncertainties of the predictions of current simulations, it is a difficult observational challenge to isolate the accreted component via photometric decomposition. The details of any decomposition are dependent on the choice of specific parameterized model components, which are not unique choices. Also, further galaxy merging would have diluted and partly erased the signatures.
4 Summary and Discussion
Huang et al. (2013) show that the surface brightness distribution of nearby Es can best be described by three photometrically distinct substructures, consisting of a compact ( 1 kpc) inner component, an intermediate-scale ( kpc) middle component, and an extended ( kpc) outer envelope. Motivated by the recently proposed two-phase formation scenario for massive galaxies, here we argue that the inner and middle components identified in nearby Es jointly comprise a structure that can be considered the evolved, local counterparts of the high- compact massive galaxies. In this scenario, these compact structures derived from highly dissipative processes at high redshifts. Dissipationless minor mergers dominate the late-time evolution of ETGs by contributing to their dramatic size growth. We identify the outer, extended component of nearby Es with this late accretion phase.
Several lines of investigation deserve further attention. In light of the photometric subcomponents uncovered in nearby Es, it is of interest to ask whether such substructure already exists in high- ETGs, and, if so, whether their gradual build up can be detected as a function of time. Current studies treat high- ETGs as single-component systems (e.g., Papovich et al. 2012; Szomoru et al. 2012), but the data may have sufficient quality to allow more complex analysis to reveal potential substructure. Of particular interest is the mechanism by which the “core” structure in massive Es (e.g. Lauer et al. 1995; Kormendy et al. 2009) arise. If scouring by binary supermassive black holes is necessary to make cores (e.g., Hopkins et al. 2009b, c), then major mergers are implicated. Similarly, if high- massive galaxies contain disky or fast-rotating structures (van der Wel et al. 2011), then they must be destroyed, presumably also through major mergers, if they are to evolve into local massive Es. How do major mergers fit into the framework of the two-phase formation scenario?
Beyond imprints on the photometric structural properties, the two-phase formation scenario also implies differences in stellar population. Recent spectroscopic studies of the faint, outer regions of nearby Es already hint that they have lower metallicity and possibly older age than the inner parts (Coccato et al. 2010; Greene et al. 2012). Photometry is vastly simpler and less expensive than spectroscopy, even if broad-band colors provide less robust constraints on stellar populations than spectra. A forthcoming paper in will make use of the multi-band CGS data to extract stellar population constraints from the color information of the subcomponents in nearby Es.
- affiliation: School of Astronomy and Space Science, Nanjing University, Nanjing 210093, China
- affiliation: The Observatories of the Carnegie Institution for Science, 813 Santa Barbara Street, Pasadena, CA 91101, USA
- affiliation: Key Laboratory of Modern Astronomy and Astrophysics, Nanjing University, Nanjing 210093, China
- affiliation: The Observatories of the Carnegie Institution for Science, 813 Santa Barbara Street, Pasadena, CA 91101, USA
- affiliation: Giant Magellan Telescope Organization, 251 South Lake Avenue, Suite 300, Pasadena, CA 91101, USA
- affiliation: Key Laboratory for Research in Galaxies and Cosmology, Shanghai Astronomical Observatory, Chinese Academy of Sciences, 80 Nandan Road, Shanghai 200030, China
- affiliation: Department of Physics and Astronomy, 4129 Frederick Reines Hall, University of California, Irvine, CA 92697-4575, USA
- slugcomment: Submitted to The Astrophysical Journal Letters.
- Bernardi, M., Roche, N., Shankar, F., & Sheth, R. K. 2011, MNRAS, 412, L6
- Bezanson, R., van Dokkum, P. G., Tal, T., et al. 2009, ApJ, 697, 1290
- Bluck, A. F. L., Conselice, C. J., Buitrago, F., et al. 2012, ApJ, 747, 34
- Buitrago, F., Trujillo, I., Conselice, C. J., et al. 2008, ApJ, 687, L61
- Cappellari, M., di Serego Alighieri, S., Cimatti, A., et al. 2009, ApJ, 704, L34
- Cimatti, A., Nipoti, C., & Cassata, P. 2012, MNRAS, 422, L62
- Coccato, L., Gerhard, O., & Arnaboldi, M. 2010, MNRAS, 407, L26
- Conselice, C. J., Bluck, A. F. L., Buitrago, F., et al. 2011, MNRAS, 413, 80
- Daddi, E., Renzini, A., Pirzkal, N., et al. 2005, ApJ, 626, 680
- Damjanov, I., Abraham, R. G., Glazebrook, K., et al. 2011, ApJ, 739, L44
- Dekel, A., Sari, R., & Ceverino, D. 2009, ApJ, 703, 785
- Fan, L., Lapi, A., De Zotti, G., & Danese, L. 2008, ApJ, 689, L101
- Gargiulo, A., Saracco, P., Longhetti, M., La Barbera, F., & Tamburri, S. 2012, MNRAS, 425, 2698
- Greene, J. E., Murphy, J. D., Comerford, J. M., Gebhardt, K., & Adams, J. J. 2012, ApJ, 750, 32
- Grogin, N. A., Kocevski, D. D., Faber, S. M., et al. 2011, ApJS, 197, 35
- Guo, Y., McIntosh, D. H., Mo, H. J., et al. 2009, MNRAS, 398, 1129
- Ho, L. C., Li, Z.-Y., Barth, A. J., Seigar, M. S., & Peng, C. Y. 2011, ApJS, 197, 21
- Hopkins, P. F., Bundy, K., Murray, N., et al. 2009a, MNRAS, 398, 898
- Hopkins, P. F., Cox, T. J., Dutta, S. N., et al. 2009b, ApJS, 181, 135
- Hopkins, P. F., Lauer, T. R., Cox, T. J., Hernquist, L., & Kormendy, J. 2009c, ApJS, 181, 486
- Huang, S., Ho, L. C., Peng, C. Y., Li, Z.-Y., & Barth, A. J. 2013, ApJ, in press (arXiv:1212.2639)
- Johansson, P. H., Naab, T., & Ostriker, J. P. 2012, ApJ, 754, 115
- Koekemoer, A. M., Faber, S. M., Ferguson, H. C., et al. 2011, ApJS, 197, 36
- Kormendy, J., Fisher, D. B., Cornell, M. E., & Bender, R. 2009, ApJS, 182, 216
- Larson, R. B. 1975, MNRAS, 173, 671
- Lackner, C. N., Cen, R., Ostriker, J. P., & Joung, M. R. 2012, MNRAS, 3426
- Lauer, T. R., Ajhar, E. A., Byun, Y.-I., et al. 1995, AJ, 110, 2622
- Li, Z.-Y., Ho, L. C., Barth, A. J., & Peng, C. Y. 2011, ApJS, 197, 22
- Merritt, D., Mikkola, S., & Szell, A. 2007, ApJ, 671, 53
- Negroponte, J., & White, S. D. M. 1983, MNRAS, 205, 1009
- Newman, A. B., Ellis, R. S., Bundy, K., & Treu, T. 2012, ApJ, 746, 162
- Nipoti, C., Treu, T., Leauthaud, A., et al. 2012, MNRAS, 422, 1714
- Onodera, M., Renzini, A., Carollo, M., et al. 2012, ApJ, 755, 26
- Oser, L., Ostriker, J. P., Naab, T., Johansson, P. H., & Burkert, A. 2010, ApJ, 725, 2312
- Papovich, C., Bassett, R., Lotz, J. M., et al. 2012, ApJ, 750, 93
- Patel, S. G., van Dokkum, P. G., Franx, M., et al. 2013, ApJ, in press (arXiv:1208.0341)
- Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2010, AJ, 139, 2097
- Prieto, M., Eliche-Moral, M. C., Balcells, M., et al. 2013, MNRAS, 428, 999
- Scarlata, C., Carollo, C. M., Lilly, S., et al. 2007, ApJS, 172, 406
- Scoville, N., Abraham, R. G., Aussel, H., et al. 2007, ApJS, 172, 38
- Sérsic, J. L. 1968, Atlas de Galaxias Australes (Córdoba: Obs. Astron., Univ. Nac. Córdoba)
- Szomoru, D., Franx, M., & van Dokkum, P. G. 2012, ApJ, 749, 121
- Tasca, L. A. M., Kneib, J.-P., Iovino, A., et al. 2009, A&A, 503, 379
- Toomre, A., & Toomre, J. 1972, ApJ, 178, 623
- Trujillo, I., Feulner, G., Goranova, Y., et al. 2006, MNRAS, 373, L36
- van der Wel, A., Rix, H.-W., Wuyts, S., et al. 2011, ApJ, 730, 38
- van Dokkum, P. G., Whitaker, K. E., Brammer, G., et al. 2010, ApJ, 709, 1018