Free Access
Issue
A&A
Volume 592, August 2016
Article Number A93
Number of page(s) 9
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/201628188
Published online 04 August 2016

© ESO, 2016

1. Introduction

Chemical patterns, such as the mass-metallicity relation, the age-metallicity relation, and metallicity gradients, store relevant information on the history of galaxy assembly. Among these patterns, observational results show that slope of the gas-phase metallicity gradients correlates with stellar mass if expressed in dex kpc-1, but that this correlation is erased if the metallicity gradients are renormalized by adopting a characteristic scale length for the disc component (e.g. Garnett et al. 1997; Sánchez et al. 2014; Ho et al. 2015; Sánchez-Menguiano et al. 2016). Similar results are found for discs formed in hydrodynamical cosmological simulations (Tissera et al. 2015) and by using analytical models (Prantzos & Boissier 2000; Ho et al. 2015). Both observational and theoretical results favour the existence of metallicity gradients, the origin of which would be linked to the joint action of gas outflows and inflows that regulate star formation activity and chemical abundances (e.g. Chiappini et al. 2001; Mollá & Díaz 2005; Gibson et al. 2013; Mott et al. 2013). An important fraction of the synthesized chemical elements, however, is locked into stellar populations (SPs) as the gas is transformed into stars. Hence, the chemical abundances of the SPs provide complementary information to understand the chemical loop between the gas-phase medium and the SPs within a context of galaxy formation (e.g. Calura et al. 2012; Tissera et al. 2012).

Most observational studies are devoted to the gas-phase abundances in discs traced by HII regions because the study of the chemical properties of the underlying SPs is more difficult. Estimations for the Milky Way and some nearby galaxies suggested the existence of negative metallicity slopes (e.g. Bland-Hawthorn & Freeman 2003; Carraro et al. 2007; Hughes et al. 1994; Kudritzki et al. 2008). The metallicity properties of the SPs in galaxies outside the Local Group have been studied using colour distributions. In these works, the colour gradients are interpreted as the result of the contribution of younger and less enriched stars to the outer parts (e.g. de Jong 1996; MacArthur et al. 2004; Muñoz-Mateos et al. 2007). Two-dimensional (2D) observations of nearby galaxies with detailed information on the age and metallicities of SPs are being gathered. As a result, more robust properties of the SPs in nearby galaxies are available (e.g. Rosales-Ortega et al. 2010). In particular, Sánchez-Blázquez et al. (2014) reported normalized metallicity gradients of SPs in the disc components for a set of face-on spirals in the CALIFA survey. Similar to the results reported for the gas-phase normalized metallicity gradients (Sánchez et al. 2014), the SP normalized metallicity slopes were not found to correlate with the global properties of galaxies such as stellar mass or luminosity. A weak trend for old SPs (>6 Gyr) to show steeper metallicity gradients than younger SPs (<2 Gyr) was reported. The metallicity and age gradients of galaxies with a variety of morphology are analysed by González Delgado et al. (2015). These authors also find a dependence of both gradients on stellar mass (or the potential well). In particular they found that spiral galaxies have negative metallicity and age profiles in agreement with an inside-out formation. For disc-dominated galaxies, they identified the steeper negative metallicity slopes to take place for intermediate stellar mass galaxies (Mollá & Díaz 2005).

The metallicity gradients of the SPs are the result of the combination of stars formed at different epochs from gas in the interstellar medium (ISM) that have different physical and chemical properties. Disentangling their contributions is a complex task, albeit crucial to build up a comprehensive understanding of how star formation, chemical production, and mixing take place in galaxies. Mergers and interactions are known to be efficient mechanisms to redistribute angular momentum and, as a consequence, they can affect the metallicity distribution of gas and stars. The action of migration was also claimed to be important in shaping the metallicity gradients, and flattening the abundance profiles of the SPs, principally of the old SPs (e.g. Friedli et al. 1994; Roškar et al. 2008; Di Matteo et al. 2013). However, more realistic simulations where the vertical velocity of the stars is taken into account suggest that strong migration might be prevented (Minchev et al. 2014).

In the current cosmological paradigm, galaxy formation is a complex process that involves the action of different physical mechanisms acting with different temporal and spatial scales. In this context, hydrodynamical simulations that include chemical evolution are powerful tools to study the relation between chemical patterns and the history of galaxy assembly (e.g. Mosconi et al. 2001; Scannapieco et al. 2006; Tissera et al. 2012; Vogelsberger et al. 2014; Schaye et al. 2015). Similarly, chemodynamical models have extensively analysed the formation of metallicity gradients in disc galaxies (e.g. Chiappini et al. 1997; Mollá et al. 1997; Ho et al. 2015). Within a hierarchical clustering scenario, different mechanisms play a role in the assembly history of galaxies, however most numerical works suggest an inside-out formation for the disc components (Scannapieco et al. 2009; Brook et al. 2012; Few et al. 2012; Minchev et al. 2014). Particularly, Pilkington et al. (2012) analysed the metallicity gradients of the young SPs as tracers of the gas-phase metallicity, finding that the efficiency of star formation as a function of radius was a main factor in determining the metallicity gradients. As expected, different supernova (SN) feedback models produced different results (Gibson et al. 2013; Snaith et al. 2013). Tissera et al. (2015) analysed the metallicity gradients of the gas phase of the same set of simulated galaxies analysed in our work. They detected a trend for the gas-phase metallicity slopes to correlate with stellar mass. These authors found a relationship between the efficiency of star formation and the slope of the gas-phase metallicity profile in agreement with the observations of Stott et al. (2014). According to detailed numerical simulations of prepared merger events, these relationships could be the result of the action of violent events in the regulation of star formation and chemical mixing (Tissera et al. 2015; and Sillero, in prep.).

