Search

High-redshift protoclusters consisting of dusty starbursts are thought to play an important role in galaxy evolution. Their dusty nature makes them bright in the FIR/submm but difficult to find in optical/NIR surveys. Radio observations are an excellent way to study these dusty starbursts, as dust is transparent in the radio and there is a tight correlation between the FIR and radio emission of a galaxy. Here, we present MeerKAT 1.28 GHz radio imaging of 3 Herschel candidate protoclusters, with a synthesised beam size of ~$7.5'' \times 6.6''$ and a central thermal noise down to $4.35~\mu$Jy/beam. Our source counts are consistent with other radio counts with no evidence of overdensities. Around 95% of the Herschel sources have 1.28 GHz IDs. Using the Herschel 250 micron primary beam size as the searching radius, we find 54.2% Herschel sources have multiple 1.28 GHz IDs. Our average FIR-radio correlation coefficient $q_{250\mu m}$ is $2.33\pm 0.26$. Adding $q_{250\mu m}$ as a new constraint, the probability of finding chance-aligned sources is reduced by a factor of ~6, but with the risk of discarding true identifications of radio-loud/quiet sources. With accurate MeerKAT positions, we cross-match our Herschel sources to optical/NIR data followed by photometric redshift estimations. By removing $z<1$ sources, the density contrasts of two of the candidate protoclusters increase, suggestive of them being real protoclusters at $z>1$. There is also potentially a $0.9

One of the greatest challenges in galaxy evolution over the last decade has been constraining the prevalence of heavily dust-obscured galaxies in the early Universe. At $z>3$, these galaxies are increasingly rare, and difficult to identify as they are interspersed among the more numerous dust-obscured galaxy population at $z=1-3$, making efforts to secure confident spectroscopic redshifts expensive, and sometimes unsuccessful. In this work, we present the Extended Mapping Obscuration to Reionization with ALMA (Ex-MORA) Survey -- a 2mm blank-field survey in the COSMOS-Web field, and the largest ever ALMA blank-field survey to-date covering 577 arcmin$^2$. Ex-MORA is an expansion of the MORA survey designed to identify primarily $z>3$ dusty, star-forming galaxies while simultaneously filtering out the more numerous $z<3$ population by leveraging the very negative $K$-correction at observed-frame 2mm. We identify 37 significant ($>$5$\sigma$) sources, 33 of which are robust thermal dust emitters. We measure a median redshift of $\langle z \rangle = 3.6^{+0.1}_{-0.2}$, with two-thirds of the sample at $z>3$, and just under half at $z>4$, demonstrating the overall success of the 2mm-selection technique. The integrated $z>3$ volume density of Ex-MORA sources is $\sim1-3\times10^{-5}$ Mpc$^{-3}$, consistent with other surveys of infrared luminous galaxies at similar epochs. We also find that techniques using rest-frame optical emission (or lack thereof) to identify $z>3$ heavily dust-obscured galaxies miss at least half of Ex-MORA galaxies. This supports the idea that the dusty galaxy population is heterogeneous, and that synergies across observatories spanning multiple energy regimes are critical to understanding their formation and evolution at $z>3$., Comment: 13 pages, 5 figures, submitted to ApJ; fully reduced mosaic will be shared upon publication

In this Roadmap, we present a vision for the future of submillimetre and millimetre astronomy in the United Kingdom over the next decade and beyond. This Roadmap has been developed in response to the recommendation of the Astronomy Advisory Panel (AAP) of the STFC in the AAP Astronomy Roadmap 2022. In order to develop our stragetic priorities and recommendations, we surveyed the UK submillimetre and millimetre community to determine their key priorities for both the near-term and long-term future of the field. We further performed detailed reviews of UK leadership in submillimetre/millimetre science and instrumentation. Our key strategic priorities are as follows: 1. The UK must be a key partner in the forthcoming AtLAST telescope, for which it is essential that the UK remains a key partner in the JCMT in the intermediate term. 2. The UK must maintain, and if possible enhance, access to ALMA and aim to lead parts of instrument development for ALMA2040. Our strategic priorities complement one another: AtLAST (a 50m single-dish telescope) and an upgraded ALMA (a large configurable interferometric array) would be in synergy, not competition, with one another. Both have identified and are working towards the same overarching science goals, and both are required in order to fully address these goals., Comment: 91 pages plus cover, 38 figures. Submitted to the Science and Technology Facilities Council, August 2024. One figure corrected (v2); new appendix with STFC Q&A; corrected SMA access statement; updated references, acronyms & author list (v3)

