The dependence of oxygen and nitrogen abundances on stellar mass from the CALIFA survey

Context. The study of the integrated properties of star-forming galaxies is central to understand their formation and evolution. Some of these properties are extensive and therefore their analysis require totally covering and spatially resolved observations. Among these properties, metallicity can be defined in spiral discs by means of integral field spectroscopy (IFS) of individual H II regions. The simultaneous analysis of the abundances of primary elements, as oxygen, and secondary, as nitrogen, also provides clues about the star formation history and the processes that shape the build-up of spiral discs. Aims. Our main aim is to analyse simultaneously O/H and N/O abundance ratios in H II regions in different radial positions of the discs in a large sample of spiral galaxies to obtain the slopes and the characteristic abundance ratios that can be related to their integrated properties. Methods. We analysed the optical spectra of individual selected H II regions extracted from a sample of 350 spiral galaxies of the CALIFA survey. We calculated total O/H abundances and, for the first time, N/O ratios using the semi-empirical routine H II-C HI-MISTRY, which, according to Perez-Montero (2014, MNRAS, 441, 2663), is consistent with the direct method and reduces the uncertainty in the O/H derivation using [N II] lines owing to the dispersion in the O/H-N/O relation. Then we performed linear fittings to the abundances as a function of the de-projected galactocentric distances. Results. The analysis of the radial distribution both for O/H and N/O in the non-interacting galaxies reveals that both average slopes are negative, but a non-negligible fraction of objects have a flat or even a positive gradient (at least 10% for O/H and 4% for N/O). The slopes normalised to the effective radius appear to have a slight dependence on the total stellar mass and the morphological type, as late low-mass objects tend to have flatter slopes. No clear relation is found, however, to explain the presence of inverted gradients in this sample, and there is no dependence between the average slopes and the presence of a bar. The relation between the resulting O/H and N/O linear fittings at the effective radius is much tighter (correlation coefficient ρ_s = 0.80) than between O/H and N/O slopes (ρ_s = 0.39) or for O/H and N/O in the individual H ii regions (ρ_s = 0.37). These O/H and N/O values at the effective radius also correlate very tightly (less than 0.03 dex of dispersion) with total luminosity and stellar mass. The relation with other integrated properties, such as star formation rate, colour, or morphology, can be understood only in light of the found relation with mass

Photoionization models of the CALIFA HII regions I. Hybrid models

Photoionization models of HIT regions require as input a description of the ionizing spectral energy distribution (SED) and of the gas distribution, in terms of ionization parameter U and chemical abundances (e.g., O/H and N/O). A strong degeneracy exists between the hardness of the SED and U, which in turn leads to high uncertainties in the determination of the other parameters, including abundances. One way to resolve the degeneracy is to fix one of the parameters using additional information. For each of the 20 000 sources of the CALIFA HIT regions catalog, a grid of photoionization models is computed assuming the ionizing SED to be described by the underlying stellar population obtained from spectral synthesis modeling. The ionizing SED is then defined as the sum of various stellar bursts of different ages and metallicities. This solves the degeneracy between the shape of the ionizing SED and U. The nebular metallicity (associated with O/H) is defined using the classical strong line method O3N2 (which gives our models the status of "hybrids"). The remaining free parameters are the abundance ratio N/O and the ionization parameter U, which are determined by looking for the model fitting [N II]/H alpha and [O III]/H beta. The models are also selected to fit [O II]/H beta. This process leads to a set of 3200 models that reproduce the three observations simultaneously. We find that the regions associated with young stellar bursts (i.e., ionized by OB stars) are affected by leaking of ionizing photons, the proportion of escaping photons having a median of 80%. The set of photoionization models satisfactorily reproduces the electron temperature derived from the [O III]lambda 4363/5007 line ratio. We determine new relations between the nebular parameters, like the ionization parameter U and the [OH]/[ III] or [S II]/[S III] line ratios. A new relation between N/O and O/H is obtained, mostly compatible with previous empirical determinations (and not with previous results obtained using photoionization models). A new relation between U and O/H is also determined. All the models are publicly available on the Mexican Million Models database 3MdB

The Califa and Hipass circular velocity function for all morphological galaxy types