In this work, we focus on the study of the metallicity gradients of the SPs in the disc components of the same simulated galaxies whose gas-phase components were analysed by Tissera et al. (2015) in the redshift range [0, 2). We use the 12 +log O/H abundance ratios and explore correlations with other physical parameters such as stellar mass and half-mass radius. We also extended the analysis of the metallicity gradients for SPs with different ages. We performed this analysis with a two-fold goal: searching for clues for galaxy formation and setting constraints on our galaxy formation model. The paper is organized as follows. Section 2 explains the main characteristics of the numerical models and the galaxy catalogue. Section 3 is dedicated to the analysis of the SP metallicity gradients. Our conclusions summarize the main findings.

thumbnail Fig. 1

Surface mass density, oxygen abundance, and ages for the SPs in a simulated discs shown for illustration purposes (only stars with ϵ< 0.5 were included). Each frame is 5reff wide and the half-mass radius of this stellar disc is reff = 11 kpc.

2. Numerical experiments

We analysed the stellar disc components of simulated galaxies selected from a cosmological simulation performed with a version of the code GADGET-3, an update of GADGET-2 (Springel 2005), optimized for massive parallel simulations of highly heterogeneous systems. This version includes treatments for metal-dependent radiative cooling, stochastic star formation (SF), chemical enrichment, and the multiphase model for the ISM and the SN feedback scheme of Scannapieco et al. (2005, 2006). The SN feedback model is able to successfully trigger galactic mass-loaded winds without introducing mass-scale parameters.

Our code considered energy feedback by type II (SNII) and type Ia (SNIa) supernovae. These processes are grafted into a multiphase model for the ISM described in detail by Scannapieco et al. (2006). Briefly, the energy released by each SN event is distributed between the cold and hot phases of the simulated ISM. The adopted ISM multiphase model faciliates the coexistence of diffuse and dense gas phases. In this model, each gas particle defines the cold and hot phases using a local entropy criteria. This scheme allows particles to decouple hydrodynamically from particular low-entropy particles if they are not part of a shock front. Each cold gas particle stores the SN energy received from nearby SN events until it fullfils the conditions to join their local hot phase. This SN feedback scheme does not include parameters that depend on the global properties of the given galaxy (e.g. the total mass and size) thus making it suitable for cosmological simulations, where systems with different masses form in a complex way. The SN energy is equally distributed between the cold and hot gas phases surrounding a stellar source.

As mentioned in the Introduction, we analysed the SPs of the simulated disc components whose gas-phase medium was studied by Tissera et al. (2015). They were selected from the so-called S230D simulation, whose initial conditions are consistent with the concordance model with ΩΛ = 0.7, Ωm = 0.3, Ωb = 0.04, a normalization of the power spectrum of σ8 = 0.9 and H0 = 100 h km s-1 Mpc-1, with h = 0.7. The simulated volume represents a box of 14 Mpc comovil side, resolved with 2 × 2303 initial particles, achieving a mass resolution of 5.9 × 106h-1M and 9.1 × 105h-1M for the dark matter and initial gas particles, respectively. The gravitational softening is 0.7 kpc. We acknowledge the fact that the initial condition represents a small volume of the Universe. It was chosen to represent a typical region with no large structure nearby. De Rossi et al. (2013) verified that the growth of the haloes was reproduced well by confronting the halo mass function with those obtained from the Millennium simulation (Fakhouri et al. 2010).

Simulation S230D used a high gas density threshold for the star formation and an energy per SN event of 0.7 × 1051 erg. This combination diminishes the over transformation of gas into stars at very high redshift. Simulation S230D injects 80 per cent of the new synthesized chemical elements into the cold phase (the remaining is pumped directly into the hot phase). Our code also includes the chemical evolution model developed by Mosconi et al. (2001). Chemical elements produced by SNII and SNIa events are estimated by adopting the chemical yields of Woosley & Weaver (1995) and Iwamoto et al. (1999), respectively. We follow 12 individual chemical elements (such as 16O and 56Fe) and that SNII yields are metallicity dependent (Mosconi et al. 2001).

2.1. The simulated galaxy sample

The simulated galaxy sample was constructed by using a friends-of-friends algorithm to identify the virialized structures. Then the substructures were selected using the SUBFIND scheme (Springel et al. 2001). To identify the disc components, we used the angular momentum content of each particle defined as ϵ = Jz/Jz,max(E), where Jz is the angular momentum component in the direction of the total angular momentum and Jz,max(E) is the maximum Jz over all of the particles at a given binding energy, E. The details are explained in Tissera et al. (2012) and Pedrosa & Tissera (2015). We considered star particles with ϵ higher than 0.5 to be part of the stellar disc. The star particles that did not satisfy this requirement were taken to belong to the spheroidal component and are not analysed in this paper. We do not distinguish between thin and thick discs. They are both included in the analysed stellar disc components. To diminish resolution issues, we restricted our study to systems with the number of star particles larger than 2000 in the disc components. With this criterium, 37 galaxies were selected at z ~ 0.

