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.

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.

Mass assembly and morphological transformations since z ∼ 3 from CANDELS

We quantify the evolution of the stellar mass functions (SMFs) of star-forming and quiescent galaxies as a function of morphology from z ∼ 3 to the present. Our sample consists of ∼50 000 galaxies in the CANDELS fields (∼880 arcmin^2), which we divide into four main morphological types, i.e. pure bulge-dominated systems, pure spiral disc-dominated, intermediate two-component bulge+disc systems and irregular disturbed galaxies. At z ∼ 2, 80 per cent of the stellar mass density of star-forming galaxies is in irregular systems. However, by z ∼ 0.5, irregular objects only dominate at stellar masses below 10^9 M_⊙. A majority of the star-forming irregulars present at z ∼ 2 undergo a gradual transformation from disturbed to normal spiral disc morphologies by z ∼ 1 without significant interruption to their star formation. Rejuvenation after a quenching event does not seem to be common except perhaps for the most massive objects, because the fraction of bulge-dominated star-forming galaxies with M^*/M_⊙ > 10^10.7 reaches 40 per cent at z < 1. Quenching implies the presence of a bulge: the abundance of massive red discs is negligible at all redshifts over 2 dex in stellar mass. However, the dominant quenching mechanism evolves. At z > 2, the SMF of quiescent galaxies above M^* is dominated by compact spheroids. Quenching at this early epoch destroys the disc and produces a compact remnant unless the star-forming progenitors at even higher redshifts are significantly more dense. At 1 < z < 2, the majority of newly quenched galaxies are discs with a significant central bulge. This suggests that mass quenching at this epoch starts from the inner parts and preserves the disc. At z < 1, the high-mass end of the passive SMF is globally in place and the evolution mostly happens at stellar masses below 10^10 M_⊙. These low-mass galaxies are compact, bulge-dominated systems, which were environmentally quenched: destruction of the disc through ram-pressure stripping is the likely process.

The Mice at play in the CALIFA survey. A case study of a gas-rich major merger between first passage and coalescence

We present optical integral field spectroscopy (IFS) observations of the Mice, a major merger between two massive (≳10^11 M_⊙) gas-rich spirals NGC 4676A and B, observed between first passage and final coalescence. The spectra provide stellar and gas kinematics, ionised gas properties, and stellar population diagnostics, over the full optical extent of both galaxies with ~1.6 kpc spatial resolution. The Mice galaxies provide a perfect case study that highlights the importance of IFS data for improving our understanding of local galaxies. The impact of first passage on the kinematics of the stars and gas has been significant, with strong bars most likely induced in both galaxies. The barred spiral NGC 4676B exhibits a strong twist in both its stellar and ionised gas disk. The edge-on disk galaxy NGC 4676A appears to be bulge free, with a strong bar causing its “boxy” light profile. On the other hand, the impact of the merger on the stellar populations has been minimal thus far. By combining the IFS data with archival multiwavelength observations we show that star formation induced by the recent close passage has not contributed significantly to the total star formation rate or stellar mass of the galaxies. Both galaxies show bicones of high ionisation gas extending along their minor axes. In NGC 4676A the high gas velocity dispersion and Seyfert-like line ratios at large scaleheight indicate a powerful outflow. Fast shocks (vs ~ 350 km s^-1) extend to ~6.6 kpc above the disk plane. The measured ram pressure (P/k = 4.8 × 10^6 K cm^-3) and mass outflow rate (~8−20 M_⊙ yr^-1) are similar to superwinds from local ultra-luminous infrared galaxies, although NGC 4676A only has a moderate infrared luminosity of 3 × 10^10 L_⊙. Energy beyond what is provided by the mechanical energy of the starburst appears to be required to drive the outflow. Finally, we compare the observations to mock kinematic and stellar population maps extracted from a hydrodynamical merger simulation. The models show little enhancement in star formation during and following first passage, in agreement with the observations. We highlight areas where IFS data could help further constrain the models.

The star formation history of CALIFA galaxies: Radial structures

