Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty star-forming galaxies. However, the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 570 deg2 of the extragalactic sky, 4 times larger than all the other Herschel extragalactic surveys combined, in five far-infrared and submillimeter bands. We describe the survey, the complementary multiwavelength data sets that will be combined with the Herschel data, and the six major science programs we are undertaking. Using new models based on a previous submillimeter survey of galaxies, we present predictions of the properties of the ATLAS sources in other wave bands

We present the first direct and unbiased measurement of the evolution of the dust mass function of galaxies over the past 5 billion years of cosmic history using data from the Science Demonstration Phase of the Herschel‐Astrophysical Terahertz Large Area Survey (Herschel‐ATLAS). The sample consists of galaxies selected at 250 m which have reliable counterparts from the Sloan Digital Sky Survey (SDSS) at z < 0.5, and contains 1867 sources. Dust masses are calculated using both a single‐temperature grey‐body model for the spectral energy distribution and also a model with multiple temperature components. The dust temperature for either model shows no trend with redshift. Splitting the sample into bins of redshift reveals a strong evolution in the dust properties of the most massive galaxies. At z= 0.4–0.5, massive galaxies had dust masses about five times larger than in the local Universe. At the same time, the dust‐to‐stellar mass ratio was about three to four times larger, and the optical depth derived from fitting the UV‐sub‐mm data with an energy balance model was also higher. This increase in the dust content of massive galaxies at high redshift is difficult to explain using standard dust evolution models and requires a rapid gas consumption time‐scale together with either a more top‐heavy initial mass function (IMF), efficient mantle growth, less dust destruction or combinations of all three. This evolution in dust mass is likely to be associated with a change in overall interstellar medium mass, and points to an enhanced supply of fuel for star formation at earlier cosmic epochs.

We propose a novel method to constrain turbulence and bulk motions in massive galaxies, galaxy groups and clusters, exploring both simulations and observations. As emerged in the recent picture of the top-down multiphase condensation, the hot gaseous halos are tightly linked to all other phases in terms of cospatiality and thermodynamics. While hot halos (∼ 10 7 K) are perturbed by subsonic turbulence, warm (∼ 10 4 K) ionized and neutral filaments condense out of the turbulent eddies. The peaks condense into cold molecular clouds (< 100 K) raining in the core via chaotic cold accretion (CCA). We show all phases are tightly linked in terms of the ensemble (wide-aperture) velocity dispersion along the line of sight. The correlation arises in complementary long-term AGN feedback simulations and high-resolution CCA runs, and is corroborated by the combined Hitomi and new Integral Field Unit measurements in Perseus cluster. The ensemble multiphase gas distributions (from UV to radio band) are characterized by substantial spectral line broadening (σ v,los ≈ 100 -200 km s −1 ) with mild line shift. On the other hand, pencil-beam detections (as HI absorption against the AGN backlight) sample the small-scale clouds displaying smaller broadening and significant line shift up to several 100 km s −1 (for those falling toward the AGN), with increased scatter due to the turbulence intermittency. We present new ensemble σ v,los of the warm Hα+[NII] gas in 72 observed cluster/group cores: the constraints are consistent with the simulations and can be used as robust proxies for the turbulent velocities, in particular for the challenging hot plasma (otherwise requiring extremely long X-ray exposures). Finally, we show the physically motivated criterion C ≡ t cool /t eddy ≈ 1 best traces the condensation extent region and presence of multiphase gas in observed clusters and groups. The ensemble method can be applied to many available spectroscopic datasets and can substantially advance our understanding of multiphase halos in light of the next-generation multiwavelength missions.

