Search

Starburst galaxies at the peak of cosmic star formation are among the most extreme starforming engines in the universe, producing stars over ~100 Myr. The star formation rates of these galaxies, which exceed 100 $M_\odot$ per year, require large reservoirs of cold molecular gas to be delivered to their cores, despite strong feedback from stars or active galactic nuclei. Starburst galaxies are therefore ideal targets to unravel the critical interplay between this feedback and the growth of a galaxy. The methylidyne cation, CH$^+$, is a most useful molecule for such studies because it cannot form in cold gas without supra-thermal energy input, so its presence highlights dissipation of mechanical energy or strong UV irradiation. Here, we report the detection of CH$^+$(J=1-0) emission and absorption lines in the spectra of six lensed starburst galaxies at redshifts z~2.5. This line has such a high critical density for excitation that it is emitted only in very dense ($>10^5$ cm$^{-3}$) gas, and is absorbed in low-density gas. We find that the CH$^+$ emission lines, which are broader than 1000 km s$^{-1}$, originate in dense shock waves powered by hot galactic winds. The CH$^+$ absorption lines reveal highly turbulent reservoirs of cool ($T\sim 100$K), low-density gas, extending far outside (>10 kpc) the starburst cores (radii <1 kpc). We show that the galactic winds sustain turbulence in the 10 kpc-scale environments of the starburst cores, processing these environments into multi-phase, gravitationally bound reservoirs. However, the mass outflow rates are found to be insufficient to balance the star formation rates. Another mass input is therefore required for these reservoirs, which could be provided by on-going mergers or cold stream accretion. Our results suggest that galactic feedback, coupled jointly to turbulence and gravity, extends the starburst phase instead of quenching it., Comment: Published as a Nature Letter

We present a sample of 80 candidate strongly lensed galaxies with flux density above 100mJy at 500{\mu}m extracted from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), over an area of 600 square degrees. Available imaging and spectroscopic data allow us to confirm the strong lensing in 20 cases and to reject it in one case. For other 8 objects the lensing scenario is strongly supported by the presence of two sources along the same line of sight with distinct photometric redshifts. The remaining objects await more follow-up observations to confirm their nature. The lenses and the background sources have median redshifts z_L = 0.6 and z_S = 2.5, respectively, and are observed out to z_L = 1.2 and z_S = 4.2. We measure the number counts of candidate lensed galaxies at 500{\mu}m and compare them with theoretical predictions, finding a good agreement for a maximum magnification of the background sources in the range 10-20. These values are consistent with the magnification factors derived from the lens modelling of individual systems. The catalogue presented here provides sub- mm bright targets for follow-up observations aimed at exploiting gravitational lensing to study with un-precedented details the morphological and dynamical properties of dusty star forming regions in z >~ 1.5 galaxies., Comment: 25 pages, 17 figures, accepted for publication in MNRAS

Until recently, only a handful of dusty, star-forming galaxies (DSFGs) were known at $z>4$, most of them significantly amplified by gravitational lensing. Here, we have increased the number of such DSFGs substantially, selecting galaxies from the uniquely wide 250-, 350- and 500-$\mu$m Herschel-ATLAS imaging survey on the basis of their extremely red far-infrared colors and faint 350- and 500-$\mu$m flux densities - ergo they are expected to be largely unlensed, luminous, rare and very distant. The addition of ground-based continuum photometry at longer wavelengths from the JCMT and APEX allows us to identify the dust peak in their SEDs, better constraining their redshifts. We select the SED templates best able to determine photometric redshifts using a sample of 69 high-redshift, lensed DSFGs, then perform checks to assess the impact of the CMB on our technique, and to quantify the systematic uncertainty associated with our photometric redshifts, $\sigma=0.14\,(1+z)$, using a sample of 25 galaxies with spectroscopic redshifts, each consistent with our color selection. For Herschel-selected ultrared galaxies with typical colors of $S_{500}/S_{250}\sim 2.2$ and $S_{500}/S_{350}\sim 1.3$ and flux densities, $S_{500}\sim 50\,$mJy, we determine a median redshift, $\hat{z}_{\rm phot}=3.66$, an interquartile redshift range, 3.30$-$4.27, with a median rest-frame 8$-$1000-$\mu$m luminosity, $\hat{L}_{\rm IR}$, of $1.3\times 10^{13}\,$L$_\odot$. A third lie at $z>4$, suggesting a space density, $\rho_{z>4}$, of $\approx 6 \times 10^{-7}\,$Mpc$^{-3}$. Our sample contains the most luminous known star-forming galaxies, and the most over-dense cluster of starbursting proto-ellipticals yet found., Comment: 22 pages; in press at The Astrophysical Journal; full-resolution PDF version available at http://www.roe.ac.uk/~rji/ivison2016.pdf