We have studied the radial structure of the stellar mass surface density (μ∗) and stellar population age as a function of the total stellar mass and morphology for a sample of 107 galaxies from the CALIFA survey. We applied the fossil record method based on spectral synthesis techniques to recover the star formation history (SFH), resolved in space and time, in spheroidal and disk dominated galaxies with masses from 10^9 to 10^12 M_⊙. We derived the half-mass radius, and we found that galaxies are on average 15% more compact in mass than in light. The ratio of half-mass radius to half-light radius (HLR) shows a dual dependence with galaxy stellar mass; it decreases with increasing mass for disk galaxies, but is almost constant in spheroidal galaxies. In terms of integrated versus spatially resolved properties, we find that the galaxy-averaged stellar population age, stellar extinction, and μ_∗ are well represented by their values at 1 HLR. Negative radial gradients of the stellar population ages are present in most of the galaxies, supporting an inside-out formation. The larger inner (≤1 HLR) age gradients occur in the most massive (10^11 M_⊙) disk galaxies that have the most prominent bulges; shallower age gradients are obtained in spheroids of similar mass. Disk and spheroidal galaxies show negative μ∗ gradients that steepen with stellar mass. In spheroidal galaxies, μ∗ saturates at a critical value (~7 × 10^2 M_⊙/pc^2 at 1 HLR) that is independent of the galaxy mass. Thus, all the massive spheroidal galaxies have similar local μ_∗ at the same distance (in HLR units) from the nucleus. The SFH of the regions beyond 1 HLR are well correlated with their local μ_∗, and follow the same relation as the galaxy-averaged age and μ_∗; this suggests that local stellar mass surface density preserves the SFH of disks. The SFH of bulges are, however, more fundamentally related to the total stellar mass, since the radial structure of the stellar age changes with galaxy mass even though all the spheroid dominated galaxies have similar radial structure in μ_∗. Thus, galaxy mass is a more fundamental property in spheroidal systems, while the local stellar mass surface density is more important in disks.

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.

Spatially resolved star formation main sequence of galaxies in the CALIFA survey

The "main sequence of galaxies"—defined in terms of the total star formation rate ψ versus the total stellar mass M *—is a well-studied tight relation that has been observed at several wavelengths and at different redshifts. All earlier studies have derived this relation from integrated properties of galaxies. We recover the same relation from an analysis of spatially resolved properties, with integral field spectroscopic (IFS) observations of 306 galaxies from the CALIFA survey. We consider the SFR surface density in units of log(M_⊙ yr^−1 Kpc^−2) and the stellar mass surface density in units of log(M_⊙ Kpc^−2) in individual spaxels that probe spatial scales of 0.5–1.5 Kpc. This local relation exhibits a high degree of correlation with small scatter (σ = 0.23 dex), irrespective of the dominant ionization source of the host galaxy or its integrated stellar mass. We highlight (i) the integrated star formation main sequence formed by galaxies whose dominant ionization process is related to star formation, for which we find a slope of 0.81 ± 0.02; (ii) for the spatially resolved relation obtained with the spaxel analysis, we find a slope of 0.72 ± 0.04; and (iii) for the integrated main sequence, we also identified a sequence formed by galaxies that are dominated by an old stellar population, which we have called the retired galaxies sequence.

Multiwavelength study of star-forming galaxies in the local Universe

We have studied in detail the properties of local active star-forming galaxies from the UCM survey, and in particular their star-formation histories. We have quantified the relative importance of the current episode of star formation in comparison to the underlying older stellar populations. We have also determined the total stellar mass function and burst mass function for the UCM sample using the M/L calculated for each galaxy. Integrating this mass function we obtained the contribution of the star-forming galaxies to the total stellar mass density of the local Universe.

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.

The merger history of massive spheroids since z ∼ 1 is size independent

Using a compilation of 379 massive (stellar mass M ≳ 10^11 M_⊙) spheroid-like galaxies from the near-infrared Palomar/DEEP-2 survey, we investigated, up to z ∼ 1, whether the presence of companions depends on the size of the host galaxy. We explored the presence of companions for mass ratios with respect to the central massive galaxy down to 1 : 10 and 1 : 100, and within projected distances of 30, 50 and 100 kpc of these objects. We found evidence that these companions are equally distributed around both compact and extended massive spheroid-like galaxies. This suggests that, at least since z ∼ 1, the merger activity in these objects is nearly homogeneous across the whole population and that the merger history is not affected by the size of the host galaxy. Our results could indicate that compact and extended massive spheroid-like galaxies are increasing in size at the same rate.

The relation between cool cluster cores and Herschel-detected star formation in brightest cluster galaxies

