Herschel Multi-tiered Extragalactic Survey Research Articles

Morphological Evolution of the Hosts of Far-infrared/Submillimeter Galaxies

Abstract We present a host morphological study of 1266 far-infrared galaxies (FIRGs) and submillimeter galaxies (SMGs) in the Cosmic Evolution Survey field using the F160W and F814W images obtained by the Hubble Space Telescope. The FIRGs and SMGs are selected from the Herschel Multi-tiered Extragalactic Survey and the SCUBA-2 Cosmology Legacy Survey, respectively. Their precise locations are based on the interferometry data from the Atacama Large Millimeter/submillimeter Array and the Very Large Array. These objects are mostly at 0.1 ≲ z ≲ 3. The SMGs can be regarded as the population at the high-redshift tail of the FIRGs. Most of our FIRGs/SMGs have a total infrared luminosity (L IR) in the regimes of luminous and ultraluminous infrared galaxies (LIRGs, L IR = 1011−12 L ⊙; ULIRGs, L IR > 1012 L ⊙). The hosts of the SMG ULIRGs, FIRG ULIRGs, and FIRG LIRGs are of sufficient numbers to allow for detailed analysis, and they are only modestly different in their stellar masses. Their morphological types are predominantly disk galaxies (type D) and irregular/interacting systems (type Irr/Int). There is a morphological transition at z ≈ 1.25 for the FIRG ULIRG hosts, above which the Irr/Int galaxies dominate and below which the D and Irr/Int galaxies have nearly the same contributions. The SMG ULIRG hosts seem to experience a similar transition. This suggests a shift in the relative importance of galaxy mergers/interactions versus secular gas accretions in “normal” disk galaxies as the possible triggering mechanisms of ULIRGs. The FIRG LIRG hosts are predominantly D galaxies over z = 0.25–1.25, where they are of sufficient statistics.

HELP: the Herschel Extragalactic Legacy Project

ABSTRACT We present the Herschel Extragalactic Legacy Project (HELP). This project collates, curates, homogenizes, and creates derived data products for most of the premium multiwavelength extragalactic data sets. The sky boundaries for the first data release cover 1270 deg2 defined by the Herschel SPIRE extragalactic survey fields; notably the Herschel Multi-tiered Extragalactic Survey (HerMES) and the Herschel Atlas survey (H-ATLAS). Here, we describe the motivation and principal elements in the design of the project. Guiding principles are transparent or ‘open’ methodologies with care for reproducibility and identification of provenance. A key element of the design focuses around the homogenization of calibration, meta data, and the provision of information required to define the selection of the data for statistical analysis. We apply probabilistic methods that extract information directly from the images at long wavelengths, exploiting the prior information available at shorter wavelengths and providing full posterior distributions rather than maximum-likelihood estimates and associated uncertainties as in traditional catalogues. With this project definition paper, we provide full access to the first data release of HELP; Data Release 1 (DR1), including a monolithic map of the largest SPIRE extragalactic field at 385 deg2 and 18 million measurements of PACS and SPIRE fluxes. We also provide tools to access and analyse the full HELP data base. This new data set includes far-infrared photometry, photometric redshifts, and derived physical properties estimated from modelling the spectral energy distributions over the full HELP sky. All the software and data presented is publicly available.

Red, redder, reddest: SCUBA-2 imaging of colour-selected Herschel sources

High-redshift, luminous, dusty star forming galaxies (DSFGs) constrain the extremity of galaxy formation theories. The most extreme are discovered through follow-up on candidates in large area surveys. Here we present 850 $\mu$m SCUBA-2 follow-up observations of 188 red DSFG candidates from the \textit{Herschel} Multi-tiered Extragalactic Survey (HerMES) Large Mode Survey, covering 274 deg$^2$. We detected 87 per cent with a signal-to-noise ratio $>$ 3 at 850~$\mu$m. We introduce a new method for incorporating the confusion noise in our spectral energy distribution fitting by sampling correlated flux density fluctuations from a confusion limited map. The new 850~$\mu$m data provide a better constraint on the photometric redshifts of the candidates, with photometric redshift errors decreasing from $\sigma_z/(1+z)\approx0.21$ to $0.15$. Comparison spectroscopic redshifts also found little bias ($\langle (z-z_{\rm spec})/(1+z_{\rm spec})\rangle = 0.08 $). The mean photometric redshift is found to be 3.6 with a dispersion of $0.4$ and we identify 21 DSFGs with a high probability of lying at $z > 4$. After simulating our selection effects we find number counts are consistent with phenomenological galaxy evolution models. There is a statistically significant excess of WISE-1 and SDSS sources near our red galaxies, giving a strong indication that lensing may explain some of the apparently extreme objects. Nevertheless, our sample should include examples of galaxies with the highest star formation rates in the Universe ($\gg10^3$ M$_\odot$yr$^{-1}$).

