Bright Submillimeter Galaxies Research Articles

The gravitationally unstable gas disk of a starburst galaxy 12 billion years ago.

Galaxies in the early Universe that are bright at submillimetre wavelengths (submillimetre-bright galaxies) are forming stars at a rate roughly 1,000 times higher than the Milky Way. A large fraction of the new stars form in the central kiloparsec of the galaxy1-3, a region that is comparable in size to the massive, quiescent galaxies found at the peak of cosmic star-formation history4 and the cores of present-day giant elliptical galaxies. The physical and kinematic properties inside these compact starburst cores are poorly understood because probing them at relevant spatial scales requires extremely high angular resolution. Here we report observations with a linear resolution of 550 parsecs of gas and dust in an unlensed, submillimetre-bright galaxy at a redshift of z = 4.3, when the Universe was less than two billion years old. We resolve the spatial and kinematic structure of the molecular gas inside the heavily dust-obscured core and show that the underlying gas disk is clumpy and rotationally supported (that is, its rotation velocity is larger than the velocity dispersion). Our analysis of the molecular gas mass per unit area suggests that the starburst disk is gravitationally unstable, which implies that the self-gravity of the gas is stronger than the differential rotation of the disk and the internal pressure due to stellar-radiation feedback. As a result of the gravitational instability in the disk, the molecular gas would be consumed by star formation on a timescale of 100 million years, which is comparable to gas depletion times in merging starburst galaxies5.

The impact of clustering and angular resolution on far-infrared and millimeter continuum observations

Follow-up observations at high-angular resolution of bright submillimeter galaxies selected from deep extragalactic surveys have shown that the single-dish sources are comprised of a blend of several galaxies. Consequently, number counts derived from low- and high-angular-resolution observations are in tension. This demonstrates the importance of resolution effects at these wavelengths and the need for realistic simulations to explore them. We built a new 2 deg2 simulation of the extragalactic sky from the far-infrared to the submillimeter. It is based on an updated version of the 2SFM (two star-formation modes) galaxy evolution model. Using global galaxy properties generated by this model, we used an abundance-matching technique to populate a dark-matter lightcone and thus simulate the clustering. We produced maps from this simulation and extracted the sources, and we show that the limited angular resolution of single-dish instruments has a strong impact on (sub)millimeter continuum observations. Taking into account these resolution effects, we are reproducing a large set of observables, as number counts and their evolution with redshift and cosmic infrared background power spectra. Our simulation consistently describes the number counts from single-dish telescopes and interferometers. In particular, at 350 and 500 μm, we find that the number counts measured by Herschel between 5 and 50 mJy are biased towards high values by a factor ~2, and that the redshift distributions are biased towards low redshifts. We also show that the clustering has an important impact on the Herschel pixel histogram used to derive number counts from P(D) analysis. We find that the brightest galaxy in the beam of a 500 μm Herschel source contributes on average to only ~60% of the Herschel flux density, but that this number will rise to ~95% for future millimeter surveys on 30 m-class telescopes (e.g., NIKA2 at IRAM). Finally, we show that the large number density of red Herschel sources found in observations but not in models might be an observational artifact caused by the combination of noise, resolution effects, and the steepness of color- and flux density distributions. Our simulation, called Simulated Infrared Dusty Extragalactic Sky (SIDES), is publicly available.

The Circumgalactic Medium of Submillimeter Galaxies. II. Unobscured QSOs within Dusty Starbursts and QSO Sightlines with Impact Parameters below 100 kpc

