EDP Sciences
Free access
Volume 579, July 2015
Article Number L2
Number of page(s) 6
Section Letters
DOI http://dx.doi.org/10.1051/0004-6361/201526354
Published online 19 June 2015

© ESO, 2015

1. Introduction

A kinematic classification of early-type galaxies (ETGs) might more closely reflect their true physical nature than photometry-based morphological classifications, which can be biased by inclination effects. In this sense, the ratio between ordered and random motion, V/σ, has long been used as a proxy for the rotational support of a galaxy: disc-like, rotation-dominated galaxies are intuitively associated with the highest V/σ, while the lowest values are expected in spheroidal-like, dispersion-dominated galaxies. With the advent of integral field spectroscopy, the SAURON and ATLAS3D teams developed an improved parameter, λRe, which divides galaxies into two groups: fast and slow rotators (FRs and SRs, respectively; Emsellem et al. 2007, 2011; Cappellari et al. 2011a). Based on this kinematic definition, most spirals and lenticulars (S0s) are classified as FRs, while ellipticals are found in both groups.

In the quest to understand the origin of FRs and SRs, a number of studies have focused on simulations. Jesseit et al. (2009) showed that λRe is a good indicator of the true stellar angular momentum in ETGs, and studied their possible merger origin, including the effect of gas. They already stated that the mass ratio of the encounter is crucial to determining the outcome as a SR or FR. This is in agreement with studies pointing to equal-mass mergers as the most likely origin of elliptical galaxies (Naab & Burkert 2003; Springel et al. 2005; Hoffman et al. 2010; Duc et al. 2011; Tsatsi et al. 2015). The results from Bois et al. (2011) also support this idea, and underline the role of spin (within 1:1 and 1:2 mergers, SRs are associated with retrograde and FRs with prograde encounters). These authors study the incidence of kinematic misalignments and kinematically decoupled cores, which are found primarily in SRs; however, they cannot explain the round SRs observed in ATLAS3D. A similar study with an exhaustive exploration of the parameter space was undertaken by Moody et al. (2014), who conclude that round SRs can emerge from the accumulation of many minor mergers, while FRs may be formed through a variety of pathways. The idea that ETGs might have been sculpted through multiple mergers is not new (Barnes 1985; Weil & Hernquist 1996; Bekki 2001; Bournaud et al. 2007), and in fact different studies have analysed the effects of minor and intermediate mergers onto early-type progenitors (see e.g. Bournaud et al. 2004; Di Matteo et al. 2009; Eliche-Moral et al. 2011; Hilz et al. 2013; Mapelli et al. 2015; Zinchenko et al. 2015). Using mock kinematic maps extracted from cosmological simulations, Naab et al. (2014) have recently confirmed that, even if the formation histories can be complex, the main results in a cosmological context are in full agreement with the conclusions drawn from idealised mergers.

As a consequence of the kinematic classification of ETGs, it has become apparent that most S0s are FRs and have kinematic properties that make them comparable to spirals (Cappellari et al. 2011b; Krajnović et al. 2011). Together with the fact that S0s span a whole range of bulge-to-disc ratios, this led Laurikainen et al. (2010), Cappellari et al. (2011b), and Kormendy & Bender (2012) to propose a sequence of S0s parallel to that of spirals (i.e. S0a – S0b – S0c). This recovers the original idea of Spitzer & Baade (1951) and van den Bergh (1976), bringing back the question of whether the classification parallelism reflects an underlying physical connection: are S0s faded spirals?

Observations show that gas stripping, e.g. from ram pressure in clusters, can effectively transform spirals into S0s (Crowl et al. 2005; Aragón-Salamanca et al. 2006; Crowl & Kenney 2008; Maltby et al. 2015). This constitutes an example of fading (any processes resulting in the suppression of star formation), and it can contribute to explaining the observational morphology-density relation in high-density regimes (Dressler 1980; Dressler et al. 1997; but see Cappellari et al. 2011).

However, we know that S0s do not preferentially inhabit the densest cluster environments; they are equally common in groups (Wilman et al. 2009; Bekki & Couch 2011), where mergers and tidal effects dominate (Mazzei et al. 2014a,b). Additionally, traces of past mergers have been observed in some S0s (e.g. Falcón-Barroso et al. 2004). Therefore, it seems compelling to study under which conditions S0s could emerge out of mergers, and their impact on the observed kinematics.

