Mass-metallicity relation explored with CALIFA I. Is there a dependence on the star-formation rate?

We studied the global and local ℳ-Z relation based on the first data available from the CALIFA survey (150 galaxies). This survey provides integral field spectroscopy of the complete optical extent of each galaxy (up to 2−3 effective radii), with a resolution high enough to separate individual H II regions and/or aggregations. About 3000 individual H II regions have been detected. The spectra cover the wavelength range between [OII]3727 and [SII]6731, with a sufficient signal-to-noise ratio to derive the oxygen abundance and star-formation rate associated with each region. In addition, we computed the integrated and spatially resolved stellar masses (and surface densities) based on SDSS photometric data. We explore the relations between the stellar mass, oxygen abundance and star-formation rate using this dataset. We derive a tight relation between the integrated stellar mass and the gas-phase abundance, with a dispersion lower than the one already reported in the literature (σ_Δlog (O/H) = 0.07 dex). Indeed, this dispersion is only slightly higher than the typical error derived for our oxygen abundances. However, we found no secondary relation with the star-formation rate other than the one induced by the primary relation of this quantity with the stellar mass. The analysis for our sample of ~3000 individual H II   regions confirms (i) a local mass-metallicity relation and (ii) the lack of a secondary relation with the star-formation rate. The same analysis was performed with similar results for the specific star-formation rate. Our results agree with the scenario in which gas recycling in galaxies, both locally and globally, is much faster than other typical timescales, such like that of gas accretion by inflow and/or metal loss due to outflows. In essence, late-type/disk-dominated galaxies seem to be in a quasi-steady situation, with a behavior similar to the one expected from an instantaneous recycling/closed-box model.

Pathways to quiescence: SHARDS view on the star formation histories of massive quiescent galaxies at 1.0 < z < 1.5

We present star formation histories (SFHs) for a sample of 104 massive (stellar mass M > 10^10 M_⊙) quiescent galaxies (MQGs) at z = 1.0–1.5 from the analysis of spectrophotometric data from the Survey for High-z Absorption Red and Dead Sources (SHARDS) and HST/WFC3 G102 and G141 surveys of the GOODS-North field, jointly with broad-band observations from ultraviolet (UV) to far-infrared (far-IR). The sample is constructed on the basis of rest-frame UVJ colours and specific star formation rates (sSFRs = SFR/Mass). The spectral energy distributions (SEDs) of each galaxy are compared to models assuming a delayed exponentially declining SFH. A Monte Carlo algorithm characterizes the degeneracies, which we are able to break taking advantage of the SHARDS data resolution, by measuring indices such as MgUV and D4000. The population of MQGs shows a duality in their properties. The sample is dominated (85 per cent) by galaxies with young mass-weighted ages, t_M t_M < 2 Gyr, short star formation time-scales, 〈τ〉 ∼ 60–200 Myr, and masses log(M/M_⊙) ∼ 10.5. There is an older population (15 per cent) with t_M t_M = 2–4 Gyr, longer star formation time-scales, 〈τ〉∼ 400 Myr, and larger masses, log(M/M_⊙) ∼ 10.7. The SFHs of our MQGs are consistent with the slope and the location of the main sequence of star-forming galaxies at z > 1.0, when our galaxies were 0.5–1.0 Gyr old. According to these SFHs, all the MQGs experienced a luminous infrared galaxy phase that lasts for ∼500 Myr, and half of them an ultraluminous infrared galaxy phase for ∼100 Myr. We find that the MQG population is almost assembled at z ∼ 1, and continues evolving passively with few additions to the population.

Episodic star formation in a group of LAEs at z=5.07

We are undertaking a search for high-redshift low-luminosity Lyman Alpha sources in the SHARDS (Survey for High-z Absorption Red and Dead Sources) survey. Among the pre-selected Lyman Alpha sources two candidates were spotted, located 3.19 arcsec apart, and tentatively at the same redshift. Here, we report on the spectroscopic confirmation with Gran Telescopio Canarias of the Lyman Alpha emission from this pair of galaxies at a confirmed spectroscopic redshifts of z=5.07. Furthermore, one of the sources is interacting/merging with another close companion that looks distorted. Based on the analysis of the spectroscopy and additional photometric data, we infer that most of the stellar mass of these objects was assembled in a burst of star formation 100 Myr ago. A more recent burst (2 Myr old) is necessary to account for the measured Lyman Alpha flux. We claim that these two galaxies are good examples of Lyman Alpha sources undergoing episodic star formation. Besides, these sources very likely constitute a group of interacting Lyman Alpha emitters (LAEs).

Large-scale clustering measurements with photometric redshifts: comparing the dark matter haloes of X-ray AGN, star-forming and passive galaxies at z ≈ 1

