R. S. Bussmann

Gravitational Lens Models Based on Submillimeter Array Imaging of Herschel-selected Strongly Lensed Sub-millimeter Galaxies at z > 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 (i.e., multiply imaged), 4 are moderately lensed (i.e., singly imaged), and the remainder require additional data to determine their lensing status. We apply a visibility-plane lens modeling technique to the SMA data to recover information about the masses of the lenses as well as the intrinsic (i.e., unlensed) sizes (r half) and far-infrared luminosities (L FIR) of the lensed submillimeter galaxies (SMGs). The sample of lenses comprises primarily isolated massive galaxies, but includes some groups and clusters as well. Several of the lenses are located at z lens > 0.7, a redshift regime that is inaccessible to lens searches based on Sloan Digital Sky Survey spectroscopy. The lensed SMGs are amplified by factors that are significantly below statistical model predictions given the 500 μm flux densities of our sample. We speculate that this may reflect a deficiency in our understanding of the intrinsic sizes and luminosities of the brightest SMGs. The lensed SMGs span nearly one decade in L FIR (median L FIR = 7.9 × 1012 L ☉) and two decades in FIR luminosity surface density (median ΣFIR = 6.0 × 1011 L ☉ kpc–2). The strong lenses in this sample and others identified via (sub-)mm surveys will provide a wealth of information regarding the astrophysics of galaxy formation and evolution over a wide range in redshift.

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 M_(dyn) of several × 10^(11) M ☉, 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 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^(14.6)-M ☉ cluster.

Blind detections of CO J = 1-0 in 11 H-ATLAS galaxies at z = 2.1-3.5 with the GBT/Zpectrometer

We report measurements of the carbon monoxide ground state rotational transition (^(12)C^(16)O J = 1-0) with the Zpectrometer ultrawideband 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 100 μm. The CO detections confirm the existence of massive gas reservoirs within these luminous dusty star-forming galaxies (DSFGs). The CO redshift distribution of the 350 μm selected galaxies is strikingly similar to the optical redshifts of 850 μm-selected submillimeter galaxies 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 that are difficult to determine 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 μ from spectroscopic data alone, yielding μ ~ 3-20. Corrected for magnification, most galaxy luminosities are consistent with an ultraluminous infrared galaxy classification, but three are candidate hyper-LIRGs with luminosities greater than 10^(13) L_☉.

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

Stellar archaeology shows that massive elliptical galaxies formed rapidly about ten billion years ago with star-formation rates of above several hundred solar masses per year. Their progenitors are probably the submillimetre bright galaxies at redshifts z greater than 2. Although the mean molecular gas mass (5 × 1010 solar masses) of the submillimetre bright galaxies can explain the formation of typical elliptical galaxies, it is inadequate to form elliptical galaxies that already have stellar masses above 2 × 1011 solar masses at z ≈ 2. Here we report multi-wavelength high-resolution observations of a rare merger of two massive submillimetre bright galaxies at z = 2.3. The system is seen to be forming stars at a rate of 2,000 solar masses per year. The star-formation efficiency is an order of magnitude greater than that of normal galaxies, so the gas reservoir will be exhausted and star formation will be quenched in only around 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 about 4 × 1011 solar masses. We conclude that gas-rich major galaxy mergers with intense star formation can form the most massive elliptical galaxies by z ≈ 1.5.

A comprehensive view of a strongly lensed Planck-associated submillimeter galaxy

We present high-resolution maps of stars, dust, and molecular gas in a strongly lensed submillimeter galaxy (SMG) at z = 3.259. HATLAS J114637.9–001132 is selected from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) as a strong lens candidate mainly based on its unusually high 500 μm flux density (~300 mJy). It is the only high-redshift Planck detection in the 130 deg^2 H-ATLAS Phase-I area. Keck Adaptive Optics images reveal a quadruply imaged galaxy in the K band while the Submillimeter Array and the Jansky Very Large Array show doubly imaged 880 μm 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 emissions 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.7 × 10^(13) L_☉; star formation rate (SFR) ~2000 M_☉ yr^(–1)), gas-rich (M_(gas)/M_(baryon) ~ 70%), young (M_(stellar)/SFR ~ 20 Myr), and short-lived (M_(gas)/SFR ~ 40 Myr) starburst. With physical properties similar to unlensed z > 2 SMGs, HATLAS J114637.9–001132 offers a detailed view of a typical SMG through a powerful cosmic microscope.