(abridged) We report rest-frame submillimeter H2O emission line observations of 11 HyLIRGs/ULIRGs at z~2-4 selected among the brightest lensed galaxies discovered in the Herschel-ATLAS. Using the IRAM NOEMA, we have detected 14 new H2O emission lines. The apparent luminosities of the H2O emission lines are $\mu L_{\rm{H_2O}} \sim 6-21 \times 10^8 L_\odot$, with velocity-integrated line fluxes ranging from 4-15 Jy km s$^{-1}$. We have also observed CO emission lines using EMIR on the IRAM 30m telescope in seven sources. The velocity widths for CO and H2O lines are found to be similar. With almost comparable integrated flux densities to those of the high-J CO line, H2O is found to be among the strongest molecular emitters in high-z Hy/ULIRGs. We also confirm our previously found correlation between luminosity of H2O ($L_{\rm{H_2O}}$) and infrared ($L_{\rm{IR}}$) that $L_{\rm{H_2O}} \sim L_{\rm{IR}}^{1.1-1.2}$, with our new detections. This correlation could be explained by a dominant role of far-infrared (FIR) pumping in the H2O excitation. Modelling reveals the FIR radiation fields have warm dust temperature $T_\rm{warm}$~45-75 K, H2O column density per unit velocity interval $N_{\rm{H_2O}}/\Delta V \gtrsim 0.3 \times 10^{15}$ cm$^{-2}$ km$^{-1}$ s and 100 $\mu$m continuum opacity $\tau_{100} > 1$ (optically thick), indicating that H2O is likely to trace highly obscured warm dense gas. However, further observations of $J\geq4$ H2O lines are needed to better constrain the continuum optical depth and other physical conditions of the molecular gas and dust. We have also detected H2O+ emission in three sources. A tight correlation between $L_{\rm{H_2O}}$ and $L_{\rm{H_2O^+}}$ has been found in galaxies from low to high redshift. The velocity-integrated flux density ratio between H2O+ and H2O suggests that cosmic rays generated by strong star formation are possibly driving the H2O+ formation., Comment: Minor changes to match the published version in A&A

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., Comment: Accepted to ApJ. 26 pages, 8 figures, 5 tables

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 ApJ

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., Comment: 27 pages, 13 figures, 7 tables, ApJ accepted

We discuss the restframe UV emission from the starbursting galaxy HFLS3 at z=6.34, discovered in Herschel/SPIRE data due to its red color in the submm wavelengths from 250-500 um. The apparent inst. SFR of HFLS3 inferred from the total FIR luminosity measured with over 15 photometric data points between 100 to 1000 um is 2900 Msun/yr. Keck/NIRC2 Ks band adaptive optics imaging data showed two potential NIR 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 HFLS3's NIR counterpart. New HST/WFC3 and ACS imaging data show both optically bright galaxies are in the foreground at z<6. A new lensing model based on HST data and mm-wave continuum emission yields a magnification of 2.2+/-0.3. The lack of multiple imaging constrains the lensing magnification to be lower than either 2.7 or 3.5 at the 95% confidence level for the two scenarios, which attribute one or two components to HFLS3 in the source plane. Correcting for gravitational lensing, the inst. SFR is 1320 Msun/yr with the 95% confidence lower limit around 830 Msun/yr. Using models for the restframe UV to FIR SED, the ave. SFR over the last 100 Myr is around 660 Msun/yr. The dust and stellar masses of HFLS3 from the same SED models are 3x10^8 Msun and ~5x10^10 Msun, respectively, with large systematic uncertainties on assumptions related to the SED model. With HST/WFC3 images we also find diffuse NIR emission about 0.5" (~3 kpc) SW of HFLS3 that remains undetected in the ACS data. The emission has a photometric redshift consistent with either z~6 or a dusty galaxy template at z~2. If at the same redshift as HFLS3 the detected diffuse emission could be part of the complex merger system that could be triggering the starburst. Alternatively, it could be part of the foreground structure at z~2.1 that is responsible for lensing of HFLS3., Comment: Submitted to ApJ, 10 pages, 6 figures

We report a highly significant ($>10\sigma$) spatial correlation between galaxies with $S_{350\mu\rm m}\ge 30\,$mJy detected in the equatorial fields of the \textsl{Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS) with estimated redshifts $\gtrsim 1.5$, and SDSS or GAMA galaxies at $0.2\le z\le 0.6$. The significance of the cross-correlation is much higher than those reported so far for samples with non-overlapping redshift distributions selected in other wavebands. Extensive, realistic simulations of clustered sub-mm galaxies amplified by foreground structures confirm that the cross-correlation is explained by weak gravitational lensing ($\mu<2$). The simulations also show that the measured amplitude and range of angular scales of the signal are larger than can be accounted for by galaxy-galaxy weak lensing. However, for scales $\lesssim 2\,$arcmin, the signal can be reproduced if SDSS/GAMA galaxies act as signposts of galaxy groups/clusters with halo masses in the range $10^{13.2}$--$10^{14.5} M_{\odot}$. The signal detected on larger scales appears to reflect the clustering of such halos., Comment: 13 pages, 14 figures, accepted for publication in MNRAS