For illustration purposes, in Fig. 1 we show face-on and edge-on projected distributions of the stellar mass, oxygen abundances, and stellar ages for one of the simulated disc components. The bulge was extracted by applying the procedure described above. From this figure, we can see the existence of spiral arms that are populated mainly by young stars. In this example, there are clear age and metallicity profiles, although there are also local inhomogenities, revealing a complex structure that varies from system to system.

The simulated galaxies analysed in this paper were also studied by Pedrosa & Tissera (2015). These authors focused on the angular momentum content of the disc and bulge components, finding that the simulated galaxies are able to reproduce observed trends reported by Romanowsky & Fall (2012). Hence, we expect the sizes of the simulated galaxies to be in agreement with observations. Tissera et al. (2015) studied the gas-phase metallicity gradients of the same galaxies we analysed in this paper in relation to the specific star formation of the galaxies. These authors reported that these galaxies reproduced observational trends at z ~ 0. In agreement with observations, Tissera et al. (2015) showed that simulated gas-phase metallicity slopes of these low stellar mass galaxies exhibit a larger variety. As a function of redshift, these authors found an increasing fraction of gas-phase metallicity profiles with negative slopes (<−0.1 dex kpc-1 in galaxies with M< 1010.5M). More massive galaxies showed a larger frequency of positive metallicity slopes for higher redshifts.

We use the CALIFA sample of Sánchez-Blázquez et al. (2014) to confront the simulated trends at z ~ 0 since it is a self-consistent sample and the stellar mass range covered is similar to our simulated mass range. For the purpose of comparison, we use the mass-weighted metallicity gradients. We note that the observed half-mass radii are defined from the surface luminosity distributions, while we use the half-mass radius of the stellar disc components. This should be taken into consideration when comparing the normalized gradients.

3. Metallicity and age gradients of the SPs

We estimated the radial abundance gradients and age profiles of the disc SP. We performed linear regression fits to the profiles in the range [0.5,1] reff, where reff is the radius that enclosed half of the stellar mass of the disc components. We find that this radial region provides stable and good fits for our profiles and is within the radial range chosen by Sánchez-Blázquez et al. (2014) for the analysis of the SP gradients in a set of face-on spiral galaxies of the CALIFA survey, which we use to compare with our simulated data. This is also the spatial range adopted to fit the gas-phase metallicity gradients by Tissera et al. (2015). The linear regressions are performed by applying a bootstrap technique. In Fig. 2 we show the metallicity and age profiles (upper and lower panels, respectively) for the same galaxy depicted in Fig. 1. The goodness of the fittings are shown with the shaded areas (95 per cent confidence bands).

thumbnail Fig. 2

Metallicity and age profiles for the SPs in the stellar disc component of the galaxy shown in Fig. 1 (thin grey lines). The thick red lines are the linear fits within the region 0.5 <r/reff< 1. The shaded areas represent the 95 per cent confidence bands of the fits.

3.1. Metallicity gradients

Figure 3 shows the metallicity slopes in units of dex kpc-1 (upper panel) and the renormalized slopes using reff (lower panel) as a function of the stellar mass of the galaxies. For comparison, we also include the observational estimations reported by Sánchez-Blázquez et al. (2014).

We define two subsamples by adopting the stellar mass limit, M< 1010.5M, similar to the analysis carried out by Tissera et al. (2015) for the gas-phase metallicity profiles, to quantify a possible dependence on the stellar mass. From Fig. 3 and Table 1, we can see that the non-normalized metallicity slopes show a very weak dependence on stellar mass. Galaxies in the low stellar mass subsample have a median metallicity slope of −0.04 ± 0.01dex kpc-1, while those in the high stellar mass subsample have −0.03 ± 0.02dex kpc-1. The median slopes and errors were estimated by applying a bootstrap resampling technique. In Table 1 we also include the standard dispersions. The most significant variation with stellar mass is the increase of the dispersion of the metallicity slopes measured in the low stellar mass galaxies.

thumbnail Fig. 3

Slopes of the non-normalized (upper panel) and normalized (lower panel) metallicity gradients of the SPs in the disc components as a function of total stellar mass of the simulated galaxies (solid, blue circles). The error bars correspond to a bootstrap error of the linear regression fits. Observations from the CALIFA survey by (Sánchez-Blázquez et al. 2014; open yellow squares) were included for comparison.

We performed similar estimations for the normalized SP metallicity gradients, (dexreff-1\hbox{${\rm dex~{\it r}_{eff}^{-1}}$}). As shown in the lower panel of from Fig. 3, there is no dependence on stellar mass and the dispersion increases significantly for all galaxies (Table 1). There are galaxies with normalized metallicity profiles with positive slopes at the two extreme of the stellar mass range, and similar behaviour is also present for the CALIFA slopes. In the case of those at the lower mass end, the fittings have a larger dispersion. The Spearman test of correlation yields no signal of correlation for both distributions: r = −0.11,p = 0.51 for the non-normalized relation and r = −0.19,p = 0.27 for the normalized relation. The gas-phase metallicity gradients of these simulated galaxies show clear trends with stellar mass as reported by Tissera et al. (2015). The weakness of the signal could be due to the action of dynamical processes or different star formation efficiencies. The SN feedback adopted for this simulation is able to regulate the transformation of gas into stars naturally in different potential wells (de Rossi et al. 2010). Interestingly, the simulated discs are able to reproduce a trend of having steeper negative metallicity gradients for masses around 1010M, as shown by Sánchez-Blázquez et al. (2014) and González Delgado et al. (2015). See also Mollá & Díaz (2005) for similar results by using analytical models. De Rossi et al. (2013) showed that there is a change in the effects of the SN feedback around this stellar mass from the combined effects of the SN strengths and the potential well of the systems. These changes affect the gas cooling time, star formation activity, and the amount of metals ejected in galactic winds. We speculate that these processes might also leave an imprint in the chemical gradients. A large statistical sample is needed to confirm these trends.