The velocity function (VF) is a fundamental observable statistic of the galaxy population that is similar to the luminosity function in importance, but much more difficult to measure. In this work we present the first directly measured circular VF that is representative between 60 < v_circ < 320 km s^-1 for galaxies of all morphological types at a given rotation velocity. For the low-mass galaxy population (60 < v_circ < 170 km s^-1), we use the HI Parkes All Sky Survey VF. For the massive galaxy population (170 < v_circ < 320 km s^-1), we use stellar circular velocities from the Calar Alto Legacy Integral Field Area Survey (CALIFA). In earlier work we obtained the measurements of circular velocity at the 80% light radius for 226 galaxies and demonstrated that the CALIFA sample can produce volume-corrected galaxy distribution functions. The CALIFA VF includes homogeneous velocity measurements of both late and early-type rotation-supported galaxies and has the crucial advantage of not missing gas-poor massive ellipticals that HI surveys are blind to. We show that both VFs can be combined in a seamless manner, as their ranges of validity overlap. The resulting observed VF is compared to VFs derived from cosmological simulations of the z = 0 galaxy population. We find that dark-matter-only simulations show a strong mismatch with the observed VF. Hydrodynamic simulations fare better, but still do not fully reproduce observations

The CALIFA survey across the Hubble sequence: Spatially resolved stellar population properties in galaxies

Various different physical processes contribute to the star formation and stellar mass assembly histories of galaxies. One important approach to understanding the significance of these different processes on galaxy evolution is the study of the stellar population content of today's galaxies in a spatially resolved manner. The aim of this paper is to characterize in detail the radial structure of stellar population properties of galaxies in the nearby universe, based on a uniquely large galaxy sample, considering the quality and coverage of the data. The sample under study was drawn from the CALIFA survey and contains 300 galaxies observed with integral field spectroscopy. These cover a wide range of Hubble types, from spheroids to spiral galaxies, while stellar masses range from M_* ∼ 10^9 to 7 x 10^11 M_⨀. We apply the fossil record method based on spectral synthesis techniques to recover the following physical properties for each spatial resolution element in our target galaxies: the stellar mass surface density (μ_*), stellar extinction (A_V), light-weighted and mass-weighted ages (_L, _M), and mass-weighted metallicity (_M). To study mean trends with overall galaxy properties, the individual radial profiles are stacked in seven bins of galaxy morphology (E, S0, Sa, Sb, Sbc, Sc, and Sd). We confirm that more massive galaxies are more compact, older, more metal rich, and less reddened by dust. Additionally, we find that these trends are preserved spatially with the radial distance to the nucleus. Deviations from these relations appear correlated with Hubble type: earlier types are more compact, older, and more metal rich for a given M-star, which is evidence that quenching is related to morphology, but not driven by mass. Negative gradients of _L are consistent with an inside-out growth of galaxies, with the largest _L gradients in Sb-Sbc galaxies. Further, the mean stellar ages of disks and bulges are correlated and with disks covering a wider range of ages, and late-type spirals hosting younger disks. However, age gradients are only mildly negative or flat beyond R∼2 HLR (half light radius), indicating that star formation is more uniformly distributed or that stellar migration is important at these distances. The gradients in stellar mass surface density depend mostly on stellar mass, in the sense that more massive galaxies are more centrally concentrated. Whatever sets the concentration indices of galaxies obviously depends less on quenching/morphology than on the depth of the potential well. There is a secondary correlation in the sense that at the same M_* early-type galaxies have steeper gradients. The μ_* gradients outside 1 HLR show no dependence on Hubble type. We find mildly negative _M gradients, which are shallower than predicted from models of galaxy evolution in isolation. In general, metallicity gradients depend on stellar mass, and less on morphology, hinting that metallicity is affected by both - the depth of the potential well and morphology/quenching. Thus, the largest _M gradients occur in Milky Way-like Sb-Sbc galaxies, and are similar to those measured above the Galactic disk. Sc spirals show flatter _M gradients, possibly indicating a larger contribution from secular evolution in disks. The galaxies from the sample have decreasing-outward stellar extinction; all spirals show similar radial profiles, independent from the stellar mass, but redder than E and S0. Overall, we conclude that quenching processes act in manners that are independent of mass, while metallicity and galaxy structure are influenced by mass-dependent processes

Stellar kinematics across the Hubble sequence in the CALIFA survey: general properties and aperture corrections