We report on deep near-infrared observations obtained with the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST) of the first five confirmed gravitational lensing events discovered by the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We succeed in disentangling the background galaxy from the lens to gain separate photometry of the two components. The HST data allow us to significantly improve on previous constraints of the mass in stars of the lensed galaxy and to perform accurate lens modelling of these systems, as described in the accompanying paper by Dye et al. We fit the spectral energy distributions of the background sources from near-IR to millimetre wavelengths and use the magnification factors estimated by Dye et al. to derive the intrinsic properties of the lensed galaxies. We find these galaxies to have star-formation rates of approximately 400 to 2000 M_sol/yr, with approximately (6-25)x10^10 M_sol of their baryonic mass already turned into stars. At these rates of star formation, all remaining molecular gas will be exhausted in less than 100 Myr, reaching a final mass in stars of a few 10^11 M_sol. These galaxies are thus proto-ellipticals caught during their major episode of star formation, and observed at the peak epoch z=1.5-3 of the cosmic star formation history of the Universe., Comment: 14 pages, submitted to MNRAS

We have determined the mass-density radial profiles of the first five strong gravitational lens systems discovered by the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We present an enhancement of the semi-linear lens inversion method of Warren & Dye which allows simultaneous reconstruction of several different wavebands and apply this to dual-band imaging of the lenses acquired with the Hubble Space Telescope. The five systems analysed here have lens redshifts which span a range, 0.22

We present Submillimeter Array (SMA) [CII] 158um and Jansky Very Large Array (JVLA) $^{12}$CO(1-0) line emission maps for the bright, lensed, submillimeter source at $z=5.2430$ behind Abell 773: HLSJ091828.6+514223 (HLS0918). We combine these measurements with previously reported line profiles, including multiple $^{12}$CO rotational transitions, [CI], water and [NII], providing some of the best constraints on the properties of the interstellar medium (ISM) in a galaxy at $z>5$. HLS0918 has a total far-infrared (FIR) luminosity L_FIR(8-1000um) = (1.6$\pm$0.1)x10^14 L_sun/mu, where the total magnification mu_total = 8.9$\pm$1.9, via a new lens model from the [CII] and continuum maps. Despite a HyLIRG luminosity, the FIR continuum shape resembles that of a local LIRG. We simultaneously fit all of the observed spectral line profiles, finding four components which correspond cleanly to discrete spatial structures identified in the maps. The two most redshifted spectral components occupy the nucleus of a massive galaxy, with a source plane separation <1 kpc. The reddest dominates the continuum map (de-magnified L_FIR = (1.1$\pm$0.2)x10^13 L_sun), and excites strong water emission in both nuclear components via a powerful FIR radiation field from the intense star formation. A third star-forming component is most likely a region of a merging companion (dV ~ 500 km/s) exhibiting generally similar gas properties. The bluest component originates from a spatially distinct region, and photo-dissociation region (PDR) analysis suggests that it is lower density, cooler and forming stars less vigorously than the other components. Strikingly, it has very strong [NII] emission which may suggest an ionized, molecular outflow. This comprehensive view of gas properties and morphology in HLS0918 previews the science possible for a large sample of high-redshift galaxies once ALMA attains full sensitivity., Comment: 28 pages, 9 figures, accepted in ApJ

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., Comment: 28 pages, 7 tables, 9 figures. Accepted to ApJ; revised following minor comments from the referee and the community at large

We report the redshift of HATLAS J132427.0+284452 (hereafter HATLAS J132427), a gravitationally lensed starburst galaxy, the first determined 'blind' by the Herschel Space Observatory. This is achieved via the detection of [C II] consistent with z = 1.68 in a far-infrared spectrum taken with the SPIRE Fourier Transform Spectrometer. We demonstrate that the [C II] redshift is secure via detections of CO J = 2 - 1 and 3 - 2 using the Combined Array for Research in Millimeter-wave Astronomy and the Institut de Radioastronomie Millimetrique's Plateau de Bure Interferometer. The intrinsic properties appear typical of high-redshift starbursts despite the high lensing-amplified fluxes, proving the ability of the FTS to probe this population with the aid of lensing. The blind detection of [C II] demonstrates the potential of the SAFARI imaging spectrometer, proposed for the much more sensitive SPICA mission, to determine redshifts of multiple dusty galaxies simultaneously without the benefit of lensing., Comment: 6 pages, 5 figures, accepted for publication in MNRAS as a Letter

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 version