Abstract We present Atacama Large Millimeter/submillimeter Array (ALMA) 870 μm observations of 29 bright Herschel sources near high-redshift QSOs. The observations confirm that 20 of the Herschel sources are submillimeter-bright galaxies (SMGs) and identify 16 new SMG−QSO pairs that are useful to studies of the circumgalactic medium (CGM) of SMGs. Eight out of the 20 SMGs are blends of multiple 870 μm sources. The angular separations for six of the Herschel-QSO pairs are less than 10″, comparable to the sizes of the Herschel beam and the ALMA primary beam. We find that four of these six “pairs” are actually QSOs hosted by SMGs. No additional submillimeter companions are detected around these QSOs, and the rest-frame ultraviolet spectra of the QSOs show no evidence of significant reddening. Black hole accretion and star formation contribute almost equally in bolometric luminosity in these galaxies. The SMGs hosting QSOs show similar source sizes, dust surface densities, and star formation rate surface densities to those of other SMGs in the sample. We find that the black holes are growing ∼3× faster than the galaxies when compared to the present-day black hole/galaxy mass ratio, suggesting a QSO duty cycle of ≲30% in SMGs at . The remaining two Herschel-detected QSOs are undetected at 870 μm, but each has an SMG “companion” only 9″ and 12″ away (71 and 95 kpc at z = 3). They could be either merging or projected pairs. If the former, they would represent a rare class of “wet−dry” mergers. If the latter, the QSOs would, for the first time, probe the CGM of SMGs at impact parameters below 100 kpc.

Blending bias impacts the host halo masses derived from a cross-correlation analysis of bright submillimetre galaxies

Placing bright sub-millimetre galaxies (SMGs) within the broader context of galaxy formation and evolution requires accurate measurements of their clustering, which can constrain the masses of their host dark matter halos. Recent work has shown that the clustering measurements of these galaxies may be affected by a `blending bias,' which results in the angular correlation function of the sources extracted from single-dish imaging surveys being boosted relative to that of the underlying galaxies. This is due to confusion introduced by the coarse angular resolution of the single-dish telescope and could lead to the inferred halo masses being significantly overestimated. We investigate the extent to which this bias affects the measurement of the correlation function of SMGs when it is derived via a cross-correlation with a more abundant galaxy population. We find that the blending bias is essentially the same as in the auto-correlation case and conclude that the best way to reduce its effects is to calculate the angular correlation function using SMGs in narrow redshift bins. Blending bias causes the inferred host halo masses of the SMGs to be overestimated by a factor of $\sim6$ when a redshift interval of $\delta z=3$ is used. However, this reduces to a factor of $\sim2$ for $\delta z=0.5$. The broadening of photometric redshift probability distributions with increasing redshift can therefore impart a mild halo `downsizing' effect onto the inferred host halo masses, though this trend is not as strong as seen in recent observational studies.

VLA AND ALMA IMAGING OF INTENSE GALAXY-WIDE STAR FORMATION IN z ∼ 2 GALAXIES

ABSTRACT We present ≃0.″4 resolution extinction-independent distributions of star formation and dust in 11 star-forming galaxies (SFGs) at z = 1.3–3.0. These galaxies are selected from sensitive blank-field surveys of the 2′ × 2′ Hubble Ultra-Deep Field at λ = 5 cm and 1.3 mm using the Karl G. Jansky Very Large Array and Atacama Large Millimeter/submillimeter Array. They have star formation rates (SFRs), stellar masses, and dust properties representative of massive main-sequence SFGs at z ∼ 2. Morphological classification performed on spatially resolved stellar mass maps indicates a mixture of disk and morphologically disturbed systems; half of the sample harbor X-ray active galactic nuclei (AGNs), thereby representing a diversity of z ∼ 2 SFGs undergoing vigorous mass assembly. We find that their intense star formation most frequently occurs at the location of stellar-mass concentration and extends over an area comparable to their stellar-mass distribution, with a median diameter of 4.2 ± 1.8 kpc. This provides direct evidence of galaxy-wide star formation in distant blank-field-selected main-sequence SFGs. The typical galactic-average SFR surface density is 2.5 M ⊙ yr−1 kpc−2, sufficiently high to drive outflows. In X-ray-selected AGN where radio emission is enhanced over the level associated with star formation, the radio excess pinpoints the AGNs, which are found to be cospatial with star formation. The median extinction-independent size of main-sequence SFGs is two times larger than those of bright submillimeter galaxies, whose SFRs are 3–8 times larger, providing a constraint on the characteristic SFR (∼300 M ⊙ yr−1) above which a significant population of more compact SFGs appears to emerge.

