Showing 1–50 of 384 results for author: Clements, D L

SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). V. Confusion-limited Submillimeter Galaxy Number Counts at 450 $μ$m and Data Release for the COSMOS Field

Authors: Zhen-Kai Gao, Chen-Fatt Lim, Wei-Hao Wang, Chian-Chou Chen, Ian Smail, Scott C. Chapman, Xian Zhong Zheng, Hyunjin Shim, Tadayuki Kodama, Yiping Ao, Siou-Yu Chang, David L. Clements, James S. Dunlop, Luis C. Ho, Yun-Hsin Hsu, Chorng-Yuan Hwang, Ho Seong Hwang, M. P. Koprowski, Douglas Scott, Stephen Serjeant, Yoshiki Toba, Sheona A. Urquhart

Abstract: We present confusion-limited SCUBA-2 450-$μ$m observations in the COSMOS-CANDELS region as part of the JCMT Large Program, SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). Our maps at 450 and 850 $μ$m cover an area of 450 arcmin$^2$. We achieved instrumental noise levels of $σ_{\mathrm{450}}=$ 0.59 mJy beam$^{-1}$ and $σ_{\mathrm{850}}=$ 0.09 mJy beam$^{-1}$ in the deepest area of each map. The co… ▽ More We present confusion-limited SCUBA-2 450-$μ$m observations in the COSMOS-CANDELS region as part of the JCMT Large Program, SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). Our maps at 450 and 850 $μ$m cover an area of 450 arcmin$^2$. We achieved instrumental noise levels of $σ_{\mathrm{450}}=$ 0.59 mJy beam$^{-1}$ and $σ_{\mathrm{850}}=$ 0.09 mJy beam$^{-1}$ in the deepest area of each map. The corresponding confusion noise levels are estimated to be 0.65 and 0.36 mJy beam$^{-1}$. Above the 4 (3.5) $σ$ threshold, we detected 360 (479) sources at 450 $μ$m and 237 (314) sources at 850 $μ$m. We derive the deepest blank-field number counts at 450 $μ$m, covering the flux-density range of 2 to 43 mJy. These are in agreement with other SCUBA-2 blank-field and lensing-cluster observations, but are lower than various model counts. We compare the counts with those in other fields and find that the field-to-field variance observed at 450 $μ$m at the $R=6^\prime$ scale is consistent with Poisson noise, so there is no evidence of strong 2-D clustering at this scale. Additionally, we derive the integrated surface brightness at 450 $μ$m down to 2.1 mJy to be $57.3^{+1.0}_{-6.2}$~Jy deg$^{-2}$, contributing to (41$\pm$4)\% of the 450-$μ$m extragalactic background light (EBL) measured by COBE and Planck. Our results suggest that the 450-$μ$m EBL may be fully resolved at $0.08^{+0.09}_{-0.08}$~mJy, which extremely deep lensing-cluster observations and next-generation submillimeter instruments with large aperture sizes may be able to achieve. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: 29 pages, 14 figures, accepted for publication in ApJ

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

A large population of strongly lensed faint submillimetre galaxies in future dark energy surveys inferred from JWST imaging

Authors: James Pearson, Stephen Serjeant, Wei-Hao Wang, Zhen-Kai Gao, Arif Babul, Scott Chapman, Chian-Chou Chen, David L. Clements, Christopher J. Conselice, James Dunlop, Lulu Fan, Luis C. Ho, Ho Seong Hwang, Maciej Koprowski, Michał Michałowski, Hyunjin Shim

Abstract: Bright galaxies at sub-millimetre wavelengths from Herschel are now well known to be predominantly strongly gravitationally lensed. The same models that successfully predicted this strongly lensed population also predict about one percent of faint $450μ$m-selected galaxies from deep James Clerk Maxwell Telescope (JCMT) surveys will also be strongly lensed. Follow-up ALMA campaigns have so far foun… ▽ More Bright galaxies at sub-millimetre wavelengths from Herschel are now well known to be predominantly strongly gravitationally lensed. The same models that successfully predicted this strongly lensed population also predict about one percent of faint $450μ$m-selected galaxies from deep James Clerk Maxwell Telescope (JCMT) surveys will also be strongly lensed. Follow-up ALMA campaigns have so far found one potential lens candidate, but without clear compelling evidence e.g. from lensing arcs. Here we report the discovery of a compelling gravitational lens system confirming the lensing population predictions, with a $z_{s} = 3.4 {\pm} 0.4$ submm source lensed by a $z_{spec} = 0.360$ foreground galaxy within the COSMOS field, identified through public JWST imaging of a $450μ$m source in the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) catalogue. These systems will typically be well within the detectable range of future wide-field surveys such as Euclid and Roman, and since sub-millimetre galaxies are predominantly very red at optical/near-infrared wavelengths, they will tend to appear in near-infrared channels only. Extrapolating to the Euclid-Wide survey, we predict tens of thousands of strongly lensed near-infrared galaxies. This will be transformative for the study of dusty star-forming galaxies at cosmic noon, but will be a contaminant population in searches for strongly lensed ultra-high-redshift galaxies in Euclid and Roman. △ Less

Submitted 9 January, 2024; v1 submitted 2 September, 2023; originally announced September 2023.

Comments: 9 pages, 5 figures, 4 tables; added figure, updated Discussion section, and minor revisions; accepted for publication in MNRAS

Observational evidence for cosmological coupling of black holes and its implications for an astrophysical source of dark energy

Authors: Duncan Farrah, Kevin S. Croker, Gregory Tarlé, Valerio Faraoni, Sara Petty, Jose Afonso, Nicolas Fernandez, Kurtis A. Nishimura, Chris Pearson, Lingyu Wang, Michael Zevin, David L Clements, Andreas Efstathiou, Evanthia Hatziminaoglou, Mark Lacy, Conor McPartland, Lura K Pitchford, Nobuyuki Sakai, Joel Weiner

Abstract: Observations have found black holes spanning ten orders of magnitude in mass across most of cosmic history. The Kerr black hole solution is however provisional as its behavior at infinity is incompatible with an expanding universe. Black hole models with realistic behavior at infinity predict that the gravitating mass of a black hole can increase with the expansion of the universe independently of… ▽ More Observations have found black holes spanning ten orders of magnitude in mass across most of cosmic history. The Kerr black hole solution is however provisional as its behavior at infinity is incompatible with an expanding universe. Black hole models with realistic behavior at infinity predict that the gravitating mass of a black hole can increase with the expansion of the universe independently of accretion or mergers, in a manner that depends on the black hole's interior solution. We test this prediction by considering the growth of supermassive black holes in elliptical galaxies over $0<z\lesssim2.5$. We find evidence for cosmologically coupled mass growth among these black holes, with zero cosmological coupling excluded at 99.98% confidence. The redshift dependence of the mass growth implies that, at $z\lesssim7$, black holes contribute an effectively constant cosmological energy density to Friedmann's equations. The continuity equation then requires that black holes contribute cosmologically as vacuum energy. We further show that black hole production from the cosmic star formation history gives the value of $Ω_Λ$ measured by Planck while being consistent with constraints from massive compact halo objects. We thus propose that stellar remnant black holes are the astrophysical origin of dark energy, explaining the onset of accelerating expansion at $z \sim 0.7$. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: 10 pages, 3 figures, published in ApJ Letters

Journal ref: ApJL 944 L31 (2023)

The SCUBA-2 Large eXtragalactic Survey: 850um map, catalogue and the bright-end number counts of the XMM-LSS field

Authors: T. K. Garratt, J. E. Geach, Y. Tamura, K. E. K. Coppin, M. Franco, Y. Ao, C. -C. Chen, C. Cheng, D. L. Clements, Y. S. Dai, H. Dannerbauer, T. R. Greve, B. Hatsukade, H. S. Hwang, L. Jiang, K. Kohno, M. P. Koprowski, M. J. Michalowski, M. Sawicki, D. Scott, H. Shim, T. T. Takeuchi, W. -H. Wang, Y. Q. Xue, C. Yang

Abstract: 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 t… ▽ More 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. △ Less

Submitted 25 January, 2023; originally announced January 2023.

Comments: Accepted for publication in MNRAS

Venus, Phosphine and the Possibility of Life

Authors: David L. Clements

Abstract: The search for life elsewhere in the universe is one of the central aims of science in the 21st century. While most of this work is aimed at planets orbiting other stars, the search for life in our own Solar System is an important part of this endeavour. Venus is often thought to have too harsh an environment for life, but it may have been a more hospitable place in the distant past. If life evolv… ▽ More The search for life elsewhere in the universe is one of the central aims of science in the 21st century. While most of this work is aimed at planets orbiting other stars, the search for life in our own Solar System is an important part of this endeavour. Venus is often thought to have too harsh an environment for life, but it may have been a more hospitable place in the distant past. If life evolved there in the past then the cloud decks of Venus are the only remaining niche where life as we know it might survive today. The discovery of the molecule phosphine, PH$_3$, in these clouds has reinvigorated research looking into the possibility of life in the clouds. In this review we examine the background to studies of the possibility of life on Venus, discuss the discovery of phosphine, review conflicting and confirming observations and analyses, and then look forward to future observations and space missions that will hopefully provide definitive answers as to the origin of phosphine on Venus and to the question of whether life might exist there. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: Invited review, accepted for publication in Contemporary Physics

A Preferential Growth Channel for Supermassive Black Holes in Elliptical Galaxies at z<2

Authors: Duncan Farrah, Sara Petty, Kevin Croker, Gregory Tarle, Michael Zevin, Evanthia Hatziminaoglou, Francesco Shankar, Lingyu Wang, David L Clements, Andreas Efstathiou, Mark Lacy, Kurtis A. Nishimura, Jose Afonso, Chris Pearson, Lura K Pitchford

