The morphology-density relationship states that dense cosmic environments such as galaxy clusters have an overabundance of quiescent elliptical galaxies, but it is unclear at which redshift this relationship is first established. We study the morphology of 4 clusters with $1.2<z<1.8$ using HST imaging and the morphology computation code statmorph. By comparing median morphology of cluster galaxies to CANDELS field galaxies using Monte Carlo analysis, we find that 2 out of 4 clusters (at z=1.19 and z=1.75) have an established morphology-density relationship with more than $3\sigma$ significance. $\sim$50% of galaxies in these clusters are bulge-dominated compared to $\sim$30% in the field, and they are significantly more compact. This result is more significant for low-mass galaxies with $\log M/M_\odot \lessapprox 10.5$, showing that low-mass galaxies are affected the most in clusters. We also find an intriguing system of two z $\approx$ 1.45 clusters at a unusually small separation 2D separation of $3'$ and 3D separation of $\approx73$ Mpc that exhibit no morphology-density relationship but have enhanced merger signatures. We conclude that the environmental mechanism responsible for the morphology-density relationship is 1) already active as early as z=1.75, 2) forms compact, bulge-dominated galaxies and 3) affects primarily low-mass galaxies. However, there is a significant degree of intracluster variance that may depend on the larger cosmological environment in which the cluster is embedded.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We present the star formation histories of 39 galaxies with high quality rest-frame optical spectra at 0.5<z<1.3 selected to have strong Balmer absorption lines and/or Balmer break, and compare to a sample of spectroscopically selected quiescent galaxies at the same redshift. Photometric selection identifies a majority of objects that have clear evidence for a recent short-lived burst of star formation within the last 1.5 Gyr, i.e. "post-starburst" galaxies, however we show that good quality continuum spectra are required to obtain physical parameters such as burst mass fraction and burst age. Dust attenuation appears to be the primary cause for misidentification of post-starburst galaxies, leading to contamination in spectroscopic samples where only the [OII] emission line is available, as well as a small fraction of objects lost from photometric samples. The 31 confirmed post-starburst galaxies have formed 40-90% of their stellar mass in the last 1-1.5 Gyr. We use the derived star formation histories to find that the post-starburst galaxies are visible photometrically for 0.5-1 Gyr. This allows us to update a previous analysis to suggest that 25-50% of the growth of the red sequence at z~1 could be caused by a starburst followed by rapid quenching. We use the inferred maximum historical star formation rates of several 100-1000 Msun/yr and updated visibility times to confirm that sub-mm galaxies are likely progenitors of post-starburst galaxies. The short quenching timescales of 100-200 Myr are consistent with cosmological hydrodynamic models in which rapid quenching is caused by the mechanical expulsion of gas due to an AGN.

We present a study of the dust, stars and atomic gas (HI) in an HI-selected sample of local galaxies (z<0.035) in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) fields. This HI-selected sample reveals a population of very high gas fraction (>80 per cent), low stellar mass sources that appear to be in the earliest stages of their evolution. We compare this sample with dust and stellar mass selected samples to study the dust and gas scaling relations over a wide range of gas fraction (proxy for evolutionary state of a galaxy). The most robust scaling relations for gas and dust are those linked to NUV-r (SSFR) and gas fraction, these do not depend on sample selection or environment. At the highest gas fractions, our additional sample shows the dust content is well below expectations from extrapolating scaling relations for more evolved sources, and dust is not a good tracer of the gas content. The specific dust mass for local galaxies peaks at a gas fraction of ~75 per cent. The atomic gas depletion time is also longer for high gas fraction galaxies, opposite to the trend found for molecular gas depletion timescale. We link this trend to the changing efficiency of conversion of HI to H2 as galaxies increase in stellar mass surface density as they evolve. Finally, we show that galaxies start out barely obscured and increase in obscuration as they evolve, yet there is no clear and simple link between obscuration and global galaxy properties.

An important aspect of quenching star formation is the removal of the cold interstellar medium (ISM; non-ionised gas and dust) from a galaxy. In addition, dust grains can be destroyed in a hot or turbulent medium. The adopted timescale of dust removal usually relies on uncertain theoretical estimates. It is tricky to track the dust removal, because usually dust is constantly replenished by consecutive generations of stars. Our objective is to measure observationally the timescale of dust removal. We here explore an approach to select galaxies which do have detectable amounts of dust and cold ISM but exhibit a low current dust production rate. Any decrease of the dust and gas content as a function of the age of such galaxies therefore must be attributed to processes governing the ISM removal. We used a sample of galaxies detected by Herschel in the far-infrared with visually assigned early-type morphology or spirals with red colours. We also obtained JCMT/SCUBA-2 observations for five of them. We discovered an exponential decline of the dust-to-stellar mass ratio with age, which we interpret as an evolutionary trend of dust removal from these galaxies. For the first time we directly measure the dust removal timescale in such galaxies to be tau=(2.5+-0.4) Gyr (the corresponding half-life time is (1.75+-0.25) Gyr). This quantity may be used in models in which it must be assumed a priori and cannot be derived. Any process which removes dust in these galaxies, such as dust grain destruction, cannot happen on shorter timescales. The timescale is comparable to the quenching timescales found in simulations for galaxies with similar stellar masses. The dust is likely of internal, not external origin. It was either formed in the past directly by supernovae, or from seeds produced by SNe and with grain growth in the ISM contributing substantially to the dust mass accumulation.

The Herschel-ATLAS unbiased survey of cold dust in the local Universe is dominated by a surprising population of very blue (FUV-K < 3.5), dust-rich galaxies with high gas fractions (f_HI = M_HI/(M*+M_HI)>0.5)). Dubbed `Blue and Dusty Gas Rich Sources' (BADGRS) they have cold diffuse dust temperatures, and the highest dust-to-stellar mass ratios of any galaxies in the local Universe. Here, we explore the molecular ISM in a representative sample of BADGRS, using very deep CO(J_up=1,2,3) observations across the central and outer disk regions. We find very low CO brightnesses (Tp=15-30 mK), despite the bright far-infrared emission and metallicities in the range 0.5<Z/Z_sun<1.0. The CO line ratios indicate a range of conditions with R_21=0.6-2.1 and R_31=0.2-1.2. Using a metallicity dependent conversion from CO luminosity to molecular gas mass we find M_H2/M_d=7-27 and Sigma_H2=0.5-6 M_sun pc^-2, around an order of magnitude lower than expected. The BADGRS have lower molecular gas depletion timescales (tau_d = 0.5 Gyr) than other local spirals, lying offset from the Kennicutt-Schmidt relation by a similar factor to Blue Compact Dwarf galaxies. The cold diffuse dust temperatures in BADGRS (13-16 K) require an interstellar radiation field 10-20 times lower than that inferred from their observed surface brightness. We speculate that the dust in these sources has either a very clumpy geometry or a very different opacity in order to explain the cold temperatures and lack of CO emission. BADGRS also have low UV attenuation for their UV colour suggestive of an SMC-type dust attenuation curve, different star formation histories or different dust/star geometry. They lie in a similar part of the IRX-beta space as z=5 galaxies and may be useful as local analogues for high gas fraction galaxies in the early Universe.

The origin of the star forming main sequence ( i.e., the relation between star formation rate and stellar mass, globally or on kpc-scales; hereafter SFMS) remains a hotly debated topic in galaxy evolution. Using the ALMA-MaNGA QUEnching and STar formation (ALMaQUEST) survey, we show that for star forming spaxels in the main sequence galaxies, the three local quantities, star-formation rate surface density (\sigsfr), stellar mass surface density (\sigsm), and the \h2~mass surface density (\sigh2), are strongly correlated with one another and form a 3D linear (in log) relation with dispersion. In addition to the two well known scaling relations, the resolved SFMS (\sigsfr~ vs. \sigsm) and the Schmidt-Kennicutt relation (\sigsfr~ vs. \sigh2; SK relation), there is a third scaling relation between \sigh2~ and \sigsm, which we refer to as the `molecular gas main sequence' (MGMS). The latter indicates that either the local gas mass traces the gravitational potential set by the local stellar mass or both quantities follow the underlying total mass distributions. The scatter of the resolved SFMS ($\sigma \sim 0.25$ dex) is the largest compared to those of the SK and MGMS relations ($\sigma \sim$ 0.2 dex). A Pearson correlation test also indicates that the SK and MGMS relations are more strongly correlated than the resolved SFMS. Our result suggests a scenario in which the resolved SFMS is the least physically fundamental and is the consequence of the combination of the SK and the MGMS relations.