Herschel surveys have found large numbers of sources with red far-IR colours, and spectral energy distributions (SEDs) rising from 250 to 500$\mu$m: 500 risers. The nature and role of these sources is not fully understood. We here present Submillimeter Array (SMA) interferometric imaging at 200 GHz of a complete sample of five 500 risers with F500 $>$ 44 mJy selected within a 4.5 square degree region of the XMMLSS field. These observations can resolve the separate components of multiple sources and allow cross identification at other wavelengths using the extensive optical-to-IR data in this field. Of our five targets we find that two are likely gravitationally lensed, two are multiple sources, and one an isolated single source. Photometric redshifts, using optical-to-IR data and far-IR/submm data, suggest they lie at redshifts $z \sim 2.5 - 3.5$. Star formation rates and stellar masses estimated from the SEDs show that the majority of our sources lie on the star-formation rate-stellar mass `main sequence', though with outliers both above and below this relation. Of particular interest is our most multiple source which consists of three submm emitters and one submm-undetected optical companion within a 7 arcsecond region, all with photometric redshifts $\sim$ 3. One of the submm emitters in this group lies above the `main sequence', while the optical companion lies well below the relation, and has an estimated stellar mass of 3.3$\pm 1.3 \times$10$^{11}$ M$_{\odot}$. We suggest this object is a forming brightest cluster galaxy (BCG) in the process of accreting actively star forming companions., Comment: Accepted for publication in MNRAS

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance., Comment: Accepted for publication in the A&A special issue`Euclid on Sky'

We present 850um imaging of the XMM-LSS field observed for 170 hours as part of the James Clerk Maxwell Telescope SCUBA-2 Large eXtragalactic Survey (S2LXS). S2LXS XMM-LSS maps an area of 9 square degrees, reaching a moderate depth of 1-sigma ~ 4 mJy/beam. This is the largest contiguous area of extragalactic sky mapped by JCMT at 850um to date. The wide area of the S2LXS XMM-LSS survey allows us to probe the ultra-bright (S_850um > 15 mJy), yet rare submillimetre population. We present the S2LXS XMM-LSS catalogue, which comprises 40 sources detected at >5-sigma significance, with deboosted flux densities in the range of 7 mJy to 48 mJy. We robustly measure the bright-end of the 850um number counts at flux densities >7 mJy, reducing the Poisson errors compared to existing measurements. The S2LXS XMM-LSS observed number counts show the characteristic upturn at bright fluxes, expected to be motivated by local sources of submillimetre emission and high-redshift strongly lensed galaxies. We find that the observed 850um number counts are best reproduced by model predictions that include either strong lensing or source blending from a 15 arcsec beam, indicating that both may make an important contribution to the observed over-abundance of bright single dish 850um selected sources. We make the S2LXS XMM-LSS 850um map and >5-sigma catalogue presented here publicly available., Comment: Accepted for publication in MNRAS

We investigate the resolved kinematics of the molecular gas, as traced by ALMA in CO (2-1), of 25 cluster member galaxies across three different clusters at a redshift of $z\sim1.6$. This is the first large-scale analysis of the molecular gas kinematics of cluster galaxies at this redshift. By separately estimating the rotation curve of the approaching and receding side of each galaxy via kinematic modeling, we quantify the difference in total circular velocity to characterize the overall kinematic asymmetry of each galaxy. 3/14 of the galaxies in our sample that we are able to model have similar degrees of asymmetry as that observed in galaxies in the field at similar redshift. However, this leaved 11/14 galaxies in our sample with significantly higher asymmetry, and some of these galaxies have degrees of asymmetry of up to $\sim$50 times higher than field galaxies observed at similar redshift. Some of these extreme cases also have one-sided tail-like morphology seen in the molecular gas, supporting a scenario of tidal and/or ram pressure interaction. Such stark differences in the kinematic asymmetry in clusters versus the field suggest the evolutionary influence of dense environments, established as being a major driver of galaxy evolution at low-redshift, is also active in the high-redshift universe., Comment: Submission pre-referee report

Protoclusters are important for studying how halo mass and stellar mass assemble in the early universe. Finding signposts of such over-dense regions is a popular method to identify protocluster candidates. Hyperluminous infrared galaxies (HLIRGs), are expected to reside in overdense regions with massive halos. We study the Mpc-scale environment of the largest HLIRG sample to date and investigate whether they predominantly live in overdense regions. We first explore the surface density of Herschel 250 $\mu$m sources around HLIRGs and compare with that around random positions. Then, we compare the spatial distribution of neighbours around HLIRGs with that around randomly selected galaxies using a deep IRAC-selected catalogue with good-quality photometric redshifts. We also use a redshift-matched quasar sample and submillimeter galaxy (SMG) sample to validate our method, as previous clustering studies have measured the host halo masses of these populations. Finally, we adopt a Friends of Friends (FOF) algorithm to seek (proto)clusters that host HLIRGs. We find that HLIRGs tend to have more bright star-forming neighbours (with 250 $\mu$m flux density >10 mJy) within 100$\arcsec$ projected radius than a random galaxy at a 3.7$\sigma$ significance. In our 3D analysis, we find relatively weak excess of IRAC-selected sources within 3 Mpc around HLIRGs compared with random galaxy neighbours, mainly influenced by photometric redshift uncertainty and survey depth. We find a more significant difference (at a 4.7$\sigma$ significance) in the number of Low Frequency Array (LOFAR)-detected neighbours in the deepest EN1 field. HLIRGs at 3 < z < 4 show stronger excess compared to HLIRGs at 2 < z < 3, consistent with cosmic downsizing. Finally, we select and present a list of 30 most promising protocluster candidates for future follow-up observations., Comment: 19 pages, 10 figures, accepted by A&A

We examine the far-IR properties of a sample of 5391 optically selected QSOs in the 0.544.7, using SPIRE data from Herschel-ATLAS. We split the sample in a grid of 74 luminosity-redshift bins and compute the average optical-infrared spectral energy distribution (SED) in each bin. By normalising an intrinsic AGN template to the AGN optical power (at 5100A) we decompose the total infrared emission (L_IR; 8-1000um) into an AGN (L_IR,AGN) and star-forming component (L_IR,SF). We find that the AGN contribution to L_IR increases as a function of AGN power which manifests as a reduction of the `far-IR bump' in the average QSO SEDs. We note that L_IR,SF does not correlate with AGN power; the mean star formation rates (SFRs) of AGN host galaxies are a function of redshift only and they range from ~6 Msun/yr at z~0 to a plateau of <200 Msun/yr at z~2.6. Our results indicate that the accuracy of far-IR emission as a proxy for SFR decreases with increasing AGN luminosity. We show that, at any given redshift, observed trends between infrared luminosity (whether monochromatic or total) and AGN power (in the optical or X-rays) can be explained by a simple model which is the sum of two components: (A) the infrared emission from star-formation, uncorrelated with AGN power and (B) the infrared emission from AGN, directly proportional to AGN power in the optical or X-rays., Comment: 20 pages, 10 figures, accepted for publication in MNRAS

We present the ultraviolet to submillimetre spectral energy distributions (SEDs) of the HERschel Ultra Luminous Infrared Galaxy Survey (HERUS) sample of 42 local ultraluminous infrared galaxies (ULIRGs) and fit them with a Markov chain Monte Carlo (MCMC) code using the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) radiative transfer models for starbursts, active galactic nucleus (AGN) tori and host galaxy. The Spitzer IRS spectroscopy data are included in the fitting. Our bayesian SED fitting method takes comparable time to popular energy balance methods but it is more physically motivated and versatile. All HERUS galaxies harbor high rates of star formation but we also find bolometrically significant AGN in all of the galaxies of the sample. We estimate the correction of the luminosities of the AGN in the ULIRGs due to the anisotropic emission of the torus and find that it could be up to a factor of $\sim10$ for nearly edge-on tori. We present a comparison of our results with the smooth torus model of Fritz et al. and the two-phase models of Siebenmorgen et al. and SKIRTOR. We find that the CYGNUS AGN torus models fit significantly better the SEDs of our sample compared to all other models. We find no evidence that strong AGN appear either at the beginning or end of a starburst episode or that starbursts and AGN affect each other. IRAS 01003-2238 and Mrk 1014 show evidence for dual AGN in their SED fits suggesting a minimum dual AGN fraction in the sample of 5%., Comment: Published in Monthly Notices of the Royal Astronomical Society. The CYGNUS models used in this paper are available at https://arc.euc.ac.cy/cygnus/

