We present Atacama Large Millimeter/sub-millimeter Array (ALMA) observations towards 27 low-redshift ($0.02< z<0.2$) star-forming galaxies taken from the Valpara\'iso ALMA/APEX Line Emission Survey (VALES). We perform stacking analyses of the $^{12}$CO($1-0$), $^{13}$CO($1-0$) and C$^{18}$O($1-0$) emission lines to explore the $L'$ ($^{12}$CO($1-0$))/$L'$($^{13}$CO($1-0$))) (hereafter $L'$($^{12}$CO)/$L'$($^{13}$CO)) and $L'$($^{13}$CO($1-0$))/$L'$(C$^{18}$O($1-0$)) (hereafter $L'$($^{13}$CO)/$L'$(C$^{18}$O) line luminosity ratio dependence as a function of different global galaxy parameters related to the star formation activity. The sample has far-IR luminosities $10^{10.1-11.9}$L$_{\odot}$ and stellar masses of $10^{9.8-10.9}$M$_{\odot}$ corresponding to typical star-forming and starburst galaxies at these redshifts. On average we find a $L'$($^{12}$CO)/$L'$($^{13}$CO) line luminosity ratio value of 16.1$\pm$2.5. Galaxies with evidences of possible merging activity tend to show higher $L'$($^{12}$CO)/$L'$($^{13}$CO) ratios by a factor of two, while variations of this order are also found in galaxy samples with higher star formation rates or star formation efficiencies. We also find an average $L'$($^{13}$CO)/$L'$(C$^{18}$O) line luminosity ratio of 2.5$\pm$0.6, which is in good agreement with those previously reported for starburst galaxies. We find that galaxy samples with high $L_{\text{IR}}$, SFR and SFE show low $L'$($^{13}$CO)/$L'$(C$^{18}$O) line luminosity ratios with high $L'$($^{12}$CO)/$L'$($^{13}$CO) line luminosity ratios, suggesting that these trends are produced by selective enrichment of massive stars in young starbursts.

On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event.

Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the external medium generates external shock waves, responsible for the afterglow emission, which lasts from days to months, and occurs over a broad energy range, from the radio to the GeV bands. The afterglow emission is generally well explained as synchrotron radiation by electrons accelerated at the external shock. Recently, an intense, long-lasting emission between 0.2 and 1 TeV was observed from the GRB 190114C. Here we present the results of our multi-frequency observational campaign of GRB~190114C, and study the evolution in time of the GRB emission across 17 orders of magnitude in energy, from $5\times10^{-6}$ up to $10^{12}$\,eV. We find that the broadband spectral energy distribution is double-peaked, with the TeV emission constituting a distinct spectral component that has power comparable to the synchrotron component. This component is associated with the afterglow, and is satisfactorily explained by inverse Compton upscattering of synchrotron photons by high-energy electrons. We find that the conditions required to account for the observed TeV component are not atypical, supporting the possibility that inverse Compton emission is commonly produced in GRBs.

Using the EMIR instrument on the IRAM 30m telescope, we conducted a spectroscopic redshift search of seven z$_{\rm phot}$ $\sim$ 4 sub-millimetre bright galaxies selected from the Herschel Bright Sources (HerBS) sample with fluxes at 500 $\mu$m greater than 80 mJy. For four sources, we obtained spectroscopic redshifts between 3.4 < z < 4.1 through the detection of multiple CO-spectral lines with J $\leq$ 3. Later, we detected low-J transitions for two of these sources with the GBT including the CO(1-0) transition. For the remaining three sources, more data are needed to determine the spectroscopic redshift unambiguously. The measured CO luminosities and line widths suggest that all these sources are gravitationally lensed. These observations demonstrate that the 2 mm window is indispensable to confirm robust spectroscopic redshifts for z < 4 sources. Finally, we present an efficient graphical method to correctly identify spectroscopic redshifts.

We analyze an extremely deep 450-$\mu$m image ($1\sigma=0.56$\,mJy\,beam$^{-1}$) of a $\simeq 300$\,arcmin$^{2}$ area in the CANDELS/COSMOS field as part of the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). We select a robust (signal-to-noise ratio $\geqslant 4$) and flux-limited ($\geqslant 4$\,mJy) sample of 164 sub-millimeter galaxies (SMGs) at 450-$\mu$m that have $K$-band counterparts in the COSMOS2015 catalog identified from radio or mid-infrared imaging. Utilizing this SMG sample and the 4705 $K$-band-selected non-SMGs that reside within the noise level $\leqslant 1$\,mJy\,beam$^{-1}$ region of the 450-$\mu$m image as a training set, we develop a machine-learning classifier using $K$-band magnitude and color-color pairs based on the thirteen-band photometry available in this field. We apply the trained machine-learning classifier to the wider COSMOS field (1.6\,deg$^{2}$) using the same COSMOS2015 catalog and identify a sample of 6182 450-$\mu$m SMG candidates with similar colors. The number density, radio and/or mid-infrared detection rates, redshift and stellar mass distributions, and the stacked 450-$\mu$m fluxes of these SMG candidates, from the S2COSMOS observations of the wide field, agree with the measurements made in the much smaller CANDELS field, supporting the effectiveness of the classifier. Using this 450-$\mu$m SMG candidate sample, we measure the two-point autocorrelation functions from $z=3$ down to $z=0.5$. We find that the 450-$\mu$m SMG candidates reside in halos with masses of $\simeq (2.0\pm0.5) \times10^{13}\,h^{-1}\,\rm M_{\odot}$ across this redshift range. We do not find evidence of downsizing that has been suggested by other recent observational studies.

Studying the nature of various types of supernovae (SNe) is important for our understanding of stellar evolution. Observations of atomic and molecular gas in the host galaxies of gamma-ray bursts (GRBs) and SNe have recently been used to learn about the nature of the explosions themselves and the star formation events during which their progenitors were born. Based on archival data for M74, which previously has not been investigated in the context of SN positions, we report the gas properties in the environment of the broad-lined type Ic (Ic-BL) SN 2002ap and the type II SNe 2003gd, 2013ej, and 2019krl. The SN 2002ap is located at the end of an off-centre, asymmetric, 55 kpc-long HI extension containing 7.5% of the total atomic gas in M74, interpreted as a signature of external gas accretion. It is the fourth known case of an explosion of a presumably massive star located close to a concentration of atomic gas (after GRBs 980425, 060505, and SN 2009bb). It is unlikely that all these associations are random (at a 3sigma significance), so the case of SN 2002ap adds to the evidence that the birth of the progenitors of type Ic-BL SNe and GRBs is connected with the accretion of atomic gas from the intergalactic medium. The HI extension could come from tidally disrupted companions of M74, or be a remnant of a galaxy or a gas cloud that accreted entirely from the intragroup medium. The other (type II) SNe in M74 are located at the outside edge of a spiral arm. This suggests that either their progenitors were born when gas was piling up there or that the SN progenitors moved away from the arm due to their orbital motions. These type II SNe do not seem to be related to gas accretion.

We measure the 850-$\mu$m source densities of 46 candidate protoclusters selected from the Planck High-z catalogue (PHz) and the Planck Catalogue of Compact Sources (PCCS) that were followed up with Herschel-SPIRE and SCUBA-2. This paper aims to search for overdensities of 850-$\mu$m sources in order to select the fields that are most likely to be genuine protoclusters. Of the 46 candidate protoclusters, 25 have significant overdensities ($>$5 times the field counts), 11 have intermediate overdensities (3--5 times the field counts) and 10 have no overdensity ($<$3 times the field counts) of 850-$\mu$m sources. We find that the enhanced number densities are unlikely to be the result of sample variance. Compared with the number counts of another sample selected from Planck's compact source catalogues, this [PHz+PCCS]-selected sample has a higher fraction of candidate protoclusters with significant overdensities, though both samples show overdensities of 850-$\mu$m sources above intermediate level. Based on the estimated star-formation rate densities (SFRDs), we suggest that both samples can efficiently select protoclusters with starbursting galaxies near the redshift at which the global field SFRD peaks ($2 < z < 3$). Based on the confirmation of overdensities found here, future follow-up observations on other PHz targets may greatly increase the number of genuine DSFG-rich clusters/protoclusters.

We construct a SCUBA-2 450-$\mu$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$\sigma$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $\mu$m with signal-to-noise ratio $>$ 4.0 and analyze the physical properties of their multi-wavelength counterparts. We find that most of the sources are at $z\lesssim3$, with a median of $z = 1.79^{+0.03}_{-0.15}$. About $35^{+32}_{-25}$% of our sources are classified as starburst galaxies based on their total star-formation rates (SFRs) and stellar masses ($M_{\ast}$). By fitting the far-infrared spectral energy distributions, we find that our 450-$\mu$m-selected sample has a wide range of dust temperatures (20 K $ \lesssim T_{\rm d} \lesssim$ 60 K), with a median of ${T}_{\rm d} = 38.3^{+0.4}_{-0.9}$ K. We do not find a redshift evolution in dust temperature for sources with $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ at $z<3$. However, we find a moderate correlation where dust temperature increases with the deviation from the SFR-$M_{\ast}$ relation. The increase in dust temperature also correlates with optical morphology, which is consistent with merger-triggered starbursts in sub-millimeter galaxies. Our galaxies do not show the tight IRX-$\beta_{\rm UV}$ correlation that has been observed in the local Universe. We construct the infrared luminosity functions of our 450-$\mu$m sources and measure their comoving SFR densities. The contribution of the $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ population to the SFR density rises dramatically from $z$ = 0 to 2 ($\propto$ ($1+z$)$^{3.9\pm1.1}$) and dominates the total SFR density at $z \gtrsim 2$.

We investigate the evolution of the gas mass fraction for galaxies in the COSMOS field using submillimetre emission from dust at 850$\mu$m. We use stacking methodologies on the 850$\mu$m S2COSMOS map to derive the gas mass fraction of galaxies out to high redshifts, 0 <= $z$ <= 5, for galaxies with stellar masses of $10^{9.5} < M_* (\rm M_{\odot}) < 10^{11.75}$. In comparison to previous literature studies we extend to higher redshifts, include more normal star-forming galaxies (on the main sequence), and also investigate the evolution of the gas mass fraction split by star-forming and passive galaxy populations. We find our stacking results broadly agree with scaling relations in the literature. We find tentative evidence for a peak in the gas mass fraction of galaxies at around $z$ ~ 2.5-3, just before the peak of the star formation history of the Universe. We find that passive galaxies are particularly devoid of gas, compared to the star-forming population. We find that even at high redshifts, high stellar mass galaxies still contain significant amounts of gas.

We obtained CO(2-1) observations of seven GRB hosts with the APEX and IRAM 30m telescopes. We analysed these data together with all other hosts with previous CO observations. We obtained detections for 3 GRB hosts (980425, 080207, and 111005A) and upper limits for the remaining 4 (031203, 060505, 060814, and 100316D). In our entire sample of 12 CO-observed GRB hosts, 3 are clearly deficient in molecular gas, even taking into account their metallicity (980425, 060814, and 080517). Four others are close to the best-fit line for other star-forming galaxies on the SFR-MH2 plot (051022, 060505, 080207, and 100316D). One host is clearly molecule rich (111005A). Finally, the data for 4 GRB hosts are not deep enough to judge whether they are molecule deficient (000418, 030329, 031203, and 090423). The median value of the molecular gas depletion time, MH2/SFR, of GRB hosts is ~0.3 dex below that of other star-forming galaxies, but this result has low statistical significance. A Kolmogorov-Smirnov test performed on MH2/SFR shows an only ~2sigma difference between GRB hosts and other galaxies. This difference can partly be explained by metallicity effects, since the significance decreases to ~1sigma for MH2/SFR versus~metallicity. We found that any molecular gas deficiency of GRB hosts has low statistical significance and that it can be attributed to their lower metallicities; and thus the sample of GRB hosts has molecular properties that are consistent with those of other galaxies, and they can be treated as representative star-forming galaxies. Given the concentration of atomic gas recently found close to GRB and supernova sites, indicating recent gas inflow, our results about the weak molecular deficiency imply that such an inflow does not enhance the SFRs significantly, or that atomic gas converts efficiently into the molecular phase, which fuels star formation.