Post-starburst galaxies, identified by their unusually strong Balmer absorption lines and weaker than average emission lines, have traditionally been selected based on their central stellar populations. Here we identify 360 galaxies with post-starburst regions from the MaNGA integral field survey and classify these galaxies into three types: 31 galaxies with central post-starburst regions (CPSB), 37 galaxies with off-center ring-like post-starburst regions (RPSB) and 292 galaxies with irregular post-starburst regions (IPSB). Focussing on the CPSB and RPSB samples, and comparing their radial gradients in D$_n$4000, H$\delta_{\rm A}$ and W(H$\alpha$) to control samples, we find that while the CPSBs have suppressed star formation throughout their bulge and disk, and clear evidence of rapid decline of star formation in the central regions, the RPSBs only show clear evidence of recently rapidly suppressed star formation in their outer regions and an ongoing central starburst. The radial profiles in mass-weighted age and stellar $v/\sigma$ indicate that CPSBs and RPSBs are not simply different evolutionary stages of the same event, rather that CPSB galaxies are caused by a significant disruptive event, while RPSB galaxies are caused by disruption of gas fuelling to the outer regions. Compared to the control samples, both CPSB and RPSB galaxies show a higher fraction of interactions/mergers, misaligned gas or bars that might be the cause of the gas inflows and subsequent quenching.

We investigate the prevalence of galactic-scale outflows in post-starburst (PSB) galaxies at high redshift ($1 < z < 1.4$), using the deep optical spectra available in the UKIDSS Ultra Deep Survey (UDS). We use a sample of $\sim40$ spectroscopically confirmed PSBs, recently identified in the UDS field, and perform a stacking analysis in order to analyse the structure of strong interstellar absorption features such as Mg ii ($\lambda2800$ Ang.). We find that for massive ($M_* > 10^{10}\rm\,M_{\odot}$) PSBs at $z > 1$, there is clear evidence for a strong blue-shifted component to the Mg ii absorption feature, indicative of high-velocity outflows ($v_{\rm out}\sim1150\pm160\rm\,km\,s^{-1}$) in the interstellar medium. We conclude that such outflows are typical in massive PSBs at this epoch, and potentially represent the residual signature of a feedback process that quenched these galaxies. Using full spectral fitting, we also obtain a typical stellar velocity dispersion $\sigma_*$ for these PSBs of $\sim200\rm\,km\,s^{-1}$, which confirms they are intrinsically massive in nature (dynamical mass $M_{\rm d}\sim10^{11}\rm\,M_{\odot}$). Given that these high-$z$ PSBs are also exceptionally compact ($r_{\rm e}\sim1$--$2\rm\,kpc$) and spheroidal (Sersic index $n\sim3$), we propose that the outflowing winds may have been launched during a recent compaction event (e.g. major merger or disc collapse) that triggered either a centralised starburst or active galactic nuclei (AGN) activity. Finally, we find no evidence for AGN signatures in the optical spectra of these PSBs, suggesting they were either quenched by stellar feedback from the starburst itself, or that if AGN feedback is responsible, the AGN episode that triggered quenching does not linger into the post-starburst phase.

A galaxy's stellar mass is one of its most fundamental properties, but it remains challenging to measure reliably. With the advent of very large optical spectroscopic surveys, efficient methods that can make use of low signal-to-noise spectra are needed. With this in mind, we created a new software package for estimating effective stellar mass-to-light ratios $\log \Upsilon^*$ that uses principal component analysis(PCA) basis set to optimize the comparison between observed spectra and a large library of stellar population synthesis models. In Paper I, we showed that a with a set of six PCA basis vectors we could faithfully represent most optical spectra from the Mapping Nearby Galaxies at APO (MaNGA) survey;and we tested the accuracy of our M/L estimates using synthetic spectra. Here, we explore sources of systematic error in our mass measurements by comparing our new measurements to data from the literature. We compare our stellar mass surface density estimates to kinematics-derived dynamical mass surface density measurements from the DiskMass Survey and find some tension between the two which could be resolved if the disk scale-heights used in the kinematic analysis were overestimated by a factor of $\sim 1.5$. We formulate an aperture-corrected stellar mass catalog for the MaNGA survey, and compare to previous stellar mass estimates based on multi-band optical photometry, finding typical discrepancies of 0.1 dex. Using the spatially resolved MaNGA data, we evaluate the impact of estimating total stellar masses from spatially unresolved spectra, and we explore how the biases that result from unresolved spectra depend upon the galaxy's dust extinction and star formation rate. Finally, we describe a SDSS Value-Added Catalog which will include both spatially resolved and total (aperture-corrected) stellar masses for MaNGA galaxies.

We present a method of fitting optical spectra of galaxies using a basis set of six vectors obtained from principal component analysis (PCA) of a library of synthetic spectra of 40000 star formation histories (SFHs). Using this library, we provide estimates of resolved effective stellar mass-to-light ratio ($\log \Upsilon^*$) for thousands of galaxies from the SDSS-IV/MaNGA integral-field spectroscopic survey. Using a testing framework built on additional synthetic SFHs, we show that the estimates of stellar mass-to-light ratio are reliable (as are their uncertainties) at a variety of signal-to-noise ratios, stellar metallicities, and dust attenuation conditions. Finally, we describe the future release of the resolved stellar mass-to-light ratios as a SDSS-IV/MaNGA Value-Added Catalog (VAC) and provide a link to the software used to conduct this analysis.

We present the HI-MaNGA programme of HI follow-up for the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. MaNGA, which is part of the Fourth phase of the Sloan Digital Sky Surveys (SDSS-IV), is in the process of obtaining integral field unit (IFU) spectroscopy for a sample of ~10,000 nearby galaxies. We give an overview of the HI 21cm radio follow-up observing plans and progress and present data for the first 331 galaxies observed in the 2016 observing season at the Robert C. Bryd Green Bank Telescope (GBT). We also provide a cross match of the current MaNGA (DR15) sample with publicly available HI data from the Arecibo Legacy Fast Arecibo L-band Feed Array (ALFALFA) survey. The addition of HI data to the MaNGA data set will strengthen the survey's ability to address several of its key science goals that relate to the gas content of galaxies, while also increasing the legacy of this survey for all extragalactic science.

We present the severe, systematic uncertainties currently facing our understanding of dust emission, which stymie our ability to truly exploit dust as a tool for studying galaxy evolution. We propose a program of study to tackle these uncertainties, describe the necessary facilities, and discuss the potential science gains that will result. This white paper was submitted to the US National Academies' Astro2020 Decadal Survey on Astronomy and Astrophysics.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

The large Integral Field Spectroscopy (IFS) surveys have allowed the classification of ionizing sources of emission lines on sub-kpc scales. In this work, we define two non-parametric parameters, quiescence (f$_{q}$) and its concentration (c$_{q}$), to quantify the strength and the spatial distribution of the quenched areas, respectively, traced by the LI(N)ER regions with low EW(H$\alpha$). With these two measurements, we classify MaNGA galaxies into inside-out and outside-in quenching types according to their locations on the f$_{q}$ vs. c$_{q}$ plane and we measure the fraction of inside-out (outside-in) quenching galaxies as a function of halo mass. We find that the fraction of galaxies showing inside-out quenching increases with halo mass, irrespective of stellar mass or galaxy type (satellites vs. centrals). In addition, high stellar mass galaxies exhibit a greater fraction of inside-out quenching compared to low stellar mass ones in all environments. In contrast, the fraction of outside-in quenching does not depend on halo mass. Our results suggest that morphological quenching may be responsible for the inside-out quenching seen in all environments. On the other hand, the flat dependence of the outside-in quenching on halo mass could be a mixed result of ram-pressure stripping and galaxy mergers. Nevertheless, at a given environment and stellar mass, the fraction of inside-out quenching is systematically greater than that of outside-in quenching, suggesting that inside-out quenching is the dominant quenching mode in all environments.

Post-starburst galaxies are typically considered to be a transition population, en route to the red sequence after a recent quenching event. Despite this, recent observations have shown that these objects typically have large reservoirs of cold molecular gas. In this paper we study the star-forming gas properties of a large sample of post-starburst galaxies selected from the cosmological, hydrodynamical EAGLE simulations. These objects resemble observed high-mass post-starburst galaxies both spectroscopically and in terms of their space density, stellar mass distribution and sizes. We find that the vast majority of simulated post-starburst galaxies have significant gas reservoirs, with star-forming gas masses of ~10$^9$ M$_{\odot}$, in good agreement with those seen in observational samples. The simulation reproduces the observed time evolution of the gas fraction of the post-starburst galaxy population, with the average galaxy losing ~90 per cent of its star-forming interstellar medium in only ~600 Myr. A variety of gas consumption/loss processes are responsible for this rapid evolution, including mergers and environmental effects, while active galactic nuclei play only a secondary role. The fast evolution in the gas fraction of post-starburst galaxies is accompanied by a clear decrease in the efficiency of star formation, due to a decrease in the dense gas fraction. We predict that forthcoming ALMA observations of the gas reservoirs of low-redshift post-starburst galaxies will show that the molecular gas is typically compact and has disturbed kinematics, reflecting the disruptive nature of many of the evolutionary pathways that build up the post-starburst galaxy population.