The Strong Gravitationally Lensed Herschel Galaxy HLock01: Optical Spectroscopy Reveals a Close Galaxy Merger with Evidence of Inflowing Gas

Abstract The submillimeter galaxy (SMG) HERMES J105751.1+573027 (hereafter HLock01) at z = 2.9574 ± 0.0001 is one of the brightest gravitationally lensed sources discovered in the Herschel Multi-tiered Extragalactic Survey. Apart from the high flux densities in the far-infrared, it is also extremely bright in the rest-frame ultraviolet (UV), with a total apparent magnitude m UV ≃ 19.7 mag. We report here deep spectroscopic observations with the Gran Telescopio Canarias of the optically bright lensed images of HLock01. Our results suggest that HLock01 is a merger system composed of the Herschel-selected SMG and an optically bright Lyman break-like galaxy (LBG), separated by only 3.3 kpc in projection. While the SMG appears very massive (M * ≃ 5 × 1011 M ⊙), with a highly extinguished stellar component (A V ≃ 4.3 ), the LBG is a young, lower-mass (M * ≃ 1 × 1010 M ⊙), but still luminous ( ) satellite galaxy. Detailed analysis of the high signal-to-noise ratio (S/N) rest-frame UV spectrum of the LBG shows complex kinematics of the gas, exhibiting both blueshifted and redshifted absorption components. While the blueshifted component is associated with strong galactic outflows from the massive stars in the LBG, as is common in most star-forming galaxies, the redshifted component may be associated with gas inflow seen along a favorable sightline to the LBG. We also find evidence of an extended gas reservoir around HLock01 at an impact parameter of 110 kpc, through the detection of C ii λλ1334 absorption in the red wing of a bright Lyα emitter at z ≃ 3.327. The data presented here highlight the power of gravitational lensing in high S/N studies to probe deeply into the physics of high-z star-forming galaxies.

HELP: star formation as a function of galaxy environment withHerschel

The Herschel Extragalactic Legacy Project (HELP) brings together a vast range of data from many astronomical observatories. Its main focus is on the Herschel data, which maps dust obscured star formation over 1300 deg$^2$. With this unprecedented combination of data sets, it is possible to investigate how the star formation vs stellar mass relation (main-sequence) of star-forming galaxies depends on environment. In this pilot study we explore this question between 0.1 < z < 3.2 using data in the COSMOS field. We estimate the local environment from a smoothed galaxy density field using the full photometric redshift probability distribution. We estimate star formation rates by stacking the SPIRE data from the Herschel Multi-tiered Extragalactic Survey (HerMES). Our analysis rules out the hypothesis that the main-sequence for star-forming systems is independent of environment at 1.5 < z < 2, while a simple model in which the mean specific star formation rate declines with increasing environmental density gives a better description. However, we cannot exclude a simple hypothesis in which the main-sequence for star-forming systems is independent of environment at z < 1.5 and z > 2. We also estimate the evolution of the star formation rate density in the COSMOS field and our results are consistent with previous measurements at z < 1.5 and z > 2 but we find a $1.4^{+0.3}_{-0.2}$ times higher peak value of the star formation rate density at $z \sim 1.9$.

SPITZERIMAGING OF STRONGLY LENSEDHERSCHEL-SELECTED DUSTY STAR-FORMING GALAXIES

