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 present the first major data release of the largest single key-project in area carried out in open time with the Herschel Space Observatory. The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) is a survey of 600 deg 2 in five photometric bands -100, 160, 250, 350 and 500 µm -with the PACS and SPIRE cameras. In this paper and a companion paper (Bourne et al. 2016) we present the survey of three fields on the celestial equator, covering a total area of 161.6 deg 2 and previously observed in the Galaxy and Mass Assembly (GAMA) spectroscopic survey. This paper describes the Herschel images and catalogues of the sources detected on the SPIRE 250 µm images. The 1σ noise for source detection, including both confusion and instrumental noise, is 7.4, 9.4 and 10.2 mJy at 250, 350 and 500 µm. Our catalogue includes 120230 sources in total, with 113995, 46209 and 11011 sources detected at > 4σ at 250, 350 and 500 µm. The catalogue contains detections at > 3σ at 100 and 160 µm for 4650 and 5685 sources, and the typical noise at these wavelengths is 44 and 49 mJy. We include estimates of the completeness of the survey and of the effects of flux bias and also describe a novel method for determining the true source counts. The H-ATLAS source counts are very similar to the source counts from the deeper HerMES survey at 250 and 350 µm, with a small difference at 500 µm. Appendix A provides a quick start in using the released datasets, including instructions and cautions on how to use them.

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 × 10 12 L ) and two decades in FIR luminosity surface density (median Σ FIR = 6.0 × 10 11 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.

Use policyThe full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-prot purposes provided that:• a full bibliographic reference is made to the original source • a link is made to the metadata record in DRO • the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders.Please consult the full DRO policy for further details.

Radio emission is a key indicator of star formation activity in galaxies, but the radio luminositystar formation relation has to date been studied almost exclusively at frequencies of 1.4 GHz or above. At lower radio frequencies, the effects of thermal radio emission are greatly reduced, and so we would expect the radio emission observed to be completely dominated by synchrotron radiation from supernova-generated cosmic rays. As part of the LOFAR Surveys Key Science project, the Herschel-ATLAS NGP field has been surveyed with LOFAR at an effective frequency of 150 MHz. We select a sample from the MPA-JHU catalogue of Sloan Digital Sky Survey galaxies in this area: the combination of Herschel, optical and mid-infrared data enable us to derive star formation rates (SFRs) for our sources using spectral energy distribution fitting, allowing a detailed study of the low-frequency radio luminosity-star formation relation in the nearby Universe. For those objects selected as star-forming galaxies (SFGs) using optical emission line diagnostics, we find a tight relationship between the 150 MHz radio luminosity (L 150 ) and SFR. Interestingly, we find that a single power-law relationship between L 150 and SFR is not a good description of all SFGs: a broken power-law model provides a better fit. This may indicate an additional mechanism for the generation of radio-emitting cosmic rays. Also, at given SFR, the radio luminosity depends on the stellar mass of the galaxy. Objects that were not classified as SFGs have higher 150-MHz radio luminosity than would be expected given their SFR, implying an important role for low-level active galactic nucleus activity.

We use the energy-balance code MAGPHYS to determine stellar and dust masses, and dust corrected star-formation rates for over 200,000 GAMA galaxies, 170,000 G10-COSMOS galaxies and 200,000 3D-HST galaxies. Our values agree well with previously reported measurements and constitute a representative and homogeneous dataset spanning a broad range in stellar mass (10 8 -10 12 M ), dust mass (10 6 -10 9 M ), and star-formation rates (0.01-100M yr −1 ), and over a broad redshift range (0.0 < z < 5.0). We combine these data to measure the cosmic star-formation history (CSFH), the stellar-mass density (SMD), and the dust-mass density (DMD) over a 12 Gyr timeline. The data mostly agree with previous estimates, where they exist, and provide a quasi-homogeneous dataset using consistent mass and star-formation estimators with consistent underlying assumptions over the full time range. As a consequence our formal errors are significantly reduced when compared to the historic literature. Integrating our cosmic star-formation history we precisely reproduce the stellar-mass density with an ISM replenishment factor of 0.50 ± 0.07, consistent with our choice of Chabrier IMF plus some modest amount of stripped stellar mass. Exploring the cosmic dust density evolution, we find a gradual increase in dust density with lookback time. We build a simple phenomenological model from the CSFH to account for the dust mass evolution, and infer two key conclusions: (1) For every unit of stellar mass which is formed 0.0065-0.004 units of dust mass is also formed; (2) Over the history of the Universe approximately 90 to 95 per cent of all dust formed has been destroyed and/or ejected.

Context. The Frontier Fields survey is a pioneering observational program aimed at collecting photometric data, both from space (Hubble Space Telescope and Spitzer Space Telescope) and from ground-based facilities (VLT Hawk-I), for six deep fields pointing at clusters of galaxies and six nearby deep parallel fields, in a wide range of passbands. The analysis of these data is a natural outcome of the Astrodeep project, an EU collaboration aimed at developing methods and tools for extragalactic photometry and creating valuable public photometric catalogues. Aims. We produce multiwavelength photometric catalogues (from B to 4.5 µm) for the first two of the Frontier Fields, Abell-2744 and MACS-J0416 (plus their parallel fields). Methods. To detect faint sources even in the central regions of the clusters, we develop a robust and repeatable procedure that uses the public codes Galapagos and Galfit to model and remove most of the light contribution from both the brightest cluster members, and the intra-cluster light. We perform the detection on the processed HST H160 image to obtain a pure H-selected sample, which is the primary catalogue that we publish. We also add a sample of sources which are undetected in the H160 image but appear on a stacked infrared image. Photometry on the other HST bands is obtained using SExtractor, again on processed images after the procedure for foreground light removal. Photometry on the Hawk-I and IRAC bands is obtained using our PSF-matching deconfusion code t-phot. A similar procedure, but without the need for the foreground light removal, is adopted for the Parallel fields.Results. The procedure of foreground light subtraction allows for the detection and the photometric measurements of ∼2500 sources per field. We deliver and release complete photometric H-detected catalogues, with the addition of the complementary sample of infrared-detected sources. All objects have multiwavelength coverage including B to H HST bands, plus K-band from Hawk-I, and 3.6−4.5 µm from Spitzer. full and detailed treatment of photometric errors is included. We perform basic sanity checks on the reliability of our results. Conclusions. The multiwavelength photometric catalogues are available publicly and are ready to be used for scientific purposes. Our procedures allows for the detection of outshone objects near the bright galaxies, which, coupled with the magnification effect of the clusters, can reveal extremely faint high redshift sources. Full analysis on photometric redshifts is presented in Paper II.

No abstract