Panchromatic observations of the best candidate HyLIRG from the widest Herschel extragalactic imaging survey have led to the discovery of at least four intrinsically luminous z=2.41 galaxies across a ~100-kpc region - a cluster of starbursting proto-ellipticals. Via sub-arcsecond interferometric imaging we have measured accurate gas and star-formation surface densities. The two brightest galaxies span ~3 kpc FWHM in submm/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 M(dyn) of several x 10^11 Msun, 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 ~>5x 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 this system is the progenitor of a B(inary)-type ~10^14.6-Msun cluster., Comment: 15 pages, 8 colour figures; accepted for publication in ApJ

We present high-spatial resolution imaging obtained with the Submillimeter Array (SMA) at 880um and the Keck Adaptive Optics (AO) system at Ks-band of a gravitationally lensed sub-millimeter galaxy (SMG) at z=4.243 discovered in the Herschel-Astrophysical Terahertz Large Area Survey. The SMA data (angular resolution ~0.6") resolve the dust emission into multiple lensed images, while the Keck AO Ks-band data (angular resolution ~0.1") resolve the lens into a pair of galaxies separated by 0.3". We present an optical spectrum of the foreground lens obtained with the Gemini-South telescope that provides a lens redshift of z_lens = 0.595 +/- 0.005. We develop and apply a new lens modeling technique in the visibility plane that shows that the SMG is magnified by a factor of mu = 4.1 +/- 0.2 and has an intrinsic infrared (IR) luminosity of L_IR = (2.1 +/- 0.2) x 10^13 Lsun. We measure a half-light radius of the background source of r_s = 4.4 +/- 0.5 kpc which implies an IR luminosity surface density of Sigma_IR = (3.4 +/- 0.9) x 10^11 Lsun kpc^-2, a value that is typical of z > 2 SMGs but significantly lower than IR luminous galaxies at z~0. The two lens galaxies are compact (r_lens ~ 0.9 kpc) early-types with Einstein radii of theta_E1 = 0.57 +/- 0.01 and theta_E2 = 0.40 +/- 0.01 that imply masses of M_lens1 = (7.4 +/- 0.5) x 10^10 Msun and M_lens2 = (3.7 +/- 0.3) x 10^10 Msun. The two lensing galaxies are likely about to undergo a dissipationless merger, and the mass and size of the resultant system should be similar to other early-type galaxies at z~0.6. This work highlights the importance of high spatial resolution imaging in developing models of strongly lensed galaxies discovered by Herschel., Comment: Accepted to the ApJ. 12 pages, 6 figures, 3 tables

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane., Comment: ApJ in press. 31 pages, 16 figures, 5 tables. This version updated to match accepted version

We report measurements of the carbon monoxide ground state rotational transition (12C16O J = 1--0) with the Zpectrometer ultra-wideband spectrometer on the 100-m diameter Green Bank Telescope. The sample comprises 11 galaxies with redshifts between z = 2.1 and 3.5 from a total sample of 24 targets identified by Herschel-ATLAS photometric colors from the SPIRE instrument. Nine of the CO measurements are new redshift determinations, substantially adding to the number of detections of galaxies with rest-frame peak submillimeter emission near 100um. The CO detections confirm the existence of massive gas reservoirs within these luminous dusty star-forming galaxies (DSFGs). The CO redshift distribution of the 350um-selected galaxies is strikingly similar to the optical redshifts of 850um-selected submillimeter galaxies (SMGs) in 2.1 < z < 3.5. Spectroscopic redshifts break a temperature-redshift degeneracy; optically thin dust models fit to the far-infrared photometry indicate characteristic dust temperatures near 34 K for most of the galaxies we detect in CO. Detections of two warmer galaxies and statistically significant nondetections hint at warmer or molecule-poor DSFGs with redshifts difficult determine from from Herschel-SPIRE photometric colors alone. Many of the galaxies identified by H-ATLAS photometry are expected to be amplified by foreground gravitational lenses. Analysis of CO linewidths and luminosities provides a method for finding approximate gravitational lens magnifications mu from spectroscopic data alone, yielding mu ~ 3--20. Corrected for magnification, most galaxy luminosities are consistent with an ultra-luminous infrared galaxy (ULIRG) classification, but three are candidate hyper-LIRGs with luminosities greater than 10^13 L_sun., Comment: ApJ, accepted. 15 pages, 7 figures, 3 tables