Using data from ALMA and near-infrared (NIR) integral field spectrographs including both SINFONI and KMOS on the VLT, we investigate the two-dimensional distributions of H$\alpha$ and rest-frame far-infrared (FIR) continuum in six submillimeter galaxies at $z\sim2$. At a similar spatial resolution ($\sim$0.5" FWHM; $\sim$4.5 kpc at $z=2$), we find that the half-light radius of H$\alpha$ is significantly larger than that of the FIR continuum in half of the sample, and on average H$\alpha$ is a median factor of $2.0\pm0.4$ larger. Having explored various ways to correct for the attenuation, we find that the attenuation-corrected H$\alpha$-based SFRs are systematically lower than the IR-based SFRs by at least a median factor of $3\pm1$, which cannot be explained by the difference in half-light radius alone. In addition, we find that in 40% of cases the total $V$-band attenuation ($A_V$) derived from energy balance modeling of the full ultraviolet(UV)-to-FIR spectral energy distributions (SEDs) is significantly higher than that derived from SED modeling using only the UV-to-NIR part of the SEDs, and the discrepancy appears to increase with increasing total infrared luminosity. Finally, considering all our findings along with the studies in the literature, we postulate that the dust distributions in SMGs, and possibly also in less IR luminous $z\sim2$ massive star-forming galaxies, can be decomposed into three main components; the diffuse dust heated by older stellar populations, the more obscured and extended young star-forming HII regions, and the heavily obscured central regions that have a low filling factor but dominate the infrared luminosity in which the majority of attenuation cannot be probed via UV-to-NIR emissions.

We analyse 870um Atacama Large Millimetre Array (ALMA) dust continuum detections of 41 canonically-selected z~3 Lyman-break galaxies (LBGs), as well as 209 ALMA-undetected LBGs, in follow-up of SCUBA-2 mapping of the UKIDSS Ultra Deep Survey (UDS) field. We find that our ALMA-bright LBGs lie significantly off the locally calibrated IRX-beta relation and tend to have relatively bluer rest-frame UV slopes (as parametrised by beta), given their high values of the 'infrared excess' (IRX=L_IR/L_UV), relative to the average 'local' IRX-beta relation. We attribute this finding in part to the young ages of the underlying stellar populations but we find that the main reason behind the unusually blue UV slopes are the relatively shallow slopes of the corresponding dust attenuation curves. We show that, when stellar masses are being established via SED fitting, it is absolutely crucial to allow the attenuation curves to vary (rather than fixing it on Calzetti-like law), where we find that the inappropriate curves may underestimate the resulting stellar masses by a factor of ~2-3x on average. In addition, we find these LBGs to have relatively high specific star-formation rates (sSFRs), dominated by the dust component, as quantified via the fraction of obscured star formation ( f_obs = SFR_IR/SFR_(UV+IR)). We conclude that the ALMA-bright LBGs are, by selection, massive galaxies undergoing a burst of a star formation (large sSFRs, driven, for example, by secular or merger processes), with a likely geometrical disconnection of the dust and stars, responsible for producing shallow dust attenuation curves.

We present sensitive 850$\mu$m imaging of the COSMOS field using 640hr of new and archival observations taken with SCUBA-2 at the East Asian Observatory's James Clerk Maxwell Telescope. The SCUBA-2 COSMOS survey (S2COSMOS) achieves a median noise level of $\sigma_{850\mu{\mathrm{m}}}$=1.2mJy/beam over an area of 1.6 sq. degree (MAIN; HST/ACS footprint), and $\sigma_{850\mu{\mathrm{m}}}$=1.7mJy/beam over an additional 1 sq. degree of supplementary (SUPP) coverage. We present a catalogue of 1020 and 127 sources detected at a significance level of >4$\sigma$ and >4.3$\sigma$ in the MAIN and SUPP regions, respectively, corresponding to a uniform 2% false-detection rate. We construct the single-dish 850$\mu$m number counts at $S_{850}$>2mJy and show that these S2COSMOS counts are in agreement with previous single-dish surveys, demonstrating that degree-scale fields are sufficient to overcome the effects of cosmic variance in the $S_{850}$=2-10mJy population. To investigate the properties of the galaxies identified by S2COSMOS sources we measure the surface density of near-infrared-selected galaxies around their positions and identify an average excess of 2.0$\pm$0.2 galaxies within a 13$''$ radius (~100kpc at $z$~2). The bulk of these galaxies represent near-infrared-selected SMGs and/or spatially-correlated sources and lie at a median photometric redshift of $z$=2.0$\pm$0.1. Finally, we perform a stacking analysis at sub-millimeter and far-infrared wavelengths of stellar-mass-selected galaxies ($M_{\star}$=10$^{10}$-10$^{12}{\rm M_{\odot}}$) from $z$=0-4, obtaining high-significance detections at 850um in all subsets (SNR=4-30), and investigate the relation between far-infrared luminosity, stellar mass, and the peak wavelength of the dust SED. The publication of this survey adds a new deep, uniform sub-millimeter layer to the wavelength coverage of this well-studied COSMOS field.

We present radio and optical afterglow observations of the TeV-bright long Gamma Ray Burst (GRB) 190114C, which was detected by the MAGIC telescope at a redshift of z=0.425. Our observations with ALMA, ATCA, and uGMRT were obtained by our low frequency observing campaign and range from ~1 to ~140 days after the burst and the optical observations were done with the 0.7-m GROWTH-India telescope up to ~25 days after the burst. Long term radio/mm observations reveal the complex nature of the afterglow, which does not conform to the predictions of the standard afterglow model. We find that the microphysical parameters of the external forward shock, representing the share of shock-created energy in the non-thermal electron population and magnetic field, are evolving with time. The kinetic energy in the blast-wave is almost an order of magnitude higher than that measured in the prompt emission.

GRB 190114C is the first GRB for which the detection of very-high energy emission up to the TeV range has been reported. It is still unclear whether environmental properties might have contributed to the production of these very high-energy photons, or if it is solely related to the released GRB emission. The relatively low redshift of the GRB (z=0.425) allows us to study the host galaxy of this event in detail, and to potentially identify idiosyncrasies that could point to progenitor characteristics or environmental properties responsible for such a unique event. We use ultraviolet, optical, infrared and submillimetre imaging and spectroscopy obtained with HST, VLT and ALMA to obtain an extensive dataset on which the analysis of the host galaxy is based. The host system is composed of a close pair of interacting galaxies (Delta v = 50 km s^-1), both of which are well-detected by ALMA in CO(3-2). The GRB occurred within the nuclear region (~170 pc from the centre) of the less massive but more star-forming galaxy of the pair. The host is more massive (log(M/M_odot)=9.3) than average GRB hosts at that redshift and the location of the GRB is rather unique. The enhanced star-formation rate was probably triggered by tidal interactions between the two galaxies. Our ALMA observations indicate that both host galaxy and companion have a high molecular gas fraction, as has been observed before in interacting galaxy pairs. The location of the GRB within the core of an interacting galaxy with an extinguished line-of-sight is indicative of a denser environment than typically observed for GRBs and could have been crucial for the generation of the very-high-energy photons that were observed.

We present a new measurement of the evolving galaxy far-IR luminosity function (LF) extending out to redshifts z~5, with resulting implications for the level of dust-obscured star-formation density in the young Universe. To achieve this we have exploited recent advances in sub-mm/mm imaging with SCUBA-2 on the James Clerk Maxwell Telescope (JCMT) and the Atacama Large Millimeter/Submillimeter Array (ALMA), which together provide unconfused imaging with sufficient dynamic range to provide meaningful coverage of the luminosity-redshift plane out to z>4. Our results support previous indications that the faint-end slope of the far-IR LF is sufficiently flat that comoving luminosity-density is dominated by bright objects (~L*). However, we find that the number-density/luminosity of such sources at high redshifts has been severely over-estimated by studies that have attempted to push the highly-confused Herschel SPIRE surveys beyond z~2. Consequently we confirm recent reports that cosmic star-formation density is dominated by UV-visible star formation at z>4. Using both direct (1/Vmax) and maximum likelihood determinations of the LF, we find that its high-redshift evolution is well characterized by continued positive luminosity evolution coupled with negative density evolution (with increasing redshift). This explains why bright sub-mm sources continue to be found at z>5, even though their integrated contribution to cosmic star-formation density at such early times is very small. The evolution of the far-IR galaxy LF thus appears similar in form to that already established for active galactic nuclei, possibly reflecting a similar dependence on the growth of galaxy mass.

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.

We present SCUBA-2 850-$\mu$m observations of 13 candidate starbursting protoclusters selected using Planck and Herschel data. The cumulative number counts of the 850-$\mu$m sources in 9/13 of these candidate protoclusters show significant overdensities compared to the field, with the probability $<$10$^{-2}$ assuming the sources are randomly distributed in the sky. Using the 250-, 350-, 500- and 850-$\mu$m flux densities, we estimate the photometric redshifts of individual SCUBA-2 sources by fitting spectral energy distribution (SED) templates with an MCMC method. The photometric redshift distribution, peaking at $2<z<3$, is consistent with that of known $z>2$ protoclusters and the peak of the cosmic star-formation rate density (SFRD). We find that the 850-$\mu$m sources in our candidate protoclusters have infrared luminosities of $L_{\mathrm{IR}}\gtrsim$10$^{12}L_{\odot}$ and star-formation rates of SFR=(500-1,500)$M_{\odot}$yr$^{-1}$. By comparing with results in the literature considering only Herschel photometry, we conclude that our 13 candidate protoclusters can be categorised into four groups: six of them being high-redshift starbursting protoclusters, one being a lower-redshift cluster/protocluster, three being protoclusters that contain lensed DSFG(s) or are rich in 850-$\mu$m sources, and three regions without significant Herschel or SCUBA-2 source overdensities. The total SFRs of the candidate protoclusters are found to be comparable or higher than those of known protoclusters, suggesting our sample contains some of the most extreme protocluster population. We infer that cross-matching Planck and Herschel data is a robust method for selecting candidate protoclusters with overdensities of 850-$\mu$m sources.

Unusual stellar explosions represent an opportunity to learn about both stellar and galaxy evolution. Mapping the atomic gas in host galaxies of such transients can lead to an understanding of the conditions triggering them. We provide resolved atomic gas observations of the host galaxy, CGCG137-068, of the unusual, poorly-understood transient AT2018cow searching for clues to understand its nature. We test whether it is consistent with a recent inflow of atomic gas from the intergalactic medium, as suggested for host galaxies of gamma-ray bursts (GRBs) and some supernovae (SNe). We observed the HI hyperfine structure line of the AT2018cow host with the Giant Metrewave Radio Telescope. There is no unusual atomic gas concentration near the position of AT2018cow. The gas distribution is much more regular than those of GRB/SN hosts. The AT2018cow host has an atomic gas mass lower by 0.24 dex than predicted from its star formation rate (SFR) and is at the lower edge of the galaxy main sequence. In the continuum we detected the emission of AT2018cow and of a star-forming region in the north-eastern part of the bar (away from AT2018cow). This region hosts a third of the galaxy's SFR. The absence of atomic gas concentration close to AT2018cow, along with a normal SFR and regular HI velocity field, sets CGCG137-068 apart from GRB/SN hosts studied in HI. The environment of AT2018cow therefore suggests that its progenitor may not have been a massive star. Our findings are consistent with an origin of the transient that does not require a connection between its progenitor and gas concentration or inflow: an exploding low-mass star, a tidal disruption event, a merger of white dwarfs, or a merger between a neutron star and a giant star. We interpret the recently reported atomic gas ring in CGCG137-068 as a result of internal processes connected with gravitational resonances caused by the bar.

We combine high-resolution ALMA and HST/CANDELS observations of 20 submillimeter galaxies (SMGs) predominantly from the AS2UDS survey at z~2 with bright rest-frame optical counterparts (Ks < 22.9) to investigate the resolved structural properties of their dust and stellar components. We derive two-dimensional stellar-mass distributions that are inferred from spatial mass-to-light ratio (M/L) corrections based on rest-frame optical colors. Due to the high central column densities of dust in our SMGs, our mass distributions likely represent a lower limit to the true central mass density. The centroid positions between the inferred stellar-mass and the dust distributions agree within 1.1 kpc, indicating an overall good spatial agreement between the two components. The majority of our sources exhibit compact dust configurations relative to the stellar component (with a median ratio of effective radii Re,dust/Re,Mstar = 0.6). This ratio does not change with specific star-formation rate (sSFR) over the factor of 30 spanned by our targets, sampling the locus of "normal" main sequence galaxies up to the starburst regime, log(sSFR/sSFRMS) > 0.5. Our results imply that massive SMGs are experiencing centrally enhanced star formation unlike typical spiral galaxies in the local Universe. The sizes and stellar densities of our SMGs are in agreement with those of the passive population at z=1.5, consistent with these systems being the descendants of z~2 SMGs.

The mechanism of dust formation in galaxies at high redshift is still unknown. Asymptotic giant branch (AGB) stars and explosions of supernovae (SNe) are possible dust producers, and non-stellar processes may substantially contribute to dust production, for example grain growth in the interstellar medium (ISM). Our aim is to determine the contribution to dust production of AGB stars and SNe in nine galaxies at z ~ 6-8.3, for which observations of dust have been recently attempted. In order to determine the origin of the observed dust we have determined dust yields per AGB star and SN required to explain the total amounts of dust in these galaxies. We find that AGB stars were not able to produce the amounts of dust observed in the galaxies in our sample. In order to explain these dust masses, SNe would have to have maximum efficiency and not destroy the dust which they formed. Therefore, the observed amounts of dust in the galaxies in the early universe were formed either by efficient supernovae or by a non-stellar mechanism, for instance the grain growth in the interstellar medium.

Using ALMA, we report high angular-resolution observations of the redshift z=3.63 galaxy, G09v1.97, one of the most luminous strongly lensed galaxies discovered by the H-ATLAS survey. We present 0"2-0"4 resolution images of the rest-frame 188 and 419$\mu$m dust continuum and the CO(6-5), H2O(211-202) and J=2 H2O+ line emission. We also report the detection of H$_2^{18}$O in this source. The dust continuum and molecular gas emission are resolved into a nearly complete ~1"5 diameter Einstein ring plus a weaker image in the center, which is caused by a special dual deflector lensing configuration. The observed line profiles of the CO, H2O and H2O+ lines are strikingly similar. In the source plane, we reconstruct the dust continuum images and the spectral cubes of the line emission at sub-kpc scales. The reconstructed dust emission in the source plane is dominated by a compact disk with an effective radius of 0.7kpc plus an overlapping extended disk with a radius twice as large. While the average magnification for the dust continuum is $\mu$~10-11, the magnification of the line emission varies 5 to 22 across different velocity components. The emission lines have similar spatial and kinematic distributions. The molecular gas and dust content reveal that G09v1.97 is a gas-rich major merger in its pre-coalescence phase. Both of the merging companions are intrinsically ULIRGs with LIR reaching $\gtrsim 4\times10^{12}L_\odot$, and the total LIR of G09v1.97 is $1.4\times10^{13}L_\odot$. The approaching southern galaxy shows no obvious kinematic structure with a semi-major half-light radius a_s=0.4kpc, while the receding galaxy resembles an a_s=1.2kpc rotating disk. The two galaxies are separated by a projected distance of 1.3kpc, bridged by weak line emission that is co-spatially located with the cold-dust-emission peak, suggesting a large amount of cold ISM in the interacting region. (abridged)

We present the catalogue and properties of sources in AS2UDS, an 870-$\mu$m continuum survey with the Atacama Large Millimetre/sub-millimetre Array (ALMA) of 716 single-dish sub-millimetre sources detected in the UKIDSS/UDS field by the SCUBA-2 Cosmology Legacy Survey. In our sensitive ALMA follow-up observations we detect 708 sub-millimetre galaxies (SMGs) at $>$\,4.3$\sigma$ significance across the $\sim$\,1-degree diameter field. We combine our precise ALMA positions with the extensive multi-wavelength coverage in the UDS field to fit the spectral energy distributions of our SMGs to derive a median redshift of $z_{\rm phot}=$\,2.61$\pm$0.09. This large sample reveals a statistically significant trend of increasing sub-millimetre flux with redshift suggestive of galaxy downsizing. 101 ALMA maps do not show a $>$\,4.3$\sigma$ SMG, but we demonstrate from stacking {\it Herschel} SPIRE observations at these positions, that the vast majority of these blank maps correspond to real single-dish sub-millimetre sources. We further show that these blank maps contain an excess of galaxies at $z_{\rm phot}=$\,1.5--4 compared to random fields, similar to the redshift range of the ALMA-detected SMGs. In addition, we combine X-ray and mid-infrared active galaxy nuclei activity (AGN) indicators to yield a likely range for the AGN fraction of 8--28\,\% in our sample. Finally, we compare the redshifts of this population of high-redshift, strongly star-forming galaxies with the inferred formation redshifts of massive, passive galaxies being found out to $z\sim$\,2, finding reasonable agreement -- in support of an evolutionary connection between these two classes of massive galaxy.

Long gamma-ray bursts mark the death of massive stars, as revealed by their association with energetic broad-lined stripped-envelope supernovae. The scarcity of nearby events and the brightness of the GRB afterglow, dominating the first days of emission, have so far prevented the study of the very early stages of the GRB-SN evolution. Here we present detailed, multi-epoch spectroscopic observations of SN 2017iuk, associated with GRB 171205A which display features at extremely high expansion velocities of $\sim$ 100,000 km s$^{-1}$ within the first day after the burst. These high-velocity components are characterized by chemical abundances different from those observed in the ejecta of SN 2017iuk at later times. Using spectral synthesis models developed for the SN 2017iuk, we explain these early features as originating not from the supernova ejecta, but from a hot cocoon generated by the energy injection of a mildly-relativistic GRB jet expanding into the medium surrounding the progenitor star. This cocoon becomes rapidly transparent and is outshone by the supernova emission which starts dominating three days after the burst. These results proves that the jet plays an important role not only in powering the GRB event but also its associated supernova.

We present measurements of the angular correlation function of sub-millimeter (sub-mm) galaxies (SMGs) identified in four out of the five fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) - GAMA-9h, GAMA-12h, GAMA-15h and NGP - with flux densities $S_{250\mu m}$>30 mJy at 250 {\mu}m. We show that galaxies selected at this wavelength trace the underlying matter distribution differently at low and high redshifts. We study the evolution of the clustering finding that at low redshifts sub-mm galaxies exhibit clustering strengths of $r_0$ $\sim$ 2 - 3 $h^{-1}$ Mpc, below z < 0.3. At high redshifts, on the other hand, we find that sub-mm galaxies are more strongly clustered with correlation lengths $r_0$ = 8.1 $\pm$ 0.5, 8.8 $\pm$ 0.8 and 13.9 $\pm$ 3.9 $h^{-1}$Mpc at z = 1 - 2, 2 - 3 and 3 - 5, respectively. We show that sub-mm galaxies across the redshift range 1 < z < 5, typically reside in dark-matter halos of mass of the order of ~ $10^{12.5}$ - $10^{13.0}$ $h^{-1} \, M_{\odot}$ and are consistent with being the progenitors of local massive elliptical galaxies that we see in the local Universe.

Galaxy interactions are often accompanied by an enhanced star formation rate (SFR). Since molecular gas is essential for star formation, it is vital to establish whether, and by how much, galaxy interactions affect the molecular gas properties. We investigate the effect of interactions on global molecular gas properties by studying a sample of 58 galaxies in pairs and 154 control galaxies. Molecular gas properties are determined from observations with the JCMT, PMO, CSO telescopes, and supplemented with data from the xCOLD GASS and JINGLE surveys at $^{12}$CO(1-0) and $^{12}$CO(2-1). The SFR, gas mass ($M_\mathrm{H_{2}}$), and gas fraction ($f_{gas}$) are all enhanced in galaxies in pairs by $\sim$ 2.5 times compared to the controls matched in redshift, mass, and effective radius, while the enhancement of star formation efficiency (SFE $\equiv$ SFR/$M_{H_{2}}$) is less than a factor of 2. We also find that the enhancements in SFR, $M_{H_{2}}$ and $f_{gas}$ increase with decreasing pair separation and are larger in systems with smaller stellar mass ratio. Conversely, the SFE is only enhanced in close pairs (separation $<$ 20 kpc) and equal-mass systems; therefore most galaxies in pairs lie in the same parameter space on the SFR-$M_{H_{2}}$ plane as controls. This is the first time that the dependence of molecular gas properties on merger configurations is probed statistically with a relatively large sample and with a carefully-selected control sample for individual galaxies. We conclude that galaxy interactions do modify the molecular gas properties, although the strength of the effect is merger configuration dependent.

Ultra-faint galaxies are hosted by small dark matter halos with shallow gravitational potential wells, hence their star formation activity is more sensitive to feedback effects. The shape of the faint-end of the high-$z$ galaxy luminosity function (LF) contains important information on star formation and its interaction with the reionization process during the Epoch of Reionization (EoR). High-$z$ galaxies with $M_{\rm UV}\gtrsim-17$ have only recently become accessible thanks to the Frontier Fields (FFs) survey combining deep {\it HST} imaging and the gravitational lensing effect. In this paper we investigate the faint-end of the LF at redshift $>$5 using the data of FFs clusters Abell 2744 (A2744), MACSJ0416.1-2403 (M0416), MACSJ0717.5+3745 (M0717) and MACSJ1149.5+2223 (M1149). We analyze both an empirical and a physically-motivated LF model to obtain constraints on a possible turn-over of LF at faint magnitudes. In the empirical model the LF drops fast when the absolute UV magnitude $M_{\rm UV}$ is much larger than a turn-over absolute UV magnitude $M_{\rm UV}^{\rm T}$. We obtain $M_{\rm UV}^{\rm T}\gtrsim-14.6 $ (15.2) at 1 (2) $\sigma$ confidence level (C.L.) for $z\sim6$. In the physically-motivated analytical model, star formation in halos with circular velocity below $v_c^*$ is fully quenched if these halos are located in ionized regions. Using updated lensing models and new additional FFs data, we re-analyze previous constraints on $v_c^*$ and $f_{\rm esc}$ presented by Castellano et al. 2016a (C16a) using a smaller dataset. We obtain new constraints on $v_c^*\lesssim 59$ km s$^{-1}$ and $f_{\rm esc}\lesssim 56\%$ (both at 2$\sigma$ C.L.) and conclude that there is no turn-over detected so far from the analyzed FFs data. Forthcoming {\it JWST} observations will be key to tight these constraints further.