Prantzos & Boissier (2000, see also Ho et al. 2015 reported a correlation of the metallicity gradients of the gas phase with reff, which got erased if the gradients were normalized by a characteristic scale length. To test the existence of such a correlation for the simulated SPs, in Fig. 4 we show the relation for the simulated discs. In this case, we find a clearer signal of correlation for the non-normalized gradients: r ~ 0.38,p ~ 0.02. The correlation signal significantly decreases for the normalized gradients: r ~ 0.031,p ~ 0.86.

Since the characteristic radius reff seems to play a significant role, it is important to show that the simulated reff are in agreement with observations. In Fig. 5 we show a function of stellar mass together with recent observations by van der Wel et al. (2014) and Sánchez-Blázquez et al. (2014). The plot shows that the mean values are consistent with observations by van der Wel et al. (2014) within a standard deviation for late-type galaxies as expected since most of our galaxies have active star formation (Tissera et al. 2015). However, at the more massive end, our simulated systems show larger values than the reported averages. The simulated reff are in very good agreement with those estimated from the CALIFA survey1.

thumbnail Fig. 4

Slopes of the non-normalized (upper panels) and normalized (lower panels) metallicity gradients of the SPs in the disc components as a function of reff (solid blue circles). The error bars correspond to the bootstrap errors of the linear regression fits.

thumbnail Fig. 5

Mean reff of the stellar disc components as a function of total stellar mass of the simulated galaxies (solid blue circles) for the observational estimations by van der Wel et al. (2014) for early-type (red inverted triangles) and late-type (blue triangles) galaxies and for observations of the disc galaxies in the CALIFA survey (Sánchez-Blázquez et al. 2014; open yellow squares).

thumbnail Fig. 6

Slopes of non-normalized (upper panel) and normalized (lower panel) age gradients of the SPs as a function of the stellar mass (left panels) and reff (right panels) of the simulated galaxies. Similarly, observations from the CALIFA survey (Sánchez-Blázquez et al. 2014; open yellow squares) were included for comparison.

3.2. Age gradients

The stellar age profiles store information on the formation history of the stellar disc. In an inside-out disc formation scenario, negative age gradients are expected since the central regions formed first and have higher density gas to feed the star formation activity. Although this might be the case on global terms, mergers and interactions as well as secular evolution could contribute to flatten or to build positive age and/or metallicity profiles as shown in previous works (e.g. Rupke et al. 2010; Perez et al. 2011). Migration could also cause old stars to populate outer regions, flattening the age and metallicity distributions. Minor mergers might also add old stars to the outer parts of the discs (Abadi et al. 2003; Tissera et al. 2012).

We fit linear regressions to the age profiles following the same procedure as in the case of the metallicity profiles. In Fig. 6, we show the age gradients as a function of the stellar mass of the simulated galaxies. On average, the mean age gradients are negative as expected for inside-out histories of formation. The slopes of the age profiles show a very mild dependence on stellar mass. In Table 1 we show the median age gradients for the two stellar-mass subsamples. The normalized age slopes show similar behaviour, again with a large dispersion. The Spearman test shows a very weak signal for the non-normalized relations (r = 0.19,p = 0.26; upper panel) and no signal at all for the normalized relation (r = −0.044,p = 0.80; lower panel).

As a function of reff, the age gradients show no clear signal of correlation. This might be caused by the intrinsic dispersion in stellar age or by the low number of galaxies. In analytical models, a parallel behaviour between metallicity and age could be expected but the relations are more complex when the non-linear assembly of structure is considered. In fact, the age profiles could be more affected by a fraction of young stars as a function of radius. In our simulation, there are high metallicity stars of all ages, although they tend to be younger than low-metallicity stars on average.

We note that there are steeper negative normalized age gradients in the simulation than in the observations reported by Sánchez-Blázquez et al. (2014). These age profiles with steeper negative slopes could be the result of a difference in the estimation of the reff and/or of the action of physical processes poorly described in our model or of differences in the galaxy sample.

Table 1

Median metallicity (Z) and age(age) gradients for the normalized (dexreff-1\hbox{${\rm dex~{\it r}_{eff}^{-1}}$}) and non-normalized (dex kpc-1) stellar profiles.

thumbnail Fig. 7

Age-metallicity relation for simulated galaxies: a) negative age gradients (upper-left galaxy: age = −0.41;∇Z = −0.34; lower-left galaxy: age = −0.40;∇Z = −0.14) and b) positive age gradients (upper-right galaxy: age = 0.32;∇Z = −0.34; lower-right galaxy: age = 0.10;∇Z = −0.14). The colours denote the distance to the centre of mass in units of the reff. Each point depicts a star particle representing a single SP. The relation between old and new SPs clearly differs and is responsible for the different age profiles in the simulations.