We use data from 1222 late-type star-forming galaxies in the SDSS IV MaNGA survey to identify regions in which the gas-phase metallicity is anomalously-low compared to expectations from the tight empirical relation between metallicity and stellar surface mass-density at a given stellar mass. We find anomalously low metallicity (ALM) gas in 10% of the star-forming spaxels, and in 25% of the galaxies in the sample. The incidence rate of ALM gas increases strongly with both global and local measures of the specific star-formation rate, and is higher in lower mass galaxies and in the outer regions of galaxies. The incidence rate is also significantly higher in morphologically disturbed galaxies. We estimate that the lifetimes of the ALM regions are a few hundred Myr. We argue that the ALM gas has been delivered to its present location by a combination of interactions, mergers, and accretion from the halo, and that this infusion of gas stimulates star-formation. Given the estimated lifetime and duty cycle of such events, we estimate that the time-averaged accretion rate of ALM gas is similar to the star-formation rate in late type galaxies over the mass-range M$_* \sim10^9$ to 10$^{10}$ M$_{\odot}$.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, one of three core programs of the fourth-generation Sloan Digital Sky Survey (SDSS-IV), is producing a massive, high-dimensional integral field spectroscopic data set. However, leveraging the MaNGA data set to address key questions about galaxy formation presents serious data-related challenges due to the combination of its spatially inter-connected measurements and sheer volume. For each galaxy, the MaNGA pipelines produce relatively large data files to preserve the spatial correlations of the spectra and measurements, but this comes at the expense of storing the data set in a coarsely-chunked manner. The coarse chunking and total volume of the data make it time-consuming to download and curate locally-stored data. Thus, accessing, querying, visually exploring, and performing statistical analyses across the whole data set at a fine-grained scale is extremely challenging using just FITS files. To overcome these challenges, we have developed \marvin: a toolkit consisting of a Python package, Application Programming Interface (API), and web application utilizing a remote database. \marvin's robust and sustainable design minimizes maintenance, while facilitating user-contributed extensions such as high level analysis code. Finally, we are in the process of abstracting out \marvin's core functionality into a separate product so that it can serve as a foundation for others to develop \marvin-like systems for new science applications.

Unresolved gas and dust observations show a surprising diversity in the amount of interstellar matter in early-type galaxies. Using ALMA observations we resolve the ISM in z$\sim$0.05 early-type galaxies. From a large sample of early-type galaxies detected in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) we selected five of the dustiest cases, with dust masses M$_d\sim$several$\times10^7$M$_\odot$, with the aim of mapping their submillimetre continuum and $^{12}$CO(2-1) line emission distributions. These observations reveal molecular gas disks. There is a lack of associated, extended continuum emission in these ALMA observations, most likely because it is resolved out or surface brightness limited, if the dust distribution is as extended as the CO gas. However, two galaxies have central continuum ALMA detections. An additional, slightly offset, continuum source is revealed in one case, which may have contributed to confusion in the Herschel fluxes. Serendipitous continuum detections further away in the ALMA field are found in another case. Large and massive rotating molecular gas disks are mapped in three of our targets, reaching a few$\times10^{9}$M$_\odot$. One of these shows evidence of kinematic deviations from a pure rotating disc. The fields of our two remaining targets contain only smaller, weak CO sources, slightly offset from the optical galaxy centres. These may be companion galaxies seen in ALMA observations, or background objects. These heterogeneous findings in a small sample of dusty early-type galaxies reveal the need for more such high spatial resolution studies, to understand statistically how dust and gas are related in early-type galaxies.

We analyse the structure of galaxies with high specific star formation rate (SSFR) in cluster and field environments in the redshift range $0.5<z<1.0$. Recent studies have shown that these galaxies are strongly depleted in dense environments due to rapid environmental quenching, giving rise to post-starburst galaxies (PSBs). We use effective radii and S\'ersic indices as tracers of galaxy structure, determined using imaging from the UKIDSS Ultra Deep Survey (UDS). We find that the high-SSFR galaxies that survive into the cluster environment have, on average, larger effective radii than those in the field. We suggest that this trend is likely to be driven by the most compact star-forming galaxies being preferentially quenched in dense environments. We also show that the PSBs in clusters have stellar masses and effective radii that are similar to the missing compact star-forming population, suggesting that these PSBs are the result of size-dependent quenching. We propose that both strong stellar feedback and the stripping of the extended halo act together to preferentially and rapidly quench the compact and low-mass star-forming systems in clusters to produce PSBs. We test this scenario using the stacked spectra of 124 high-SSFR galaxies, showing that more compact galaxies are more likely to host outflows. We conclude that a combination of environmental and secular processes is the most likely explanation for the appearance of PSBs in galaxy clusters.

Measuring molecular gas mass is vital for understanding the evolution of galaxies at high redshifts (z$\geq$1). Most measurements rely on CO as a tracer, but dependences on metallicity, dynamics and surface density lead to systematic uncertainties in high-z galaxies, where these physical properties are difficult to observe, and where the physical environments can differ systematically from those at z=0. Dust continuum emission provides a potential alternative assuming a known dust/gas ratio, but this must be calibrated on a direct gas tracer at z$\geq$1. In this paper we consider the [CI] 492-GHz emission line, which has been shown to trace molecular gas closely throughout Galactic clouds and has the advantages of being optically thin in typical conditions (unlike CO), and being observable at accessible frequencies at high redshifts (in contrast to the low-excitation lines of CO). We use the Atacama Large Millimetre/submillimetre Array (ALMA) to measure [CI], CO(4-3) and dust emission in a representative sample of star-forming galaxies at z=1, and combine these data with multi-wavelength spectral energy distributions to study relationships between dust and gas components of galaxies. We uncover a strong [CI]-dust correlation, suggesting that both trace similar phases of the gas. By incorporating other samples from the literature, we show that this correlation persists over a wide range of luminosities and redshifts up to z$\sim$4. Finally we explore the implications of our results as an independent test of literature calibrations for dust as a tracer of gas mass, and for predicting the CI abundance.

JINGLE is a new JCMT legacy survey designed to systematically study the cold interstellar medium of galaxies in the local Universe. As part of the survey we perform 850um continuum measurements with SCUBA-2 for a representative sample of 193 Herschel-selected galaxies with M*>10^9Msun, as well as integrated CO(2-1) line fluxes with RxA3m for a subset of 90 of these galaxies. The sample is selected from fields covered by the Herschel-ATLAS survey that are also targeted by the MaNGA optical integral-field spectroscopic survey. The new JCMT observations combined with the multi-wavelength ancillary data will allow for the robust characterization of the properties of dust in the nearby Universe, and the benchmarking of scaling relations between dust, gas, and global galaxy properties. In this paper we give an overview of the survey objectives and details about the sample selection and JCMT observations, present a consistent 30 band UV-to-FIR photometric catalog with derived properties, and introduce the JINGLE Main Data Release (MDR). Science highlights include the non-linearity of the relation between 850um luminosity and CO line luminosity, and the serendipitous discovery of candidate z>6 galaxies.

A key task of observational extragalactic astronomy is to determine where -- within galaxies of diverse masses and morphologies -- stellar mass growth occurs, how it depends on galaxy properties and what processes regulate star formation. Using spectroscopic indices derived from the stellar continuum at $\sim 4000$\AA, we determine the spatially resolved star-formation histories of 980000 spaxels in 2404 galaxies in the SDSS-IV MaNGA IFU survey. We examine the spatial distribution of star-forming, quiescent, green valley, starburst and post-starburst spaxels as a function of stellar mass and morphology to see where and in what types of galaxy star formation is occurring. The spatial distribution of star-formation is dependent primarily on stellar mass, with a noticeable change in the distribution at \mstar$>10^{10}$\msun. Galaxies above this mass have an increasing fraction of regions that are forming stars with increasing radius, whereas lower mass galaxies have a constant fraction of star forming regions with radius. Our findings support a picture of inside-out growth and quenching at high masses. We find that morphology (measured via concentration) correlates with the fraction of star-forming spaxels, but not with their radial distribution. We find (post-)starburst regions are more common outside of the galaxy centre, are preferentially found in asymmetric galaxies, and have lower gas-phase metallicity than other regions, consistent with interactions triggering starbursts and driving low metallicity gas into regions at $<1.5R_e$.