We present Atacama Large Millimeter/submillimeter Array (ALMA) resolved observations of molecular gas in galaxies up to $z=0.35$ to characterise the role of global galactic dynamics on the global interstellar medium (ISM) properties. These observations consist of a sub-sample of 39 galaxies taken from the Valpara\'iso ALMA Line Emission Survey (VALES). From the CO($J=1-0)$ emission line, we quantify the kinematic parameters by modelling the velocity fields. We find that the IR luminosity increases with the rotational to dispersion velocity ratio ($V_{\rm rot}/\sigma_v$, corrected for inclination). We find a dependence between $V_{\rm rot}/\sigma_v$ and the [CII]/IR ratio, suggesting that the so-called `[CII] deficit' is related to the dynamical state of the galaxies. We find that global pressure support is needed to reconcile the dynamical mass estimates with the stellar masses in our systems with low $V_{\rm rot}/\sigma_v$ values. The star formation rate (SFR) is weakly correlated with the molecular gas fraction ($f_{\rm H_2}$) in our sample, suggesting that the release of gravitational energy from cold gas may not be the main energy source of the turbulent motions seen in the VALES galaxies. By defining a proxy of the `star formation efficiency' parameter as the SFR divided by the CO luminosity (SFE$'\equiv$ SFR/L$'_{\rm CO}$), we find a constant SFE$'$ per crossing time ($t_{\rm cross}$). We suggest that $t_{\rm cross}$ may be the controlling timescale in which the star formation occurs in dusty $z\sim0.03-0.35$ galaxies.

We present structural parameters and morphological properties of faint 450-um selected submillimeter galaxies (SMGs) from the JCMT Large Program, STUDIES, in the COSMOS-CANDELS region. Their properties are compared to an 850um selected and a matched star-forming samples. We investigate stellar structures of 169 faint 450-um sources (S450=2.8-29.6mJy; S/N>4) at z<3 using HST near-infrared observations. Based on our spectral energy distribution fitting, half of such faint SMGs (LIR=10^11.65+-0.98Lsun) lie above the star-formation rate (SFR)/stellar mass plane. The size-mass relation shows that these SMGs are generally similar to less-luminous star-forming galaxies selected by NUV-r vs. r-J colors. Because of the intrinsic luminosity of the sample, their rest-frame optical emission is less extended than the 850um sources (S850>2mJy), and more extended than the star-forming galaxies in the same redshift range. For the stellar mass and SFR matched sample at z~=1 and z~=2, the size differences are marginal between faint SMGs and the matched galaxies. Moreover, faint SMGs have similar Sersic indices and projected axis ratios as star-forming galaxies with the same stellar mass and SFR. Both SMGs and the matched galaxies show high fractions (~70%) of disturbed features at z~=2, and the fractions depend on the SFRs. These suggest that their star formation activity is related to galaxy merging, and the stellar structures of SMGs are similar to those of star-forming galaxies. We show that the depths of submillimeter surveys are approaching the lower luminosity end of star-forming galaxies, allowing us to detect galaxies on the main sequence.

The potential similarity of the powering mechanisms of relativistic SNe and GRBs allowed us to make a prediction that relativistic SNe are born in environments similar to those of GRBs, that is, ones which are rich in atomic gas. Here we embark on testing this hypothesis by analysing the properties of the host galaxy NGC 3278 of the relativistic SN 2009bb. This is the first time the atomic gas properties of a relativistic SN host are provided and the first time resolved 21 cm-hydrogen-line (HI) information is provided for a host of an SN of any type in the context of the SN position. We obtained radio observations with ATCA covering the HI line, and optical integral field unit spectroscopy observations with MUSE. The atomic gas distribution of NGC 3278 is not centred on the optical galaxy centre, but instead around a third of atomic gas resides in the region close to the SN position. This galaxy has a few times lower atomic and molecular gas masses than predicted from its SFR. SN 2009bb exploded close to the region with the highest SFR density and the lowest age (~5.5 Myr, corresponding to the initial mass of the progenitor star ~36 Mo). As for GRB hosts, the gas properties of NGC 3278 are consistent with a recent inflow of gas from the intergalactic medium, which explains the concentration of atomic gas close to the SN position and the enhanced SFR. Super-solar metallicity at the position of the SN (unlike for most GRBs) may mean that relativistic explosions signal a recent inflow of gas (and subsequent star formation), and their type (GRBs or SNe) is determined either i) by the metallicity of the inflowing gas, so that metal-poor gas results in a GRB explosion and metal-rich gas results in a relativistic SN explosion without an accompanying GRB, or ii) by the efficiency of gas mixing, or iii) by the type of the galaxy.

We describe the application of the supervised machine-learning algorithms to identify the likely multi-wavelength counterparts to submillimeter sources detected in panoramic, single-dish submillimeter surveys. As a training set, we employ a sample of 695 ($S_{\rm 870\mu m}$ >1 mJy) submillimeter galaxies (SMGs) with precise identifications from the ALMA follow-up of the SCUBA-2 Cosmology Legacy Survey's UKIDSS-UDS field (AS2UDS). We show that radio emission, near-/mid-infrared colors, photometric redshift, and absolute $H$-band magnitude are effective predictors that can distinguish SMGs from submillimeter-faint field galaxies. Our combined radio+machine-learning method is able to successfully recover $\sim$85 percent of ALMA-identified SMGs which are detected in at least three bands from the ultraviolet to radio. We confirm the robustness of our method by dividing our training set into independent subsets and using these for training and testing respectively, as well as applying our method to an independent sample of $\sim$100 ALMA-identified SMGs from the ALMA/LABOCA ECDF-South Survey (ALESS). To further test our methodology, we stack the 870$\mu$m ALMA maps at the positions of those $K$-band galaxies that are classified as SMG counterparts by the machine-learning but do not have a $>$4.3$ \sigma$ ALMA detection. The median peak flux density of these galaxies is $S_{\rm 870\mu m}=(0.61\pm0.03)$ mJy, demonstrating that our method can recover faint and/or diffuse SMGs even when they are below the detection threshold of our ALMA observations. In future, we will apply this method to samples drawn from panoramic single-dish submillimeter surveys which currently lack interferometric follow-up observations, to address science questions which can only be tackled with large, statistical samples of SMGs.

We report the detection of the radio afterglow of a long gamma-ray burst (GRB) 111005A at 5-345 GHz, including the very long baseline interferometry observations with the positional error of 0.2 mas. The afterglow position is coincident with the disk of a galaxy ESO 580-49 at z= 0.01326 (~1" from its center), which makes GRB 111005A the second closest GRB known to date, after GRB 980425. The radio afterglow of GRB 111005A was an order of magnitude less luminous than those of local low-luminosity GRBs, and obviously than those of cosmological GRBs. The radio flux was approximately constant and then experienced an unusually rapid decay a month after the GRB explosion. Similarly to only two other GRBs, we did not find the associated supernovae (SN), despite deep near- and mid-infrared observations 1-9 days after the GRB explosion, reaching ~20 times fainter than other SNe associated with GRBs. Moreover, we measured twice solar metallicity for the GRB location. The low gamma-ray and radio luminosities, rapid decay, lack of a SN, and super-solar metallicity suggest that GRB 111005A represents a different rare class of GRBs than typical core-collapse events. We modelled the spectral energy distribution of the GRB 111005A host finding that it is a dwarf, moderately star-forming galaxy, similar to the host of GRB 980425. The existence of two local GRBs in such galaxies is still consistent with the hypothesis that the GRB rate is proportional to the cosmic star formation rate (SFR) density, but suggests that the GRB rate is biased towards low SFRs. Using the far-infrared detection of ESO 580-49, we conclude that the hosts of both GRBs 111005A and 980425 exhibit lower dust content than what would be expected from their stellar masses and optical colours.

We report the first results of AS2UDS: an 870 $\mu$m continuum survey with the Atacama Large Millimeter/Submillimeter Array (ALMA) of a total area of $\sim$ 50 arcmin$^2$ comprising a complete sample of 716 submillimeter sources drawn from the SCUBA-2 Cosmology Legacy Survey (S2CLS) map of the UKIDSS/UDS field. The S2CLS parent sample covers a 0.96 degree$^2$ field at $\sigma_{850}=0.90\pm0.05$ mJy beam$^{-1}$. Our deep, high-resolution ALMA observations with $\sigma_{\rm 870}\sim$ 0.25 mJy and a 0.15"--0.30" FWHM synthesized beam, provide precise locations for 695 submillimetre galaxies (SMGs) responsible for the submillimeter emission corresponding to 606 sources in the low resolution, single-dish map. We measure the number counts of SMGs brighter than $S_{\rm 870}\geq$ 4 mJy, free from the effects of blending and show that the normalisation of the counts falls by 28 $\pm$ 2% in comparison to the SCUBA-2 parent sample, but that the shape remains unchanged. We determine that 44$^{+16}_{-14}$% of the brighter single-dish sources with $S_{850}\geq$ 9 mJy consist of a blend of two or more ALMA-detectable SMGs brighter than $S_{870}\sim$ 1 mJy (corresponding to a galaxy with a total-infrared luminosity of $L_{IR}\geq$ 10$^{12}$ L$_\odot$), in comparison to 28 $\pm$ 2% for the single-dish sources at $S_{\rm 850}\geq$ 5 mJy. Using the 46 single-dish submillimeter sources that contain two or more ALMA-detected SMGs with photometric redshifts, we show that there is a significant statistical excess of pairs of SMGs with similar redshifts ($<$1% probability of occurring by chance), suggesting that at least 30% of these blends arise from physically associated pairs of SMGs.

Dust-enshrouded, starbursting, submillimeter galaxies (SMGs) at $z \geq 3$ have been proposed as progenitors of $z \geq 2$ compact quiescent galaxies (cQGs). To test this connection, we present a detailed spatially resolved study of the stars, dust and stellar mass in a sample of six submillimeter-bright starburst galaxies at $z \sim 4.5$. The stellar UV emission probed by HST is extended, irregular and shows evidence of multiple components. Informed by HST, we deblend Spitzer/IRAC data at rest-frame optical finding that the systems are undergoing minor mergers, with a typical stellar mass ratio of 1:6.5. The FIR dust continuum emission traced by ALMA locates the bulk of star formation in extremely compact regions (median $r_{\rm{e}} = 0.70 \pm 0.29$ kpc) and it is in all cases associated with the most massive component of the mergers (median $\log (M_{*}/M_{\odot}) = 10.49 \pm 0.32$). We compare spatially resolved UV slope ($\beta$) maps with the FIR dust continuum to study the infrared excess ($\rm{IRX} = L_{\rm{IR}}/L_{\rm{UV}}$)-$\beta$ relation. The SMGs display systematically higher $\rm{IRX}$ values than expected from the nominal trend, demonstrating that the FIR and UV emissions are spatially disconnected. Finally, we show that the SMGs fall on the mass-size plane at smaller stellar masses and sizes than cQGs at $z = 2$. Taking into account the expected evolution in stellar mass and size between $z = 4.5$ and $z = 2$ due to the ongoing starburst and mergers with minor companions, this is in agreement with a direct evolutionary connection between the two populations.

Interferometric observations have demonstrated that a significant fraction of single-dish submillimetre (submm) sources are blends of multiple submm galaxies (SMGs), but the nature of this multiplicity, i.e. whether the galaxies are physically associated or chance projections, has not been determined. We performed spectroscopy of 11 SMGs in six multi-component submm sources, obtaining spectroscopic redshifts for nine of them. For an additional two component SMGs, we detected continuum emission but no obvious features. We supplement our observed sources with four sources from the literature. This sample allows us to statistically constrain the physical nature of single-dish submm source multiplicity for the first time. In three [3/7, or 43 (-33/+39) per cent at 95% confidence] of the single-dish sources for which the nature of the blending is unambiguous, the components for which spectroscopic redshifts are available are physically associated, whereas 4/7 [57 (-39/+33) per cent] have at least one unassociated component. When components whose spectra exhibit continuum but no features and for which the photometric redshift is significantly different from the spectroscopic redshift of the other component are also considered, 6/9 [67 (-37/+26) per cent] of the single-dish sources are comprised of at least one unassociated component. The nature of the multiplicity of one single-dish source is ambiguous. We conclude that physically associated systems and chance projections both contribute to the multi-component single-dish submm source population. This result contradicts the conventional wisdom that bright submm sources are solely a result of merger-induced starbursts, as blending of unassociated galaxies is also important.

We present a multi-wavelength analysis of galaxies selected at 450 and 850um from the deepest SCUBA-2 observations in the Extended Groth Strip (EGS) field, which have an average depth of sigma_450=1.9 and sigma_850=0.46 mJy/beam over ~70 sq. arcmin. The final sample comprises 95 sources: 56 (59%) are detected at both wavelengths, 31 (33%) are detected only at 850um, and 8 (8%) are detected only at 450um. We identify counterparts for 75% of the whole sample. The redshift distributions of the 450 and 850um samples peak at different redshifts with median values of z=1.66 +/- 0.18 and z=2.30 +/- 0.20, respectively. However, the two populations have similar IR luminosities, SFRs, and stellar masses, with mean values of (1.5 +/- 0.2) x 10^12 L_sun, (150 +/- 20) M_sun/yr, and (9.0 +/- 0.6) x 10^10 M_sun, respectively. This places most of our sources (>85%) on the high-mass end of the `main-sequence' of star-forming galaxies. Exploring the IR excess vs UV-slope (IRX-beta) relation we find that the most luminous galaxies are consistent with the Meurer law, while the less luminous galaxies lie below this relation. Using the results of a two-dimensional modelling of the HST H_160-band imaging, we derive a median Sersic index of n=1.4 +0.3 -0.1 and a median half-light radius of R_1/2 = 4.8 +/ 0.4 kpc. Based on a visual-like classification in the same band, we find that the dominant component for most of the galaxies at all redshifts is a disk-like structure, although there is a transition from irregular disks to disks with a spheroidal component at z~1.4, which morphologically supports the scenario of SMGs as progenitors of massive elliptical galaxies.

With the advent of wide-area submillimeter surveys, a large number of high-redshift gravitationally lensed dusty star-forming galaxies (DSFGs) has been revealed. Due to the simplicity of the selection criteria for candidate lensed sources in such surveys, identified as those with $S_{500\mu m} > 100$ mJy, uncertainties associated with the modelling of the selection function are expunged. The combination of these attributes makes submillimeter surveys ideal for the study of strong lens statistics. We carried out a pilot study of the lensing statistics of submillimetre-selected sources by making observations with the Atacama Large Millimetre Array (ALMA) of a sample of strongly-lensed sources selected from surveys carried out with the Herschel Space Observatory. We attempted to reproduce the distribution of image separations for the lensed sources using a halo mass function taken from a numerical simulation which contains both dark matter and baryons. We used three different density distributions, one based on analytical fits to the halos formed in the EAGLE simulation and two density distributions (Singular Isothermal Sphere (SIS) and SISSA) that have been used before in lensing studies. We found that we could reproduce the observed distribution with all three density distributions, as long as we imposed an upper mass transition of $\sim$$10^{13} M_{\odot}$ for the SIS and SISSA models, above which we assumed that the density distribution could be represented by an NFW profile. We show that we would need a sample of $\sim$500 lensed sources to distinguish between the density distributions, which is practical given the predicted number of lensed sources in the Herschel surveys.