The CALIFA team has recently raised further doubts about S0s as faded spirals through a new diagnostic diagram: the λRe-concentration plane (van de Ven et al., in prep.). The population of late-type spirals (Sb, Sc, Sd) shows a clear incompatibility with S0s when both angular momentum (λRe) and concentration (R90/R50) are simultaneously taken into account. Provided that simple fading is not expected to significantly change the angular momentum of the galaxy, this contradicts the idea that most S0s are faded spirals. Here, we study whether major mergers can explain those differences using N-body dissipative simulations from the GalMer database. We will analyse the encounters that end up in relaxed S0-like remnants, showing that the induced changes in ellipticity, stellar angular momentum, and concentration are in agreement with the CALIFA observations.

2. Binary merger models

GalMer1 is a public database of binary N-body merger simulations that sample a wide range of mass ratios, morphological types, and orbital characteristics. The progenitor galaxies are modelled using a spherical non-rotating dark-matter halo, an optional disc and bulge, with a spatial resolution of 0.28 kpc and ~ 105 particles per galaxy. The simulations use a TreeSPH code (Semelin & Combes 2002), and take into account the effects of gas and star formation, with total simulation times in the range 34 Gyr, which typically implies ~1 Gyr of relaxation after full merger (see Chilingarian et al. 2010, for more details).

2.1. Lenticular remnants

We consider the major mergers (mass ratios 1:1 to 1:3) involving all possible combinations of two spiral progenitors (Sa, Sb, or Sd) that give rise to realistic, dynamically-relaxed S0-like remnants, based on quantitative criteria that impose structural, kinematic, SFR, and gas-content parameters typical of lenticular galaxies (Eliche-Moral et al., in prep.). In addition, we have performed a visual morphological classification to identify which remnants would have been classified as S0-like by observers (i.e. as disc galaxies without noticeable spiral arms). To this end, we have simulated photometric images of the resulting galaxies in several broad bands (B, V, R, I, and K), mimicking typical conditions of current observational surveys. We use a mass-to-light ratio which considers the stellar mass, age, and metallicity of each simulation particle according to Bruzual & Charlot (2003), with a Chabrier IMF and the Padova evolutionary tracks; for details, see Borlaff et al. (2014) and Querejeta et al. (2015). This provides us with a final sample of 67 S0-like remnants, which will be the focus of this Letter (see Table 1).

2.2. Stellar angular momentum, ellipticity, and concentration

We assume that the simulated merger remnants are observed at the median distance of the 300 CALIFA galaxies in Falcón-Barroso et al. (in prep.): Dmedian = 67 Mpc. We also consider the same spatial resolution and field of view as the survey, PSFFWHM = 1.6′′, Rmax = 35′′ (Sánchez et al. 2012).

We calculate λe according to Emsellem et al. (2011), (1)where Fi, Ri, Vi, and σi are the flux (using an appropriate mass-to-light ratio), radius, velocity, and velocity dispersion measured within each spatial bin out to the effective radius of the galaxy (RRe). The ellipticity at the effective radius εe is obtained from interpolation of the curve of growth of the luminosity-weighted ellipticities within increasingly larger isophotes.

The Petrosian concentration parameter R90/R50 It is defined as the ratio of the radii enclosing 90% and 50% of the Petrosian flux, measured on the 1D azimuthally averaged radial surface brightness profile in the simulated SDSS r band. According to Gadotti (2009), there is an equivalence between R90/R50 and bulge-to-total ratio: R90/R50 = 1.93 + 2.02(B + bar) /T, which we add on top of our plot for reference. This relation is valid for (B + bar) /T < 0.6, but then flattens for higher values; thus, the upper axis in the right panel of Fig. 1 must be interpreted with caution for (B + bar) /T > 0.6.

thumbnail Fig. 1

Stellar angular momentum (λRe) plotted against the ellipticity (εe, left) and concentration (R90/R50, right) for the GalMer simulations (original models and merger remnants), in comparison to CALIFA galaxies. All the parameters correspond to an edge-on view. The dotted line in the left plot represents the division between fast and slow rotators, and the top axis in the right plot is from Gadotti (2009), valid only up to (B + bar) /T < 0.6. See legend at the right for the symbols and colour-coding.

Open with DEXTER

3. Results

3.1. Angular momentum versus ellipticity

We compare the distribution of the remnants in the λReεe plane with the real galaxies from the CALIFA survey. The spiral progenitor models lie in the region of FRs, as expected.