We present an analysis of the structure of post-starburst (PSB) galaxies in the redshift range $0.5 < z < 2$, using a photometrically-selected sample identified in the Ultra Deep Survey (UDS) field. We examine the structure of $\sim80$ of these transient galaxies using radial light $\mu(r)$ profiles obtained from CANDELS $\textit{Hubble Space Telescope}$ near-infrared/optical imaging, and compare to a large sample of $\sim2000$ passive and star-forming galaxies. For each population, we determine their typical structural properties (effective radius $r_{\rm e}$, S\'ersic index $n$) and find significant differences in PSB structure at different epochs. At high redshift ($z > 1$), PSBs are typically massive ($M_* > 10^{10}\rm\,M_{\odot}$), very compact and exhibit high S\'ersic indices, with structures that differ significantly from their star-forming progenitors but are similar to massive passive galaxies. In contrast, at lower redshift ($0.5 < z < 1$), PSBs are generally of low mass ($M_* < 10^{10}\rm\,M_{\odot}$) and exhibit compact but less concentrated profiles (i.e. lower S\'ersic indices), with structures similar to low-mass passive discs. Furthermore, for both epochs we find remarkably consistent PSB structure across the optical/near-infrared wavebands (which largely trace different stellar populations), suggesting that any preceding starburst and/or quenching in PSBs was not strongly centralized. Taken together, these results imply that PSBs at $z > 1$ have been recently quenched during a major disruptive event (e.g. merger or protogalactic collapse) which formed a compact remnant, while at $z < 1$ an alternative less disruptive process is primarily responsible. Our results suggest that high-$z$ PSBs are an intrinsically different population to those at lower redshifts, and indicate different quenching routes are active at different epochs.

Post-starburst galaxies can be identified via the presence of prominent Hydrogen Balmer absorption lines in their spectra. We present a comprehensive study of the origin of strong Balmer lines in a volume-limited sample of 189 galaxies with $0.01<z<0.05$, $\log(\mbox{M}_{\star}/\mbox{M}_{\odot})>9.5$ and projected axis ratio $b/a>0.32$. We explore their structural properties, environments, emission lines and star formation histories, and compare them to control samples of star-forming and quiescent galaxies, and simulated galaxy mergers. Excluding contaminants, in which the strong Balmer lines are most likely caused by dust-star geometry, we find evidence for three different pathways through the post-starburst phase, with most events occurring in intermediate-density environments: (1) a significant disruptive event, such as a gas-rich major merger, causing a starburst and growth of a spheroidal component, followed by quenching of the star formation (70% of post-starburst galaxies at $9.5<\log(\mbox{M}_{\star}/\mbox{M}_{\odot})<10.5$ and 60% at $\log(\mbox{M}_{\star}/\mbox{M}_{\odot})>10.5$); (2) at $9.5<\log(\mbox{M}_{\star}/\mbox{M}_{\odot})<10.5$, stochastic star formation in blue-sequence galaxies, causing a weak burst and subsequent return to the blue sequence (30%); (3) at $\log(\mbox{M}_{\star}/\mbox{M}_{\odot})>10.5$, cyclic evolution of quiescent galaxies which gradually move towards the high-mass end of the red sequence through weak starbursts, possibly as a result of a merger with a smaller gas-rich companion (40%). Our analysis suggests that AGN are `on' for $50%$ of the duration of the post-starburst phase, meaning that traditional samples of post-starburst galaxies with strict emission line cuts will be at least $50%$ incomplete due to the exclusion of narrow-line AGN.

We investigate the relationship between the quenching of star formation and the structural transformation of massive galaxies, using a large sample of photometrically-selected post-starburst galaxies in the UKIDSS UDS field. We find that post-starburst galaxies at high-redshift ($z>1$) show high S\'ersic indices, significantly higher than those of active star-forming galaxies, but with a distribution that is indistinguishable from the old quiescent population. We conclude that the morphological transformation occurs before (or during) the quenching of star formation. Recently quenched galaxies are also the most compact; we find evidence that massive post-starburst galaxies (M$_{\ast}> 10^{10.5} ~$M$_{\odot}$) at high redshift ($z>1$) are on average smaller than comparable quiescent galaxies at the same epoch. Our findings are consistent with a scenario in which massive passive galaxies are formed from three distinct phases: (1) gas-rich dissipative collapse to very high densities, forming the proto-spheroid; (2) rapid quenching of star formation, to create the "red nugget" with post-starburst features; (3) a gradual growth in size as the population ages, perhaps as a result of minor mergers.

Gas-rich minor mergers contribute significantly to the gas reservoir of early-type galaxies (ETGs) at low redshift, yet the star formation efficiency (SFE; the star formation rate divided by the molecular gas mass) appears to be strongly suppressed following some of these events, in contrast to the more well-known merger-driven starbursts. We present observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of six ETGs, which have each recently undergone a gas-rich minor merger, as evidenced by their disturbed stellar morphologies. These galaxies were selected because they exhibit extremely low SFEs. We use the resolving power of ALMA to study the morphology and kinematics of the molecular gas. The majority of our galaxies exhibit spatial and kinematical irregularities, such as detached gas clouds, warps, and other asymmetries. These asymmetries support the interpretation that the suppression of the SFE is caused by dynamical effects stabilizing the gas against gravitational collapse. Through kinematic modelling we derive high velocity dispersions and Toomre Q stability parameters for the gas, but caution that such measurements in edge-on galaxies suffer from degeneracies. We estimate merger ages to be about 100~Myr based on the observed disturbances in the gas distribution. Furthermore, we determine that these galaxies lie, on average, two orders of magnitude below the Kennicutt-Schmidt relation for star-forming galaxies as well as below the relation for relaxed ETGs. We discuss potential dynamical processes responsible for this strong suppression of star formation surface density at fixed molecular gas surface density.

One key problem in astrophysics is understanding how and why galaxies switch off their star formation, building the quiescent population that we observe in the local Universe. From the GAMA and VIPERS surveys, we use spectroscopic indices to select quiescent and candidate transition galaxies. We identify potentially rapidly transitioning post-starburst galaxies, and slower transitioning green-valley galaxies. Over the last 8 Gyrs the quiescent population has grown more slowly in number density at high masses (M$_*>10^{11}$M$_\odot$) than at intermediate masses (M$_*>10^{10.6}$M$_\odot$). There is evolution in both the post-starburst and green valley stellar mass functions, consistent with higher mass galaxies quenching at earlier cosmic times. At intermediate masses (M$_*>10^{10.6}$M$_\odot$) we find a green valley transition timescale of 2.6 Gyr. Alternatively, at $z\sim0.7$ the entire growth rate could be explained by fast-quenching post-starburst galaxies, with a visibility timescale of 0.5 Gyr. At lower redshift, the number density of post-starbursts is so low that an unphysically short visibility window would be required for them to contribute significantly to the quiescent population growth. The importance of the fast-quenching route may rapidly diminish at $z<1$. However, at high masses (M$_*>10^{11}$M$_\odot$), there is tension between the large number of candidate transition galaxies compared to the slow growth of the quiescent population. This could be resolved if not all high mass post-starburst and green-valley galaxies are transitioning from star-forming to quiescent, for example if they rejuvenate out of the quiescent population following the accretion of gas and triggering of star formation, or if they fail to completely quench their star formation.

The Herschel Space Observatory has revealed a very different galaxyscape from that shown by optical surveys which presents a challenge for galaxy-evolution models. The Herschel surveys reveal (1) that there was rapid galaxy evolution in the very recent past and (2) that galaxies lie on a a single Galaxy Sequence (GS) rather than a star-forming `main sequence' and a separate region of `passive' or `red-and-dead' galaxies. The form of the GS is now clearer because far-infrared surveys such as the Herschel ATLAS pick up a population of optically-red star-forming galaxies that would have been classified as passive using most optical criteria. The space-density of this population is at least as high as the traditional star-forming population. By stacking spectra of H-ATLAS galaxies over the redshift range 0.001 < z < 0.4, we show that the galaxies responsible for the rapid low-redshift evolution have high stellar masses, high star-formation rates but, even several billion years in the past, old stellar populations - they are thus likely to be relatively recent ancestors of early-type galaxies in the Universe today. The form of the GS is inconsistent with rapid quenching models and neither the analytic bathtub model nor the hydrodynamical EAGLE simulation can reproduce the rapid cosmic evolution. We propose a new gentler model of galaxy evolution that can explain the new Herschel results and other key properties of the galaxy population.

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 Msol), dust mass (10^6---10^9 Msol), and star-formation rates (0.01---100 Msol per yr), 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.

The process of planet formation offers a rich source of dust production via grain growth in protostellar discs, and via grinding of larger bodies in debris disc systems. Chemical evolution models, designed to follow the build up of metals and dust in galaxies, do not currently account for planet formation. We consider the possibility that the apparent under-prediction of dust mass in high redshift galaxies by chemical evolution models could be in part, due to these models neglecting this process, specifically due to their assumption that a large fraction of the dust mass is removed from the interstellar medium during star formation (so-called astration). By adding a planet formation phase into galaxy chemical evolution, we demonstrate that the dust budget crisis can be partially ameliorated by a factor of 1.3-1.5 only if a) circumstellar discs prevent a large fraction of the dust mass entering the star during its birth, and b) that dust mass is preferentially liberated via jets, winds and outflows rather than accreted into planetary-mass bodies.

We derive the low redshift galaxy stellar mass function (GSMF), inclusive of dust corrections, for the equatorial Galaxy And Mass Assembly (GAMA) dataset covering 180 deg$^2$. We construct the mass function using a density-corrected maximum volume method, using masses corrected for the impact of optically thick and thin dust. We explore the galactic bivariate brightness plane ($M_\star-\mu$), demonstrating that surface brightness effects do not systematically bias our mass function measurement above 10$^{7.5}$ M$_{\odot}$. The galaxy distribution in the $M-\mu$-plane appears well bounded, indicating that no substantial population of massive but diffuse or highly compact galaxies are systematically missed due to the GAMA selection criteria. The GSMF is {fit with} a double Schechter function, with $\mathcal M^\star=10^{10.78\pm0.01\pm0.20}M_\odot$, $\phi^\star_1=(2.93\pm0.40)\times10^{-3}h_{70}^3$Mpc$^{-3}$, $\alpha_1=-0.62\pm0.03\pm0.15$, $\phi^\star_2=(0.63\pm0.10)\times10^{-3}h_{70}^3$Mpc$^{-3}$, and $\alpha_2=-1.50\pm0.01\pm0.15$. We find the equivalent faint end slope as previously estimated using the GAMA-I sample, although we find a higher value of $\mathcal M^\star$. Using the full GAMA-II sample, we are able to fit the mass function to masses as low as $10^{7.5}$ $M_\odot$, and assess limits to $10^{6.5}$ $M_\odot$. Combining GAMA-II with data from G10-COSMOS we are able to comment qualitatively on the shape of the GSMF down to masses as low as $10^{6}$ $M_\odot$. Beyond the well known upturn seen in the GSMF at $10^{9.5}$ the distribution appears to maintain a single power-law slope from $10^9$ to $10^{6.5}$. We calculate the stellar mass density parameter given our best-estimate GSMF, finding $\Omega_\star= 1.66^{+0.24}_{-0.23}\pm0.97 h^{-1}_{70} \times 10^{-3}$, inclusive of random and systematic uncertainties.

Active galactic nuclei (AGN), particularly the most luminous AGN, are commonly assumed to be triggered through major mergers, however observational evidence for this scenario is mixed. To investigate any influence of galaxy mergers on AGN triggering and luminosities through cosmic time, we present a sample of 106 luminous X-ray selected type 1 AGN from the COSMOS survey. These AGN occupy a large redshift range (0.5 < z < 2.2) and two orders of magnitude in X-ray luminosity ($\sim$10$^{43}$ - 10$^{45}$ erg s$^{-1}$). AGN hosts are carefully mass and redshift matched to 486 control galaxies. A novel technique for identifying and quantifying merger features in galaxies is developed, subtracting GALFIT galaxy models and quantifying the residuals. Comparison to visual classification confirms this measure reliably picks out disturbance features in galaxies. No enhancement of merger features with increasing AGN luminosity is found with this metric, or by visual inspection. We analyse the redshift evolution of AGN associated with galaxy mergers and find no merger enhancement in lower redshift bins. Contrarily, in the highest redshift bin (z$\sim$2) AGN are $\sim$4 times more likely to be in galaxies exhibiting evidence of morphological disturbance compared to control galaxies, at 99% confidence level ($\sim$2.4$\sigma$) from visual inspection. Since only $\sim$15% of these AGN are found to be in morphologically disturbed galaxies, it is implied that major mergers at high redshift make a noticeable but subdominant contribution to AGN fuelling. At low redshifts other processes dominate and mergers become a less significant triggering mechanism.

Galaxy interactions are thought to be one of the main triggers of Active Galactic Nuclei (AGN), especially at high luminosities, where the accreted gas mass during the AGN lifetime is substantial. Evidence for a connection between mergers and AGN, however, remains mixed. Possible triggering mechanisms remain particularly poorly understood for luminous AGN, which are thought to require triggering by major mergers, rather than secular processes. We analyse the host galaxies of a sample of 20 optically and X-ray selected luminous AGN (log($L_{bol}$ [erg/s]) $>$ 45) at z $\sim$ 0.6 using HST WFC3 data in the F160W/H band. 15/20 sources have resolved host galaxies. We create a control sample of mock AGN by matching the AGN host galaxies to a control sample of non-AGN galaxies. Visual signs of disturbances are found in about 25% of sources in both the AGN hosts and control galaxies. Using both visual classification and quantitative morphology measures, we show that the levels of disturbance are not enhanced when compared to a matched control sample. We find no signs that major mergers play a dominant role in triggering AGN at high luminosities, suggesting that minor mergers and secular processes dominate AGN triggering up to the highest AGN luminosities. The upper limit on the enhanced fraction of major mergers is $\leqslant$20%. While major mergers might increase the incidence of (luminous AGN), they are not the prevalent triggering mechanism in the population of unobscured AGN.

We present the evolution in the number density and stellar mass functions of photometrically selected post-starburst galaxies in the UKIDSS Deep Survey (UDS), with redshifts of 0.5<z<2 and stellar masses logM>10. We find that this transitionary species of galaxy is rare at all redshifts, contributing ~5% of the total population at z~2, to <1% by z~0.5. By comparing the mass functions of quiescent galaxies to post-starburst galaxies at three cosmic epochs, we show that rapid quenching of star formation can account for 100% of quiescent galaxy formation, if the post-starburst spectral features are visible for ~250Myr. The flattening of the low mass end of the quiescent galaxy stellar mass function seen at z~1 can be entirely explained by the addition of rapidly quenched galaxies. Only if a significant fraction of post-starburst galaxies have features that are visible for longer than 250Myr, or they acquire new gas and return to the star-forming sequence, can there be significant growth of the red sequence from a slower quenching route. The shape of the mass function of these transitory post-starburst galaxies resembles that of quiescent galaxies at z~2, with a preferred stellar mass of logM~10.6, but evolves steadily to resemble that of star-forming galaxies at z<1. This leads us to propose a dual origin for post-starburst galaxies: (1) at z>2 they are exclusively massive galaxies that have formed the bulk of their stars during a rapid assembly period, followed by complete quenching of further star formation, (2) at z<1 they are caused by the rapid quenching of gas-rich star-forming galaxies, independent of stellar mass, possibly due to environment and/or gas-rich major mergers.

This paper is the second in a pair of articles presenting data release 1 (DR1) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), the largest single open-time key project carried out with the Herschel Space Observatory. The H-ATLAS is a wide-area imaging survey carried out in five photometric bands at 100, 160, 250, 350 and 500$\mu$m covering a total area of 600deg$^2$. In this paper we describe the identification of optical counterparts to submillimetre sources in DR1, comprising an area of 161 deg$^2$ over three equatorial fields of roughly 12$^\circ$x4.5$^\circ$ centred at 9$^h$, 12$^h$ and 14.5$^h$ respectively. Of all the H-ATLAS fields, the equatorial regions benefit from the greatest overlap with current multi-wavelength surveys spanning ultraviolet (UV) to mid-infrared regimes, as well as extensive spectroscopic coverage. We use a likelihood-ratio technique to identify SDSS counterparts at r<22.4 for 250-$\mu$m-selected sources detected at $\geq$ 4$\sigma$ ($\approx$28mJy). We find `reliable' counterparts (reliability R$\geq$0.8) for 44,835 sources (39 per cent), with an estimated completeness of 73.0 per cent and contamination rate of 4.7 per cent. Using redshifts and multi-wavelength photometry from GAMA and other public catalogues, we show that H-ATLAS-selected galaxies at $z<0.5$ span a wide range of optical colours, total infrared (IR) luminosities, and IR/UV ratios, with no strong disposition towards mid-IR-classified AGN in comparison with optical selection. The data described herein, together with all maps and catalogues described in the companion paper (Valiante et al. 2016), are available from the H-ATLAS website at www.h-atlas.org.

Despite decades of study, we still do not fully understand why some massive galaxies abruptly switch off their star formation in the early Universe, and what causes their rapid transition to the red sequence. Post-starburst galaxies provide a rare opportunity to study this transition phase, but few have currently been spectroscopically identified at high redshift ($z>1$). In this paper we present the spectroscopic verification of a new photometric technique to identify post-starbursts in high-redshift surveys. The method classifies the broad-band optical-near--infrared spectral energy distributions (SEDs) of galaxies using three spectral shape parameters (super-colours), derived from a principal component analysis of model SEDs. When applied to the multiwavelength photometric data in the UKIDSS Ultra Deep Survey (UDS), this technique identified over 900 candidate post-starbursts at redshifts $0.5<z<2.0$. In this study we present deep optical spectroscopy for a subset of these galaxies, in order to confirm their post-starburst nature. Where a spectroscopic assessment was possible, we find the majority (19/24 galaxies; ~80 per cent) exhibit the strong Balmer absorption (H $\delta$ equivalent width $W_{\lambda}$ >5 Ang.) and Balmer break, characteristic of post-starburst galaxies. We conclude that photometric methods can be used to select large samples of recently-quenched galaxies in the distant Universe.

We present a new morphological indicator designed for automated recognition of galaxies with faint asymmetric tidal features suggestive of an ongoing or past merger. We use this new indicator, together with preexisting diagnostics of galaxy structure to study the role of galaxy mergers in inducing (post-)starburst spectral signatures in local galaxies, and investigate whether (post-)starburst galaxies play a role in the build up of the `red sequence'. Our morphological and structural analysis of an evolutionary sample of 335 (post-)starburst galaxies in the SDSS DR7 with starburst ages 0<tSB<0.6 Gyr, shows that 45% of galaxies with young starbursts (tSB<0.1 Gyr) show signatures of an ongoing or past merger. This fraction declines with starburst age, and we find a good agreement between automated and visual classifications. The majority of the oldest (post-)starburst galaxies in our sample (tSB~0.6Gyr) have structural properties characteristic of early-type disks and are not as highly concentrated as the fully quenched galaxies commonly found on the `red sequence' in the present day Universe. This suggests that, if (post-)starburst galaxies are a transition phase between active star-formation and quiescence, they do not attain the structure of presently quenched galaxies within the first 0.6 Gyr after the starburst.

We combine Herschel/SPIRE sub-millimeter (submm) observations with existing multi-wavelength data to investigate the characteristics of low redshift, optically red galaxies detected in submm bands. We select a sample of galaxies in the redshift range 0.01$\leq$z$\leq$0.2, having >5$\sigma$ detections in the SPIRE 250 micron submm waveband. Sources are then divided into two sub-samples of $red$ and $blue$ galaxies, based on their UV-optical colours. Galaxies in the $red$ sample account for $\approx$4.2 per cent of the total number of sources with stellar masses M$_{*}\gtrsim$10$^{10}$ Solar-mass. Following visual classification of the $red$ galaxies, we find that $\gtrsim$30 per cent of them are early-type galaxies and $\gtrsim$40 per cent are spirals. The colour of the $red$-spiral galaxies could be the result of their highly inclined orientation and/or a strong contribution of the old stellar population. It is found that irrespective of their morphological types, $red$ and $blue$ sources occupy environments with more or less similar densities (i.e., the $\Sigma_5$ parameter). From the analysis of the spectral energy distributions (SEDs) of galaxies in our samples based on MAGPHYS, we find that galaxies in the $red$ sample (of any morphological type) have dust masses similar to those in the $blue$ sample (i.e. normal spiral/star-forming systems). However, in comparison to the $red$-spirals and in particular $blue$ systems, $red$-ellipticals have lower mean dust-to-stellar mass ratios. Besides galaxies in the $red$-elliptical sample have much lower mean star-formation/specific-star-formation rates in contrast to their counterparts in the $blue$ sample. Our results support a scenario where dust in early-type systems is likely to be of an external origin.

We present the GAMA Panchromatic Data Release (PDR) constituting over 230deg$^2$ of imaging with photometry in 21 bands extending from the far-UV to the far-IR. These data complement our spectroscopic campaign of over 300k galaxies, and are compiled from observations with a variety of facilities including: GALEX, SDSS, VISTA, WISE, and Herschel, with the GAMA regions currently being surveyed by VST and scheduled for observations by ASKAP. These data are processed to a common astrometric solution, from which photometry is derived for 221,373 galaxies with r<19.8 mag. Online tools are provided to access and download data cutouts, or the full mosaics of the GAMA regions in each band. We focus, in particular, on the reduction and analysis of the VISTA VIKING data, and compare to earlier datasets (i.e., 2MASS and UKIDSS) before combining the data and examining its integrity. Having derived the 21-band photometric catalogue we proceed to fit the data using the energy balance code MAGPHYS. These measurements are then used to obtain the first fully empirical measurement of the 0.1-500$\mu$m energy output of the Universe. Exploring the Cosmic Spectral Energy Distribution (CSED) across three time-intervals (0.3-1.1Gyr, 1.1-1.8~Gyr and 1.8---2.4~Gyr), we find that the Universe is currently generating $(1.5 \pm 0.3) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$, down from $(2.5 \pm 0.2) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$ 2.3~Gyr ago. More importantly, we identify significant and smooth evolution in the integrated photon escape fraction at all wavelengths, with the UV escape fraction increasing from 27(18)% at z=0.18 in NUV(FUV) to 34(23)% at z=0.06. The GAMA PDR will allow for detailed studies of the energy production and outputs of individual systems, sub-populations, and representative galaxy samples at $z<0.5$. The GAMA PDR can be found at: http://gama-psi.icrar.org/

In this work we present IRAM-30m telescope observations of a sample of bulge-dominated galaxies with large dust lanes, which have had a recent minor merger. We find these galaxies are very gas rich, with H2 masses between 4x10^8 and 2x10^10 Msun. We use these molecular gas masses, combined with atomic gas masses from an accompanying paper, to calculate gas-to-dust and gas-to-stellar mass ratios. The gas-to-dust ratios of our sample objects vary widely (between ~50 and 750), suggesting many objects have low gas-phase metallicities, and thus that the gas has been accreted through a recent merger with a lower mass companion. We calculate the implied minor companion masses and gas fractions, finding a median predicted stellar mass ratio of ~40:1. The minor companion likely had masses between ~10^7 - 10^10 Msun. The implied merger mass ratios are consistent with the expectation for low redshift gas-rich mergers from simulations. We then go on to present evidence that (no matter which star-formation rate indicator is used) our sample objects have very low star-formation efficiencies (star-formation rate per unit gas mass), lower even than the early-type galaxies from ATLAS3D which already show a suppression. This suggests that minor mergers can actually suppress star-formation activity. We discuss mechanisms that could cause such a suppression, include dynamical effects induced by the minor merger.

We present the properties of the first 250 $\mu$m blind sample of nearby galaxies (15 < D < 46 Mpc) containing 42 objects from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). Herschel's sensitivity probes the faint end of the dust luminosity function for the first time, spanning a range of stellar mass (7.4 < log$_{10}$ M$_{\star}$ < 11.3 M$_{\odot}$), star formation activity (-11.8 < log$_{10}$ SSFR < -8.9 yr$^{-1}$), gas fraction (3-96 per cent), and colour (0.6 < FUV-Ks < 7.0 mag). The median cold dust temperature is 14.6 K, colder than in the Herschel Reference Survey (18.5 K) and Planck Early Release Compact Source Catalogue 17.7 K. The mean dust-to-stellar mass ratio in our sample is higher than these surveys by factors of 3.7 and 1.8, with a dust mass volume density of (3.7 $\pm$ 0.7) x 10$^{5}$ M$_{\odot}$ Mpc$^{-3}$. Counter-intuitively, we find that the more dust rich a galaxy, the lower its UV attenuation. Over half of our dust-selected sample are very blue in FUV-Ks colour, with irregular and/or highly flocculent morphology, these galaxies account for only 6 per cent of the sample's stellar mass but contain over 35 per cent of the dust mass. They are the most actively star forming galaxies in the sample, with the highest gas fractions and lowest UV attenuation. They also appear to be in an early stage of converting their gas into stars, providing valuable insights into the chemical evolution of young galaxies.