The discovery of gas-rich, dusty starbursts in luminous reddened quasars atz∼ 2.5 with ALMA

We present ALMA observations of cold dust and molecular gas in four high-luminosity, heavily reddened (A$_{\rm{V}} \sim 2.5-6$ mag) Type 1 quasars at $z\sim2.5$ with virial M$_{\rm{BH}} \sim 10^{10}$M$_\odot$, to test whether dusty, massive quasars represent the evolutionary link between submillimetre bright galaxies (SMGs) and unobscured quasars. All four quasars are detected in both the dust continuum and in the $^{12}$CO(3-2) line. The mean dust mass is 6$\times$10$^{8}$M$_\odot$ assuming a typical high redshift quasar spectral energy distribution (T=41K, $\beta$=1.95 or T=47K, $\beta$=1.6). The implied star formation rates are very high - $\gtrsim$1000 M$_\odot$ yr$^{-1}$ in all cases. Gas masses estimated from the CO line luminosities cover $\sim$1-5$\times10^{10}$($\alpha_{\rm{CO}} / 0.8$)M$_\odot$ and the gas depletion timescales are very short - $\sim5-20$Myr. A range of gas-to-dust ratios is observed in the sample. We resolve the molecular gas in one quasar - ULASJ2315$+$0143 ($z=2.561$) - which shows a strong velocity gradient over $\sim$20 kpc. The velocity field is consistent with a rotationally supported gas disk but other scenarios, e.g. mergers, cannot be ruled out at the current resolution of these data. In another quasar - ULASJ1234+0907 ($z=2.503$) - we detected molecular line emission from two millimetre bright galaxies within 200 kpc of the quasar, suggesting that this quasar resides in a significant over-density. The high detection rate of both cold dust and molecular gas in these sources, suggests that reddened quasars could correspond to an early phase in massive galaxy formation associated with large gas reservoirs and significant star formation.

FAINT SUBMILLIMETER GALAXIES IDENTIFIED THROUGH THEIR OPTICAL/NEAR-INFRARED COLORS. I. SPATIAL CLUSTERING AND HALO MASSES

ABSTRACT The properties of submillimeter galaxies (SMGs) that are fainter than the confusion limit of blank-field single-dish surveys ( 2 mJy) are poorly constrained. Using a newly developed color selection technique, Optical-Infrared Triple Color (OIRTC), that has been shown to successfully select such faint SMGs, we identify a sample of 2938 OIRTC-selected galaxies, dubbed Triple Color Galaxies (TCGs), in the UKIDSS-UDS field. We show that these galaxies have a median 850 μm flux of mJy (equivalent to a star formation rate SFR yr−1 based on spectral energy distribution fitting), representing the first large sample of faint SMGs that bridges the gap between bright SMGs and normal star-forming galaxies in S 850 and L IR. We assess the basic properties of TCGs and their relationship with other galaxy populations at . We measure the two-point autocorrelation function for this population and derive a typical halo mass of log10( ) , , and at , 2–3, and 3–5, respectively. Together with the bright SMGs ( mJy) and a comparison sample of less far-infrared luminous star-forming galaxies, we find a lack of dependence between spatial clustering and S 850 (or SFR), suggesting that the difference between these populations may lie in their local galactic environment. Lastly, on the scale of at we find a tentative enhancement of the clustering of TCGs over the comparison star-forming galaxies, suggesting that some faint SMGs are physically associated pairs, perhaps reflecting a merging origin in their triggering.