We combine multi-wavelength data in the AEGIS-XD and C-COSMOS surveys to measure the typical dark matter halo mass of X-ray selected active galactic nuclei (AGN) [L_X(2–10 keV) > 10^42 erg s^− 1] in comparison with far-infrared selected star-forming galaxies detected in the Herschel/PEP survey (PACS Evolutionary Probe; L_IR > 10^11 L_⊙) and quiescent systems at z ≈ 1. We develop a novel method to measure the clustering of extragalactic populations that uses photometric redshift probability distribution functions in addition to any spectroscopy. This is advantageous in that all sources in the sample are used in the clustering analysis, not just the subset with secure spectroscopy. The method works best for large samples. The loss of accuracy because of the lack of spectroscopy is balanced by increasing the number of sources used to measure the clustering. We find that X-ray AGN, far-infrared selected star-forming galaxies and passive systems in the redshift interval 0.6 < z < 1.4 are found in haloes of similar mass, log M_DMH/(M_⊙ h^−1) ≈ 13.0. We argue that this is because the galaxies in all three samples (AGN, star-forming, passive) have similar stellar mass distributions, approximated by the J-band luminosity. Therefore, all galaxies that can potentially host X-ray AGN, because they have stellar masses in the appropriate range, live in dark matter haloes of log M_DMH/(M_⊙ h^−1) ≈ 13.0 independent of their star formation rates. This suggests that the stellar mass of X-ray AGN hosts is driving the observed clustering properties of this population. We also speculate that trends between AGN properties (e.g. luminosity, level of obscuration) and large-scale environment may be related to differences in the stellar mass of the host galaxies.

Testing diagnostics of nuclear activity and star formation in galaxies at z > 1

We present some of the first science data with the new Keck/MOSFIRE instrument to test the effectiveness of different AGN/SF diagnostics at z ~ 1.5. MOSFIRE spectra were obtained in three H-band multi-slit masks in the GOODS-S field, resulting in 2 hr exposures of 36 emission-line galaxies. We compare X-ray data with the traditional emission-line ratio diagnostics and the alternative mass-excitation and color-excitation diagrams, combining new MOSFIRE infrared data with previous HST/WFC3 infrared spectra (from the 3D-HST survey) and multiwavelength photometry. We demonstrate that a high [O III]/Hβ ratio is insufficient as an active galactic nucleus (AGN) indicator at z > 1. For the four X-ray-detected galaxies, the classic diagnostics ([O III]/Hβ versus [N II]/Hα and [S II]/Hα) remain consistent with X-ray AGN/SF classification. The X-ray data also suggest that "composite" galaxies (with intermediate AGN/SF classification) host bona fide AGNs. Nearly ~2/3 of the z ~ 1.5 emission-line galaxies have nuclear activity detected by either X-rays or the classic diagnostics. Compared to the X-ray and line ratio classifications, the mass-excitation method remains effective at z > 1, but we show that the color-excitation method requires a new calibration to successfully identify AGNs at these redshifts.

Measuring star formation in high-z massive galaxies: a mid-infrared to submillimetre study of the GOODS NICMOS Survey sample

We present measurements of the mean mid-infrared to submillimetre flux densities of massive (M_*≳ 10^11 M_⊙) galaxies at redshifts 1.7 < z < 2.9, obtained by stacking positions of known objects taken from the GOODS NICMOS Survey (GNS) catalogue on maps at 24 μm (Spitzer/MIPS); 70, 100 and 160 μm (Herschel/PACS); 250, 350 and 500 μm (BLAST); and 870 μm (LABOCA). A modified blackbody spectrum fit to the stacked flux densities indicates a median [interquartile] star formation rate (SFR) of SFR = 63[48, 81] M_⊙ yr^−1. We note that not properly accounting for correlations between bands when fitting stacked data can significantly bias the result. The galaxies are divided into two groups, disc-like and spheroid-like, according to their Sérsic indices, n. We find evidence that most of the star formation is occurring in n≤ 2 (disc-like) galaxies, with median [interquartile] SFR = 122[100, 150] M_⊙ yr^−1, while there are indications that the n > 2 (spheroid-like) population may be forming stars at a median [interquartile] SFR = 14[9, 20] M_⊙ yr^−1, if at all. Finally, we show that star formation is a plausible mechanism for size evolution in this population as a whole, but find only marginal evidence that it is what drives the expansion of the spheroid-like galaxies.

Star formation in a stellar mass-selected sample of galaxies to z=3 from the GOODS-NICMOS Survey