The Herschel Space Observatory has had a tremendous impact on the study of extragalactic dust. Specifically, early-type galaxies (ETG) have been the focus of several studies. In this paper we combine results from two Herschel studies - a Virgo cluster study HeViCS and a broader, low-redshift H-ATLAS/GAMA study - and contrast the dust and associated properties for similar mass galaxies. This comparison is motivated by differences in results exhibited between multiple Herschel studies of early-type galaxies. A comparison between consistent modified blackbody derived dust mass is carried out, revealing strong differences between the two samples in both dust mass and dust-to-stellar mass ratio. In particular, the HeViCS sample lacks massive ETG with as high a specific dust content as found in H-ATLAS. This is most likely connected with the difference in environment for the two samples. We calculate nearest neighbour environment densities in a consistent way, showing that H-ATLAS ETG occupy sparser regions of the local Universe, whereas HeViCS ETG occupy dense regions. This is also true for ETG that are not Herschel-detected but are in the Virgo and GAMA parent samples. Spectral energy distributions are fit to the panchromatic data. From these we find that in H-ATLAS the specific star formation rate anticorrelates with stellar mass and reaches values as high as in our Galaxy. On the other hand HeViCS ETG appear to have little star formation. Based on the trends found here, H-ATLAS ETG are thought to have more extended star formation histories and a younger stellar population than HeViCS ETG.