We perform lens modelling and source reconstruction of Submillimeter Array (SMA) data for a sample of 12 strongly lensed galaxies selected at 500$\mu$m in the Herschel Astrophysical Terahertz Large Area Survey H-ATLAS. A previous analysis of the same dataset used a single S\`ersic profile to model the light distribution of each background galaxy. Here we model the source brightness distribution with an adaptive pixel scale scheme, extended to work in the Fourier visibility space of interferometry. We also present new SMA observations for seven other candidate lensed galaxies from the H-ATLAS sample. Our derived lens model parameters are in general consistent with previous findings. However, our estimated magnification factors, ranging from 3 to 10, are lower. The discrepancies are observed in particular where the reconstructed source hints at the presence of multiple knots of emission. We define an effective radius of the reconstructed sources based on the area in the source plane where emission is detected above 5$\sigma$. We also fit the reconstructed source surface brightness with an elliptical Gaussian model. We derive a median value $r_{eff}\,\sim 1.77\,$kpc and a median Gaussian full width at half maximum $\sim1.47\,$kpc. After correction for magnification, our sources have intrinsic star formation rates SFR$\,\sim900-3500\,M_{\odot}yr^{-1}$, resulting in a median star formation rate surface density $\Sigma_{SFR}\sim132\,M_{\odot}$ yr$^{-1}$ kpc$^{-2}$ (or $\sim 218\,M_{\odot}$ yr$^{-1}$ kpc$^{-2}$ for the Gaussian fit). This is consistent with what observed for other star forming galaxies at similar redshifts, and is significantly below the Eddington limit for a radiation pressure regulated starburst.

We use a sample of 4178 Lyman break galaxies (LBGs) at z = 3, 4 and 5 in the UKIRT Infrared Deep Sky Survey (UKIDSS) Ultra Deep Survey (UDS) field to investigate the relationship between the observed slope of the stellar continuum emission in the ultraviolet, {\beta}, and the thermal dust emission, as quantified via the so-called 'infrared excess' (IRX = LIR/LUV). Through a stacking analysis we directly measure the 850-{\mu}m flux density of LBGs in our deep (0.9mJy) James Clerk Maxwell Telescope (JCMT) SCUBA-2 850-{\mu}m map, as well as deep public Herschel/SPIRE 250-, 350- and 500-{\mu}m imaging. We establish functional forms for the IRX-{\beta} relation to z ~ 5, confirming that there is no significant redshift evolution of the relation and that the resulting average IRX-{\beta} curve is consistent with a Calzetti-like attenuation law. We compare our results with recent work in the literature, finding that discrepancies in the slope of the IRX-{\beta} relation are driven by biases in the methodology used to determine the ultraviolet slopes. Consistent results are found when IRX-{\beta} is evaluated by stacking in bins of stellar mass, M, and we argue that the near-linear IRX-M relationship is a better proxy for correcting observed UV luminosities to total star formation rates, provided an accurate handle on M can be had, and also gives clues as to the physical driver of the role of dust-obscured star formation in high-redshift galaxies.

By determining the nature of all the Planck compact sources within 808.4 deg^2 of large Herschel surveys, we have identified 27 candidate proto-clusters of dusty star forming galaxies (DSFGs) that are at least 3{\sigma} overdense in either 250, 350 or 500 $\mu$mm sources. We find roughly half of all the Planck compact sources are resolved by Herschel into multiple discrete objects, with the other half remaining unresolved by Herschel. We find a significant difference between versions of the Planck catalogues, with earlier releases hosting a larger fraction of candidate proto-clusters and Galactic Cirrus than later releases, which we ascribe to a difference in the filters used in the creation of the three catalogues. We find a surface density of DSFG candidate proto-clusters of (3.3 $\pm$ 0.7) x 10^(-2) sources deg^(-2), in good agreement with previous similar studies. We find that a Planck colour selection of S_{857}/S_{545} < 2 works well to select candidate proto-clusters, but can miss proto-clusters at z < 2. The Herschel colours of individual candidate proto-cluster members indicate our candidate proto-clusters all likely all lie at z > 1. Our candidate proto-clusters are a factor of 5 times brighter at 353 GHz than expected from simulations, even in the most conservative estimates. Further observations are needed to confirm whether these candidate proto-clusters are physical clusters, multiple proto-clusters along the line of sight, or chance alignments of unassociated sources.

Binary neutron-star mergers (BNSMs) are among the most readily detectable gravitational-wave (GW) sources with LIGO. They are also thought to produce short $\gamma$-ray bursts (SGRBs), and kilonovae that are powered by r-process nuclei. Detecting these phenomena simultaneously would provide an unprecedented view of the physics during and after the merger of two compact objects. Such a Rosetta Stone event was detected by LIGO/Virgo on 17 August 2017 at a distance of $\sim 44$ Mpc. We monitored the position of the BNSM with ALMA at 338.5 GHz and GMRT at 1.4 GHz, from 1.4 to 44 days after the merger. Our observations rule out any afterglow more luminous than $3\times 10^{26}~{\rm erg\,s}^{-1}\,{\rm Hz}^{-1}$ in these bands, probing $>$2--4 dex fainter than previous SGRB limits. We match these limits, in conjunction with public data announcing the appearance of X-ray and radio emission in the weeks after the GW event, to templates of off-axis afterglows. Our broadband modeling suggests that GW170817 was accompanied by a SGRB and that the GRB jet, powered by $E_{\rm AG,\,iso}\sim10^{50}$~erg, had a half-opening angle of $\sim20^\circ$, and was misaligned by $\sim41^\circ$ from our line of sight. The data are also consistent with a more collimated jet: $E_{\rm AG,\,iso}\sim10^{51}$~erg, $\theta_{1/2,\,\rm jet}\sim5^\circ$, $\theta_{\rm obs}\sim17^\circ$. This is the most conclusive detection of an off-axis GRB afterglow and the first associated with a BNSM-GW event to date. Assuming a uniform top-hat jet, we use the viewing angle estimates to infer the initial bulk Lorentz factor and true energy release of the burst.

(abridged) We present the IRAM-30m observations of multiple-J CO and CI line emission in a sample of redshift ~2-4 Herschel-ATLAS SMGs. A non-negligible effect of differential lensing is found for the CO emission lines, which could have caused significant underestimations of the linewidths, hence of the dynamical masses. The CO SLEDs are found to be similar to those of the local starburst-dominated ULIRGs and of the previously studied SMGs. After correcting for lensing amplification, we derived the global properties of the bulk of molecular gas in the SMGs using non-LTE radiative transfer modelling. The gas thermal pressure is found to be correlated with star formation efficiency. Further decomposing the CO SLEDs into two excitation components, we find a low-excitation component, which is less correlated with star formation, and a high-excitation one which is tightly related to the on-going star-forming activity. Additionally, tight linear correlations between the FIR and CO line luminosities have been confirmed for the $J \ge 5$ CO lines, implying that these CO lines are good tracers of star formation. The [CI](2-1) lines follow the tight linear correlation between the luminosities of the [CI](2-1) and the CO(1-0) line found in local starbursts, indicating that CI lines could serve as good total molecular gas mass tracers for high-redshift SMGs. The total mass of the molecular gas reservoir, $(1-30) \times 10^{10} M_\odot$, suggests a typical molecular gas depletion time ~20-100 Myr and a gas to dust mass ratio ${\delta}_{\rm GDR}$~30-100. The ratio between CO line luminosity and the dust mass appears to be slowly increasing with redshift for the SMGs, which need to be further confirmed. Finally, through comparing the linewidth of CO and H2O lines, we find that they agree well in almost all our SMGs, confirming that the emitting regions are co-spatially located.

We present physical properties [redshifts (z), star-formation rates (SFRs) and stellar masses (Mstar)] of bright (S850>4mJy) submm galaxies in the ~2deg2 COSMOS and UDS fields selected with SCUBA-2/JCMT. We complete the galaxy identification process for all (~2000) S/N>3.5 850-um sources, but focus our scientific analysis on a high-quality sub-sample of 651 S/N>4 sources with complete multi-wavelength coverage including 1.1-mm imaging. We check the reliability of our identifications, and the robustness of the SCUBA-2 fluxes by revisiting the recent ALMA follow-up of 29 sources in our sample. Considering >4mJy ALMA sources, our identification method has a completeness of ~86 per cent with a reliability of ~92 per cent, and only ~15-20 per cent of sources are significantly affected by multiplicity (when a secondary component contributes >1/3 of the primary source flux). The impact of source blending on the 850-um source counts as determined with SCUBA-2 is modest; scaling the single-dish fluxes by ~0.9 reproduces the ALMA source counts. For our final SCUBA-2 sample we find median z=2.40+0.10-0.04, SFR=287+-6Moyr-1, and log(Mstar/Mo)=11.12+-0.02 (the latter for 349/651 sources with optical identifications). These properties clearly locate bright submm galaxies on the high-mass end of the 'main sequence' of star-forming galaxies out to z~6, suggesting that major mergers are not a dominant driver of the high-redshift submm-selected population. Their number densities are also consistent with the evolving galaxy stellar mass function. Hence, the submm galaxy population is as expected, albeit reproducing the evolution of the main sequence of star-forming galaxies remains a challenge for theoretical models/simulations.

We present stacked average far-infrared spectra of a sample of 197 dusty, star-forming galaxies (DSFGs) at $0.005 < z < 4$ using close to 90% of the SPIRE Fourier Transform Spectrometer (FTS) extragalactic data archive from the Herschel Space Observatory based on 3.5 years of science operations. These spectra explore an observed-frame $\rm 447\,GHz-1568\,GHz$ ($\rm 191\,\mu m-671\,\mu m$) frequency (wavelength) range allowing us to observe the main atomic and molecular lines emitted by gas in the interstellar medium. The sample is sub-divided into five redshift bins at $0.005 < z < 0.05$, $0.05 < z < 0.2$, $0.2 < z < 0.5$, $0.8 < z <2$, and $2 < z < 4$. To study the dependence of observed spectral lines on total infrared luminosity, the sources in a subset of the redshift bins are stacked in luminosity bins. These stacked spectra are used to determine the average properties of the interstellar medium and dense molecular gas properties of DSFGs, in particular, the fine-structure line ([CII] 158 $\mu$m and [OI] 63 $\mu$m) luminosity ratios, and the line to far-IR luminosity ratios are used to model the gas density and radiation field strength in the photodissociation regions (PDRs). For the low-redshift sample, we additionally present the average spectral line energy distributions (SLED) of CO and $\rm{H_2O}$ rotational transitions and also consider PDR conditions based on a combination of [CI] 370 $\mu$m and 609 $\mu$m and $\rm CO (7-6)$ lines. For the high-z ($0.8 < z < 4$) sample PDR models suggest a molecular gas distribution in the presence of a radiation field that is at least a factor of 10$^3$ larger than the Milky-Way and with a neutral gas density of roughly 10$^3$ to 10$^5$ cm$^{-3}$. The corresponding PDR models for the low-z sample suggest a UV radiation field and gas density comparable to those at high-z.

We present high-spatial-resolution ($\sim 0.12''$ or $\approx 800 \, {\rm pc}$ at $z = 4.5$) ALMA $870\,\mu$m dust continuum observations of a sample of 44 ultrared dusty star-forming galaxies (DSFGs) selected from the H-ATLAS and HerMES far-infrared surveys because of their red colors from 250 to 500 $\mu$m: $S_{500} / S_{250} > 1.5$ and $S_{500} / S_{350} > 1.0$. With photometric redshifts in the range $z \sim 4$-6, our sample includes the most luminous starbursting systems in the early Universe known so far, with total obscured star-formation rates (SFRs) of up to $\sim 4,500 \, M_\odot \, {\rm yr}^{-1}$, as well as a population of lensed, less intrinsically luminous sources. The lower limit on the number of ultrared DSFGs at 870 $\mu$m (with flux densities measured from the ALMA maps and thus not affected by source confusion) derived in this work is in reasonable agreement with models of galaxy evolution, whereas there have been reports of conflicts at 500 $\mu$m (where flux densities are derived from SPIRE). Ultrared DSFGs have a variety of morphologies (from relatively extended disks with smooth radial profiles, to compact sources, both isolated and interacting) and an average size, $\theta_{\rm FWHM}$, of $1.46 \pm 0.41\, {\rm kpc}$, considerably smaller than the values reported in previous work for less-luminous DSFGs at lower redshifts. The size and the estimated gas-depletion times of our sources are compatible with their being the progenitors of the most massive, compact, red-and-dead galaxies at $z \sim 2$-3, and ultimately of local ultra-massive elliptical galaxies or massive galaxy clusters. We are witnessing the birth of the high-mass tail of the red sequence of galaxies.

We present the Herschel Bright Sources (HerBS) sample, a sample of bright, high-redshift Herschel sources detected in the 616.4 square degree H-ATLAS survey. The HerBS sample contains 209 galaxies, selected with a 500 {\mu}m flux density greater than 80 mJy and an estimated redshift greater than 2. The sample consists of a combination of HyLIRGs and lensed ULIRGs during the epoch of peak cosmic star formation. In this paper, we present SCUBA-2 observations at 850 ${\mu}$m of 189 galaxies of the HerBS sample, 152 of these sources were detected. We fit a spectral template to the Herschel-SPIRE and 850 ${\mu}$m SCUBA-2 flux densities of 22 sources with spectroscopically determined redshifts, using a two-component modified blackbody spectrum as a template. We find a cold- and hot-dust temperature of 21.29 K and 45.80 K, a cold-to-hot dust mass ratio of 26.62 and a $\beta$ of 1.83. The poor quality of the fit suggests that the sample of galaxies is too diverse to be explained by our simple model. Comparison of our sample to a galaxy evolution model indicates that the fraction of lenses is high. Out of the 152 SCUBA-2 detected galaxies, the model predicts 128.4 $\pm$ 2.1 of those galaxies to be lensed (84.5%). The SPIRE 500 ${\mu}$m flux suggests that out of all 209 HerBS sources, we expect 158.1 $\pm$ 1.7 lensed sources, giving a total lensing fraction of 76 per cent.

Star formation rate density (SFRD) has not been constant throughout the history of the Universe. The rate at which stars form greatly affects the evolution of the Universe, but the factors which drive SFRD evolution remain uncertain. There must be sufficient amount of gas to fuel the star formation, either as a reservoir within a galaxy, or as inflow from the intergalactic medium (IGM). This work explores how the gas accretion rate onto galaxies over time has affected star formation rate. We propose a novel method of measuring cosmic gas accretion rate. This involves comparing the comoving densities of available Hi and H2 gas and the densities of existing stars at different redshifts. We constrained gas accretion until z = 5, and we found that the gas accretion rate density (GARD) is relatively constant in the range from z = 5 to z = 0. This constancy in the GARD is not reflected by the SFRD, which declines significantly between z = 1.0 and z = 0. This work suggests that the decline is not due to a reduction in GARD.

We have combined the wide-area Herschel-ATLAS far-IR survey with spectroscopic redshifts from GAMA and SDSS to define a sample of 21 low--redshift ($z_{\rm spec} < 0.5$) analogs of submm galaxies (SMGs). These have been selected because their dust temperatures and total IR luminosities are similar to those for the classical high-redshift SMG population. As well as presenting the sample, in this paper we report $^{12}$CO(2-1) and $^{12}$CO(1-0) observations of 16 low-redshift analogs of SMGs taken with the IRAM-30m telescope. We have obtained that low-redshift analogs of SMGs represent a very diverse population, similar to what has been found for high-redshift SMGs. A large variety in the molecular gas excitation or $^{12}$CO(2-1)/$^{12}$CO(1-0) line ratio is seen, meaning that extrapolations from $J \geq 2$ CO lines can result in very uncertain molecular gas mass determinations. Our sources with $^{12}$CO(1-0) detections follow the dust--gas correlation found in previous work at different redshifts and luminosities. The molecular gas mass of low-redshift SMGs has an average value of $M_{\rm H_2} \sim 1.6 \times 10^{10}\,M_\odot$ and will be consumed in $\sim 100 \, {\rm Myr}$ . We also find a wide range of molecular gas fractions, with the highest values being compatible with those found in high-redshift SMGs with $^{12}$CO(1-0) detections, which are only the most luminous. Low-redshift SMGs offer a unique opportunity to study the properties of extreme star formation in a detail not possible at higher redshifts.

We present the results of combined deep Keck/NIRC2, HST/WFC3 near-infrared and Herschel far infrared observations of an extremely star forming dusty lensed galaxy identified from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS J133542.9+300401). The galaxy is gravitationally lensed by a massive WISE identified galaxy cluster at $z\sim1$. The lensed galaxy is spectroscopically confirmed at $z=2.685$ from detection of $\rm {CO (1 \rightarrow 0)}$ by GBT and from detection of $\rm {CO (3 \rightarrow 2)}$ obtained with CARMA. We use the combined spectroscopic and imaging observations to construct a detailed lens model of the background dusty star-forming galaxy (DSFG) which allows us to study the source plane properties of the target. The best-fit lens model provide magnification of $\mu_{\rm star}=2.10\pm0.11$ and $\mu_{\rm dust}=2.02\pm0.06$ for the stellar and dust components respectively. Multi-band data yields a magnification corrected star formation rate of $1900(\pm200)\,M_{\odot}{\rm yr^{-1}}$ and stellar mass of $6.8_{-2.7}^{+0.9}\times10^{11}\,M_{\odot}$ consistent with a main sequence of star formation at $z\sim2.6$. The CO observations yield a molecular gas mass of $8.3(\pm1.0)\times10^{10}\,M_{\odot}$, similar to the most massive star-forming galaxies, which together with the high star-formation efficiency are responsible for the intense observed star formation rates. The lensed DSFG has a very short gas depletion time scale of $\sim40$ Myr. The high stellar mass and small gas fractions observed indicate that the lensed DSFG likely has already formed most of its stellar mass and could be a progenitor of the most massive elliptical galaxies found in the local Universe.

We present the calibration between the dust continuum luminosity and interstellar gas content obtained from the Valpara\'{i}so ALMA Line Emission Survey (VALES) sample of 67 main-sequence star-forming galaxies at 0.02<$z$<0.35. We use CO(1-0) observations from the Atacama Large Millimetre/submillimetre Array (ALMA) to trace the molecular gas mass, $M_{\mathrm{H}_{2}}$, and estimate the rest-frame monochromatic luminosity at 850 $\mu$m, $L_{\nu_{850}}$, by extrapolating the dust continuum from MAGPHYS modelling of the far-ultraviolet to submillimetre spectral energy distribution sampled by the Galaxy And Mass Assembly (GAMA) survey. Adopting $\alpha_{\rm CO}$ = 6.5 (K km s$^{-1}$ pc$^{2}$)$^{-1}$, the average ratio of $L_{\nu_{850}}/M_{\mathrm{H}_{2}}$ = (6.4$\pm$1.4)$\times10^{19}$ erg s$^{-1}$ Hz$^{-1}$ $\mathrm{M}_{\odot}^{-1}$, in excellent agreement with literature values. We obtain a linear fit of $\log_{10}$ ($M_{\mathrm{H}_{2}}/\mathrm{M}_{\odot}$) = (0.92$\pm$0.02) $\log_{10}$ ($L_{\nu_{850}}$/erg s$^{-1}$ Hz$^{-1}$)-(17.31$\pm$0.59). We provide relations between $L_{\nu_{850}}$, $M_{\mathrm{H}_{2}}$ and $M_{\mathrm{ISM}}$ when combining the VALES and literature samples, and adopting a Galactic $\alpha_{\rm CO}$ value.

Submillimeter galaxies (SMGs) at $z\gtrsim1$ are luminous in the far-infrared and have star-formation rates, SFR, of hundreds to thousands of solar masses per year. However, it is unclear whether they are true analogs of local ULIRGs or whether the mode of their star formation is more similar to that in local disk galaxies. We target these questions by using Herschel-PACS to examine the conditions in the interstellar medium (ISM) in far-infrared luminous SMGs at z~1-4. We present 70-160 micron photometry and spectroscopy of the [OIV]26 micron, [FeII]26 micron, [SIII]33 micron, [SiII]34 micron, [OIII]52 micron, [NIII]57 micron, and [OI]63 micron fine-structure lines and the S(0) and S(1) hydrogen rotational lines in 13 lensed SMGs identified by their brightness in early Herschel data. Most of the 13 targets are not individually spectroscopically detected and we instead focus on stacking these spectra with observations of an additional 32 SMGs from the \herschel\ archive -- representing a complete compilation of PACS spectroscopy of SMGs. We detect [OI]63 micron, [SiII]34 micron, and [NIII]57 micron at >3sigma in the stacked spectra, determining that the average strengths of these lines relative to the far-IR continuum are $(0.36\pm0.12)\times10^{-3}$, $(0.84\pm0.17)\times10^{-3}$, and $(0.27\pm0.10)\times10^{-3}$, respectively. Using the [OIII]52/[NIII]57 emission line ratio we show that SMGs have average gas-phase metallicities $\gtrsim Z_{\rm sun}$. By using PDR modelling and combining the new spectral measurements with integrated far-infrared fluxes and existing [CII]158 micron data we show that SMGs have average gas densities, n, of $\sim10^{1-3}{\rm cm^{-3}}$ and FUV field strengths, $G_0\sim10^{2.2-4.5}$ (in Habing units: $1.6\times10^{-3}{\rm erg~cm^{-2}~s^{-1}}$), consistent with both local ULIRGs and lower luminosity star-forming galaxies.

We use new Band-3 CO(1-0) observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA) to study the physical conditions in the interstellar gas of a sample of 27 dusty main-sequence star-forming galaxies at 0.03<$z$<0.2 present in the Valpara\'iso ALMA Line Emission Survey (VALES). The sample is drawn from far-IR bright galaxies over $\sim$160 deg$^{2}$ in the Herschel Astrophysical Terahertz Large Area Survey (HATLAS), which is covered by Herschel [CII] 158 $\mu$m spectroscopy and far-infrared (FIR) photometry. The [CII] and CO lines are both detected at >5$\sigma$ in 26 sources. We find an average [CII] to CO(1-0) luminosity ratio of 3500$\pm$1200 for our sample that is consistent with previous studies. Using the [CII], CO and FIR measurements as diagnostics of the physical conditions of the interstellar medium, we compare these observations to the predictions of a photodissociation region (PDR) model to determine the gas density, surface temperature, pressure, and the strength of the incident far-ultraviolet (FUV) radiation field, $G_{0}$, normalised to the Habing Field. The majority of our sample exhibit hydrogen densities of 4 < $\log n/\mathrm{cm}^{3}$ < 5.5 and experience an incident FUV radiation field with strengths of 2 < $\log G_0$ < 3 when adopting standard adjustments. A comparison to galaxy samples at different redshifts indicates that the average strength of the FUV radiation field appears constant up to redshift $z\sim$6.4, yet the neutral gas density increases with redshift by a factor of $\sim$100, that persists regardless of various adjustments to our observable quantities. This evolution could provide an explanation for the observed evolution of the star formation rate density with cosmic time, yet could arise from a combination of observational biases when using different suites of emission lines as diagnostic tracers of PDR gas.

We present a sample of 80 candidate strongly lensed galaxies with flux density above 100mJy at 500{\mu}m extracted from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), over an area of 600 square degrees. Available imaging and spectroscopic data allow us to confirm the strong lensing in 20 cases and to reject it in one case. For other 8 objects the lensing scenario is strongly supported by the presence of two sources along the same line of sight with distinct photometric redshifts. The remaining objects await more follow-up observations to confirm their nature. The lenses and the background sources have median redshifts z_L = 0.6 and z_S = 2.5, respectively, and are observed out to z_L = 1.2 and z_S = 4.2. We measure the number counts of candidate lensed galaxies at 500{\mu}m and compare them with theoretical predictions, finding a good agreement for a maximum magnification of the background sources in the range 10-20. These values are consistent with the magnification factors derived from the lens modelling of individual systems. The catalogue presented here provides sub- mm bright targets for follow-up observations aimed at exploiting gravitational lensing to study with un-precedented details the morphological and dynamical properties of dusty star forming regions in z >~ 1.5 galaxies.

The purported spiral host galaxy of GRB 020819B at z=0.41 has been seminal in establishing our view of the diversity of long-duration gamma-ray burst environments: optical spectroscopy of this host provided evidence that GRBs can form even at high metallicities, while millimetric observations suggested that GRBs may preferentially form in regions with minimal molecular gas. We report new observations from VLT (MUSE and X-shooter) which demonstrate that the purported host is an unrelated foreground galaxy. The probable radio afterglow is coincident with a compact, highly star-forming, dusty galaxy at z=1.9621. The revised redshift naturally explains the apparent nondetection of CO(3-2) line emission at the afterglow site from ALMA. There is no evidence that molecular gas properties in GRB host galaxies are unusual, and limited evidence that GRBs can form readily at super-Solar metallicity.

We present deep observations at 450 um and 850 um in the Extended Groth Strip field taken with the SCUBA-2 camera mounted on the James Clerk Maxwell Telescope as part of the deep SCUBA-2 Cosmology Legacy Survey (S2CLS), achieving a central instrumental depth of $\sigma_{450}=1.2$ mJy/beam and $\sigma_{850}=0.2$ mJy/beam. We detect 57 sources at 450 um and 90 at 850 um with S/N > 3.5 over ~70 sq. arcmin. From these detections we derive the number counts at flux densities $S_{450}>4.0$ mJy and $S_{850}>0.9$ mJy, which represent the deepest number counts at these wavelengths derived using directly extracted sources from only blank-field observations with a single-dish telescope. Our measurements smoothly connect the gap between previous shallower blank-field single-dish observations and deep interferometric ALMA results. We estimate the contribution of our SCUBA-2 detected galaxies to the cosmic infrared background (CIB), as well as the contribution of 24 um-selected galaxies through a stacking technique, which add a total of $0.26\pm0.03$ and $0.07\pm0.01$ MJy/sr, at 450 um and 850 um, respectively. These surface brightnesses correspond to $60\pm20$ and $50\pm20$ per cent of the total CIB measurements, where the errors are dominated by those of the total CIB. Using the photometric redshifts of the 24 um-selected sample and the redshift distributions of the submillimetre galaxies, we find that the redshift distribution of the recovered CIB is different at each wavelength, with a peak at $z\sim1$ for 450 um and at $z\sim2$ for 850um, consistent with previous observations and theoretical models.

Submillimetre galaxies (SMGs) are among the most luminous dusty galaxies in the Universe, but their true nature remains unclear; are SMGs the progenitors of the massive elliptical galaxies we see in the local Universe, or are they just a short-lived phase among more typical star-forming galaxies? To explore this problem further, we investigate the clustering of SMGs identified in the SCUBA-2 Cosmology Legacy Survey. We use a catalogue of submillimetre ($850\mu$m) source identifications derived using a combination of radio counterparts and colour/IR selection to analyse a sample of 914 SMGs in the UKIDSS Ultra Deep Survey (UDS), making this the largest high redshift sample of these galaxies to date. Using angular cross-correlation techniques, we estimate the halo masses for this large sample of SMGs and compare them with passive and star-forming galaxies selected in the same field. We find that SMGs, on average, occupy high-mass dark matter halos (M$_{\text{halo}} >10^{13}$M$_{\odot}$) at redshifts $z > 2.5$, consistent with being the progenitors of massive quiescent galaxies in present-day galaxy clusters. We also find evidence of downsizing, in which SMG activity shifts to lower mass halos at lower redshifts. In terms of their clustering and halo masses, SMGs appear to be consistent with other star-forming galaxies at a given redshift.

(abridged) We report rest-frame submillimeter H2O emission line observations of 11 HyLIRGs/ULIRGs at z~2-4 selected among the brightest lensed galaxies discovered in the Herschel-ATLAS. Using the IRAM NOEMA, we have detected 14 new H2O emission lines. The apparent luminosities of the H2O emission lines are $\mu L_{\rm{H_2O}} \sim 6-21 \times 10^8 L_\odot$, with velocity-integrated line fluxes ranging from 4-15 Jy km s$^{-1}$. We have also observed CO emission lines using EMIR on the IRAM 30m telescope in seven sources. The velocity widths for CO and H2O lines are found to be similar. With almost comparable integrated flux densities to those of the high-J CO line, H2O is found to be among the strongest molecular emitters in high-z Hy/ULIRGs. We also confirm our previously found correlation between luminosity of H2O ($L_{\rm{H_2O}}$) and infrared ($L_{\rm{IR}}$) that $L_{\rm{H_2O}} \sim L_{\rm{IR}}^{1.1-1.2}$, with our new detections. This correlation could be explained by a dominant role of far-infrared (FIR) pumping in the H2O excitation. Modelling reveals the FIR radiation fields have warm dust temperature $T_\rm{warm}$~45-75 K, H2O column density per unit velocity interval $N_{\rm{H_2O}}/\Delta V \gtrsim 0.3 \times 10^{15}$ cm$^{-2}$ km$^{-1}$ s and 100 $\mu$m continuum opacity $\tau_{100} > 1$ (optically thick), indicating that H2O is likely to trace highly obscured warm dense gas. However, further observations of $J\geq4$ H2O lines are needed to better constrain the continuum optical depth and other physical conditions of the molecular gas and dust. We have also detected H2O+ emission in three sources. A tight correlation between $L_{\rm{H_2O}}$ and $L_{\rm{H_2O^+}}$ has been found in galaxies from low to high redshift. The velocity-integrated flux density ratio between H2O+ and H2O suggests that cosmic rays generated by strong star formation are possibly driving the H2O+ formation.

We present a source-plane reconstruction of a ${\it Herschel}$ and ${\it Planck}$-detected gravitationally-lensed dusty star-forming galaxy (DSFG) at $z=1.68$ using {\it Hubble}, Sub-millimeter Array (SMA), and Keck observations. The background sub-millimeter galaxy (SMG) is strongly lensed by a foreground galaxy cluster at $z=0.997$ and appears as an arc of length $\sim 15^{\prime \prime}$ in the optical images. The continuum dust emission, as seen by SMA, is limited to a single knot within this arc. We present a lens model with source plane reconstructions at several wavelengths to show the difference in magnification between the stars and dust, and highlight the importance of a multi-wavelength lens models for studies involving lensed DSFGs. We estimate the physical properties of the galaxy by fitting the flux densities to model SEDs leading to a magnification-corrected star formation rate of $390 \pm 60$ M$_{\odot}$ yr$^{-1}$ and a stellar mass of $1.1 \pm 0.4\times 10^{11}$ M$_{\odot}$. These values are consistent with high-redshift massive galaxies that have formed most of their stars already. The estimated gas-to-baryon fraction, molecular gas surface density, and SFR surface density have values of $0.43 \pm 0.13$, $350 \pm 200$ M$_{\odot}$ pc$^{-2}$, and $\sim 12 \pm 7~$M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$, respectively. The ratio of star formation rate surface density to molecular gas surface density puts this among the most star-forming systems, similar to other measured SMGs and local ULIRGS.

We have recently suggested that gas accretion can be studied using host galaxies of gamma-ray bursts (GRBs). We obtained the first ever far-infrared (FIR) line observations of a GRB host, namely Herschel/PACS resolved [CII] 158 um and [OI] 63 um spectroscopy, as well as APEX CO(2-1) and ALMA CO(1-0) observations of the GRB 980425 host. It has elevated [CII]/FIR and [OI]/FIR ratios and higher values of star formation rate (SFR) derived from line ([CII], [OI], Ha) than from continuum (UV, IR, radio) indicators. [CII] emission exhibits a normal morphology, peaking at the galaxy center, whereas [OI] is concentrated close to the GRB position and the nearby Wolf-Rayet region. The high [OI] flux indicates high radiation field and gas density. The [CII]/CO luminosity ratio of the GRB 980425 host is close to the highest values found for local star-forming galaxies. Its CO-derived molecular gas mass is low given its SFR and metallicity, but the [CII]-derived molecular gas mass is close to the expected value. The [OI] and HI concentrations, and the high radiation field and density are consistent with the hypothesis of a very recent (at most a few tens of Myr ago) inflow of atomic gas triggering star formation. Dust has not had time to build up (explaining high line-to-continuum ratios). Such a recent enhancement of star-formation would indeed manifest itself in high SFR_line/SFR_continuum ratios, because the line indicators are sensitive only to recent (<10 Myr) activity, whereas the continuum indicators measure the SFR averaged over much longer periods (~100 Myr). Other GRB hosts exhibit a mean SFR_line/SFR_continuum of 1.74+-0.32. This is consistent with a very recent enhancement of star formation being common among GRB hosts, so galaxies which have recently experienced inflow of gas may preferentially host stars exploding as GRBs. Hence GRB hosts may be used to investigate recent gas accretion.

Extinction curves, especially those in the Milky Way (MW), the Large Magellanic Cloud (LMC), and the Small Magellanic Cloud (SMC), have provided us with a clue to the dust properties in the nearby Universe. We examine whether or not these extinction curves can be explained by well known dust evolution processes. We treat the dust production in stellar ejecta, destruction in supernova shocks, dust growth by accretion and coagulation, and dust disruption by shattering. To make a survey of the large parameter space possible, we simplify the treatment of the grain size distribution evolution by adopting the `two-size approximation', in which we divide the grain population into small ($\lesssim 0.03~\mu$m) and large ($\gtrsim 0.03~\mu$m) grains. It is confirmed that the MW extinction curve can be reproduced in reasonable ranges for the time-scale of the above processes with a silicate-graphite mixture. This indicates that the MW extinction curve is a natural consequence of the dust evolution through the above processes. We also find that the same models fail to reproduce the SMC/LMC extinction curves. Nevertheless, this failure can be remedied by giving higher supernova destruction rates for small carbonaceous dust and considering amorphous carbon for carbonaceous dust; these modification fall in fact in line with previous studies. Therefore, we conclude that the current dust evolution scenario composed of the aforementioned processes is successful in explaining the extinction curves. All the extinction curves favor efficient interstellar processing of dust, especially, strong grain growth by accretion and coagulation.