The inside-out disc formation scenario seems to be the common formation path for our discs in agreement with observations. However, there are a few of them that show positive age profiles (less than 10 per cent). We found that these discs have higher percentage of recent newborn stars within reff, which, as expected, are responsible for the flattening or changing in the slope of the age profiles. However, we found that not all of these discs show positive metallicity slopes. In previous works, the positive metallicity slopes of the gas-phase components were shown to appear as a result of the effects of mergers or interactions that drove gas infalls and triggered star formation. These gas inflows are expected to have lower metallicity than the gas present in the central regions (e.g. Perez et al. 2006, 2011; Rupke et al. 2010; Tissera et al. 2015). If these inflows took place very recently, the gas-phase disc would be already enriched (Tissera et al. 2015) and the newborn stars would be young but not so metal poor with respect to the underlying SPs. Hence the age profiles would reflect the recent star formation activity in the central regions more efficiently than the metallicity gradients. Our simulations suggest that there might be a variety of situations, depending on the metallicity of the infalling gas (i.e. the slope of the gas-phase metallicity gradients) and the relative fraction of young-to-old stars present in the central regions.

As an illustration, in Fig. 7 we show the age-metallicity relation for two galaxy discs with positive (right panel) and negative (left panel) age gradients that have negative metallicity slopes. From this figure we can see that the galaxies with positive age slopes have a larger fraction of young stars within reff. Those with negative age slopes have young stars formed mainly in the outer parts of the discs or more evenly distributed along the discs. We expect positive age slopes if there is available gas to form stars in the central regions and this might happen if the gas is transported inwards by tidal torques produced during interactions with nearby satellites or galaxies.

Table 2

Median metallicity gradients, both normalized to the effective radius of each disc (in dexreff-1\hbox{${\rm dex~{\it r}_{eff}^{-1}}$}) and non-normalized (in dex kpc-1).

We also find that old stars are present at all radii. In particular, in Fig. 7 the presence of old stars in the outer regions can be clearly seen. This could be ascribed in part to migration processes (e.g. Spitoni et al. 2015) or bar diffusion (Brunetti et al. 2011). It should be also mentioned that the accretion of satellite galaxies could also contribute with old stars in the outer parts (Abadi et al. 2003; Tissera et al. 2012).

Hence, it might be possible to have galaxies with positive age profiles that have experienced secular evolution recently but their metallicity gradients still show negative slopes. This situation adds noise to observations and may prevent us from drawing clear conclusions in relation to the presence of bars and their impact on metallicity slopes (Sánchez-Blázquez et al. 2014).

3.3. Metallicity gradients for old and young stars

Sánchez-Blázquez et al. (2014) reported a trend, albeit weak, for the old SP to have normalized metallicity profiles with steeper negative slopes than those determined by young stars. We estimated the metallicity profiles for stars older than 6 Gyr and younger than 2 Gyr for both normalized and non-normalized profiles to evaluate this trend in our simulated galaxies. Figure 8 shows the estimated metallicity gradients. The median values have a mild trend for the old populations to have steeper metallicity gradients (−0.06 ± 0.01 dex kpc-1) compared to the young populations (−0.04 ± 0.02 dex kpc-1). A stronger trend is found for the normalized metallicity slopes although with larger bootstrap errors and dispersions (Table 2).

For both age subsamples, we estimated the slopes of the stellar metallicity profiles as a function of the galaxy stellar mass as shown in Fig. 9. The non-normalized gradients show a very weak correlation with stellar mass for the young SPs The slopes of the normalized metallicity profiles show no correlations for the old and young SPs with stellar mass (Table 2). The normalization of the metallicity profiles by the corresponding reff erases any dependence on the stellar mass that might have been imprinted on the metallicity profiles of the young SP. The Spearman tests provide no clear signal of correlations. Hence, the trends estimated above should be taken as indicative and be confirmed with a larger sample. Both the non-normalized and normalized metallicity gradients of the old SPs show a clear U-shape distribution (the young populations are noiser). From these trends it seems that the old populations are primarily responsible for this behaviour in the simulation (Figs. 3 and 4).

The gas-phase metallicity profiles show a clearer dependence on stellar mass as reported by Tissera et al. (2015), however, the SPs are not able to preserve this correlation in our simulations and, even for the younger stars, the signal is very weak. This suggests the action of mechanisms that contribute to mix up stars such as mergers/interactions and/or migration, as previously mentioned.

thumbnail Fig. 8

Slopes of the normalized metallicity profiles of the old SP (>6 Gyr) vs. those of the young SP (<2 Gyr) for the simulated disc components.

thumbnail Fig. 9

Normalized (left panels) and non-normalized (right panels) metallicity slopes as a function of stellar mass separated into two subpopulations: old (>6 Gyr; upper panels) and young <2 Gyr; lower panels).

4. Conclusions

We studied the metallicity and age gradients of the stellar discs z ~ 0 in a hierarchical clustering scenario with the aim of analysing whether they reproduce observational trends and how they depend on stellar mass. The metallicity gradients reflect the enrichment of the ISM at the time of star formation and, hence, can store information that is relevant to understand the star formation process, the enrichment cycle, and the effects of galaxy assembly.