We present the rest-frame optical spectral energy distribution and stellar masses of six Herschel- selected gravitationally lensed dusty, star-forming galaxies (DSFGs) at 1 < z < 3. These galaxies were first identified with Herschel/SPIRE imaging data from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). The targets were observed with Spitzer/IRAC at 3.6 and 4.5um. Due to the spatial resolution of the IRAC observations at the level of 2 arcseconds, the lensing features of a background DSFG in the near-infrared are blended with the flux from the foreground lensing galaxy in the IRAC imaging data. We make use of higher resolution Hubble/WFC3 or Keck/NIRC2 Adaptive Optics imaging data to fit light profiles of the foreground lensing galaxy (or galaxies) as a way to model the foreground components, in order to successfully disentangle the foreground lens and background source flux densities in the IRAC images. The flux density measurements at 3.6 and 4.5um, once combined with Hubble/WFC3 and Keck/NIRC2 data, provide important constraints on the rest-frame optical spectral energy distribution of the Herschel-selected lensed DSFGs. We model the combined UV- to millimeter-wavelength SEDs to establish the stellar mass, dust mass, star-formation rate, visual extinction, and other parameters for each of these Herschel-selected DSFGs. These systems have inferred stellar masses in the range 8 x 10^10 to 4 x 10^11 Msun and star-formation rates of around 100 Msun yr-1. This puts these lensed sub-millimeter systems well above the SFR-M* relation observed for normal star-forming galaxies at similar redshifts. The high values of SFR inferred for these systems are consistent with a major merger-driven scenario for star formation.

HerMES: ALMA IMAGING OFHERSCHEL-SELECTED DUSTY STAR-FORMING GALAXIES

The Herschel Multi-tiered Extragalactic Survey (HerMES) has identified large numbers of dusty star-forming galaxies (DSFGs) over a wide range in redshift. A detailed understanding of these DSFGs is hampered by the limited spatial resolution of Herschel. We present 870um 0.45" resolution imaging from the Atacama Large Millimeter/submillimeter Array (ALMA) of 29 HerMES DSFGs with far-infrared (FIR) flux densities in between the brightest of sources found by Herschel and fainter DSFGs found in ground-based sub-millimeter (sub-mm) surveys. We identify 62 sources down to the 5-sigma point-source sensitivity limit in our ALMA sample (sigma~0.2mJy), of which 6 are strongly lensed (showing multiple images) and 36 experience significant amplification (mu>1.1). To characterize the properties of the ALMA sources, we introduce and make use of uvmcmcfit, a publicly available Markov chain Monte Carlo analysis tool for interferometric observations of lensed galaxies. Our lens models tentatively favor intrinsic number counts for DSFGs with a steep fall off above 8mJy at 880um. Nearly 70% of the Herschel sources comprise multiple ALMA counterparts, consistent with previous research indicating that the multiplicity rate is high in bright sub-mm sources. Our ALMA sources are located significantly closer to each other than expected based on results from theoretical models as well as fainter DSFGs identified in the LABOCA ECDFS Submillimeter Survey. The high multiplicity rate and low projected separations argue in favor of interactions and mergers driving the prodigious emission from the brightest DSFGs as well as the sharp downturn above S_880=8mJy.

HERMES: CURRENT COSMIC INFRARED BACKGROUND ESTIMATES CAN BE EXPLAINED BY KNOWN GALAXIES AND THEIR FAINT COMPANIONS AT z &lt; 4

We report contributions to cosmic infrared background (CIB) intensities originating from known galaxies and their faint companions at submillimeter wavelengths. Using the publicly-available UltraVISTA catalog, and maps at 250, 350, and 500 {\mu}m from the \emph{Herschel} Multi-tiered Extragalactic Survey (HerMES), we perform a novel measurement that exploits the fact that uncatalogued sources may bias stacked flux densities --- particularly if the resolution of the image is poor --- and intentionally smooth the images before stacking and summing intensities. By smoothing the maps we are capturing the contribution of faint (undetected in K_S ~ 23.4) sources that are physically associated, or correlated, with the detected sources. We find that the cumulative CIB increases with increased smoothing, reaching 9.82 +- 0.78, 5.77 +- 0.43, and 2.32 +- 0.19$\, \rm nW m^{-2} sr^{-1}$ at 250, 350, and 500 {\mu}m at 300 arcsec FWHM. This corresponds to a fraction of the fiducial CIB of 0.94 +- 0.23, 1.07 +- 0.31, and 0.97 +- 0.26 at 250, 350, and 500 {\mu}m, where the uncertainties are dominated by those of the absolute CIB. We then propose, with a simple model combining parametric descriptions for stacked flux densities and stellar mass functions, that emission from galaxies with log(M/Msun) > 8.5 can account for the most of the measured total intensities, and argue against contributions from extended, diffuse emission. Finally, we discuss prospects for future survey instruments to improve the estimates of the absolute CIB levels, and observe any potentially remaining emission at z > 4.

