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.

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.

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.

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.

Star formation efficiency in galaxy clusters

Context. The luminous material in clusters of galaxies exists in two forms: the visible galaxies and the X-ray emitting intra-cluster medium. The hot intra-cluster gas is the major observed baryonic component of clusters, about six times more massive than the stellar component. The mass contained within visible galaxies is approximately 3% of the dynamical mass. Aims. Our aim was to analyze both baryonic components, combining X-ray and optical data of a sample of five galaxy clusters (Abell 496, 1689, 2050, 2631 and 2667), within the redshift range 0.03 < z < 0.3. We determined the contribution of stars in galaxies and the intra-cluster medium to the total baryon budget. Methods. We used public XMM-Newton data to determine the gas mass and to obtain the X-ray substructures. Using the optical counterparts from SDSS or CFHT we determined the stellar contribution. Results. We examine the relative contribution of galaxies, intra-cluster light and intra-cluster medium to baryon budget in clusters through the stellar-to-gas mass ratio, estimated with recent data. We find that the stellar-to-gas mass ratio within r(500) (the radius within which the mean cluster density exceeds the critical density by a factor of 500), is anti-correlated with the ICM temperature, which range from 24% to 6% while the temperature ranges from 4.0 to 8.3 keV. This indicates that less massive cold clusters are more prolific star forming environments than massive hot clusters.

Recovery of the star formation history of the LMC from the VISTA survey of the Magellanic system

The VISTA near infrared survey of the Magellanic System (VMC) will provide deep YJK(s) photometry reaching stars in the oldest turn-off point throughout the Magellanic Clouds (MCs). As part of the preparation for the survey, we aim to access the accuracy in the star formation history (SFH) that can be expected from VMC data, in particular for the Large Magellanic Cloud (LMC). To this aim, we first simulate VMC images containing not only the LMC stellar populations but also the foreground Milky Way (MW) stars and background galaxies. The simulations cover the whole range of density of LMC field stars. We then perform aperture photometry over these simulated images, access the expected levels of photometric errors and incompleteness, and apply the classical technique of SFH-recovery based on the reconstruction of colour-magnitude diagrams (CMD) via the minimisation of a chi-squared-like statistics. We verify that the foreground MW stars are accurately recovered by the minimisation algorithms, whereas the background galaxies can be largely eliminated from the CMD analysis due to their particular colours and morphologies. We then evaluate the expected errors in the recovered star formation rate as a function of stellar age, SFR(t), starting from models with a known age-metallicity relation (AMR). It turns out that, for a given sky area, the random errors for ages older than similar to 0.4 Gyr seem to be independent of the crowding. This can be explained by a counterbalancing effect between the loss of stars from a decrease in the completeness and the gain of stars from an increase in the stellar density. For a spatial resolution of similar to 0.1 deg(2), the random errors in SFR(t) will be below 20% for this wide range of ages. On the other hand, due to the lower stellar statistics for stars younger than similar to 0.4 Gyr, the outer LMC regions will require larger areas to achieve the same level of accuracy in the SFR( t). If we consider the AMR as unknown, the SFH-recovery algorithm is able to accurately recover the input AMR, at the price of an increase of random errors in the SFR(t) by a factor of about 2.5. Experiments of SFH-recovery performed for varying distance modulus and reddening indicate that these parameters can be determined with (relative) accuracies of Delta(m-M)(0) similar to 0.02 mag and Delta E(B-V) similar to 0.01 mag, for each individual field over the LMC. The propagation of these errors in the SFR(t) implies systematic errors below 30%. This level of accuracy in the SFR(t) can reveal significant imprints in the dynamical evolution of this unique and nearby stellar system, as well as possible signatures of the past interaction between the MCs and the MW.

The Coordinated Radio and Infrared Survey for High-mass Star Formation. II. Source Catalog

The CORNISH project is the highest resolution radio continuum survey of the Galactic plane to date. It is the 5 GHz radio continuum part of a series of multi-wavelength surveys that focus on the northern GLIMPSE region (10° < l < 65°), observed by the Spitzer satellite in the mid-infrared. Observations with the Very Large Array in B and BnA configurations have yielded a 1.''5 resolution Stokes I map with a root mean square noise level better than 0.4 mJy beam 1. Here we describe the data-processing methods and data characteristics, and present a new, uniform catalog of compact radio emission. This includes an implementation of automatic deconvolution that provides much more reliable imaging than standard CLEANing. A rigorous investigation of the noise characteristics and reliability of source detection has been carried out. We show that the survey is optimized to detect emission on size scales up to 14'' and for unresolved sources the catalog is more than 90% complete at a flux density of 3.9 mJy. We have detected 3062 sources above a 7σ detection limit and present their ensemble properties. The catalog is highly reliable away from regions containing poorly sampled extended emission, which comprise less than 2% of the survey area. Imaging problems have been mitigated by down-weighting the shortest spacings and potential artifacts flagged via a rigorous manual inspection with reference to the Spitzer infrared data. We present images of the most common source types found: H II regions, planetary nebulae, and radio galaxies. The CORNISH data and catalog are available online at http://cornish.leeds.ac.uk.

Local star formation triggered by supernova shocks in magnetized diffuse neutral clouds