The non-normalized metallicity gradients of the SPs show a very weak dependence on stellar mass of the galaxies, which is erased when normalized by reff. There is a large variety of normalized metallicity gradients. The results are in global agreement with those reported by the CALIFA survey (Sánchez-Blázquez et al. 2014), considering the differences in the definition of the scale lengths (Sect. 3). Our simulated discs show the steeper metallicity gradients for galaxies with stellar masses around 1010M in agreement with observational results (e.g. González Delgado et al. 2015).

The simulated reff of the stellar discs as a function of stellar mass are consistent with observational estimations by van der Wel et al. (2014) and by Sánchez-Blázquez et al. (2014) within a standard deviation. Our simulated stellar discs are formed by conserving the specific angular momentum content (Pedrosa & Tissera 2015) with mean values in agreement with observations (Romanowsky & Fall 2012). These are important achievements for our simulations considering that the free parameters were not fixed to reproduce these mean values.

We also find a correlation between the metallicity gradients and the reff that is erased if the profiles are normalized by this characteristic scale length. As pointed out by previous works (Prantzos & Boissier 2000), this suggests that the build up of the metallicity profiles is strongly linked with the formation of discs. Other processes such as galaxy mergers, secular evolution, migration, and accretion of small satellites could perturbate the discs and mix the gas and SPs, thereby producing changes in the relation. As reported in previous works, we expect mergers to be the one of the most important mechanisms responsible for driving inflows that contribute to produce positive metallicity and/or age profiles (Rupke et al. 2010; Perez et al. 2006, 2011). In fact, the metallicity gradients may change back to negative but they tend to remain with values above the original values. Also, the ejection of enriched material from the central regions by SN outflows could be re-accreted in the external regions of the discs later on, and contribute to flattening the metallicity profiles (Perez et al. 2011).

The simulated stellar discs show negative age profiles that, on average, are consistent with an inside-out history of formation. There are some discs with positive (or close to flat) age gradients that can be linked to the recent formation of stars in the central region r/reff< 1. These cases could be related to galaxies with a recent history of galaxy interactions that have driven gas inflows. The metallicity profiles are weakly modified since recent gas infalls would tend to transport already enriched material. Conversely the age profiles are strongly modified by these new stars since the dynamical range covered by age is much larger. Our results suggest that, in some cases, relating the presence of a bar responsible for gas inflow with the slope of metallicity profile might not be so direct.

The SPs show metallicity profiles with different slopes when they are separated according to age. The metallicity gradients of young (<2 Gyr) stars tend to be slightly flatter than those of the old (>6 Gyr) stars. These trends are in agreement with results from CALIFA survey and the Milky Way (Maciel et al. 2003).

Overall, our findings suggest that the underlying process that is responsible for determining the age and metallicity profiles is the inside-out formation of discs with specific angular momentum conservation. In spite of this, dynamical processes that can disturb the stellar distributions, together with those that affect the transformation of gas into stars, play an important role.


1

We also fitted exponential discs to the mass surface density and calculated the reff from the disc scale lengths. These values agreed with those estimated directly from the mass distributions.

Acknowledgments

Support for P.B.T. is provided by the Ministry of Economy, Development, and Toursim’s Millennium Science Initiative through grant IC120009, awarded to The Millennium Institute of Astrophysics, MAS. R.E.G.M. acknowledges support from Ciência sem Fronteiras (CNPq, Brazil). This work was partially supported by PICT 2011-0959 (ANCyT, Argentina) and PIP 2012-0396 (Conicet, Argentina). PBT acknowledges partial support from the Regular Grant UNAB 2014, Nucleo UNAB 2015 (DI-677-15/N) of Universidad Andres Bello and Fondecyt 1150334 (Conicyt).