LENS MODELS OFHERSCHEL-SELECTED GALAXIES FROM HIGH-RESOLUTION NEAR-IR OBSERVATIONS

We present Keck-Adaptive Optics and Hubble Space Telescope high resolution near-infrared (IR) imaging for 500 um-bright candidate lensing systems identified by the Herschel Multi-tiered Extra-galactic Survey (HerMES) and Herschel Astrophysical Terahertz Survey (H-ATLAS). Out of 87 candidates with near-IR imaging, 15 (~17%) display clear near-IR lensing morphologies. We present near-IR lens models to reconstruct and recover basic rest-frame optical morphological properties of the background galaxies from 12 new systems. Sources with the largest near-IR magnification factors also tend to be the most compact, consistent with the size bias predicted from simulations and pre- vious lensing models for sub-millimeter galaxies. For four new sources that also have high-resolution sub-mm maps, we test for differential lensing between the stellar and dust components and find that the 880 um magnification factor (u_880) is ~1.5 times higher than the near-IR magnification factor (u_NIR), on average. We also find that the stellar emission is ~2 times more extended in size than dust. The rest-frame optical properties of our sample of Herschel-selected lensed SMGs are consistent with those of unlensed SMGs, which suggests that the two populations are similar.

HerMES: point source catalogues from Herschel-SPIRE observations II★