Dust-obscured galaxies (DOGs) are a subset of high-redshift (z \approx 2) optically-faint ultra-luminous infrared galaxies (ULIRGs, e.g. L_{IR} > 10^{12} Lsun). We present new far-infrared photometry, at 250, 350, and 500 um (observed-frame), from the Herschel Space Telescope for a large sample of 113 DOGs with spectroscopically measured redshifts. Approximately 60% of the sample are detected in the far-IR, confirming their high IR luminosities, which range from 10^{11.6} Lsun < L_{IR} (8-1000 um) <10^{13.6} Lsun. 90% of the Herschel detected DOGs in this sample are ULIRGs and 30% have L_{IR} > 10^{13} Lsun. The rest-frame near-IR (1 - 3 um) SEDs of the Herschel detected DOGs are predictors of their SEDs at longer wavelengths. DOGs with "power-law" SEDs in the rest-frame near-IR show observed-frame 250/24 um flux density ratios similar to the QSO-like local ULIRG, Mrk 231. DOGs with a stellar "bump" in their rest-frame near-IR show observed-frame 250/24 um flux density ratios similar to local star-bursting ULIRGs like NGC 6240. For the Herschel detected DOGs, accurate estimates (within \approx 25%) of total IR luminosity can be predicted from their rest-frame mid-IR data alone (e.g. from Spitzer observed-frame 24 um luminosities). Herschel detected DOGs tend to have a high ratio of infrared luminosity to rest-frame 8 um luminosity (the IR8= L_{IR}(8-1000 um)/v L_{v}(8 um) parameter of Elbaz et al. 2011). Instead of lying on the z=1-2 "infrared main-sequence" of star forming galaxies (like typical LIRGs and ULIRGs at those epochs) the DOGs, especially large fractions of the bump sources, tend to lie in the starburst sequence. While, Herschel detected DOGs are similar to scaled up versions of local ULIRGs in terms of 250/24 um flux density ratio, and IR8, they tend to have cooler far-IR dust temperatures (20-40 K for DOGs vs. 40-50 K for local ULIRGs). Abridged., Comment: 24 pages, 14 figures, 3 tables, accepted for publication in the Astronomical Journal

We present high-resolution maps of stars, dust, and molecular gas in a strongly lensed submillimeter galaxy (SMG) at z = 3.259. HATLAS12--00 is selected from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) as a strong lens candidate mainly based on its unusually high 500um flux density (~300 mJy). It is the only high-redshift Planck detection in the 130 deg^2 H-ATLAS Phase 1 area. Keck Adaptive Optics images reveal a quadruply imaged galaxy in the K-band while the Submillimeter Array and the Extended Very Large Array show doubly imaged 880um and CO(1-0) sources, indicating differentiated distributions of the various components in the galaxy. In the source plane, the stars reside in three major kpc-scale clumps extended over ~1.6 kpc, the dust in a compact (~1 kpc) region ~3 kpc north of the stars, and the cold molecular gas in an extended (~7 kpc) disk ~5 kpc northeast of the stars. The emission from the stars, dust, and gas are magnified by ~17, 8, and 7 times, respectively, by four lensing galaxies at z ~ 1. Intrinsically, the lensed galaxy is a warm (T_dust ~ 40-65 K), hyper-luminous (L_IR ~ 1.7e13 Lsun; SFR ~ 2000 Msun/yr), gas-rich (M_gas/M_baryon ~ 70%), young (M_stellar/SFR ~ 20 Myr), and short-lived (M_gas/SFR ~ 40 Myr) starburst, without a significant active galactic nucleus. With physical properties similar to unlensed z > 2 SMGs, HATLAS12--00 offers a detailed view of a typical SMG through a powerful cosmic microscope., Comment: ApJ accepted version. Minor revisions. 12 pages, 4 figures, emulateapj style

We use spitzer-IRAC data to identify near-infrared counterparts to submillimeter galaxies detected with Herschel-SPIRE at 250um in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). Using a likelihood ratio analysis we identify 146 reliable IRAC counterparts to 123 SPIRE sources out of the 159. We find that, compared to the field population, the SPIRE counterparts occupy a distinct region of 3.6 and 4.5um color-magnitude space, and we use this property to identify a further 23 counterparts to 13 SPIRE sources. The IRAC identification rate of 86% is significantly higher than those that have been demonstrated with wide-field ground-based optical and near-IR imaging of Herschel fields. We estimate a false identification rate of 3.6%, corresponding to 4 to 5 sources. Among the 73 counterparts that are undetected in SDSS, 57 have both 3.6 and 4.5um coverage. Of these 43 have [3.6] - [4.5]> 0 indicating that they are likely to be at z > 1.4. Thus, ~ 40% of identified SPIRE galaxies are likely to be high redshift (z > 1.4) sources. We discuss the statistical properties of the IRAC-identified SPIRE galaxy sample including far-IR luminosities, dust temperatures, star-formation rates, and stellar masses. The majority of our detected galaxies have 10^10 to 10^11 L_sun total IR luminosities and are not intense starbursting galaxies as those found at z ~ 2, but they have a factor of 2 to 3 above average specific star-formation rates compared to near-IR selected galaxy samples., Comment: 16 pages, 12 figures, two tables. Accepted for publication in ApJ