Abstract: The assembly of stellar and supermassive black hole (SMBH) mass in elliptical galaxies since $z\sim1$ can help to diagnose the origins of locally-observed correlations between SMBH mass and stellar mass. We therefore construct three samples of elliptical galaxies, one at $z\sim0$ and two at $0.7\lesssim z \lesssim2.5$, and quantify their relative positions in the $M_{BH}-M_*$ plane. Using a Bayesi… ▽ More The assembly of stellar and supermassive black hole (SMBH) mass in elliptical galaxies since $z\sim1$ can help to diagnose the origins of locally-observed correlations between SMBH mass and stellar mass. We therefore construct three samples of elliptical galaxies, one at $z\sim0$ and two at $0.7\lesssim z \lesssim2.5$, and quantify their relative positions in the $M_{BH}-M_*$ plane. Using a Bayesian analysis framework, we find evidence for translational offsets in both stellar mass and SMBH mass between the local sample and both higher redshift samples. The offsets in stellar mass are small, and consistent with measurement bias, but the offsets in SMBH mass are much larger, reaching a factor of seven between $z\sim1$ and $z\sim0$. The magnitude of the SMBH offset may also depend on redshift, reaching a factor of $\sim20$ at $z\sim 2$. The result is robust against variation in the high and low redshift samples and changes in the analysis approach. The magnitude and redshift evolution of the offset are challenging to explain in terms of selection and measurement biases. We conclude that either there is a physical mechanism that preferentially grows SMBHs in elliptical galaxies at $z\lesssim 2$, or that selection and measurement biases are both underestimated, and depend on redshift. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: ApJ, accepted

Recovering Phosphine in Venus' Atmosphere from SOFIA Observations

Authors: Jane S. Greaves, Janusz J. Petkowski, Anita M. S. Richards, Clara Sousa-Silva, Sara Seager, David L. Clements

Abstract: Searches for phosphine in Venus' atmosphere have sparked a debate. Cordiner et al. 2022 analyse spectra from the Stratospheric Observatory For Infrared Astronomy (SOFIA) and infer <0.8 ppb of PH3. We noticed that some spectral artefacts arose from non-essential calibration-load signals. By-passing these signals allows simpler post-processing and a 5.7σ candidate detection, suggesting approx. 3 ppb… ▽ More Searches for phosphine in Venus' atmosphere have sparked a debate. Cordiner et al. 2022 analyse spectra from the Stratospheric Observatory For Infrared Astronomy (SOFIA) and infer <0.8 ppb of PH3. We noticed that some spectral artefacts arose from non-essential calibration-load signals. By-passing these signals allows simpler post-processing and a 5.7σ candidate detection, suggesting approx. 3 ppb of PH3 above the clouds. Compiling six phosphine results hints at an inverted abundance trend: decreasing above the clouds but rising again in the mesosphere from some unexplained source. However, no such extra source is needed if phosphine is undergoing destruction by sunlight (photolysis), to a similar degree as on Earth. Low phosphine values/limits are found where the viewed part of the super-rotating Venusian atmosphere had passed through sunlight, while high values are from views moving into sunlight. We suggest Venusian phosphine is indeed present, and so merits further work on models of its origins. △ Less

Submitted 2 March, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: v2 is an updated and revised version of the paper

Journal ref: Geophysical Research Letters, 50, e2023GL103539

A large-scale kinematic study of molecular gas in high-z cluster galaxies: Evidence for high levels of kinematic asymmetry

Authors: W. J. Cramer, A. G. Noble, K. Massingill, J. Cairns, D. L. Clements, M. C. Cooper, R. Demarco, J. Matharu, M. McDonald, A. Muzzin, J. Nantais, G. Rudnick, H. Übler, E. van Kampen, T. M. A. Webb, G. Wilson, H. K. C. Yee

Abstract: 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 kinemati… ▽ More 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. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: Submission pre-referee report

The star-formation rates of QSOs

Authors: M. Symeonidis, N. Maddox, M. J. Jarvis, M. J. Michalowski, P. Andreani, D. L. Clements, G. De Zotti, S. Duivenvoorden, J. Gonzalez-Nuevo, E. Ibar, R. J. Ivison, L. Leeuw, M. J. Page, R. Shirley, M. W. L. Smith, M. Vaccari

Abstract: We examine the far-IR properties of a sample of 5391 optically selected QSOs in the 0.5<z<2.65 redshift range down to log[nuLnu,2500 (erg/s)]>44.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 powe… ▽ More We examine the far-IR properties of a sample of 5391 optically selected QSOs in the 0.5<z<2.65 redshift range down to log[nuLnu,2500 (erg/s)]>44.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. △ Less

Submitted 23 May, 2022; originally announced May 2022.

Comments: 20 pages, 10 figures, accepted for publication in MNRAS

Stellar and black hole assembly in z<0.3 infrared-luminous mergers: intermittent starbursts vs. super-Eddington accretion

Authors: Duncan Farrah, Andreas Efstathiou, Jose Afonso, Jeronimo Bernard-Salas, Joe Cairns, David L Clements, Kevin Croker, Evanthia Hatziminaoglou, Maya Joyce, Mark Lacy, Vianney Lebouteiller, Alix Lieblich, Carol Lonsdale, Seb Oliver, Chris Pearson, Sara Petty, Lura K Pitchford, Dimitra Rigopoulou, Michael Rowan-Robinson, Jack Runburg, Henrik Spoon, Aprajita Verma, Lingyu Wang

Abstract: We study stellar and black hole mass assembly in a sample of 42 infrared-luminous galaxy mergers at z<0.3 by combining results from radiative transfer modelling with archival measures of molecular gas and black hole mass. The ratios of stellar mass, molecular gas mass, and black hole mass to each other are consistent with those of massive gas-rich galaxies at z<0.3. The advanced mergers may show i… ▽ More We study stellar and black hole mass assembly in a sample of 42 infrared-luminous galaxy mergers at z<0.3 by combining results from radiative transfer modelling with archival measures of molecular gas and black hole mass. The ratios of stellar mass, molecular gas mass, and black hole mass to each other are consistent with those of massive gas-rich galaxies at z<0.3. The advanced mergers may show increased black hole mass to stellar mass ratios, consistent with the transition from AGN to ellipticals and implying substantial black hole mass growth over the course of the merger. Star formation rates are enhanced relative to the local main sequence, by factors of ~100 in the starburst and ~1.8 in the host. The starburst star formation rates appear distinct to star formation in the main sequence at all redshifts up to at least z~5. Starbursts may prefer late-stage mergers, but are observed at any merger stage. We do not find evidence that the starbursts in these low-redshift systems substantially increase the total stellar mass, with a soft upper limit on the stellar mass increase from starburst activity of about a factor of two. In contrast, 12 objects show evidence for super-Eddington accretion, associated with late-stage mergers, suggesting that many AGN in infrared-luminous mergers go through a super-Eddington phase. The super-Eddington phase may increase black hole mass by up to an order of magnitude at an accretion efficiency of 42+/-33% over a period of 44+/-22Myr. Our results imply that super-Eddington accretion is an important black hole growth channel in infrared-luminous galaxies at all redshifts. △ Less

Submitted 29 April, 2022; originally announced May 2022.

Comments: MNRAS Accepted

A new look at local ultraluminous infrared galaxies: the atlas and radiative transfer models of their complex physics

Authors: A. Efstathiou, D. Farrah, J. Afonso, D. L. Clements, E. González-Alfonso, M. Lacy, S. Oliver, V. Papadopoulou Lesta, C. Pearson, D. Rigopoulou, M. Rowan-Robinson, H. W. W. Spoon, A. Verma, L. Wang