SXDF-ALMA 2 arcmin2 deep survey: Resolving and characterizing the infrared extragalactic background light down to 0.5 mJy

Abstract We present a multiwavelength analysis of five submillimeter sources (S1.1mm = 0.54–2.02 mJy) that were detected during our 1.1 mm deep continuum survey in the Subaru/XMM-Newton Deep Survey Field (SXDF)-UDS-CANDELS field (2 arcmin2, 1σ = 0.055 mJy beam−1) using the Atacama Large Millimeter/submillimeter Array (ALMA). The two brightest sources correspond to a known single-dish (AzTEC) selected bright submillimeter galaxy (SMG), whereas the remaining three are faint SMGs newly uncovered by ALMA. If we exclude the two brightest sources, the contribution of the ALMA-detected faint SMGs to the infrared extragalactic background light is estimated to be $\sim 4.1^{+5.4}_{-3.0}\:$Jy deg−2, which corresponds to $\sim 16^{+22}_{-12}\%$ of the infrared extragalactic background light. This suggests that their contribution to the infrared extragalactic background light is as large as that of bright SMGs. We identified multiwavelength counterparts of the five ALMA sources. One of the sources (SXDF-ALMA3) is extremely faint in the optical to near-infrared region despite its infrared luminosity ($L_\mathrm{IR}\simeq 1\times 10^{12}\,L_{\odot}$ or SFR ≃ 100 M⊙ yr−1). By fitting the spectral energy distributions at the optical-to-near-infrared wavelengths of the remaining four ALMA sources, we obtained the photometric redshifts (zphoto) and stellar masses (M*): zphoto ≃ 1.3–2.5, M* ≃ (3.5–9.5) × 1010 M⊙. We also derived their star formation rates (SFRs) and specific SFRs as ≃30–200 M⊙ yr−1 and ≃0.8–2 Gyr−1, respectively. These values imply that they are main sequence star-forming galaxies.

CANDIDATE GRAVITATIONALLY LENSED DUSTY STAR-FORMING GALAXIES IN THE HERSCHEL WIDE AREA SURVEYS*

ABSTRACT We present a list of candidate gravitationally lensed dusty star-forming galaxies (DSFGs) from the HerMES Large Mode Survey and the Herschel Stripe 82 Survey. Together, these partially overlapping surveys cover 372 deg2 on the sky. After removing local spiral galaxies and known radio-loud blazars, our candidate list of lensed DSFGs is composed of 77 sources with 500 μm flux densities (S 500) greater than 100 mJy. Such sources are dusty starburst galaxies similar to the first bright sub-millimeter galaxies (SMGs) discovered with SCUBA. We expect a large fraction of this list to be strongly lensed, with a small fraction made up of bright SMG–SMG mergers that appear as hyper-luminous infrared galaxies ( L IR > 10 13 L ⊙ ). Thirteen of the 77 candidates have spectroscopic redshifts from CO spectroscopy with ground-based interferometers, putting them at z > 1 and well above the redshift of the foreground lensing galaxies. The surface density of our sample is 0.21 ± 0.03 deg−2. We present follow-up imaging of a few of the candidates to confirm their lensing nature. The sample presented here is an ideal tool for higher-resolution imaging and spectroscopic observations to understand the detailed properties of starburst phenomena in distant galaxies.

Early science with the Large Millimeter Telescope: observations of extremely luminous high-zsources identified byPlanck