[C II] AND 12CO(1-0) EMISSION MAPS IN HLSJ091828.6+514223: A STRONGLY LENSED INTERACTING SYSTEM AT z = 5.24

We present Submillimeter Array [C II] 158 μm and Karl G. Jansky Very Large Array ^(12)CO(1-0) line emission maps for the bright, lensed, submillimeter source at z = 5.2430 behind A 773: HLSJ091828.6+514223 (HLS0918). We combine these measurements with previously reported line profiles, including multiple ^(12)CO rotational transitions, [C I], water, and [N II], providing some of the best constraints on the properties of the interstellar medium in a galaxy at z > 5. HLS0918 has a total far-infrared (FIR) luminosity L _(FIR(8–1000 μm)) = (1.6 ± 0.1) × 10^(14) L_☉ μ^(–1), where the total magnification μ_(total) = 8.9 ± 1.9, via a new lens model from the [C II] 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 that 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 (demagnified L_(FIR, component) = (1.1 ± 0.2) × 10^(13) L_☉) 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 (ΔV ~ 500 km s^(–1)) exhibiting generally similar gas properties. The bluest component originates from a spatially distinct region and photodissociation region analysis suggests that it is lower density, cooler, and forming stars less vigorously than the other components. Strikingly, it has very strong [N II] 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.

Herschel-ATLAS: modelling the first strong gravitational lenses

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. We present an enhancement of the semilinear 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 ≤ z ≤ 0.94. Our findings are consistent with other studies by concluding that: (1) the logarithmic slope of the total mass density profile steepens with decreasing redshift; (2) the slope is positively correlated with the average total projected mass density of the lens contained within half the effective radius and negatively correlated with the effective radius; (3) the fraction of dark matter contained within half the effective radius increases with increasing effective radius and increases with redshift.

H2O emission in high-z ultra-luminous infrared galaxies

Using the IRAM Plateau de Bure interferometer (PdBI), we report the detection of water vapor in six new lensed ultra-luminous starburst galaxies at high redshift, discovered in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). The sources are detected either in the 202−111 or 211−202 H2O emission lines with integrated line fluxes ranging from 1.8 to 14 Jy km s-1. The corresponding apparent luminosities are μLH2O ~ 3−12 × 108 L⊙, where μ is the lensing magnification factor (3  100 mJy) discovered in the Herschel surveys. Correcting the luminosities for amplification, using existing lensing models, LH2O is found to have a strong dependence on the infrared luminosity, varying as ~LIR1.2. This relation, which needs to be confirmed with better statistics, may indicate a role of radiative (infrared) excitation of the H2O lines, and implies that high-z galaxies with LIR ≳ 1013 L⊙ tend to be very strong emitters in water vapor, that have no equivalent in the local universe.

Spitzer-IRAC identification of Herschel-Atlas SPIRE sources

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.

Herschel -ATLAS/GAMA: SDSS cross-correlation induced by weak lensing

We report a highly significant (>10σ) spatial correlation between galaxies with S350 μm ≥ 30 mJy detected in the equatorial fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) with estimated redshifts ≳ 1.5, and Sloan Digital Sky Survey (SDSS) or Galaxy And Mass Assembly (GAMA) galaxies at 0.2 ≤ z ≤ 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 can be explained by weak gravitational lensing (μ < 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 ≲ 2 arcmin, the signal can be reproduced if SDSS/GAMA galaxies act as signposts of galaxy groups/clusters with halo masses in the range 1013.2–1014.5 M⊙. The signal detected on larger scales appears to reflect the clustering of such haloes.