Naturally, for a meaningful comparison with observational data, one has to consider projection effects: the maximum λRe is attained for an edge-on view, and it is closer to the real angular momentum (Jesseit et al. 2009); for this reason, we compare our results with the CALIFA measurements corrected for inclination. The deprojection of λRe to edge-on follows van de Ven et al. (in prep.), assuming δ = 0.5 for the global anisotropy (the mean value of the distribution proposed by Cappellari et al. 2007), but we check that varying it by ± 1σ has a <3% effect on the resulting deprojected values. When deprojecting εe, due to observational uncertainties in the inclination, the result is unfeasible for a few cases, in which we assume the maximum εintr = 1.

The left panel of Fig. 1 shows that major mergers can produce both fast- and slow-rotating S0s (9 out the 67 S0-like remnants are SRs). While only one SR S0s is present in the first CALIFA kinematic sample, in ATLAS3D, 13 out of the 36 SRs were S0-like (Cappellari et al. 2011a). We also find that, in agreement with previous studies, mergers tend to reduce both the stellar angular momentum λRe and ellipticity εe by an amount that varies largely depending on the mass ratio, gas fractions, and orbital characteristics (Jesseit et al. 2009; Bois et al. 2011; Tapia et al. 2014; Naab et al. 2014). In this diagram, real spirals (especially Sb and Sc) have a very limited overlap with S0s. Remarkably, the S0-like remnants that we find in these simulations are in excellent agreement with the real S0s, and far from the progenitor spirals that they originate from. This is a first hint that a good fraction of S0s could have originated through a mechanism involving mergers.

3.2. Angular momentum versus concentration

As mentioned in Sect. 1, the differences between spirals and S0s become more striking in the λRe-concentration plane introduced by van de Ven et al. (in prep.). The right panel of Fig. 1 compares the distribution of our remnant S0s with that of real galaxies from CALIFA in λRe versus R90/R50, showing that S0s resulting from major mergers are consistent with real lenticulars. In this plane, Sb and Sc galaxies clearly cluster towards the top-left corner; Sa galaxies show a larger scatter, extending down diagonally towards the bottom-right, where S0 and E galaxies are located. The number of Sd galaxies in this first CALIFA kinematic sample is small, but they also show little overlap with S0s.

The novel and main result presented in this Letter is that major mergers can transform spiral progenitors on the region of high λRe and low concentration into realistic S0 systems of lower λRe and higher concentration, in agreement with the CALIFA observations. The bulk shift between spirals and lenticulars is significant both observationally and according to these simulations. This confirms that major mergers are a plausible mechanism to explain this observed change, which is incompatible with simple fading.

It might seem paradoxical that the progenitor Sa galaxy has a lower R90/R50 concentration than the Sb galaxy, but this should not be seen as unrealistic: in the CALIFA sample, many Sa indeed show lower R90/R50 than a great deal of Sb galaxies. Even if B/T and concentration are often invoked as a proxy for the morphological type, there is not a one-to-one relation between them (see e.g. Graham & Worley 2008).

Within the S0 merger remnants there are some outliers, especially in terms of concentration, R90/R50, because for these remnants the total simulated time after coalescence is short (typically 0.5 Gyr after full merger). Since we consider a mass-to-light ratio (M/L) that varies according to the age and metallicity of the stellar particles in the simulation, when star formation takes place in the centre (which is very frequent in mergers, see Barnes & Hernquist 1991), we end up having an unusually high light concentration due to young stars, which reflects on the R90/R50 parameter (see also Naab & Trujillo 2006; Hopkins et al. 2009a, on the evolution of surface density profiles for elliptical galaxies resulting from simulated mergers). In fact, the most extreme concentration outliers (R90/R50> 4.5) correspond to the encounters involving the highest initial gas fractions (20%), and the highest fractions of merger-triggered star formation (new stars 15% of total stellar mass). The variation of M/L in simulations is often overlooked (assuming a constant value), and this is another important result: if we want to measure concentrations for wet merger remnants, the effects of young stellar populations on the light distribution cannot be ignored, as already pointed out by Jesseit et al. (2009). After an extra ~1 Gyr of relaxation, the central young populations will decrease their r-band luminosity by at least a factor of 2 (see Querejeta et al. 2015), and this would probably shift the points to the left and fully overlap the observational values. Additionally, it must be emphasised that most remnants have reasonable concentrations even after the short relaxation periods considered (~90% with R90/R50< 4), fully compatible with the observed values for S0s.