We present the multi-wavelength counterparts of 850-$\mu$m selected submillimetre sources over a 2-deg$^2$ field centred on the North Ecliptic Pole. In order to overcome the large beam size (15 arcsec) of the 850-$\mu$m images, deep optical to near-infrared (NIR) photometric data and arcsecond-resolution 20-cm images are used to identify counterparts of submillimetre sources. Among 647 sources, we identify 514 reliable counterparts for 449 sources (69 per cent in number), based either on probabilities of chance associations calculated from positional offsets or offsets combined with the optical-to-NIR colours. In the radio imaging, the fraction of 850-$\mu$m sources having multiple counterparts is 7 per cent. The photometric redshift, infrared luminosity, stellar mass, star-formation rate (SFR), and the AGN contribution to the total infrared luminosity of the identified counterparts are investigated through spectral energy distribution fitting. The SMGs are infrared-luminous galaxies at an average $\langle z\rangle=2.5$ with $\mathrm{log}_{10} (L_\mathrm{IR}/\mathrm{L}_\odot)=11.5-13.5$, with a mean stellar mass of $\mathrm{log}_{10} (M_\mathrm{star}/\mathrm{M}_\odot)=10.90$ and SFR of $\mathrm{log}_{10} (\mathrm{SFR/M_\odot\,yr^{-1}})=2.34$. The SMGs show twice as large SFR as galaxies on the star-forming main sequence, and about 40 per cent of the SMGs are classified as objects with bursty star formation. At $z\ge4$, the contribution of AGN luminosity to total luminosity for most SMGs is larger than 30 per cent. The FIR-to-radio correlation coefficient of SMGs is consistent with that of main-sequence galaxies at $z\simeq2$., Comment: 24 pages, 16 figures. Accepted for publication in MNRAS

We present SCUBA-2 and SMA follow-up observations of four candidate high redshift Dusty Star-Forming Galaxies, selected as sources with rising SEDs in the 250, 350 and 500$\mu$m Herschel SPIRE bands. Previous SMA observations showed no counterparts to these sources, but in our deeper sub-mm observations we detect counterparts to all four 500$\mu$m risers, with three resolving into multiple systems. For these three multiple systems, the SMA 345GHz ($\approx 870\mu$m) observations recover $123 \pm 73\%$, $60 \pm 15\%$ and $19 \pm 4\%$ respectively of the integrated 850$\mu$m flux density from SCUBA-2, indicating that there may be additional sources below our SMA detection limit making up a dense, protocluster core. The fourth 500$\mu$m riser was observed at a lower frequency and so we cannot make a similar comparison. We estimate photometric redshifts based on FIR/sub-mm colours, finding that 3/4 likely lie at $z \geq 2$. This fits with the interpretation that the 500$\mu$m riser selection criterion selects both intrinsically red, individual galaxies at $z > 4$, and multiple systems at more moderate redshifts, artificially reddened by the effects of blending. We use the SCUBA-2 850$\mu$m maps to investigate the environments of these 500$\mu$m risers. By constructing cumulative number counts and estimating photometric redshifts for surrounding SCUBA-2 detections, we find that one of our 500$\mu$m risers could plausibly reside in a $z \geq 2$ protocluster. We infer that bright 500$\mu$m risers with faint 850$\mu$m flux densities are typically multiple systems at $z \geq 2$ that may reside in overdensities of bright sub-mm galaxies., Comment: 20 pages, 13 figures, 3 tables. Submitted to MNRAS

The luminosity function (LF) of active galactic nuclei (AGN) probes the history of supermassive black hole assembly and growth across cosmic time. To mitigate selection biases, we present a consistent analysis of the AGN LFs derived for both X-ray and mid-infrared (MIR) selected AGN in the XMM-Large Scale Structure (XMM-LSS) field. There are 4268 AGN used to construct the MIR luminosity function (IRLF) and 3427 AGN used to construct the X-ray luminosity function (XLF), providing the largest census of the AGN population out to $z=4$ in both bands with significant reduction in uncertainties. We are able for the first time to see the knee of the IRLF at $z>2$ and observe a flattening of the faint-end slope as redshift increases. The bolometric luminosity density, a proxy for the cosmic black hole accretion history, computed from our LFs shows a peak at $z\approx2.25$ consistent with recent estimates of the peak in the star formation rate density (SFRD). However, at earlier epochs, the AGN luminosity density is flatter than the SFRD. If confirmed, this result suggests that the build up of black hole mass outpaces the growth of stellar mass in high mass systems at $z\gtrsim 2.5$. This is consistent with observations of redshift $z\sim 6$ quasars which lie above the local $M-\sigma$ relationship. The luminosity density derived from the IRLF is higher than that from the XLF at all redshifts. This is consistent with the dominant role of obscured AGN activity in the cosmic growth of supermassive black holes., Comment: 17 pages, 14 figures, 1 appendix