We present the evolution of dust and molecular gas properties in a sample of 11 $z\sim0.03$ starburst to post-starburst (PSB) galaxies selected to span an age sequence from ongoing starburst to 1 Gyr after the starburst ended. All PSBs harbour significant molecular gas and dust reservoirs and residual star formation, indicating that complete quenching of the starburst due to exhaustion or expulsion of gas has not occurred during this timespan. As the starburst ages, we observe a clear decrease in the star-formation efficiency, molecular gas and SFR surface density, and effective dust temperature, from levels coincident with starburst galaxies to those of normal star-forming galaxies. These trends are consistent with a natural decrease in the SFR following consumption of molecular gas by the starburst, and corresponding decrease in the interstellar radiation field strength as the starburst ages. The gas and dust contents of the PSBs are coincident with those of star-forming galaxies and molecular gas-rich early-type galaxies, and are not consistent with galaxies on the red-sequence. We find no evidence that the global gas reservoir is expelled by stellar winds or AGN feedback. Our results show that although a strong starburst in a low-redshift galaxy may cause the galaxy to ultimately have a lower specific SFR and be of an earlier morphological type, the galaxy will remain in the "green valley" for an extended time. Multiple such episodes may be needed to complete migration of the galaxy from the blue- to red-sequence.

We report an unexpected variation in the positional offset distributions between Herschel-ATLAS sub-millimetre (submm) sources and their optical associations, depending on both 250-{\mu}m signal-to-noise ratio and 250/350-{\mu}m colour. We show that redder and brighter submm sources have optical associations with a broader distribution of positional offsets than would be expected if these offsets were due to random positional errors in the source extraction. The observation can be explained by two possible effects: either red submm sources trace a more clustered population than blue ones, and their positional errors are increased by confusion; or red submm sources are generally at high redshifts and are frequently associated with low-redshift lensing structures which are identified as false counterparts. We perform various analyses of the data, including the multiplicity of optical associations, the redshift and magnitude distributions in H-ATLAS in comparison to HerMES, and simulations of weak lensing, and we conclude that the effects are most likely to be explained by widespread weak lensing of Herschel-SPIRE sources by foreground structures. This has important consequences for counterpart identification and derived redshift distributions and luminosity functions of submm surveys.

