Search

The discovery of phosphine in Venus' atmosphere provides lessons for the search for life. The detection has survived all challenges and has acquired independent support from archival data from PVP. The presence of phosphine in Venus' oxidising environment is perplexing, and comprehensive studies rule out all known abiotic sources. More data is needed to understand the origin of phosphine, leading to JCMT-Venus, a long term atmospheric monitoring programme. This can find how phosphine varies in relation to other species providing clues to its origin. We present the latest JCMT-Venus results. The discovery and subsequent papers were explicit that they did not constitute evidence for life, only of phosphine. Media and public reaction to the discovery and its implications provide lessons for future life searches, as does the reaction of the scientific community. How this was handled by the team, media, and general public will be reviewed., Comment: Accepted for publication in Proceedings of the KAVLI-IAU SYMPOSIUM (IAUS 387): (Toward) Discovery of Life Beyond Earth and its Impact

We present confusion-limited SCUBA-2 450-$\mu$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 $\mu$m cover an area of 450 arcmin$^2$. We achieved instrumental noise levels of $\sigma_{\mathrm{450}}=$ 0.59 mJy beam$^{-1}$ and $\sigma_{\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) $\sigma$ threshold, we detected 360 (479) sources at 450 $\mu$m and 237 (314) sources at 850 $\mu$m. We derive the deepest blank-field number counts at 450 $\mu$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 $\mu$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 $\mu$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-$\mu$m extragalactic background light (EBL) measured by COBE and Planck. Our results suggest that the 450-$\mu$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., Comment: 29 pages, 14 figures, accepted for publication in ApJ

We have carried out HST snapshot observations at 1.1 $\mu$m of 281 candidate strongly lensed galaxies identified in the wide-area extragalactic surveys conducted with the Herschel space observatory. Our candidates comprise systems with flux densities at $500\,\mu$m$ S_{500}\geq 80$ mJy. We model and subtract the surface brightness distribution for 130 systems, where we identify a candidate for the foreground lens candidate. After combining visual inspection, archival high-resolution observations, and lens subtraction, we divide the systems into different classes according to their lensing likelihood. We confirm 65 systems to be lensed. Of these, 30 are new discoveries. We successfully perform lens modelling and source reconstruction on 23 systems, where the foreground lenses are isolated galaxies and the background sources are detected in the HST images. All the systems are successfully modelled as a singular isothermal ellipsoid. The Einstein radii of the lenses and the magnifications of the background sources are consistent with previous studies. However, the background source circularised radii (between 0.34 kpc and 1.30 kpc) are $\sim$3 times smaller than the ones measured in the sub-mm/mm for a similarly selected and partially overlapping sample. We compare our lenses with those in the SLACS survey, confirming that our lens-independent selection is more effective at picking up fainter and diffuse galaxies and group lenses. This sample represents the first step towards characterising the near-IR properties and stellar masses of the gravitationally lensed dusty star-forming galaxies., Comment: 57 pages, 18 figures, 11 tables. Accepted for publication in MNRAS

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{\mu}$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{\mu}$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., Comment: 9 pages, 5 figures, 4 tables; added figure, updated Discussion section, and minor revisions; accepted for publication in MNRAS

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

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., Comment: Invited review, accepted for publication in Contemporary Physics

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., Comment: ApJ, accepted

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{\sigma} 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., Comment: v2 is an updated and revised version of the paper

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., Comment: MNRAS Accepted

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

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$\sigma$ 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 $32$. Analysis of MORA sources' spectral energy distributions hint at steeper empirically-measured dust emissivity indices than typical literature studies, with $\langle\beta\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., Comment: 42 pages, 14 figures, accepted for publication in ApJ

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., Comment: 26 pages, 19 figures. Submitted to ApJS June 2021, Accepted September 2021

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., Comment: v3 is an proofread version of the paper that will appear in print in MNRAS

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{\sigma} confidence using Venus bandpass self-calibration, and two simpler approaches are shown to yield example 4.5-4.8{\sigma} 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., Comment: Revised and updated re-analysis of the phosphine spectra in Venus' clouds; Submitted to Nature Astronomy "Matters Arising" on Apr/19/2021

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{\sigma} 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., Comment: 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

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{\sigma} 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., Comment: 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

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., Comment: 20 pages, 11 figures, 6 tables; to appear in the Astrophysical Journal (accepted November 29, 2020)

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.

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

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., Comment: 10 pages 3 figures. EAO Submillimetre Futures White Paper Series, 2019

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~\mu$m wavelength range is successfully detected with $>5\sigma$ 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 $\mu$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., Comment: 26 pages, 12 figures, Accepted for publication in ApJ

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$\mu$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 $\tau_{\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)., Comment: 24 pages, 8 figures, accepted for publication in ApJ

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., Comment: MNRAS in press, 31 pages