We present the stellar kinematic maps of a large sample of galaxies from the integral-field spectroscopic survey CALIFA. The sample comprises 300 galaxies displaying a wide range of morphologies across the Hubble sequence, from ellipticals to late-type spirals. This dataset allows us to homogeneously extract stellar kinematics up to several effective radii. In this paper, we describe the level of completeness of this subset of galaxies with respect to the full CALIFA sample, as well as the virtues and limitations of the kinematic extraction compared to other well-known integral-field surveys. In addition, we provide averaged integrated velocity dispersion radial profiles for different galaxy types, which are particularly useful to apply aperture corrections for single aperture measurements or poorly resolved stellar kinematics of high-redshift sources. The work presented in this paper sets the basis for the study of more general properties of galaxies that will be explored in subsequent papers of the survey

No direct coupling between bending of galaxy disc stellar age and light profiles

We study the stellar properties of 44 face-on spiral galaxies from the Calar Alto Legacy Integral Field Area survey via full spectrum fitting techniques. We compare the age profiles with the surface brightness distribution in order to highlight differences between profile types (type I, exponential profile; and II, down-bending profile). We observe an upturn ('U-shape') in the age profiles for 17 out of these 44 galaxies with reliable stellar information up to their outer parts. This 'U-shape' is not a unique feature for type II galaxies but can be observed in type I as well. These findings suggest that the mechanisms shaping the surface brightness and stellar population distributions are not directly coupled. This upturn in age is only observable in the light-weighted profiles while it flattens out in the mass-weighted profiles. Given recent results on the outer parts of nearby systems and the results presented in this Letter, one of the most plausible explanations for the age upturn is an early formation of the entire disc(~10 Gyr ago) followed by an inside-out quenching of the star formation

Nearby supernova host galaxies from the CALIFA survey II. Supernova environmental metallicity

The metallicity of a supernova progenitor, together with its mass, is one of the main parameters that can rule their outcome. We present the second study of nearby supernova (SN) host galaxies (0.005 10 dex) by targeted searches. We also found no evidence that the metallicity at the SN location differs from the average metallicity at the galactocentric distance of the SNe. By extending our SN sample with published metallicities at the SN location, we are able to study the metallicity distributions for all SN subtypes split into SN discovered in targeted and untargeted searches. We confirm a bias toward higher host masses and metallicities in the targeted searches. Combining data from targeted and untargeted searches we found a sequence from higher to lower local metallicity: SN Ia, Ic, and II show the highest metallicity, which is significantly higher than SN Ib, IIb, and Ic-BL. Our results support the picture of SN Ib resulting from binary progenitors and, at least part of, SN Ic being the result of single massive stars stripped of their outer layers by metallicity driven winds. We studied several proxies of the local metallicity frequently used in the literature and found that the total host metallicity allows for the estimation of the metallicity at the SN location with an accuracy better than 0.08 dex and very small bias. In addition, weak AGNs that cannot be seen in the total spectrum may weakly bias (by 0.04 dex) the metallicity estimate from the galaxy integrated spectrum

Two-dimensional multi-component photometric decomposition of CALIFA galaxies

We present a two-dimensional multi-component photometric decomposition of 404 galaxies from the Calar Alto Legacy Integral Field Area data release 3 (CALIFA-DR3). They represent all possible galaxies with no clear signs of interaction and not strongly inclined in the final CALIFA data release. Galaxies are modelled in the g, r, and i Sloan Digital Sky Survey (SDSS) images including, when appropriate, a nuclear point source, bulge, bar, and an exponential or broken disc component. We use a human-supervised approach to determine the optimal number of structures to be included in the fit. The dataset, including the photometric parameters of the CALIFA sample, is released together with statistical errors and a visual analysis of the quality of each fit. The analysis of the photometric components reveals a clear segregation of the structural composition of galaxies with stellar mass. At high masses (log (M-*/M-circle dot) > 11), the galaxy population is dominated by galaxies modelled with a single Sersic or a bulge+disc with a bulge-to-total (B/T) luminosity ratio B/T > 0.2. At intermediate masses (9.5 < log (M-*/M-circle dot) < 11), galaxies described with bulge+disc but B/T < 0.2 are preponderant, whereas, at the low mass end (log (M-*/M-circle dot) < 9.5), the prevailing population is constituted by galaxies modelled with either pure discs or nuclear point sources+discs (i.e., no discernible bulge). We obtain that 57% of the volume corrected sample of disc galaxies in the CALIFA sample host a bar. This bar fraction shows a significant drop with increasing galaxy mass in the range 9.5 < log (M-*/M-circle dot) < 11.5. The analyses of the extended multi-component radial profile result in a volume-corrected distribution of 62%, 28%, and 10% for the so-called Type I (pure exponential), Type II (down-bending), and Type III (up-bending) disc profiles, respectively. These fractions are in discordance with previous findings. We argue that the different methodologies used to detect the breaks are the main cause for these differences