We present the first results from an on-going survey to characterize the circumgalactic medium (CGM) of the massive high-redshift galaxies detected as submillimeter galaxies (SMGs). We constructed a parent sample of 163 SMG-QSO pairs with separations less than $\sim$36 arcsec by cross-matching far-infrared-selected galaxies from Herschel with spectroscopically confirmed QSOs. The Herschel sources were selected to match the properties of SMGs. We determined the sub-arcsecond positions of six Herschel sources with the Very Large Array and obtained secure redshift identification for three of those with near-infrared spectroscopy. The QSO sightlines probe transverse proper distances of 112, 157, and 198 kpc at foreground redshifts of 2.043, 2.515, and 2.184, respectively, which are comparable to the virial radius of the $\sim10^{13}$ Msun halos expected to host SMGs. High-quality absorption-line spectroscopy of the QSOs reveals systematically strong HI Lyman-alpha absorption around all three SMGs, with rest-frame equivalent widths of $\sim2-3$ \AA. However, none of the three absorbers exhibits compelling evidence for optically thick HI gas or metal absorption, in contrast to the dominance of strong neutral absorbers in the CGM of luminous $z \sim 2$ QSOs. The low covering factor of optically thick HI gas around SMGs tentatively indicates that SMGs may not have as prominent cool gas reservoirs in their halos as the co-eval QSOs and that they may inhabit less massive halos than previously thought.

Improving the capabilities of detecting faint X-ray sources is fundamental to increase the statistics on faint high-z AGN and star-forming galaxies. We performed a simultaneous Maximum Likelihood PSF fit in the [0.5-2] keV and [2-7] keV energy bands of the 4 Ms{\em Chandra} Deep Field South (CDFS) data at the position of the 34930 CANDELS H-band selected galaxies. For each detected source we provide X-ray photometry and optical counterpart validation. We validated this technique by means of a raytracing simulation. We detected a total of 698 X-ray point-sources with a likelihood $\mathcal{L}$$>$4.98 (i.e. $>$2.7$\sigma$). We show that the prior knowledge of a deep sample of Optical-NIR galaxies leads to a significant increase of the detection of faint (i.e. $\sim$10$^{-17}$ cgs in the [0.5-2] keV band) sources with respect to "blind" X-ray detections. By including previous X-ray catalogs, this work increases the total number of X-ray sources detected in the 4 Ms CDFS, CANDELS area to 793, which represents the largest sample of extremely faint X-ray sources assembled to date. Our results suggest that a large fraction of the optical counterparts of our X-ray sources determined by likelihood ratio actually coincides with the priors used for the source detection. Most of the new detected sources are likely star-forming galaxies or faint absorbed AGN. We identified a few sources sources with putative photometric redshift z$>$4. Despite the low number statistics and the uncertainties on the photo-z, this sample significantly increases the number of X--ray selected candidate high-z AGN.

We exploit a sample of ultra-faint high-redshift galaxies (demagnified HST $H_{160}$ magnitude $>30$) in the Frontier Fields clusters A2744 and M0416 to constrain a theoretical model for the UV luminosity function (LF) in the presence of photoionization feedback. The objects have been selected on the basis of accurate photometric redshifts computed from multi-band photometry including 7 HST bands and deep $K_s$ and IRAC observations. Magnification is computed on an object-by-object basis from all available lensing models of the two clusters. We take into account source detection completeness as a function of luminosity and size, magnification effects and systematics in the lens modeling of the clusters under investigation. We find that our sample of high-$z$ galaxies constrain the cut-off halo circular velocity below which star-formation is suppressed by photo-ionization feedback to $v_c^{\rm cut} < 50$ km s$^{-1}$. This circular velocity corresponds to a halo mass of $\approx5.6\times10^9~M_\odot$ and $\approx2.3\times10^9~M_\odot$ at $z=5$ and 10 respectively: higher mass halos can thus sustain continuous star formation activity without being quenched by external ionizing flux. More stringent constraints are prevented by the uncertainty in the modeling of the cluster lens, as embodied by systematic differences among the lens models available.

We present multiwavelength photometric catalogues (HST, Spitzer and Hawk-I K band) for the first two of the Frontier Fields, Abell2744 and MACSJ0416 (plus their parallel fields). 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 as well as the ICL. We perform the detection on the HST H160 processed image to obtain a pure H-selected sample. We also add a sample of sources which are undetected in the H160 image but appear on a stacked infrared image. Photometry in the other HST bands is obtained using SExtractor, performed again on residual images after the Galfit procedure for foreground light removal. Photometry on the Hawk-I and IRAC bands has been 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. The procedure 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 a 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 {\mu}m from Spitzer. Full and detailed treatment of photometric errors is included. We perform basic sanity checks on the reliability of our results. The multiwavelength catalogues are publicly available and are ready to be used for scientific purposes. Our procedures allows for the detection of outshined 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 a companion Paper II. [abridged]

We present the first public release of photometric redshifts, galaxy rest-frame properties and associated magnification values in the cluster and parallel pointings of the first two Frontier Fields, Abell-2744 and MACS-J0416. We exploit a multi-wavelength catalogue ranging from HST to ground-based K and Spitzer IRAC which is specifically designed to enable detection and measurement of accurate fluxes in crowded cluster regions. The multi-band information is used to derive photometric redshifts and physical properties of sources detected either in the H-band image alone or from a stack of four WFC3 bands. To minimize systematics median photometric redshifts are assembled from six different approaches to photo-z estimates. Their reliability is assessed through a comparison with available spectroscopic samples. State of the art lensing models are used to derive magnification values on an object-by-object basis by taking into account sources positions and redshifts. We show that photometric redshifts reach a remarkable ~3-5% accuracy. After accounting for magnification the H band number counts are found in agreement at bright magnitudes with number counts from the CANDELS fields, while extending the presently available samples to galaxies intrinsically as faint as H160~32-33 thanks to strong gravitational lensing. The Frontier Fields allow to probe the galaxy stellar mass distribution at 0.5-1.5 dex lower masses, depending on magnification, with respect to extragalactic wide fields, including sources at Mstar~ 10^7-10^8 Msun at z>5. Similarly, they allow the detection of objects with intrinsic SFRs>1dex lower than in the CANDELS fields reaching 0.1-1 Msun/yr at z~6-10. [abridged]

Exploiting the sensitivity and spatial resolution of the Atacama Large Millimeter/submillimeter Array (ALMA), we have studied the morphology and the physical scale of the interstellar medium - both gas and dust - in SGP38326, an unlensed pair of interacting starbursts at $z= 4.425$. SGP38326 is the most luminous star bursting system known at $z > 4$ with an IR-derived ${\rm SFR \sim 4300 \,} M_\odot \, {\rm yr}^{-1}$. SGP38326 also contains a molecular gas reservoir among the most massive ever found in the early Universe, and it is the likely progenitor of a massive, red-and-dead elliptical galaxy at $z \sim 3$. Probing scales of $\sim 0.1"$ or $\sim 800 \, {\rm pc}$ we find that the smooth distribution of the continuum emission from cool dust grains contrasts with the more irregular morphology of the gas, as traced by the [CII] fine structure emission. The gas is also extended over larger physical scales than the dust. The velocity information provided by the resolved [CII] emission reveals that the dynamics of the two components of SGP38326 are compatible with disk-like, ordered rotation, but also reveals an ISM which is turbulent and unstable. Our observations support a scenario where at least a subset of the most distant extreme starbursts are highly dissipative mergers of gas-rich galaxies.

We introduce the Swift Gamma-Ray Burst Host Galaxy Legacy Survey ("SHOALS"), a multi-observatory high-redshift galaxy survey targeting the largest unbiased sample of long-duration gamma-ray burst hosts yet assembled (119 in total). We describe the motivations of the survey and the development of our selection criteria, including an assessment of the impact of various observability metrics on the success rate of afterglow-based redshift measurement. We briefly outline our host-galaxy observational program, consisting of deep Spitzer/IRAC imaging of every field supplemented by similarly-deep, multi-color optical/NIR photometry, plus spectroscopy of events without pre-existing redshifts. Our optimized selection cuts combined with host-galaxy follow-up have so far enabled redshift measurements for 110 targets (92%) and placed upper limits on all but one of the remainder. About 20% of GRBs in the sample are heavily dust-obscured, and at most 2% originate from z>5.5. Using this sample we estimate the redshift-dependent GRB rate density, showing it to peak at z~2.5 and fall by about an order of magnitude towards low (z=0) redshift, while declining more gradually towards high (z~7) redshift. This behavior is consistent with a progenitor whose formation efficiency varies modestly over cosmic history. Our survey will permit the most detailed examination to date of the connection between the GRB host population and general star-forming galaxies, directly measure evolution in the host population over cosmic time and discern its causes, and provide new constraints on the fraction of cosmic star-formation occurring in undetectable galaxies at all redshifts.

We present rest-frame NIR luminosities and stellar masses for a large and uniformly-selected population of GRB host galaxies using deep Spitzer Space Telescope imaging of 119 targets from the Swift GRB Host Galaxy Legacy Survey spanning 0.03 < z < 6.3, and determine the effects of galaxy evolution and chemical enrichment on the mass distribution of the GRB host population across cosmic history. We find strong evolution in the host luminosity distribution between z~0.5 (median absolute NIR AB magnitude ~ -18.5, corresponding to M* ~ 3x10^8 M_sun and z~1.5), but negligible variation between z~1.5 and z~5 (median magnitude ~ -21.2, corresponding to M* ~ 5x10^9 M_sun). Dust-obscured GRBs dominate the massive host population but are only rarely seen associated with low-mass hosts, indicating that massive star-forming galaxies are universally and (to some extent) homogeneously dusty at high-redshift while low-mass star-forming galaxies retain little dust in their ISM. Comparing our luminosity distributions to field surveys and measurements of the high-z mass-metallicity relation, our results have good consistency with a model in which the GRB rate per unit star-formation is constant in galaxies with gas-phase metallicity below approximately the Solar value but heavily suppressed in more metal-rich environments. This model also naturally explains the previously-reported "excess" in the GRB rate beyond z>2; metals stifle GRB production in most galaxies at z<1.5 but have only minor impact at higher redshifts. The metallicity threshold we infer is much higher than predicted by single-star models and favors a binary progenitor. Our observations also constrain the fraction of cosmic star-formation in low-mass galaxies undetectable to Spitzer to be a small minority at most redshifts (~10% at z~2, ~25% at z~3, and ~50% at z=3.5-6.0).

We present multiwavelength identifications for the counterparts of 1088 submillimeter sources detected at 850$\mu$m in the SCUBA-2 Cosmology Legacy Survey study of the UKIDSS-UDS field. By utilising an ALMA pilot study on a subset of our bright SCUBA-2 sample as a training set, along with the deep optical-near-infrared data available in this field, we develop a novel technique, Optical-IR Triple Color (OIRTC), using $z-K$, $K-[3.6]$, $[3.6]-[4.5]$ colors to select the candidate submillimeter galaxy (SMG) counterparts. By combining radio identification and the OIRTC technique, we find counterpart candidates for 80% of the Class = 1 $\geq4\,\sigma$ SCUBA-2 sample, defined as those that are covered by both radio and OIR imaging and the base sample for our scientific analyses. Based on the ALMA training set, we expect the accuracy of these identifications to be $82\pm20$%, with a completeness of $69\pm16$%, essentially as accurate as the traditional $p$-value technique but with higher completeness. We find that the fraction of SCUBA-2 sources having candidate counterparts is lower for fainter 850$\mu$m sources, and we argue that for follow-up observations sensitive to SMGs with $S_{850}\gtrsim 1$ mJy across the whole ALMA beam, the fraction with multiple counterparts is likely to be $>40$% for SCUBA-2 sources at $S_{850} \gtrsim 4$ mJy. We find that the photometric redshift distribution for the SMGs is well fit by a lognormal distribution, with a median redshift of $z=2.3\pm0.1$. After accounting for the sources without any radio and/or OIRTC counterpart, we estimate the median redshift to be $z=2.6\pm0.1$ for SMGs with $S_{850} >1$ mJy. We also use this new large sample to study the clustering of SMGs and the the far-infrared properties of the unidentified submillimeter sources by stacking their Herschel SPIRE far-infrared emission.

We introduce a new color-selection technique to identify high-redshift, massive galaxies that are systematically missed by Lyman-break selection. The new selection is based on the H_{160} and IRAC 4.5um bands, specifically H - [4.5] > 2.25 mag. These galaxies, dubbed "HIEROs", include two major populations that can be separated with an additional J - H color. The populations are massive and dusty star-forming galaxies at z > 3 (JH-blue) and extremely dusty galaxies at z < 3 (JH-red). The 350 arcmin^2 of the GOODS-N and GOODS-S fields with the deepest HST/WFC3 and IRAC data contain 285 HIEROs down to [4.5] < 24 mag. We focus here primarily on JH-blue (z > 3) HIEROs, which have a median photometric redshift z ~4.4 and stellar massM_{*}~10^{10.6} Msun, and are much fainter in the rest-frame UV than similarly massive Lyman-break galaxies (LBGs). Their star formation rates (SFRs) reaches ~240 Msun yr^{-1} leading to a specific SFR, sSFR ~4.2 Gyr^{-1}, suggesting that the sSFRs for massive galaxies continue to grow at z > 2 but at a lower growth rate than from z=0 to z=2. With a median half-light radius of 2 kpc, including ~20% as compact as quiescent galaxies at similar redshifts, JH-blue HIEROs represent perfect star-forming progenitors of the most massive (M_{*} > 10^{11.2} Msun) compact quiescent galaxies at z ~ 3 and have the right number density. HIEROs make up ~60% of all galaxies with M_{*} > 10^{10.5} Msun identified at z > 3 from their photometric redshifts. This is five times more than LBGs with nearly no overlap between the two populations. While HIEROs make up 15-25% of the total SFR density at z ~ 4-5, they completely dominate the SFR density taking place in M_{*} >10^{10.5} Msun galaxies, and are therefore crucial to understanding the very early phase of massive galaxy formation.

Cosmic dust is an important component of the Universe, and its origin, especially at high redshifts, is still unknown. I present a simple but powerful method of assessing whether dust observed in a given galaxy could in principle have been formed by asymptotic giant branch (AGB) stars or supernovae (SNe). Using this method I show that for most of the galaxies with detected dust emission between z=4 and z=7.5 (1.5-0.7 billion years after the Big Bang) AGB stars are not numerous and efficient enough to be responsible for the measured dust masses. Supernovae could account for most of the dust, but only if all of them had efficiencies close to the maximal theoretically allowed value. This suggests that a different mechanism is responsible for dust production at high redshifts, and the most likely possibility is the grain growth in the interstellar medium.