We present the first results from the 2mm Mapping Obscuration to Reionization (MORA) survey, the largest ALMA contiguous blank-field survey to-date with a total area of 184 sq. arcmin and the only at 2mm to search for dusty star-forming galaxies (DSFGs). We use the 13 sources detected above 5sigma to estimate the first ALMA galaxy number counts at this wavelength. These number counts are then combined with the state-of-the-art galaxy number counts at 1.2mm and 3mm and with a backward evolution model to place constraints on the evolution of the IR luminosity function and dust-obscured star formation in the last 13 billion years. Our results suggest a steep redshift evolution on the space density of DSFGs and confirm the flattening of the IR luminosity function at faint luminosities, with a slope of $\alpha_{LF} = -0.42^{+0.02}_{-0.04}$. We conclude that the dust-obscured component, which peaks at z=2-2.5, has dominated the cosmic history of star formation for the past ~12 billion years, back to z~4. At z=5, the dust-obscured star formation is estimated to be ~35% of the total star formation rate density and decreases to 25%-20% at z=6-7, implying a minor contribution of dust-enshrouded star formation in the first billion years of the Universe. With the dust-obscured star formation history constrained up to the end of the epoch of reionization, our results provide a benchmark to test galaxy formation models, to study the galaxy mass assembly history, and to understand the dust and metal enrichment of the Universe at early times., Comment: Accepted for publication in The Astrophysical Journal

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will observe several Deep Drilling Fields (DDFs) to a greater depth and with a more rapid cadence than the main survey. In this paper, we describe the ``DeepDrill'' survey, which used the Spitzer Space Telescope Infrared Array Camera (IRAC) to observe three of the four currently defined DDFs in two bands, centered on 3.6 $\mu$m and 4.5 $\mu$m. These observations expand the area which was covered by an earlier set of observations in these three fields by the Spitzer Extragalactic Representative Volume Survey (SERVS). The combined DeepDrill and SERVS data cover the footprints of the LSST DDFs in the Extended Chandra Deep Field-South field (ECDFS), the ELAIS-S1 field (ES1), and the XMM Large-Scale Structure Survey field (XMM-LSS). The observations reach an approximate $5\sigma$ point-source depth of 2 $\mu$Jy (corresponding to an AB magnitude of 23.1; sufficient to detect a 10$^{11} M_{\odot}$ galaxy out to $z\approx 5$) in each of the two bands over a total area of $\approx 29\,$deg$^2$. The dual-band catalogues contain a total of 2.35 million sources. In this paper we describe the observations and data products from the survey, and an overview of the properties of galaxies in the survey. We compare the source counts to predictions from the SHARK semi-analytic model of galaxy formation. We also identify a population of sources with extremely red ([3.6]$-$[4.5] $>1.2$) colours which we show mostly consists of highly-obscured active galactic nuclei., Comment: 20 pages, 9 figures; MNRAS in press

We analyse the physical properties of 121 SNR $\geq$ 5 sub-millimetre galaxies (SMGs) from the STUDIES 450-$\mu$m survey. We model their UV-to-radio spectral energy distributions using MAGPHYS+photo-$z$ and compare the results to similar modelling of 850-$\mu$m-selected SMG sample from AS2UDS, to understand the fundamental physical differences between the two populations at the observed depths. The redshift distribution of the 450-$\mu$m sample has a median of $z$ = 1.85 $\pm$ 0.12 and can be described by strong evolution of the far-infrared luminosity function. The fainter 450-$\mu$m sample has $\sim$14 times higher space density than the brighter 850-$\mu$m sample at $z$ $\lesssim$2, and a comparable space density at $z$ = 2-3, before rapidly declining, suggesting LIRGs are the main obscured population at $z$ $\sim$ 1-2, while ULIRGs dominate at higher redshifts. We construct rest-frame $\sim$ 180-$\mu$m-selected and dust-mass-matched samples at $z$ = 1-2 and $z$ = 3-4 from the 450-$\mu$m and 850-$\mu$m samples, respectively, to probe the evolution of a uniform sample of galaxies spanning the cosmic noon era. Using far-infrared luminosity, dust masses and an optically-thick dust model, we suggest that higher-redshift sources have higher dust densities due to inferred dust continuum sizes which are roughly half of those for the lower-redshift population at a given dust mass, leading to higher dust attenuation. We track the evolution in the cosmic dust mass density and suggest that the dust content of galaxies is governed by a combination of both the variation of gas content and dust destruction timescale., Comment: 20 pages, 7 figures, MNRAS accepted

We present SMA observations at resolutions from 0.35 to 3 arcseconds of a sample of 34 candidate high redshift dusty star forming galaxies (DSFGs). These sources were selected from the HerMES Herschel survey catalogues to have SEDs rising from 250 to 350 to 500$\mu$m, a population termed 500-risers. We detect counterparts to 24 of these sources, with four having two counterparts. We conclude that the remaining ten sources that lack detected counterparts are likely to have three or more associated sources which blend together to produce the observed Herschel source. We examine the role of lensing, which is predicted to dominate the brightest (F500 $>$ 60mJy) half of our sample. We find that while lensing plays a role, at least 35% of the bright sources are likely to be multiple sources rather than the result of lensing. At fainter fluxes we find a blending rate comparable to, or greater than, the predicted 40%. We determine far-IR luminosities and star formation rates for the non-multiple sources in our sample and conclude that, in the absence of strong lensing, our 500-risers are very luminous systems with L$_{FIR} > 10^{13}$L$_{\odot}$ and star formation rates $> 1000$M$_{\odot}$/yr., Comment: Accepted for publication in MNRAS

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., Comment: 20 Pages, 8 Figures + 7 Appendix Figures + 1 animated figure. Resubmitted in order to reflect the version as printed in MNRAS

Scaling laws of dust, HI gas and metal mass with stellar mass, specific star formation rate and metallicity are crucial to our understanding of the buildup of galaxies through their enrichment with metals and dust. In this work, we analyse how the dust and metal content varies with specific gas mass ($M_{\text{HI}}$/$M_{\star}$) across a diverse sample of 423 nearby galaxies. The observed trends are interpreted with a set of Dust and Element evolUtion modelS (DEUS) - incluidng stellar dust production, grain growth, and dust destruction - within a Bayesian framework to enable a rigorous search of the multi-dimensional parameter space. We find that these scaling laws for galaxies with $-1.0\lesssim \log M_{\text{HI}}$/$M_{\star}\lesssim0$ can be reproduced using closed-box models with high fractions (37-89$\%$) of supernova dust surviving a reverse shock, relatively low grain growth efficiencies ($\epsilon$=30-40), and long dus lifetimes (1-2\,Gyr). The models have present-day dust masses with similar contributions from stellar sources (50-80\,$\%$) and grain growth (20-50\,$\%$). Over the entire lifetime of these galaxies, the contribution from stardust ($>$90\,$\%$) outweighs the fraction of dust grown in the interstellar medium ($<$10$\%$). Our results provide an alternative for the chemical evolution models that require extremely low supernova dust production efficiencies and short grain growth timescales to reproduce local scaling laws, and could help solving the conundrum on whether or not grains can grow efficiently in the interstellar medium., Comment: 56 pages, 30 figures, accepted for publication in MNRAS

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., Comment: 8 pages, 1 figure, 1 table, accepted for publication in MNRAS

