Lens Models of Herschel-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 µm bright candidate lensing systems identified by the Herschel Multi-tiered Extragalactic Survey and Herschel Astrophysical Terahertz Large Area Survey. 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 previous lensing models for submillimeter galaxies (SMGs). For four new sources that also have high-resolution submillimeter maps, we test for differential lensing between the stellar and dust components and find that the 880 µm magnification factor (µ_(880)) is ~1.5 times higher than the near-IR magnification factor (µ_(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

The composite nature of Dust-Obscured Galaxies (DOGs) at z∼2-3 in the COSMOS field: I. A far-infrared view

Dust-Obscured galaxies (DOGs) are bright 24μm-selected sources with extreme obscuration at optical wavelengths. They are typically characterized by a rising power-law continuum of hot dust (TD ∼ 200-1000K) in the near-IR indicating that their mid-IR luminosity is dominated by an active galactic nucleus (AGN). DOGs with a fainter 24 μm flux display a stellar bump in the near-IR and their mid-IR luminosity appears to be mainly powered by dusty star formation. Alternatively, it may be that the mid-IR emission arising from AGN activity is dominant but the torus is sufficiently opaque to make the near-IR emission from the AGN negligible with respect to the emission from the host component. In an effort to characterize the astrophysical nature of the processes responsible for the IR emission in DOGs, this paper exploits Herschel data (PACS + SPIRE) on a sample of 95 DOGs within the COSMOS field. We derive a wealth of far-IR properties (e.g., total IR luminosities; mid-to-far IR colours; dust temperatures and masses) based on SED fitting. Of particular interest are the 24 μm-bright DOGs (F 24μm >1mJy). They present bluer far-IR/mid-IR colours than the rest of the sample, unveiling the potential presence of an AGN. The AGN contribution to the total 8-1000μm flux increases as a function of the rest-frame 8μm-luminosity irrespective of the redshift. This confirms that faint DOGs (L 8μm < 10 12 Lo ) are dominated by star-formation while brighter DOGs show a larger contribution from an AGN

Constraining the Lyα escape fraction with far-infrared observations of Lyα emitters

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 f esc (Lyα) at these redshifts. Thus, if the LAEs f esc (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

Environment of the submillimeter-bright massive starburst HFLS3 at z∼z\sim6.34

We describe the search for Lyman-break galaxies (LBGs) near the sub-millimeter bright starburst galaxy HFLS3 at $z$$=$6.34 and a study on the environment of this massive galaxy during the end of reionization.We performed two independent selections of LBGs on images obtained with the \textit{Gran Telescopio Canarias} (GTC) and the \textit{Hubble Space Telescope} (HST) by combining non-detections in bands blueward of the Lyman-break and color selection. A total of 10 objects fulfilling the LBG selection criteria at $z$$>$5.5 were selected over the 4.54 and 55.5 arcmin$^2$ covered by our HST and GTC images, respectively. The photometric redshift, UV luminosity, and the star-formation rate of these sources were estimated with models of their spectral energy distribution. These $z$$\sim$6 candidates have physical properties and number densities in agreement with previous results. The UV luminosity function at $z$$\sim$6 and a Voronoi tessellation analysis of this field shows no strong evidence for an overdensity of relatively bright objects (m$_{F105W}$$<$25.9) associated with \textit{HFLS3}. However, the over-density parameter deduced from this field and the surface density of objects can not excluded definitively the LBG over-density hypothesis. Moreover we identified three faint objects at less than three arcseconds from \textit{HFLS3} with color consistent with those expected for $z$$\sim$6 galaxies. Deeper data are needed to confirm their redshifts and to study their association with \textit{HFLS3} and the galaxy merger that may be responsible for the massive starburst.Comment: 14 pages, 12 figures, accepted for publication in Ap

HerMES: the rest-frame UV emission and a lensing model for the z= 6.34 luminous dusty starburst galaxy HFLS3

We discuss the rest-frame ultraviolet emission from the starbursting galaxy HFLS3 at a redshift of 6.34. The galaxy was discovered in Herschel/SPIRE data due to its red color in the submillimeter wavelengths from 250 to 500 μm. Keck/NIRC2 K s -band adaptive optics imaging data showed two potential near-IR counterparts near HFLS3. Previously, the northern galaxy was taken to be in the foreground at z = 2.1, while the southern galaxy was assumed to be HFLS3's near-IR counterpart. The recently acquired Hubble/WFC3 and Advanced Camera for Surveys (ACS) imaging data show conclusively that both optically bright galaxies are in the foreground at z < 6. A new lensing model based on the Hubble imaging data and the millimeter-wave continuum emission yields a magnification factor of 2.2 ± 0.3, with a 95% confidence upper limit on the magnification of 3.5. When corrected for lensing, the instantaneous star formation rate is 1320 M ☉ yr–1, with the 95% confidence lower limit around 830 M ☉ yr–1. The dust and stellar masses of HFLS3 from the same spectral energy distribution (SED) models are at the level of 3 × 108 M ☉ and ~5 × 1010 M ☉, respectively, with large systematic uncertainties on assumptions related to the SED model. With Hubble/WFC3 images, we also find diffuse near-IR emission about 0.5 arcsec (~3 kpc) to the southwest of HFLS3 that remains undetected in the ACS imaging data. The emission has a photometric redshift consistent with either z ~ 6 or a dusty galaxy template at z ~ 2

Spitzer imaging of strongly lensed Herschel-selected dusty star-forming galaxies

We present the rest-frame optical spectral energy distribution (SED) 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.5 μm. Due to the spatial resolution of the IRAC observations at the level of 2″, 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.5 μm, once combined with Hubble/WFC3 and Keck/NIRC2 data, provide important constraints on the rest-frame optical SED 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 × 1010–4 × 1011 Me and star formation rates of around 100 Me yr−1 . This puts these lensed submillimeter 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

The rapid assembly of an elliptical galaxy of 400 billion solar masses at a redshift of 2.3

Stellar archeology shows that massive elliptical galaxies today formed rapidly about ten billion years ago with star formation rates above several hundreds solar masses per year (M_sun/yr). Their progenitors are likely the sub-millimeter-bright galaxies (SMGs) at redshifts (z) greater than 2. While SMGs' mean molecular gas mass of 5x10^10 M_sun can explain the formation of typical elliptical galaxies, it is inadequate to form ellipticals that already have stellar masses above 2x10^11 M_sun at z ~ 2. Here we report multi-wavelength high-resolution observations of a rare merger of two massive SMGs at z = 2.3. The system is currently forming stars at a tremendous rate of 2,000 M_sun/yr. With a star formation efficiency an order-of-magnitude greater than that of normal galaxies, it will quench the star formation by exhausting the gas reservoir in only ~200 million years. At a projected separation of 19 kiloparsecs, the two massive starbursts are about to merge and form a passive elliptical galaxy with a stellar mass of ~4x10^11 M_sun. Our observations show that gas-rich major galaxy mergers, concurrent with intense star formation, can form the most massive elliptical galaxies by z ~ 1.5.Comment: Appearing in Nature online on May 22 and in print on May 30. Submitted here is the accepted version (including the Supplementary Information), see nature.com for the final versio

ALMA observations of lensed Herschel sources: testing the dark matter halo paradigm

With the advent of wide-area submillimetre surveys, a large number of high-redshift gravitationally lensed dusty star-forming galaxies have been revealed. Because of the simplicity of the selection criteria for candidate lensed sources in such surveys, identified as those with S500 μm > 100 mJy, uncertainties associated with the modelling of the selection function are expunged. The combination of these attributes makes submillimetre surveys ideal for the study of strong lens statistics. We carried out a pilot study of the lensing statistics of submillimetre-selected sources by making observations with the Atacama Large Millimeter Array (ALMA) of a sample of strongly lensed sources selected from surveys carried out with the Herschel Space Observatory. We attempted to reproduce the distribution of image separations for the lensed sources using a halo mass function taken from a numerical simulation that contains both dark matter and baryons. We used three different density distributions, one based on analytical fits to the haloes formed in the EAGLE simulation and two density distributions [Singular Isothermal Sphere (SIS) and SISSA] that have been used before in lensing studies. We found that we could reproduce the observed distribution with all three density distributions, as long as we imposed an upper mass transition of ∼1013 M⊙ for the SIS and SISSA models, above which we assumed that the density distribution could be represented by a Navarro–Frenk–White profile. We show that we would need a sample of ∼500 lensed sources to distinguish between the density distributions, which is practical given the predicted number of lensed sources in the Herschel surveys

Herschel-ATLAS: A Binary HyLIRG Pinpointing a Cluster of Starbursting Protoellipticals

Panchromatic observations of the best candidate hyperluminous infrared galaxies from the widest Herschel extragalactic imaging survey have led to the discovery of at least four intrinsically luminous z = 2.41 galaxies across an 100 kpc region-a cluster of starbursting protoellipticals. Via subarcsecond interferometric imaging we have measured accurate gas and star formation surface densities. The two brightest galaxies span 3 kpc FWHM in submillimeter/radio continuum and CO J = 4-3, and double that in CO J = 1-0. The broad CO line is due partly to the multitude of constituent galaxies and partly to large rotational velocities in two counter-rotating gas disks-a scenario predicted to lead to the most intense starbursts, which will therefore come in pairs. The disks have Mdyn of several 10(sup 11) solar Mass, and gas fractions of 40%. Velocity dispersions are modest so the disks are unstable, potentially on scales commensurate with their radii: these galaxies are undergoing extreme bursts of star formation, not confined to their nuclei, at close to the Eddington limit. Their specific star formation rates place them greater than or approx. equal to 5 above the main sequence, which supposedly comprises large gas disks like these. Their high star formation efficiencies are difficult to reconcile with a simple volumetric star formation law. N-body and dark matter simulations suggest that this system is the progenitor of a B(inary)-type 10(sup 14.6) -solar mass cluster

Multi-wavelength lens construction of a Planck and Herschel-detected star-bursting galaxy

We present a source-plane reconstruction of a Herschel and Planck-detected gravitationally lensed dusty star-forming galaxy (DSFG) at z = 1.68 using Hubble, Submillimeter Array (SMA), and Keck observations. The background submillimeter galaxy (SMG) is strongly lensed by a foreground galaxy cluster at z = 0.997 and appears as an arc with a length of ∼15″ in the optical images. The continuum dust emission, as seen by SMA, is limited to a single knot within this arc. We present a lens model with source-plane reconstructions at several wavelengths to show the difference in magnification between the stars and dust, and highlight the importance of multi-wavelength lens models for studies involving lensed DSFGs. We estimate the physical properties of the galaxy by fitting the flux densities to model spectral energy distributions leading to a magnification-corrected starformation rate (SFR) of 390 ± 60 M yr−1 and a stellar mass of 1.1 ± 0.4 10 x 11 M. These values are consistent with high-redshift massive galaxies that have formed most of their stars already. The estimated gas-to-baryon fraction, molecular gas surface density, and SFR surface density have values of 0.43 ± 0.13, 350 ± 200 M pc−2, and ~ 12 7 M yr−1 kpc−2, respectively. The ratio of SFR surface density to molecular gas surface density puts this among the most star-forming systems, similar to other measured SMGs and local ULIRGs