We present a study of the star-forming properties of a stellar mass-selected sample of galaxies in the GOODS (Great Observatories Origins Deep Survey) NICMOS Survey (GNS), based on deep Hubble Space Telescope (HST) imaging of the GOODS North and South fields. Using a stellar mass-selected sample, combined with HST/ACS and Spitzer data to measure both ultraviolet (UV) and infrared-derived star formation rates (SFRs), we investigate the star forming properties of a complete sample of ∼1300 galaxies down to log M_*= 9.5 at redshifts 1.5 < z < 3. Eight per cent of the sample is made up of massive galaxies with M_*≥ 10^11 M_⊙. We derive optical colours, dust extinctions and UV and infrared SFR to determine how the SFR changes as a function of both stellar mass and time. Our results show that SFR increases at higher stellar mass such that massive galaxies nearly double their stellar mass from star formation alone over the redshift range studied, but the average value of SFR for a given stellar mass remains constant over this ∼2 Gyr period. Furthermore, we find no strong evolution in the SFR for our sample as a function of mass over our redshift range of interest; in particular we do not find a decline in the SFR among massive galaxies, as is seen at z < 1. The most massive galaxies in our sample (log M_*≥ 11) have high average SFRs with values SFR_UV, corr= 103 ± 75 M_⊙ yr^−1, and yet exhibit red rest-frame (U−B) colours at all redshifts. We conclude that the majority of these red high-redshift massive galaxies are red due to dust extinction. We find that A_2800 increases with stellar mass, and show that between 45 and 85 per cent of massive galaxies harbour dusty star formation. These results show that even just a few Gyr after the first galaxies appear, there are strong relations between the global physical properties of galaxies, driven by stellar mass or another underlying feature of galaxies strongly related to the stellar mass.

Deep Herschel view of obscured star formation in the Bullet cluster

We use deep, five band (100–500 μm) data from the Herschel Lensing Survey (HLS) to fully constrain the obscured star formation rate, SFR_FIR, of galaxies in the Bullet cluster (z = 0.296), and a smaller background system (z = 0.35) in the same field. Herschel detects 23 Bullet cluster members with a total SFR_FIR = 144±14 M_⨀ yr^-1. On average, the background system contains brighter far-infrared (FIR) galaxies, with ~50% higher SFR_FIR (21 galaxies; 207± 9 M_⨀ yr^-1). SFRs extrapolated from 24 μm flux via recent templates (SFR_24 µm) agree well with SFRFIR for ~60% of the cluster galaxies. In the remaining ~40%, SFR_24 µm underestimates SFR_FIR due to a significant excess in observed S_100/S_24 (rest frame S_75/S_18) compared to templates of the same FIR luminosity.

MID-IR luminosities and UV/optical star formation rates at z < 1.4

Ultraviolet (UV) nonionizing continuum and mid-infrared (IR) emission constitute the basis of two widely used star formation (SF) indicators at intermediate and high redshifts. We study 2430 galaxies with z < 1.4 in the Extended Groth Strip with deep MIPS 24 μm observations from FIDEL, spectroscopy from DEEP2, and UV, optical, and near-IR photometry from the AEGIS. The data are coupled with dust-reddened stellar population models and Bayesian spectral energy distribution (SED) fitting to estimate dust-corrected star formation rates (SFRs). In order to probe the dust heating from stellar populations of various ages, the derived SFRs were averaged over various timescales—from 100 Myr for "current" SFR (corresponding to young stars) to 1-3 Gyr for long-timescale SFRs (corresponding to the light-weighted age of the dominant stellar populations). These SED-based UV/optical SFRs are compared to total IR luminosities extrapolated from 24 μm observations, corresponding to 10-18 μm rest frame. The total IR luminosities are in the range of normal star-forming galaxies and luminous IR galaxies (10^10-10^12 L_☉). We show that the IR luminosity can be estimated from the UV and optical photometry to within a factor of 2, implying that most z < 1.4 galaxies are not optically thick. We find that for the blue, actively star-forming galaxies the correlation between the IR luminosity and the UV/optical SFR shows a decrease in scatter when going from shorter to longer SFR-averaging timescales. We interpret this as the greater role of intermediate age stellar populations in heating the dust than what is typically assumed. Equivalently, we observe that the IR luminosity is better correlated with dust-corrected optical luminosity than with dust-corrected UV light. We find that this holds over the entire redshift range. Many so-called green valley galaxies are simply dust-obscured actively star-forming galaxies. However, there exist 24 μm detected galaxies, some with L_IR>10^11 L_☉, yet with little current SF. For them a reasonable amount of dust absorption of stellar light (but presumably higher than in nearby early-type galaxies) is sufficient to produce the observed levels of IR, which includes a large contribution from intermediate and old stellar populations. In our sample, which contains very few ultraluminous IR galaxies, optical and X-ray active galactic nuclei do not contribute on average more than ~50% to the mid-IR luminosity, and we see no evidence for a large population of "IR excess" galaxies.

