We quantify the fraction of the cosmic infrared background (CIB) that originates from galaxies identified in the UV/optical/near-infrared by stacking 81,250 (∼35.7 arcmin-2) K-selected sources (K AB < 24.0) split according to their rest-frame U-V versus V-J colors into 72,216 star-forming and 9034 quiescent galaxies, on maps from Spitzer/MIPS (24 μm), Herschel/PACS (100, 160 μm), Herschel/SPIRE (250, 350, 500 μm), and AzTEC (1100 μm). The fraction of the CIB resolved by our catalog is (69% ± 15%) at 24 μm, (78% ± 17%) at 70 μm, (58% ± 13%) at 100 μm, (78% ± 18%) at 160 μm, (80% ± 17%) at 250 μm, (69% ± 14%) at 350 μm, (65% ± 12%) at 500 μm, and (45% ± 8%) at 1100 μm. Of that total, about 95% originates from star-forming galaxies, while the remaining 5% is from apparently quiescent galaxies. The CIB at λ ≲ 200 μm appears to be sourced predominantly from galaxies at z ≲ 1, while at λ ≳ 200 μm the bulk originates from 1 ≲ z ≲ 2. Galaxies with stellar masses log(M/M ) = 9.5-11 are responsible for the majority of the CIB, with those in the log(M/M ) = 9.5-10 bin contributing mostly at λ < 250 μm, and those in the log(M/M ) = 10-11 bin dominating at λ > 350 μm. The contribution from galaxies in the log(M/M ) = 9.0-9.5 (lowest) and log(M/M ) = 11.0-12.0 (highest) stellar-mass bins contribute the least - both of order 5% - although the highest stellar-mass bin is a significant contributor to the luminosity density at z ≳ 2. The luminosities of the galaxies responsible for the CIB shifts from combinations of "normal" and luminous infrared galaxies (LIRGs) at λ ≲ 160 μm, to LIRGs at 160 ≲ λ ≲ 500 μm, to finally LIRGs and ultra-luminous infrared galaxies at λ ≳ 500 μm. Stacking analyses were performed using SIMSTACK, a novel algorithm designed to account for possible biases in the stacked flux density due to clustering. It is made available to the public at www.astro.caltech.edu/∼viero/viero-homepage/ toolbox.html. © 2013. The American Astronomical Society. All rights reserved..

The potential for Planck to detect clusters of dusty, star-forming galaxies at z > 1 is tested by examining the Herschel-SPIRE images of Planck Early Release Compact Source Catalog sources lying in fields observed by the Herschel Multitiered Extragalactic Survey. Of the 16 Planck sources that lie in the ~90 sq. deg. examined, we find that 12 are associated with single bright Herschel sources. The remaining four are associated with overdensities of Herschel sources, making them candidate clusters of dusty, star-forming galaxies. We use complementary optical/near-IR data for these 'clumps' to test this idea, and find evidence for the presence of galaxy clusters in all four cases. We use photometric redshiftsand red sequence galaxies to estimate the redshifts of these clusters, finding that they range from 0.8 to 2.3. These redshifts imply that the Herschel sources in these clusters, which contribute to the detected Planck flux, are forming stars very rapidly, with typical total cluster star formation rates > 1000M⊙ yr-1. The high-redshift clusters discovered inthese observations are used to constrain the epoch of cluster galaxy formation, finding that the galaxies in our clusters are 1-1.5 Gyr old at z ~ 1-2. Prospects for the discovery of further clusters of dusty galaxies are discussed, using not only all sky Planck surveys, but also deeper, smaller area, Herschel surveys. © 2014 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society.