References

  1. Abadi, M. G., Navarro, J. F., Steinmetz, M., & Eke, V. R. 2003, ApJ, 591, 499 [NASA ADS] [CrossRef] [Google Scholar]
  2. Bland-Hawthorn, J., & Freeman, K. C. 2003, in Star Formation Through Time, eds. E. Perez, R. M. Gonzalez Delgado, & G. Tenorio-Tagle, ASP Conf. Ser., 297, 457 [Google Scholar]
  3. Brook, C. B., Stinson, G. S., Gibson, B. K., et al. 2012, MNRAS, 426, 690 [NASA ADS] [CrossRef] [Google Scholar]
  4. Brunetti, M., Chiappini, C., & Pfenniger, D. 2011, A&A, 534, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Calura, F., Gibson, B. K., Michel-Dansac, L., et al. 2012, MNRAS, 427, 1401 [NASA ADS] [CrossRef] [Google Scholar]
  6. Carraro, G., Geisler, D., Villanova, S., Frinchaboy, P. M., & Majewski, S. R. 2007, A&A, 476, 217 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Chiappini, C., Matteucci, F., & Gratton, R. 1997, ApJ, 477, 765 [NASA ADS] [CrossRef] [Google Scholar]
  8. Chiappini, C., Matteucci, F., & Romano, D. 2001, ApJ, 554, 1044 [NASA ADS] [CrossRef] [Google Scholar]
  9. de Jong, R. S. 1996, A&S, 313, 45 [Google Scholar]
  10. de Rossi, M. E., Tissera, P. B., & Pedrosa, S. E. 2010, A&A, 519, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. De Rossi, M. E., Avila-Reese, V., Tissera, P. B., González-Samaniego, A., & Pedrosa, S. E. 2013, MNRAS, 435, 2736 [NASA ADS] [CrossRef] [Google Scholar]
  12. Di Matteo, P., Haywood, M., Combes, F., Semelin, B., & Snaith, O. N. 2013, A&A, 553, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Fakhouri, O., Ma, C.-P., & Boylan-Kolchin, M. 2010, MNRAS, 406, 2267 [NASA ADS] [CrossRef] [Google Scholar]
  14. Few, C. G., Gibson, B. K., Courty, S., et al. 2012, A&A, 547, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Friedli, D., Benz, W., & Kennicutt, R. 1994, ApJ, 430, L105 [NASA ADS] [CrossRef] [Google Scholar]
  16. Garnett, D. R., Shields, G. A., Skillman, E. D., Sagan, S. P., & Dufour, R. J. 1997, ApJ, 489, 63 [NASA ADS] [CrossRef] [Google Scholar]
  17. Gibson, B. K., Pilkington, K., Brook, C. B., Stinson, G. S., & Bailin, J. 2013, A&A, 554, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. González Delgado, R. M., García-Benito, R., Pérez, E., et al. 2015, A&A, 581, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Ho, I.-T., Kudritzki, R.-P., Kewley, L. J., et al. 2015, MNRAS, 448, 2030 [NASA ADS] [CrossRef] [Google Scholar]
  20. Hughes, S. M. G., Stetson, P. B., Turner, A., et al. 1994, ApJ, 428, 143 [NASA ADS] [CrossRef] [Google Scholar]
  21. Iwamoto, K., Brachwitz, F., Nomoto, K., et al. 1999, ApJS, 125, 439 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  22. Kudritzki, R.-P., Urbaneja, M. A., Bresolin, F., et al. 2008, ApJ, 681, 269 [NASA ADS] [CrossRef] [Google Scholar]
  23. MacArthur, L. A., Courteau, S., Bell, E., & Holtzman, J. A. 2004, ApJS, 152, 175 [NASA ADS] [CrossRef] [Google Scholar]
  24. Maciel, W. J., Costa, R. D. D., & Uchida, M. M. M. 2003, A&A, 397, 667 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Minchev, I., Chiappini, C., & Martig, M. 2014, A&A, 572, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Mollá, M., & Díaz, A. I. 2005, MNRAS, 358, 521 [NASA ADS] [CrossRef] [Google Scholar]
  27. Mollá, M., Ferrini, F., & Díaz, A. I. 1997, ApJ, 475, 519 [NASA ADS] [CrossRef] [Google Scholar]
  28. Mosconi, M. B., Tissera, P. B., Lambas, D. G., & Cora, S. A. 2001, MNRAS, 325, 34 [NASA ADS] [CrossRef] [Google Scholar]
  29. Mott, A., Spitoni, E., & Matteucci, F. 2013, MNRAS, 435, 2918 [NASA ADS] [CrossRef] [Google Scholar]
  30. Muñoz-Mateos, J. C., Gil de Paz, A., Boissier, S., et al. 2007, ApJ, 658, 1006 [NASA ADS] [CrossRef] [Google Scholar]
  31. Pedrosa, S., & Tissera, P. 2015, Boletin de la Asociaciòn Argentina de Astronomia, 57, 14 [NASA ADS] [Google Scholar]
  32. Perez, J., Michel-Dansac, L., & Tissera, P. B. 2011, MNRAS, 417, 580 [NASA ADS] [CrossRef] [Google Scholar]
  33. Perez, M. J., Tissera, P. B., Scannapieco, C., Lambas, D. G., & de Rossi, M. E. 2006, A&A, 459, 361 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Pilkington, K., Gibson, B. K., Brook, C. B., et al. 2012, MNRAS, 425, 969 [NASA ADS] [CrossRef] [Google Scholar]
  35. Prantzos, N., & Boissier, S. 2000, MNRAS, 313, 338 [NASA ADS] [CrossRef] [Google Scholar]
  36. Romanowsky, A. J., & Fall, S. M. 2012, ApJS, 203, 17 [NASA ADS] [CrossRef] [Google Scholar]
  37. Rosales-Ortega, F. F., Kennicutt, R. C., Sánchez, S. F., et al. 2010, MNRAS, 405, 735 [NASA ADS] [Google Scholar]
  38. Roškar, R., Debattista, V. P., Stinson, G. S., et al. 2008, ApJ, 675, L65 [NASA ADS] [CrossRef] [Google Scholar]
  39. Rupke, D. S. N., Kewley, L. J., & Chien, L. 2010, ApJ, 723, 1255 [NASA ADS] [CrossRef] [Google Scholar]
  40. Sánchez, S. F., Rosales-Ortega, F. F., Iglesias-Páramo, J., et al. 2014, A&A, 563, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Sánchez-Blázquez, P., Rosales-Ortega, F. F., Méndez-Abreu, J., et al. 2014, A&A, 570, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Sánchez-Menguiano, L., Sánchez, S. F., Pérez, I., et al. 2016, A&A, 587, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Scannapieco, C., Tissera, P. B., White, S. D. M., & Springel, V. 2005, MNRAS, 364, 552 [NASA ADS] [CrossRef] [Google Scholar]
  44. Scannapieco, C., Tissera, P. B., White, S. D. M., & Springel, V. 2006, MNRAS, 371, 1125 [NASA ADS] [CrossRef] [Google Scholar]
  45. Scannapieco, C., White, S. D. M., Springel, V., & Tissera, P. B. 2009, MNRAS, 396, 696 [NASA ADS] [CrossRef] [Google Scholar]
  46. Schaye, J., Crain, R. A., Bower, R. G., et al. 2015, MNRAS, 446, 521 [Google Scholar]
  47. Snaith, O. N., Tissera, P. B., Pedrosa, S., de Rossi, M. E., & Vílchez, J. M. 2013, Asociacion Argentina de Astronomia La Plata Argentina Book Series, 4, 129 [NASA ADS] [Google Scholar]
  48. Spitoni, E., Romano, D., Matteucci, F., & Ciotti, L. 2015, ApJ, 802, 129 [NASA ADS] [CrossRef] [Google Scholar]
  49. Springel, V. 2005, MNRAS, 364, 1105 [NASA ADS] [CrossRef] [Google Scholar]
  50. Springel, V., Yoshida, N., & White, S. D. M. 2001, New Astron., 6, 79 [Google Scholar]
  51. Stott, J. P., Sobral, D., Swinbank, A. M., et al. 2014, MNRAS, 443, 2695 [NASA ADS] [CrossRef] [Google Scholar]
  52. Tissera, P. B., White, S. D. M., & Scannapieco, C. 2012, MNRAS, 420, 255 [NASA ADS] [CrossRef] [Google Scholar]
  53. Tissera, P. B., Pedrosa, S. E., Sillero, E., & Vilchez, J. M. 2015, MNRAS, 456, 2982 [Google Scholar]
  54. van der Wel, A., Franx, M., van Dokkum, P. G., et al. 2014, ApJ, 788, 28 [NASA ADS] [CrossRef] [Google Scholar]
  55. Vogelsberger, M., Genel, S., Springel, V., et al. 2014, MNRAS, 444, 1518 [NASA ADS] [CrossRef] [Google Scholar]
  56. Woosley, S. E., & Weaver, T. A. 1995, ApJS, 101, 181 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Median metallicity (Z) and age(age) gradients for the normalized (dexreff-1\hbox{${\rm dex~{\it r}_{eff}^{-1}}$}) and non-normalized (dex kpc-1) stellar profiles.