Quality early STEM experiences provide a critical foundation for learning about science, technology, engineering, and mathematics (STEM) in ways that facilitate later learning. Early childhood educators have a unique opportunity to benefit all students by bringing the STEM disciplines into greater balance with literacy, reducing early inequities in access and opportunity to learn STEM, and developing skills and understandings that students can build on as they progress through school and into civic life and the workforce. This brief draws on research and development supported by the National Science Foundation to highlight important considerations about STEM educational experiences for young children and professional learning for educators who provide those experiences.

The cosmic infrared background (CIB) provides a fundamental observational constraint on the star-formation history of galaxies over cosmic history. We estimate the contribution to the CIB from catalogued galaxies in the COSMOS field by using a novel map fitting technique on the \textit{Herschel} SPIRE maps. Prior galaxy positions are obtained using detections over a large range in wavelengths in the $K_{\rm s}$--3\,GHz range. Our method simultaneously fits the galaxies, the system foreground, and the leakage of flux from galaxies located in masked areas and corrects for an "over-fitting" effect not previously accounted for in stacking methods. We explore the contribution to the CIB as a function of galaxy survey wavelength and depth. We find high contributions to the CIB with the deep $r$ ($m_{\rm AB} \le 26.5$), $K_{\rm s}$ ($m_{\rm AB} \le 24.0$) and 3.6\,$\mu$m ($m_{\rm AB} \le 25.5$) catalogues. We combine these three deep catalogues and find a total CIB contributions of 10.5 $\pm$ 1.6, 6.7 $\pm$ 1.5 and 3.1 $\pm$ 0.7\,nWm$^{-2}$sr$^{-1}$ at 250, 350 and 500\,$\mu$m, respectively. Our CIB estimates are consistent with recent phenomenological models, prior based SPIRE number counts and with (though more precise than) the diffuse total measured by FIRAS. Our results raise the interesting prospect that the CIB contribution at $\lambda \le 500\,\mu$m from known galaxies has converged. Future large-area surveys like those with the Large Synoptic Survey Telescope are therefore likely to resolve a substantial fraction of the population responsible for the CIB at 250\,$\mu$m $\leq \lambda \leq$ 500\,$\mu$m., Comment: 16 pages, 14 figures

Dusty star-forming galaxies (DSFGs) detected at $z > 4$ provide important examples of the first generations of massive galaxies. However, few examples with spectroscopic confirmation are currently known, with Hershel struggling to detect significant numbers of $z > 6$ DSFGs. NGP6_D1 is a bright 850 $\mu m$ source (12.3 $\pm$ 2.5 mJy) with no counterparts at shorter wavelengths (a SPIRE dropout). Interferometric observations confirm it is a single source, with no evidence for any optical or NIR emission, or nearby likely foreground lensing sources. No $>3\sigma$ detected lines are seen in both LMT RSR and IRAM 30m EMIR spectra of NGP6_D1 across 32 $GHz$ of bandwidth despite reaching detection limits of $\sim 1 mJy/500 km~s^{-1}$, so the redshift remains unknown. Template fitting suggests that NGP6_D1 is most likely between $z = 5.8$ and 8.3. SED analysis finds that NGP6_D1 is a ULIRG, with a dust mass $\sim 10^8$ - $10^9$ $M_{\odot}$ and a SFR of $\sim$ 500 $M_{\odot}~yr^{-1}$. We place upper limits on the gas mass of NGP6_D1 of $M_{H2}$ $ < (1.1~\pm~3.5) \times 10^{11}$ $M_{\odot}$, consistent with a gas-to-dust ratio of $\sim$ 100 - 1000. We discuss the nature of NGP6_D1 in the context of the broader submm population, and find that comparable SPIRE dropouts account for $\sim$ 20% of all SCUBA-2 detected sources, but with a similar flux density distribution to the general population., Comment: Accepoted for publication in MNRAS

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 $22$ 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., Comment: 30 pages, 8 figures, 5 tables, accepted for publication in MNRAS

We present the discovery of HLock01-LAB, a luminous and large Lya nebula at z=3.326. Medium-band imaging and long-slit spectroscopic observations with the Gran Telescopio Canarias reveal extended emission in the Lya 1215\AA, CIV1550\AA, and HeII 1640\AA lines over ~100kpc, and a total luminosity L(Lya)=(6.4+/-0.1)x10^44 erg s^-1. HLock01-LAB presents an elongated morphology aligned with two faint radio sources contained within the central ~8kpc of the nebula. The radio structures are consistent to be faint radio jets or lobes of a central galaxy, whose spectrum shows nebular emission characteristic of a type-II active galactic nucleus (AGN). The continuum emission of the AGN at short wavelengths is, however, likely dominated by stellar emission of the host galaxy, for which we derive a stellar mass M* = 2.3x10^11 Msun. The detection of extended emission in CIV and CIII] indicates that the gas within the nebula is not primordial. Feedback may have enriched the halo at at least 50 kpc from the nuclear region. Using rest-frame UV emission-line diagnostics, we find that the gas in the nebula is likely heated by the AGN. Nevertheless, at the center of the nebula we find extreme emission line ratios of Lya/CIV~60 and Lya/HeII~80, one of the highest values measured to date, and well above the standard values of photoionization models (Lya/HeII~30 for case B photoionization). Our data suggest that jet-induced shocks are likely responsible for the increase of the electron temperature and, thus, the observed Lya enhancement in the center of the nebula. This scenario is further supported by the presence of radio structures and perturbed kinematics in this region. The large Lya luminosity in HLock01-LAB is likely due to a combination of AGN photoionization and jet-induced shocks, highlighting the diversity of sources of energy powering Lya nebulae. [abridged], Comment: 13 pages, 8 figures, 3 tables. Accepted for publication in A&A. Abstract abridged for arXiv submission