[Abridged] Aims: This work focuses on one lensed system, HATLAS J142935.3-002836 (H1429-0028), selected in the Herschel-ATLAS field. Gathering a rich, multi-wavelength dataset, we aim to confirm the lensing hypothesis and model the background source's morphology and dynamics, as well as to provide a full physical characterisation. Methods: Multi-wavelength high-resolution data is utilised to assess the nature of the system. A lensing-analysis algorithm which simultaneously fits different wavebands is adopted to characterise the lens. The background galaxy dynamical information is studied by reconstructing the 3-D source-plane of the ALMA CO(J:4-3) transition. Near-IR imaging from HST and Keck-AO allows to constrain rest-frame optical photometry independently for the foreground and background systems. Physical parameters (such as stellar and dust masses) are estimated via modelling of the spectral energy distribution taking into account source blending, foreground obscuration, and differential magnification. Results: The system comprises a foreground edge-on disk galaxy (at z_sp=0.218) with an almost complete Einstein ring around it. The background source (at z_sp=1.027) is magnified by a factor of ~8-10 depending on wavelength. It is comprised of two components and a tens of kpc long tidal tail resembling the Antennae merger. As a whole, the system is a massive stellar system (1.32[-0.41,+0.63] x1E11 Mo) forming stars at a rate of 394+-90 Mo/yr, and has a significant gas reservoir M_ISM = 4.6+-1.7 x1E10 Mo. Its depletion time due to star formation alone is thus expected to be tau_SF=M_ISM/SFR=117+-51 Myr. The dynamical mass of one of the components is estimated to be 5.8+-1.7 x1E10 Mo, and, together with the photometric total mass estimate, it implies that H1429-0028 is a major merger system (1:2.8[-1.5,+1.8]).

We present a comparison of the physical properties of a rest-frame $250\mu$m selected sample of massive, dusty galaxies from $0<z<5.3$. Our sample comprises 29 high-redshift submillimetre galaxies (SMGs) from the literature, and 843 dusty galaxies at $z<0.5$ from the Herschel-ATLAS, selected to have a similar stellar mass to the SMGs. The $z>1$ SMGs have an average SFR of $390^{+80}_{-70}\,$M$_\odot$yr$^{-1}$ which is 120 times that of the low-redshift sample matched in stellar mass to the SMGs (SFR$=3.3\pm{0.2}$ M$_\odot$yr$^{-1}$). The SMGs harbour a substantial mass of dust ($1.2^{+0.3}_{-0.2}\times{10}^9\,$M$_\odot$), compared to $(1.6\pm0.1)\times{10}^8\,$M$_\odot$ for low-redshift dusty galaxies. At low redshifts the dust luminosity is dominated by the diffuse ISM, whereas a large fraction of the dust luminosity in SMGs originates from star-forming regions. At the same dust mass SMGs are offset towards a higher SFR compared to the low-redshift H-ATLAS galaxies. This is not only due to the higher gas fraction in SMGs but also because they are undergoing a more efficient mode of star formation, which is consistent with their bursty star-formation histories. The offset in SFR between SMGs and low-redshift galaxies is similar to that found in CO studies, suggesting that dust mass is as good a tracer of molecular gas as CO.