We present 8.5 arcsec resolution 1.1mm continuum imaging and CO spectroscopic redshift measurements of eight extremely bright submillimetre galaxies identified from the Planck and Herschel surveys, taken with the Large Millimeter Telescope's AzTEC and Redshift Search Receiver instruments. We compiled a candidate list of high redshift galaxies by cross-correlating the Planck Surveyor mission's highest frequency channel (857 GHz, FWHM = 4.5 arcmin) with the archival Herschel Spectral and Photometric Imaging Receiver (SPIRE) imaging data, and requiring the presence of a unique, single Herschel counterpart within the 150 arcsec search radius of the Planck source positions with 350 micron flux density larger than 100 mJy, excluding known blazars and foreground galaxies. All eight candidate objects observed are detected in 1.1mm continuum by AzTEC bolometer camera, and at least one CO line is detected in all cases with a spectroscopic redshift between 1.3 < z(CO) < 3.3. Their infrared spectral energy distributions mapped using the Herschel and AzTEC photometry are consistent with cold dust emission with characteristic temperature between $T_d$ = 43 K and 84 K. With apparent infrared luminosity of up to L(IR) = $3\times10^{14} \mu^{-1} L_\odot$, they are some of the most luminous galaxies ever found (with yet unknown gravitational magnification factor $\mu$). The analysis of their spectral energy distributions (SEDs) suggests that star formation is powering the bulk of their extremely large IR luminosities. Derived molecular gas masses of $M_{H2}=(0.6-7.8)\times 10^{11} M_\odot$ (for $\mu$~10) also make them some of the most gas-rich high redshift galaxies ever detected.

The formation of submillimetre-bright galaxies from gas infall over a billion years.

Submillimetre-bright galaxies at high redshift are the most luminous, heavily star-forming galaxies in the Universe and are characterized by prodigious emission in the far-infrared, with a flux of at least five millijanskys at a wavelength of 850 micrometres. They reside in haloes with masses about 10(13) times that of the Sun, have low gas fractions compared to main-sequence disks at a comparable redshift, trace complex environments and are not easily observable at optical wavelengths. Their physical origin remains unclear. Simulations have been able to form galaxies with the requisite luminosities, but have otherwise been unable to simultaneously match the stellar masses, star formation rates, gas fractions and environments. Here we report a cosmological hydrodynamic galaxy formation simulation that is able to form a submillimetre galaxy that simultaneously satisfies the broad range of observed physical constraints. We find that groups of galaxies residing in massive dark matter haloes have increasing rates of star formation that peak at collective rates of about 500-1,000 solar masses per year at redshifts of two to three, by which time the interstellar medium is sufficiently enriched with metals that the region may be observed as a submillimetre-selected system. The intense star formation rates are fuelled in part by the infall of a reservoir gas supply enabled by stellar feedback at earlier times, not through major mergers. With a lifetime of nearly a billion years, our simulations show that the submillimetre-bright phase of high-redshift galaxies is prolonged and associated with significant mass buildup in early-Universe proto-clusters, and that many submillimetre-bright galaxies are composed of numerous unresolved components (for which there is some observational evidence).

OPTICAL–INFRARED PROPERTIES OF FAINT 1.3 mm SOURCES DETECTED WITH ALMA

We report optical-infrared (IR) properties of faint 1.3 mm sources (S_1.3mm = 0.2-1.0 mJy) detected with the Atacama Large Millimeter/submillimeter Array (ALMA) in the Subaru/XMM-Newton Deep Survey (SXDS) field. We searched for optical/IR counterparts of 8 ALMA-detected sources (>=4.0 sigma, the sum of the probability of spurious source contamination is ~1) in a K-band source catalog. Four ALMA sources have K-band counterpart candidates within a 0.4" radius. Comparison between ALMA-detected and undetected K-band sources in the same observing fields shows that ALMA-detected sources tend to be brighter, more massive, and more actively forming stars. While many of the ALMA-identified submillimeter-bright galaxies (SMGs) in previous studies lie above the sequence of star-forming galaxies in stellar mass--star-formation rate plane, our ALMA sources are located in the sequence, suggesting that the ALMA-detected faint sources are more like `normal' star-forming galaxies rather than `classical' SMGs. We found a region where multiple ALMA sources and K-band sources reside in a narrow photometric redshift range (z ~ 1.3-1.6) within a radius of 5" (42 kpc if we assume z = 1.45). This is possibly a pre-merging system and we may be witnessing the early phase of formation of a massive elliptical galaxy.