The largest Herschel extragalactic surveys, H-ATLAS and HerMES, have selected a sample of "ultrared" dusty, star-forming galaxies (DSFGs) with rising SPIRE flux densities ($S_{500} > S_{350} > S_{250}$; so-called "500 $\mu$m-risers") as an efficient way for identifying DSFGs at higher redshift ($z > 4$). In this paper, we present a large Spitzer follow-up program of 300 Herschel ultrared DSFGs. We have obtained high-resolution ALMA, NOEMA, and SMA data for 63 of them, which allow us to securely identify the Spitzer/IRAC counterparts and classify them as gravitationally lensed or unlensed. Within the 63 ultrared sources with high-resolution data, $\sim$65% appear to be unlensed, and $\sim$27% are resolved into multiple components. We focus on analyzing the unlensed sample by directly performing multi-wavelength spectral energy distribution (SED) modeling to derive their physical properties and compare with the more numerous $z \sim 2$ DSFG population. The ultrared sample has a median redshift of 3.3, stellar mass of 3.7 $\times$ 10$^{11}$ $M_{\odot}$, star formation rate (SFR) of 730 $M_{\odot}$yr$^{-1}$, total dust luminosity of 9.0 $\times$ 10$^{12}$ $L_{\odot}$, dust mass of 2.8 $\times$ 10$^9$ $M_{\odot}$, and V-band extinction of 4.0, which are all higher than those of the ALESS DSFGs. Based on the space density, SFR density, and stellar mass density estimates, we conclude that our ultrared sample cannot account for the majority of the star-forming progenitors of the massive, quiescent galaxies found in infrared surveys. Our sample contains the rarer, intrinsically most dusty, luminous and massive galaxies in the early universe that will help us understand the physical drivers of extreme star formation., Comment: 33 pages, 10 figures; accepted for publication in ApJS. The catalog tables will be available on ApJS and VizieR; authors' version is available now upon request

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., Comment: Accepted to ApJ. 46 pages, 14 figures, 5 tables

We present the SCUBA-2 850 ${\mu}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 ${\mu}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 ${\mu}m$ maps, statistically accounting for the contamination by CO(J=3-2) in the 850 ${\mu}m$ band. Of our initial sample of 193 galaxies, 191 are detected at 250 ${\mu}m$ with a $\geq$ 5${\sigma}$ significance. In the SCUBA-2 850 ${\mu}m$ band we detect 126 galaxies with $\geq$ 3${\sigma}$ 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 $(\beta)$. Instead, our new 850 ${\mu}m$ data suggest either that a large fraction of our objects require $\beta < 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 $\beta$ 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., Comment: Accepted by MNRAS; data available at http://www.star.ucl.ac.uk/JINGLE/

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., Comment: Science White paper submitted to the Astro2020 Decadal Survey

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$\alpha$ luminosity of $\log(L_{\rm Ly\alpha}/[\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}$., Comment: 17 pages, 11 figures, Accepted for Publication in ApJ

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., Comment: Accepted for publication in the Journal of the British Interplanetary Society

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., Comment: MNRAS in press, 25 pages

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., Comment: Accepted for publication in PASJ, 11 pages, 8 figures

The galaxies found in optical surveys fall in two distinct regions of a diagram of optical colour versus absolute magnitude: the red sequence and the blue cloud with the green valley in between. We show that the galaxies found in a submillimetre survey have almost the opposite distribution in this diagram, forming a `green mountain'. We show that these distinctive distributions follow naturally from a single, continuous, curved Galaxy Sequence in a diagram of specific star-formation rate versus stellar mass without there being the need for a separate star-forming galaxy Main Sequence and region of passive galaxies. The cause of the red sequence and the blue cloud is the geometric mapping between stellar mass/specific star-formation rate and absolute magnitude/colour, which distorts a continuous Galaxy Sequence in the diagram of intrinsic properties into a bimodal distribution in the diagram of observed properties. The cause of the green mountain is Malmquist bias in the submillimetre waveband, with submillimetre surveys tending to select galaxies on the curve of the Galaxy Sequence, which have the highest ratios of submillimetre-to-optical luminosity. This effect, working in reverse, causes galaxies on the curve of the Galaxy Sequence to be underrepresented in optical samples, deepening the green valley. The green valley is therefore not evidence (1) for there being two distinct populations of galaxies, (2) for galaxies in this region evolving more quickly than galaxies in the blue cloud and the red sequence, (c) for rapid quenching processes in the galaxy population., Comment: Accepted by Monthly Notices of the Royal Astronomical Society

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., Comment: Accepted for publication in the PASJ special issue on AKARI. 3 Figures

The submillimeter galaxy (SMG) HERMES J105751.1+573027 (hereafter HLock01) at z = 2.9574 +/- 0.0001 is one of the brightest gravitationally lensed sources discovered in the Herschel Multi-tiered Extragalactic Survey. Apart from the high flux densities in the far-infrared, it is also extremely bright in the rest-frame ultraviolet (UV), with a total apparent magnitude m_UV = 19.7 mag. We report here deep spectroscopic observations with the Gran Telescopio Canarias of the optically bright lensed images of HLock01. Our results suggest that HLock01 is a merger system composed of the Herschel-selected SMG and an optically bright Lyman break-like galaxy (LBG), separated by only 3.3 kpc in projection. While the SMG appears very massive (M* = 5x10^11 Msun), with a highly extinguished stellar component (A_V = 4.3), the LBG is a young, lower-mass (M* = 1x10^10 Msun), but still luminous (10xL*_UV) satellite galaxy. Detailed analysis of the high signal-to-noise (S/N) rest-frame UV spectrum of the LBG shows complex kinematics of the gas, exhibiting both blueshifted and redshifted absorption components. While the blueshifted component is associated with strong galactic outflows from the massive stars in the LBG, as is common in most star-forming galaxies, the redshifted component may be associated with gas inflow seen along a favorable sightline to the LBG. We also find evidence of an extended gas reservoir around HLock01 at an impact parameter of 110 kpc, through the detection of C II 1334$\AA$ absorption in the red wing of a bright Ly-alpha emitter at z = 3.327. The data presented here highlight the power of gravitational lensing in high S/N studies to probe deeply into the physics of high-z star forming galaxies., Comment: 24 pages, 11 figures. Accepted for publication in ApJ