Discovery of "Warm Dust" galaxies in clusters at z ~ 0.3: evidence for stripping of cool dust in the dense environment?

Using far-infrared imaging from the "Herschel Lensing Survey," we derive dust properties of spectroscopically confirmed cluster member galaxies within two massive systems at z ~ 0.3: the merging Bullet Cluster and the more relaxed MS2137.3-2353. Most star-forming cluster sources (~90%) have characteristic dust temperatures similar to local field galaxies of comparable infrared (IR) luminosity (T_dust ~ 30 K). Several sub-luminous infrared galaxy (LIRG; L_IR < 10^11 L_☉) Bullet Cluster members are much warmer (T_dust > 37 K) with far-infrared spectral energy distribution (SED) shapes resembling LIRG-type local templates. X-ray and mid-infrared data suggest that obscured active galactic nuclei do not contribute significantly to the infrared flux of these "warm dust" galaxies. Sources of comparable IR luminosity and dust temperature are not observed in the relaxed cluster MS2137, although the significance is too low to speculate on an origin involving recent cluster merging. "Warm dust" galaxies are, however, statistically rarer in field samples (>3σ), indicating that the responsible mechanism may relate to the dense environment. The spatial distribution of these sources is similar to the whole far-infrared bright population, i.e., preferentially located in the cluster periphery, although the galaxy hosts tend toward lower stellar masses (M_* < 10^10 M_☉). We propose dust stripping and heating processes which could be responsible for the unusually warm characteristic dust temperatures. A normal star-forming galaxy would need 30%-50% of its dust removed (preferentially stripped from the outer reaches, where dust is typically cooler) to recover an SED similar to a "warm dust" galaxy. These progenitors would not require a higher IR luminosity or dust mass than the currently observed normal star-forming population.

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.

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

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.

Exploring the star forming regions in vigorous star forming galaxies at Z=0.84

We analyze the properties of star forming regions in a sample of star forming galaxies at z = 0.84. Star forming regions are extracted from B band ACS-HST images. Previously we have substracted a model of the galaxy, fitting a bulged+disk model to the whole galaxy. Special care has been taken masking the star forming regions in the model fitting procedure, yielding more reliable results. We present here the properties of these star forming regions.

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.

The AGN, star-forming, and morphological properties of luminous ir-bright/optically-faint galaxies

We present the active galactic nucleus (AGN), star-forming, and morphological properties of a sample of 13 MIR-luminous (∫_24 700 μJy) IR-bright/optically-faint galaxies (IRBGs, ∫_24/f_R≲ 1000). While these z ∼ 2 sources were drawn from deep Chandra fields with >200 ks X-ray coverage, only seven are formally detected in the X-ray and four lack X-ray emission at even the 2σ level. Spitzer InfraRed Spectrograph (IRS) spectra, however, confirm that all of the sources are AGN-dominated in the mid-IR, although half have detectable polycyclic aromatic hydrocarbon (PAH) emission responsible for ∼25% of their mid-infrared flux density. When combined with other samples, this indicates that at least 30%–40% of luminous IRBGs have star formation rates in the ultraluminous infrared galaxy (ULIRG) range (∼100–2000 M_⨀ yr^−1). X-ray hardness ratios and MIR to X-ray luminosity ratios indicate that all members of the sample contain heavily X-ray obscured AGNs, 80% of which are candidates to be Compton thick. Furthermore, the mean X-ray luminosity of the sample, log L_2–10 keV(erg s^−1) ∼44.6, indicates that these IRBGs are Type 2 QSOs, at least from the X-ray perspective. While those sources most heavily obscured in the X-ray are also those most likely to display strong silicate absorption in the mid-IR, silicate absorption does not always accompany X-ray obscuration. Finally, ∼70% of the IRBGs are merger candidates, a rate consistent with that of sub-mm galaxies (SMGs), although SMGs appear to be physically larger than IRBGs. These characteristics are consistent with the proposal that these objects represent a later, AGN-dominated, and more relaxed evolutionary stage following soon after the star-formation-dominated one represented by the SMGs.

The far-infrared/submillimeter properties of galaxies located behind the Bullet cluster

The Herschel Lensing Survey (HLS) takes advantage of gravitational lensing by massive galaxy clusters to sample a population of high-redshift galaxies which are too faint to be detected above the confusion limit of current far-infrared/submillimeter telescopes. Measurements from 100-500 μm bracket the peaks of the far-infrared spectral energy distributions of these galaxies, characterizing their infrared luminosities and star formation rates. We introduce initial results from our science demonstration phase observations, directed toward the Bullet cluster (1E0657-56). By combining our observations with LABOCA 870 μm and AzTEC 1.1 mm data we fully constrain the spectral energy distributions of 19 MIPS 24 μm-selected galaxies which are located behind the cluster. We find that their colors are best fit using templates based on local galaxies with systematically lower infrared luminosities. This suggests that our sources are not like local ultra-luminous infrared galaxies in which vigorous star formation is contained in a compact highly dust-obscured region. Instead, they appear to be scaled up versions of lower luminosity local galaxies with star formation occurring on larger physical scales.

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.