The Herschel Multi-tiered Extragalactic Survey (HerMES) is the largest Guaranteed Time Key Programme on the Herschel Space Observatory. With a wedding cake survey strategy, it consists of nested fields with varying depth and area totalling ~380 deg^2. In this paper, we present deep point source catalogues extracted from Herschel-SPIRE observations of all HerMES fields, except for the later addition of the 270 deg^2 HeLMS field. These catalogues constitute the second Data Release (DR2) made in October 2013. A subset of these catalogues, which consists of bright sources extracted from Herschel-SPIRE observations completed by May 1, 2010 (covering ~ 74 deg^2) were released earlier in the first extensive Data Release (DR1) in March 2012. Two different methods are used to generate the point source catalogues, the SUSSEXtractor (SXT) point source extractor used in two earlier data releases (EDR and EDR2) and a new source detection and photometry method. The latter combines an iterative source detection algorithm, StarFinder (SF), and a De-blended SPIRE Photometry (DESPHOT) algorithm. We use end-to-end Herschel-SPIRE simulations with realistic number counts and clustering properties to characterise basic properties of the point source catalogues, such as the completeness, reliability, photometric and positional accuracy. Over 500, 000 catalogue entries in HerMES fields (except HeLMS) are released to the public through the HeDAM website (http://hedam.oamp.fr/herMES).

CONSTRAINING THE Lyα ESCAPE FRACTION WITH FAR-INFRARED OBSERVATIONS OF Lyα EMITTERS

We study the far-infrared properties of 498 Lyman Alpha Emitters (LAEs) at z=2.8, 3.1 and 4.5 in the Extended Chandra Deep Field-South, using 250, 350 and 500 micron data from the Herschel Multi-tiered Extragalactic Survey (HerMES) and 870 micron data from the LABOCA ECDFS Submillimeter Survey (LESS). 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\sigma$ depths of ~0.1 to 0.4 mJy. The LAEs are also undetected at $\ge3\sigma$ in the stacks, although a $2.5\sigma$ signal is observed at 870 micron for the z=2.8 sources. We consider a wide range of far-infrared spectral energy distributions (SEDs), including a 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$\alpha$ and UV emission to determine lower limits on the LAEs Ly$\alpha$ escape fraction ($f_{\rm esc}($Ly$\alpha$)). For the Sd SED template, the inferred LAEs $f_{\rm esc}($Ly$\alpha$) are $\gtrsim30%$ ($1\sigma$) at z=2.8, 3.1 and 4.5, which are all significantly higher than the global $f_{\rm esc}($Ly$\alpha$) at these redshifts. Thus, if the LAEs $f_{\rm esc}($Ly$\alpha$) 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_{\rm esc}($Ly$\alpha$) of ~10 to 20% ($1\sigma$), all of which are slightly higher than the global evolution of $f_{\rm esc}($Ly$\alpha$) but consistent with it at the 2 to 3$\sigma$ level.

Herschel Multitiered Extragalactic Survey: clusters of dusty galaxies uncovered by Herschel★ and Planck†

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 redshifts and 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 >1000 M⊙ yr−1. The high-redshift clusters discovered in these 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.

HerMES: dust attenuation and star formation activity in ultraviolet-selected samples from z∼ 4 to ∼ 1.5★

We study the link between observed ultraviolet (UV) luminosity, stellar mass and dust attenuation within rest-frame UV-selected samples at z ∼ 4, ∼ 3 and ∼1.5. We measure by stacking at 250, 350 and 500 μm in the Herschel/Spectral and Photometric Imaging Receiver images from the Herschel Multi-Tiered Extragalactic Survey (HerMES) program the average infrared luminosity as a function of stellar mass and UV luminosity. We find that dust attenuation is mostly correlated with stellar mass. There is also a secondary dependence with UV luminosity: at a given UV luminosity, dust attenuation increases with stellar mass, while at a given stellar mass it decreases with UV luminosity. We provide new empirical recipes to correct for dust attenuation given the observed UV luminosity and the stellar mass. Our results also enable us to put new constraints on the average relation between star formation rate (SFR) and stellar mass at z ∼ 4, ∼3 and ∼1.5. The SFR–stellar mass relations are well described by power laws (⁠|${\rm SFR}\propto M_*^{0.7}$|⁠), with the amplitudes being similar at z ∼ 4 and ∼3, and decreasing by a factor of 4 at z ∼ 1.5 at a given stellar mass. We further investigate the evolution with redshift of the specific SFR. Our results are in the upper range of previous measurements, in particular at z ∼ 3, and are consistent with a plateau at 3 < z < 4. Current model predictions (either analytic, semi-analytic or hydrodynamic) are inconsistent with these values, as they yield lower predictions than the observations in the redshift range we explore. We use these results to discuss the star formation histories of galaxies in the framework of the main sequence of star-forming galaxies. Our results suggest that galaxies at high redshift (2.5 < z < 4) stay around 1 Gyr on the main sequence. With decreasing redshift, this time increases such that z = 1 main-sequence galaxies with |$10^{8}<M_*/\rm {M}_{{\odot }}<10^{10}$| stay on the main sequence until z = 0.

GRAVITATIONAL LENS MODELS BASED ON SUBMILLIMETER ARRAY IMAGING OFHERSCHEL-SELECTED STRONGLY LENSED SUB-MILLIMETER GALAXIES ATz&gt; 1.5

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 (multiply imaged), 4 are moderately lensed (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 500um 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.9x10^12 L_sun) and two decades in FIR luminosity surface density (median Sigma_FIR=6.0x10^11 L_sun 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.

AGN and Star Formation in HerMES-IRS sources

AbstractOne of the remaining open issues in the context of the analysis of Active Galactic Nuclei (AGN) is the evidence that nuclear gravitational accretion is often accompanied by a concurrent starburst activity. We developed a spectral energy distribution (SED) fitting technique to derive simultaneously the physical properties of active galaxies and coexisting starbursts making the best use of Spitzer and Herschel IR observations. We apply the SED fitting procedure to a large sample of extragalactic sources representing the HerMES (Herschel/Multi-tiered Extragalactic Survey) population with IRS spectra with a plethora of multi-wavelength data in order to study the impact of a possible presence of an AGN on the host galaxy's properties. We analyze the star formation rate (SFR) in conncetion to the presence of an AGN and compared the properties of the hot (AGN) and cold (starburst) dust component. Our findings are consistent with no evidence for the presence of an AGN affecting the star formation processes of the host galaxies.

Herschel reveals the obscured star formation in HiZELS H emitters at z = 1.47

We describe the far-infrared (far-IR; rest-frame 8–1000-μm) properties of a sample of 443 Hα-selected star-forming galaxies in the Cosmic Evolution Survey (COSMOS) and Ultra Deep Survey (UDS) fields detected by the High-redshift Emission Line Survey (HiZELS) imaging survey. Sources are identified using narrow-band filters in combination with broad-band photometry to uniformly select Hα (and [O ii] if available) emitters in a narrow redshift slice at z = 1.47 ± 0.02. We use a stacking approach in Spitzer-MIPS mid-IR, Herschel-PACS/SPIRE far-IR [from the PACS Evolutionary Prove (PEP) and Herschel Multi-tiered Extragalactic Survey (HerMES)] and AzTEC mm-wave images to describe their typical far-IR properties. We find that HiZELS galaxies with observed Hα luminosities of L(Hα)obs ≈ 108.1-9.1 L⊙ ( ≈ 1041.7-42.7 erg s−1) have bolometric far-IR luminosities of typical luminous IR galaxies, |$L(8{\rm -}1000\,\mu {\rm m})\approx 10^{11.41^{+0.04}_{-0.06}} $| L⊙. Combining the Hα and far-IR luminosities, we derive median star formation rates (SFRs) of SFRHα, FIR = 32 ± 5 M⊙ yr−1 and Hα extinctions of AHα = 1.0 ± 0.2 mag. Perhaps surprisingly, little difference is seen in typical HiZELS extinction levels compared to local star-forming galaxies. We confirm previous empirical stellar mass (M*) to AHα relations and the little or no evolution up to z = 1.47. For HiZELS galaxies (or similar samples) we provide an empirical parametrization of the SFR as a function of rest-frame (u − z) colours and 3.6-μm photometry – a useful proxy to aid in the absence of far-IR detections in high-z galaxies. We find that the observed Hα luminosity is a dominant SFR tracer when rest-frame (u − z) colours are ≲0.9 mag or when Spitzer-3.6-μm photometry is fainter than 22 mag (Vega) or when stellar masses are lower than 109.7 M⊙. We do not find any correlation between the [O ii]/Hα and far-IR luminosity, suggesting that this emission line ratio does not trace the extinction of the most obscured star-forming regions, especially in massive galaxies where these dominate. The luminosity-limited HiZELS sample tends to lie above of the so-called main sequence for star-forming galaxies, especially at low stellar masses, indicating high star formation efficiencies in these galaxies. This work has implications for SFR indicators and suggests that obscured star formation is linked to the assembly of stellar mass, with deeper potential wells in massive galaxies providing dense, heavily obscured environments in which stars can form rapidly.

The roles of star formation and AGN activity of IRS sources in the HerMES fields

In this work, we explore the impact of the presence of an active galactic nucleus (AGN) on the mid- and far-infrared (IR) properties of galaxies as well as the effects of simultaneous AGN and starburst activity in the same galaxies. To do this, we apply a multicomponent, multiband spectral synthesis technique to a sample of 250 μm selected galaxies of the Herschel Multi-tiered Extragalactic Survey (HerMES), with Infrared Spectrograph (IRS) spectra available for all galaxies. Our results confirm that the inclusion of the IRS spectra plays a crucial role in the spectral analysis of galaxies with an AGN component improving the selection of the best-fitting hot dust (torus) model. We find a correlation between the obscured star formation rate, SFRIR, derived from the IR luminosity of the starburst component, and SFRPAH, derived from the luminosity of the PAH features, LPAH, with SFRFIR taking higher values than SFRPAH. The correlation is different for AGN- and starburst-dominated objects. The ratio of LPAH to that of the starburst component, LPAH/LSB, is almost constant for AGN-dominated objects but decreases with increasing LSB for starburst-dominated objects. SFRFIR increases with the accretion luminosity, Lacc, with the increase less prominent for the very brightest, unobscured AGN-dominated sources. We find no correlation between the masses of the hot (AGN-heated) and cold (starburst-heated) dust components. We interpret this as a non-constant fraction of gas driven by the gravitational effects to the AGN while the starburst is ongoing. We also find no evidence of the AGN affecting the temperature of the cold dust component, though this conclusion is mostly based on objects with a non-dominant AGN component. We conclude that our findings do not provide evidence that the presence of AGN affects the star formation process in the host galaxy, but rather that the two phenomena occur simultaneously over a wide range of luminosities.

The Herschel★ PEP/HerMES luminosity function – I. Probing the evolution of PACS selected Galaxies to z ≃ 4

We exploit the deep and extended far-IR data sets (at 70, 100 and 160 μm) of the Herschel Guaranteed Time Observation (GTO) PACS Evolutionary Probe (PEP) Survey, in combination with the Herschel Multi-tiered Extragalactic Survey data at 250, 350 and 500 μm, to derive the evolution of the rest-frame 35-, 60-, 90- and total infrared (IR) luminosity functions (LFs) up to z ∼ 4. We detect very strong luminosity evolution for the total IR LF (LIR ∝ (1 + z)3.55 ± 0.10 up to z ∼ 2, and ∝ (1 + z)1.62 ± 0.51 at 2 < z ≲ 4) combined with a density evolution (∝(1 + z)−0.57 ± 0.22 up to z ∼ 1 and ∝ (1 + z)−3.92 ± 0.34 at 1 < z ≲ 4). In agreement with previous findings, the IR luminosity density (ρIR) increases steeply to z ∼ 1, then flattens between z ∼ 1 and z ∼ 3 to decrease at z ≳ 3. Galaxies with different spectral energy distributions, masses and specific star formation rates (SFRs) evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to ρIR. Galaxies occupying the well-established SFR–stellar mass main sequence (MS) are found to dominate both the total IR LF and ρIR at all redshifts, with the contribution from off-MS sources (≥0.6 dex above MS) being nearly constant (∼20 per cent of the total ρIR) and showing no significant signs of increase with increasing z over the whole 0.8 < z < 2.2 range. Sources with mass in the range 10 ≤ log(M/M⊙) ≤ 11 are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z (≳2) LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the Super Massive Black Holes (SMBH) grows and is obscured by dust (possibly triggered by a major merging event), is followed by an AGN-dominated phase, then evolving towards a local elliptical. On the other hand, moderately star-forming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

Connecting stellar mass and star-formation rate to dark matter halo mass out to z ∼ 2

We have constructed an extended halo model (EHM) which relates the total stellar mass and star-formation rate (SFR) to halo mass (M_h). An empirical relation between the distribution functions of total stellar mass of galaxies and host halo mass, tuned to match the spatial density of galaxies over 0<z<2 and the clustering properties at z~0, is extended to include two different scenarios describing the variation of SFR on M_h. We also present new measurements of the redshift evolution of the average SFR for star-forming galaxies of different stellar mass up to z=2, using data from the Herschel Multi-tiered Extragalactic Survey (HerMES) for infrared-bright galaxies. Combining the EHM with the halo accretion histories from numerical simulations, we trace the stellar mass growth and star-formation history in halos spanning a range of masses. We find that: (1) The intensity of the star-forming activity in halos in the probed mass range has steadily decreased from z~2 to 0; (2) At a given epoch, halos in the mass range between a few times 10^{11} M_Sun and a few times 10^{12} M_Sun are the most efficient at hosting star formation; (3) The peak of SFR density shifts to lower mass halos over time; (4) Galaxies that are forming stars most actively at z~2 evolve into quiescent galaxies in today's group environments, strongly supporting previous claims that the most powerful starbursts at z~2 are progenitors of today's elliptical galaxies.

Inferring the mass of submillimetre galaxies by exploiting their gravitational magnification of background galaxies

Dust emission at sub-millimetre wavelengths allows us to trace the early phases of star formation in the Universe. In order to understand the physical processes involved in this mode of star formation, it is essential to gain knowledge about the dark matter structures - most importantly their masses - that sub-millimetre galaxies live in. Here we use the magnification effect of gravitational lensing to determine the average mass and dust content of sub-millimetre galaxies with 250mu flux densities of S_250>15mJy selected using data from the Herschel Multi-tiered Extragalactic Survey. The positions of hundreds of sub-millimetre foreground lenses are cross-correlated with the positions of background Lyman-break galaxies at z~3-5 selected using optical data from the Canada-France Hawaii Telescope Legacy Survey. We detect a cross-correlation signal at the 7-sigma level over a sky area of one square degree, with ~80% of this signal being due to magnification, whereas the remaining ~20% comes from dust extinction. Adopting some simple assumptions for the dark matter and dust profiles and the redshift distribution enables us to estimate the average mass of the halos hosting the sub-millimetre galaxies to be log(M_200/M_sun)=13.17+0.05-0.08(stat.) and their average dust mass fraction (at radii of >10kpc) to be M_dust/M_200~6x10^-5. This supports the picture that sub-millimetre galaxies are dusty, forming stars at a high rate, reside in massive group-sized halos, and are a crucial phase in the assembly and evolution of structure in the Universe.