ALMA 870$\mu$m continuum imaging has uncovered a population of blends of multiple dusty star-forming galaxies (DSFGs) in sources originally detected with the Herschel Space Observatory. However, their pairwise separations are much smaller that what is found by ALMA follow-up of other single-dish surveys or expected from theoretical simulations. Using ALMA and VLA, we have targeted three of these systems to confirm whether the multiple 870$\mu$m continuum sources lie at the same redshift, successfully detecting $^{12}$CO($J = 3$-2) and $^{12}$CO($J = 1$-0) lines and being able to confirm that in the three cases all the multiple DSFGs are likely physically associated within the same structure. Therefore, we report the discovery of two new gas-rich dusty protocluster cores (HELAISS02, $z = 2.171 \pm 0.004$; HXMM20, $z = 2.602 \pm 0.002$). The third target is located in the well known COSMOS overdensity at $z = 2.51$ (named CL J1001+0220 in the literature), for which we do not find any new secure CO(1-0) detection, although some of its members show only tentative detections and require further confirmation. From the gas, dust, and stellar properties of the two new protocluster cores, we find very large molecular gas fractions yet low stellar masses, pushing the sources above the main sequence, while not enhancing their star formation efficiency. We suggest that the sources might be newly formed galaxies migrating to the main sequence. The properties of the three systems compared to each other and to field galaxies may suggest a different evolutionary stage between systems., Comment: Accepted by ApJ. 24 pages, 9 figures

We present the science case for mapping several thousand galaxy (proto)clusters at z=1-10 with a large aperture single dish sub-mm facility, producing a high-redshift counterpart to local large surveys of rich clusters like the well-studied Abell catalogue. Principal goals of a large survey of distant clusters are the evolution of galaxy clusters over cosmic time and the impact of environment on the evolution and formation of galaxies. To make a big leap forward in this emerging research field, the community would benefit from a large-format, wide-band, direct-detection spectrometer (e.g., based on MKID technology), covering a wide field of ~1 square degree and a frequency coverage from 70 to 700 GHz., Comment: Science white paper submitted to the Astro2020 Decadal Survey

With the recent discovery of a dozen dusty star-forming galaxies and around 30 quasars at z>5 that are hyper-luminous in the infrared ($\mu$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $\mu$ is a lensing magnification factor), the possibility has opened up for SPICA, the proposed ESA M5 mid-/far-infrared mission, to extend its spectroscopic studies toward the epoch of reionization and beyond. In this paper, we examine the feasibility and scientific potential of such observations with SPICA's far-infrared spectrometer SAFARI, which will probe a spectral range (35-230 $\mu$m) that will be unexplored by ALMA and JWST. Our simulations show that SAFARI is capable of delivering good-quality spectra for hyper-luminous infrared galaxies (HyLIRGs) at z=5-10, allowing us to sample spectral features in the rest-frame mid-infrared and to investigate a host of key scientific issues, such as the relative importance of star formation versus AGN, the hardness of the radiation field, the level of chemical enrichment, and the properties of the molecular gas. From a broader perspective, SAFARI offers the potential to open up a new frontier in the study of the early Universe, providing access to uniquely powerful spectral features for probing first-generation objects, such as the key cooling lines of low-metallicity or metal-free forming galaxies (fine-structure and H2 lines) and emission features of solid compounds freshly synthesized by Population III supernovae. Ultimately, SAFARI's ability to explore the high-redshift Universe will be determined by the availability of sufficiently bright targets (whether intrinsically luminous or gravitationally lensed). With its launch expected around 2030, SPICA is ideally positioned to take full advantage of upcoming wide-field surveys such as LSST, SKA, Euclid, and WFIRST, which are likely to provide extraordinary targets for SAFARI., Comment: 24 pages, 8 figures, accepted for publication in PASA (SPICA special issue)

Tidal disruption events (TDEs) are transient flares produced when a star is ripped apart by the gravitational field of a supermassive black hole (SMBH). We have observed a transient source in the western nucleus of the merging galaxy pair Arp 299 that radiated >1.5x10^52 erg in the infrared and radio, but was not luminous at optical or X-ray wavelengths. We interpret this as a TDE with much of its emission re-radiated at infrared wavelengths by dust. Efficient reprocessing by dense gas and dust may explain the difference between theoretical predictions and observed luminosities of TDEs. The radio observations resolve an expanding and decelerating jet, probing the jet formation and evolution around a SMBH., Comment: Published online by Science on 14 Jun 2018

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., Comment: 14 pages, 7 figures, accepted for publication in MNRAS

High-redshift, luminous, dusty star forming galaxies (DSFGs) constrain the extremity of galaxy formation theories. The most extreme are discovered through follow-up on candidates in large area surveys. Here we present 850 $\mu$m SCUBA-2 follow-up observations of 188 red DSFG candidates from the \textit{Herschel} Multi-tiered Extragalactic Survey (HerMES) Large Mode Survey, covering 274 deg$^2$. We detected 87 per cent with a signal-to-noise ratio $>$ 3 at 850~$\mu$m. We introduce a new method for incorporating the confusion noise in our spectral energy distribution fitting by sampling correlated flux density fluctuations from a confusion limited map. The new 850~$\mu$m data provide a better constraint on the photometric redshifts of the candidates, with photometric redshift errors decreasing from $\sigma_z/(1+z)\approx0.21$ to $0.15$. Comparison spectroscopic redshifts also found little bias ($\langle (z-z_{\rm spec})/(1+z_{\rm spec})\rangle = 0.08 $). The mean photometric redshift is found to be 3.6 with a dispersion of $0.4$ and we identify 21 DSFGs with a high probability of lying at $z > 4$. After simulating our selection effects we find number counts are consistent with phenomenological galaxy evolution models. There is a statistically significant excess of WISE-1 and SDSS sources near our red galaxies, giving a strong indication that lensing may explain some of the apparently extreme objects. Nevertheless, our sample should include examples of galaxies with the highest star formation rates in the Universe ($\gg10^3$ M$_\odot$yr$^{-1}$)., Comment: 22 pages, 14 figures and 3 tables

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., Comment: Resubmitted to MNRAS after request for minor revision

We present the Herschel-SPIRE photometric atlas for a complete flux limited sample of 43 local Ultraluminous Infrared Galaxies (ULIRGs), selected at 60${\mu}$m by IRAS, as part of the HERschel ULIRG Survey (HERUS). Photometry observations were obtained using the SPIRE instrument at 250, 350 and 500${\mu}$m. We describe these observations, present the results, and combine the new observations with data from IRAS to examine the far-IR spectral energy distributions (SEDs) of these sources. We fit the observed SEDs of HERUS objects with a simple parameterised modified black body model where temperature and emissivity $\beta$ are free parameters. We compare the fitted values to those of non-ULIRG local galaxies, and find, in agreement with earlier results, that HERUS ULIRGs have warmer dust (median temperature T = 37.9+/-4.7 K compared to 21.3+/-3.4 K) but a similar $\beta$ distribution (median $\beta$ = 1.7 compared to 1.8) to the Herschel reference sample (HRS, Cortese et al., 2014) galaxies. Dust masses are found to be in the range of 10^7.5 to 10^9 Msun significantly higher than that of Herschel Reference Sample (HRS) sources. We compare our results for local ULIRGs with higher redshift samples selected at 250 and 850${\mu}$m. These latter sources generally have cooler dust and/or redder 100-to-250 ${\mu}$m colours than our 60${\mu}$m-selected ULIRGs. We show that this difference may in part be the result of the sources being selected at different wavelengths rather than being a simple indication of rapid evolution in the properties of the population., Comment: Accepted for publication in MNRAS