We apply a chemical evolution model to investigate the sources and evolution of dust in a sample of 26 high-redshift ($z>1$) submillimetre galaxies (SMGs) from the literature, with complete photometry from ultraviolet to the submillimetre. We show that dust produced only by low-intermediate mass stars falls a factor 240 short of the observed dust masses of SMGs, the well-known `dust-budget crisis'. Adding an extra source of dust from supernovae can account for the dust mass in 19 per cent of the SMG sample. Even after accounting for dust produced by supernovae the remaining deficit in the dust mass budget provides support for higher supernova yields, substantial grain growth in the interstellar medium or a top-heavy IMF. Including efficient destruction of dust by supernova shocks increases the tension between our model and observed SMG dust masses. The models which best reproduce the physical properties of SMGs have a rapid build-up of dust from both stellar and interstellar sources and minimal dust destruction. Alternatively, invoking a top-heavy IMF or significant changes in the dust grain properties can solve the dust budget crisis only if dust is produced by both low mass stars and supernovae and is not efficiently destroyed by supernova shocks.

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.

We present an analysis of CO molecular gas tracers in a sample of 500{\mu}m-selected Herschel-ATLAS galaxies at z<0.05 (cz<14990km/s). Using 22-500{\mu}m photometry from WISE, IRAS and Herschel, with HI data from the literature, we investigate correlations between warm and cold dust, and tracers of the gas in different phases. The correlation between global CO(3-2) line fluxes and FIR-submillimetre fluxes weakens with increasing IR wavelength ({\lambda}>60{\mu}m), as a result of colder dust being less strongly associated with dense gas. Conversely, CO(2-1) and HI line fluxes both appear to be better correlated with longer wavelengths, suggesting that cold dust is more strongly associated with diffuse atomic and molecular gas phases, consistent with it being at least partially heated by radiation from old stellar populations. The increased scatter at long wavelengths implies that submillimetre fluxes are a poorer tracer of SFR. Fluxes at 22 and 60{\mu}m are also better correlated with diffuse gas tracers than dense CO(3-2), probably due to very-small-grain emission in the diffuse interstellar medium, which is not correlated with SFR. The FIR/CO luminosity ratio and the dust mass/CO luminosity ratio both decrease with increasing luminosity, as a result of either correlations between mass and metallicity (changing CO/H2) or between CO luminosity and excitation [changing CO(3-2)/CO(1-0)].

We use multi-wavelength (0.12 - 500 micron) photometry from Herschel-ATLAS, WISE, UKIDSS, SDSS and GALEX, to study 23 nearby spheroidal galaxies with prominent dust lanes (DLSGs). DLSGs are considered to be remnants of recent minor mergers, making them ideal laboratories for studying both the interstellar medium (ISM) of spheroids and minor-merger-driven star formation in the nearby Universe. The DLSGs exhibit star formation rates (SFRs) between 0.01 and 10 MSun yr^-1, with a median of 0.26 MSun yr^-1 (a factor of 3.5 greater than the average SG). The median dust mass, dust-to-stellar mass ratio and dust temperature in these galaxies are around 10^7.6 MSun yr^-1, ~0.05% and ~19.5 K respectively. The dust masses are at least a factor of 50 greater than that expected from stellar mass loss and, like the SFRs, show no correlation with galaxy luminosity, suggesting that both the ISM and the star formation have external drivers. Adopting literature gas-to-dust ratios and star formation histories derived from fits to the panchromatic photometry, we estimate that the median current and initial gas-to-stellar mass ratios in these systems are ~4% and ~7% respectively. If, as indicated by recent work, minor mergers that drive star formation in spheroids with (NUV-r)>3.8 (the colour range of our DLSGs) have stellar mass ratios between 1:6 and 1:10, then the satellite gas fractions are likely >50%.

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.

We present two large, nearby (0.013$\le$z$\le$0.06) samples of Early-Type Galaxies (ETGs): a visually classified sample of 220 ETGs, created using source-matched data from the Galaxy and Mass Assembly (GAMA) database with FIR/sub-mm detections from $Herschel$-ATLAS; and a visually classified sample of 551 ETGs which are undetected with $Herschel$-ATLAS. Active galactic nuclei (AGN) are removed from our samples using optical emission line diagnostics. These samples are scrutinised to determine characteristics of sub-mm detected versus undetected ETGs. We find similarities in the stellar mass distributions of the two ETG samples but testing other properties uncovers significant differences. The sub-mm detected sample is shown to have lower concentration and S\'ersic indices than those in the undetected sample - a result which may be linked to the presence of dust in the former. Optical and UV-optical colours are also shown to be much bluer, indicating that the dust is linked with recent star formation. The intrinsic effective radii are on average 1.5 times larger for the sub-mm detected ETGs. Surface densities and groups data from the GAMA database are examined for the two samples, leading to the conclusion that dusty ETGs inhabit sparser environments than non-dusty ETGs in the nearby universe, although environments of the brightest ETGs are shown to differ the least. Modified Planck functions are fit to the H-ATLAS detected PACS and SPIRE fluxes for ETGs with sub-mm flux densities of at least 3$\sigma$ in the 350$\mu$m SPIRE band, giving a resultant mean cold dust temperature of T$_{d}$=22.1K, with a range of 9-30K. The corresponding mean dust mass is 1.8$\times10^{7}$M$_{\odot}$, with a range of (0.08-35.0)$\times10^{7}$M$_{\odot}$. The dust masses calculated from these fits, normalised by stellar mass, are shown to increase with decreasing stellar mass and bluer colours.

We present a pan-chromatic analysis of an unprecedented sample of 1402 250 micron-selected galaxies at z < 0.5 (mean z = 0.24) from the Herschel-ATLAS survey. We complement our Herschel 100-500 micron data with UV-K-band photometry from the Galaxy And Mass Assembly (GAMA) survey and apply the MAGPHYS energy-balance technique to produce pan-chromatic SEDs for a representative sample of 250 micron selected galaxies spanning the most recent 5 Gyr of cosmic history. We derive estimates of physical parameters, including star formation rates, stellar masses, dust masses and infrared luminosities. The typical H-ATLAS galaxy at z < 0.5 has a far-infrared luminosity in the range 10^10 - 10^12 Lsolar (SFR: 1-50 Msolar/yr) thus is broadly representative of normal star forming galaxies over this redshift range. We show that 250 micron-selected galaxies contain a larger mass of dust at a given infra-red luminosity or star formation rate than previous samples selected at 60 micron from IRAS. We derive typical SEDs for H-ATLAS galaxies, and show that the emergent SED shape is most sensitive to specific star formation rate. The optical-UV SEDs also become more reddened due to dust at higher redshifts. Our template SEDs are significantly cooler than existing infra-red templates. They may therefore be most appropriate for inferring total IR luminosities from moderate redshift submillimetre selected samples and for inclusion in models of the lower redshift submillimetre galaxy populations.

While the selection of strongly lensed galaxies with 500{\mu}m flux density S_500>100 mJy has proven to be rather straightforward (Negrello et al. 2010), for many applications it is important to analyze samples larger than the ones obtained when confining ourselves to such a bright limit. Moreover, only by probing to fainter flux densities is possible to exploit strong lensing to investigate the bulk of the high-z star-forming galaxy population. We describe HALOS (the Herschel-ATLAS Lensed Objects Selection), a method for efficiently selecting fainter candidate strongly lensed galaxies, reaching a surface density of ~1.5-2 deg^-2, i.e. a factor of about 4 to 6 higher than that at the 100 mJy flux limit. HALOS will allow the selection of up to ~1000 candidate strongly lensed galaxies (with amplifications \mu>2) over the full H-ATLAS survey area. Applying HALOS to the H-ATLAS Science Demonstration Phase field (~14.4 deg^2) we find 31 candidate strongly lensed galaxies, whose candidate lenses are identified in the VIKING near-infrared catalog. Using the available information on candidate sources and candidate lenses we tentatively estimate a ~72% purity of the sample. The redshift distribution of the candidate lensed sources is close to that reported for most previous surveys for lensed galaxies, while that of candidate lenses extends to substantially higher redshifts than found in the other surveys. The counts of candidate strongly lensed galaxies are also in good agreement with model predictions (Lapi et al. 2011). Even though a key ingredient of the method is the deep near-infrared VIKING photometry, we show that H-ATLAS data alone allow the selection of a similarly deep sample of candidate strongly lensed galaxies with an efficiency close to 50%; a slightly lower surface density (~1.45 deg^-2) can be reached with a ~70% efficiency.

We present a simple, physically-motivated model to interpret consistently the emission from galaxies at ultraviolet, optical and infrared wavelengths. We combine this model with a Bayesian method to obtain robust statistical constraints on key parameters describing the stellar content, star formation activity and dust content of galaxies. Our model is now publicly available via a user-friendly code package, MAGPHYS at www.iap.fr/magphys. We present an application of this model to interpret a sample of ~1400 local (z<0.5) galaxies from the H-ATLAS survey. We find that, for these galaxies, the diffuse interstellar medium, powered mainly by stars older than 10 Myr, accounts for about half the total infrared luminosity. We discuss the implications of this result to the use of star formation rate indicators based on total infrared luminosity.

We present the dust properties and star-formation histories of local submillimetre-selected galaxies in Herschel-ATLAS, classified by optical morphology. The early-type galaxies (ETGs) that are detected contain as much dust as typical spirals, and form a unique sample that has been blindly selected at submillimetre wavelengths. Comparing H-ATLAS galaxies to a control sample of optically selected galaxies, we find 5.5% of luminous ETGs are detected in H-ATLAS. The H-ATLAS ETGs contain a significant mass of cold dust: the mean dust mass is 5.5x10^7 Msun, with individual galaxies ranging from 9x10^5-4x10^8 Msun. This is comparable to that of spirals in our sample, and is an order of magnitude more dust than that found for the control ETGs, which have a median dust mass inferred from stacking of (0.8-4.0)x10^6 Msun. The ETGs detected in H-ATLAS have bluer NUV-r colours, higher specific star-formation rates and younger stellar populations than ETGs which are optically selected, and may be transitioning from the blue cloud to the red sequence. We also find that H-ATLAS and control ETGs inhabit similar low-density environments. We conclude that the dust in H-ATLAS and control ETGs cannot be solely from stellar sources, and a large contribution from dust formed in the ISM or external sources is required. Alternatively, dust destruction may not be as efficient as predicted. We also explore the properties of the most passive spiral galaxies in our sample with SSFR<10^-11/yr. We find these passive spirals have lower dust-to-stellar mass ratios, higher stellar masses and older stellar population ages than normal spirals. The passive spirals inhabit low density environments similar to those of the normal spiral galaxies in our sample. This shows that the processes which turn spirals passive do not occur solely in the intermediate density environments of group and cluster outskirts. (Abridged)

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-ATLAS. The sample consists of galaxies selected at 250{\mu}m which have reliable counterparts from 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 using 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 3-4 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 timescale together with either a more top-heavy 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 ISM mass, and points to an enhanced supply of fuel for star formation at earlier cosmic epochs.