Determining star formation rates for infrared galaxies

We show that measures of star formation rates (SFRs) for infrared galaxies using either single-band 24 μm or extinction-corrected Paα luminosities are consistent in the total infrared luminosity = L(TIR) ~ 10^10 L_☉ range. MIPS 24 μm photometry can yield SFRs accurately from this luminosity upward: SFR(M_☉ yr^–1) = 7.8 × 10^–10 L(24 μm, L_☉) from L(TIR) = 5× 10^9 L_☉ to 10^11 L_☉ and SFR = 7.8 × 10^–10 L(24 μm, L_☉)(7.76 × 10^–11 L(24))^0.048 for higher L(TIR). For galaxies with L(TIR) ≥ 10^10 L_☉, these new expressions should provide SFRs to within 0.2 dex. For L(TIR) ≥ 10^11 L_☉, we find that the SFR of infrared galaxies is significantly underestimated using extinction-corrected Paα (and presumably using any other optical or near-infrared recombination lines). As a part of this work, we constructed spectral energy distribution templates for eleven luminous and ultraluminous purely star forming infrared galaxies and over the spectral range 0.4 μm to 30 cm. We use these templates and the SINGS data to construct average templates from 5 μm to 30 cm for infrared galaxies with L(TIR) = 5× 10^9 to 10^13 L_☉. All of these templates are made available online.

Near-infrared and star-forming properties of local luminous infrared galaxies

We use Hubble Space Telescope (HST) NICMOS continuum and Paα observations to study the near-infrared and star formation properties of a representative sample of 30 local (d ~ 35-75 Mpc) luminous infrared galaxies (LIRGs, infrared [8-1000 μm] luminosities of log L_IR = 11-11.9 L_☉). The data provide spatial resolutions of 25-50 pc and cover the central ~3.3-7.1 kpc regions of these galaxies. About half of the LIRGs show compact (~1-2 kpc) Paα emission with a high surface brightness in the form of nuclear emission, rings, and minispirals. The rest of the sample show Paα emission along the disk and the spiral arms extending over scales of 3-7 kpc and larger. About half of the sample contains H II regions with Hα luminosities significantly higher than those observed in normal galaxies. There is a linear empirical relationship between the mid-IR 24 μm and hydrogen recombination (extinction-corrected Paα) luminosity for these LIRGs, and the H II regions in the central part of M51. This relation holds over more than four decades in luminosity, suggesting that the mid-IR emission is a good tracer of the star formation rate (SFR). Analogous to the widely used relation between the SFR and total IR luminosity of R. Kennicutt, we derive an empirical calibration of the SFR in terms of the monochromatic 24 μm luminosity that can be used for luminous, dusty galaxies.

Are long gamma-ray bursts biased tracers of star formation? Clues from the host galaxies of the Swift/BAT6 complete sample of bright LGRBs II. Star formation rates and metallicities at z < 1

Aims. Long gamma-ray bursts (LGRBs) are associated with the deaths of massive stars and might therefore be a potentially powerful tool for tracing cosmic star formation. However, especially at low redshifts (z< 1.5) LGRBs seem to prefer particular types of environment. Our aim is to study the host galaxies of a complete sample of bright LGRBs to investigate the effect of the environment on GRB formation. Methods. We studied host galaxy spectra of the Swift/BAT6 complete sample of 14 z< 1 bright LGRBs. We used the detected nebular emission lines to measure the dust extinction, star formation rate (SFR), and nebular metallicity (Z) of the hosts and supplemented the data set with previously measured stellar masses M_*. The distributions of the obtained properties and their interrelations (e.g. mass-metallicity and SFR-M_* relations) are compared to samples of field star-forming galaxies. Results. We find that LGRB hosts at z< 1 have on average lower SFRs than if they were direct star formation tracers. By directly comparing metallicity distributions of LGRB hosts and star-forming galaxies, we find a good match between the two populations up to 12 +log (O/H)~8.4−8.5, after which the paucity of metal-rich LGRB hosts becomes apparent. The LGRB host galaxies of our complete sample are consistent with the mass-metallicity relation at similar mean redshift and stellar masses. The cutoff against high metallicities (and high masses) can explain the low SFR values of LGRB hosts. We find a hint of an increased incidence of starburst galaxies in the Swift/BAT6 z< 1 sample with respect to that of a field star-forming population. Given that the SFRs are low on average, the latter is ascribed to low stellar masses. Nevertheless, the limits on the completeness and metallicity availability of current surveys, coupled with the limited number of LGRB host galaxies, prevents us from investigating more quantitatively whether the starburst incidence is such as expected after taking into account the high-metallicity aversion of LGRB host galaxies.

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.