The dependence on initial conditions is complex, but we have found a clear trend for the highest reduction in angular momentum to be associated with retrograde encounters, in agreement with Bois et al. (2011). In fact, the emergence of SRs is almost only associated with retrograde mergers (8 out of 9, see Table 1). We have observed some degree of dependence with the pericentre distance, with a preference for short pericentres in the formation of SRs, but the dependence is milder than the correlation with the spin-orbit coupling. Within the range of major mergers considered, the mass ratio 1:1 to 1:3 does not seem to affect the outcome significantly. Interestingly, the change in λRe does not show any systematic dependence on the gas content or measured SFRs. In terms of concentration, however, there is some relation with the gas content and newborn stars, as commented above: the highest concentrations are preferentially associated with the most gas-rich encounters.

4. Discussion

As we have seen, the systematic offset between spirals and S0s in the λRe-concentration plane reported by the CALIFA team can be attributed to the effect of major mergers. A number of studies have already made it clear that disc-like remnants can emerge from mergers, even those with equal mass-ratios (Barnes 2002; Springel & Hernquist 2005; Robertson et al. 2004, 2006; Hopkins et al. 2009b; Moster et al. 2011; but see the doubts raised by Bournaud et al. 2011 about the need to include supersonic turbulence in such simulations of gas-rich mergers). In this sense, in Querejeta et al. (2015) we proved that not only can discs survive major mergers, but also that they can produce remnants with disc-bulge coupling in perfect agreement with observations of S0s (e.g. Laurikainen et al. 2010). Moreover, we have found that the discs of our S0-like remnants tend to exhibit antitruncations in agreement with those in real S0s (even reproducing their tight scaling relations, see Borlaff et al. 2014). If we combine these results with the angular momentum and concentration changes presented in this Letter, major mergers stand out as a plausible transformation mechanism to explain the origin of many S0s.

The role of mergers with other mass ratios (intermediate and minor mergers) and the combined effect of multiple mergers remains to be quantified. Based on the dry minor mergers considered by Tapia et al. (2014), and the studies by Bois et al. (2011), Moody et al. (2014) and Naab et al. (2014), we would expect intermediate and minor mergers to induce more modest changes on λRe and concentrations, but a series of mergers (which is very probable in a cosmological context) could add up and produce a net effect similar to the one that we have presented here. A major merger origin of a significant fraction of S0s would also be in agreement with the fact that higher gas fractions are expected in the early Universe, making it more likely to obtain disc-like remnants out of major mergers (keeping in mind that the gas fractions of our progenitors are representative of present-day spirals).

The fact that our results are compatible with a merger origin of S0s should not, of course, be overinterpreted. We are claiming it as a plausible mechanism, but not necessarily the only one, to explain the observed offset between most spirals and S0s in the λRe-concentration plane. In fact, some spirals and S0s clearly overlap (especially Sa galaxies); in these cases, simple fading cannot be ruled out. Remarkably, from their analysis of pseudobulges, Vaghmare et al. (2015) have recently pointed out that gas stripping alone is a viable process to transform early-type spirals into S0s, but not for late-type spirals, which agrees with the overlap found by CALIFA. Thus, it is not reasonable to claim that all lenticular galaxies are the by-products of mergers, but their relative relevance is probably higher in the transformation starting from late-type spirals. Additionally, the role of merging in “pre-processing” galaxies in filaments before falling into a cluster and in “post-processing” them during their infall should also be considered (see Vijayaraghavan & Ricker 2013; Head et al. 2014).


GalMer project: http://galmer.obspm.fr


The authors would like to acknowledge the referee for a helpful report, and the GalMer team for creating such a powerful tool, with especial thanks to Paola Di Matteo for her kind support with the database. M.Q., G.v.d.V. and J.F.B. acknowledge financial support to the DAGAL network from the People Programme (Marie Curie Actions) of the European Union’s Seventh Framework Programme FP7/2007-2013/ under REA grant agreement number PITN-GA-2011-289313. Supported by the Spanish Ministry of Economy and Competitiveness (MINECO) under projects AYA2012-31277 and AYA2013-48226-C3-1-P. J.M.A. acknowledges support from the European Research Council Starting Grant (SEDmorph; P.I. V. Wild).


Online material

Table 1

Properties of the simulated progenitors and remnants.

All Tables

Table 1

Properties of the simulated progenitors and remnants.

All Figures

thumbnail Fig. 1

Stellar angular momentum (λRe) plotted against the ellipticity (εe, left) and concentration (R90/R50, right) for the GalMer simulations (original models and merger remnants), in comparison to CALIFA galaxies. All the parameters correspond to an edge-on view. The dotted line in the left plot represents the division between fast and slow rotators, and the top axis in the right plot is from Gadotti (2009), valid only up to (B + bar) /T < 0.6. See legend at the right for the symbols and colour-coding.

Open with DEXTER
In the text