Two millimeter observations of the MACS J1149.6+2223 cluster have detected a source that was consistent with the location of the lensed MACS1149-JD galaxy at z=9.6. A positive identification would have rendered this galaxy as the youngest dust forming galaxy in the universe. Follow up observation with the AzTEC 1.1 mm camera and the IRAM NOrthern Extended Millimeter Array (NOEMA) at 1.3 mm have not confirmed this association. In this paper we show that the NOEMA observations associate the 2 mm source with [PCB2012] 2882 ([PCB2012] 2882 is the NED-searchable name for this source.), source number 2882 in the Hubble Space Telescope ( HST) Cluster Lensing and Supernova (CLASH) catalog of MACS J1149.6+2223. This source, hereafter referred to as CLASH 2882, is a gravitationally lensed spiral galaxy at z=0.99. We combine the GISMO 2 mm and NOEMA 1.3 mm fluxes with other (rest frame) UV to far-IR observations to construct the full spectral energy distribution (SED) of this galaxy, and derive its star formation history, and stellar and interstellar dust content. The current star formation rate of the galaxy is 54/mu Msun yr-1, and its dust mass is about 5 10^7/mu Msun, where mu is the lensing magnification factor for this source, which has a mean value of 2.7. The inferred dust mass is higher than the maximum dust mass that can be produced by core collapse supernovae (CCSN) and evolved AGB stars. As with many other star forming galaxies, most of the dust mass in CLASH 2882 must have been accreted in the dense phases of the ISM.

Recent observations with the GISMO (Goddard-IRAM Superconducting 2 Millimeter Observer) 2 mm camera revealed a detection 8 arcsec away from the lensed galaxy MACS1149-JD1 at z=9.6. Within the 17.5 arcsec FWHM GISMO beam, this detection is consistent with the position of the high-redshift galaxy and therefore, if confirmed, this object could be claimed to be the youngest galaxy producing significant quantities of dust. We present higher resolution (8.5 arcsec) observations of this system taken with the AzTEC 1.1 mm camera mounted on the Large Millimeter Telescope Alfonso Serrano. Dust continuum emission at the position of MACS1149-JD1 is not detected with an r.m.s. of 0.17 mJy/beam. However, we find a detection ~ 11 arcsec away from MACS1149-JD1, still within the GISMO beam which is consistent with an association to the GISMO source. Combining the AzTEC and GISMO photometry, together with Herschel ancillary data, we derive a z_phot= 0.7-1.6 for the dusty galaxy. We conclude therefore that the GISMO and AzTEC detections are not associated with MACS1149-JD1. From the non-detection of MACS1149-JD1 we derive the following (3 \sigma) upper limits corrected for gravitational lensing magnification and for cosmic microwave background effects: dust mass < 1.6 x 10^7 M_sun, IR luminosity < 8 x 10^10 L_sun, star formation rate < 14 M_sun/yr, and UV attenuation < 2.7 mag. These limits are comparable to those derived for other high-redshift galaxies from deep Atacama Large Millimeter/submillimeter Array (ALMA) observations.

We report the detection and mapping of atomic hydrogen in HI 21cm emission from ESO 184-G82, the host galaxy of the gamma ray burst 980425. This is the first instance where HI in emission has been detected from a galaxy hosting a gamma ray burst. ESO 184-G82 is an isolated galaxy and contains a Wolf-Rayet region close to the location of the gamma ray burst and the associated supernova, SN 1998bw. This is one of the most luminous HII regions identified in the local Universe, with a very high inferred density of star formation. The HI 21cm observations reveal a high HI mass for the galaxy, twice as large as the stellar mass. The spatial and velocity distribution of the HI 21cm emission reveals a disturbed rotating gas disk, which suggests that the galaxy has undergone a recent minor merger that disrupted its rotation. We find that the Wolf-Rayet region and the gamma ray burst are both located in the highest HI column density region of the galaxy. We speculate that the merger event has resulted in shock compression of the gas, triggering extreme star formation activity, and resulting in the formation of both the Wolf-Rayet region and the gamma ray burst. The high HI column density environment of the GRB is consistent with the high HI column densities seen in absorption in the host galaxies of high redshift gamma ray bursts.

Long gamma-ray bursts (GRBs), among the most energetic events in the Universe, are explosions of massive and short-lived stars, so they pinpoint locations of recent star formation. However, several GRB host galaxies have recently been found to be deficient in molecular gas (H2), believed to be the fuel of star formation. Moreover, optical spectroscopy of GRB afterglows implies that the molecular phase constitutes only a small fraction of the gas along the GRB line-of-sight. Here we report the first ever 21 cm line observations of GRB host galaxies, using the Australia Telescope Compact Array, implying high levels of atomic hydrogen (HI), which suggests that the connection between atomic gas and star formation is stronger than previously thought, with star formation being potentially directly fuelled by atomic gas (or with very efficient HI-to-H2 conversion and rapid exhaustion of molecular gas), as has been theoretically shown to be possible. This can happen in low metallicity gas near the onset of star formation, because cooling of gas (necessary for star formation) is faster than the HI-to-H2 conversion. Indeed, large atomic gas reservoirs, together with low molecular gas masses, stellar and dust masses are consistent with GRB hosts being preferentially galaxies which have very recently started a star formation episode after accreting metal-poor gas from the intergalactic medium. This provides a natural route for forming GRBs in low-metallicity environments. The gas inflow scenario is also consistent with the existence of the companion HI object with no optical counterpart ~19 kpc from the GRB 060505 host, and with the fact that the HI centroids of the GRB 980425 and 060505 hosts do not coincide with optical centres of these galaxies, but are located close to the GRB positions.

Interstellar dust in galaxies can be traced either through its extinction effects on the star light, or through its thermal emission at infrared wavelengths. Recent radiative transfer studies of several nearby edge-on galaxies have found an apparent inconsistency in the dust energy balance: the radiative transfer models that successfully explain the optical extinction underestimate the observed fluxes by an average factor of three. We investigate the dust energy balance for IC4225 and NGC5166, two edge-on spiral galaxies observed by the Herschel Space Observatory in the frame of the H-ATLAS survey. We start from models which were constrained from optical data and extend them to construct the entire spectral energy distribution of our galaxies. These predicted values are subsequently compared to the observed far-infrared fluxes. We find that including a young stellar population in the modelling is necessary as it plays a non-negligible part in the heating of the dust grains. While the modelling approach for both galaxies is nearly identical, we find two very different results. As is often seen in other edge-on spiral galaxies, the far-infrared emission of our radiative transfer model of IC4225 underestimates the observed fluxes by a factor of about three. For NGC5166 on the other hand, we find that both the predicted spectral energy distribution as well as the simulated images match the observations particularly well. We explore possible reasons for this difference and conclude that it is unlikely that one single mechanism is the cause of the dust energy balance problem in spiral galaxies. We discuss the different approaches that can be considered in order to get a conclusive answer on the origin this discrepancy.

Dust plays an important role in our understanding of the Universe, but it is not obvious yet how the dust in the distant universe was formed. I derived the dust yields per asymptotic giant branch (AGB) star and per supernova (SN) required to explain dust masses of galaxies at z=6.3-7.5 (680-850 million years after the Big Bang) for which dust emission has been detected (HFLS3 at z=6.34, ULAS J1120+0641 at z=7.085, and A1689-zD1 at z=7.5), or unsuccessfully searched for. I found very high required yields, implying that AGB stars could not contribute substantially to dust production at these redshifts, and that SNe could explain these dust masses, but only if they do not destroy most of the dust they form (which is unlikely given the upper limits on the SN dust yields derived for galaxies where dust is not detected). This suggests that the grain growth in the interstellar medium is likely required at these early epochs.

We present Hubble Space Telescope (HST) WFC3 imaging and grism spectroscopy observations of the {\it Herschel}-selected gravitationally-lensed starburst galaxy HATLASJ1429-0028. The lensing system consists of an edge-on foreground disk galaxy at $z=0.218$ with a nearly complete Einstein ring of the infrared luminous galaxy at $z=1.027$. The WFC3 spectroscopy with G102 and G141 grisms, covering the wavelength range of 0.8 to 1.7 $\mu$m, resulted in detections of H$\alpha$+[NII], H$\beta$, [SII], and [OIII] for the background galaxy from which we measure line fluxes and ratios. The Balmer line ratio H$\alpha$/H$\beta$ of $7.5 \pm 4.4$, when corrected for [NII], results in an extinction for the starburst galaxy of $E(B-V)=0.8 \pm 0.5$. The $H\alpha$ based star-formation rate, when corrected for extinction, is $60 \pm 50$ M$_{\odot}$ yr$^{-1}$, lower than the instantaneous star-formation rate of 390 $\pm$ 90 M$_{\odot}$ yr$^{-1}$ from the total IR luminosity. We also compare the nebular line ratios of HATLASJ1429-0028 with other star-forming and sub-mm bright galaxies. The nebular line ratios are consistent with an intrinsic ultra-luminous infrared galaxy with no evidence for excitation by an active galactic nucleus (AGN). We estimate the metallicity, $12 + log(O/H)$, of HATLASJ1429-0028 to be 8.49 $\pm$ 0.16. Such a low value is below the average relations for stellar mass vs. metallicity of galaxies at $z \sim 1$ for a galaxy with stellar mass of $\sim 2 \times 10^{11}$ M$_{\odot}$. The combination of high stellar mass, lack of AGN indicators, low metallicity, and the high star-formation rate of HATLASJ1429-0028 suggest that this galaxy is currently undergoing a rapid formation.

We explore the behaviour of [CII]-157.74um forbidden fine-structure line observed in a sample of 28 galaxies selected from ~50deg^2 of the H-ATLAS survey. The sample is restricted to galaxies with flux densities higher than S_160um>150mJy and optical spectra from the GAMA survey at 0.02<z<0.2. Far-IR spectra centred on this redshifted line were taken with the PACS instrument on-board the Herschel Space Observatory. The galaxies span 10<log(L_IR/Lo)<12 (where L_IR=L_IR[8-1000um]) and 7.3<log(L_[CII]/Lo)<9.3, covering a variety of optical galaxy morphologies. The sample exhibits the so-called [CII] deficit at high IR luminosities, i.e. L_[CII]/L_IR (hereafter [CII]/IR) decreases at high L_IR. We find significant differences between those galaxies presenting [CII]/IR>2.5x10^-3 with respect to those showing lower ratios. In particular, those with high ratios tend to have: (1) L_IR<10^11Lo; (2) cold dust temperatures, T_d<30K; (3) disk-like morphologies in r-band images; (4) a WISE colour 0.5<S_12um/S_22um<1.0; (5) low surface brightness Sigma_IR~10^8-9 Lo kpc^-2, (6) and specific star-formation rates of sSFR~0.05-3 Gyr^-1. We suggest that the strength of the far-UV radiation fields (<G_O>) is main parameter responsible for controlling the [CII]/IR ratio. It is possible that relatively high <G_O> creates a positively charged dust grain distribution, impeding an efficient photo-electric extraction of electrons from these grains to then collisionally excite carbon atoms. Within the brighter IR population, 11<log(L_IR/Lo)<12, the low [CII]/IR ratio is unlikely to be modified by [CII] self absorption or controlled by the presence of a moderately luminous AGN (identified via the BPT diagram).

Candidates for the modest galaxies that formed most of the stars in the early universe, at redshifts $z > 7$, have been found in large numbers with extremely deep restframe-UV imaging. But it has proved difficult for existing spectrographs to characterise them in the UV. The detailed properties of these galaxies could be measured from dust and cool gas emission at far-infrared wavelengths if the galaxies have become sufficiently enriched in dust and metals. So far, however, the most distant UV-selected galaxy detected in dust emission is only at $z = 3.25$, and recent results have cast doubt on whether dust and molecules can be found in typical galaxies at this early epoch. Here we report thermal dust emission from an archetypal early universe star-forming galaxy, A1689-zD1. We detect its stellar continuum in spectroscopy and determine its redshift to be $z = 7.5\pm0.2$ from a spectroscopic detection of the Ly{\alpha} break. A1689-zD1 is representative of the star-forming population during reionisation, with a total star-formation rate of about 12M$_\odot$ yr$^{-1}$. The galaxy is highly evolved: it has a large stellar mass, and is heavily enriched in dust, with a dust-to-gas ratio close to that of the Milky Way. Dusty, evolved galaxies are thus present among the fainter star-forming population at $z > 7$, in spite of the very short time since they first appeared.

Luminous infrared galaxies and submillimeter galaxies contribute significantly to stellar mass assembly and the frequency of GRBs in these systems provides an important test of the connection between the gamma-ray burst rate and that of overall cosmic star-formation. We present sensitive 3 GHz radio observations using the Karl G. Jansky Very Large Array of 32 uniformly-selected GRB host galaxies spanning a redshift range from 0 < z < 2.5, providing the first fully dust- and sample-unbiased measurement of the fraction of GRBs originating from the Universe's most bolometrically luminous galaxies. Four galaxies are detected, with inferred radio star-formation rates ranging between 50-300 Msun/yr. Three of the four detections correspond to events consistent with being optically-obscured "dark" bursts. Our overall detection fraction implies that between 9-23% of GRBs between 0.5 < z < 2.5 occur in galaxies with S_3GHz > 10 uJy, corresponding to SFR >50 Msun/yr at z~1 or >250 Msun/yr at z~2. Similar galaxies contribute approximately 10-30% of all cosmic star-formation, so our results are consistent with a GRB rate which is not strongly biased with respect to the total star-formation rate of a galaxy. However, all four radio-detected hosts have modest stellar masses (~few x 10^10 Msun), significantly lower than IR/submillimeter-selected field galaxies of similar luminosities. We suggest that the GRB rate may be suppressed in metal-rich environments but independently enhanced in intense starbursts, producing a strong efficiency dependence on mass but little net dependence on bulk galaxy star-formation rate.

Gamma-ray bursts (GRBs) are the most energetic phenomena in the Universe; believed to result from the collapse and subsequent explosion of massive stars. Even though it has profound consequences for our understanding of their nature and selection biases, little is known about the dust properties of the galaxies hosting GRBs. We present analysis of the far-infrared properties of an unbiased sample of 20 \textit{BeppoSAX} and \textit{Swift} GRB host galaxies (at an average redshift of $z\,=\,3.1$) located in the {\it Herschel} Astrophysical Terahertz Large Area Survey, the {\it Herschel} Virgo Cluster Survey, the {\it Herschel} Fornax Cluster Survey, the {\it Herschel} Stripe 82 Survey and the {\it Herschel} Multi-tiered Extragalactic Survey, totalling $880$ deg$^2$, or $\sim 3$\% of the sky in total. Our sample selection is serendipitous, based only on whether the X-ray position of a GRB lies within a large-scale {\it Herschel} survey -- therefore our sample can be considered completely unbiased. Using deep data at wavelengths of 100\,--\,500$\,\mu$m, we tentatively detected 1 out of 20 GRB hosts located in these fields. We constrain their dust masses and star formation rates (SFRs), and discuss these in the context of recent measurements of submillimetre galaxies and ultraluminous infrared galaxies. The average far-infrared flux of our sample gives an upper limit on SFR of $<114\,{\rm M}\odot\,\mbox{yr}^{-1}$. The detection rate of GRB hosts is consistent with that predicted assuming that GRBs trace the cosmic SFR density in an unbiased way, i.e. that the fraction of GRB hosts with $\mbox{SFR}>500\,{\rm M}\odot\,\mbox{yr}^{-1}$ is consistent with the contribution of such luminous galaxies to the cosmic star formation density.

Long-duration gamma-ray bursts (LGRBs) are produced by the collapse of very massive stars. Due to the short lifetime of their progenitors, LGRBs pinpoint star-forming galaxies. We present here a multi-band search for the host galaxy of the long dark GRB 050219A within the enhanced Swift/XRT error circle. We used spectroscopic observations acquired with VLT/X-shooter to determine the redshift and star-formation rate of the putative host galaxy. We compared the results with the optical/IR spectral energy distribution obtained with different facilities. Surprisingly, the host galaxy is a old and quiescent early-type galaxy at z = 0.211 characterised by an unprecedentedly low specific star-formation rate. It is the first LGRB host to be also an early-type post-starburst galaxy. This is further evidence that GRBs can explode in all kind of galaxies, with the only requirement being an episode of star-formation.

Submillimetre (submm) galaxies are among the most rapidly star-forming and most massive high-redshift galaxies; thus, their properties provide important constraints on galaxy evolution models. However, there is still a debate about their stellar masses and their nature in the context of the general galaxy population. To test the reliability of their stellar mass determinations, we used a sample of simulated submm galaxies for which we derived stellar masses via spectral energy distribution (SED) modelling (with Grasil, Magphys, Hyperz and LePhare) adopting various star formation histories (SFHs). We found that the assumption of SFHs with two independent components leads to the most accurate stellar masses. Exponentially declining SFHs (tau) lead to lower masses (albeit still consistent with the true values), while the assumption of single-burst SFHs results in a significant mass underestimation. Thus, we conclude that studies based on the higher masses inferred from fitting the SEDs of real submm galaxies with double SFHs are most likely to be correct, implying that submm galaxies lie on the high-mass end of the main sequence of star-forming galaxies. This conclusion appears robust to assumptions of whether or not submm galaxies are driven by major mergers, since the suite of simulated galaxies modelled here contains examples of both merging and isolated galaxies. We identified discrepancies between the true and inferred stellar ages (rather than the dust attenuation) as the primary determinant of the success/failure of the mass recovery. Regardless of the choice of SFH, the SED-derived stellar masses exhibit a factor of ~2 scatter around the true value; this scatter is an inherent limitation of the SED modelling due to simplified assumptions. Finally, we found that the contribution of active galactic nuclei does not have any significant impact on the derived stellar masses.