The Spitzer Space Telescope has identified a population of ultra-luminous infrared galaxies (ULIRGs) at z ~ 2 that may play an important role in the evolution of massive galaxies. We measure the stellar masses of two populations of Spitzer-selected ULIRGs, both of which have extremely red R-[24] colors (dust-obscured galaxies, or DOGs) and compare our results with sub-millimeter selected galaxies (SMGs). One set of 39 DOGs has a local maximum in their mid-IR spectral energy distribution (SED) at rest-frame 1.6um associated with stellar emission ("bump DOGs"), while the other set of 51 DOGs has a power-law dominated mid-IR SED with spectral features typical of obscured AGN ("power-law DOGs"). We use stellar population synthesis models applied self-consistently to broad-band photometry in the rest-frame ultra-violet, optical, and near-infrared of each of these populations and test a variety of stellar population synthesis codes, star-formation histories (SFHs), and initial mass functions (IMFs). Assuming a simple stellar population SFH and a Chabrier IMF, we find that the median and inner quartile stellar masses of SMGs, bump DOGs and power-law DOGs are given by log(M_*/M_sun) = 10.42_-0.36^+0.42, 10.62_-0.32^+0.36, and 10.71_-0.34^+0.40, respectively. Implementing more complicated SFHs with multiple age components increases these mass estimates by up to 0.5 dex. Our stellar mass estimates are consistent with physical mechanisms for the origin of z~2 ULIRGs that result in high star-formation rates for a given stellar mass. Such mechanisms are usually driven by a major merger of two gas-rich systems, rather than smooth accretion of gas and small satellites., Comment: 17 pages, 7 figures, 3 tables. Plus figures showing SEDs and best-fit synthesized stellar population model. Accepted to the ApJ

We present Hubble Space Telescope (HST) imaging of 22 ultra-luminous infrared galaxies (ULIRGs) at z~2 with extremely red R-[24] colors (called dust-obscured galaxies, or DOGs) which have a local maximum in their spectral energy distribution (SED) at rest-frame 1.6um associated with stellar emission. These sources, which we call "bump DOGs", have star-formation rates of 400-4000 Msun/yr and have redshifts derived from mid-IR spectra which show strong polycyclic aromatic hydrocarbon emission --- a sign of vigorous on-going star-formation. Using a uniform morphological analysis, we look for quantifiable differences between bump DOGs, power-law DOGs (Spitzer-selected ULIRGs with mid-IR SEDs dominated by a power-law and spectral features that are more typical of obscured active galactic nuclei than starbursts), sub-millimeter selected galaxies (SMGs), and other less-reddened ULIRGs from the Spitzer extragalactic First Look Survey (XFLS). Bump DOGs are larger than power-law DOGs (median Petrosian radius of 8.4 +/- 2.7 kpc vs. 5.5 +/- 2.3 kpc) and exhibit more diffuse and irregular morphologies (median M_20 of -1.08 +/- 0.05 vs. -1.48 +/- 0.05). These trends are qualitatively consistent with expectations from simulations of major mergers in which merging systems during the peak star-formation rate period evolve from M_20 = -1.0 to M_20 = -1.7. Less obscured ULIRGs (i.e., non-DOGs) tend to have more regular, centrally peaked, single-object morphologies rather than diffuse and irregular morphologies. This distinction in morphologies may imply that less obscured ULIRGs sample the merger near the end of the peak star-formation rate period. Alternatively, it may indicate that the intense star-formation in these less-obscured ULIRGs is not the result of a recent major merger., Comment: Accepted to ApJ; 22 pages, 8 Figures, 7 Tables

We have obtained high spatial resolution Keck OSIRIS integral field spectroscopy of four z~1.5 ultra-luminous infrared galaxies that exhibit broad H-alpha emission lines indicative of strong AGN activity. The observations were made with the Keck laser guide star adaptive optics system giving a spatial resolution of 0.1", or <1 kpc at these redshifts. These high spatial resolution observations help to spatially separate the extended narrow-line regions --- possibly powered by star formation --- from the nuclear regions, which may be powered by both star formation and AGN activity. There is no evidence for extended, rotating gas disks in these four galaxies. Assuming dust correction factors as high as A(H-alpha)=4.8 mag, the observations suggest lower limits on the black hole masses of (1 - 9) x 10^8 solar masses, and star formation rates <100 solar masses per year. The black hole masses and star formation rates of the sample galaxies appear low in comparison to other high-z galaxies with similar host luminosities. We explore possible explanations for these observations including, host galaxy fading, black hole growth, and the shut down of star formation., Comment: Accepted for publication in the Astronomical Journal. 12 pages, 6 figures, 5 tables