We study the far-infrared properties of 498 Lyα emitters (LAEs) at z = 2.8, 3.1, and 4.5 in the Extended Chandra Deep Field-South, using 250, 350, and 500 μm data from the Herschel Multi-tiered Extragalactic Survey and 870 μm data from the LABOCA ECDFS Submillimeter Survey. None of the 126, 280, or 92 LAEs at z = 2.8, 3.1, and 4.5, respectively, are individually detected in the far-infrared data. We use stacking to probe the average emission to deeper flux limits, reaching 1σ depths of ̃0.1 to 0.4 mJy. The LAEs are also undetected at ≥3σ in the stacks, although a 2.5σ signal is observed at 870 μm for the z = 2.8 sources. We consider a wide range of far-infrared spectral energy distributions (SEDs), including an M82 and an Sd galaxy template, to determine upper limits on the far-infrared luminosities and far-infrared-derived star formation rates of the LAEs. These star formation rates are then combined with those inferred from the Lyα and UV emission to determine lower limits on the LAEs' Lyα escape fraction (f esc(Lyα)). For the Sd SED template, the inferred LAEs f esc(Lyα) are ≳ 30% (1σ) at z = 2.8, 3.1, and 4.5, which are all significantly higher than the global fesc(Lyα) at these redshifts. Thus, if the LAEs fesc(Lyα) follows the global evolution, then they have warmer far-infrared SEDs than the Sd galaxy template. The average and M82 SEDs produce lower limits on the LAE f esc(Lyα) of ̃10%-20% (1σ), all of which are slightly higher than the global evolution of f esc(Lyα), but consistent with it at the 2σ-3σ level. © 2014. The American Astronomical Society. All rights reserved.

Strong gravitational lenses are now being routinely discovered in wide-field surveys at (sub-)millimeter wavelengths. We present Submillimeter Array (SMA) high-spatial resolution imaging and Gemini-South and Multiple Mirror Telescope optical spectroscopy of strong lens candidates discovered in the two widest extragalactic surveys conducted by the Herschel Space Observatory: the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). From a sample of 30 Herschel sources with S 500 > 100 mJy, 21 are strongly lensed (i.e., multiply imaged), 4 are moderately lensed (i.e., singly imaged), and the remainder require additional data to determine their lensing status. We apply a visibility-plane lens modeling technique to the SMA data to recover information about the masses of the lenses as well as the intrinsic (i.e., unlensed) sizes (r half) and far-infrared luminosities (L FIR) of the lensed submillimeter galaxies (SMGs). The sample of lenses comprises primarily isolated massive galaxies, but includes some groups and clusters as well. Several of the lenses are located at z lens > 0.7, a redshift regime that is inaccessible to lens searches based on Sloan Digital Sky Survey spectroscopy. The lensed SMGs are amplified by factors that are significantly below statistical model predictions given the 500 μm flux densities of our sample. We speculate that this may reflect a deficiency in our understanding of the intrinsic sizes and luminosities of the brightest SMGs. The lensed SMGs span nearly one decade in L FIR (median L FIR = 7.9 × 10 12 L ) and two decades in FIR luminosity surface density (median ΣFIR = 6.0 × 1011 L kpc-2). The strong lenses in this sample and others identified via (sub-)mm surveys will provide a wealth of information regarding the astrophysics of galaxy formation and evolution over a wide range in redshift. © 2013. The American Astronomical Society. All rights reserved..

We present a method for selecting z > 4 dusty, star-forming galaxies (DSFGs) using Herschel/Spectral and Photometric Imaging Receiver 250/350/500 μm flux densities to search for red sources. We apply this method to 21 deg2 of data from the HerMES survey to produce a catalog of 38 high-z candidates. Follow-up of the first five of these sources confirms that this method is efficient at selecting high-z DSFGs, with 4/5 at z = 4.3-6.3 (and the remaining source at z = 3.4), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 μm) and in single-band surveys, shows that our method is much more efficient at selecting high-z DSFGs, in the sense that a much larger fraction are at z > 3. Correcting for the selection completeness and purity, we find that the number of bright (S 500 μm ≥ 30 mJy), red Herschel sources is 3.3 ± 0.8 deg-2. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-z DSFGs is similar to that at z ∼ 2, rest-frame UV based studies may be missing a significant component of the star formation density at z = 4-6, even after correction for extinction. © 2014. The American Astronomical Society. All rights reserved.