Exploiting the H-ATLAS Science Demonstration Phase (SDP) survey data, we have determined the luminosity functions (LFs) at rest-frame wavelengths of 100 and 250 µm and at several redshifts z 1, for bright sub-mm galaxies with star formation rates (SFR) 100 M ⊙ yr −1 . We find that the evolution of the comoving LF is strong up to z ≈ 2.5, and slows down at higher redshifts. From the LFs and the information on halo masses inferred from clustering analysis, we derived an average relation between SFR and halo mass (and its scatter). We also infer that the timescale of the main episode of dust-enshrouded star formation in massive halos (M H 3 × 10 12 M ⊙ ) amounts to ∼ 7 × 10 8 yr. Given the SFRs, which are in the range 10 2 − 10 3 M ⊙ yr −1 , this timescale implies final stellar masses of order of 10 11 − 10 12 M ⊙ . The corresponding stellar mass function matches the observed mass function of passively evolving galaxies at z 1. The comparison of the statistics for sub-mm and UV selected galaxies suggests that the dust-free, UV bright phase, is 10 2 times shorter than the sub-mm bright phase, implying that the dust must form soon after the onset of star formation. Using a single reference Spectral Energy Distribution (SED; the one of the z ≈ 2.3 galaxy SMM J2135-0102), our simple physical model is able to reproduce not only the LFs at different redshifts > 1 but also the counts at wavelengths ranging from 250 µm to ≈ 1 mm. Owing to the steepness of the counts and their relatively broad frequency range, this result suggests that the dispersion of sub-mm SEDs of z > 1 galaxies around the reference one is rather small.

We report measurements of the carbon monoxide ground state rotational transition ( 12 C 16 O 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 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 (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 µ from spectroscopic data alone, yielding µ ∼ 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 ⊙ .

An ALMA observation of the early-type galaxy NGC 5044, which resides at the center of an X-ray bright group with a moderate cooling flow, detected 24 molecular structures within the central 2.5 kpc. The masses of the molecular structures vary from 3 × 10 5 M to 10 7 M and the CO(2-1) linewidths vary from 15 to 65 km s −1 . Given the large CO(2-1) linewidths, the observed structures are likely giant molecular associations (GMAs) and not individual molecular clouds (GMCs). Only a few of the GMAs are spatially resolved and the average density of these GMAs yields a GMC volume filling factor of about 15%. The masses of the resolved GMAs are insufficient for them to be gravitationally bound, however, the most massive GMA does contain a less massive component with a linewidth of 5.5 km s −1 (typical of an individual virialized GMC). We also show that the GMAs cannot be pressure confined by the hot gas. Given the CO(2-1) linewidths of the GMAs (i.e., the velocity dispersion of the embedded GMCs) they should disperse on a timescale of about 12 Myr. No disk-like molecular structures are detected and all indications suggest that the molecular gas follows ballistic trajectories after condensing out of the thermally unstable hot gas. The 230 GHz luminosity of the central continuum source is 500 times greater than its low frequency radio luminosity and probably reflects a recent accretion event. The spectrum of the central continuum source also exhibits an absorption feature with a linewidth typical of an individual GMC and an infalling velocity of 250 km s −1 .

A fundamental gap in the current understanding of galaxies concerns the thermodynamical evolution of ordinary, baryonic matter. On the one hand, radiative emission drastically decreases the thermal energy content of the interstellar plasma (ISM), inducing a slow cooling flow towards the centre. On the other hand, the active galactic nucleus (AGN) struggles to prevent the runaway cooling catastrophe, injecting huge amount of energy into the ISM. The present study intends to investigate thoroughly the role of mechanical AGN feedback in (isolated or massive) elliptical galaxies, extending and completing the mass range of tested cosmic environments. Our previously successful feedback models in galaxy clusters and groups demonstrated that AGN outflows, self‐regulated by cold gas accretion, are able to quench the cooling flow properly without destroying the cool core. Via three‐dimensional hydrodynamic simulations (flash 3.3), also including stellar evolution, we show that massive mechanical AGN outflows can indeed solve the cooling‐flow problem for the entire life of the galaxy, at the same time reproducing typical observational features and constraints such as buoyant underdense bubbles, elliptical shock cocoons, sonic ripples, dredge‐up of metals, subsonic turbulence and extended filamentary or nuclear cold gas. In order to avoid overheating and totally emptying the isolated galaxy, the frequent mechanical AGN feedback should be less powerful and efficient (ε∼ 10−4) compared with the heating required for more massive and bound ellipticals surrounded by the intragroup medium (ε∼ 10−3).