We have identified a sample of cool field brown dwarf candidates using IRAC data from the Spitzer Deep, Wide-Field Survey (SDWFS). The candidates were selected from 400,000 SDWFS sources with [4.5] <= 18.5 mag and required to have [3.6]-[4.5] >= 1.5 and [4.5] - [8.0] <= 2.0 on the Vega system. The first color requirement selects objects redder than all but a handful of presently known brown dwarfs with spectral classes later than T7, while the second eliminates 14 probable reddened AGN. Optical detection of 4 of the remaining 18 sources implies they are likely also AGN, leaving 14 brown dwarf candidates. For two of the brightest candidates (SDWFS J143524.44+335334.6 and SDWFS J143222.82+323746.5), the spectral energy distributions including near-infrared detections suggest a spectral class of ~ T8. The proper motion is < 0.25 "/yr, consistent with expectations for a luminosity inferred distance of >70 pc. The reddest brown dwarf candidate (SDWFS J143356.62+351849.2) has [3.6] - [4.5]=2.24 and H - [4.5] > 5.7, redder than any published brown dwarf in these colors, and may be the first example of the elusive Y-dwarf spectral class. Models from Burrows et al. (2003) predict larger numbers of cool brown dwarfs should be found for a Chabrier (2003) mass function. Suppressing the model [4.5] flux by a factor of two, as indicated by previous work, brings the Burrows models and observations into reasonable agreement. The recently launched Wide-field Infrared Survey Explorer (WISE) will probe a volume ~40x larger and should find hundreds of brown dwarfs cooler than T7., Comment: 13 pages, 6 figures, accepted for publication in the June 2010 issue of The Astronomical Journal

We present a detailed analysis of the relation between infrared luminosity and molecular line luminosity, for a variety of molecular transitions, using a sample of 34 nearby galaxies spanning a broad range of infrared luminosities (10^{10} < L_{IR} < 10^{12.5} L_sun). We show that the power-law index of the relation is sensitive to the critical density of the molecular gas tracer used, and that the dominant driver in observed molecular line ratios in galaxies is the gas density. As most nearby ultraluminous infrared galaxies (ULIRGs) exhibit strong signatures of active galactic nuclei (AGN) in their center, we revisit previous claims questioning the reliability of HCN as a probe of the dense gas responsible for star formation in the presence of AGN. We find that the enhanced HCN(1-0)/CO(1-0) luminosity ratio observed in ULIRGs can be successfully reproduced using numerical models with fixed chemical abundances and without AGN-induced chemistry effects. We extend this analysis to a total of ten molecular line ratios by combining the following transitions: CO(1-0), HCO+(1-0), HCO+(3-2), HCN(1-0), and HCN(3-2). Our results suggest that AGNs reside in systems with higher dense gas fraction, and that chemistry or other effects associated with their hard radiation field may not dominate (NGC 1068 is one exception). Galaxy merger could be the underlying cause of increased dense gas fraction and the evolutionary stage of such mergers may be another determinant of the HCN/CO luminosity ratio., Comment: 20 pages, 10 figures. Published in The Astrophysical Journal (ApJ 707, 1217). Fixed a reference

We present SHARC-II 350um imaging of twelve 24um-bright (F_24um > 0.8 mJy) Dust-Obscured Galaxies (DOGs) and CARMA 1mm imaging of a subset of 2 DOGs, all selected from the Bootes field of the NOAO Deep Wide-Field Survey. Detections of 4 DOGs at 350um imply IR luminosities which are consistent within a factor of 2 of expectations based on a warm dust spectral energy distribution (SED) scaled to the observed 24um flux density. The 350um upper limits for the 8 non-detected DOGs are consistent with both Mrk231 and M82 (warm dust SEDs), but exclude cold dust (Arp220) SEDs. The two DOGs targeted at 1mm were not detected in our CARMA observations, placing strong constraints on the dust temperature: T_dust > 35-60 K. Assuming these dust properties apply to the entire sample, we find dust masses of ~3x10^8 M_sun. In comparison to other dusty z ~ 2 galaxy populations such as sub-millimeter galaxies (SMGs) and other Spitzer-selected high-redshift sources, this sample of DOGs has higher IR luminosities (2x10^13 L_sun vs. 6x10^12 L_sun for the other galaxy populations), warmer dust temperatures (>35-60 K vs. ~30 K), and lower inferred dust masses (3x10^8 M_sun vs. 3x10^9 M_sun). Herschel and SCUBA-2 surveys should be able to detect hundreds of these power-law dominated DOGs. We use HST and Spitzer/IRAC data to estimate stellar masses of these sources and find that the stellar to gas mass ratio may be higher in our 24um-bright sample of DOGs than in SMGs and other Spitzer-selected sources. Although larger sample sizes are needed to provide a definitive conclusion, the data are consistent with an evolutionary trend in which the formation of massive galaxies at z~2 involves a sub-millimeter bright, cold-dust and star-formation dominated phase followed by a 24um-bright, warm-dust and AGN-dominated phase., Comment: 16 pages, 7 figures, 6 tables; accepted to the ApJ