Abstract: 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 (A… ▽ More 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%. △ Less

Submitted 29 April, 2022; originally announced May 2022.

Comments: 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/

Journal ref: 2022MNRAS.512.5183E

The Nature of 500 Micron Risers II: Sub-mm Faint Dusty Star Forming Galaxies

Authors: J. Cairns, D. L. Clements, J. Greenslade, G. Petitpas, T. Cheng, Y. Ding, A. Parmar, I. Pérez-Fournon, D. Riechers

Abstract: 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$μ$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$μ$m risers, with three resolving into multiple systems. Fo… ▽ More 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$μ$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$μ$m risers, with three resolving into multiple systems. For these three multiple systems, the SMA 345GHz ($\approx 870μ$m) observations recover $123 \pm 73\%$, $60 \pm 15\%$ and $19 \pm 4\%$ respectively of the integrated 850$μ$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$μ$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$μ$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$μ$m maps to investigate the environments of these 500$μ$m risers. By constructing cumulative number counts and estimating photometric redshifts for surrounding SCUBA-2 detections, we find that one of our 500$μ$m risers could plausibly reside in a $z \geq 2$ protocluster. We infer that bright 500$μ$m risers with faint 850$μ$m flux densities are typically multiple systems at $z \geq 2$ that may reside in overdensities of bright sub-mm galaxies. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: 20 pages, 13 figures, 3 tables. Submitted to MNRAS

Consistent analysis of the AGN LF in X-ray and MIR in the XMM-LSS field

Authors: Jack Runburg, Duncan Farrah, Anna Sajina, Mark Lacy, Jenna Lidua, Evanthia Hatziminaoglou, W. N. Brandt, Chien-Ting J. Chen, Kristina Nyland, Raphael Shirley, D. L. Clements, Lura K. Pitchford

Abstract: 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 (… ▽ More 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-σ$ 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. △ Less

Submitted 30 November, 2021; originally announced December 2021.

Comments: 17 pages, 14 figures, 1 appendix

Characterization of Two 2mm-detected Optically-Obscured Dusty Star-Forming Galaxies

Authors: Sinclaire M. Manning, Caitlin M. Casey, Jorge A. Zavala, Georgios E. Magdis, Patrick M. Drew, Jaclyn B. Champagne, Manuel Aravena, Matthieu Béthermin, David L. Clements, Steven L. Finkelstein, Seiji Fujimoto, Christopher C. Hayward, Jacqueline A. Hodge, Olivier Ilbert, Jeyhan S. Kartaltepe, Kirsten K. Knudsen, Anton M. Koekemoer, Allison W. S. Man, David B. Sanders, Kartik Sheth, Justin S. Spilker, Johannes Staguhn, Margherita Talia, Ezequiel Treister, Min S. Yun

Abstract: The 2mm Mapping Obscuration to Reionization with ALMA (MORA) Survey was designed to detect high redshift ($z\gtrsim4$), massive, dusty star-forming galaxies (DSFGs). Here we present two, likely high redshift sources, identified in the survey whose physical characteristics are consistent with a class of optical/near-infrared (OIR) invisible DSFGs found elsewhere in the literature. We first perform… ▽ More The 2mm Mapping Obscuration to Reionization with ALMA (MORA) Survey was designed to detect high redshift ($z\gtrsim4$), massive, dusty star-forming galaxies (DSFGs). Here we present two, likely high redshift sources, identified in the survey whose physical characteristics are consistent with a class of optical/near-infrared (OIR) invisible DSFGs found elsewhere in the literature. We first perform a rigorous analysis of all available photometric data to fit spectral energy distributions and estimate redshifts before deriving physical properties based on our findings. Our results suggest the two galaxies, called MORA-5 and MORA-9, represent two extremes of the "OIR-dark" class of DSFGs. MORA-5 ($z_{\rm phot}=4.3^{+1.5}_{-1.3}$) is a significantly more active starburst with a star-formation rate of 830$^{+340}_{-190}$M$_\odot$yr$^{-1}$ compared to MORA-9 ($z_{\rm phot}=4.3^{+1.3}_{-1.0}$) whose star-formation rate is a modest 200$^{+250}_{-60}$M$_\odot$yr$^{-1}$. Based on the stellar masses (M$_{\star}\approx10^{10-11}$M$_\odot$), space density ($n\sim(5\pm2)\times10^{-6}$Mpc$^{-3}$, which incorporates two other spectroscopically confirmed OIR-dark DSFGs in the MORA sample at $z=4.6$ and $z=5.9$), and gas depletion timescales ($<1$Gyr) of these sources, we find evidence supporting the theory that OIR-dark DSFGs are the progenitors of recently discovered $3<z<4$ massive quiescent galaxies. △ Less

Submitted 3 November, 2021; originally announced November 2021.

Comments: 18 pages, 7 figures, accepted for publication in ApJ

Mapping Obscuration to Reionization with ALMA (MORA): 2mm Efficiently Selects the Highest-Redshift Obscured Galaxies

Authors: Caitlin M. Casey, Jorge A. Zavala, Sinclaire M. Manning, Manuel Aravena, Matthieu Béthermin, Karina I. Caputi, Jaclyn B. Champagne, David L. Clements, Patrick Drew, Steven L. Finkelstein, Seiji Fujimoto, Christopher C. Hayward, Anton M. Koekemoer, Vasily Kokorev, Claudia del P. Lagos, Arianna S. Long, Georgios E. Magdis, Allison W. S. Man, Ikki Mitsuhashi, Gergö Popping, Justin Spilker, Johannes Staguhn, Margherita Talia, Sune Toft, Ezequiel Treister , et al. (2 additional authors not shown)

Abstract: We present the characteristics of 2mm-selected sources from the largest Atacama Large Millimeter and submillimeter Array (ALMA) blank-field contiguous survey conducted to-date, the Mapping Obscuration to Reionization with ALMA (MORA) survey covering 184arcmin$^2$ at 2mm. Twelve of the thirteen detections above 5$σ$ are attributed to emission from galaxies, eleven of which are dominated by cold dus… ▽ More We present the characteristics of 2mm-selected sources from the largest Atacama Large Millimeter and submillimeter Array (ALMA) blank-field contiguous survey conducted to-date, the Mapping Obscuration to Reionization with ALMA (MORA) survey covering 184arcmin$^2$ at 2mm. Twelve of the thirteen detections above 5$σ$ are attributed to emission from galaxies, eleven of which are dominated by cold dust emission. These sources have a median redshift of $\langle z_{\rm 2mm}\rangle=3.6^{+0.4}_{-0.3}$ primarily based on optical/near-infrared (OIR) photometric redshifts with some spectroscopic redshifts, with 77$\pm$11% of sources at $z>3$ and 38$\pm$12% of sources at $z>4$. This implies that 2mm selection is an efficient method for identifying the highest redshift dusty star-forming galaxies (DSFGs). Lower redshift DSFGs ($z<3$) are far more numerous than those at $z>3$ yet likely to drop out at 2mm. MORA shows that DSFGs with star-formation rates in excess of 300M$_\odot$ yr$^{-1}$ and relative rarity of $\sim$10$^{-5}$ Mpc$^{-3}$ contribute $\sim$30% to the integrated star-formation rate density between $3<z<6$. The volume density of 2mm-selected DSFGs is consistent with predictions from some cosmological simulations and is similar to the volume density of their hypothesized descendants: massive, quiescent galaxies at $z>2$. Analysis of MORA sources' spectral energy distributions hint at steeper empirically-measured dust emissivity indices than typical literature studies, with $\langleβ\rangle=2.2^{+0.5}_{-0.4}$. The MORA survey represents an important step in taking census of obscured star-formation in the Universe's first few billion years, but larger area 2mm surveys are needed to more fully characterize this rare population and push to the detection of the Universe's first dusty galaxies. △ Less

Submitted 13 October, 2021; originally announced October 2021.

Comments: 42 pages, 14 figures, accepted for publication in ApJ

The HASHTAG project: The First Submillimeter Images of the Andromeda Galaxy from the Ground

Authors: Matthew W. L. Smith, Stephen A. Eales, Thomas G. Williams, Bumhyun Lee, Zongnan Li, Pauline Barmby, Martin Bureau, Scott Chapman, Brian S. Cho, Aeree Chung, Eun Jung Chung, Hui-Hsuan Chung, Christopher J. R. Clark, David L. Clements, Timothy A. Davis, Ilse De Looze, David J. Eden, Gayathri Athikkat-Eknath, George P. Ford, Yu Gao, Walter Gear, Haley L. Gomez, Richard de Grijs, Jinhua He, Luis C. Ho , et al. (24 additional authors not shown)

Abstract: Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing… ▽ More Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing the galaxy's large-scale structure. We have developed a technique for producing high-resolution submillimeter images of galaxies of large angular size by using the telescope on the ground to determine the small-scale structure (the large Fourier components) and a space telescope (Herschel or Planck) to determine the large-scale structure (the small Fourier components). Using this technique, we are carrying out the HARP and SCUBA-2 High Resolution Terahertz Andromeda Galaxy Survey (HASHTAG), an international Large Program on the James Clerk Maxwell Telescope, with one aim being to produce the first high-fidelity high-resolution submillimeter images of Andromeda. In this paper, we describe the survey, the method we have developed for combining the space-based and ground-based data, and present the first HASHTAG images of Andromeda at 450 and 850um. We also have created a method to predict the CO(J=3-2) line flux across M31, which contaminates the 850um band. We find that while normally the contamination is below our sensitivity limit, the contamination can be significant (up to 28%) in a few of the brightest regions of the 10 kpc ring. We therefore also provide images with the predicted line emission removed. △ Less

Submitted 30 September, 2021; originally announced October 2021.

Comments: 26 pages, 19 figures. Submitted to ApJS June 2021, Accepted September 2021

Low levels of sulphur dioxide contamination of Venusian phosphine spectra

Authors: Jane S. Greaves, Paul B. Rimmer, Anita M. S. Richards, Janusz J. Petkowski, William Bains, Sukrit Ranjan, Sara Seager, David L. Clements, Clara Sousa Silva, Helen J. Fraser

Abstract: New analysis is presented of the 1.1 mm wavelength absorption lines in Venus' atmosphere that suggested the presence of phosphine. We retrieve a sulphur dioxide observation from the JCMT archive that was simultaneous within a few days of the PH3 1-0 spectrum obtained in June 2017, and demonstrate via a radiative transfer calculation that contamination of PH3 by SO2 was ~10 per cent. We also presen… ▽ More New analysis is presented of the 1.1 mm wavelength absorption lines in Venus' atmosphere that suggested the presence of phosphine. We retrieve a sulphur dioxide observation from the JCMT archive that was simultaneous within a few days of the PH3 1-0 spectrum obtained in June 2017, and demonstrate via a radiative transfer calculation that contamination of PH3 by SO2 was ~10 per cent. We also present ALMA 2019 spectra of PH3 1-0 and an SO2 transition acquired simultaneously, and infer that SO2 line-contamination was ~2 percent (for the least-noisy half of the planetary disc). The contamination-subtracted ALMA and JCMT spectra (of 6-8 sigma confidence) are now consistent with similar absorption-depths at the two epochs. The two values span -1.9(+/-0.2) 10-4 of the continuum signal (which was re-estimated for ALMA), albeit for differing planetary areas. This suggests that the abundance attributed to phosphine in Venus' atmosphere was broadly similar in 2017 and 2019. △ Less

Submitted 21 June, 2022; v1 submitted 18 August, 2021; originally announced August 2021.

Comments: v3 is an proofread version of the paper that will appear in print in MNRAS

Journal ref: MNRAS 2022

Recovery of Spectra of Phosphine in Venus' Clouds

Authors: Jane S. Greaves, Anita M. S. Richards, William Bains, Paul B. Rimmer, David L. Clements, Sara Seager, Janusz J. Petkowski, Clara Sousa-Silva, Sukrit Ranjan, Helen J. Fraser

Abstract: We recover PH3 in the atmosphere of Venus in data taken with ALMA, using three different calibration methods. The whole-planet signal is recovered with 5.4σ confidence using Venus bandpass self-calibration, and two simpler approaches are shown to yield example 4.5-4.8σ detections of the equatorial belt. Non-recovery by Villanueva et al. is attributable to (a) including areas of the planet with hig… ▽ More We recover PH3 in the atmosphere of Venus in data taken with ALMA, using three different calibration methods. The whole-planet signal is recovered with 5.4σ confidence using Venus bandpass self-calibration, and two simpler approaches are shown to yield example 4.5-4.8σ detections of the equatorial belt. Non-recovery by Villanueva et al. is attributable to (a) including areas of the planet with high spectral-artefacts and (b) retaining all antenna baselines which raises the noise by a factor ~2.5. We release a data-processing script that enables our whole-planet result to be reproduced. The JCMT detection of PH3 remains robust, with the alternative SO2 attribution proposed by Villanueva et al. appearing inconsistent both in line-velocity and with millimetre-wavelength SO2 monitoring. SO2 contamination of the ALMA PH3-line is minimal. Net abundances for PH3, in the gas column above ~55 km, are up to ~20 ppb planet-wide with JCMT, and ~7 ppb with ALMA (but with signal-loss possible on scales approaching planetary size). Derived abundances will differ if PH3 occupies restricted altitudes - molecules in the clouds will contribute significantly less absorption at line-centre than equivalent numbers of mesospheric molecules - but in the latter zone, PH3 lifetime is expected to be short. Given we recover phosphine, we suggest possible solutions (requiring substantial further testing): a small collisional broadening coefficient could give narrow lines from lower altitude, or a high eddy diffusion coefficient could allow molecules to survive longer at higher altitudes. Alternatively, PH3 could be actively produced by an unknown mechanism in the mesosphere, but this would need to be in addition to cloud-level PH3 detected retrospectively by Pioneer-Venus. △ Less

Submitted 19 April, 2021; originally announced April 2021.

Comments: Revised and updated re-analysis of the phosphine spectra in Venus' clouds; Submitted to Nature Astronomy "Matters Arising" on Apr/19/2021

Journal ref: Nature Astronomy 5.7 (2021): 636-639

The Evolution of the IR Luminosity Function and Dust-obscured Star Formation in the Last 13 Billion Years

Authors: J. A. Zavala, C. M. Casey, S. M. Manning, M. Aravena, M. Bethermin, K. I. Caputi, D. L. Clements, E. da Cunha, P. Drew, S. L. Finkelstein, S. Fujimoto, C. Hayward, J. Hodge, J. S. Kartaltepe, K. Knudsen, A. M. Koekemoer, A. S. Long, G. E. Magdis, A. W. S. Man, G. Popping, D. Sanders, N. Scoville, K. Sheth, J. Staguhn, S. Toft , et al. (3 additional authors not shown)

Abstract: 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 comb… ▽ More 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 $α_{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. △ Less

Submitted 12 January, 2021; originally announced January 2021.

Comments: Accepted for publication in The Astrophysical Journal

On the Robustness of Phosphine Signatures in Venus' Clouds

Authors: Jane S. Greaves, William Bains, Janusz J. Petkowski, Sara Seager, Clara Sousa-Silva, Sukrit Ranjan, David L. Clements, Paul B. Rimmer, Helen J. Fraser, Steve Mairs, Malcolm J. Currie

Abstract: We published spectra of phosphine molecules in Venus' clouds, following open-science principles in releasing data and scripts (with community input leading to ALMA re-processing, now benefiting multiple projects). Some misconceptions about de-trending of spectral baselines have also emerged, which we address here. Using the JCMT PH3-discovery data, we show that mathematically-correct polynomial fi… ▽ More We published spectra of phosphine molecules in Venus' clouds, following open-science principles in releasing data and scripts (with community input leading to ALMA re-processing, now benefiting multiple projects). Some misconceptions about de-trending of spectral baselines have also emerged, which we address here. Using the JCMT PH3-discovery data, we show that mathematically-correct polynomial fitting of periodic ripples does not lead to "fake lines" (probability < ~1%). We then show that the ripples can be characterised in a non-subjective manner via Fourier transforms. A 20 ppb PH3 feature is ~5σ compared to the JCMT baseline-uncertainty, and is distinctive as a narrow perturber of the periodic ripple pattern. The structure of the FT-derived baseline also shows that polynomial fitting, if unguided, can amplify artefacts and so artificially reduce significance of real lines. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Comments: Response to: Snellen I. et al. Astron. Astroph., in press, arXiv:2010:09761 (2020); Villanueva G. et al. Nat. Ast. Matters Arising; arXiv:2010.14305 (2020); Thompson M. Mon. Not. Roy. Astron. Soc., in press, arXiv:2010.15118 (2020); Submitted to Nature Astronomy "Matters Arising" on Dec/10/2020

Journal ref: Nature Astronomy 5.7 (2021): 726-728

A Spitzer survey of Deep Drilling Fields to be targeted by the Vera C. Rubin Observatory Legacy Survey of Space and Time

Authors: M. Lacy, J. A. Surace, D. Farrah, K. Nyland, J. Afonso, W. N. Brandt, D. L. Clements, C. D. P. Lagos, C. Maraston, J. Pforr, A. Sajina, M. Sako, M. Vaccari, G. Wilson, D. R. Ballantyne, W. A. Barkhouse, R. Brunner, R. Cane, T. E. Clarke, M. Cooper, A. Cooray, G. Covone, C. D'Andrea, A. E. Evrard, H. C. Ferguson , et al. (38 additional authors not shown)

Abstract: 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 $μ$m and… ▽ More 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 $μ$m and 4.5 $μ$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σ$ point-source depth of 2 $μ$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. △ Less

Submitted 30 November, 2020; originally announced November 2020.

Comments: 20 pages, 9 figures; MNRAS in press

Re-analysis of Phosphine in Venus' Clouds

Authors: Jane S. Greaves, Anita M. S. Richards, William Bains, Paul B. Rimmer, David L. Clements, Sara Seager, Janusz J. Petkowski, Clara Sousa-Silva, Sukrit Ranjan, Helen J. Fraser

Abstract: We first respond to two points raised by Villanueva et al. We show the JCMT discovery spectrum of PH3 can not be re-attributed to SO2, as the line width is larger than observed for SO2 features, and the required abundance would be an extreme outlier. The JCMT spectrum is also consistent with our simple model, constant PH3-abundance with altitude, with no discrepancy in line profile (within data li… ▽ More We first respond to two points raised by Villanueva et al. We show the JCMT discovery spectrum of PH3 can not be re-attributed to SO2, as the line width is larger than observed for SO2 features, and the required abundance would be an extreme outlier. The JCMT spectrum is also consistent with our simple model, constant PH3-abundance with altitude, with no discrepancy in line profile (within data limits); reconciliation with a full photochemical model is the subject of future work. Section 2 presents initial results from re-processed ALMA data. Villanueva et al. noted an issue with bandpass calibration. They have worked on a partially re-processed subset of the ALMA data, so we note where their conclusions, and those of Greaves et al., are now superseded. To summarise: we recover PH3 in Venus' atmosphere with ALMA (~5σ confidence). Localised abundance appears to peak at ~5-10 parts-per-billion (ppb), with suggestions of spatial variation. Advanced data-products suggest a planet-averaged PH3 abundance ~1-4 ppb, lower than from the earlier ALMA processing (which indicated 7+ ppb). The ALMA data are reconcilable with the JCMT detection (~20 ppb) if there is order-of-magnitude temporal variation; more advanced processing of the JCMT data is underway to check methods. Independent PH3 measurements suggest possible altitude dependence (under ~5 ppb at 60+ km, up to ~100 ppb at 50+ km; see Section 2: Conclusions.). Given that both ALMA and JCMT were working at the limit of observatory capabilities, new spectra should be obtained. The ALMA data in-hand are no longer limited by calibration, but spectral ripples still exist, probably due to size and brightness of Venus in relation to the primary beam. Further, spatial ripples are present, potentially reducing significance of real narrow spectral features. △ Less

Submitted 10 December, 2020; v1 submitted 16 November, 2020; originally announced November 2020.

Comments: Response to the Matters Arising: "No phosphine in the atmosphere of Venus" by G. Villanueva et al., preprint at arXiv:2010.14305; v2 contains updated values and additional technical details

Journal ref: Nature Astronomy 5.7 (2021): 636-639

Rise of the Titans: Gas Excitation and Feedback in a Binary Hyper-Luminous Dusty Starburst Galaxy at z~6

Authors: Dominik A. Riechers, Hooshang Nayyeri, Denis Burgarella, Bjorn H. C. Emonts, David L. Clements, Asantha Cooray, Rob J. Ivison, Seb Oliver, Ismael Perez-Fournon, Dimitra Rigopoulou, Douglas Scott

Abstract: We report new observations toward the hyper-luminous dusty starbursting major merger ADFS-27 (z=5.655), using ATCA and ALMA. We detect CO 2-1, 8-7, 9-8, 10-9 and H2O(321-221) emission, and a P-Cygni-shaped OH+(11-01) absorption/emission feature. We also tentatively detect H2O(321-312) and OH+(12-01) emission and CH+(1-0) absorption. We find a total cold molecular mass of M_gas = (2.1+/-0.2) x 10^1… ▽ More We report new observations toward the hyper-luminous dusty starbursting major merger ADFS-27 (z=5.655), using ATCA and ALMA. We detect CO 2-1, 8-7, 9-8, 10-9 and H2O(321-221) emission, and a P-Cygni-shaped OH+(11-01) absorption/emission feature. We also tentatively detect H2O(321-312) and OH+(12-01) emission and CH+(1-0) absorption. We find a total cold molecular mass of M_gas = (2.1+/-0.2) x 10^11 (alpha_CO/1.0) Msun. We also find that the excitation of the star-forming gas is overall moderate for a z>5 dusty starburst, which is consistent with its moderate dust temperature. A high density, high kinetic temperature gas component embedded in the gas reservoir is required to fully explain the CO line ladder. This component is likely associated with the "maximum starburst" nuclei in the two merging galaxies, which are separated by only (140+/-13) km/s along the line of sight and 9.0 kpc in projection. The kinematic structure of both components is consistent with galaxy disks, but this interpretation remains limited by the spatial resolution of the current data. The OH+ features are only detected towards the northern component, which is also the one that is more enshrouded in dust and thus remains undetected up to 1.6 um even in our sensitive new HST/WFC3 imaging. The absorption component of the OH+ line is blueshifted and peaks near the CO and continuum emission peak while the emission is redshifted and peaks offset by 1.7 kpc from the CO and continuum emission peak, suggesting that the gas is associated with a massive molecular outflow from the intensely star-forming nucleus that supplies 125 Msun/yr of enriched gas to its halo. △ Less

Submitted 30 November, 2020; v1 submitted 28 October, 2020; originally announced October 2020.

Comments: 20 pages, 11 figures, 6 tables; to appear in the Astrophysical Journal (accepted November 29, 2020)

Phosphine Gas in the Cloud Decks of Venus

Authors: Jane S. Greaves, Anita M. S. Richards, William Bains, Paul B. Rimmer, Hideo Sagawa, David L. Clements, Sara Seager, Janusz J. Petkowski, Clara Sousa-Silva, Sukrit Ranjan, Emily Drabek-Maunder, Helen J. Fraser, Annabel Cartwright, Ingo Mueller-Wodarg, Zhuchang Zhan, Per Friberg, Iain Coulson, E'lisa Lee, Jim Hoge

Abstract: Measurements of trace-gases in planetary atmospheres help us explore chemical conditions different to those on Earth. Our nearest neighbor, Venus, has cloud decks that are temperate but hyper-acidic. We report the apparent presence of phosphine (PH3) gas in Venusian atmosphere, where any phosphorus should be in oxidized forms. Single-line millimeter-waveband spectral detections (quality up to ~15… ▽ More Measurements of trace-gases in planetary atmospheres help us explore chemical conditions different to those on Earth. Our nearest neighbor, Venus, has cloud decks that are temperate but hyper-acidic. We report the apparent presence of phosphine (PH3) gas in Venusian atmosphere, where any phosphorus should be in oxidized forms. Single-line millimeter-waveband spectral detections (quality up to ~15 sigma) from the JCMT and ALMA telescopes have no other plausible identification. Atmospheric PH3 at ~20 parts-per-billion abundance is inferred. The presence of phosphine is unexplained after exhaustive study of steady-state chemistry and photochemical pathways, with no currently-known abiotic production routes in Venusian atmosphere, clouds, surface and subsurface, or from lightning, volcanic or meteoritic delivery. Phosphine could originate from unknown photochemistry or geochemistry, or, by analogy with biological production of phosphine on Earth, from the presence of life. Other PH3 spectral features should be sought, while in-situ cloud and surface sampling could examine sources of this gas. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Journal ref: Greaves, J. S., Richards, A. M., Bains, W., Rimmer, P. B., Sagawa, H., Clements, D. L., ... & Hoge, J. (2020). Phosphine gas in the cloud decks of Venus. Nature Astronomy, 1-10

The Nature of 500 micron Risers I: SMA Observations

Authors: J. Greenslade, D. L. Clements, G. Petitpas, V. Asboth, A. Conley, I. Pérez-Fournon, D. Riechers

Abstract: 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$μ$m, a population termed 500-risers. We detect counterparts to 24 of these sources, with four having two counterparts. We conclude that th… ▽ More 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$μ$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. △ Less

Submitted 12 June, 2020; originally announced June 2020.

Comments: Accepted for publication in MNRAS

IRAM 30m-EMIR Redshift Search of z = 3-4 Lensed Dusty Starbursts selected from the HerBS sample

Authors: T. J. L. C. Bakx, H. Dannerbauer, D. Frayer, S. A. Eales, I. Pérez-Fournon, Z. -Y. Cai, D. L. Clements, G. De Zotti, J. González-Nuevo, R. J. Ivison, A. Lapi, M. J. Michałowski, M. Negrello, S. Serjeant, M. W. L. Smith, P. Temi, S. Urquhart, P. van der Werf

Abstract: 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 $μ$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… ▽ More 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 $μ$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. △ Less

Submitted 14 July, 2020; v1 submitted 10 June, 2020; originally announced June 2020.

Comments: 20 Pages, 8 Figures + 7 Appendix Figures + 1 animated figure. Resubmitted in order to reflect the version as printed in MNRAS

JINGLE -- IV. Dust, HI gas and metal scaling laws in the local Universe

Authors: I. De Looze, I. Lamperti, A. Saintonge, M. Relano, M. W. L. Smith, C. J. R. Clark, C. D. Wilson, M. Decleir, A. P. Jones, R. C. Kennicutt, G. Accurso, E. Brinks, M. Bureau, P. Cigan, D. L. Clements, P. De Vis, L Fanciullo, Y. Gao, W. K. Gear, L. C. Ho, H. S. Hwang, M. J. Michalowski, J. C. Lee, C. Li, L. Lin , et al. (7 additional authors not shown)

Abstract: 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 a… ▽ More 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 ($ε$=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. △ Less

Submitted 2 June, 2020; originally announced June 2020.

Comments: 56 pages, 30 figures, accepted for publication in MNRAS

SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) IV: Spatial clustering and halo masses of 450-$μ$m-selected sub-millimeter galaxies

Authors: Chen-Fatt Lim, Chian-Chou Chen, Ian Smail, Wei-Hao Wang, Wei-Leong Tee, Yen-Ting Lin, Douglas Scott, Yoshiki Toba, Yu-Yen Chang, YiPing Ao, Arif Babul, Andy Bunker, Scott C. Chapman, David L Clements, Christopher J. Conselice, Yu Gao, Thomas R. Greve, Luis C. Ho, Sungwook E. Hong, Ho Seong Hwang, Maciej Koprowski, Michał J. Michałowski, Hyunjin Shim, Xinwen Shu, James M. Simpson

Abstract: We analyze an extremely deep 450-$μ$m image ($1σ=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-$μ$m that have $K$-band counterparts in the COSMOS… ▽ More We analyze an extremely deep 450-$μ$m image ($1σ=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-$μ$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-$μ$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-$μ$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-$μ$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-$μ$m SMG candidate sample, we measure the two-point autocorrelation functions from $z=3$ down to $z=0.5$. We find that the 450-$μ$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. △ Less

Submitted 2 June, 2020; v1 submitted 28 April, 2020; originally announced April 2020.

Comments: ApJ accepted

SCUBA-2 overdensities associated with candidate protoclusters selected from Planck data

Authors: T. Cheng, D. L. Clements, J. Greenslade, J. Cairns, P. Andreani, M. Bremer, L. Conversi, A. Cooray, H. Dannerbauer, G. De Zotti, S. Eales, J. González-Nuevo, E. Ibar, L. Leeuw, J. Ma, M. J. Michałowski, H. Nayyeri, D. A. Riechers, D. Scott, P. Temi, M. Vaccari, I. Valtchanov, E. van Kampen, L. Wang

Abstract: We measure the 850-$μ$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-$μ$m sources in order to select the fields that are most likely to be genuine protoclusters. Of the 46 candidate protocluste… ▽ More We measure the 850-$μ$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-$μ$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-$μ$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-$μ$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. △ Less

Submitted 18 April, 2020; originally announced April 2020.

Comments: 8 pages, 1 figure, 1 table, accepted for publication in MNRAS

Submillimeter Galaxy studies in the next decade: EAO Submillimetre Futures White Paper Series, 2019

Authors: Ran Wang, Wei-Hao Wang, David L. Clements, Haojing Yan, Yiping Ao

Abstract: Over the last two decades, the Submillimetre Common-User Bolometer Array (SCUBA) and SCUBA-2 on the James Clerk Maxwell Telescope (JCMT) achieved gread success in discovering the population of dusty starburst galaxies in the early universe. The SCUBA-2 surveys at 450 micron and 850 micron set important constraints on the obscured star formation over cosmic time, and in combination of deep optical… ▽ More Over the last two decades, the Submillimetre Common-User Bolometer Array (SCUBA) and SCUBA-2 on the James Clerk Maxwell Telescope (JCMT) achieved gread success in discovering the population of dusty starburst galaxies in the early universe. The SCUBA-2 surveys at 450 micron and 850 micron set important constraints on the obscured star formation over cosmic time, and in combination of deep optical and near-IR data, allows the study of protoclusters and structure formation. However, the current submillimeter (submm) surveys by JCMT are still limited by area of sky coverage (confusion limit mapping of only a few deg^2), which prevent a systematic study of large samples of the obscured galaxy population. In this white paper, we review the studies of the submm galaxies with current submillimeter/millimeter (submm/mm) observations, and discuss the important science with the new submm instruments in the next decade. In particular, with a 10 times faster mapping speed of the new camera, we will expect deep 850 micron surveys over 10 to 100 times larger sky area to i) largely increase the sample size of submm detections toward the highest redshift, ii) improve our knowledge of galaxy and structure formation in the early universe. △ Less

Submitted 21 January, 2020; originally announced January 2020.

Comments: 10 pages 3 figures. EAO Submillimetre Futures White Paper Series, 2019

Planck far-infrared detection of Hyper Suprime-Cam protoclusters at $\bf z\sim4$: hidden AGN and star formation activity

Authors: Mariko Kubo, Jun Toshikawa, Nobunari Kashikawa, Yi-Kuan Chiang, Roderik Overzier, Hisakazu Uchiyama, David L. Clements, David M. Alexander, Yuichi Matsuda, Tadayuki Kodama, Yoshiaki Ono, Tomotsugu Goto, Tai-An Cheng, Kei Ito

Abstract: We perform a stacking analysis of {\it Planck}, {\it AKARI}, Infrared Astronomical Satellite ($IRAS$), Wide-field Infrared Survey Eplorer ($WISE$), and {\it Herschel} images of the largest number of (candidate) protoclusters at $z\sim3.8$ selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Stacking the images of the $179$ candidate protoclusters, the combined infrared (IR) emis… ▽ More We perform a stacking analysis of {\it Planck}, {\it AKARI}, Infrared Astronomical Satellite ($IRAS$), Wide-field Infrared Survey Eplorer ($WISE$), and {\it Herschel} images of the largest number of (candidate) protoclusters at $z\sim3.8$ selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Stacking the images of the $179$ candidate protoclusters, the combined infrared (IR) emission of the protocluster galaxies in the observed $12-850~μ$m wavelength range is successfully detected with $>5σ$ significance (at $Planck$). This is the first time that the average IR spectral energy distribution (SED) of a protocluster has been constrained at $z\sim4$. The observed IR SEDs of the protoclusters exhibit significant excess emission in the mid-IR compared to that expected from typical star-forming galaxies (SFGs). They are reproduced well using SED models of intense starburst galaxies with warm/hot dust heated by young stars, or by a population of active galactic nuclei (AGN)/SFG composites. For the pure star-forming model, a total IR (from 8 to 1000 $μ$m) luminosity of $19.3_{-4.2}^{+0.6}\times10^{13}~L_{\odot}$ and a star formation rate (SFR) of $16.3_{-7.8}^{+1.0}\times10^3~M_{\odot}$ yr$^{-1}$ are found whereas for the AGN/SFG composite model, $5.1_{-2.5}^{+2.5}\times10^{13}~L_{\odot}$ and $2.1^{+6.3}_{-1.7}\times10^3~M_{\odot}$ yr$^{-1}$ are found. Uncertainty remaining in the total SFRs; however, the IR luminosities of the most massive protoclusters are likely to continue increasing up to $z\sim4$. Meanwhile, no significant IR flux excess is observed around optically selected QSOs at similar redshifts, which confirms previous results. Our results suggest that the $z\sim4$ protoclusters trace dense, intensely star-forming environments that may also host obscured AGNs missed by the selection in the optical. △ Less

Submitted 21 November, 2019; originally announced November 2019.

Comments: 26 pages, 12 figures, Accepted for publication in ApJ

Have we seen all the galaxies that comprise the cosmic infrared background at 250\,$μ$m $\le λ\le$ 500\,$μ$m?

Authors: S. Duivenvoorden, S. Oliver, M. Bethermin, D. L. Clements, G. De Zotti, A. Efstathiou, D. Farrah, P. D. Hurley, R. J. Ivison, G. Lagache, D. Scott, R. Shirley, L. Wang, M. Zemcov

Abstract: 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 t… ▽ More 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\,$μ$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\,$μ$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 $λ\le 500\,μ$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\,$μ$m $\leq λ\leq$ 500\,$μ$m. △ Less

Submitted 4 November, 2019; originally announced November 2019.

Comments: 16 pages, 14 figures

Physical Characterization of an Unlensed Dusty Star-Forming Galaxy at $z=5.85$

Authors: Caitlin M. Casey, Jorge A. Zavala, Manuel Aravena, Matthieu Bethermin, Karina I. Caputi, Jaclyn B. Champagne, David L. Clements, Elisabete Da Cunha, Patrick Drew, Steven L. Finkelstein, Christopher C. Hayward, Jeyhan S. Kartaltepe, Kirsten Knudsen, Anton M. Koekemoer, Georgios E. Magdis, Allison Man, Sinclaire M. Manning, Nick Z. Scoville, Kartik Sheth, Justin Spilker, Johannes Staguhn, Margherita Talia, Yoshiaki Taniguchi, Sune Toft, Ezequiel Treister , et al. (1 additional authors not shown)

Abstract: We present a physical characterization of MMJ100026.36+021527.9 (a.k.a. ``MAMBO-9''), a dusty star-forming galaxy (DSFG) at $z=5.850\pm0.001$. This is the highest redshift unlensed DSFG (and fourth most distant overall) found to-date, and is the first source identified in a new 2mm blank-field map in the COSMOS field. Though identified in prior samples of DSFGs at 850$μ$m-1.2mm with unknown redshi… ▽ More We present a physical characterization of MMJ100026.36+021527.9 (a.k.a. ``MAMBO-9''), a dusty star-forming galaxy (DSFG) at $z=5.850\pm0.001$. This is the highest redshift unlensed DSFG (and fourth most distant overall) found to-date, and is the first source identified in a new 2mm blank-field map in the COSMOS field. Though identified in prior samples of DSFGs at 850$μ$m-1.2mm with unknown redshift, the detection at 2mm prompted further follow-up as it indicated a much higher probability that the source was likely to sit at $z>4$. Deep observations from the Atacama Large Millimeter and submillimeter Array (ALMA) presented here confirm the redshift through the secure detection of $^{12}$CO($J\!=$6$\rightarrow$5) and p-H$_{2}$O(2$_{1,1}\!\rightarrow$2$_{0,2}$). MAMBO-9 is comprised of a pair of galaxies separated by 6kpc with corresponding star-formation rates of 590M$_\odot$yr$^{-1}$ and 220M$_\odot$yr$^{-1}$ total molecular hydrogen gas mass of (1.7$\pm$0.4)$\times10^{11}$M$_\odot$, dust mass of (1.3$\pm$0.3)$\times10^{9}$M$_\odot$ and stellar mass of (3.2$^{+1.0}_{-1.5}$)$\times10^{9}$M$_\odot$. The total halo mass, (3.3$\pm$0.8)$\times10^{12}$M$_\odot$, is predicted to exceed $>10^{15}$M$_\odot$ by $z=0$. The system is undergoing a merger-driven starburst which will increase the stellar mass of the system tenfold in $τ_{\rm depl}=40-80$Myr, converting its large molecular gas reservoir (gas fraction of 96$^{+1}_{-2}$%) into stars. MAMBO-9 evaded firm spectroscopic identification for a decade, following a pattern that has emerged for some of the highest redshift DSFGs found. And yet, the systematic identification of unlensed DSFGs like MAMBO-9 is key to measuring the global contribution of obscured star-formation to the star-formation rate density at $z>4$, the formation of the first massive galaxies, and the formation of interstellar dust at early times ($<$1Gyr). △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: 24 pages, 8 figures, accepted for publication in ApJ

A SCUBA-2 Selected Herschel-SPIRE Dropout and the Nature of this Population

Authors: J. Greenslade, E. Aguilar, D. L. Clements, H. Dannerbauer, T. Cheng, G. Petitpas, C. Yang, H. Messias, I. Oteo, D. Farrah, M. J. Michalowski, I. Perez Fournon, I. Aretxaga, M. S. Yun, S. Eales, L. Dunne, A. Cooray, P. Andreani, D. H. Hughes, M. Velazquez, D. Sanchez-Arguelles, N. Ponthieu

Abstract: 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 $μm$ source (12.3 $\pm$ 2.5 mJy) with no counterparts at shorter wavelengths (a SPIRE dropout). Inte… ▽ More 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 $μ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σ$ 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. △ Less

Submitted 8 October, 2019; originally announced October 2019.

Comments: Accepoted for publication in MNRAS

SCUBA-2 observations of candidate starbursting protoclusters selected by Planck and Herschel-SPIRE

Authors: T. Cheng, D. L. Clements, J. Greenslade, J. Cairns, P. Andreani, M. Bremer, L. Conversi, A. Cooray, H. Dannerbauer, G. De Zotti, S. Eales, J. González-Nuevo, E. Ibar, L. Leeuw, J. Ma, M. J. Michałowski, H. Nayyeri, D. A. Riechers, D. Scott, P. Temi, M. Vaccari, I. Valtchanov, E. van Kampen, L. Wang

Abstract: We present SCUBA-2 850-$μ$m observations of 13 candidate starbursting protoclusters selected using Planck and Herschel data. The cumulative number counts of the 850-$μ$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-… ▽ More We present SCUBA-2 850-$μ$m observations of 13 candidate starbursting protoclusters selected using Planck and Herschel data. The cumulative number counts of the 850-$μ$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-$μ$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-$μ$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-$μ$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-$μ$m sources. △ Less

Submitted 27 September, 2019; v1 submitted 19 September, 2019; originally announced September 2019.

Comments: 30 pages, 8 figures, 5 tables, accepted for publication in MNRAS

JINGLE V: Dust properties of nearby galaxies derived from hierarchical Bayesian SED fitting

Authors: Isabella Lamperti, Amélie Saintonge, Ilse De Looze, Gioacchino Accurso, Christopher J. R. Clark, Matthew W. L. Smith, Christine D. Wilson, Elias Brinks, Toby Brown, Martin Bureau, David L. Clements, Stephen Eales, David H. W. Glass, Ho Seong Hwang, Jong Chul Lee, Lihwai Lin, Michal J. Michalowski, Mark Sargent, Thomas G. Williams, Ting Xiao, Chentao Yang

Abstract: We study the dust properties of 192 nearby galaxies from the JINGLE survey using photometric data in the 22-850micron range. We derive the total dust mass, temperature T and emissivity index beta of the galaxies through the fitting of their spectral energy distribution (SED) using a single modified black-body model (SMBB). We apply a hierarchical Bayesian approach that reduces the known degeneracy… ▽ More We study the dust properties of 192 nearby galaxies from the JINGLE survey using photometric data in the 22-850micron range. We derive the total dust mass, temperature T and emissivity index beta of the galaxies through the fitting of their spectral energy distribution (SED) using a single modified black-body model (SMBB). We apply a hierarchical Bayesian approach that reduces the known degeneracy between T and beta. Applying the hierarchical approach, the strength of the T-beta anti-correlation is reduced from a Pearson correlation coefficient R=-0.79 to R=-0.52. For the JINGLE galaxies we measure dust temperatures in the range 17-30 K and dust emissivity indices beta in the range 0.6-2.2. We compare the SMBB model with the broken emissivity modified black-body (BMBB) and the two modified black-bodies (TMBB) models. The results derived with the SMBB and TMBB are in good agreement, thus applying the SMBB, which comes with fewer free parameters, does not penalize the measurement of the cold dust properties in the JINGLE sample. We investigate the relation between T and beta and other global galaxy properties in the JINGLE and Herschel Reference Survey (HRS) sample. We find that beta correlates with the stellar mass surface density (R=0.62) and anti-correlates with the HI mass fraction (M(HI)/M*, R=-0.65), whereas the dust temperature correlates strongly with the SFR normalized by the dust mass (R=0.73). These relations can be used to estimate T and beta in galaxies with insufficient photometric data available to measure them directly through SED fitting. △ Less

Submitted 11 September, 2019; originally announced September 2019.

Comments: MNRAS in press, 31 pages

Discovery of a giant and luminous Lya+CIV+HeII nebula at z=3.326 with extreme emission line ratios

Authors: R. Marques-Chaves, I. Pérez-Fournon, M. Villar-Martín, R. Gavazzi, D. Riechers, D. Rigopoulou, J. Wardlow, A. Cabrera-Lavers, D. L. Clements, L. Colina, A. Cooray, D. Farrah, R. J. Ivison, C. Jiménez-Ángel, P. Martínez-Navajas, H. Nayyeri, S. Oliver, A. Omont, D. Scott, Y. Shu

Abstract: 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Å, CIV1550Å, and HeII 1640Å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 s… ▽ More 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Å, CIV1550Å, and HeII 1640Å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] △ Less

Submitted 19 July, 2019; originally announced July 2019.

Comments: 13 pages, 8 figures, 3 tables. Accepted for publication in A&A. Abstract abridged for arXiv submission

Journal ref: A&A 629, A23 (2019)

Large Molecular Gas Reservoirs in Star Forming Cluster Galaxies

Authors: Joseph Cairns, Andra Stroe, Carlos De Breuck, Tony Mroczkowski, David L. Clements

Abstract: We present CO(2-1) observations of 72 galaxies in the nearby, disturbed Antlia galaxy cluster with the Atacama Pathfinder Experiment (APEX) telescope. The galaxies in our sample are selected to span a wide range of stellar masses ($10^{8}M_{\odot}\lesssim M_{\star} \lesssim 10^{10}M_{\odot}$) and star formation rates ($0.0005M_{\odot}\text{yr}^{-1}<\text{SFR}<0.3M_{\odot}\text{yr}^{-1}$). Reaching… ▽ More We present CO(2-1) observations of 72 galaxies in the nearby, disturbed Antlia galaxy cluster with the Atacama Pathfinder Experiment (APEX) telescope. The galaxies in our sample are selected to span a wide range of stellar masses ($10^{8}M_{\odot}\lesssim M_{\star} \lesssim 10^{10}M_{\odot}$) and star formation rates ($0.0005M_{\odot}\text{yr}^{-1}<\text{SFR}<0.3M_{\odot}\text{yr}^{-1}$). Reaching a depth of $23\text{mJy}$ in $50\text{km}\text{s}^{-1}$ channels, we report a total CO detection rate of $37.5\%$ and a CO detection rate of $86\%$ for sources within 1 dex of the main sequence. We compare our sample with a similar sample of galaxies in the field, finding that, for a fixed stellar mass and SFR, galaxies in the Antlia cluster have comparable molecular gas reservoirs to field galaxies. We find that $\sim41\%$ (11/27) of our CO detections display non-Gaussian CO(2-1) emission line profiles, and a number of these sources display evidence of quenching in their optical images. We also find that the majority of our sample lie either just below, or far below the main sequence of field galaxies, further hinting at potential ongoing quenching. We conclude that the Antlia cluster represents an intermediate environment between fields and dense clusters, where the gentler intracluster medium (ICM) allows the cluster members to retain their reservoirs of molecular gas, but in which the disturbed ICM is just beginning to influence the member galaxies, resulting in high SFRs and possible ongoing quenching. △ Less

Submitted 28 August, 2019; v1 submitted 17 July, 2019; originally announced July 2019.

Comments: Accepted to ApJ. 46 pages, 14 figures, 5 tables

Journal ref: 2019 ApJ, 882, 132

Confirming Herschel candidate protoclusters from ALMA/VLA CO observations

Authors: C. Gómez-Guijarro, D. A. Riechers, R. Pavesi, G. E. Magdis, T. K. D. Leung, F. Valentino, S. Toft, M. Aravena, S. C. Chapman, D. L. Clements, H. Dannerbauer, S. J. Oliver, I. Pérez-Fournon, I. Valtchanov

Abstract: ALMA 870$μ$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 sys… ▽ More ALMA 870$μ$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$μ$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. △ Less

Submitted 21 May, 2019; originally announced May 2019.

Comments: Accepted by ApJ. 24 pages, 9 figures

Journal ref: The Astrophysical Journal, Volume 872, Issue 2, article id. 117, 21 pp. (2019)

JINGLE, a JCMT legacy survey of dust and gas for galaxy evolution studies: II. SCUBA-2 850 μm data reduction and dust flux density catalogues

Authors: Matthew W. L. Smith, Christopher J. R. Clark, Ilse De Looze, Isabella Lamperti, Amélie Saintonge, Christine D. Wilson, Gioacchino Accurso, Elias Brinks, Martin Bureau, Eun Jung Chung, Phillip J. Cigan, David L. Clements, Thavisha Dharmawardena, Lapo Fanciullo, Yang Gao, Yu Gao, Walter K. Gear, Haley L. Gomez, Joshua Greenslade, Ho Seong Hwang, Francisca Kemper, Jong Chul Lee, Cheng Li, Lihwai Lin, Lijie Liu , et al. (11 additional authors not shown)

Abstract: We present the SCUBA-2 850 $μm$ component of JINGLE, the new JCMT large survey for dust and gas in nearby galaxies, which with 193 galaxies is the largest targeted survey of nearby galaxies at 850 $μm$. We provide details of our SCUBA-2 data reduction pipeline, optimised for slightly extended sources, and including a calibration model adjusted to match conventions used in other far-infrared data.… ▽ More We present the SCUBA-2 850 $μm$ component of JINGLE, the new JCMT large survey for dust and gas in nearby galaxies, which with 193 galaxies is the largest targeted survey of nearby galaxies at 850 $μm$. We provide details of our SCUBA-2 data reduction pipeline, optimised for slightly extended sources, and including a calibration model adjusted to match conventions used in other far-infrared data. We measure total integrated fluxes for the entire JINGLE sample in 10 infrared/submillimetre bands, including all WISE, Herschel-PACS, Herschel-SPIRE and SCUBA-2 850 $μm$ maps, statistically accounting for the contamination by CO(J=3-2) in the 850 $μm$ band. Of our initial sample of 193 galaxies, 191 are detected at 250 $μm$ with a $\geq$ 5$σ$ significance. In the SCUBA-2 850 $μm$ band we detect 126 galaxies with $\geq$ 3$σ$ significance. The distribution of the JINGLE galaxies in far-infrared/sub-millimetre colour-colour plots reveals that the sample is not well fit by single modified-blackbody models that assume a single dust-emissivity index $(β)$. Instead, our new 850 $μm$ data suggest either that a large fraction of our objects require $β< 1.5$, or that a model allowing for an excess of sub-mm emission (e.g., a broken dust emissivity law, or a very cold dust component 10 K) is required. We provide relations to convert far-infrared colours to dust temperature and $β$ for JINGLE-like galaxies. For JINGLE the FIR colours correlate more strongly with star-formation rate surface-density rather than the stellar surface-density, suggesting heating of dust is greater due to younger rather than older stellar-populations, consistent with the low proportion of early-type galaxies in the sample. △ Less

Submitted 23 April, 2019; originally announced April 2019.

Comments: Accepted by MNRAS; data available at http://www.star.ucl.ac.uk/JINGLE/

Mapping Galaxy Clusters in the Distant Universe

Authors: H. Dannerbauer, E. van Kampen, J. Afonso, P. Andreani, F. Arrigoni Battaia, F. Bertoldi, C. Casey, C. -C. Chen, D. L. Clements, C. De Breuck, B. Frye, J. Geach, K. Harrington, M. Hayashi, S. Jin, P. Klaassen, K. Kohno, M. D. Lehnert, I. Matute, T. Mroczkowski, A. Noble, C. Pappalardo, Y. Tamura, J. Zavala

Abstract: 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 t… ▽ More 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. △ Less

Submitted 14 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

The case for a 'sub-millimeter SDSS': a 3D map of galaxy evolution to z~10

Authors: James E. Geach, Manda Banerji, Frank Bertoldi, Matthieu Bethermin, Caitlin M. Casey, Chian-Chou Chen, David L. Clements, Claudia Cicone, Francoise Combes, Christopher Conselice, Asantha Cooray, Kristen Coppin, Emanuele Daddi, Helmut Dannerbauer, Romeel Dave, Matthew Doherty, James S. Dunlop, Alastair Edge, Duncan Farrah, Maximilien Franco, Gary Fuller, Tracy Garratt, Walter Gear, Thomas R. Greve, Evanthia Hatziminaoglou , et al. (31 additional authors not shown)

Abstract: The Sloan Digital Sky Survey (SDSS) was revolutionary because of the extraordinary breadth and ambition of its optical imaging and spectroscopy. We argue that a 'sub-millimeter SDSS' - a sensitive large-area imaging+spectroscopic survey in the sub-mm window - will revolutionize our understanding of galaxy evolution in the early Universe. By detecting the thermal dust continuum emission and atomic… ▽ More The Sloan Digital Sky Survey (SDSS) was revolutionary because of the extraordinary breadth and ambition of its optical imaging and spectroscopy. We argue that a 'sub-millimeter SDSS' - a sensitive large-area imaging+spectroscopic survey in the sub-mm window - will revolutionize our understanding of galaxy evolution in the early Universe. By detecting the thermal dust continuum emission and atomic and molecular line emission of galaxies out to z~10 it will be possible to measure the redshifts, star formation rates, dust and gas content of hundreds of thousands of high-z galaxies down to ~L*. Many of these galaxies will have counterparts visible in the deep optical imaging of the Large Synoptic Survey Telescope. This 3D map of galaxy evolution will span the peak epoch of galaxy formation all the way back to cosmic dawn, measuring the co-evolution of the star formation rate density and molecular gas content of galaxies, tracking the production of metals and charting the growth of large-scale structure. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Science White paper submitted to the Astro2020 Decadal Survey

SILVERRUSH. VIII. Spectroscopic Identifications of Early Large Scale Structures with Protoclusters Over 200 Mpc at z~6-7: Strong Associations of Dusty Star-Forming Galaxies

Authors: Yuichi Harikane, Masami Ouchi, Yoshiaki Ono, Seiji Fujimoto, Darko Donevski, Takatoshi Shibuya, Andreas L. Faisst, Tomotsugu Goto, Bunyo Hatsukade, Nobunari Kashikawa, Kotaro Kohno, Takuya Hashimoto, Ryo Higuchi, Akio K. Inoue, Yen-Ting Lin, Crystal L. Martin, Roderik Overzier, Ian Smail, Jun Toshikawa, Hideki Umehata, Yiping Ao, Scott Chapman, David L. Clements, Myungshin Im, Yipeng Jing , et al. (10 additional authors not shown)

Abstract: We have obtained three-dimensional maps of the universe in $\sim200\times200\times80$ comoving Mpc$^3$ (cMpc$^3$) volumes each at $z=5.7$ and $6.6$ based on a spectroscopic sample of 179 galaxies that achieves $\gtrsim80$\% completeness down to the Ly$α$ luminosity of $\log(L_{\rm Lyα}/[\mathrm{erg\ s^{-1}}])=43.0$, based on our Keck and Gemini observations and the literature. The maps reveal fila… ▽ More We have obtained three-dimensional maps of the universe in $\sim200\times200\times80$ comoving Mpc$^3$ (cMpc$^3$) volumes each at $z=5.7$ and $6.6$ based on a spectroscopic sample of 179 galaxies that achieves $\gtrsim80$\% completeness down to the Ly$α$ luminosity of $\log(L_{\rm Lyα}/[\mathrm{erg\ s^{-1}}])=43.0$, based on our Keck and Gemini observations and the literature. The maps reveal filamentary large-scale structures and two remarkable overdensities made out of at least 44 and 12 galaxies at $z=5.692$ (z57OD) and $z=6.585$ (z66OD), respectively, making z66OD the most distant overdensity spectroscopically confirmed to date with $>10$ spectroscopically confirmed galaxies. We compare spatial distributions of submillimeter galaxies at $z\simeq 4-6$ with our $z=5.7$ galaxies forming the large-scale structures, and detect a $99.97\%$ signal of cross correlation, indicative of a clear coincidence of dusty star-forming galaxy and dust unobscured galaxy formation at this early epoch. The galaxies in z57OD and z66OD are actively forming stars with star formation rates (SFRs) $\gtrsim5$ times higher than the main sequence, and particularly the SFR density in z57OD is 10 times higher than the cosmic average at the redshift (a.k.a. the Madau-Lilly plot). Comparisons with numerical simulations suggest that z57OD and z66OD are protoclusters that are progenitors of the present-day clusters with halo masses of $\sim10^{14}\ \mathrm{M_\odot}$. △ Less

Submitted 24 June, 2019; v1 submitted 25 February, 2019; originally announced February 2019.

Comments: 17 pages, 11 figures, Accepted for Publication in ApJ

The ISM Properties and Gas Kinematics of a Redshift 3 Massive Dusty Star-forming Galaxy

Authors: T. K. Daisy Leung, Dominik A. Riechers, Andrew J. Baker, Dave L. Clements, Asantha Cooray, Christopher C. Hayward, R. J. Ivison, Roberto Neri, Alain Omont, Ismael Perez-Fournon, Douglas Scott, Julie L. Wardlow

Abstract: We present CO(J= 1-0; 3-2; 5-4; 10-9) and 1.2-kpc resolution [CII] line observations of the dusty star-forming galaxy (SFG) HXMM05 -- carried out with the Karl G. Jansky Very Large Array, the Combined Array for Research in Millimeter-wave Astronomy, the Plateau de Bure Interferometer, and the Atacama Large Millimeter/submillimeter Array, measuring an unambiguous redshift of z = 2.9850+/-0.0009. We… ▽ More We present CO(J= 1-0; 3-2; 5-4; 10-9) and 1.2-kpc resolution [CII] line observations of the dusty star-forming galaxy (SFG) HXMM05 -- carried out with the Karl G. Jansky Very Large Array, the Combined Array for Research in Millimeter-wave Astronomy, the Plateau de Bure Interferometer, and the Atacama Large Millimeter/submillimeter Array, measuring an unambiguous redshift of z = 2.9850+/-0.0009. We find that HXMM05 is a hyper-luminous infrared galaxy (LIR=(4+/-1)x10^13 Lsun) with a total molecular gas mass of (2.1+/-0.7)x10^11 (alpha_CO/0.8) Msun. The CO(J=1-0) and [CII] emission are extended over ~9 kpc in diameter, and the CO line FWHM exceeds 1100 km s^-1. The [CII] emission shows a monotonic velocity gradient consistent with a disk, with a maximum rotation velocity of v_c = 616+/-100 km s^-1 and a dynamical mass of (7.7+/-3.1)x10^11 Msun. We find a star formation rate (SFR) of 2900^750_-595 Msun yr^-1. HXMM05 is thus among the most intensely star-forming galaxies known at high redshift. Photo-dissociation region modeling suggests physical conditions similar to nearby SFGs, showing extended star formation, which is consistent with our finding that the gas and dust emission are co-spatial. Its molecular gas excitation resembles the local major merger Arp 220. The broad CO and [CII] lines and a pair of compact dust nuclei suggest the presence of a late-stage major merger at the center of the extended disk, again reminiscent of Arp 220. The observed gas kinematics and conditions together with the presence of a companion and the pair of nuclei suggest that HXMM05 is experiencing multiple mergers as a part of the evolution. △ Less

Submitted 14 December, 2018; originally announced December 2018.

Comments: 32 pages, 19 figures, 7 tables. Accepted to the ApJ

Life Before Fermi - Back to the Solar System

Authors: David L Clements

Abstract: The existence of intelligent, interstellar traveling and colonising life is a key assumption behind the Fermi Paradox. Until recently, detecting signs of life elsewhere has been so technically challenging as to seem almost impossible. However, new observational insights and other developments mean that signs of life elsewhere might realistically be uncovered in the next decade or two. We here revi… ▽ More The existence of intelligent, interstellar traveling and colonising life is a key assumption behind the Fermi Paradox. Until recently, detecting signs of life elsewhere has been so technically challenging as to seem almost impossible. However, new observational insights and other developments mean that signs of life elsewhere might realistically be uncovered in the next decade or two. We here review what are believed to be the basic requirements for life, the history of life on Earth, and then apply this knowledge to potential sites for life in our own Solar System. We conclude that the necessities of life - liquid water and sources of energy - are in fact quite common in the Solar System, but most potential sites are beneath the icy surfaces of gas giant moons. If this is the case elsewhere in the Galaxy, life may be quite common but, even if intelligence develops, is essentially sealed in a finite environment, unable to communicate with the outside world. △ Less

Submitted 15 November, 2018; originally announced November 2018.

Comments: Accepted for publication in the Journal of the British Interplanetary Society

Probing the High-Redshift Universe with SPICA: Toward the Epoch of Reionization and Beyond

Authors: E. Egami, S. Gallerani, R. Schneider, A. Pallottini, L. Vallini, E. Sobacchi, A. Ferrara, S. Bianchi, M. Bocchio, S. Marassi, L. Armus, L. Spinoglio, A. W. Blain, M. Bradford, D. L. Clements, H. Dannerbauer, J. A. Fernández-Ontiveros, E. González-Alfonso, M. J. Griffin, C. Gruppioni, H. Kaneda, K. Kohno, S. C. Madden, H. Matsuhara, P. Najarro , et al. (14 additional authors not shown)

Abstract: 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 ($μ$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $μ$ 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… ▽ More 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 ($μ$$L_{\rm IR}>10^{13}$ L$_{\odot}$, where $μ$ 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 $μ$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. △ Less

Submitted 25 September, 2018; v1 submitted 20 September, 2018; originally announced September 2018.

Comments: 24 pages, 8 figures, accepted for publication in PASA (SPICA special issue)

JINGLE, a JCMT legacy survey of dust and gas for galaxy evolution studies: I. Survey overview and first results

Authors: Amelie Saintonge, Christine D. Wilson, Ting Xiao, Lihwai Lin, Ho Seong Hwang, Tomoka Tosaki, Martin Bureau, Phillip J. Cigan, Christopher J. R. Clark, David L. Clements, Ilse De Looze, Thavisha Dharmawardena, Yang Gao, Walter K. Gear, Joshua Greenslade, Isabella Lamperti, Jong Chul Lee, Cheng Li, Michal J. Michalowski, Angus Mok, Hsi-An Pan, Anne E. Sansom, Mark Sargent, Matthew W. L. Smith, Thomas Williams , et al. (66 additional authors not shown)

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

Submitted 19 September, 2018; originally announced September 2018.

Comments: MNRAS in press, 25 pages

The Herschel-PACS North Ecliptic Pole Survey

Authors: Chris Pearson, Laia Barrufet, Maria Del Carmen Campos Varillas, Stephen Serjeant, David L Clements, Tomotsugu Goto, Myungshin Im, Woong-Seob Jeong, Seong Jin Kim, Hideo Matsuhara, Chris Sedgwick, Ivan Valtchanov

Abstract: A detailed analysis of Herschel-PACS observations at the North Ecliptic Pole is presented. High quality maps, covering an area of 0.44 square degrees, are produced and then used to derive potential candidate source lists. A rigorous quality control pipeline has been used to create final legacy catalogues in the PACS Green 100 micron and Red 160 micron bands, containing 1384 and 630 sources respect… ▽ More A detailed analysis of Herschel-PACS observations at the North Ecliptic Pole is presented. High quality maps, covering an area of 0.44 square degrees, are produced and then used to derive potential candidate source lists. A rigorous quality control pipeline has been used to create final legacy catalogues in the PACS Green 100 micron and Red 160 micron bands, containing 1384 and 630 sources respectively. These catalogues reach to more than twice the depth of the current archival Herschel/PACS Point Source Catalogue, detecting 400 and 270 more sources in the short and long wavelength bands respectively. Galaxy source counts are constructed that extend down to flux densities of 6mJy and 19mJy (50% completeness) in the Green 100 micron and Red 160 micron bands respectively. These source counts are consistent with previously published PACS number counts in other fields across the sky. The source counts are then compared with a galaxy evolution model identifying a population of luminous infrared galaxies as responsible for the bulk of the galaxy evolution over the flux range (5-100mJy) spanned by the observed counts, contributing approximate fractions of 50% and 60% to the cosmic infrared background (CIRB) at 100 microns and 160 microns respectively. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: Accepted for publication in PASJ, 11 pages, 8 figures

AKARI and IRAS: From Beam Corrections to SEDs

Authors: David L. Clements, Michael Rowan-Robinson, Chris Pearson, Jose Afonso, Vianney Labouteiller, Duncan Farrah, Andreas Efstathiou, Josh Greenslade, Lingyu Wang

Abstract: There is significant scientific value to be gained from combining AKARI fluxes with data at other far-IR wavelengths from the IRAS and Herschel missions. To be able to do this we must ensure that there are no systematic differences between the datasets that need to be corrected before the fluxes are compatible with each other. One such systematic effect identified in the BSCv1 data is the issue of… ▽ More There is significant scientific value to be gained from combining AKARI fluxes with data at other far-IR wavelengths from the IRAS and Herschel missions. To be able to do this we must ensure that there are no systematic differences between the datasets that need to be corrected before the fluxes are compatible with each other. One such systematic effect identified in the BSCv1 data is the issue of beam corrections. We determine these for the BSCv2 data by correlating ratios of appropriate IRAS and AKARI bands with the difference in 2MASS J band extended and point source magnitudes for sources cross matched between the IRAS FSC, AKARI BSCv2 and 2MASS catalogs. We find significant correlations (p<< 10^-13) indicating that beam corrections are necessary in the 65 and 90 micron bands. We then use these corrected fluxes to supplement existing data in spectral energy distribution (SED) fits for ultraluminous infrared galaxies in the HERUS survey. The addition of AKARI fluxes makes little difference to the results of simple (T,beta) fits to the SEDs of these sources, though there is a general decrease in reduced chi^2 values. The utility of the extra AKARI data, however, is in allowing physically more realistic SED models with more parameters to be fit to the data. We also extend our analysis of beam correction issues in the AKARI data by examining the Herschel Reference Sample, which have Herschel photometry from 100 to 500 microns and which are more spatially extended than the HERUS ULIRGs. Thirty-four of the HRS sources have good Herschel SEDs and matching data from AKARI. This investigation finds that our simple 2MASS-based beam correction scheme is inadequate for these larger and more complex sources. There are also indications that additional beam corrections at 140 and 160 microns are needed for these sources, extended on scales >1 arcminute. △ Less

Submitted 10 August, 2018; originally announced August 2018.

Comments: Accepted for publication in the PASJ special issue on AKARI. 3 Figures