Table 2

Median metallicity gradients, both normalized to the effective radius of each disc (in dexreff-1\hbox{${\rm dex~{\it r}_{eff}^{-1}}$}) and non-normalized (in dex kpc-1).

All Figures

thumbnail Fig. 1

Surface mass density, oxygen abundance, and ages for the SPs in a simulated discs shown for illustration purposes (only stars with ϵ< 0.5 were included). Each frame is 5reff wide and the half-mass radius of this stellar disc is reff = 11 kpc.

In the text
thumbnail Fig. 2

Metallicity and age profiles for the SPs in the stellar disc component of the galaxy shown in Fig. 1 (thin grey lines). The thick red lines are the linear fits within the region 0.5 <r/reff< 1. The shaded areas represent the 95 per cent confidence bands of the fits.

In the text
thumbnail Fig. 3

Slopes of the non-normalized (upper panel) and normalized (lower panel) metallicity gradients of the SPs in the disc components as a function of total stellar mass of the simulated galaxies (solid, blue circles). The error bars correspond to a bootstrap error of the linear regression fits. Observations from the CALIFA survey by (Sánchez-Blázquez et al. 2014; open yellow squares) were included for comparison.

In the text
thumbnail Fig. 4

Slopes of the non-normalized (upper panels) and normalized (lower panels) metallicity gradients of the SPs in the disc components as a function of reff (solid blue circles). The error bars correspond to the bootstrap errors of the linear regression fits.

In the text
thumbnail Fig. 5

Mean reff of the stellar disc components as a function of total stellar mass of the simulated galaxies (solid blue circles) for the observational estimations by van der Wel et al. (2014) for early-type (red inverted triangles) and late-type (blue triangles) galaxies and for observations of the disc galaxies in the CALIFA survey (Sánchez-Blázquez et al. 2014; open yellow squares).

In the text
thumbnail Fig. 6

Slopes of non-normalized (upper panel) and normalized (lower panel) age gradients of the SPs as a function of the stellar mass (left panels) and reff (right panels) of the simulated galaxies. Similarly, observations from the CALIFA survey (Sánchez-Blázquez et al. 2014; open yellow squares) were included for comparison.

In the text
thumbnail Fig. 7

Age-metallicity relation for simulated galaxies: a) negative age gradients (upper-left galaxy: age = −0.41;∇Z = −0.34; lower-left galaxy: age = −0.40;∇Z = −0.14) and b) positive age gradients (upper-right galaxy: age = 0.32;∇Z = −0.34; lower-right galaxy: age = 0.10;∇Z = −0.14). The colours denote the distance to the centre of mass in units of the reff. Each point depicts a star particle representing a single SP. The relation between old and new SPs clearly differs and is responsible for the different age profiles in the simulations.

In the text
thumbnail Fig. 8

Slopes of the normalized metallicity profiles of the old SP (>6 Gyr) vs. those of the young SP (<2 Gyr) for the simulated disc components.

In the text
thumbnail Fig. 9

Normalized (left panels) and non-normalized (right panels) metallicity slopes as a function of stellar mass separated into two subpopulations: old (>6 Gyr; upper panels) and young <2 Gyr; lower panels).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.