We present far-infrared (FIR) analysis of 68 brightest cluster galaxies (BCGs) at 0.08 < z < 1.0. Deriving total infrared luminosities directly from Spitzer and Herschel photometry spanning the peak of the dust component (24-500 μm), we calculate the obscured star formation rate (SFR). 22^+6.2 _–5.3% of the BCGs are detected in the far-infrared, with SFR = 1-150 M ☉ yr^–1. The infrared luminosity is highly correlated with cluster X-ray gas cooling times for cool-core clusters (gas cooling time <1 Gyr), strongly suggesting that the star formation in these BCGs is influenced by the cluster-scale cooling process. The occurrence of the molecular gas tracing Hα emission is also correlated with obscured star formation. For all but the most luminous BCGs (L_TIR > 2 × 10^11 L_☉), only a small (≤0.4 mag) reddening correction is required for SFR(Hα) to agree with SFR_FIR. The relatively low Hα extinction (dust obscuration), compared to values reported for the general star-forming population, lends further weight to an alternate (external) origin for the cold gas. Finally, we use a stacking analysis of non-cool-core clusters to show that the majority of the fuel for star formation in the FIR-bright BCGs is unlikely to originate from normal stellar mass loss.

The host galaxies and black holes of typical z~0.5-1.4 AGNs

We study the stellar and star formation properties of the host galaxies of 58 X-ray-selected AGNs in the GOODS portion of the Chandra Deep Field South (CDF-S) region at z ~ 0.5-1.4. The AGNs are selected such that their rest-frame UV to near-infrared spectral energy distributions (SEDs) are dominated by stellar emission; i.e., they show a prominent 1.6 μm bump, thus minimizing the AGN emission "contamination." This AGN population comprises approximately 50% of the X-ray-selected AGNs at these redshifts. We find that AGNs reside in the most massive galaxies at the redshifts probed here. Their characteristic stellar masses (M_* ~ 7.8 × 10^10 and M_* ~ 1.2 × 10^11 M_☉ at median redshifts of 0.67 and 1.07, respectively) appear to be representative of the X-ray-selected AGN population at these redshifts and are intermediate between those of local type 2 AGNs and high-redshift (z ~ 2) AGNs. The inferred black hole masses (M_BH ~ 2 × 10^8 M_☉) of typical AGNs are similar to those of optically identified quasars at similar redshifts. Since the AGNs in our sample are much less luminous (L_2–10 keV < 10^44 erg s^−1) than quasars, typical AGNs have low Eddington ratios (η ~ 0.01-0.001). This suggests that, at least at intermediate redshifts, the cosmic AGN "downsizing" is due to both a decrease in the characteristic stellar mass of typical host galaxies and less efficient accretion. Finally, there is no strong evidence in AGN host galaxies for either highly suppressed star formation (expected if AGNs played a role in quenching star formation) or elevated star formation when compared to mass-selected (i.e., IRAC-selected) galaxies of similar stellar masses and redshifts.

The star formation and extinction coevolution of UV-selected galaxies over 0.05 < z < 1.2

We use a new stacking technique to obtain mean mid-IR and far-IR to far-UV flux ratios over the rest-frame near-UV, near-IR color-magnitude diagram. We employ COMBO-17 redshifts and COMBO-17 optical, GALEX far- and near-UV, and Spitzer IRAC and MIPS mid-IR photometry. This technique permits us to probe the infrared excess (IRX), the ratio of far-IR to far-UV luminosity, and the specific star formation rate (SSFR) and their coevolution over 2 orders of magnitude of stellar mass and over redshift 0.1 < z < 1.2. We find that the SSFR and the characteristic mass (Script M_0) above which the SSFR drops increase with redshift (downsizing). At any given epoch, the IRX is an increasing function of mass up to Script M_0. Above this mass the IRX falls, suggesting gas exhaustion. In a given mass bin below Script M_0, the IRX increases with time in a fashion consistent with enrichment. We interpret these trends using a simple model with a Schmidt-Kennicutt law and extinction that tracks gas density and enrichment. We find that the average IRX and SSFR follow a galaxy age parameter ξ, which is determined mainly by the galaxy mass and time since formation. We conclude that blue-sequence galaxies have properties which show simple, systematic trends with mass and time such as the steady buildup of heavy elements in the interstellar media of evolving galaxies and the exhaustion of gas in galaxies that are evolving off the blue sequence. The IRX represents a tool for selecting galaxies at various stages of evolution.

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.