We present a photometric catalog for Spitzer Space Telescope warm mission observations of the North Ecliptic Pole (NEP; centered at $\rm R.A.=18^h00^m00^s$, $\rm Decl.=66^d33^m38^s.552$). The observations are conducted with IRAC in 3.6 $\mu$m and 4.5 $\mu$m bands over an area of 7.04 deg$^2$ reaching 1$\sigma$ depths of 1.29 $\mu$Jy and 0.79 $\mu$Jy in the 3.6 $\mu$m and 4.5 $\mu$m bands respectively. The photometric catalog contains 380,858 sources with 3.6 $\mu$m and 4.5 $\mu$m band photometry over the full-depth NEP mosaic. Point source completeness simulations show that the catalog is 80% complete down to 19.7 AB. The accompanying catalog can be utilized in constraining the physical properties of extra-galactic objects, studying the AGN population, measuring the infrared colors of stellar objects, and studying the extra-galactic infrared background light., Comment: 10 pages, 11 figures and 3 tables. Accepted to the ApJS

The SPICA mid and far-infrared telescope will address fundamental issues in our understanding of star formation and ISM physics in galaxies. A particular hallmark of SPICA is the outstanding sensitivity enabled by the cold telescope, optimized detectors, and wide instantaneous bandwidth throughout the mid- and far-infrared. The spectroscopic, imaging and polarimetric observations that SPICA will be able to collect will help in clarifying the complex physical mechanisms which underlie the baryon cycle of galaxies. In particular: (i) The access to a large suite of atomic and ionic fine-structure lines for large samples of galaxies will shed light on the origin of the observed spread in star formation rates within and between galaxies. (ii) Observations of HD rotational lines (out to $\sim$10 Mpc) and fine structure lines such as [CII] 158 $\mu$m (out to $\sim$100 Mpc) will clarify the main reservoirs of interstellar matter in galaxies, including phases where CO does not emit. (iii) Far-infrared spectroscopy of dust and ice features will address uncertainties in the mass and composition of dust in galaxies, and the contributions of supernovae to the interstellar dust budget will be quantified by photometry and monitoring of supernova remnants in nearby galaxies. (iv) Observations of far-infrared cooling lines such as [OI] 63 $\mu$m from star-forming molecular clouds in our Galaxy will evaluate the importance of shocks to dissipate turbulent energy. The paper concludes with requirements for the telescope and instruments, and recommendations for the observing strategy., Comment: Accepted for publication in PASA; 19 pages, 8 figures, 1 table

We present images obtained with LABOCA on the APEX telescope of a sample of 22 galaxies selected via their red Herschel SPIRE 250-, 350- and $500\textrm{-}\mu\textrm{m}$ colors. We aim to see if these luminous, rare and distant galaxies are signposting dense regions in the early Universe. Our $870\textrm{-}\mu\textrm{m}$ survey covers an area of $\approx0.8\,\textrm{deg}^2$ down to an average r.m.s. of $3.9\,\textrm{mJy beam}^{-1}$, with our five deepest maps going $\approx2\times$ deeper still. We catalog 86 DSFGs around our 'signposts', detected above a significance of $3.5\sigma$. This implies a $100\pm30\%$ over-density of $S_{870}>8.5\,\textrm{mJy}$ DSFGs, excluding our signposts, when comparing our number counts to those in 'blank fields'. Thus, we are $99.93\%$ confident that our signposts are pinpointing over-dense regions in the Universe, and $\approx95\%$ confident that these regions are over-dense by a factor of at least $\ge1.5\times$. Using template SEDs and SPIRE/LABOCA photometry we derive a median photometric redshift of $z=3.2\pm0.2$ for our signposts, with an interquartile range of $z=2.8\textrm{-}3.6$. We constrain the DSFGs likely responsible for this over-density to within $|\Delta z|\le0.65$ of their respective signposts. These 'associated' DSFGs are radially distributed within $1.6\pm0.5\,\textrm{Mpc}$ of their signposts, have median SFRs of $\approx(1.0\pm0.2)\times10^3\,M_{\odot}\,\textrm{yr}^{-1}$ (for a Salpeter stellar IMF) and median gas reservoirs of $\sim1.7\times10^{11}\,M_{\odot}$. These candidate proto-clusters have average total SFRs of at least $\approx (2.3\pm0.5)\times10^3\,M_{\odot}\,\textrm{yr}^{-1}$ and space densities of $\sim9\times10^{-7}\,\textrm{Mpc}^{-3}$, consistent with the idea that their constituents may evolve to become massive ETGs in the centers of the rich galaxy clusters we see today., Comment: 24 pages, 13 figures, 5 tables, submitted for publication in ApJ