A millimetre-wave redshift search for the unlensed HyLIRG, HS1700.850.1

We report the redshift of an unlensed, highly obscured submillimetre galaxy (SMG), HS1700.850.1, the brightest SMG (S850um =19.1 mJy) detected in the JCMT/SCUBA-2 Baryonic Structure Survey, based on the detection of its CO line emission. Using the IRAM PdBI-WIDEX with 3.6GHz band width, we serendipitously detect an emission line at 150.6 GHz. From a search over 14.5 GHz in the 3-mm and 2-mm atmospheric windows, we confirm the identification of this line as CO(5-4) at z = 2.816, meaning that it does not reside in the z~2.30 proto-cluster in this field. Measurement of the 870um source size (<0.85") from the Sub-Millimeter Array (SMA) confirms a compact emission in a S870um =14.5mJy, LIR~10^13 Lsun component, suggesting an Eddington-limited starburst. We use the double-peaked CO line profile measurements along with the SMA size constraints to study the gas dynamics of a HyLIRG, estimating the gas and dynamical masses of HS1700.850.1. While HS1700.850.1 is one of the most extreme galaxies known in the Universe, we find that it occupies a relative void in the Lyman-Break Galaxy distribution in this field. Comparison with other extreme objects at similar epochs (HyLIRG Quasars), and cosmological simulations, suggests such an anti-bias of bright SMGs could be relatively common, with the brightest SMGs rarely occupying the most overdense regions at z=2-4.

Early Science with the Large Millimeter Telescope: observations of dust continuum and CO emission lines of cluster-lensed submillimetre galaxies atz=2.0–4.7

We present Early Science observations with the Large Millimeter Telescope, AzTEC 1.1 mm continuum images and wide bandwidth spectra (73–111 GHz) acquired with the Redshift Search Receiver, towards four bright lensed submillimetre galaxies identified through the Herschel Lensing Survey-snapshot and the Submillimetre Common-User Bolometer Array-2 Cluster Snapshot Survey. This pilot project studies the star formation history and the physical properties of the molecular gas and dust content of the highest redshift galaxies identified through the benefits of gravitational magnification. We robustly detect dust continuum emission for the full sample and CO emission lines for three of the targets. We find that one source shows spectroscopic multiplicity and is a blend of three galaxies at different redshifts (z = 2.040, 3.252, and 4.680), reminiscent of previous high-resolution imaging follow-up of unlensed submillimetre galaxies, but with a completely different search method, that confirm recent theoretical predictions of physically unassociated blended galaxies. Identifying the detected lines as 12CO (Jup = 2–5) we derive spectroscopic redshifts, molecular gas masses, and dust masses from the continuum emission. The mean H2 gas mass of the full sample is (2.0 ± 0.2) × 1011 M⊙/μ, and the mean dust mass is (2.0 ± 0.2) × 109 M⊙/μ, where μ ≈ 2–5 is the expected lens amplification. Using these independent estimations we infer a gas-to-dust ratio of δGDR ≈ 55–75, in agreement with other measurements of submillimetre galaxies. Our magnified high-luminosity galaxies fall on the same locus as other high-redshift submillimetre galaxies, extending the |$L^{\prime }_{\rm CO}$|–LFIR correlation observed for local luminous and ultraluminous infrared galaxies to higher far-infrared and CO luminosities.

Jansky VLA S-band view of Hα emitters (HAEs) associated with a protocluster 4C23.56 at z = 2.5

AbstractWe present recent results on Karl Jansky Very Large Array (JVLA) deep S-band (2-4 GHz) observation towards a protocluster 4C23.56 at redshift z ∼ 2.5. The protocluster 4C23.56 is known to have a significant over density (∼ 5 times) of star-burst galaxies selected to be Hα line-bright by a Subaru narrow band imaging. Now we have found 25 HAEs associated with the protocluster. These starburst HAEs are likely to become massive ellipticals at z = 0 in a cluster. Various other galaxy populations also reside in this field and the fact makes the field very unique as a tool to understand galaxy formation in a over dense region. Subsequent deep 1100-μm continuum surveys by the ASTE 10-m dish have discovered that several submillimeter bright galaxies (SMGs) coincide with HAEs, suggesting HAEs undergoing dusty starbursts. As star formation rates (SFRs) of HAEs might have been underestimated, we use radio being resistant to dust extinction. We investigate the correlation between SFR1.4 GHz and SFRHα for radio index α = 0.8 to see if the correlation holds for the sources and to check the number of dusty star forming galaxies. Our final results will allow us to evaluate quantitatively how the galaxy formation channel may be different under the condition of over-densities.

ULASJ1234+0907: the reddest type 1 quasar at z = 2.5 revealed in the X-ray and far-infrared

Abstract We present Herschel and XMM–Newton observations of ULASJ1234+0907 (z = 2.503), the reddest broad-line type 1 quasar currently known with (i − K)AB &amp;gt; 7.1. Herschel observations indicate that the quasar host is a hyperluminous infrared galaxy with a total infrared luminosity of log10(LIR/L⊙) = 13.90 ± 0.02. A greybody fit gives a dust temperature of Td = 60 ± 3 K assuming an emissivity index of β = 1.5, considerably higher than in submillimeter bright galaxies observed at similar redshifts. The star formation rate is estimated to be &amp;gt;2000 M⊙ yr−1 even accounting for a significant contribution from an active galactic nucleus (AGN) component to the total infrared luminosity or requiring that only the far-infrared luminosity is powered by a starburst. XMM–Newton observations constrain the hard X-ray luminosity to be L2-10 keV = 1.3 × 1045 erg s−1, putting ULASJ1234+0907 among the brightest X-ray quasars known. Through very deep optical and near-infrared imaging of the field at subarcsecond seeing, we demonstrate that despite its extreme luminosity, it is highly unlikely that ULASJ1234+0907 is being lensed. We measure a neutral hydrogen column density of NH = 9.0 × 1021 cm−2 corresponding to AV ∼ 6. The observed properties of ULASJ1234+0907 – high luminosity and Eddington ratio, broad lines, moderate column densities and significant infrared emission from reprocessed dust – are similar to those predicted by galaxy formation simulations for the AGN blowout phase. The high Eddington ratio, combined with the presence of significant amounts of dust, is expected to drive strong outflows due to the effects of radiation pressure on dust.

HIGH-RESOLUTION SPECTROSCOPIC IMAGING OF CO IN Az= 4.05 PROTO-CLUSTER

We present a study of the formation of clustered, massive galaxies at large look-back times via spectroscopic imaging of CO in the unique GN20 proto-cluster at z = 4.05. Existing observations show that this is a dense concentration of gas-rich, very active star forming galaxies, including multiple bright submillimeter galaxies (SMGs). Using deep, high-resolution VLA CO(2-1) observations, we image the molecular gas with a resolution of ~1 kpc just 1.6 Gyr after the Big Bang. The SMGs GN20.2a and GN20.2b have deconvolved sizes of ~5 kpc X 3 kpc and ~8 kpc X 5 kpc (Gaussian FWHM) in CO(2-1), respectively, and we measure gas surface densities up to ~12,700/1,700X(sin i) (\alpha_CO/0.8) M_sun/pc^2 for GN20.2a/GN20.2b in the highest-resolution maps. Dynamical mass estimates allow us to constrain the CO-to-H_2 conversion factor to \alpha_CO = 1.7+/-0.8 M_sun (K km s^{-1} pc^2)^-1 for GN20.2a and \alpha_CO = 1.1+/-^{1.5}_{1.1} M_sun (K km/s pc^2)^-1 for GN20.2b. We measure significant offsets (0.5"-1") between the CO and optical emission, indicating either dust obscuration on scales of tens of kpc or that the emission originates from distinct galaxies. CO spectral line energy distributions imply physical conditions comparable to other SMGs and reveal further evidence that GN20.2a and GN20.2b are in different merging stages. We carry out a targeted search for CO emission from the 14 known B-band Lyman break galaxies (LBGs) in the field, tentatively detecting CO in a previously-undetected LBG and placing 3sigma upper limits on the CO luminosities of those that may lie within our bandpass. A blind search for emission line sources down to a 5sigma limiting CO luminosity of L'_CO(2-1) = 8 X 10^9 K km/s pc^2 and covering \Delta z = 0.0273 (~20 comoving Mpc) produces no other strong contenders associated with the proto-cluster.

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 × 10(10) 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 × 10(11) 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 × 10(11) 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.

Submillimetre galaxies in cosmological hydrodynamic simulations: source number counts and the spatial clustering

We use large cosmological smoothed particle hydrodynamics simulations to study the formation and evolution of submillimetre galaxies (SMGs). In our previous work, we studied the statistical properties of ultraviolet-selected star-forming galaxies at high redshifts. We populate the same cosmological simulations with SMGs by calculating the reprocessing of stellar light by dust grains into far-infrared to millimetre wavebands in a self-consistent manner. We generate light-cone outputs to compare directly the statistical properties of the simulated SMGs with available observations. Our model reproduces the submillimetre source number counts and the clustering amplitude. We show that bright SMGs with flux S > 1 mJy reside in haloes with masses of ∼1013 M⊙ and have stellar masses greater than 1011 M⊙. The angular cross-correlation between the SMGs and Lyman α emitters is significantly weaker than that between the SMGs and Lyman-break galaxies. The cross-correlation is also weaker than the autocorrelation of the SMGs. The redshift distribution of the SMGs shows a broad peak at z ∼ 2, where bright SMGs contribute significantly to the global cosmic star formation rate density. Our model predicts that there are hundreds of SMGs with S > 0.1 mJy at z > 5 per 1 deg2 field. Such SMGs can be detected by ALMA.

How to distinguish starbursts and quiescently star-forming galaxies: the ‘bimodal’ submillimetre galaxy population as a case study

In recent work (arXiv:1101.0002) we have suggested that the high-redshift (z ~ 2-4) bright submillimetre galaxy (SMG) population is heterogeneous, with major mergers contributing both at early stages, where quiescently star-forming discs are blended into one submm source (`galaxy-pair SMGs'), and late stages, where mutual tidal torques drive gas inflows and cause strong starbursts. Here we combine hydrodynamic simulations of major mergers with 3-D dust radiative transfer calculations to determine observational diagnostics that can distinguish between quiescently star-forming SMGs and starburst SMGs via integrated data alone. We fit the far-IR SEDs of the simulated galaxies with the optically thin single-temperature modified blackbody, the full form of the single-temperature modified blackbody, and a power-law temperature-distribution model. The effective dust temperature, T_dust, and power-law index of the dust emissivity in the far-IR, \beta, derived can significantly depend on the fitting form used, and the intrinsic \beta\ of the dust is not recovered. However, for all forms used here, there is a T_dust above which almost all simulated galaxies are starbursts, so a T_dust cut is very effective at selecting starbursts. Simulated merger-induced starbursts also have higher L_IR/M_gas and L_IR/L_FUV than quiescently star-forming galaxies and lie above the star formation rate-stellar mass relation. These diagnostics can be used to test our claim that the SMG population is heterogeneous and to observationally determine what star formation mode dominates a given galaxy population. We comment on applicability of these diagnostics to ULIRGs that would not be selected as SMGs. These `hot-dust ULIRGs' are typically starburst galaxies lower in mass than SMGs, but they can also simply be SMGs observed from a different viewing angle.