We present statistically significant detections at 850um of the Lyman Break Galaxy (LBG) population at z=3, 4, and 5 using data from the Submillimetre Common User Bolometer Array 2 (SCUBA-2) Cosmology Legacy Survey (S2CLS) in the United Kingdom Infrared Deep Sky Survey Ultra Deep Survey (UKIDSS-UDS) field. We employ a stacking technique to probe beneath the survey limit to measure the average 850um flux density of LBGs at z=3, 4, and 5 with typical ultraviolet luminosities of L(1700A)~10^29 erg/s/Hz. We measure 850um flux densities of (0.25 +/- 0.03, (0.41 +/- 0.06), and (0.88 +/- 0.23) mJy respectively, and find that they contribute at most 20 per cent to the cosmic far-infrared background at 850um. Fitting an appropriate range of spectral energy distributions to the z=3, 4, and 5 LBG stacked 24-850um fluxes, we derive infrared (IR) luminosities of L(8-1000um)~3.2, 5.5, and 11.0x10^11 Lsun (corresponding to star formation rates of ~50-200 Msun/yr) respectively. We find that the evolution in the IR luminosity density of LBGs is broadly consistent with model predictions for the expected contribution of luminous IR galaxy (LIRG) to ultraluminous IR galaxy (ULIRG) type systems at these epochs. We also see a strong positive correlation between stellar mass and IR luminosity. Our data are consistent with the main sequence of star formation showing little or no evolution from z=3 to 5. We have also confirmed that, for a fixed mass, the reddest LBGs (UV slope Beta -> 0) are indeed redder due to dust extinction, with SFR(IR)/SFR(UV) increasing by approximately an order of magnitude over -2<Beta<0 such that SFR(IR)/SFR(UV)~20 for the reddest LBGs. Furthermore, the most massive LBGs also tend to have higher obscured-to-unobscured ratio, hinting at a variation in the obscuration properties across the mass range.

We report a highly significant ($>10\sigma$) spatial correlation between galaxies with $S_{350\mu\rm m}\ge 30\,$mJy detected in the equatorial fields of the \textsl{Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS) with estimated redshifts $\gtrsim 1.5$, and SDSS or GAMA galaxies at $0.2\le z\le 0.6$. The significance of the cross-correlation is much higher than those reported so far for samples with non-overlapping redshift distributions selected in other wavebands. Extensive, realistic simulations of clustered sub-mm galaxies amplified by foreground structures confirm that the cross-correlation is explained by weak gravitational lensing ($\mu<2$). The simulations also show that the measured amplitude and range of angular scales of the signal are larger than can be accounted for by galaxy-galaxy weak lensing. However, for scales $\lesssim 2\,$arcmin, the signal can be reproduced if SDSS/GAMA galaxies act as signposts of galaxy groups/clusters with halo masses in the range $10^{13.2}$--$10^{14.5} M_{\odot}$. The signal detected on larger scales appears to reflect the clustering of such halos.

Gamma-ray bursts (GRBs) have been proposed as a tool for studying star formation in the Universe, so it is crucial to investigate whether their host galaxies and immediate environments are in any way special compared with other star-forming galaxies. Here we present spatially resolved maps of dust emission of the host galaxy of the closest known GRB 980425 at z=0.0085 using our new high-resolution observations from Herschel, Atacama Pathfinder Experiment (APEX), Atacama Large Millimeter Array (ALMA) and Australia Telescope Compact Array (ATCA). We modelled the spectral energy distributions of the host and of the star-forming region displaying the Wolf-Rayet signatures in the spectrum (WR region), located 800 pc from the GRB position. The host is characterised by low dust content and a high fraction of UV-visible star formation, similar to other dwarf galaxies. These galaxies are abundant in the local universe, so it is not surprising to find a GRB in one of them, assuming the correspondence between the GRB rate and star formation. The WR region contributes substantially to the host emission at the far-infrared, millimetre, and radio wavelengths and we propose that this is a consequence of its high gas density. If dense environments are also found close to the positions of other GRBs, then the ISM density should also be considered, along with metallicity, an important factor influencing whether a given stellar population can produce a GRB.

Until recently, dust emission has been detected in very few host galaxies of gamma-ray bursts (GRBHs). With Herschel, we have now observed 17 GRBHs up to redshift z~3 and detected seven of them at infrared (IR) wavelengths. This relatively high detection rate (41%) may be due to the composition of our sample which at a median redshift of 1.1 is dominated by the hosts of dark GRBs. Although the numbers are small, statistics suggest that dark GRBs are more likely to be detected in the IR than their optically-bright counterparts. Combining our IR data with optical, near-infrared, and radio data from our own datasets and from the literature, we have constructed spectral energy distributions (SEDs) which span up to 6 orders of magnitude in wavelength. By fitting the SEDs, we have obtained stellar masses, dust masses, star-formation rate (SFR), and extinctions for our sample galaxies. We find that GRBHs are galaxies that tend to have a high specfic SFR (sSFR), and like other star-forming galaxies, their ratios of dust-to-stellar mass are well correlated with sSFR. We incorporate our Herschel sample into a larger compilation of GRBHs, and compare this combined sample to SFR-weighted median stellar masses of the widest, deepest galaxy survey to date. This is done in order to establish whether or not GRBs can be used as an unbiased tracer of cosmic comoving SFR density (SFRD) in the universe. In contrast with previous results, this comparison shows that GRBHs are medium-sized galaxies with relatively high sSFRs; stellar masses and sSFRs of GRBHs as a function of redshift are similar to what is expected for star-forming galaxy populations at similar redshifts. We conclude that there is no strong evidence that GRBs are biased tracers of SFRD; thus they should be able to reliably probe the SFRD to early epochs.

There is a remarkably tight relation between the observationally inferred dust masses and star-formation rates (SFRs) of SDSS galaxies, Mdust $\propto$ SFR$^{1.11}$ (Da Cunha et al. 2010). Here we extend the Mdust-SFR relation to the high end and show that it bends over at very large SFRs (i.e., dust masses are lower than predicted for a given SFR). We identify several distinct evolutionary processes in the diagram: (1) A star-bursting phase in which dust builds up rapidly at early times. The maximum attainable dust mass in this process is the cause of the bend-over of the relation. A high dust-formation efficiency, a bottom-light initial mass function, and negligible supernova shock dust destruction are required to produce sufficiently high dust masses. (2) A quiescent star-forming phase in which the subsequent parallel decline in dust mass and SFR gives rise to the Mdust-SFR relation, through astration and dust destruction. The dust-to-gas ratio is approximately constant along the relation. We show that the power-law slope of the Mdust-SFR relation is inversely proportional to the global Schmidt-Kennicutt law exponent (i.e., $\sim 0.9$) in simple chemical evolution models. (3) A quenching phase which causes star formation to drop while the dust mass stays roughly constant or drops proportionally. Combined with merging, these processes, as well as the range in total baryonic mass, give rise to a complex population of the diagram which adds significant scatter to the original Mdust-SFR relation. (4) At very high redshifts, a population of galaxies located significantly below the local relation is predicted.

Upon its completion the Herschel ATLAS (H-ATLAS) will be the largest submillimetre survey to date, detecting close to half-a-million sources. It will only be possible to measure spectroscopic redshifts for a small fraction of these sources. However, if the rest-frame spectral energy distribution (SED) of a typical H-ATLAS source is known, this SED and the observed Herschel fluxes can be used to estimate the redshifts of the H-ATLAS sources without spectroscopic redshifts. In this paper, we use a subset of 40 H-ATLAS sources with previously measured redshifts in the range 0.5<z<4.2 to derive a suitable average template for high redshift H-ATLAS sources. We find that a template with two dust components T_c = 23.9 K, T_h = 46.9 K and ratio of mass of cold dust to mass of warm dust of 30.1) provides a good fit to the rest-frame fluxes of the sources in our calibration sample. We use a jackknife technique to estimate the accuracy of the redshifts estimated with this template, finding a root mean square of Delta z/(1+z) = 0.26. For sources for which there is prior information that they lie at z > 1 we estimate that the rms of Delta z/(1+z) = 0.12. We have used this template to estimate the redshift distribution for the sources detected in the H-ATLAS equatorial fields, finding a bimodal distribution with a mean redshift of 1.2, 1.9 and 2.5 for 250, 350 and 500 um selected sources respectively. \end{abstract}

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 investigate the multi-wavelength properties of a sample of 450-\mu m selected sources from the SCUBA-2 Cosmology Legacy Survey (S2CLS). A total of 69 sources were identified above 4\sigma\ in deep SCUBA-2 450-\mu m observations overlapping the UDS and COSMOS fields and covering 210 sq. arcmin to a typical depth of \sigma 450=1.5 mJy. Reliable cross identification are found for 58 sources (84 per cent) in Spitzer and Hubble Space Telescope WFC3/IR data. The photometric redshift distribution (dN/dz) of 450\mu m-selected sources is presented, showing a broad peak in the redshift range 1<z<3, and a median of z=1.4. Combining the SCUBA-2 photometry with Herschel SPIRE data from HerMES, the submm spectral energy distribution (SED) is examined via the use of modified blackbody fits, yielding aggregate values for the IR luminosity, dust temperature and emissivity of <LIR>=10^12 +/- 0.8 L_sol, <T_D>=42 +/- 11 K and <\beta_D>=1.6 +/- 0.5, respectively. The relationship between these SED parameters and the physical properties of galaxies is investigated, revealing correlations between T_D and LIR and between \beta_D and both stellar mass and effective radius. The connection between star formation rate and stellar mass is explored, with 24 per cent of 450 \mu m sources found to be ``star-bursts'', i.e. displaying anomalously high specific SFRs. However, both the number density and observed properties of these ``star-burst'' galaxies are found consistent with the population of normal star-forming galaxies.

We report the discovery of a well-defined correlation between B-band face-on central optical depth due to dust, \tau^f_B, and the stellar mass surface density, \mu_{*}, of nearby (z < 0.13) spiral galaxies: log(\tau^f_B) = 1.12(+-0.11)log(\mu_{*}/M_sol kpc^2)-8.6(+-0.8). This relation was derived from a sample of spiral galaxies taken from the Galaxy and Mass Assembly (GAMA) survey and detected in the FIR/submm in the Herschel-ATLAS survey. Using a quantitative analysis of the NUV attenuation-inclination relation for complete samples of GAMA spirals categorized according to \mu_{*} we demonstrate that this correlation can be used to statistically correct for dust attenuation purely on the basis of optical photometry and S'ersic-profile morphological fits. Considered together with previously established empirical relationships between stellar mass, metallicity and gas mass, the near linearity and high constant of proportionality of the \tau^f_B-\mu_{*} relation disfavors a stellar origin for the bulk of refractory grains in spiral galaxies, instead being consistent with the existence of a ubiquitous and very rapid mechanism for the growth of dust in the ISM. We use the \tau^f_B-\mu_{*} relation in conjunction with the radiation transfer model for spiral galaxies of Popescu & Tuffs (2011) to derive intrinsic scaling relations between specific star formation rate (sSFR), stellar mass, and \mu_{*}, in which the attenuation of the UV light used to measure the SFR is corrected on an object-to-object basis. A marked reduction in scatter in these relations is achieved which is demonstrably due to correction of both the inclination-dependent and face-on components of attenuation. Our results are consistent with a picture of spiral galaxies in which most of the submm emission originates from grains residing in translucent structures, exposed to UV in the diffuse interstellar radiation field.

The larger number of models of asteroid shapes and their rotational states derived by the lightcurve inversion give us better insight into both the nature of individual objects and the whole asteroid population. With a larger statistical sample we can study the physical properties of asteroid populations, such as main-belt asteroids or individual asteroid families, in more detail. Shape models can also be used in combination with other types of observational data (IR, adaptive optics images, stellar occultations), e.g., to determine sizes and thermal properties. We use all available photometric data of asteroids to derive their physical models by the lightcurve inversion method and compare the observed pole latitude distributions of all asteroids with known convex shape models with the simulated pole latitude distributions. We used classical dense photometric lightcurves from several sources and sparse-in-time photometry from the U.S. Naval Observatory in Flagstaff, Catalina Sky Survey, and La Palma surveys (IAU codes 689, 703, 950) in the lightcurve inversion method to determine asteroid convex models and their rotational states. We also extended a simple dynamical model for the spin evolution of asteroids used in our previous paper. We present 119 new asteroid models derived from combined dense and sparse-in-time photometry. We discuss the reliability of asteroid shape models derived only from Catalina Sky Survey data (IAU code 703) and present 20 such models. By using different values for a scaling parameter cYORP (corresponds to the magnitude of the YORP momentum) in the dynamical model for the spin evolution and by comparing synthetics and observed pole-latitude distributions, we were able to constrain the typical values of the cYORP parameter as between 0.05 and 0.6.

GRB-selected galaxies are broadly known to be faint, blue, young, star-forming dwarf galaxies. This insight, however, is based in part on heterogeneous samples of optically selected, lower-redshift galaxies. To study the statistical properties of GRB-selected galaxies we here introduce The Optically Unbiased GRB Host (TOUGH) complete sample of 69 X-ray selected Swift GRB host galaxies spanning the redshift range 0.03-6.30 and summarise the first results of a large observational survey of these galaxies.

Long-duration gamma-ray bursts (GRBs) are powerful tracers of star-forming galaxies. We have defined a homogeneous subsample of 69 Swift GRB-selected galaxies spanning a very wide redshift range. Special attention has been devoted to making the sample optically unbiased through simple and well-defined selection criteria based on the high-energy properties of the bursts and their positions on the sky. Thanks to our extensive follow-up observations, this sample has now achieved a comparatively high degree of redshift completeness, and thus provides a legacy sample, useful for statistical studies of GRBs and their host galaxies. In this paper we present the survey design and summarize the results of our observing program conducted at the ESO Very Large Telescope (VLT) aimed at obtaining the most basic properties of galaxies in this sample, including a catalog of R and Ks magnitudes and redshifts. We detect the host galaxies for 80 % of the GRBs in the sample, although only 42 % have Ks-band detections, which confirms that GRB-selected host galaxies are generally blue. The sample is not uniformly blue, however, with two extremely red objects detected. Moreover, galaxies hosting GRBs with no optical/NIR afterglows, whose identification therefore relies on X-ray localizations, are significantly brighter and redder than those with an optical/NIR afterglow. Our spectroscopic campaign has resulted in 77 % now having redshift measurements, with a median redshift of 2.14 +- 0.18. TOUGH alone includes 17 detected z > 2 Swift GRB host galaxies suitable for individual and statistical studies. Seven hosts have detections of the Ly-alpha emission line and we can exclude an early indication that Ly-alpha emission is ubiquitous among GRB hosts, but confirm that Ly-alpha is stronger in GRB-selected galaxies than in flux-limited samples of Lyman break galaxies.

The AzTEC 1.1 mm survey of the SCUBA HAlf Degree Extragalactic Survey (SHADES) fields is the largest (0.7 deg2) blank-field millimetre-wavelength survey undertaken to date at a resolution of ~18" and a depth of ~1 mJy. We have used the deep optical-to-radio multi-wavelength data in the SHADES Lockman Hole East and SXDF/UDS fields to obtain galaxy identifications for ~64% (~80% with tentative identifications) of the 148 AzTEC-SHADES 1.1 mm sources, exploiting deep radio and 24 um data complemented by methods based on 8 um flux-density and red optical-infrared (i-K) colour. This unusually high identification rate can be attributed to the relatively bright millimetre-wavelength flux-density threshold, combined with the relatively deep supporting multi-frequency data. We have further exploited the optical-mid-infrared-radio data to derive a ~60% (~75% with tentative identifications) complete redshift distribution for the AzTEC-SHADES sources, yielding a median redshift of z~2.2, with a high-redshift tail extending to at least z~4. Despite the larger area probed by the AzTEC survey relative to the original SCUBA SHADES imaging, the redshift distribution of the AzTEC sources is consistent with that displayed by the SCUBA sources, and reinforces tentative evidence that the redshift distribution of mm/sub-mm sources in the Lockman Hole field is significantly different from that found in the SXDF/UDS field. Comparison with simulated surveys of similar scale extracted from semi-analytic models based on the Millennium simulation indicates that this is as expected if the mm/sub-mm sources are massive (M>10^11 Mo) star-forming galaxies tracing large-scale structures over scales of 10-20 Mpc. This confirms the importance of surveys covering several square degrees (as now underway with SCUBA2) to obtain representative samples of bright (sub)mm-selected galaxies.

We use spitzer-IRAC data to identify near-infrared counterparts to submillimeter galaxies detected with Herschel-SPIRE at 250um in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). Using a likelihood ratio analysis we identify 146 reliable IRAC counterparts to 123 SPIRE sources out of the 159. We find that, compared to the field population, the SPIRE counterparts occupy a distinct region of 3.6 and 4.5um color-magnitude space, and we use this property to identify a further 23 counterparts to 13 SPIRE sources. The IRAC identification rate of 86% is significantly higher than those that have been demonstrated with wide-field ground-based optical and near-IR imaging of Herschel fields. We estimate a false identification rate of 3.6%, corresponding to 4 to 5 sources. Among the 73 counterparts that are undetected in SDSS, 57 have both 3.6 and 4.5um coverage. Of these 43 have [3.6] - [4.5]> 0 indicating that they are likely to be at z > 1.4. Thus, ~ 40% of identified SPIRE galaxies are likely to be high redshift (z > 1.4) sources. We discuss the statistical properties of the IRAC-identified SPIRE galaxy sample including far-IR luminosities, dust temperatures, star-formation rates, and stellar masses. The majority of our detected galaxies have 10^10 to 10^11 L_sun total IR luminosities and are not intense starbursting galaxies as those found at z ~ 2, but they have a factor of 2 to 3 above average specific star-formation rates compared to near-IR selected galaxy samples.

The objective of this paper is to determine the level of obscured star formation activity and dust attenuation in a sample of gamma-ray burst (GRB) hosts; and to test the hypothesis that GRB hosts have properties consistent with those of the general star-forming galaxy populations. We present a radio continuum survey of all z<1 GRB hosts in The Optically Unbiased GRB Host (TOUGH) sample supplemented with radio data for all (mostly pre-Swift) GRB-SN hosts discovered before October 2006. We present new radio data for 22 objects and have obtained a detection for three of them (GRB 980425, 021211, 031203; none in the TOUGH sample), increasing the number of radio-detected GRB hosts from two to five. The star formation rate (SFR) for the GRB 021211 host of ~825 Mo yr^-1, the highest ever reported for a GRB host, places it in the category of ultraluminous infrared galaxies. We found that at least 63% of GRB hosts have SFR < 100 Mo yr^-1 and at most 8% can have SFR > 500 Mo yr^-1. For the undetected hosts the mean radio flux (<35 uJy 3sigma) corresponds to an average SFR < 15 Mo yr^-1. Moreover, ~88% of the z<1 GRB hosts have ultraviolet dust attenuation A_UV < 6.7 mag (visual attenuation A_V < 3 mag). Hence we did not find evidence for large dust obscuration in a majority of GRB hosts. Finally, we found that the distributions of SFRs and A_UV of GRB hosts are consistent with those of Lyman break galaxies, Halpha emitters at similar redshifts and of galaxies from cosmological simulations. The similarity of the GRB population with other star-forming galaxies is consistent with the hypothesis that GRBs, a least at z<1, trace a large fraction of all star formation, and are therefore less biased indicators than once thought.

Establishing the stellar masses (M*), and hence specific star-formation rates (sSFRs) of submillimetre galaxies (SMGs) is crucial for determining their role in the cosmic galaxy/star formation. However, there is as yet no consensus over the typical M* of SMGs. Specifically, even for the same set of SMGs, the reported average M* have ranged over an order of magnitude, from ~5x10^10 Mo to ~5x10^11 Mo. Here we study how different methods of analysis can lead to such widely varying results. We find that, contrary to recent claims in the literature, potential contamination of IRAC 3-8 um photometry from hot dust associated with an active nucleus is not the origin of the published discrepancies in derived M*. Instead, we expose in detail how inferred M* depends on assumptions made in the photometric fitting, and quantify the individual and cumulative effects of different choices of initial mass function, different brands of evolutionary synthesis models, and different forms of assumed star-formation history. We review current observational evidence for and against these alternatives as well as clues from the hydrodynamical simulations, and conclude that, for the most justifiable choices of these model inputs, the average M* of SMGs is ~2x10^11 Mo. We also confirm that this number is perfectly reasonable in the light of the latest measurements of their dynamical masses, and the evolving M* function of the overall galaxy population. M* of this order imply that the average sSFR of SMGs is comparable to that of other star-forming galaxies at z>2, at 2-3 Gyr^-1. This supports the view that, while rare outliers may be found at any M*, most SMGs simply form the top end of the main-sequence of star-forming galaxies at these redshifts. Conversely, this argues strongly against the viewpoint that SMGs are extreme pathological objects, of little relevance in the cosmic history of star-formation.

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.

Gamma-ray burst (GRB) host galaxies have been studied extensively in optical photometry and spectroscopy. Here we present the first mid-infrared spectrum of a GRB host, HG031203. It is one of the nearest GRB hosts at z=0.1055, allowing both low and high-resolution spectroscopy with Spitzer-IRS. Medium resolution UV-to-K-band spectroscopy with the X-shooter spectrograph on the VLT is also presented, along with Spitzer IRAC and MIPS photometry, as well as radio and sub-mm observations. These data allow us to construct a UV-to-radio spectral energy distribution with almost complete spectroscopic coverage from 0.3-35 micron of a GRB host galaxy for the first time, potentially valuable as a template for future model comparisons. The IRS spectra show strong, high-ionisation fine structure line emission indicative of a hard radiation field in the galaxy, suggestive of strong ongoing star-formation and a very young stellar population. The selection of HG031203 via the presence of a GRB suggests that it might be a useful analogue of very young star-forming galaxies in the early universe, and hints that local BCDs may be used as more reliable analogues of star-formation in the early universe than typical local starbursts. We look at the current debate on the ages of the dominant stellar populations in z~7 and z~8 galaxies in this context. The nebular line emission is so strong in HG031203, that at z~7, it can reproduce the spectral energy distributions of z-band dropout galaxies with elevated IRAC 3.6 and 4.5 micron fluxes without the need to invoke a 4000A break.