We present high spatial resolution optical and near-infrared imaging obtained using the ACS, WFPC2 and NICMOS cameras aboard the Hubble Space Telescope of 31 24um--bright z~2 Dust Obscured Galaxies (DOGs) identified in the Bootes Field of the NOAO Deep Wide-Field Survey. Although this subset of DOGs have mid-IR spectral energy distributions dominated by a power-law component suggestive of an AGN, all but one of the galaxies are spatially extended and not dominated by an unresolved component at rest-frame UV or optical wavelengths. The observed V-H and I-H colors of the extended components are 0.2-3 magnitudes redder than normal star-forming galaxies. All but 1 have axial ratios >0.3, making it unlikely that DOGs are composed of an edge-on star-forming disk. We model the spatially extended component of the surface brightness distributions of the DOGs with a Sersic profile and find effective radii of 1-6 kpc. This sample of DOGs is smaller than most sub-millimeter galaxies (SMGs), but larger than quiescent high-redshift galaxies. Non-parametric measures (Gini and M20) of DOG morphologies suggest that these galaxies are more dynamically relaxed than local ULIRGs. We estimate lower limits to the stellar masses of DOGs based on the rest-frame optical photometry and find that these range from ~10^(9-11) M_sun. If major mergers are the progenitors of DOGs, then these observations suggest that DOGs may represent a post-merger evolutionary stage., Comment: 23 pages, 9 figures, 6 tables, accepted to ApJ; lower limits on stellar mass revised upwards by factor of (1+z)

We present observations made with the 10m Heinrich Hertz Sub-Millimeter Telescope of HCN (3-2) emission from a sample of 30 nearby galaxies ranging in infrared luminosity from 10^10 - 10^12.5 L_sun and HCN (3-2) luminosity from 10^6 - 10^9 K km s^-1 pc^2. We examine the correlation between the infrared luminosity and HCN (3-2) luminosity and find that the best fit linear regression has a slope (in log-log space) of 0.74+/-0.12. Including recently published data from Gracia-Carpio et al. tightens the constraints on the best-fit slope to 0.79+/-0.09. This slope below unity suggests that the HCN (3-2) molecular line luminosity is not linearly tracing the amount of dense gas. Our results are consistent with predictions from recent theoretical models that find slopes below unity when the line luminosity depends upon the average gas density with a power-law index greater than a Kennicutt-Schmidt index of 1.5., Comment: 4 pages, 1 table, 2 figures; accepted to the ApJL

We present a 5'x5' integrated intensity map of 12CO (J=3-2) emission from the rho-Ophiuchi cloud core that traces low-luminosity outflow emission from two protostars: Elias 29 and, most likely, LFAM 26. The morphology of the outflow from Elias 29 is bipolar and has a curved axis that traces the S-shaped symmetry seen in H_2 emission. The outflow from LFAM 26 is a new detection and oriented in the east/west direction near the plane of the sky with most of the blue-shifted emission being absorbed by intervening clouds. The outflow axis of this object also appears to intersect a knot of H_2 emission previously attributed to Elias 29. LFAM 26 is a low luminosity source (L_bol = 0.06 L_sun) which, in combination with the observed outflow, makes it a candidate Very Low Luminosity Object (VeLLO). We derive lower limits to the gas column densities and energetics for both outflows. The mechanical luminosities for Elias 29 and LFAM 26 are 6.4 and 10.3 x 10^{-3} L_sun, respectively., Comment: 9 pages, 2 figures, accepted to the ApJL

We report measurements of the millimeter wavelength brightness fluctuations produced by the atmosphere above the South Pole made with the Arcminute Cosmology Bolometer Array Receiver (ACBAR). The data span the 2002 Austral winter during which ACBAR was mounted on the Viper telescope at the South Pole. We recover the atmospheric signal in the presence of instrument noise by calculating the correlation between signals from distinct elements of the ACBAR bolometer array. With this method, it is possible to measure atmospheric brightness fluctuations with high SNR even under the most stable atmospheric conditions. The observed atmospheric signal is characterized by the parameters of the Komolgorov-Taylor (KT) model, which are the amplitude and power law exponent describing the atmospheric power spectrum, and the two components of the wind angular velocity at the time of the observation. The KT model is typically a good description of the observed fluctuations, and fits to the data produce values of the Komolgorov exponent that are consistent with theoretical expectations. By combining the wind angular velocity results with measurements of the wind linear velocity, we find that the altitude of the observed atmospheric fluctuations is consistent with the distribution of water vapor determined from radiosonde data. For data corresponding to frequency passbands centered on 150, 219, and 274 GHz, we obtain median fluctuation power amplitudes of [10, 38, 74] mK^{2} rad^{-5/3} in Rayleigh-Jeans temperature units. Comparing with previous work, we find that these median amplitudes are approximately an order of magnitude smaller than those found at the South Pole during the Austral summer and at least 30 times lower than found at the ALMA site in the Atacama desert., Comment: 13 pages, 15 figures, submitted to ApJ, vertical margins fixed