In this work, considering the impact of a supernova remnant (SNR) with a neutral magnetized cloud we derived analytically a set of conditions that are favourable for driving gravitational instability in the cloud and thus star formation. Using these conditions, we have built diagrams of the SNR radius, R(SNR), versus the initial cloud density, n(c), that constrain a domain in the parameter space where star formation is allowed. This work is an extension to previous study performed without considering magnetic fields (Melioli et al. 2006, hereafter Paper I). The diagrams are also tested with fully three-dimensional MHD radiative cooling simulations involving a SNR and a self-gravitating cloud and we find that the numerical analysis is consistent with the results predicted by the diagrams. While the inclusion of a homogeneous magnetic field approximately perpendicular to the impact velocity of the SNR with an intensity similar to 1 mu G within the cloud results only a small shrinking of the star formation zone in the diagram relative to that without magnetic field, a larger magnetic field (similar to 10 mu G) causes a significant shrinking, as expected. Though derived from simple analytical considerations these diagrams provide a useful tool for identifying sites where star formation could be triggered by the impact of a supernova blast wave. Applications of them to a few regions of our own Galaxy (e.g. the large CO shell in the direction of Cassiopeia, and the Edge Cloud 2 in the direction of the Scorpious constellation) have revealed that star formation in those sites could have been triggered by shock waves from SNRs for specific values of the initial neutral cloud density and the SNR radius. Finally, we have evaluated the effective star formation efficiency for this sort of interaction and found that it is generally smaller than the observed values in our own Galaxy (SFE similar to 0.01-0.3). This result is consistent with previous work in the literature and also suggests that the mechanism presently investigated, though very powerful to drive structure formation, supersonic turbulence and eventually, local star formation, does not seem to be sufficient to drive global star formation in normal star-forming galaxies, not even when the magnetic field in the neutral clouds is neglected.

Molecular gas in the galaxy M83 : Its distribution, kinematics, and relation to star formation

The barred spiral galaxy M83 (NGC5236) has been observed in the 12CO J=1–0 and J=2–1 millimetre lines with the Swedish-ESO Submillimetre Telescope (SEST). The sizes of the CO maps are 100×100, and they cover the entire optical disk. The CO emission is strongly peaked toward the nucleus. The molecular spiral arms are clearly resolved and can be traced for about 360º. The total molecular gas mass is comparable to the total Hi mass, but H2 dominates in the optical disk. Iso-velocity maps show the signature of an inclined, rotating disk, but also the effects of streaming motions along the spiral arms. The dynamical mass is determined and compared to the gas mass. The pattern speed is determined from the residual velocity pattern, and the locations of various resonances are discussed. The molecular gas velocity dispersion is determined, and a trend of decreasing dispersion with increasing galactocentric radius is found. A total gas (H2+Hi+He) mass surface density map is presented, and compared to the critical density for star formation of an isothermal gaseous disk. The star formation rate (SFR) in the disk is estimated using data from various star formation tracers. The different SFR estimates agree well when corrections for extinctions, based on the total gas mass map, are made. The radial SFR distribution shows features that can be associated with kinematic resonances. We also find an increased star formation efficiency in the spiral arms. Different Schmidt laws are fitted to the data. The star formation properties of the nuclear region, based on high angular resolution HST data, are also discussed.

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.

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).

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.

Role of galaxy mergers in cosmic star formation history

We present a morphology study of intermediate-redshift (0.2 < z < 1.2) luminous infrared galaxies (LIRGs) and general field galaxies in the GOODS fields using a revised asymmetry measurement method optimized for deep fields. By taking careful account of the importance of the underlying sky-background structures, our new method does not suffer from systematic bias and offers small uncertainties. By redshifting local LIRGs and low-redshift GOODS galaxies to different higher redshifts, we have found that the redshift dependence of the galaxy asymmetry due to surface-brightness dimming is a function of the asymmetry itself, with larger corrections for more asymmetric objects. By applying redshift-, infrared (IR)-luminosity- and optical-brightness-dependent asymmetry corrections, we have found that intermediate-redshift LIRGs generally show highly asymmetric morphologies, with implied merger fractions ~50% up to z = 1.2, although they are slightly more symmetric than local LIRGs. For general field galaxies, we find an almost constant relatively high merger fraction (20%-30%). The B-band luminosity functions (LFs) of galaxy mergers are derived at different redshifts up to z = 1.2 and confirm the weak evolution of the merger fraction after breaking the luminosity-density degeneracy. The IR LFs of galaxy mergers are also derived, indicating a larger merger fraction at higher IR luminosity. The integral of the merger IR LFs indicates a dramatic evolution of the merger-induced IR energy density [(1 + z)^~(5-6)], and that galaxy mergers start to dominate the cosmic IR energy density at z greater than or ~ 1.

On the metallicity dependence of the 24 μm luminosity as a star formation tracer

We investigate the use of the rest-frame 24 μm luminosity as an indicator of the star formation rate (SFR) in galaxies with different metallicities by comparing it to the (extinction-corrected) Hα luminosity. We carry out this analysis in two steps: First, we compare the emission from H (II) regions in different galaxies with metallicities between 12 + and 8.9. We find that the 24 μm and the extinction-corrected Hα luminosities from individual H (II) log (O/H) = 8.1 regions follow the same correlation for all galaxies, independent of their metallicity. Second, the role of metallicity is explored further for the integrated luminosity in a sample of galaxies with metallicities in the range of 12 +. For this sample we compare the 24 μm and Hα luminosities integrated over the entire galaxies log (O/ H) = 7.2-9.1 and find a lack of the 24 μm emission for a given Hα luminosity for low-metallicity objects, likely reflecting a low dust content. These results suggest that the 24 μm luminosity is a good metallicity-independent tracer for the SFR in individual H (II) regions. On the other hand, metallicity has to be taken into account when using the 24 μm luminosity as a tracer for the SFR of entire galaxies.