The CALIFA Survey : Exploring the Oxygen Abundance in the Local Universe

We present here a review of the latest results on the spatially-resolved analysis of the stellar populations and ionized gas of disk-dominated galaxies based on Calar Alto Legacy Integral Field Area (CALIFA) data. CALIFA is an ongoing integral field spectroscopy (IFS) survey of galaxies in the Local Universe (0.005 < z < 0.03) that has already obtained spectroscopic information up to ∼2.5 re with a spatial resolution better than ∼1 kpc for a total number of more than 600 galaxies of different morphological types, covering the color-magnitude diagram up to MR < −18 mag. With nearly 2000 spectra obtained for each galaxy, CALIFA offers one of the best IFUdatasets to study the star Galaxies 2015, 3 165 formation histories and chemical enrichment of galaxies. In this article, we focus on the main results from the analysis of the oxygen abundances based on the study of ionized gas in H II regions and individual spaxels and their relation to the global properties of galaxies, using an updated/revised dataset with more galaxies and ionized regions. In summary, we have confirmed previous published results indicating that: (1) the M-Z relation does not present a secondary relation to the star formation rate, when the abundance is measured at the effective radius; (2) the oxygen abundance presents a strong correlation with the stellar surface density (Σ-Z relation); (3) the oxygen abundance profiles present three well-defined regimes: (i) an overall negative radial gradient between 0.5 and 2 re, with a characteristic slope of αO/H ∼−0.1 dex/re; (ii) a universal flattening beyond >2 re; and (iii) an inner drop at <0.5 re that depends on mass; (4) the presence of bending in the surface brightness profile of disk galaxies is not clearly related to either the change in the shape of the oxygen abundance profile or the properties of the underlying stellar population. All of these results indicate that disk galaxies present an overall inside-out growth, with chemical enrichment and stellar mass growth tightly correlated and dominated by local processes and limited effects of radial mixing or global outflows. However, clear deviations are shown with respect to this simple scenario, which affect the abundance profiles in both the innermost and outermost regions of galaxies.


Introduction
We present here a brief summary of the status of the on-going Calar Alto Legacy Integral Field Area (CALIFA) survey with an emphasis on the results that have been recently published.In particular, we make a summary of the most relevant results found regarding the properties of H II regions discovered using this survey, focused on two main results: (i) the existence of a characteristic gradient of oxygen abundance irrespective of the properties of the galaxies; and (ii) the lack of a secondary correlation with the star formation in the mass-metallicity (oxygen abundance) relation.