GRBs generate an afterglow emission that can be detected from radio to X-rays during days, or even weeks after the initial explosion. The peak of this emission crosses the mm/submm range during the first hours to days, making their study in this range crucial for constraining the models. Observations have been limited until now due to the low sensitivity of the observatories in this range. We present observations of 10 GRB afterglows obtained from APEX and SMA, as well as the first detection of a GRB with ALMA, and put them into context with all the observations that have been published until now in the spectral range that will be covered by ALMA. The catalogue of mm/submm observations collected here is the largest to date and is composed of 102 GRBs, of which 88 had afterglow observations, whereas the rest are host galaxy searches. With our programmes, we contributed with data of 11 GRBs and the discovery of 2 submm counterparts. In total, the full sample, including data from the literature, has 22 afterglow detections with redshift ranging from 0.168 to 8.2. GRBs have been detected in mm/submm wavelengths with peak luminosities spanning 2.5 orders of magnitude, the most luminous reaching 10^33erg s^-1 Hz^-1. We observe a correlation between the X-ray brightness at 0.5 days and the mm/submm peak brightness. Finally we give a rough estimate of the distribution of peak flux densities of GRB afterglows, based on the current mm/submm sample. Observations in the mm/submm bands have been shown to be crucial for our understanding of the physics of GRBs, but have until now been limited by the sensitivity of the observatories. With the start of the operations at ALMA, the sensitivity will be increased by more than an order of magnitude. Our estimates predict that, once completed, ALMA will detect up to 98% of the afterglows if observed during the passage of the peak synchrotron emission.

To discriminate between different dust formation processes is a key issue in order to understand its properties. We analysed six submillimeter galaxies at redshifts 4<z<5 and nine quasars at 5<z<6.4. We estimated their dust masses from their (sub)millimeter emission and their stellar masses from the spectral energy distribution modelling or from the dynamical and gas masses obtained from the CO line detections. We calculated the dust yields per AGB star and per SN required to explain these dust masses and concluded that AGB stars are not efficient enough to form dust in the majority of these galaxies. SN could be responsible for dust production, but only if dust destruction in the SN shocks is not taken into account. Otherwise even SNe are not efficient enough, which advocates for some other dust production mechanism. We present the hypothesis that grain growth in the interstellar medium is responsible for bulk of the dust mass accumulation in these galaxies.

Q0151+048 is a physical QSO pair at z ~ 1.929 with a separation of 3.3 arcsec on the sky. In the spectrum of Q0151+048A (qA), a DLA is observed at a higher redshift. We have previously detected the host galaxies of both QSOs, as well as a Lya blob. We performed low-resolution spectroscopy with the slit aligned with the extended emission. We also observed the system using the medium-resolution VLT/X-shooter spectrograph and the slit aligned with the two QSOs. We measure systemic redshifts of zem(A)=1.92924{\pm}0.00036 and zem(B)=1.92863{\pm}0.00042 from the H{\beta} and H{\alpha} emission lines, respectively. We estimate the masses of the black holes of the two QSOs to be 10^9.33 M{\odot} and 10^8.38 M{\odot} for qA and qB, respectively. From this we infered the mass of the dark matter halos hosting the two QSOs: 10^13.74 M{\odot} and 10^13.13 M{\odot} for qA and qB, respectively. We observe a velocity gradient along the major axis of the Lya blob consistent with the rotation curve of a large disk galaxy, but it may also be caused by gas inflow or outflow. We detect residual continuum in the DLA trough which we interpret as emission from the host galaxy of qA. The derived H0 column density of the DLA is log NH0 = 20.34 {\pm} 0.02. Metal column densities results in an overall metallicity of 0.01 Z{\odot}. We detect CII* which allows us to make a physical model of the DLA cloud. From the systemic redshifts of the QSOs, we conclude that the Lya blob is associated with qA rather than with the DLA. The DLA must be located in front of both the Lya blob and qA at a distance larger than 30 kpc. The two QSOs accrete at normal eddington ratios. The DM halo of this double quasar will grow to the mass of our local super-cluster at z=0. We point out that those objects therefore form an ideal laboratory to study the physical interactions in a z=2 pre-cursor of our local super-cluster.

Aims: To gain better insight on the physics of stripped-envelope core-collapse supernovae through studying their environments. Methods: We obtained low-resolution optical spectroscopy with the New Technology Telescope (+ EFOSC2) at the locations of 20 Type Ib/c supernovae. We measure the flux of emission lines in the stellar-continuum-subtracted spectra from which local metallicities are computed. For the supernova regions we estimate both the mean stellar age, interpreting the stellar absorption with population synthesis models, and the age of the youngest stellar populations using the H-alpha equivalent width as an age indicator. These estimates are compared with the lifetimes of single massive stars. Results: Based on our sample, we detect a tentative indication that Type Ic supernovae might explode in environments that are more metal-rich than those of Type Ib supernovae (average difference of 0.08 dex), but this is not a statistically significant result. The lower limits placed on the ages of the supernova birthplaces are overall young, although there are several cases where these appear older than what is expected for the evolution of single stars more massive than 25-30 M_{sun}. This is only true, however, assuming that the supernova progenitors were born during an instantaneous (not continuous) episode of star formation. Conclusions: These results do not conclusively favor any of the two evolutionary paths (single or binary) leading to stripped supernovae. We do note a fraction of events for which binary evolution is more likely, due to their associated age limits. The fact, however, that the supernova environments contain areas of recent (< 15 Myr) star formation and that the environmental metallicities are, at least, not against the single evolutionary scenario, suggest that this channel is also broadly consistent with the observations.

We use the science demonstration field data of the Herschel-ATLAS to study how star formation, traced by the far-infrared Herschel data, is related to both the accretion luminosity and redshift of quasars selected from the Sloan Digital Sky Survey and the 2SLAQ survey. By developing a maximum likelihood estimator to investigate the presence of correlations between the far-infrared and optical luminosities we find evidence that the star-formation in quasar hosts is correlated with both redshift and quasar accretion luminosity. Assuming a relationship of the form L_IR \propto L_QSO^{\theta} (1 + z)^{\zeta}, we find {\theta} = 0.22 +/- 0.08 and {\zeta} = 1.6 +/- 0.4, although there is substantial additional uncertainty in {\zeta} of order +/- 1, due to uncertainties in the host galaxy dust temperature. We find evidence for a large intrinsic dispersion in the redshift dependence, but no evidence for intrinsic dispersion in the correlation between L_QSO and L_IR, suggesting that the latter may be due to a direct physical connection between star formation and black hole accretion. This is consistent with the idea that both the quasar activity and star formation are dependent on the same reservoir of cold gas, so that they are both affected by the influx of cold gas during mergers or heating of gas via feedback processes.

We present physical properties of two submillimeter selected gravitationally lensed sources, identified in the Herschel Astrophysical Terahertz Large Area Survey. These submillimeter galaxies (SMGs) have flux densities > 100 mJy at 500 um, but are not visible in existing optical imaging. We fit light profiles to each component of the lensing systems in Spitzer IRAC 3.6 and 4.5 um data and successfully disentangle the foreground lens from the background source in each case, providing important constraints on the spectral energy distributions (SEDs) of the background SMG at rest-frame optical-near-infrared wavelengths. The SED fits show that these two SMGs have high dust obscuration with Av ~4 to 5 and star formation rates of ~100 M_sun/yr. They have low gas fractions and low dynamical masses compared to 850 um selected galaxies.

We analyse Spitzer images of 30 long-duration gamma-ray burst (GRB) host galaxies. We estimate their total stellar masses (M*) based on the rest-frame K-band luminosities (L_Krest) and constrain their star formation rates (SFRs, not corrected for dust extinction) based on the rest-frame UV continua. Further, we compute a mean M*/L_Krest = 0.45 Msun/Lsun. We find that the hosts are low M*, star-forming systems. The median M* in our sample (<M*> = 10^9.7 Msun) is lower than that of "field" galaxies (e.g., Gemini Deep Deep Survey). The range spanned by M* is 10^7 Msun < M* < 10^11 Msun, while the range spanned by the dust-uncorrected UV SFR is 10^-2 Msun yr^-1 < SFR < 10 Msun yr^-1. There is no evidence for intrinsic evolution in the distribution of M* with redshift. We show that extinction by dust must be present in at least 25% of the GRB hosts in our sample and suggest that this is a way to reconcile our finding of a relatively lower UV-based, specific SFR (PHI = SFR/M*) with previous claims that GRBs have some of the highest PHI values. We also examine the effect that the inability to resolve the star-forming regions in the hosts has on PHI.

We investigate whether stellar dust sources i.e. asymptotic giant branch (AGB) stars and supernovae (SNe) can account for dust detected in 5<z<6.5 quasars (QSOs). We calculate the required dust yields per AGB star and per SN using the dust masses of QSOs inferred from their millimeter emission and stellar masses approximated as the difference between the dynamical and the H_2 gas masses of these objects. We find that AGB stars are not efficient enough to form dust in the majority of the z>5 QSOs, whereas SNe may be able to account for dust in some QSOs. However, they require very high dust yields even for a top-heavy initial mass function. This suggests additional non-stellar dust formation mechanism e.g. significant dust grain growth in the interstellar medium of at least three out of nine z>5 QSOs. SNe (but not AGB stars) may deliver enough heavy elements to fuel this growth.

We have determined the luminosity function of 250um-selected galaxies detected in the ~14 sq.deg science demonstration region of the Herschel-ATLAS project out to a redshift of z=0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250um luminosity density out to z=0.2 where it is 3.6+1.4-0.9 times higher than the local value.

We present a derivation of the star formation rate per comoving volume of quasar host galaxies, derived from stacking analyses of far-infrared to mm-wave photometry of quasars with redshifts 0<z<6 and absolute I-band magnitudes -22>I_AB>-32. We use the science demonstration observations of the first ~16 deg^2 from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) in which there are 240 quasars from the Sloan Digital Sky Survey (SDSS) and a further 171 from the 2dF-SDSS LRG and QSO (2SLAQ) survey. We supplement this data with a compilation of data from IRAS, ISO, Spitzer, SCUBA and MAMBO. H-ATLAS alone statistically detects the quasars in its survey area at >5sigma at 250, 350 and 500um. From the compilation as a whole we find striking evidence of downsizing in quasar host galaxy formation: low-luminosity quasars with absolute magnitudes in the range -22>I_AB>-24 have a comoving star formation rate (derived from 100um rest-frame luminosities) peaking between redshifts of 1 and 2, while high-luminosity quasars with I_AB<-26 have a maximum contribution to the star formation density at z~3. The volume-averaged star formation rate of -22>I_AB>-24 quasars evolves as (1+z)^{2.3 +/- 0.7} at z<2, but the evolution at higher luminosities is much faster reaching (1+z)^{10 +/- 1} at -26>I_AB>-28. We tentatively interpret this as a combination of a declining major merger rate with time and gas consumption reducing fuel for both black hole accretion and star formation.

The existence of submillimeter-selected galaxies (SMGs) at redshifts z>4 has recently been confirmed. Using simultaneously all the available data from UV to radio we have modeled the spectral energy distributions of the six known spectroscopically confirmed SMGs at z>4. We find that their star formation rates (average ~2500 MSun yr^{-1}), stellar (~3.6x10^{11} MSun) and dust (~6.7x10^{8} MSun) masses, extinction (A_V~2.2 mag), and gas-to-dust ratios (~60) are within the ranges for 1.7<z<3.6 SMGs. Our analysis suggests that infrared-to-radio luminosity ratios of SMGs do not change up to redshift ~5 and are lower by a factor of ~2.1 than the value corresponding to the local IR-radio correlation. However, we also find dissimilarities between z>4 and lower-redshift SMGs. Those at z>4 tend to be among the most star-forming, least massive and hottest (~60 K) SMGs and exhibit the highest fraction of stellar mass formed in the ongoing starburst (~45%). This indicates that at z>4 we see earlier stages of evolution of submillimeter-bright galaxies. Using the derived properties for z>4 SMGs we investigate the origin of dust at epochs less than 1.5 Gyr after the big bang. This is significant to our understanding of the evolution of the early universe. For three z>4 SMGs, asymptotic giant branch stars could be the dominant dust producers. However, for the remaining three only supernovae (SNe) are efficient and fast enough to be responsible for dust production, though requiring a very high dust yield per SN (0.15-0.65 MSun). The required dust yields are lower if a top-heavy initial mass function or significant dust growth in the interstellar medium is assumed. We estimate lower limits of the contribution of SMGs to the cosmic star formation and stellar mass densities at z~4-5 to be ~4% and ~1%, respectively.

Context: GRB afterglows are excellent probes of gas and dust in star-forming galaxies at all epochs. It has been posited that dust in the early Universe must be different from dust at lower z. To date two reports directly support this contention, one of which is based on the spectral shape of GRB 050904 at z = 6.295. Aims: We reinvestigate the afterglow to understand dust at high z. We address the claimed evidence for unusual (SN-origin) dust in its host galaxy by simultaneously examining the X-ray and optical/NIR spectrophotometric data. Methods: We derive the intrinsic SED of the afterglow at 0.47, 1.25 and 3.4 days, by re-reducing the Swift X-ray data, the 1.25 days FORS2 z-Gunn photometric data, the spectroscopic and z'-band photometric data at ~3 days from the Subaru telescope, as well as the critical UKIRT Z-band photometry at 0.47 days, upon which the claim of dust detection largely relies. Results: We find no evidence of dust extinction in the SED. We compute flux densities at lambda_rest = 1250 AA directly from the observed counts at all epochs. In the earliest epoch, 0.47 days, the Z-band suppression is found to be smaller (0.3 +- 0.2 mag) than previously reported and statistically insignificant (<1.5 sigma). Furthermore we find that the photometry of this band is unstable and difficult to calibrate. Conclusions: From the afterglow SED we demonstrate that there is no evidence for dust extinction -- the SED at all times can be reproduced without dust, and at 1.25 days in particular, significant extinction can be excluded, with A(3000 AA) < 0.27 mag at 95% confidence using the SN-type extinction curve. We conclude that there is no evidence of any extinction in the afterglow of GRB 050904 and that the presence of SN-origin dust in the host of GRB 050904 must be viewed skeptically. [abridged]

The nature of galaxies selected at submillimeter wavelengths (SMGs, S_850 > 3 mJy), some of the bolometrically most luminous objects at high redshifts, is still elusive. In particular their star formation histories and source of emission are not accurately constrained. In this paper we introduce a new approach to analyse the SMG data. Namely, we present the first self-consistent UV-to-radio spectral energy distribution fits of 76 SMGs with spectroscopic redshifts using all photometric datapoints from ultraviolet to radio simultaneously. We find that they are highly star-forming (median star formation rate 713 MSun yr^-1 for SMGs at z>0.5), moderately dust-obscured (median A_V~2 mag), hosting significant stellar populations (median stellar mass 3.7x10^11 MSun) of which only a minor part has been formed in the ongoing starburst episode. This implies that in the past, SMGs experienced either another starburst episode or merger with several galaxies. The properties of SMGs suggest that they are progenitors of present-day elliptical galaxies. We find that these bright SMGs contribute significantly to the cosmic star formation rate density (~20%) and stellar mass density (~30-50%) at redshifts 2-4. Using number counts at low fluxes we find that as much as 80% of the cosmic star formation at these redshifts took place in SMGs brighter than 0.1 mJy. We find evidence that a linear infrared-radio correlation holds for SMGs in an unchanged form up to redshift of 3.6, though its normalization is offset from the local relation by a factor of ~2.1 towards higher radio luminosities. We present a compilation of photometry data of SMGs and determinations of cosmic SFR and stellar mass densities.

We present an analysis of the spectral energy distribution (SED) of the galaxy ESO 184-G82, the host of the closest known long gamma-ray burst (GRB) 980425 and its associated supernova (SN) 1998bw. We use our observations obtained at the Australia Telescope Compact Array (the third >3 sigma radio detection of a GRB host) as well as archival infrared and ultraviolet (UV) observations to estimate its star formation state. We find that ESO 184-G82 has a UV star formation rate (SFR) and stellar mass consistent with the population of cosmological GRB hosts and of local dwarf galaxies. However, it has a higher specific SFR (per unit stellar mass) than luminous spiral galaxies. The mass of ESO 184-G82 is dominated by an older stellar population in contrast to the majority of GRB hosts. The Wolf-Rayet region ~800 pc from the SN site experienced a starburst episode during which the majority of its stellar population was built up. Unlike that of the entire galaxy, its SED is similar to those of cosmological submillimeter/radio-bright GRB hosts with hot dust content. These findings add to the picture that in general, the environments of GRBs on 1-3 kpc scales are associated with high specific SFR and hot dust.

SN 2008D was discovered while following up an unusually bright X-ray transient in the nearby spiral galaxy NGC 2770. We present early optical spectra (obtained 1.75 days after the X-ray transient) which allowed the first identification of the object as a supernova at redshift z = 0.007. These spectra were acquired during the initial declining phase of the light curve, likely produced in the stellar envelope cooling after shock breakout, and rarely observed. They exhibit a rather flat spectral energy distribution with broad undulations, and a strong, W-shaped feature with minima at 3980 and 4190 AA (rest frame). We also present extensive spectroscopy and photometry of the supernova during the subsequent photospheric phase. Unlike supernovae associated with gamma-ray bursts, SN 2008D displayed prominent He features and is therefore of Type Ib.

We present optical to radio spectral energy distribution fitting of the host galaxies of four long gamma-ray bursts: 980703, 000210, 000418 and 010222, which were detected at submillimetre and/or radio wavelengths. We find that only very young starburst galaxy models are consistent with the data having both blue optical colors and a pronounced submm emission. For each host we are able to construct a model consistent with the short- and long-wavelength parts of the spectra. We find galaxy ages ranging from 0.09 to 2.0 Gyrs and star formation rates ranging from 138 to 380 MSun/yr.

We present detailed fits of the spectral energy distributions (SEDs) of four submillimeter (submm) galaxies selected by the presence of a gamma-ray burst (GRB) event (GRBs 980703, 000210, 000418 and 010222). These faint ~3 mJy submm emitters at redshift ~1 are characterized by an unusual combination of long- and short-wavelength properties, namely enhanced submm and/or radio emission combined with optical faintness and blue colors. We exclude an active galactic nucleus as the source of long-wavelength emission. From the SED fits we conclude that the four galaxies are young (ages <2 Gyr), highly starforming (star formation rates ~150 MSun/yr), low-mass (stellar masses ~10^10 MSun) and dusty (dust masses ~3x10^8 MSun). Their high dust temperatures (Td>45 K) indicate that GRB host galaxies are hotter, younger, and less massive counterparts to submm-selected galaxies detected so far. Future facilities like Herschel, JCMT/SCUBA-2 and ALMA will test this hypothesis enabling measurement of dust temperatures of fainter GRB-selected galaxies.

Using detailed spectral energy distribution fits we present evidence that submillimeter- and radio-bright gamma-ray burst host galaxies are hotter counterparts to submillimeter galaxies. This hypothesis makes them of special interest since hotter submm galaxies are difficult to find and are believed to contribute significantly to the star formation history of the Universe.

We analyse 4.5, 8 and 24 um band Spitzer images of six gamma ray burst host galaxies at redshifts close to 1. We constrain their star formation rates (SFR) based on the entire available spectral energy distribution rather than the 24 um band only. Further, we estimate their stellar masses (M*) based on rest frame K band luminosities. Our sample spans a wide range of galaxy properties: derived SFRs range from less than 10 to a few hundred solar masses per year; values of M* range from 10^9 to 10^10 Mo with a median of 5.6 x 10^9 Mo. Comparing the specific star formation rate (PHI = SFR/M*) of our sample as a function of M* to other representative types of galaxies (distant red galaxies, Ly-alpha emitters, Lyman break galaxies, submillimeter galaxies and z ~ 2 galaxies from the Great Observatories Origins Deep Survey-North field), we find that gamma ray burst hosts are among those with the highest PHI.