CALIFA, the Calar Alto Legacy Integral Field Area Survey IV. Third public data release

This paper describes the third public data release (DR3) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. Science-grade quality data for 667 galaxies are made public, including the 200 galaxies of the second public data release (DR2). Data were obtained with the integral field spectrograph PMAS/PPak mounted on the 3.5 m telescope at the Calar Alto Observatory. Three different spectral setups are available: i) a low-resolution V500 setup covering the wavelength range 3745-7500 angstrom (4240-7140 angstrom unvignetted) with a spectral resolution of 6.0 angstrom (FWHM) for 646 galaxies, ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 angstrom (3650-4620 angstrom unvignetted) with a spectral resolution of 2.3 angstrom (FWHM) for 484 galaxies, and iii) the combination of the cubes from both setups (called COMBO) with a spectral resolution of 6.0 angstrom and a wavelength range between 3700-7500 angstrom (3700-7140 angstrom unvignetted) for 446 galaxies. The Main Sample, selected and observed according to the CALIFA survey strategy covers a redshift range between 0.005 and 0.03, spans the color-magnitude diagram and probes a wide range of stellar masses, ionization conditions, and morphological types. The Extension Sample covers several types of galaxies that are rare in the overall galaxy population and are therefore not numerous or absent in the CALIFA Main Sample. All the cubes in the data release were processed using the latest pipeline, which includes improved versions of the calibration frames and an even further improved image reconstruction quality. In total, the third data release contains 1576 datacubes, including-1.5 million independent spectra

Aperture effects on the oxygen abundance determinations from Califa data

This paper aims to provide aperture corrections for emission lines in a sample of spiral galaxies from the Calar Alto Legacy Integral Field Area Survey (CALIFA) database. In particular, we explore the behavior of the log([O III] λ5007/Hβ)/([N II] λ6583/Hα) (O3N2) and log[N II] lambda 6583/Hα (N2) flux ratios since they are closely connected to different empirical calibrations of the oxygen abundances in star-forming galaxies. We compute the median growth curves of Hα, Hα/Hβ, O3N2, and N-2 up to 2.5R(50) and 1.5 disk R-eff. These distances cover most of the optical spatial extent of the CALIFA galaxies. The growth curves simulate the effect of observing galaxies through apertures of varying radii. We split these growth curves by morphological types and stellar masses to check if there is any dependence on these properties. The median growth curve of the Hα flux shows a monotonous increase with radius with no strong dependence on galaxy inclination, morphological type, and stellar mass. The median growth curve of the Hα/HβH ratio monotonically decreases from the center toward larger radii, showing for small apertures a maximum value of ≈10% larger than the integrated one. It does not show any dependence on inclination, morphological type, and stellar mass. The median growth curve of N-2 shows a similar behavior, decreasing from the center toward larger radii. No strong dependence is seen on the inclination, morphological type, and stellar mass. Finally, the median growth curve of O3N2 increases monotonically with radius, and it does not show dependence on the inclination. However, at small radii it shows systematically higher values for galaxies of earlier morphological types and for high stellar mass galaxies. Applying our aperture corrections to a sample of galaxies from the SDSS survey at 0.02 ≤ z ≤ 0.3 shows that the average difference between fiber-based and aperture-corrected oxygen abundances, for different galaxy stellar mass and redshift ranges, reaches typically to ≈11%, depending on the abundance calibration used. This average difference is found to be systematically biased, though still within the typical uncertainties of oxygen abundances derived from empirical calibrations. Caution must be exercised when using observations of galaxies for small radii (e.g., below 0.5 R_eff) given the high dispersion shown around the median growth curves. Thus, the application of these median aperture corrections to derive abundances for individual galaxies is not recommended when their fluxes come from radii much smaller than either R_50 or R_eff