CALIFA: A Brief Introduction
The Calar Alto Legacy Integral Field Area (CALIFA) survey [1] is an ongoing large project of the Centro Astronómico Hispano-Alemán at the Calar Alto observatory to obtain spatially-resolved spectra for 600 local (0.005 < z < 0.03) galaxies by means of integral field spectroscopy (IFS).CALIFA observations started in June 2010 with the Potsdam Multi-Aperture Spectrograph (PMAS), mounted to the 3.5-m telescope, utilizing the large (74 × 64 ) hexagonal field-of-view (FoV) offered by the PPak fiber bundle [2,3].PPak was created for the Disk Mass Survey [4].Each galaxy is observed using two different setups, an intermediate spectral resolution one (V 1200, R∼1650) that covers the blue range of the optical wavelength range (3700-4700 Å) and a low-resolution one (V 500, R∼850) that covers the first order of the optical wavelength range (3750-7500 Å).A diameter-selected sample of 939 galaxies was drawn from the seventh data release of the Sloan Digital Sky Survey (SDSS), which is described in [5].From this mother sample, the 600 target galaxies are randomly selected.
Combining the techniques of imaging and spectroscopy through optical IFS provides a more comprehensive view of individual galaxy properties than any traditional survey.CALIFA-like observations were collected during the feasibility studies [6,7] and the PPak IFS Nearby Galaxy Survey (PINGS, [8]), a predecessor of this survey.The first results based on those datasets already explored their information content (e.g., [8,9]).
Compared to other IFS surveys, CALIFA offers a unique combination of: (i) a sample covering a wide range of morphological types in a wide range of masses, sampling the color-magnitude diagram for M R > − 18 mag; (ii) a large FoV, that guarantees covering the entire optical extension of the galaxies up to 2.5 r e for 80% of the sample; and (iii) an accurate spatial sampling, with a typical spatial resolution of ∼1 kpc for the entire sample, which allows one to spatially resolve the spectroscopic properties of most relevant structures in galaxies (spiral arms, bars, bulges, H II regions).The penalty for a better spatial sampling of the galaxies is the somehow limited number of galaxies in the survey (e.g., with respect to MaNGA and SAMI ).In terms of spectral resolution, only the blue wavelength range is sampled with a similar spectral resolution as these other two surveys.
As a legacy survey, one of the main goals of the CALIFA collaboration is to grant public access of the fully-reduced data cubes, in a sequence of data releases (DR).In November 2012, we delivered our first DR [10], comprising 200 data cubes corresponding to 100 objects (http://califa.caha.es/DR1/).After almost two years, and a major improvement in the data reduction, we presented our second data release [11], comprising 400 data cubes corresponding to 200 objects (http://califa.caha.es/DR2/), the first of October 2014.The third and final DR will be presented in spring 2016, and it will comprise more than 600 objects (1200 data cubes).

CALIFA: Main Science Results
The data products that can be derived from the IFU datasets obtained by the CALIFA survey comprise information on the stellar populations, ionized gas, mass distribution and stellar and gas kinematics.Similar data products are derived for any of the indicated projects: Atlas3D [12], MaNGA [13] or SAMI [14].In summary, they comprise a panoramic view of the spatially-resolved spectroscopic properties of these galaxies.
Different science goals have been already addressed using this information: (i) New techniques have been developed to understand the spatially-resolved star formation histories (SFH) of galaxies [15,16].We found solid evidence that mass assembly in typical galaxies happens from inside-out [17].The SFH and chemical enrichment of bulges and early-type galaxies are fundamentally related to the total stellar mass, while for disk galaxies, it is more related to the local stellar mass density [18,19]; negative age gradients indicate that the quenching is progressing outwards in massive galaxies [18], and age and metallicity gradients suggest that galaxy bars have not significantly altered the SFH of spirals [20]; finally, we explore the spatially-resolved stellar populations across the Hubble sequence [19] and how mergers influence the assembly of blue elliptical galaxies [21].(ii) We explore the origin of the low intensity, LINER -like, ionized gas in galaxies.These regions are clearly not related to star formation activity or to AGN activity [22].They are most probably related to post-AGB ionization in many cases [23,24] and [25].(iii) We explore the aperture and resolution effects on the data.CALIFA provides a unique tool to understand the aperture and resolution effects in larger single-fiber (like SDSS) and IFS surveys (like MaNGA, SAMI).We explored the effects of the dilution of the signal in different gas and stellar population properties [26] and proposed a new empirical aperture correction for the SDSS data [27].(iv) CALIFA is the first IFU survey that allows gas and stellar kinematic studies for all morphologies with enough spectroscopic resolution to study: (1) the kinematics of the ionized gas [28]; (2) the effects of bars in the kinematics of galaxies [29,30]; (3) the effects of the interaction stage on the kinematic signatures [31]; (4) and measure the bar pattern speeds in late-type galaxies [32].(v) We extend the measurements of the angular momentum of galaxies to previously-unexplored ranges of morphology and ellipticity [33], and propose a new dynamical classification scheme for galaxies [34].(vi) We explore in detail the effects of galaxy interaction in the enhancement of the star formation rate and the ignition of galactic outflows [35,36].(vii) We study the nature of the progenitors of SNe [37].(viii) Finally, we explore the star formation indicators for extended objects and the capability of using Hα as a good SFR tracer [38].Further results based on the analysis of the H II regions will be discussed in the next section.

Results of Our Studies of the H II Regions
The project to derive the properties of the H II regions in galaxies in the Local Universe using IFS data was initiated based on the data from the PINGS survey [39].This survey acquired IFS mosaic data for a dozen medium-sized nearby galaxies.In [8,40], we studied in detail the ionized gas and H II regions of the largest galaxy in the sample (NGC 628).The main results of these studies are included in the contribution by Rosales-Ortega in this volume.We then continued the acquisition of IFS data for a larger sample of visually-classified face-on spiral galaxies [6], as part of the feasibility studies for the CALIFA survey [1].The spatially-resolved properties of a typical galaxy in this sample, UGC9837, were presented by [7].
In [41], we presented a new method to detect, segregate and extract the main spectroscopic properties of H II regions from IFS data (HIIEXPLORER).Figure 1, left-panel, illustrates how this method works, showing the Hα intensity map of NGC 6155 as observed by CALIFA, together with the segmentation map of the ionized regions detected.A preliminary catalog of ∼2600 H II regions and aggregations extracted from 38 face-on (i < 45 • ) spiral galaxies compiled from the PINGS and CALIFA feasibility studies was presented.We found a new local scaling relation between the stellar mass density and oxygen abundance, the so-called Σ-Z relation [9].
The same catalog allows us to explore the galactocentric radial gradient of the oxygen abundance [41].Figure 1, right-panel, shows a typical radial abundance gradient for the galaxies within the CALIFA sample, derived using the O3N2 calibrator by 12 + log(O/H) = 8.73 − 0.32 log(O3N2) [42].We confirmed that up to ∼2 disk effective radii (r e ), there is a negative gradient of the oxygen abundance in all of the spiral galaxies analyzed.The gradient presents a very similar slope for all of the galaxies (∼−0.12 ± 0.11 dex/r e ), when the radial distances are measured in units of the disk effective radii.Beyond ∼2 disk effective radii, our data show evidence of a flattening in the abundance for those galaxies that cover that galactocentric distance, as we will show later, consistent with several other spectroscopic explorations, based mostly on a few objects, e.g., [43].
In [44], we used the updated CALIFA catalog of ∼7000 H II regions, extracted from 306 individual galaxies, to study the radial oxygen abundance gradient up to 3-4 disk effective radii.We confirmed that the abundance gradients present a common slope up to ∼2 effective radii, with a distribution compatible with being produced by random fluctuations, for all galaxies when normalized to the disk effective radius of α O/H = −0.1 dex/r e (using the O3N2 calibrator by [42]).No significant differences are found on the basis of the morphological type, presence or absence of bars, absolute magnitude and/or stellar mass.The only clear deviation from the common slope is seen in galaxies with evidence of interaction or undergoing a merging process.For these galaxies, the gradient is significantly flatter, in agreement with [45].To further explore these issues, we applied the same procedure to the IFS data provided by the CALIFA survey.In [46], we presented the first results based on the catalog of H II regions extracted from an enlarged sample of galaxies (∼100).We studied the dependence of the M-Z relation on the star formation rate.We found no secondary relation different from the one induced by the well-known relation between star formation and mass, contrary to what was claimed by other authors [47,48] based on single-aperture spectroscopic data (SDSS).Although the reason for the discrepancy is still not clear, we postulated that simple aperture bias, like the one presented in previous datasets, may induce the reported secondary relation.Figure 2 presents an updated version of these results, including the last list of galaxies analyzed until January 2015 (271 galaxies from the CALIFA sample together with 31 galaxies from the CALIFA pilot studies).This sample of galaxies comprises all galaxies observed up to November 2014 (492 objects), for which we have at least five H II regions/aggregations within a range between 0.5 and 2.1 r e of galactocentric distances, thus those galaxies for which we have a well-defined oxygen abundance at the effective radius as defined in [46].In total, the current sample is three-times larger than the one presented in our previous studies, with a much better coverage of the mass range.
The left panel of Figure 2 shows the M-Z relation found for these galaxies, with the color code indicating the integrated SFR for each galaxy.Notably, the stronger gradient in SFR is along the stellar mass, as expected for star-forming galaxies.Once subtracting the best fitted function to the M-Z relation, the residuals of the abundance do not present any evident secondary relation to the SFR (right panel).Thus, the results presented in [46] are confirmed with a sample of galaxies enlarged by almost a factor of two.We tested the results using different calibrators for the oxygen abundance, including the O3N2 and N2 of both [42,49], the ONS by [50] and the counter-method by [51], without significant qualitative results.We also confirmed the local Σ-Z relation unveiled by [9], with a larger statistical sample of H II regions (∼20,000).This nebular Σ-Z relation is flatter than the one derived for the average stellar populations [52], but both agree for the younger stars, as expected if the most recent stars are born from the chemically-enriched ISM. Figure 3 shows the most recently-derived Σ-Z relation using the ∼9000 H II regions detected in the 492 CALIFA galaxies described before.The dispersion around the average relation is ∼0.08 dex, just slightly larger than the one reported for the global M-Z relation, in agreement with the results presented in [8,46].In addition, we derived the same relation using a spaxel-by-spaxel selection of the star-forming regions, presented in [53], that comprise a total of nearly one million regions.Figure 3, right panel, shows the derived distribution using this dataset, which presents general good agreement with the Σ-Z relation derived using the H II regions.The main difference is that the relation in the right panel seems to be slightly steeper than the relation in the left panel.This is induced by an aperture effect, since the figure on the left corresponds to a smoothing of the figure on the right for apertures of the typical size of the H II regions (3 -4 ).
Figure 3. Left panel: distribution of the oxygen abundance along the stellar mass surface density for the ∼9000 H II regions detected for the galaxies presented in Figure 2. The image and contours show the density of points along the distribution, with the first contour encircling 95% of the regions and the percentage decreasing by 20% for each successive contour.The yellow solid points indicate the average abundance for a set of successive bins of 0.2 dex in stellar mass density, with the error bars indicating the standard deviation of the abundance distribution for each bin (with a typical value of ∼0.08 dex).Right panel: this shows the same distribution of the ∼1 million spaxels extracted from the same CALIFA galaxies for which ionization conditions is dominated by star formation, following the criteria outlined in [46].Image and contours present the same density distribution.
These results agree with the main conclusion of our previous study [41], where a limited sample of 38 face-on spiral galaxies was analyzed, using similar methods.In a companion article [20], we analyze the radial gradient of stellar metallicity, where we have found consistent results.New studies seem to confirm the existence of a characteristic or benchmark slope in the abundance gradient of disk galaxies [54].
Beyond ∼2 r e , the data confirms our previous result of the presence of a flattening in the abundance in most galaxies with detected H II regions at this radial distance.A change in the slope of the radial gradient of oxygen abundance has been reported by several authors, e.g., [8,43,55], although with less significant numbers.On the other hand, in the inner regions, the abundance gradient seems to present a drop or flattening for some galaxies, although this does not seem to be a general property of the abundance gradients.
All of these results were recently confirmed by a study performed both using the H II regions and the ionized gas spaxel-by-spaxel Sánchez-Menguiano et al. [53].In that study a detailed analysis was performed of the overall shape of the abundance gradient, restricting our sample to the face-on spiral galaxies within the CALIFA survey (122 objects).The total sample of individual abundance values comprises 185,000 spaxels and 7100 H II regions/aggregations, distributed across the optical extension of these galaxies.It was found that the abundance profiles present three well-defined ranges with different behaviors: (1) within ∼0.5-2 r e , the abundance profile presents a negative gradient with a common slope ranging between −0.07 dex/r e and −0.14 dex/r e depending on the calibrator adopted; (2) beyond >2 r e , all galaxies seem to present a flattening in the abundance gradient, irrespective of their properties (morphological type, mass, luminosity); and (3) within <0.5 r e , some galaxies present a drop/flattening in the abundance distribution, deviating from the overall negative gradient, which seems to depend on the stellar mass of the galaxy (being more intense at higher masses).We should note here that the mixing of ionization between star formation and diffuse emission does not significantly affect the result due to the selection criteria, which impose a minimum EW (Hα) of 6Å for the selected spaxels (and H II aggregations), well above the typical EW(Hα) < 3Å of the diffuse ionization.That selection criteria implies that for the same continuum level, the diffuse ionization contributes always less than 50% to the total flux and typically less than 10%.We repeated the analysis restricting the selection to those spaxels with an EW(Hα) > 20Å without any significant variation in the results.
Figure 4 illustrates these results, showing the common abundance gradient (left panel), very similar to the one reported in [44], and the different shapes in the abundance profile at different masses (right panel), where the three regimes described before are evident.For this particular case, the O3N2 calibrator proposed by [49] was adopted (12 + log(O/H) = 8.50 − 0.22 log(O3N2)).However, we tested the qualitative results using not only this calibrator, but the O3N2 presented by [42], the ONS by [50] and the PYQZ [56].
We would like to enter into the discussion of the accuracy of the different calibrators; however, in all of our studies, we prefer those calibrators anchored to oxygen abundances derived using the direct method.Those calibrators per se avoid the intrinsic problem in any photo-ionization model-based calibrator, that is the double dependence of each line ratio with both the abundance and the ionization strength (plus many other physical parameters, like the relative abundance of nitrogen to oxygen, the corresponding sulfur to oxygen ratio, dust attenuation, geometry of the nebula, shape of the spectrum of the ionization source, and so on).The abundance calibrators based on the direct method provide a direct relation between the line ratios and the abundances, at the expense of a larger nominal dispersion in those relations and once assumed that the sample of H II regions used for the calibration cover a similar range of ionization conditions as the nebulae to be measured.Thus, the main advantage of this method is that the possible degeneracies between the abundance and the ionization parameter are intrinsically included, since both parameters present a well-known correlation, e.g., [57], a fact that has been well known since the early studies of the ionization conditions in H II regions [58].Only the abundance derivations based on photoionization models should either impose somehow that correlation, e.g., [59], to derive reliable results or restrict the ionization source to reproduce that correlation (e.g., PYQZ), a problem that is avoided using calibrators based on the direct method, like O3N2 or ONS.
The second criticism is that the calibrators using the direct method hardly derive oxygen abundances above solar, in contrast with the metallicity-derived for the stellar populations using different methods, e.g., [19,52], not only in the Local Universe, but at larger redshifts.In the naive perception that chemical enrichment follows an instant recycling and a closed-box model, gas and stellar metallicity should be tightly linked, and they can only increase with time; therefore, the oxygen abundances derived using the direct method should be biases, since they cannot reproduce the supra-solar metallicities observed in some stellar populations.This criticism, while following the common sense, is basically wrong.It is true that the first generations of stars enrich the ISM fast and that gas-phase metallicity grows quickly.However, most of the chemical enrichment models predict a decrease of the gas-phase metallicity as metal-poor gas (or even pristine gas) is accreted by the galaxy, in particular in the last 8 Gy, when the galaxy evolution has been dominated by quiescent processes.Since stars are metal traps, they exhibit an abundance of gas when they were formed, and since most of the stellar continuum in the optical range is dominated by main-sequence stars, they are basically showing the average metallicity of the gas ∼4 Gy ago, not the current one that could well be lower than that of the average stellar population.[46]: i.e., the abundance gradient for each galaxy was subtracted by the predicted value corresponding to its mass as derived using the M-Z relation, and then, the average abundance for all the galaxies was added, to scale the gradient to an average value.The outermost contour encircles 95% of the total number of spaxels, decreasing 20% in each consecutive contour.Finally, for most of our studies, we adopt the O3N2 calibrator instead of more refined ones, like the ONS and the counter methods described before, for the following reasons: (i) it uses two pairs of line ratios that are very near in wavelength, and therefore, it is basically independent of the dust attenuation; (ii) it does not use just one pair of line ratios, like N2, which could be more sensitive to the relative abundance between nitrogen and oxygen; (iii) it does not use [OII], which is an emission line at the edge of the spectral regime covered by CALIFA (and most of the spectroscopic surveys at low-z), which normally present the largest problems in the spectrophotometric calibration; in addition, it is a non-resolved doublet (to our resolution), and therefore, the derivation of its flux is in many cases less accurate than the other two; (iv) it does not use the [SII] doublet that suffers from the same problems of [OII], being at the edge of the covered spectral range and being a doublet (this time resolved), but in addition, it can be affected by the telluric absorptions (at ∼6860Å) and the water vapor vibrational lines of the atmosphere.In summary, although the O3N2 calibrator many present a nominal dispersion of ∼0.08 dex, being less accurate than the nominal one presented by other calibrators, since it uses a lower number of emission lines, stronger ones and those less affected by spectrophotometric inaccuracies and dust attenuation issues, the final accuracy of the derived abundance used to be much better than the one provided by more complex calibrators.
The shape of the profile of the abundance gradient is directly related to the chemical enrichment processes and the secular evolution in galaxies.In a recent study [60], we investigate the relation between the oxygen abundance gradients in the outer disk and the presence of breaks in the surface brightness profiles and broad-band colors, for a large sample comprising 350 CALIFA galaxies.Galaxies were classified on the basis of the presence of a bending in their surface brightness profile in the three usual types (Type I, with no bending, Type II with a down-bending and Type III with an up-bending).For Type II, the presence of a U-shape in the color was confirmed, in agreement with the location of the bending.However, this reddening seems to be correlated with the flattening in the abundance gradient only for the low mass galaxies.Mass seems to affect more the shape of the abundance profile than the presence of a bending (in agreement with the results shown in Figure 4).Only for Type III galaxies, there seems to be a correlation between the change in the color and the abundance gradient, which could be explained as a downsizing effect on the population of Type III galaxies that enhances the inside-out growth in these galaxies.On the other hand, recent results by [61], based on the analysis of the star formation history of spiral galaxies with different surface brightness profile types, extracted from the CALIFA sample, found that the U-shape in the age distribution is found irrespective of the profile type.This seems to be the first clear evidence of a decoupling of the mechanisms that generate the bending in the profile and the evolution of the stellar populations.

Evidence for an Inside-Out Growth
Along this article, we describe the different results we have obtained in our studies of H II regions using IFS.Our main results can be summarized as follows: • The updated analysis presented here does not reveal any evidence of a secondary relation between the M-Z relation with the star formation rate, different from the well-known relation between the SFR and the stellar mass for star-forming galaxies.This result is based on a sample three-times larger than the one analyzed by [46] and with a much wider and better coverage in masses.• The gas abundance is strongly correlated with the stellar surface mass density (Σ-Z relation), following a correlation similar to the M-Z relation, that expands over three orders of magnitude in stellar surface density.Therefore, the stellar mass growth and the oxygen abundance enrichment are tightly related.Since the first grows inside-out, the second follows a similar trend.• The gas abundance presents a radial profile with three well-defined regimes: (1) An overall negative radial gradient, with more metal-rich H II regions in the inner areas and more metal-poor ones in the outer ones between 0.5 and 2 r e .Independently of the galaxy properties, all galaxies present a common slope for this gradient (when the galactocentric distance is normalized to the effective radius of the disk) with a value of α O/H = −0.07dex/r e and −0.14 dex/r e (depending on the oxygen abundance calibrator adopted).( 2) Beyond >2 r e , all galaxies seem to present a flattening in the abundance gradient irrespective of their properties (morphological type, mass, luminosity).(3) Within <0.5 r e , some galaxies present a drop/flattening in the abundance distribution, deviating from the overall negative gradient, which seems to depend on the stellar mass of the galaxy (being more intense at higher masses).• The presence of a bending in the surface brightness profile of disk galaxies is not clearly related to either the change in the shape of the oxygen abundance profile or the properties of the underlying stellar population.Only for Type III disk galaxies is it possible to describe a trend between some of these properties.
The existence of a universal radial decrease in the oxygen abundance has been already reported in many previous studies e.g., [8,41,43,55,[62][63][64].This observational property is compatible with our current understanding of the formation and evolution of spiral galaxies, e.g., [65] and the references therein.Gas accretion brings gas into the galactic center, where it first reaches the required density to ignite star formation.Thus, the inner regions are populated by older stars; they have undergone a faster gas reprocessing, and galaxies experience an inside-out mass growth [66,67].Both the extinction-corrected color gradients in nearby galaxies [68] and the weak dependence of the mass-size relation with redshift [69][70][71] support an inside-out scenario for the evolution of disks.Recent results, based on the analysis of the star formation history of CALIFA data, found undisputed evidence of the inside-out growth of the stellar mass in galaxies [17], at least for galaxies more massive than 10 10 M .These results are also supported by the radial distribution of the stellar ages found for the same dataset [18].
However, the described characteristic slope for the abundance gradient, independent of many of the properties of the galaxies, was only recently proposed [41].This result imposes a more severe restriction on our current understanding of how disk galaxies grow.In essence, it agrees with the recently-proposed Σ-Z relation [9] that links the gas abundance with the mass density of the underlying stellar population.
It describes how the stellar mass and the gas abundance, both fundamental products of the star formation history, grow consistently in disk galaxies, from the center to the outer parts.Together with the M-Z relation, they indicate that more massive galaxies (that trace the strongest potential well) form before and faster, accumulating more stellar mass and more metals.The presence of a common gradient in the abundance indicates that all disk-dominated galaxies of the same disk effective radius (hence, same disk mass) build up their metal content at a given normalized radius with similar efficiency.
The common slope suggests that the chemical evolution of galaxies is very similar in all disk galaxies, being compatible with a closed-box model.The classical closed-box model considers that each radial bin of a galaxy comprises primordial gas from which stars are born, live all their life-time and die in situ, according to some SFR and IMF prescriptions [72].Therefore, if the amount of primordial gas is proportional to the depth of the potential well and the efficiency of the SFR is the same for all galaxies, both the stellar mass and the enrichment would be just proportional to the time, for a given halo mass.
The lack of dependence of the M-Z and Σ-Z relations on the SFR (or the SFR surface density) also supports this scenario.Under this assumption, all galaxies should have a universal gradient of their oxygen abundance with its zero-point being proportional to the total mass.However, it is well known that the closed-box model cannot predict the right fraction of metal-poor stars with respect to the observed metallicity distribution of nearby long-lived stars in the Milky Way, e.g., [73].A more realistic model overcomes this problem by allowing the disk of galaxies to form via continuous accretion gas, driven by the gravitational force.This accretion can be compensated or even overcome for certain galaxies and certain periods by supernova explosions, e.g., [74].However, the outflow of gas is not expected to feature in the history of most spiral galaxies and is usually neglected in the models [73].This modified model is consistent with the described common radial gradient if the local gas recycling is faster than other timescales involved [75] and if the radial inflow is similar for those radial bins with the same stellar mass.
Despite this evidence, other authors consider that galactic winds play a fundamental role in the shaping of spiral galaxies, e.g., [76], and may be an essential ingredient in their formation.Galactic winds are normally associated with strong star formation events, but there is also evidence for such winds at more moderate SFR.They seem to be a ubiquitous event in high redshift galaxies, e.g., [77,78], and some theoretical work emphasizes that they are a fundamental ingredient in the evolution of galaxies, e.g., [79,80].In particular, they may play a role in the enrichment of the intergalactic medium, although in most cases, they do not seem to escape the dark matter halo around the galaxy.However, our results indicate that gas flows seem to have a less relevant effect than anticipated in the shape of the chemical gradient in galaxies, i.e., their net effect is not evident.
Our results do not rule out the existence of radial gas motions and metal mixing; they, however, add new constraints on their net effect on the chemical abundance patterns of disk galaxies.Contrary to isolated disk galaxies, interacting/merging galaxies show a clearly flatter metallicity distribution, in agreement with [81,82].This indicates that galaxy interactions can result in an effective mixing of metals.The absence of a clear difference in the radial abundance slopes of barred and unbarred galaxies suggests that bars: (i) do not significantly enhance the efficiency of metal mixing in galaxies; or (ii) produce a change roughly proportional in the gas abundance and stellar mass distribution, which compensate each other when normalized to the disk effective radius; or (iii) are of a temporary nature, and their lifetime is shorter than the timescale for chemical abundance mixing in galaxies.
In general, the properties of the ionized gas in late-type galaxies are consistent with a quiescent evolution, where gas recycling is faster than other time scales involved [75].This would imply that galaxies seem to behave locally in a similar manner rather than globally, dominated by a radial mass distribution following the potential well of the matter, with an inside-out growth that is regulated by gas inflow and local downsizing star formation.Therefore, the dominant parameter that defines the amount of metals is the stellar mass, since both parameters are the consequence of an (almost) closed-box star formation process.
The deviations from the simple monotonic decrease in the radial profile of the abundance gradient and the tight correlation between the stellar mass density and the oxygen abundance may indicate that the pure inside-out growth and pseudo-closed-box model are indeed valid only at first order.Processes that transfer gas across the galaxy, that allow the inflow of more metal-poor gas or gas recycling within the halo and quenching and/or a slow decrease of the star formation activity in the inner regions of the more massive galaxies should be fundamental ingredients in any realistic scenario that tries to describe the observed properties summarized in the current report.

Figure 1 .
Figure 1.Left panel: integral field spectroscopy (IFS)-based Hα maps, in units of 10 −16 erg s −1 cm −2 arcsec −2 , derived for NGC 6155 (color image), together with the detected H II regions shown as black segmented contours; right panel: radial distribution of the oxygen abundance derived for the individual H II regions with abundance errors below 0.15 dex, as a function of the deprojected distance (i.e., corrected for inclination), normalized to the effective radius, for the same galaxy.The size of the circles is proportional to the Hα intensity.

Figure 2 .
Figure2.Left panel: distribution of the oxygen abundances at the effective radii as a function of the integrated stellar masses for the Calar Alto Legacy Integral Field Area (CALIFA) galaxies (271, circles), together with those from the CALIFA feasibility studies(31, squares).Right panel: distribution of the differential oxygen abundances with respect to the solid-line shown in the left panel (i.e., the dependence on the stellar mass), as a function of the integrated SFRfor the CALIFA galaxies.

Figure 4 .
Figure 4. Left panel: radial density distribution of the oxygen abundance after re-scaling the oxygen abundances of each galaxy following the M-Z relation derived in[46]: i.e., the abundance gradient for each galaxy was subtracted by the predicted value corresponding to its mass as derived using the M-Z relation, and then, the average abundance for all the galaxies was added, to scale the gradient to an average value.The outermost contour encircles 95% of the total number of spaxels, decreasing 20% in each consecutive contour.The yellow diamonds represent the mean oxygen abundance values, with the error bars indicating the corresponding standard deviations, for bins of 0.25 r e .The solid-black line represents the error-weighted linear fit derived for those mean values within the range between 0.5 and 2.0 r e and the dashed-white line represents the linear relation corresponding to the characteristic values of the zero-points and slopes for the complete sample.Right panel: mean oxygen abundance radial profiles derived considering galaxies in four different bins according to their integrated stellar mass.The limits of the bins were chosen to ensure a similar number of elements in each bin: log(M/M ) < 10.1, blue diamonds; 10.1 < log(M/M ) < 10.35, red squares; 10.35 < log(M/M ) < 10.65, yellow dots; log(M/M ) > 10.65, purple triangles.Dashed vertical lines delimit the three different regimes in the oxygen abundance profiles.
Figure 4. Left panel: radial density distribution of the oxygen abundance after re-scaling the oxygen abundances of each galaxy following the M-Z relation derived in[46]: i.e., the abundance gradient for each galaxy was subtracted by the predicted value corresponding to its mass as derived using the M-Z relation, and then, the average abundance for all the galaxies was added, to scale the gradient to an average value.The outermost contour encircles 95% of the total number of spaxels, decreasing 20% in each consecutive contour.The yellow diamonds represent the mean oxygen abundance values, with the error bars indicating the corresponding standard deviations, for bins of 0.25 r e .The solid-black line represents the error-weighted linear fit derived for those mean values within the range between 0.5 and 2.0 r e and the dashed-white line represents the linear relation corresponding to the characteristic values of the zero-points and slopes for the complete sample.Right panel: mean oxygen abundance radial profiles derived considering galaxies in four different bins according to their integrated stellar mass.The limits of the bins were chosen to ensure a similar number of elements in each bin: log(M/M ) < 10.1, blue diamonds; 10.1 < log(M/M ) < 10.35, red squares; 10.35 < log(M/M ) < 10.65, yellow dots; log(M/M ) > 10.65, purple triangles.Dashed vertical lines delimit the three different regimes in the oxygen abundance profiles.