Our current knowledge of star formation and accretion luminosity at high-redshift (z>3-4), as well as the possible connections between them, relies mostly on observations in the rest-frame ultraviolet (UV), which are strongly affected by dust obscuration. Due to the lack of sensitivity of past and current infrared (IR) instrumentation, so far it has not been possible to get a glimpse into the early phases of the dust-obscured Universe. Among the next generation of IR observatories, SPICA, observing in the 12-350 micron range, will be the only facility that can enable us to make the required leap forward in understanding the obscured star-formation rate and black-hole accretion rate densities (SFRD and BHARD, respectively) with respect to what Spitzer and Herschel achieved in the mid- and far-IR at z<3. In particular, SPICA will have the unique ability to trace the evolution of the obscured SFRD and BHARD over cosmic time, from the peak of their activity back to the reionisation epoch (i.e., 3

IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features., Comment: Paper accepted for publication on PASA on 21st September 2017, as part of the SPICA Special Issue

Saturn's moon Enceladus has vents emerging from a sub-surface ocean, offering unique probes into the liquid environment. These vents drain into the larger neutral torus in orbit around Saturn. We present a methanol (CH3OH) detection observed with IRAM 30-m from 2008 along the line-of-sight through Saturn's E-ring. Additionally, we also present supporting observations from the Herschel public archive of water (ortho-H2O; 1669.9 GHz) from 2012 at a similar elongation and line-of-sight. The CH3OH 5(1,1)-4(1,1) transition was detected at 5.9 sigma confidence. The line has 0.43 km/s width and is offset by +8.1 km/s in the moon's reference frame. Radiative transfer models allow for gas cloud dimensions from 1750 km up to the telescope beam diameter ~73000 km. Taking into account the CH3OH lifetime against solar photodissociation and the redshifted line velocity, there are two possible explanations for the CH3OH emission: methanol is primarily a secondary product of chemical interactions within the neutral torus that (1) spreads outward throughout the E-ring or (2) originates from a compact, confined gas cloud lagging Enceladus by several km/s. We find either scenario to be consistent with significant redshifted H2O emission (4 sigma) measured from the Herschel public archive. The measured CH3OH:H2O abundance (> 0.5 per cent) significantly exceeds the observed abundance in the direct vicinity of the vents (~0.01 per cent), suggesting CH3OH is likely chemically processed within the gas cloud with methane (CH4) as its parent species., Comment: 16 pages, 4 figures, accepted for publication in the International Journal of Astrobiology (IJA)

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., Comment: Submitted to ApJ. Comments welcome

We report the identification of an extreme proto-cluster of galaxies in the early Universe whose core (nicknamed Distant Red Core, DRC) is formed by at least ten dusty star-forming galaxies (DSFGs), confirmed to lie at $z_{\rm spec} = 4.002$ via detection of [CI](1-0), $^{12}$CO(6-5), $^{12}$CO(4-3), $^{12}$CO(2-1) and ${\rm H_2O} (2_{11} - 2_{02})$ emission lines, detected using ALMA and ATCA. The spectroscopically-confirmed components of the proto-cluster are distributed over a ${\rm 260\, kpc \times 310\, kpc}$ region and have a collective obscured star-formation rate (SFR) of $\sim 6500 \, M_\odot \, {\rm yr}^{-1}$, considerably higher than has been seen before in any proto-cluster of galaxies or over-densities of DSFGs at $z \gtrsim 4$. Most of the star formation is taking place in luminous DSFGs since no Ly$\alpha$ emitters are detected in the proto-cluster core, apart from a Ly$\alpha$ blob located next to one of the DRC dusty components and extending over $60\,{\rm kpc}$. The total obscured SFR of the proto-cluster could rise to ${\rm SFR} \sim 14,400 \, M_\odot \, {\rm yr}^{-1}$ if all the members of an over-density of bright DSFGs discovered around DRC in a wide-field LABOCA 870-$\mu$m image are part of the same structure. The total halo mass of DRC could be as high as $\sim 4.4 \times 10^{13}\,M_\odot$ and could be the progenitor of a Coma-like cluster at $z = 0$. The relatively short gas-depletion times of the DRC components suggest either the presence of a mechanism able to trigger extreme star formation simultaneously in galaxies spread over a few hundred kpc or the presence of gas flows from the cosmic web able to sustain star formation over several hundred million years., Comment: Submitted to ApJ. Minor updates added, including a change of the source name. Comments welcome

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.

We present 1.3- and/or 3-mm continuum images and 3-mm spectral scans, obtained using NOEMA and ALMA, of 21 distant, dusty, star-forming galaxies (DSFGs). Our sample is a subset of the galaxies selected by Ivison et al. (2016) on the basis of their extremely red far-infrared (far-IR) colours and low {\it Herschel} flux densities; most are thus expected to be unlensed, extraordinarily luminous starbursts at $z \gtrsim 4$, modulo the considerable cross-section to gravitational lensing implied by their redshift. We observed 17 of these galaxies with NOEMA and four with ALMA, scanning through the 3-mm atmospheric window. We have obtained secure redshifts for seven galaxies via detection of multiple CO lines, one of them a lensed system at $z=6.027$ (two others are also found to be lensed); a single emission line was detected in another four galaxies, one of which has been shown elsewhere to lie at $z=4.002$. Where we find no spectroscopic redshifts, the galaxies are generally less luminous by 0.3-0.4 dex, which goes some way to explaining our failure to detect line emission. We show that this sample contains amongst the most luminous known star-forming galaxies. Due to their extreme star-formation activity, these galaxies will consume their molecular gas in $\lesssim 100$ Myr, despite their high molecular gas masses, and are therefore plausible progenitors of the massive, `red-and-dead' elliptical galaxies at $z \approx 3$., Comment: 15 pages, 10 figures; accepted for publication in MNRAS

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., Comment: Submitted to ApJ. Comments welcome

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission., Comment: 79 pages, 14 figures

Observational investigations of the abundance of massive precursors of local galaxy clusters ("proto-clusters") allow us to test the growth of density perturbations, to constrain cosmological parameters that control it, to test the theory of non-linear collapse and how the galaxy formation takes place in dense environments. The Planck collaboration has recently published a catalogue of >~ 2000 cold extra-galactic sub-millimeter sources, i.e. with colours indicative of z >~ 2, almost all of which appear to be over-densities of star-forming galaxies. They are thus considered as proto-cluster candidates. Their number densities (or their flux densities) are far in excess of expectations from the standard scenario for the evolution of large-scale structure. Simulations based on a physically motivated galaxy evolution model show that essentially all cold peaks brighter than S_{545GHz} = 500 mJy found in Planck maps after having removed the Galactic dust emission can be interpreted as positive Poisson fluctuations of the number of high-z dusty proto-clusters within the same Planck beam, rather then being individual clumps of physically bound galaxies. This conclusion does not change if an empirical fit to the luminosity function of dusty galaxies is used instead of the physical model. The simulations accurately reproduce the statistic of the Planck detections and yield distributions of sizes and ellipticities in qualitative agreement with observations. The redshift distribution of the brightest proto-clusters contributing to the cold peaks has a broad maximum at 1.5 <~ z <~ 3. Therefore follow-up of Planck proto-cluster candidates will provide key information on the high-z evolution of large scale structure., Comment: 10 pages, 8 figures, accepted for publication in MNRAS