Showing 1–50 of 154 results for author: Hodge, J

A photo-$z$ cautionary tale: Redshift confirmation of COSBO-7 at $z=2.625$

Authors: Shuowen Jin, Nikolaj B. Sillassen, Jacqueline Hodge, Georgios E. Magdis, Caitlin Casey, Francesca Rizzo, Anton M. Koekemoer, Francesco Valentino, Vasily Kokorev, Benjamin Magnelli, Raphael Gobat, Steven Gillman, Maximilien Franco, Andreas Faisst, Jeyhan Kartaltepe, Eva Schinnerer, Sune Toft, Hiddo S. B. Algera, Santosh Harish, Minju Lee, Daizhong Liu, Marko Shuntov, Margherita Talia, Aswin Vijayan

Abstract: Photometric redshifts are widely used in studies of dusty star-forming galaxies (DSFGs), but catastrophic photo-$z$ failure can undermine all redshift-dependent results. Here we report the spectroscopic redshift confirmation of COSBO-7, a strongly lensed DSFG in the COSMOS-PRIMER field. Recently, using 10 bands of JWST NIRCam and MIRI imaging data on COSBO-7, Ling et al. (2024) reported a photomet… ▽ More Photometric redshifts are widely used in studies of dusty star-forming galaxies (DSFGs), but catastrophic photo-$z$ failure can undermine all redshift-dependent results. Here we report the spectroscopic redshift confirmation of COSBO-7, a strongly lensed DSFG in the COSMOS-PRIMER field. Recently, using 10 bands of JWST NIRCam and MIRI imaging data on COSBO-7, Ling et al. (2024) reported a photometric redshift solution of $z\gtrsim7.0$, favored by four independent spectral energy distribution (SED) fitting codes, and therefore providing an appealing candidate of the most distant massive DSFG. This photo-$z$ solution was also supported by a single line detection in ALMA Band 3 consistent with CO(7-6) at $z=7.46$. However, our new ALMA observations robustly detect two lines in Band 6 identified as CO(7-6) and [CI](2-1) at $z_{\rm spec}=2.625$, and thus the Band 3 line as CO(3-2). The three robust line detections decidedly place COSBO-7 at $z=2.625$, refuting the photo-$z$ solution. We derive physical parameters by fitting NIR-to-mm photometry and lens modeling, revealing that COSBO-7 is a main-sequence galaxy. We examine possible reasons for this photo-$z$ failure and attribute it to (1) the likely underestimation of photometric uncertainty at 0.9$μ$m, and (2) the lack of photometry at wavelengths beyond 20$μ$m. Notably, we recover a bona-fide $z_{\rm phot}\sim 2.3$ by including the existing MIPS $24μ$m photometry, demonstrating the critical importance of mid-infrared data for photo-$z$ robustness. This work highlights a common challenge in modeling SEDs of DSFGs, cautioning against the reliability of photometric redshifts, as well as pseudo-spectroscopic redshifts based on single line detection. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: Submitted to A&A Letter

REBELS-25: Discovery of a dynamically cold disc galaxy at z = 7.31

Authors: Lucie E. Rowland, Jacqueline Hodge, Rychard Bouwens, Pavel Mancera Piña, Alexander Hygate, Hiddo Algera, Manuel Aravena, Rebecca Bowler, Elisabete da Cunha, Pratika Dayal, Andrea Ferrara, Thomas Herard-Demanche, Hanae Inami, Ivana van Leeuwen, Ilse de Looze, Pascal Oesch, Andrea Pallottini, Siân Phillips, Matus Rybak, Sander Schouws, Renske Smit, Laura Sommovigo, Mauro Stefanon, Paul van der Werf

Abstract: We present high resolution ($\sim0.14$" = 710 pc) ALMA [CII] 158$μ$m and dust continuum follow-up observations of REBELS-25, a [CII]-luminous ($L_{\mathrm{[CII]}}=(1.7\pm0.2)\times 10^9 \mathrm{L_{\odot}}$) galaxy at redshift $z=7.3065\pm0.0001$. These high resolution, high signal-to-noise observations allow us to study the sub-kpc morphology and kinematics of this massive (… ▽ More We present high resolution ($\sim0.14$" = 710 pc) ALMA [CII] 158$μ$m and dust continuum follow-up observations of REBELS-25, a [CII]-luminous ($L_{\mathrm{[CII]}}=(1.7\pm0.2)\times 10^9 \mathrm{L_{\odot}}$) galaxy at redshift $z=7.3065\pm0.0001$. These high resolution, high signal-to-noise observations allow us to study the sub-kpc morphology and kinematics of this massive ($M_* = 8^{+4}_{-2} \times 10^9 \mathrm{M_{\odot}}$) star-forming (SFR$_{\mathrm{UV+IR}} = 199^{+101}_{-63} \mathrm{M_{\odot}} \mathrm{yr}^{-1}$) galaxy in the Epoch of Reionisation. By modelling the kinematics with $^{\mathrm{3D}}$BAROLO, we find it has a low velocity dispersion ($\barσ = 33 \pm 9$ km s$^{-1}$) and a high ratio of ordered-to-random motion ($V_{\mathrm{rot, ~max}}/\barσ = 11 ^{+8}_{-4}$), indicating that REBELS-25 is a dynamically cold disc. Additionally, we find that the [CII] distribution is well fit by a near-exponential disc model, with a Sérsic index, $n$, of $1.3 \pm 0.2$, and we see tentative evidence of more complex non-axisymmetric structures suggestive of a bar in the [CII] and dust continuum emission. By comparing to other high spatial resolution cold gas kinematic studies, we find that dynamically cold discs seem to be more common in the high redshift Universe than expected based on prevailing galaxy formation theories, which typically predict more turbulent and dispersion-dominated galaxies in the early Universe as an outcome of merger activity, gas accretion and more intense feedback. This higher degree of rotational support seems instead to be consistent with recent cosmological simulations that have highlighted the contrast between cold and warm ionised gas tracers, particularly for massive galaxies. We therefore show that dynamically settled disc galaxies can form as early as 700 Myr after the Big Bang. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: Submitted to MNRAS

A Comparative Study of the Ground State Transitions of CO and [C I] as Molecular Gas Tracers at High Redshift

Authors: Marta Frias Castillo, Matus Rybak, Jacqueline A. Hodge, Paul Van der Werk, Ian Smail, Joshua Butterworth, Jasper Jansen, Theodoros Topkaras, Chian-Chou Chen, Scott C. Chapman, Axel Weiss, Hiddo Algera, Jack E. Birkin, Elisabete da Cunha, Jianhang Chen, Helmut Dannerbauer, E. F. Jiménez-Andrade, Soh Ikarashi, Cheng-Lin Liao, Eric J. Murphy, A. M. Swinbank, Fabian Walter, Gabriela Calistro Rivera, R. J. Ivison, Claudia del P. Lagos

Abstract: The CO(1--0) and [\ion{C}{1}](1--0) emission lines are well-established tracers of cold molecular gas mass in local galaxies. At high redshift, where the interstellar medium (ISM) is likely to be denser, there have been limited direct comparisons of both ground state transitions. Here we present a study of CO(1--0) and [\ion{C}{1}](1--0) emission in a sample of 20 unlensed dusty, star-forming gala… ▽ More The CO(1--0) and [\ion{C}{1}](1--0) emission lines are well-established tracers of cold molecular gas mass in local galaxies. At high redshift, where the interstellar medium (ISM) is likely to be denser, there have been limited direct comparisons of both ground state transitions. Here we present a study of CO(1--0) and [\ion{C}{1}](1--0) emission in a sample of 20 unlensed dusty, star-forming galaxies at $z=2-5$. The CO(1--0)/[\ion{C}{1}](1--0) ratio is constant up to at least $z=5$, supporting the use of [CI](1-0) as a gas mass tracer. PDR modelling of the available data indicates a median H$_2$ density of log$(n~[$cm$^{-3}])=4.7\pm0.2$, and UV radiation field log$(G_{\mathrm{UV}} [G$_0$])=3.2\pm0.2$. We use the CO(1--0), [\ion{C}{1}](1--0) and 3mm dust continuum measurements to cross--calibrate the respective gas mass conversion factors, finding no dependence of these factors on either redshift or infrared luminosity. Assuming a variable CO conversion factor then implies [\ion{C}{1}] and dust conversion factors that differ from canonically assumed values but are consistent with the solar/super-solar metallicities expected for our sources. Radiative transfer modelling shows that the warmer CMB at high redshift can significantly affect the [\ion{C}{1}] as well as CO emission, which can change the derived molecular gas masses by up to 70\% for the coldest kinetic gas temperatures expected. Nevertheless, we show that the magnitude of the effect on the ratio of the tracers is within the known scatter of the $L'_\mathrm{CO}-L'_\mathrm{[CI]}$ relation. Further determining the absolute decrease of individual line intensities will require well-sampled spectral line energy distributions (SLEDs) to model the gas excitation conditions in more detail. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Optical and soft X-ray light-curve analysis during the 2022 eruption of U Scorpii: structural changes in the accretion disk

Authors: Katsuki Muraoka, Naoto Kojiguchi, Junpei Ito, Daisaku Nogami, Taichi Kato, Yusuke Tampo, Kenta Taguchi, Keisuke Isogai, Teofilo Arranz, John Blackwell, David Blane, Stephen M. Brincat, Graeme Coates, Walter Cooney, Shawn Dvorak, Charles Galdies, Daniel Glomski, Franz-Josef Hambsch, Barbara Harris, John Hodge, Jose L. Hernández-Verdejo, Marco Iozzi, Hiroshi Itoh, Seiichiro Kiyota, Darrell Lee , et al. (30 additional authors not shown)

Abstract: We present our optical photometric observations of the 2022 eruption of the recurrent nova U Scorpii (U Sco) using 49,152 data points over 70 d following the optical peak. We have also analyzed its soft X-ray (0.3--1 keV) light curve by the Neil Gehrels Swift Observatory. During the 2022 eruption, the optical plateau stage started 13.8--15.0 d and ended 23.8--25.0 d after the optical peak. The sof… ▽ More We present our optical photometric observations of the 2022 eruption of the recurrent nova U Scorpii (U Sco) using 49,152 data points over 70 d following the optical peak. We have also analyzed its soft X-ray (0.3--1 keV) light curve by the Neil Gehrels Swift Observatory. During the 2022 eruption, the optical plateau stage started 13.8--15.0 d and ended 23.8--25.0 d after the optical peak. The soft X-ray stage started 14.6--15.3 d and ended 38.7--39.5 d after the optical peak. Both stages started later and had shorter durations, and the soft X-ray light curve peaked earlier and was less luminous compared to those during the U Sco 2010 eruption. These points suggest that there were differences in the envelope mass between the different cycles of the nova eruption. Furthermore, we have analyzed the optical eclipses during the 2022 eruption. The primary eclipse was first observed 10.4--11.6 d after the optical peak, earlier than the beginning of the optical plateau stage. This sequence of events can be explained by the receding ejecta photosphere associated with the expanding nova ejecta. We have determined the ingress and egress phases of the primary eclipses and estimated the outer radius of the optical light source centered at the white dwarf (WD). During the optical plateau stage, the source radius remained $\sim$1.2 times larger than the Roche volume radius of the primary WD, being close to the L1 point. When the optical plateau stage ended, the source radius drastically shrank to the tidal truncation radius within a few orbital periods. This previously unresolved phenomenon can be interpreted as a structural change in U Sco where the temporarily expanded accretion disk due to the nova wind returned to a steady state. △ Less

Submitted 13 February, 2024; v1 submitted 9 February, 2024; originally announced February 2024.

Comments: 16 pages, 7 figures, 7 tables, accepted for publication in PASJ; doi:10.1093/pasj/psae010

MSC Class: 85-11

At the end of cosmic noon: Short gas depletion times in unobscured quasars at $z \sim$ 1

Authors: M. Frias Castillo, M. Rybak, J. Hodge, P. van der Werf, L. J. Abbo, F. J. Ballieux, S. Ward, C. Harrison, G. Calistro Rivera, J. P. McKean, H. R. Stacey

Abstract: Unobscured quasars (QSOs) are predicted to be the final stage in the evolutionary sequence from gas-rich mergers to gas-depleted, quenched galaxies. Studies of this population, however, find a high incidence of far-infrared-luminous sources -suggesting significant dust-obscured star formation-but direct observations of the cold molecular gas fuelling this star formation are still necessary. We pre… ▽ More Unobscured quasars (QSOs) are predicted to be the final stage in the evolutionary sequence from gas-rich mergers to gas-depleted, quenched galaxies. Studies of this population, however, find a high incidence of far-infrared-luminous sources -suggesting significant dust-obscured star formation-but direct observations of the cold molecular gas fuelling this star formation are still necessary. We present a NOEMA study of CO(2-1) emission, tracing the cold molecular gas, in ten lensed z=1-1.5 unobscured QSOs. We detected CO(2-1) in seven of our targets, four of which also show continuum emission (λ_rest = 1.3mm). After subtracting the foreground galaxy contribution to the photometry, spectral energy distribution fitting yielded stellar masses of 10^9-11 M_\odot, with star formation rates of 25-160 M_\odot yr^-1 for the host galaxies. These QSOs have lower $L'_\mathrm{CO}$ than star-forming galaxies with the same L_IR, and show depletion times spanning a large range (50-900 Myr), but with a median of just 90 Myr. We find molecular gas masses in the range 2-40 x 10^9(alpha_CO/4) M_\odot, which suggest gas fractions above ~50% for most of the targets. Despite the presence of an unobscured QSO, the host galaxies are able to retain significant amounts of cold gas. However, with a median depletion time of ~90 Myr, the intense burst of star formation taking place in these targets will quickly deplete their molecular gas reservoirs in the absence of gas replenishment, resulting in a quiescent host galaxy. The non-detected QSOs are three of the four radio-loud QSOs in the sample, and their properties indicate that they are likely already transitioning into quiescence. Recent cosmological simulations tend to overestimate the depletion times expected for these z~1 QSO-host galaxies, which is likely linked to their difficulty producing starbursts across the general high-redshift galaxy population. △ Less

Submitted 12 January, 2024; v1 submitted 9 January, 2024; originally announced January 2024.

Comments: 20 pages. Accepted for publication in A&A

The kinematics of massive high-redshift dusty star-forming galaxies

Authors: A. Amvrosiadis, J. L. Wardlow, J. E. Birkin, I. Smail, A. M. Swinbank, J. Nightingale, F. Bertoldi, W. N. Brandt, C. M. Casey, S. C. Chapman, C. -C. Chen, P. Cox, E. da Cunha, H. Dannerbauer, U. Dudzevičiūtė, B. Gullberg, J. A. Hodge, K. K. Knudsen, K. Menten, F. Walter, P. van der Werf

Abstract: We present a new method for modelling the kinematics of galaxies from interferometric observations by performing the optimization of the kinematic model parameters directly in visibility-space instead of the conventional approach of fitting velocity fields produced with the CLEAN algorithm in real-space. We demonstrate our method on ALMA observations of $^{12}$CO (2$-$1), (3$-$2) or (4$-$3) emissi… ▽ More We present a new method for modelling the kinematics of galaxies from interferometric observations by performing the optimization of the kinematic model parameters directly in visibility-space instead of the conventional approach of fitting velocity fields produced with the CLEAN algorithm in real-space. We demonstrate our method on ALMA observations of $^{12}$CO (2$-$1), (3$-$2) or (4$-$3) emission lines from an initial sample of 30 massive 850$μ$m-selected dusty star-forming galaxies with far-infrared luminosities $\gtrsim$$\,10^{12}\,$L$_{\odot}$ in the redshift range $z \sim\,$1.2$-$4.7. Using the results from our modelling analysis for the 12 sources with the highest signal-to-noise emission lines and disk-like kinematics, we conclude the following: (i) Our sample prefers a CO-to-$H_2$ conversion factor, of $α_{\rm CO} = 0.92 \pm 0.36$; (ii) These far-infrared luminous galaxies follow a similar Tully$-$Fisher relation between the circularized velocity, $V_{\rm circ}$, and baryonic mass, $M_{\rm b}$, as more typical star-forming samples at high redshift, but extend this relation to much higher masses $-$ showing that these are some of the most massive disk-like galaxies in the Universe; (iii) Finally, we demonstrate support for an evolutionary link between massive high-redshift dusty star-forming galaxies and the formation of local early-type galaxies using the both the distributions of the baryonic and kinematic masses of these two populations on the $M_{\rm b}\,-\,σ$ plane and their relative space densities. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Metal and dust evolution in ALMA REBELS galaxies: insights for future JWST observations

Authors: Marco Palla, Ilse De Looze, Monica Relaño, Stefan van der Giessen, Pratika Dayal, Andrea Ferrara, Raffaella Schneider, Luca Graziani, Hiddo S. B. Algera, Manuel Aravena, Rebecca A. A. Bowler, Alexander P. S. Hygate, Hanae Inami, Ivana van Leeuwen, Rychard Bouwens, Jacqueline Hodge, Renske Smit, Mauro Stefanon, Paul van der Werf

Abstract: ALMA observations revealed the presence of significant amounts of dust in the first Gyr of Cosmic time. However, the metal and dust buildup picture remains very uncertain due to the lack of constraints on metallicity. JWST has started to reveal the metal content of high-redshift targets, which may lead to firmer constraints on high-redshift dusty galaxies evolution. In this work, we use detailed c… ▽ More ALMA observations revealed the presence of significant amounts of dust in the first Gyr of Cosmic time. However, the metal and dust buildup picture remains very uncertain due to the lack of constraints on metallicity. JWST has started to reveal the metal content of high-redshift targets, which may lead to firmer constraints on high-redshift dusty galaxies evolution. In this work, we use detailed chemical and dust evolution models to explore the evolution of galaxies within the ALMA REBELS survey, testing different metallicity scenarios that could be inferred from JWST observations. In the models, we track the buildup of stellar mass by using non-parametric SFHs for REBELS galaxies. Different scenarios for metal and dust evolution are simulated by allowing different prescriptions for gas flows and dust processes. The model outputs are compared with measured dust scaling relations, by employing metallicity-dependent calibrations for the gas mass based on the [CII]158micron line. Independently of the galaxies metal content, we found no need for extreme dust prescriptions to explain the dust masses revealed by ALMA. However, different levels of metal enrichment will lead to different dominant dust production mechanisms, with stardust production dominant over other ISM dust processes only in the metal-poor case. This points out how metallicity measurements from JWST will significantly improve our understanding of the dust buildup in high-redshift galaxies. We also show that models struggle to reproduce observables such as dust-to-gas and dust-to-stellar ratios simultaneously, possibly indicating an overestimation of the gas mass through current calibrations, especially at high metallicities. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 16 pages + appendices, 9 Figures, 1 Table. Resubmitted to MNRAS after moderate revision

TESS photometry of the nova eruption in V606 Vul: asymmetric photosphere and multiple ejections?

Authors: Kirill V. Sokolovsky, Elias Aydi, Konstantin Malanchev, Colin J. Burke, Koji Mukai, J. L. Sokoloski, Brian D. Metzger, Kirill E. Atapin, Aleksandre A. Belinski, Yu-Ching Chen, Laura Chomiuk, Pavol A. Dubovsky, Claude-Andre Faucher-Giguere, Rebekah A. Hounsell, Natalia P. Ikonnikova, Vsevolod Yu. Lander, Junyao Li, Justin D. Linford, Amy J. Mioduszewski, Isabella Molina, Ulisse Munari, Sergey A. Potanin, Robert M. Quimby, Michael P. Rupen, Simone Scaringi , et al. (46 additional authors not shown)

Abstract: Lightcurves of many classical novae deviate from the canonical "fast rise - smooth decline" pattern and display complex variability behavior. We present the first TESS-space-photometry-based investigation of this phenomenon. We use TESS Sector 41 full-frame images to extract a lightcurve of the slow Galactic nova V606 Vul that erupted nine days prior to the start of the TESS observations. The ligh… ▽ More Lightcurves of many classical novae deviate from the canonical "fast rise - smooth decline" pattern and display complex variability behavior. We present the first TESS-space-photometry-based investigation of this phenomenon. We use TESS Sector 41 full-frame images to extract a lightcurve of the slow Galactic nova V606 Vul that erupted nine days prior to the start of the TESS observations. The lightcurve covers the first of two major peaks of V606 Vul that was reached 19 days after the start of the eruption. The nova reached its brightest visual magnitude V=9.9 in its second peak 64 days after the eruption onset, following the completion of Sector 41 observations. To increase the confidence level of the extracted lightcurve, we performed the analysis using four different codes implementing the aperture photometry (Lightkurve, VaST) and image subtraction (TESSreduce, tequila_shots) and find good agreement between them. We performed ground-based photometric and spectroscopic monitoring to complement the TESS data. The TESS lightcurve reveals two features: periodic variations (0.12771d, 0.01mag average peak-to-peak amplitude) that disappeared when the source was within 1mag of peak optical brightness and a series of isolated mini-flares (with peak-to-peak amplitudes of up to 0.5mag) appearing at seemingly random times. We interpret the periodic variations as the result of azimuthal asymmetry of the photosphere engulfing the nova-hosting binary that was distorted by and rotating with the binary. Whereas we use spectra to associate the two major peaks in the nova lightcurve with distinct episodes of mass ejection, the origin of mini-flares remains elusive. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 31 pages, 10 figures, submitted to ApJ; comments welcome

The ALMA REBELS survey: obscured star formation in massive Lyman-break galaxies at z = 4-8 revealed by the IRX-$β$ and $M_{\star}$ relations

Authors: R. A. A. Bowler, H. Inami, L. Sommovigo, R. Smit, H. S. B. Algera, M. Aravena, L. Barrufet, R. Bouwens, E. da Cunha, F. Cullen, P. Dayal, I. de Looze, J. S. Dunlop, Y. Fudamoto, V. Mauerhofer, R. J. McLure, M. Stefanon, R. Schneider, A. Ferrara, L. Graziani, J. A. Hodge, T. Nanayakkara, M. Palla, S. Schouws, D. P. Stark , et al. (1 additional authors not shown)

Abstract: We investigate the degree of dust obscured star formation in 49 massive (${\rm log}_{10}(M_{\star}/{\rm M}_{\odot})>9$) Lyman-break galaxies (LBGs) at $z = 6.5$-$8$ observed as part of the ALMA Reionization Era Bright Emission Line Survey (REBELS) large program. By creating deep stacks of the photometric data and the REBELS ALMA measurements we determine the average rest-frame UV, optical and far-… ▽ More We investigate the degree of dust obscured star formation in 49 massive (${\rm log}_{10}(M_{\star}/{\rm M}_{\odot})>9$) Lyman-break galaxies (LBGs) at $z = 6.5$-$8$ observed as part of the ALMA Reionization Era Bright Emission Line Survey (REBELS) large program. By creating deep stacks of the photometric data and the REBELS ALMA measurements we determine the average rest-frame UV, optical and far-infrared (FIR) properties which reveal a significant fraction ($f_{\rm obs} = 0.4$-$0.7$) of obscured star formation, consistent with previous studies. From measurements of the rest-frame UV slope, we find that the brightest LBGs at these redshifts show bluer ($β\simeq -2.2$) colours than expected from an extrapolation of the colour-magnitude relation found at fainter magnitudes. Assuming a modified blackbody spectral-energy distribution (SED) in the FIR (with dust temperature of $T_{\rm d} = 46\,{\rm K}$ and $β_{\rm d} = 2.0$), we find that the REBELS sources are in agreement with the local ''Calzetti-like'' starburst Infrared-excess (IRX)-$β$ relation. By reanalysing the data available for 108 galaxies at $z \simeq 4$-$6$ from the ALPINE ALMA large program using a consistent methodology and assumed FIR SED, we show that from $z \simeq 4$-$8$, massive galaxies selected in the rest-frame UV have no appreciable evolution in their derived IRX-$β$ relation. When comparing the IRX-$M_{\star}$ relation derived from the combined ALPINE and REBELS sample to relations established at $z < 4$, we find a deficit in the IRX, indicating that at $z > 4$ the proportion of obscured star formation is lower by a factor of $\gtrsim 3$ at a given a $M_{\star}$. Our IRX-$β$ results are in good agreement with the high-redshift predictions of simulations and semi-analytic models for $z \simeq 7$ galaxies with similar stellar masses and SFRs. △ Less

Submitted 28 November, 2023; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: 18 pages, 10 figures, 2 tables (plus 1 figure and 2 tables in the appendix). Updated to match MNRAS accepted version after minor corrections

Bright beacons? ALMA non-detection of a supposedly bright [OI] 63-um line in a redshift-6 dusty galaxy

Authors: M. Rybak, L. Lemsom, A. Lundgren, J. Zavala, J. A. Hodge, C. de Breuck, C. M. Casey, R. Decarli, K. Torstensson, J. L. Wardlow, P. P. van der Werf

Abstract: We report a non-detection of the [OI] 63-um emission line from the z = 6.03 galaxy G09.83808 using ALMA Band 9 observations, refuting the previously claimed detection with APEX by (Rybak et al. 2020); the new upper limit on the [OI] 63-um flux is almost 20-times lower. [OI] 63-um line could be a powerful tracer of neutral gas in the Epoch of Reionisation: yet our null result shows that detecting [… ▽ More We report a non-detection of the [OI] 63-um emission line from the z = 6.03 galaxy G09.83808 using ALMA Band 9 observations, refuting the previously claimed detection with APEX by (Rybak et al. 2020); the new upper limit on the [OI] 63-um flux is almost 20-times lower. [OI] 63-um line could be a powerful tracer of neutral gas in the Epoch of Reionisation: yet our null result shows that detecting [OI] 63-um from z$\geq$6 galaxies is more challenging than previously hypothesised. △ Less

Submitted 22 September, 2023; originally announced September 2023.

Journal ref: RNAAS, 7, 188 (2023)

PRUSSIC II -- ALMA imaging of dense-gas tracers in SDP.81: Evidence for low mechanical heating and a sub-solar metallicity in a z=3.04 dusty galaxy

Authors: M. Rybak, J. van Marrewijk, J. A. Hodge, P. Andreani, G. Calistro Rivera, L. Graziani, J. P. McKean, S. Viti, P. P. van der Werf

Abstract: We present deep ALMA Band 3 observations of the HCN, HCO+, and HNC (4-3) emission in SDP.81, a well-studied z = 3.042 strongly lensed galaxy. These lines trace the high-density gas, which remains almost entirely unexplored in z$\geq$1 galaxies. Additionally, these dense-gas tracers are potentially powerful diagnostics of the mechanical heating of the interstellar medium. While the HCN(4-3) and HNC… ▽ More We present deep ALMA Band 3 observations of the HCN, HCO+, and HNC (4-3) emission in SDP.81, a well-studied z = 3.042 strongly lensed galaxy. These lines trace the high-density gas, which remains almost entirely unexplored in z$\geq$1 galaxies. Additionally, these dense-gas tracers are potentially powerful diagnostics of the mechanical heating of the interstellar medium. While the HCN(4-3) and HNC(4-3) lines are not detected, the HCO+(4-3) emission is clearly detected and resolved. This is the third detection of this line in a high-redshift star-forming galaxy. We find an unusually high HCO+/HCN intensity ratio of $\geq$2.2. Based on the photodissociation region modelling, the most likely explanation for the elevated HCO+/HCN ratio is that SDP.81 has low mechanical heating - less than 10% of the total energy budget - and a sub-solar metallicity, Z=0.5 Z$_\odot$. While such conditions might not be representative of the general population of high-redshift dusty galaxies, lower-than-solar metallicity might have a significant impact on gas masses inferred from CO observations. In addition, we report the detection of CO(0-1) absorption from the foreground lensing galaxy and CO(1-0) emission from a massive companion to the lensing galaxy, approximately 50 kpc to the southeast. △ Less

Submitted 25 September, 2023; v1 submitted 5 August, 2023; originally announced August 2023.

Comments: A&A accepted, corrected version

Journal ref: A&A 679, A119 (2023)

The ALMA REBELS Survey: Discovery of a massive, highly star-forming and morphologically complex ULIRG at $z =7.31$

Authors: A. P. S. Hygate, J. A. Hodge, E. da Cunha, M. Rybak, S. Schouws, H. Inami, M. Stefanon, L. Graziani, R. Schneider, P. Dayal, R. J. Bouwens, R. Smit, R. A. A. Bowler, R. Endsley, V. Gonzalez, P. A. Oesch, D. P. Stark, H. S. B. Algera, M. Aravena, L. Barrufet, A. Ferrara, Y. Fudamoto, J. H. A, I. De Looze, T. Nanayakkara , et al. (5 additional authors not shown)

Abstract: We present Atacama Large Millimeter/Submillimeter Array (ALMA) [CII] and $\sim158$ $\rmμm$ continuum observations of REBELS-25, a massive, morphologically complex ultra-luminous infrared galaxy (ULIRG; $L_{\rm IR}=1.5^{+0.8}_{-0.5}\times10^{12}$ L$_\odot$) at $z=7.31$, spectroscopically confirmed by the Reionization Era Bright Emission Line Survey (REBELS) ALMA Large Programme. REBELS-25 has a sig… ▽ More We present Atacama Large Millimeter/Submillimeter Array (ALMA) [CII] and $\sim158$ $\rmμm$ continuum observations of REBELS-25, a massive, morphologically complex ultra-luminous infrared galaxy (ULIRG; $L_{\rm IR}=1.5^{+0.8}_{-0.5}\times10^{12}$ L$_\odot$) at $z=7.31$, spectroscopically confirmed by the Reionization Era Bright Emission Line Survey (REBELS) ALMA Large Programme. REBELS-25 has a significant stellar mass of $M_{*}=8^{+4}_{-2}\times10^{9}$ M$_\odot$. From dust-continuum and ultraviolet observations, we determine a total obscured + unobscured star formation rate of SFR $=199^{+101}_{-63}$ M$_\odot$ yr$^{-1}$. This is about four times the SFR estimated from an extrapolated main-sequence. We also infer a [CII]-based molecular gas mass of $M_{\rm H_2}=5.1^{+5.1}_{-2.6}\times10^{10}$ $M_\odot$, implying a molecular gas depletion time of $ t_{\rm depl, H_2}=0.3^{+0.3}_{-0.2}$ Gyr. We observe a [CII] velocity gradient consistent with disc rotation, but given the current resolution we cannot rule out a more complex velocity structure such as a merger. The spectrum exhibits excess [CII] emission at large positive velocities ($\sim500$ km s$^{-1}$), which we interpret as either a merging companion or an outflow. In the outflow scenario, we derive a lower limit of the mass outflow rate of 200 M$_\odot$ yr$^{-1}$, which is consistent with expectations for a star formation-driven outflow. Given its large stellar mass, SFR and molecular gas reservoir $\sim700$ Myr after the Big Bang, we explore the future evolution of REBELS-25. Considering a simple, conservative model assuming an exponentially declining star formation history, constant star formation efficiency, and no additional gas inflow, we find that REBELS-25 has the potential to evolve into a galaxy consistent with the properties of high-mass quiescent galaxies recently observed at $z\sim4$. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: Accepted for publication in MNRAS. 21 pages, 8 figures

VLA Legacy Survey of Molecular Gas in Massive Star-forming Galaxies at High Redshift

Authors: Marta Frias Castillo, Jacqueline Hodge, Matus Rybak, Paul van der Werf, Ian Smail, Jack Birkin, Chian-Chou Chen, Scott Chapman, Ryley Hill, Claudia del P. Lagos, Cheng-Lin Liao, Elisabete da Cunha, Gabriela Calistro Rivera, Jianhang Chen, Eric Jimenez Andrade, Eric Murphy, Douglas Scott, Mark Swinbank, Fabian Walter, Rob Ivison, Helmut Dannerbauer

Abstract: We present initial results of an ongoing survey with the Karl G. Jansky Very Large Array targeting the CO($J$ = 1-0) transition in a sample of 30 submillimeter-selected, dusty star-forming galaxies at $z =$ 2-5 with existing mid--$J$ CO detections from ALMA and NOEMA, of which 17 have been fully observed. We detect CO(1-0) emission in 11 targets, along with three tentative ($\sim$1.5-2$σ$) detecti… ▽ More We present initial results of an ongoing survey with the Karl G. Jansky Very Large Array targeting the CO($J$ = 1-0) transition in a sample of 30 submillimeter-selected, dusty star-forming galaxies at $z =$ 2-5 with existing mid--$J$ CO detections from ALMA and NOEMA, of which 17 have been fully observed. We detect CO(1-0) emission in 11 targets, along with three tentative ($\sim$1.5-2$σ$) detections; three galaxies are undetected. Our results yield total molecular gas masses of 6-23$\times$10$^{10}$ ($α_\mathrm{CO}$/1) M$_\odot$, with gas mass fractions, $f_\mathrm{gas}$=$M_\mathrm{mol}$/($M_*$+$M_\mathrm{mol}$), of 0.1-0.8 and a median depletion time of (140$\pm$70) Myr. We find median CO excitation ratios of $r_{31}$ = 0.75$\pm$0.39 and $r_{41}$ = 0.63$\pm$0.44, with a significant scatter. We find no significant correlation between the excitation ratio and a number of key parameters such as redshift, CO(1-0) line width or $Σ_\mathrm{SFR}$. We only find a tentative positive correlation between $r_{41}$ and the star-forming efficiency, but we are limited by our small sample size. Finally, we compare our results to predictions from the SHARK semi-analytical model, finding a good agreement between the molecular gas masses, depletion times and gas fractions of our sources and their SHARK counterparts. Our results highlight the heterogeneous nature of the most massive star-forming galaxies at high-redshift, and the importance of CO(1--0) observations to robustly constrain their total molecular gas content and ISM properties. △ Less

Submitted 10 February, 2023; v1 submitted 7 February, 2023; originally announced February 2023.

Comments: 10 pages, 15 figures. Accepted for publication in ApJ

The ALMA REBELS Survey: The Dust-obscured Cosmic Star Formation Rate Density at Redshift 7

Authors: Hiddo Algera, Hanae Inami, Pascal Oesch, Laura Sommovigo, Rychard Bouwens, Michael Topping, Sander Schouws, Mauro Stefanon, Daniel Stark, Manuel Aravena, Laia Barrufet, Elisabete da Cunha, Pratika Dayal, Ryan Endsley, Andrea Ferrara, Yoshinobu Fudamoto, Valentino Gonzalez, Luca Graziani, Jacqueline Hodge, Alexander Hygate, Ilse de Looze, Themiya Nanayakkara, Rafaella Schneider, Paul van der Werf

Abstract: Cosmic dust is an essential component shaping both the evolution of galaxies and their observational signatures. How quickly dust builds up in the early Universe remains an open question that requires deep observations at (sub-)millimeter wavelengths to resolve. Here we use Atacama Large Millimeter Array observations of 45 galaxies from the Reionization Era Bright Emission Line Survey (REBELS) and… ▽ More Cosmic dust is an essential component shaping both the evolution of galaxies and their observational signatures. How quickly dust builds up in the early Universe remains an open question that requires deep observations at (sub-)millimeter wavelengths to resolve. Here we use Atacama Large Millimeter Array observations of 45 galaxies from the Reionization Era Bright Emission Line Survey (REBELS) and its pilot programs, designed to target [CII] and dust emission in UV-selected galaxies at $z\sim7$, to investigate the dust content of high-redshift galaxies through a stacking analysis. We find that the typical fraction of obscured star formation $f_\mathrm{obs} = \mathrm{SFR}_\mathrm{IR} / \mathrm{SFR}_\mathrm{UV + IR}$ depends on stellar mass, similar to what is observed at lower redshift, and ranges from $f_\mathrm{obs} \approx 0.3 - 0.6$ for galaxies with $\log_{10}\left(M_\star / M_\odot\right) = 9.4 - 10.4$. We further adopt the $z\sim7$ stellar mass function from the literature to extract the obscured cosmic star formation rate density (SFRD) from the REBELS survey. Our results suggest only a modest decrease in the SFRD between $3\lesssim z \lesssim 7$, with dust-obscured star formation still contributing $\sim30\%$ at $z\sim7$. While we extensively discuss potential caveats, our analysis highlights the continued importance of dust-obscured star formation even well into the epoch of reionization. △ Less

Submitted 9 November, 2022; v1 submitted 17 August, 2022; originally announced August 2022.

Comments: MNRAS Accepted (01 Nov 2022)

PRUSSIC I - a JVLA survey of HCN/HCO+/HNC (1-0) emission in z$\sim$3 dusty galaxies: Low dense-gas fractions in high-redshift star-forming galaxies

Authors: M. Rybak, J. A. Hodge, T. R. Greve, D. Riechers, I. Lamperti, J. van Marrewijk, F. Walter, J. Wagg, P. P. van der Werf

Abstract: Dusty star-forming galaxies (DSFGs) at redshift z$\geq$1 are among the most vigorously star-forming galaxies in the Universe. However, their dense ($\geq$10$^5$ cm$^{-3}$ ) gas phase - typically traced by HCN(1-0) - remains almost entirely unexplored: only two DSFGs have been detected in HCN(1-0) to date. We present results of a JVLA survey of the J=1-0 transition of HCN, HCO+, and HNC(1-0) in six… ▽ More Dusty star-forming galaxies (DSFGs) at redshift z$\geq$1 are among the most vigorously star-forming galaxies in the Universe. However, their dense ($\geq$10$^5$ cm$^{-3}$ ) gas phase - typically traced by HCN(1-0) - remains almost entirely unexplored: only two DSFGs have been detected in HCN(1-0) to date. We present results of a JVLA survey of the J=1-0 transition of HCN, HCO+, and HNC(1-0) in six strongly lensed DSFGs at z = 2.5 - 3.3, effectively doubling the number of DSFGs with deep observations of these lines. We detect HCN(1-0) emission in one source (J1202+5354, 4.4$σ$), with a tentative HCO+ (1-0) detection in another (J1609+6045, 3.3$σ$). Spectral stacking yields strict upper limits on the HCN/FIR ($\leq$3.6$\times$10$^{-4}$) and HCN/CO(1-0) ratios ($\leq$0.045). The inferred HCN/FIR ratios (a proxy for the star-formation efficiency) are consistent with those in z$\sim$0 FIR-luminous starbursts. However, the HCN/CO ratios - a proxy for the dense-gas fraction - are a factor of a few lower than suggested by the two previous DSFG detections. Our results imply that most DSFGs have low dense-gas fractions. A comparison with Krumholz & Thompson (2007) models of star-forming galaxies indicates that the bulk of gas in DSFGs is at lower densities ($\approx$10$^2$ cm$^{-3}$ ), similar to "normal" star-forming galaxies, rather than ultraluminous starbursts. △ Less

Submitted 23 August, 2022; v1 submitted 14 July, 2022; originally announced July 2022.

Comments: Submitted to A&A 28th April 2022, accepted 19th July 2022

Journal ref: A&A 667, A70 (2022)

Dual constraints with ALMA: new [O III] 88 ${\rm μ}$m and dust-continuum observations reveal the ISM conditions of luminous LBGs at $z \sim 7$

Authors: Joris Witstok, Renske Smit, Roberto Maiolino, Nimisha Kumari, Manuel Aravena, Leindert Boogaard, Rychard Bouwens, Stefano Carniani, Jacqueline A. Hodge, Gareth C. Jones, Mauro Stefanon, Paul van der Werf, Sander Schouws

Abstract: We present new [O III] 88 ${\rm μ}$m observations of five bright $z \sim 7$ Lyman-break galaxies spectroscopically confirmed by ALMA through the [C II] 158 ${\rm μ}$m line, unlike recent [O III] detections where Lyman-${\rm α}$ was used. This nearly doubles the sample of Epoch of Reionisation galaxies with robust ($5 σ$) detections of [C II] and [O III]. We perform a multi-wavelength comparison wi… ▽ More We present new [O III] 88 ${\rm μ}$m observations of five bright $z \sim 7$ Lyman-break galaxies spectroscopically confirmed by ALMA through the [C II] 158 ${\rm μ}$m line, unlike recent [O III] detections where Lyman-${\rm α}$ was used. This nearly doubles the sample of Epoch of Reionisation galaxies with robust ($5 σ$) detections of [C II] and [O III]. We perform a multi-wavelength comparison with new deep HST images of the rest-frame UV, whose compact morphology aligns well with [O III] tracing ionised gas. By contrast, we find more spatially extended [C II] emission likely produced in neutral gas, as indicated by a [N II] 205 ${\rm μ}$m non-detection in one source. We find a positive correlation between the equivalent width of the optical [O III] and H${\rm β}$ lines and the [O III]/[C II] ratio, as seen in local metal-poor dwarf galaxies. Cloudy models of a nebula of typical density harbouring a young stellar population with a high ionisation parameter appear to adequately reproduce the far-infrared lines. Surprisingly, however, our models fail to reproduce the strength of [O III] 88 ${\rm μ}$m, unless we assume an ${\rm α}$/Fe enhancement and a near-solar nebular oxygen abundance. On spatially resolved scales, we find [O III]/[C II] shows a tentative anti-correlation with infrared excess, $L_{\rm IR}/L_{\rm UV}$, also seen on global scales in the local Universe. Finally, we introduce the far-infrared spectral energy distribution fitting code MERCURIUS to show that dust-continuum measurements of one source appear to favour a low dust temperature coupled with a high dust mass. This implies a high stellar metallicity yield and may point towards the need of dust production or grain-growth mechanisms beyond supernovae. △ Less

Submitted 5 July, 2022; v1 submitted 30 June, 2022; originally announced July 2022.

Comments: 23 pages, 11 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

The ALMA REBELS Survey: Average [CII] $158\,\rm{μm}$ sizes of Star-Forming Galaxies from $z\sim 7$ to $z\sim 4$

Authors: Y. Fudamoto, R. Smit, R. A. A. Bowler, P. A. Oesch, R. Bouwens, M. Stefanon, H. Inami, R. Endsley, V. Gonzalez, S. Schouws, D. Stark, H. S. B. Algera, M. Aravena, L. Barrufet, E. da Cunha, P. Dayal, A. Ferrara, L. Graziani, J. A. Hodge, A. P. S. Hygate, A. K. Inoue, T. Nanayakkara, A. Pallottini, E. Pizzati, R. Schneider , et al. (14 additional authors not shown)

Abstract: We present the average [CII] $158\,\rm{μm}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{μm}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective ra… ▽ More We present the average [CII] $158\,\rm{μm}$ emission line sizes of UV-bright star-forming galaxies at $z\sim7$. Our results are derived from a stacking analysis of [CII] $158\,\rm{μm}$ emission lines and dust continua observed by ALMA, taking advantage of the large program Reionization Era Bright Emission Line Survey (REBELS). We find that the average [CII] emission at $z\sim7$ has an effective radius $r_e$ of $2.2\pm0.2\,\rm{kpc}$. It is $\gtrsim2\times$ larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for $z\lesssim6$ galaxies. Additionally, we compared the average [CII] size with $4<z<6$ galaxies observed by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). By analysing [CII] sizes of $4<z<6$ galaxies in two redshift bins, we find an average [CII] size of $r_{\rm e}=2.2\pm0.2\,\rm{kpc}$ and $r_{\rm e}=2.5\pm0.2\,\rm{kpc}$ for $z\sim5.5$ and $z\sim4.5$ galaxies, respectively. These measurements show that star-forming galaxies, on average, show no evolution in the size of the [CII] $158\,{\rm μm}$ emitting regions at redshift between $z\sim7$ and $z\sim4$. This finding suggest that the star-forming galaxies could be morphologically dominated by gas over a wide redshift range. △ Less

Submitted 15 June, 2022; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: 9pages, 4 figures, submitted to ApJ; minor corrections were made for typos in the author list and the title

An ultra-deep multi-band VLA survey of the faint radio sky (COSMOS-XS): New constraints on the cosmic star formation history

Authors: D. van der Vlugt, J. A. Hodge, H. S. B. Algera, I. Smail, S. K. Leslie, J. F. Radcliffe, D. A. Riechers, H. Röttgering

Abstract: We make use of ultra-deep 3 GHz Karl G. Jansky Very Large Array observations of the COSMOS field from the multi-band COSMOS-XS survey to infer radio luminosity functions (LFs) of star-forming galaxies (SFGs). Using $\sim$1300 SFGs with redshifts out to $z\sim4.6$, and fixing the faint and bright end shape of the radio LF to the local values, we find a strong redshift trend that can be fitted by pu… ▽ More We make use of ultra-deep 3 GHz Karl G. Jansky Very Large Array observations of the COSMOS field from the multi-band COSMOS-XS survey to infer radio luminosity functions (LFs) of star-forming galaxies (SFGs). Using $\sim$1300 SFGs with redshifts out to $z\sim4.6$, and fixing the faint and bright end shape of the radio LF to the local values, we find a strong redshift trend that can be fitted by pure luminosity evolution with the luminosity parameter given by $α_L \propto (3.40 \pm 0.11) - (0.48 \pm 0.06)z$. We then combine the ultra-deep COSMOS-XS data-set with the shallower VLA-COSMOS $\mathrm{3\,GHz}$ large project data-set over the wider COSMOS field in order to fit for joint density+luminosity evolution, finding evidence for significant density evolution. By comparing the radio LFs to the observed far-infrared (FIR) and ultraviolet (UV) LFs, we find evidence of a significant underestimation of the UV LF by $21.6\%\, \pm \, 14.3 \, \%$ at high redshift ($3.3\,<\,z\,<\,4.6$, integrated down to $0.03\,L^{\star}_{z=3}$). We derive the cosmic star formation rate density (SFRD) by integrating the fitted radio LFs and find that the SFRD rises up to $z\,\sim\,1.8$ and then declines more rapidly than previous radio-based estimates. A direct comparison between the radio SFRD and a recent UV-based SFRD, where we integrate both LFs down to a consistent limit ($0.038\,L^{\star}_{z=3}$), reveals that the discrepancy between the radio and UV LFs translates to a significant ($\sim$1 dex) discrepancy in the derived SFRD at $z>3$, even assuming the latest dust corrections and without accounting for optically dark sources. △ Less

Submitted 8 April, 2022; originally announced April 2022.

Comments: Submitted to ApJ; 27 pages, 13 figures, 4 tables

Kiloparsec-scale imaging of the CO(1-0)-traced cold molecular gas reservoir in a z~3.4 submillimeter galaxy

Authors: Marta Frias Castillo, Matus Rybak, Jacqueline Hodge, Paul van der Werf, Dominik A. Riechers, Daniel Vieira, Gabriela Calistro Rivera, Laura N. Martinez-Ramirez, Fabian Walter, Erwin de Blok, Desika Narayanan, Jeff Wagg

Abstract: We present a high-resolution study of the cold molecular gas as traced by CO(1-0) in the unlensed z$\sim$3.4 submillimeter galaxy SMM J13120+4242, using multi-configuration observations with the Karl G. Jansky Very Large Array (JVLA). The gas reservoir, imaged on 0.39" ($\sim$3 kpc) scales, is resolved into two components separated by $\sim$11 kpc with a total extent of 16 $\sim$3 kpc. Despite the… ▽ More We present a high-resolution study of the cold molecular gas as traced by CO(1-0) in the unlensed z$\sim$3.4 submillimeter galaxy SMM J13120+4242, using multi-configuration observations with the Karl G. Jansky Very Large Array (JVLA). The gas reservoir, imaged on 0.39" ($\sim$3 kpc) scales, is resolved into two components separated by $\sim$11 kpc with a total extent of 16 $\sim$3 kpc. Despite the large spatial extent of the reservoir, the observations show a CO(1-0) FWHM linewidth of only 267 $\pm$ 64 km s$^{-1}$. We derive a revised line luminosity of L'$_\mathrm{CO(1-0)}$ = (10 $\pm$ 3) $\times$ 10$^{10}$ K km s$^{-1}$ pc$^2$ and a molecular gas mass of M$_\mathrm{gas}$ = (13 $\pm$ 3) $\times$ 10$^{10}$ ($α_\mathrm{CO}$/1) M$_{\odot}$. Despite the presence of a velocity gradient (consistent with previous resolved CO(6-5) imaging), the CO(1-0) imaging shows evidence for significant turbulent motions which are preventing the gas from fully settling into a disk. The system likely represents a merger in an advanced stage. Although the dynamical mass is highly uncertain, we use it to place an upper limit on the CO-to-H$_2$ mass conversion factor $α_\mathrm{CO}$ of 1.4. We revisit the SED fitting, finding that this galaxy lies on the very massive end of the main sequence at z = 3.4. Based on the low gas fraction, short gas depletion time and evidence for a central AGN, we propose that SMM J13120 is in a rapid transitional phase between a merger-driven starburst and an unobscured quasar. The case of SMM J13120 highlights the how mergers may drive important physical changes in galaxies without pushing them off the main sequence. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: 19 pages, 12 figures, accepted for publication in ApJ

The ALMA REBELS Survey: Dust Continuum Detections at z > 6.5

Authors: Hanae Inami, Hiddo S. B. Algera, Sander Schouws, Laura Sommovigo, Rychard Bouwens, Renske Smit, Mauro Stefanon, Rebecca A. A. Bowler, Ryan Endsley, Andrea Ferrara, Pascal Oesch, Daniel Stark, Manuel Aravena, Laia Barrufet, Elisabete da Cunha, Pratika Dayal, Ilse De Looze, Yoshinobu Fudamoto, Valentino Gonzalez, Luca Graziani, Jacqueline A. Hodge, Alexander P. S. Hygate, Themiya Nanayakkara, Andrea Pallottini, Dominik A. Riechers , et al. (3 additional authors not shown)

Abstract: We report 18 dust continuum detections ($\geq 3.3σ$) at $\sim88{\rm μm}$ and $158{\rm μm}$ out of 49 ultraviolet(UV)-bright galaxies ($M_{\rm UV} < -21.3$ mag) at $z>6.5$, observed by the Cycle-7 ALMA Large Program, REBELS and its pilot programs. This has more than tripled the number of dust continuum detections known at $z>6.5$. Out of these 18 detections, 12 are reported for the first time as pa… ▽ More We report 18 dust continuum detections ($\geq 3.3σ$) at $\sim88{\rm μm}$ and $158{\rm μm}$ out of 49 ultraviolet(UV)-bright galaxies ($M_{\rm UV} < -21.3$ mag) at $z>6.5$, observed by the Cycle-7 ALMA Large Program, REBELS and its pilot programs. This has more than tripled the number of dust continuum detections known at $z>6.5$. Out of these 18 detections, 12 are reported for the first time as part of REBELS. In addition, 15 of the dust continuum detected galaxies also show a [CII]$_{\rm 158{\rm μm}}$ emission line, providing us with accurate redshifts. We anticipate more line emission detections from six targets (including three continuum detected targets) where observations are still ongoing. The dust continuum detected sources in our sample tend to have a redder UV spectral slope than the ones without a dust continuum detection. We estimate that all of the sources have an infrared (IR) luminosity ($L_{\rm IR}$) in a range of $3-8 \times 10^{11} L_\odot$, except for one with $L_{\rm IR} = 1.5^{+0.8}_{-0.5} \times 10^{12}\,L_{\odot}$. Their fraction of obscured star formation is significant at $\gtrsim 50\%$. Some of the dust continuum detected galaxies show spatial offsets ($\sim 0.5-1.5''$) between the rest-UV and far-IR emission peaks. These separations appear to have an increasing trend against an indicator that suggests spatially decoupled phases of obscured and unobscured star formation. REBELS offers the best available statistical constraints on obscured star formation in UV-bright, massive galaxies at $z > 6.5$. △ Less

Submitted 26 September, 2022; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: 17 pages, 9 figures, 4 tables, Accepted for publication in MNRAS

Journal ref: MNRAS, 515, 3126 (2022)

Spectroscopic confirmation of a gravitationally lensed Lyman break galaxy at z$_{[CII]}$ = 6.827 using NOEMA

Authors: S. J. Molyneux, R. Smit, D. Schaerer, R. J. Bouwens, L. Bradley, J. A. Hodge, S. N. Longmore, S. Schouws, P. van der Werf, A. Zitrin, S. Phillips

Abstract: We present the spectroscopic confirmation of the brightest known gravitationally lensed Lyman break galaxy in the Epoch of Reionisation, A1703-zD1, through the detection of [C II] at a redshift of z = 6.8269 +/- 0.0004. This source was selected behind the strong lensing cluster Abell 1703, with an intrinsic L$_{UV}$ ~ L$^*$$_{z=7}$ luminosity and a very blue Spitzer/IRAC [3.6]-[4.5] colour, implyi… ▽ More We present the spectroscopic confirmation of the brightest known gravitationally lensed Lyman break galaxy in the Epoch of Reionisation, A1703-zD1, through the detection of [C II] at a redshift of z = 6.8269 +/- 0.0004. This source was selected behind the strong lensing cluster Abell 1703, with an intrinsic L$_{UV}$ ~ L$^*$$_{z=7}$ luminosity and a very blue Spitzer/IRAC [3.6]-[4.5] colour, implying high equivalent width line emission of [O III]+H$β$. [C II] is reliably detected at 6.1$σ$ co-spatial with the rest-frame UV counterpart, showing similar spatial extent. Correcting for the lensing magnification, the [C II] luminosity in A1703-zD1 is broadly consistent with the local L$_{[CII]}$ - SFR relation. We find a clear velocity gradient of 103 +/- 22 km/s across the source which possibly indicates rotation or an ongoing merger. We furthermore present spectral scans with no detected [C II] above 4.6$σ$ in two unlensed Lyman break galaxies in the EGS-CANDELS field at z ~ 6.6 - 6.9. This is the first time that NOEMA has been successfully used to observe [C II] in a 'normal' star-forming galaxy at z > 6, and our results demonstrate its capability to complement ALMA in confirming galaxies in the Epoch of Reionisation. △ Less

Submitted 28 February, 2022; originally announced March 2022.

Comments: 9 pages, 15 figures. Accepted for publication in MNRAS

The ALMA REBELS Survey: the dust content of $z \sim 7$ Lyman Break Galaxies

Authors: P. Dayal, A. Ferrara, L. Sommovigo, R. Bouwens, P. A. Oesch, R. Smit, V. Gonzalez, S. Schouws, M. Stefanon, C. Kobayashi, J. Bremer, H. S. B. Algera, M. Aravena, R. A. A. Bowler, E. da Cunha, Y. Fudamoto, L. Graziani, J. Hodge, H. Inami, I. De Looze, A. Pallottini, D. Riechers, R. Schneider, D. Stark, R. Endsley

Abstract: We include a fully coupled treatment of metal and dust enrichment into the Delphi semi-analytic model of galaxy formation to explain the dust content of 13 Lyman Break Galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at $z\simeq 7$. We find that the galaxy dust mass, $M_d$, is regulated by the combination of SNII dust production, astration, shock destructi… ▽ More We include a fully coupled treatment of metal and dust enrichment into the Delphi semi-analytic model of galaxy formation to explain the dust content of 13 Lyman Break Galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at $z\simeq 7$. We find that the galaxy dust mass, $M_d$, is regulated by the combination of SNII dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard timescale $τ_0= 30$ Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of $\sim 0.07-0.1\%$ and a UV-to-total star formation rate relation such that $log (ψ_{\rm UV}) = -0.05 ~[log (ψ)]^{2} + 0.86 ~log(ψ) -0.05$ (implying that 55-80\% of the star formation is obscured) for REBELS galaxies with stellar mass $M_* = 10^{9-10} M_\odot$. This relation reconciles the intrinsic UV luminosity of LBGs with their observed luminosity function at $z=7$. However, 2 out of the 13 systems show dust-to-stellar mass ratios ($\sim 0.94-1.1\%$) that are up to $18\times$ larger than expected from the fiducial relation. Due to the physical coupling between dust and metal enrichment, even decreasing $τ_0$ to very low values (0.3 Myr) only increases the dust-to-stellar mass ratio by a factor $ \sim 2$. Given that grain growth is not a viable explanation for such high observed ratios of the dust-to-stellar mass, we propose alternative solutions. △ Less

Submitted 22 February, 2022; originally announced February 2022.

Comments: Accepted to MNRAS

The ALMA REBELS Survey. Epoch of Reionization giants: properties of dusty galaxies at $z \approx 7$

Authors: A. Ferrara, L. Sommovigo, P. Dayal, A. Pallottini, R. J. Bouwens, V. Gonzalez, H. Inami, R. Smit, R. A. A. Bowler, R. Endsley, P. Oesch, S. Schouws, D. Stark, M. Stefanon, M. Aravena, E. da Cunha, I. De Looze, Y. Fudamoto, L. Graziani, J. Hodge, D. Riechers, R. Schneider, H. S. B. Algera, L. Barrufet, A. P. S. Hygate , et al. (5 additional authors not shown)

Abstract: We analyse FIR dust continuum measurements for 14 galaxies ($z\approx 7$) in the ALMA REBELS LP to derive their physical properties. Our model uses three input data: (a) the UV spectral slope, $β$, (b) the observed UV continuum flux at $1500$A, $F_{\rm UV}$, (c) the observed continuum flux at $\approx 158μ$m, $F_{158}$, and considers Milky Way (MW) and SMC extinction curves, along with different d… ▽ More We analyse FIR dust continuum measurements for 14 galaxies ($z\approx 7$) in the ALMA REBELS LP to derive their physical properties. Our model uses three input data: (a) the UV spectral slope, $β$, (b) the observed UV continuum flux at $1500$A, $F_{\rm UV}$, (c) the observed continuum flux at $\approx 158μ$m, $F_{158}$, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have (28-90.5)% of their star formation obscured; the total (UV+IR) star formation rates are in the range $31.5 < {\rm SFR}/ (M_\odot {\rm yr}^{-1}) < 129.5$. The sample-averaged dust mass and temperature are $(1.3\pm 1.1)\times 10^7 M_\odot$ and $52 \pm 11$ K, respectively. In some galaxies dust is abundant (REBELS-14, $M'_d \approx 3.4 \times 10^7 M_\odot$), or hot (REBELS-18, $T'_d \approx 67$ K). The dust distribution is compact ($<0.3$ kpc for 70% of the galaxies). The dust yield per supernova is $0.1 \le y_d/M_\odot \le 3.3$, with 70% of the galaxies requiring $y_d < 0.25 M_\odot$. Three galaxies (REBELS-12, 14, 39) require $y_d > 1 M_\odot$. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in [CII] nicely follow the local $L_{\rm CII}-$SFR relation, and are approximately located on the Kennicutt-Schmidt relation. The sample-averaged gas depletion time is of $0.11\, y_P^{-2}$ Gyr, where $y_P$ is the ratio of the gas-to-stellar distribution radius. For some systems a solution simultaneously matching the observed ($β, F_{\rm UV}, F_{158}$) values cannot be found. This occurs when the index $I_m = (F_{158}/F_{\rm UV})/(β-β_{\rm int})$, where $β_{\rm int}$ is the intrinsic UV slope, exceeds $I_m^*\approx 1120$ for a MW curve. For these objects we argue that the FIR and UV emitting regions are not co-spatial, questioning the use of the IRX-$β$ relation. △ Less

Submitted 15 February, 2022; originally announced February 2022.

Comments: 22 pages, 7 figures, accepted by MNRAS. Comments welcome

ALMA as a Redshift Machine: Using [CII] to Efficiently Confirm Galaxies in the Epoch of Reionization

Authors: Sander Schouws, Rychard Bouwens, Renske Smit, Jacqueline Hodge, Mauro Stefanon, Joris Witstok, Juliëtte Hilhorst, Ivo Labbe, Hiddo Algera, Leindert Boogaard, Michael Maseda, Pascal Oesch, Huub Röttgering, Paul van der Werf

Abstract: The [CII]$_{158μm}$ line has long been proposed as a promising line to spectroscopically confirm galaxies in the epoch of reionization. In this paper we present the results of new ALMA observations spectral scanning for [CII] in six particularly luminous Lyman Break Galaxies at $z\sim7$. The six sources were drawn from a sample of bright $z\sim7$ galaxies identified using the wide-area optical, ne… ▽ More The [CII]$_{158μm}$ line has long been proposed as a promising line to spectroscopically confirm galaxies in the epoch of reionization. In this paper we present the results of new ALMA observations spectral scanning for [CII] in six particularly luminous Lyman Break Galaxies at $z\sim7$. The six sources were drawn from a sample of bright $z\sim7$ galaxies identified using the wide-area optical, near-IR, and Spitzer/IRAC data over the COSMOS/UltraVISTA field and were targeted on the basis of tight constraints on their redshifts from their IRAC [3.6]-[4.5] colors. We detect significant ($>9σ$) [CII] lines in three of our six targets ($50\%$) co-spatial with the rest-$UV$ emission from the ground/space-based near-IR imaging. The luminosities of the [CII] lines lie in the range $5.6$ to $8.8\times10^{8}L_{\odot}$, consistent with the local [CII]-SFR relation. Meanwhile, their [CII]/$L_{IR}\sim1-3\times10^{-3}$ ratios are slightly elevated compared to local (U)LIRGS. This could be due to lower dust-to-gas or dust-to-metal ratios. We also find that our sources display a large kinematic diversity, with one source showing signs of rotation, one source a likely major merger and one dispersion dominated source that might contain a bright star-forming clump. Our results highlight the effectiveness of spectral scans with ALMA in spectroscopically confirming luminous galaxies in the epoch of reionization, something that is being be applied on a significantly larger sample in the on-going REBELS large program. △ Less

Submitted 8 February, 2022; originally announced February 2022.

Comments: 16 pages, 7 figures, 3 tables. Submitted to ApJ

The REBELS ALMA Survey: cosmic dust temperature evolution out to z $\sim$ 7

Authors: L. Sommovigo, A. Ferrara, A. Pallottini, P. Dayal, R. J. Bouwens, R. Smit, E. da Cunha, I. De Looze, R. A. A. Bowler, J. Hodge, H. Inami, P. Oesch, R. Endsley, V. Gonzalez, S. Schouws, D. Stark, M. Stefanon, M. Aravena, L. Graziani, D. Riechers, R. Schneider, P. van der Werf, H. Algera, L. Barrufet, Y. Fudamoto , et al. (5 additional authors not shown)

Abstract: ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature $T_d$ remains mostly unconstrained due to the few available FIR continuum data at redshift $z>5$. This introduces large uncertainties in several properties of high-$z$ galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star format… ▽ More ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature $T_d$ remains mostly unconstrained due to the few available FIR continuum data at redshift $z>5$. This introduces large uncertainties in several properties of high-$z$ galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star formation. Using a new method based on simultaneous [CII] 158$μ$m line and underlying dust continuum measurements, we derive $T_ d$ in the continuum and [CII] detected $z\approx 7$ galaxies in the ALMA Large Project REBELS sample. We find $39\ \mathrm{K} < T_d < 58\ \mathrm{K}$, and dust masses in the narrow range $M_d = (0.9-3.6)\times 10^7 M_{\odot}$. These results allow us to extend for the first time the reported $T_d(z)$ relation into the Epoch of Reionization. We produce a new physical model that explains the increasing $T_ d(z)$ trend with the decrease of gas depletion time, $t_{dep}=M_g/\mathrm{SFR}$, induced by the higher cosmological accretion rate at early times; this hypothesis yields $T_d \propto (1+z)^{0.4}$. The model also explains the observed $T_d$ scatter at a fixed redshift. We find that dust is warmer in obscured sources, as a larger obscuration results in more efficient dust heating. For UV-transparent (obscured) galaxies, $T_d$ only depends on the gas column density (metallicity), $T_d \propto N_H^{1/6}$ ($T_d \propto Z^{-1/6}$). REBELS galaxies are on average relatively transparent, with effective gas column densities around $N_H \simeq (0.03-1)\times 10^{21} \mathrm{cm}^{-2}$. We predict that other high-$z$ galaxies (e.g. MACS0416-Y1, A2744-YD4), with estimated $T_d \gg 60$ K, are significantly obscured, low-metallicity systems. In fact $T_d$ is higher in metal-poor systems due to their smaller dust content, which for fixed $L_{ IR}$ results in warmer temperatures. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: Accepted for publication in MNRAS

Probing star formation and ISM properties using galaxy disk inclination III: Evolution in dust opacity and clumpiness between redshift 0.0 < z < 0.7 constrained from UV to NIR

Authors: S. A. van der Giessen, S. K. Leslie, B. Groves, J. A. Hodge, C. C. Popescu, M. T. Sargent, E. Schinnerer, R. J. Tuffs

Abstract: (Abridged) In this paper, we use the Tuffs et al. attenuation - inclination models in ultraviolet (UV), optical, and near-infrared (NIR) bands to investigate the average global dust properties in galaxies as a function of stellar mass $M_{*}$, stellar mass surface density $μ_{*}$, star-formation rate $SFR$, specific star-formation rate $sSFR$, star-formation main-sequence offset $dMS$, and star-fo… ▽ More (Abridged) In this paper, we use the Tuffs et al. attenuation - inclination models in ultraviolet (UV), optical, and near-infrared (NIR) bands to investigate the average global dust properties in galaxies as a function of stellar mass $M_{*}$, stellar mass surface density $μ_{*}$, star-formation rate $SFR$, specific star-formation rate $sSFR$, star-formation main-sequence offset $dMS$, and star-formation rate surface density $Σ_{SFR}$ at redshifts $z \sim 0$ and $z \sim 0.7$. We use star-forming galaxies from SDSS ($\sim$ 20000) and GAMA ($\sim$ 2000) to form our low-z sample at $0.04 < z < 0.1$ and star-forming galaxies from COSMOS ($\sim$ 2000) for the sample at $0.6 <z < 0.8$. We find that galaxies at $z \sim 0.7$ have higher optical depth $τ_{B}^{f}$ and clumpiness $F$ than galaxies at $z \sim 0$. The increase in $F$ hints that the stars of $z \sim 0.7$ galaxies are less likely to escape their birth cloud, which might indicate that the birth clouds are larger. We also found that $τ_{B}^{f}$ increases with $M_{*}$ and $μ_{*}$independent of sample and therefore redshift. We found no clear trends in $τ_{B}^{f}$ or $F$ with $SFR$, which could imply that the dust mass distribution is independent of $SFR$. In turn, this would imply that the balance of dust formation and destruction is independent of the $SFR$. Based on an analysis of the inclination-dependence of the Balmer decrement, we find that reproducing the Balmer line emission requires not only a completely optically thick dust component associated with star forming regions, as in the standard Tuffs et al. model, but an extra component of optically thin dust within the birth clouds. This new component implies the existence of dust inside HII regions that attenuates the Balmer emission before it escapes through gaps in the birth cloud and we find it is more important in high-mass galaxies. △ Less

Submitted 13 June, 2022; v1 submitted 25 January, 2022; originally announced January 2022.

Comments: 25 pages, 24 figures, published in A&A

Journal ref: A&A 662, A26 (2022)

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

COLDz: Probing Cosmic Star Formation With Radio Free-free Emission

Authors: Hiddo S. B. Algera, Jacqueline A. Hodge, Dominik A. Riechers, Sarah K. Leslie, Ian Smail, Manuel Aravena, Elisabete da Cunha, Emanuele Daddi, Roberto Decarli, Mark Dickinson, Hansung B. Gim, Lucia Guaita, Benjamin Magnelli, Eric J. Murphy, Riccardo Pavesi, Mark T. Sargent, Chelsea E. Sharon, Jeff Wagg, Fabian Walter, Min Yun

Abstract: Radio free-free emission is considered to be one of the most reliable tracers of star formation in galaxies. However, as it constitutes the faintest part of the radio spectrum -- being roughly an order of magnitude less luminous than radio synchrotron emission at the GHz frequencies typically targeted in radio surveys -- the usage of free-free emission as a star formation rate tracer has mostly re… ▽ More Radio free-free emission is considered to be one of the most reliable tracers of star formation in galaxies. However, as it constitutes the faintest part of the radio spectrum -- being roughly an order of magnitude less luminous than radio synchrotron emission at the GHz frequencies typically targeted in radio surveys -- the usage of free-free emission as a star formation rate tracer has mostly remained limited to the local Universe. Here we perform a multi-frequency radio stacking analysis using deep Karl G. Jansky Very Large Array observations at 1.4, 3, 5, 10 and 34 GHz in the COSMOS and GOODS-North fields to probe free-free emission in typical galaxies at the peak of cosmic star formation. We find that $z \sim 0.5 - 3$ star-forming galaxies exhibit radio emission at rest-frame frequencies of $\sim 65 - 90$ GHz that is $\sim 1.5 - 2\times$ fainter than would be expected from a simple combination of free-free and synchrotron emission, as in the prototypical starburst galaxy M82. We interpret this as a deficit in high-frequency synchrotron emission, while the level of free-free emission is as expected from M82. We additionally provide the first constraints on the cosmic star formation history using free-free emission at $0.5 \lesssim z \lesssim 3$, which are in good agreement with more established tracers at high redshift. In the future, deep multi-frequency radio surveys will be crucial in order to accurately determine the shape of the radio spectrum of faint star-forming galaxies, and to further establish radio free-free emission as a tracer of high-redshift star formation. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: Accepted for publication in ApJ; 22 pages + appendices, 11 figures in main text

Normal, Dust-Obscured Galaxies in the Epoch of Reionization

Authors: Y. Fudamoto, P. A. Oesch, S. Schouws, M. Stefanon, R. Smit, R. J. Bouwens, R. A. A. Bowler, R. Endsley, V. Gonzalez, H. Inami, I. Labbe, D. Stark, M. Aravena, L. Barrufet, E. da Cunha, P. Dayal, A. Ferrara, L. Graziani, J. Hodge, A. Hutter, Y. Li, I. De Looze, T. Nanayakkara, A. Pallottini, D. Riechers , et al. (4 additional authors not shown)

Abstract: Over the past decades, rest-frame ultraviolet (UV) observations have provided large samples of UV luminous galaxies at redshift (z) greater than 6, during the so-called epoch of reionization. While a few of these UV identified galaxies revealed significant dust reservoirs, very heavily dust-obscured sources at these early times have remained elusive. They are limited to a rare population of extrem… ▽ More Over the past decades, rest-frame ultraviolet (UV) observations have provided large samples of UV luminous galaxies at redshift (z) greater than 6, during the so-called epoch of reionization. While a few of these UV identified galaxies revealed significant dust reservoirs, very heavily dust-obscured sources at these early times have remained elusive. They are limited to a rare population of extreme starburst galaxies, and companions of rare quasars. These studies conclude that the contribution of dust-obscured galaxies to the cosmic star formation rate density at $z>6$ is sub-dominant. Recent ALMA and Spitzer observations have identified a more abundant, less extreme population of obscured galaxies at $z=3-6$. However, this population has not been confirmed in the reionization epoch so far. Here, we report the discovery of two dust-obscured star forming galaxies at $z=6.6813\pm0.0005$ and $z=7.3521\pm0.0005$. These objects are not detected in existing rest-frame UV data, and were only discovered through their far-infrared [CII] lines and dust continuum emission as companions to typical UV-luminous galaxies at the same redshift. The two galaxies exhibit lower infrared luminosities and star-formation rates than extreme starbursts, in line with typical star-forming galaxies at $z\sim7$. This population of heavily dust-obscured galaxies appears to contribute 10-25 per cent to the $z>6$ cosmic star formation rate density. △ Less

Submitted 21 September, 2021; originally announced September 2021.

Comments: 16 pages, 7 figures, 4 tables, published in Nature

ALMA reveals large overdensity and strong clustering of galaxies in quasar environments at $z\sim4$

Authors: Cristina García-Vergara, Matus Rybak, Jacqueline Hodge, Joseph F. Hennawi, Roberto Decarli, Jorge González-López, Fabrizio Arrigoni-Battaia, Manuel Aravena, Emanuele P. Farina

Abstract: We present an Atacama Large Millimeter/submillimeter Array (ALMA) survey of CO(4--3) line emitting galaxies in 17 quasar fields at $z\sim4$ aimed to perform the first systematic search of dusty galaxies in high$-z$ quasar environments. Our blind search of galaxies around the quasars results in 5 CO emitters with S/N$\geq5.6$ within a projected radius of $R\lesssim1.5\,h^{-1}\,$cMpc and a velocity… ▽ More We present an Atacama Large Millimeter/submillimeter Array (ALMA) survey of CO(4--3) line emitting galaxies in 17 quasar fields at $z\sim4$ aimed to perform the first systematic search of dusty galaxies in high$-z$ quasar environments. Our blind search of galaxies around the quasars results in 5 CO emitters with S/N$\geq5.6$ within a projected radius of $R\lesssim1.5\,h^{-1}\,$cMpc and a velocity range of $\rmΔ\,v=\pm1000\,$km\,s$^{-1}$ around the quasar. In blank fields, we expect to detect only 0.28 CO emitters within the same volume, implying a total overdensity of $17.6^{+11.9}_{-7.6}$ in our fields, and indicating that quasars trace massive structures in the early universe. We quantify this overdensity by measuring the small-scale clustering of CO emitters around quasars, resulting in a cross-correlation length of $r_{\rm 0,QG}=8.37^{+2.42}_{-2.04}\,h^{-1}\,$cMpc, assuming a fixed slope $γ=1.8$. This contradicts the reported mild overdensities (x1.4) of Lyman alpha emitters (LAEs) in the same fields at scales of $R\lesssim7\,h^{-1}\,$cMpc which is well described by a cross-correlation length 3 times lower than that measured for CO emitters. We discuss some possibilities to explain this discrepancy, including low star formation efficiency, and excess of dust in galaxies around quasars. Finally, we constrain, for the first time, the clustering of CO emitters at $z\sim4$, finding an auto-correlation length of $r_{\rm 0,CO}=3.14\pm1.71\,h^{-1}\,$cMpc (with $γ=1.8$). Our work together with the previous study of LAEs around quasars traces simultaneously the clustering properties of both optical and dusty galaxy populations in quasars fields, stressing the importance of multi-wavelength studies, and highlighting important questions about galaxy properties in high$-z$ dense environments. △ Less

Submitted 20 September, 2021; originally announced September 2021.

Comments: 17 pages, 8 figures, submitted to the ApJ

Reionization Era Bright Emission Line Survey: Selection and Characterization of Luminous Interstellar Medium Reservoirs in the z>6.5 Universe

Authors: R. J. Bouwens, R. Smit, S. Schouws, M. Stefanon, R. Bowler, R. Endsley, V. Gonzalez, H. Inami, D. Stark, P. Oesch, J. Hodge, M. Aravena, E. da Cunha, P. Dayal, I. de Looze, A. Ferrara, Y. Fudamoto, L. Graziani, C. Li, T. Nanayakkara, A. Pallotini, R. Schneider, L. Sommovigo, M. Topping, P. van der Werf , et al. (6 additional authors not shown)

Abstract: The Reionization Era Bright Emission Line Survey (REBELS) is a cycle-7 ALMA Large Program (LP) that is identifying and performing a first characterization of many of the most luminous star-forming galaxies known in the z>6.5 universe. REBELS is providing this probe by systematically scanning 40 of the brightest UV-selected galaxies identified over a 7-deg**2 area for bright 158-micron [CII] and 88… ▽ More The Reionization Era Bright Emission Line Survey (REBELS) is a cycle-7 ALMA Large Program (LP) that is identifying and performing a first characterization of many of the most luminous star-forming galaxies known in the z>6.5 universe. REBELS is providing this probe by systematically scanning 40 of the brightest UV-selected galaxies identified over a 7-deg**2 area for bright 158-micron [CII] and 88-micron [OIII] lines and dust-continuum emission. Selection of the 40 REBELS targets was done by combining our own and other photometric selections, each of which is subject to extensive vetting using three completely independent sets of photometry and template-fitting codes. Building on the observational strategy deployed in two pilot programs, we are increasing the number of massive interstellar medium (ISM) reservoirs known at z>6.5 by ~4-5x to >30. In this manuscript, we motivate the observational strategy deployed in the REBELS program and present initial results. Based on the 60.6 hours of ALMA observations taken in the first year of the program (November 2019 to January 2020), 18 highly significant >~7sigma [CII] lines have already been discovered, the bulk of which (13/18) also show >~3.3 sigma dust-continuum emission. These newly discovered lines more than triple the number of bright ISM-cooling lines known in the z>6.5 universe, such that the number of ALMA-derived redshifts at z>6.5 already rival Lya redshift discoveries. An analysis of the completeness of our search results vs. star formation rate (SFR) suggests an ~79% efficiency in scanning for [CII] when the SFR(UV+IR) is in excess of 28 M_sol/yr. These new LP results further demonstrate ALMA's efficiency as a "redshift machine", particularly in the epoch of reionization. △ Less

Submitted 2 March, 2022; v1 submitted 25 June, 2021; originally announced June 2021.

Comments: 30 pages, 4 tables, 18 figures, accepted for publication in ApJ

Measurements of the dust properties in z~1-3 SMGs with ALMA

Authors: E. da Cunha, J. A. Hodge, C. M. Casey, H. S. B. Algera, M. Kaasinen, I. Smail, F. Walter, W. N. Brandt, H. Dannerbauer, R. Decarli, B. A. Groves, K. K. Knudsen, A. M. Swinbank, A. Weiss, P. van der Werf, J. A. Zavala

Abstract: We present Atacama Large Millimetre Array (ALMA) 2mm continuum observations of a complete and unbiased sample of 99 870micron-selected sub-millimeter galaxies (SMGs) in the Extended Chandra Deep Field South (ALESS). Our observations of each SMG reach average sensitivities of 53 microJy/beam. We measure the flux densities for 70 sources, for which we obtain a typical 870micron-to-2mm flux ratio of… ▽ More We present Atacama Large Millimetre Array (ALMA) 2mm continuum observations of a complete and unbiased sample of 99 870micron-selected sub-millimeter galaxies (SMGs) in the Extended Chandra Deep Field South (ALESS). Our observations of each SMG reach average sensitivities of 53 microJy/beam. We measure the flux densities for 70 sources, for which we obtain a typical 870micron-to-2mm flux ratio of 14 +/- 5. We do not find a redshift dependence of this flux ratio, which would be expected if the dust emission properties of our SMGs were the same at all redshifts. By combining our ALMA measurements with existing Herschel/SPIRE observations, we construct a (biased) subset of 27 galaxies for which the cool dust emission is sufficiently well sampled to obtain precise constraints on their dust properties using simple isothermal models. Thanks to our new 2mm observations, the dust emissivity index is well-constrained and robust against different dust opacity assumptions. The median dust emissivity index of our SMGs is $β\simeq1.9\pm0.4$, consistent with the emissivity index of dust in the Milky Way and other local and high-redshift galaxies, as well as classical dust grain model predictions. We also find a negative correlation between the dust temperature and $β$, similar to low-redshift observational and theoretical studies. Our results indicate that $β\simeq2$ in high-redshift dusty star-forming galaxies, implying little evolution in dust grain properties between our SMGs and local dusty galaxy samples, and suggesting these high-mass and high-metallicity galaxies have dust reservoirs driven by grain growth in their ISM. △ Less

Submitted 16 June, 2021; originally announced June 2021.

Comments: 31 pages, 20 figures, accepted for publication in the Astrophysical Journal

Significant Dust-Obscured Star Formation in Luminous Lyman-break Galaxies at $z$$\sim$$7$$-$$8$

Authors: Sander Schouws, Mauro Stefanon, Rychard J. Bouwens, Renske Smit, Jacqueline A. Hodge, Ivo Labbé, Hiddo S. Algera, Leindert Boogaard, Stefano Carniani, Yoshi Fudamoto, Benne W. Holwerda, Garth D. Illingworth, Roberto Maiolino, Michael V. Maseda, Pascal A. Oesch, Paul P. van der Werf

Abstract: We make use of ALMA continuum observations of $15$ luminous Lyman-break galaxies at $z$$\sim$$7$$-$$8$ to probe their dust-obscured star-formation. These observations are sensitive enough to probe to obscured SFRs of $20$ $M_{\odot}$$/$$yr$ ($3σ$). Six of the targeted galaxies show significant ($\geq$$3$$σ$) dust continuum detections, more than doubling the number of known dust-detected galaxies a… ▽ More We make use of ALMA continuum observations of $15$ luminous Lyman-break galaxies at $z$$\sim$$7$$-$$8$ to probe their dust-obscured star-formation. These observations are sensitive enough to probe to obscured SFRs of $20$ $M_{\odot}$$/$$yr$ ($3σ$). Six of the targeted galaxies show significant ($\geq$$3$$σ$) dust continuum detections, more than doubling the number of known dust-detected galaxies at $z$$>$$6.5$. Their IR luminosities range from $2.7$$\times$$10^{11}$ $L_{\odot}$ to $1.1$$\times$$10^{12}$ $L_{\odot}$, equivalent to obscured SFRs of $20$ to $105$ $M_{\odot}$$/$$yr$. We use our results to quantify the correlation of the infrared excess IRX on the UV-continuum slope $β_{UV}$ and stellar mass. Our results are most consistent with an SMC attenuation curve for intrinsic $UV$-slopes $β_{UV,intr}$ of $-2.63$ and most consistent with an attenuation curve in-between SMC and Calzetti for $β_{UV,intr}$ slopes of $-2.23$, assuming a dust temperature $T_d$ of $50$ K. Our fiducial IRX-stellar mass results at $z$$\sim$$7$$-$$8$ are consistent with marginal evolution from $z$$\sim$$0$. We then show how both results depend on $T_d$. For our six dust-detected sources, we estimate their dust masses and find that they are consistent with dust production from SNe if the dust destruction is low ($<$$90$%). Finally we determine the contribution of dust-obscured star formation to the star formation rate density for $UV$ luminous ($<$$-$$21.5$ mag: $\gtrsim$$1.7$$L_{UV} ^*$) $z$$\sim$$7$$-$$8$ galaxies, finding that the total SFR density at $z$$\sim$$7$ and $z$$\sim$$8$ from bright galaxies is $0.18_{-0.10}^{+0.08}$ dex and $0.20_{-0.09}^{+0.05}$ dex higher, respectively, i.e. $\sim$$\frac{1}{3}$ of the star formation in $\gtrsim$$1.7$$L_{UV} ^*$ galaxies at $z$$\sim$$7$$-$$8$ is obscured by dust. △ Less

Submitted 25 May, 2021; originally announced May 2021.

Comments: 21 pages, 12 figures, 2 tables. Submitted to ApJ

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

Ultra-faint [CII] emission in a redshift = 2 gravitationally-lensed metal-poor dwarf galaxy

Authors: M. Rybak, E. da Cunha, B. Groves, J. A. Hodge, M. Aravena, M. Maseda, L. Boogaard, D. Berg, S. Charlot, R. Decarli, D. K. Erb, E. Nelson, C. Pacifici, K. B. Schmidt, F. Walter, A. van der Wel

Abstract: Extreme emission-line galaxies (EELGs) at redshift z=1-2 provide a unique view of metal-poor, starburst sources that are the likely drivers of the cosmic reionization at z$\geq6$. However, the molecular gas reservoirs of EELGs - the fuel for their intense star-formation - remain beyond the reach of current facilities. We present ALMA [CII] and PdBI CO(2-1) observations of a z=1.8, strongly lensed… ▽ More Extreme emission-line galaxies (EELGs) at redshift z=1-2 provide a unique view of metal-poor, starburst sources that are the likely drivers of the cosmic reionization at z$\geq6$. However, the molecular gas reservoirs of EELGs - the fuel for their intense star-formation - remain beyond the reach of current facilities. We present ALMA [CII] and PdBI CO(2-1) observations of a z=1.8, strongly lensed EELG SL2S 0217, a bright Lyman-$α$ emitter with a metallicity 0.05 $Z_\odot$. We obtain a tentative (3-4$σ$) detection of the [CII] line and set an upper limit on the [CII]/SFR ratio of $\leq1\times10^6$ $L_\odot$/($M_\odot$ yr$^{-1}$), based on the synthesized images and visibility-plane analysis. The CO(2-1) emission is not detected. Photoionization modelling indicates that up to 80% of the [CII] emission originates from neutral or molecular gas, although we can not rule out that the gas is fully ionized. The very faint [CII] emission is in line with both nearby metal-poor dwarfs and high-redshift Lyman-$α$ emitters, and predictions from hydrodynamical simulations. However, the [CII] line is 30$\times$ fainter than predicted by the De Looze et al. [CII]-SFR relation for local dwarfs, illustrating the danger of extrapolating locally-calibrated relations to high-redshift, metal-poor galaxies. △ Less

Submitted 28 March, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Comments: Corrected to match the published version. We thank Y. Harikane and I. Smail for their comments on the previous arXiv version

COLDz: Deep 34 GHz Continuum Observations and Free-free Emission in High-redshift Star-forming Galaxies

Authors: H. S. B. Algera, J. A. Hodge, D. Riechers, E. J. Murphy, R. Pavesi, M. Aravena, E. Daddi, R. Decarli, M. Dickinson, M. Sargent, C. E. Sharon, J. Wagg

Abstract: The high-frequency radio sky has historically remained largely unexplored due to the typical faintness of sources in this regime, and the modest survey speed compared to observations at lower frequencies. However, high-frequency radio surveys present an invaluable tracer of high-redshift star-formation, as they directly target the faint radio free-free emission. We present deep continuum observati… ▽ More The high-frequency radio sky has historically remained largely unexplored due to the typical faintness of sources in this regime, and the modest survey speed compared to observations at lower frequencies. However, high-frequency radio surveys present an invaluable tracer of high-redshift star-formation, as they directly target the faint radio free-free emission. We present deep continuum observations at 34 GHz in the COSMOS and GOODS-North fields from the Karl G. Jansky Very Large Array (VLA), as part of the COLDz survey. The deep COSMOS mosaic spans $\sim10\text{arcmin}^2$ down to $σ=1.3μ\text{Jy beam}^{-1}$, while the wider GOODS-N observations cover $\sim50\text{arcmin}^2$ to $σ=5.3μ\text{Jy beam}^{-1}$. We present the deepest 34 GHz radio number counts to date, with five and thirteen continuum detections in COSMOS and GOODS-N, respectively. Nine galaxies show 34 GHz continuum emission that is originating from star-formation, although for two sources this is likely due to thermal emission from dust. Utilizing deep ancillary radio data at 1.4, 3, 5 and 10 GHz, we decompose the spectra of the remaining seven star-forming galaxies into their synchrotron and thermal free-free components, finding typical thermal fractions and synchrotron spectral indices comparable to those observed in local star-forming galaxies. Using calibrations from the literature, we determine free-free star-formation rates (SFRs), and show that these are in agreement with SFRs from spectral energy distribution fitting and the far-infrared/radio correlation. Our observations place strong direct constraints on the high-frequency radio emission in typical galaxies at high-redshift, and provide some of the first insight in what is set to become a key area of study with future radio facilities as the Square Kilometer Array Phase 1 and next-generation VLA. △ Less

Submitted 15 December, 2020; originally announced December 2020.

Comments: Submitted to ApJ; 22 pages + appendices; 13 figures in main text

An ALMA survey of the S2CLS UDS field: Optically invisible submillimetre galaxies

Authors: Ian Smail, U. Dudzevičiūtė, S. M. Stach, O. Almaini, J. E. Birkin, S. C. Chapman, Chian-Chou Chen, J. E. Geach, B. Gullberg, J. A. Hodge, S. Ikarashi, R. J. Ivison, D. Scott, Chris Simpson, A. M. Swinbank, A. P. Thomson, F. Walter, J. L. Wardlow, P. van der Werf

Abstract: We analyse a robust sample of 30 near-infrared-faint (K>25.3, 5 sigma) submillimetre galaxies selected across a 0.96 deg^2 field, to investigate their properties and the cause of their lack of detectable optical/near-infrared emission. Our analysis exploits precise identifications based on ALMA 870um continuum imaging, combined with the very deep near-infrared imaging from the UKIDSS-UDS survey. W… ▽ More We analyse a robust sample of 30 near-infrared-faint (K>25.3, 5 sigma) submillimetre galaxies selected across a 0.96 deg^2 field, to investigate their properties and the cause of their lack of detectable optical/near-infrared emission. Our analysis exploits precise identifications based on ALMA 870um continuum imaging, combined with the very deep near-infrared imaging from the UKIDSS-UDS survey. We estimate that K>25.3 submillimetre galaxies represent 15+/-2 per cent of the total population brighter than S870=3.6mJy, with an expected surface density of ~450/deg^2 above S870>1mJy. As such they pose a source of contamination in surveys for both high-redshift "quiescent" galaxies and very-high-redshift Lyman-break galaxies. We show that these K-faint submillimetre galaxies are simply the tail of the broader submillimetre population, with comparable dust and stellar masses to K<25.3 mag submillimetre galaxies, but lying at significantly higher redshifts (z=3.44+/-0.06 versus z=2.36+/-0.11) and having higher dust attenuation (Av=5.2+/-0.3 versus Av=2.9+/-0.1). We investigate the origin of the strong dust attenuation and find indications that these K-faint galaxies have smaller dust continuum sizes than the K<25.3 galaxies, as measured by ALMA, which suggests their high attenuation is related to their compact sizes. We find a correlation of dust attenuation with star-formation rate surface density (Sigma_SFR), with the K-faint submillimetre galaxies representing the higher-Sigma_SFR and highest-Av galaxies. The concentrated, intense star-formation activity in these systems is likely to be associated with the formation of spheroids in compact galaxies at high redshifts, but as a result of their high obscuration these are completely missed in UV, optical and even near-infrared surveys. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: 9 pages, 4 figures. Submitted to MNRAS, comments welcome

The clustering of submillimeter galaxies detected with ALMA

Authors: Cristina Garcia-Vergara, Jacqueline Hodge, Joseph F. Hennawi, Axel Weiss, Julie Wardlow, Adam D. Myers, Ryan Hickox

Abstract: Previous studies measuring the clustering of submillimeter galaxies (SMGs) have based their measurements on single-dish detected sources, finding evidence for strong clustering. However, ALMA has revealed that, due to the coarse angular resolution of these instruments, single-dish sources can be comprised of multiple sources. This implies that the clustering inferred from single-dish surveys may b… ▽ More Previous studies measuring the clustering of submillimeter galaxies (SMGs) have based their measurements on single-dish detected sources, finding evidence for strong clustering. However, ALMA has revealed that, due to the coarse angular resolution of these instruments, single-dish sources can be comprised of multiple sources. This implies that the clustering inferred from single-dish surveys may be overestimated. Here, we measure the clustering of SMGs based on the ALESS survey, an ALMA follow-up of sources previously identified in the LABOCA ECDFS Submillimeter Survey (LESS). We present a method to measure the clustering of ALMA sources that have been previously identified using single-dish telescopes, based on forward modeling both the single-dish and the ALMA observations. We constrain upper limits for the median mass of halos hosting SMGs at $1<z<3$, finding $M_{\rm{halo}}\leq2.4\times10^{12}\,\rm\,M_{\odot}$ for SMGs with flux densities $S_{870}\geq4.0\,$mJy, which is at least $3.8^{+3.8}_{-2.6}$ times lower than the mass inferred based on the clustering of the LESS sources alone. This suggests that the strength of SMG clustering based on single-dish observations was overestimated and therefore SMGs might be hosted by dark matter halos less massive than has previously been estimated. By extrapolating our models down to flux densities of $S_{870}\geq1.2\,\rm\,$mJy, we find that such SMGs inhabit halos with median mass $M_{\rm{halo}}\leq3.2\times10^{11}\,\rm\,M_{\odot}$. We conclude that only the brightest ($S_{870}\gtrsim5-6\,\rm\,$mJy) SMGs would trace massive structures at $z\sim2$ and only SMGs with $S_{870}\gtrsim6\,\rm\,$mJy may be connected to massive local elliptical galaxies, quasars at intermediate redshifts and high-redshift star-forming galaxies, whereas fainter SMGs are unlikely linked to these populations. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: 20 pages, 14 figures, Accepted for publication in ApJ

A Multi-wavelength Analysis of the Faint Radio Sky (COSMOS-XS): the Nature of the Ultra-faint Radio Population

Authors: H. S. B. Algera, D. Van der Vlugt, J. A. Hodge, I. Smail, M. Novak, J. F. Radcliffe, D. A. Riechers, H. Röttgering, V. Smolčić, F. Walter

Abstract: Ultra-deep radio surveys are an invaluable probe of dust-obscured star formation, but require a clear understanding of the relative contribution from radio AGN to be used to their fullest potential. We study the composition of the $μ$Jy radio population detected in the Karl G. Jansky Very Large Array COSMOS-XS survey based on a sample of 1540 sources detected at 3 GHz over an area of… ▽ More Ultra-deep radio surveys are an invaluable probe of dust-obscured star formation, but require a clear understanding of the relative contribution from radio AGN to be used to their fullest potential. We study the composition of the $μ$Jy radio population detected in the Karl G. Jansky Very Large Array COSMOS-XS survey based on a sample of 1540 sources detected at 3 GHz over an area of $\sim350\text{arcmin}^2$. This ultra-deep survey consists of a single pointing in the well-studied COSMOS field at both 3 and 10 GHz and reaches RMS-sensitivities of $0.53$ and $0.41μ$Jy beam$^{-1}$, respectively. We find multi-wavelength counterparts for $97\%$ of radio sources, based on a combination of near-UV/optical to sub-mm data, and through a stacking analysis at optical/near-infrared wavelengths we further show that the sources lacking such counterparts are likely to be high-redshift in nature (typical $z\sim4-5$). Utilizing the multi-wavelength data over COSMOS, we identify AGN through a variety of diagnostics and find these to make up $23.2\pm1.3\%$ of our sample, with the remainder constituting uncontaminated star-forming galaxies. However, more than half of the AGN exhibit radio emission consistent with originating from star-formation, with only $8.8\pm0.8\%$ of radio sources showing a clear excess in radio luminosity. At flux densities of $\sim30μ$Jy at 3 GHz, the fraction of star-formation powered sources reaches $\sim90\%$, and this fraction is consistent with unity at even lower flux densities. Overall, our findings imply that ultra-deep radio surveys such as COSMOS-XS constitute a highly effective means of obtaining clean samples of star-formation powered radio sources. △ Less

Submitted 28 September, 2020; originally announced September 2020.

Comments: Accepted for publication in ApJ; 23 pages + appendices; 15 figures in main text. This is a companion paper to Van der Vlugt et al., on arXiv today

An Ultra-deep Multi-band VLA Survey of the Faint Radio Sky (COSMOS-XS): Source Catalog and Number Counts

Authors: D. Van der Vlugt, H. S. B. Algera, J. A. Hodge, M. Novak, J. F. Radcliffe, D. A. Riechers, H. Röttgering, V. Smolčić, F. Walter

Abstract: We present ultra-deep, matched-resolution Karl G. Jansky Very Large Array (VLA) observations at 10 and $3$ GHz in the COSMOS field: the COSMOS-XS survey. The final 10 and $3$ GHz images cover $\sim16\rm{arcmin}^{2}$ and $\sim180\rm{arcmin}^{2}$ and reach median rms values of $0.41μ\rm{Jy\,beam}^{-1}$ and $0.53μ\rm{Jy\,beam}^{-1}$, respectively. Both images have an angular resolution of… ▽ More We present ultra-deep, matched-resolution Karl G. Jansky Very Large Array (VLA) observations at 10 and $3$ GHz in the COSMOS field: the COSMOS-XS survey. The final 10 and $3$ GHz images cover $\sim16\rm{arcmin}^{2}$ and $\sim180\rm{arcmin}^{2}$ and reach median rms values of $0.41μ\rm{Jy\,beam}^{-1}$ and $0.53μ\rm{Jy\,beam}^{-1}$, respectively. Both images have an angular resolution of $\sim 2.0''$. To fully account for the spectral shape and resolution variations across the broad bands, we image all data with a multi-scale, multi-frequency synthesis algorithm. We present source catalogs for the 10 and $3$ GHz image with 91 and 1498 sources, respectively, above a peak brightness threshold of $5σ$. We present source counts with completeness corrections included that are computed via Monte Carlo simulations. Our corrected radio counts at $3$ GHz with direct detections down to $\sim2.8μ$Jy are consistent within the uncertainties with other results at 3 and 1.4 GHz, but extend to fainter flux densities than previous direct detections. The ultra-faint $3$ GHz number counts are found to exceed the counts predicted by the semi-empirical radio sky simulations developed in the framework of the SKA Simulated Skies project, consistent with previous P(D) analyses. Our measured source counts suggest a steeper luminosity function evolution for these faint star-forming sources. The semi-empirical Tiered Radio Extragalactic Continuum Simulation (T-RECS) predicts this steeper evolution and is in better agreement with our results. The $10$ GHz radio number counts also agree with the counts predicted by the T-RECS simulation within the expected variations from cosmic variance. In summary, the multi-band, matched-resolution COSMOS-XS survey in the well-studied COSMOS field provides a high-resolution view of the ultra-faint radio sky that can help guide next generation radio facilities. △ Less

Submitted 28 September, 2020; originally announced September 2020.

Comments: Resubmitted to ApJ after incorporating suggestions from the referee; 20 pages, 17 figures. This is a companion paper to Algera et al., on arXiv today

The Evolution of the Baryons Associated with Galaxies Averaged over Cosmic Time and Space

Authors: Fabian Walter, Chris Carilli, Marcel Neeleman, Roberto Decarli, Gergo Popping, Rachel S. Somerville, Manuel Aravena, Frank Bertoldi, Leindert Boogaard, Pierre Cox, Elisabete da Cunha, Benjamin Magnelli, Danail Obreschkow, Dominik Riechers, Hans-Walter Rix, Ian Smail, Axel Weiss, Roberto J. Assef, Franz Bauer, Rychard Bouwens, Thierry Contini, Paulo C. Cortes, Emanuele Daddi, Tanio Diaz-Santo, Jorge Gonzalez-Lopez , et al. (8 additional authors not shown)

Abstract: We combine the recent determination of the evolution of the cosmic density of molecular gas (H_2) using deep, volumetric surveys, with previous estimates of the cosmic density of stellar mass, star formation rate and atomic gas (HI), to constrain the evolution of baryons associated with galaxies averaged over cosmic time and space. The cosmic HI and H_2 densities are roughly equal at z~1.5. The H_… ▽ More We combine the recent determination of the evolution of the cosmic density of molecular gas (H_2) using deep, volumetric surveys, with previous estimates of the cosmic density of stellar mass, star formation rate and atomic gas (HI), to constrain the evolution of baryons associated with galaxies averaged over cosmic time and space. The cosmic HI and H_2 densities are roughly equal at z~1.5. The H_2 density then decreases by a factor 6^{+3}_{-2} to today's value, whereas the HI density stays approximately constant. The stellar mass density is increasing continuously with time and surpasses that of the total gas density (HI and H_2) at redshift z~1.5. The growth in stellar mass cannot be accounted for by the decrease in cosmic H_2 density, necessitating significant accretion of additional gas onto galaxies. With the new H_2 constraints, we postulate and put observational constraints on a two step gas accretion process: (i) a net infall of ionized gas from the intergalactic/circumgalactic medium to refuel the extended HI reservoirs, and (ii) a net inflow of HI and subsequent conversion to H_2 in the galaxy centers. Both the infall and inflow rate densities have decreased by almost an order of magnitude since z~2. Assuming that the current trends continue, the cosmic molecular gas density will further decrease by about a factor of two over the next 5 Gyr, the stellar mass will increase by approximately 10%, and cosmic star formation activity will decline steadily toward zero, as the gas infall and accretion shut down. △ Less

Submitted 23 September, 2020; originally announced September 2020.

Comments: Accepted for publication in The Astrophysical Journal

The ALMA Spectroscopic Survey in the HUDF: Multi-band constraints on line luminosity functions and the cosmic density of molecular gas

Authors: Roberto Decarli, Manuel Aravena, Leindert Boogaard, Chris Carilli, Jorge González-López, Fabian Walter, Paulo C. Cortes, Pierre Cox, Elisabete da Cunha, Emanuele Daddi, Tanio Díaz-Santos, Jacqueline A. Hodge, Hanae Inami, Marcel Neeleman, Mladen Novak, Pascal Oesch, Gergö Popping, Dominik Riechers, Ian Smail, Bade Uzgil, Paul van der Werf, Jeff Wagg, Axel Weiss

Abstract: We present a CO and atomic fine-structure line luminosity function analysis using the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). ASPECS consists of two spatially-overlapping mosaics that cover the entire ALMA 3mm and 1.2mm bands. We combine the results of a line candidate search of the 1.2mm data cube with those previously obtained from the 3mm cube. Our analysis shows that… ▽ More We present a CO and atomic fine-structure line luminosity function analysis using the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). ASPECS consists of two spatially-overlapping mosaics that cover the entire ALMA 3mm and 1.2mm bands. We combine the results of a line candidate search of the 1.2mm data cube with those previously obtained from the 3mm cube. Our analysis shows that $\sim$80% of the line flux observed at 3mm arises from CO(2-1) or CO(3-2) emitters at $z$=1-3 (`cosmic noon'). At 1.2mm, more than half of the line flux arises from intermediate-J CO transitions ($J_{\rm up}$=3-6); $\sim12$% from neutral carbon lines; and $< 1$% from singly-ionized carbon, [CII]. This implies that future [CII] intensity mapping surveys in the epoch of reionization will need to account for a highly significant CO foreground. The CO luminosity functions probed at 1.2mm show a decrease in the number density at a given line luminosity (in units of $L'$) at increasing $J_{\rm up}$ and redshift. Comparisons between the CO luminosity functions for different CO transitions at a fixed redshift reveal sub-thermal conditions on average in galaxies up to $z\sim 4$. In addition, the comparison of the CO luminosity functions for the same transition at different redshifts reveals that the evolution is not driven by excitation. The cosmic density of molecular gas in galaxies, $ρ_{\rm H2}$, shows a redshift evolution with an increase from high redshift up to $z\sim1.5$ followed by a factor $\sim 6$ drop down to the present day. This is in qualitative agreement with the evolution of the cosmic star-formation rate density, suggesting that the molecular gas depletion time is approximately constant with redshift, after averaging over the star-forming galaxy population. △ Less

Submitted 22 September, 2020; originally announced September 2020.

Comments: 23 pages, 11 figures, 5 tables. Accepted for publication in the Astrophysical Journal

An ALMA Survey of the SCUBA-2 Cosmology Legacy Survey UKIDSS/UDS Field: The Far-infrared/Radio correlation for High-redshift Dusty Star-forming Galaxies

Authors: H. S. B. Algera, I. Smail, U. Dudzevičiūtė, A. M. Swinbank, S. Stach, J. A. Hodge, A. P. Thomson, O. Almaini, V. Arumugam, A. W. Blain, G. Calistro-Rivera, S. C. Chapman, C. -C Chen, E. da Cunha, D. Farrah, S. Leslie, D. Scott, D. Van der Vlugt, J. L. Wardlow, P. Van der Werf

Abstract: We study the radio properties of 706 sub-millimeter galaxies (SMGs) selected at 870$μ$m with the Atacama Large Millimeter Array from the SCUBA-2 Cosmology Legacy Survey map of the Ultra Deep Survey field. We detect 273 SMGs at $>4σ$ in deep Karl G. Jansky Very Large Array 1.4 GHz observations, of which a subset of 45 SMGs are additionally detected in 610 MHz Giant Metre-Wave Radio Telescope imagin… ▽ More We study the radio properties of 706 sub-millimeter galaxies (SMGs) selected at 870$μ$m with the Atacama Large Millimeter Array from the SCUBA-2 Cosmology Legacy Survey map of the Ultra Deep Survey field. We detect 273 SMGs at $>4σ$ in deep Karl G. Jansky Very Large Array 1.4 GHz observations, of which a subset of 45 SMGs are additionally detected in 610 MHz Giant Metre-Wave Radio Telescope imaging. We quantify the far-infrared/radio correlation through parameter $q_\text{IR}$, defined as the logarithmic ratio of the far-infrared and radio luminosity, and include the radio-undetected SMGs through a stacking analysis. We determine a median $q_\text{IR} = 2.20\pm0.03$ for the full sample, independent of redshift, which places these $z\sim2.5$ dusty star-forming galaxies $0.44\pm0.04$ dex below the local correlation for both normal star-forming galaxies and local ultra-luminous infrared galaxies (ULIRGs). Both the lack of redshift-evolution and the offset from the local correlation are likely the result of the different physical conditions in high-redshift starburst galaxies, compared to local star-forming sources. We explain the offset through a combination of strong magnetic fields ($B\gtrsim0.2$mG), high interstellar medium (ISM) densities and additional radio emission generated by secondary cosmic rays. While local ULIRGs are likely to have similar magnetic field strengths, we find that their compactness, in combination with a higher ISM density compared to SMGs, naturally explains why local and high-redshift dusty star-forming galaxies follow a different far-infrared/radio correlation. Overall, our findings paint SMGs as a homogeneous population of galaxies, as illustrated by their tight and non-evolving far-infrared/radio correlation. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: Accepted for publication in ApJ; 18 pages + appendices; 7 figures in main text

An ALMA/NOEMA survey of the molecular gas properties of high-redshift star-forming galaxies

Authors: Jack E. Birkin, Axel Weiss, J. L. Wardlow, Ian Smail, A. M. Swinbank, U. Dudzevičiūtė, Fang Xia An, Y. Ao, S. C. Chapman, Chian-Chou Chen, E. da Cunha, H. Dannerbauer, B. Gullberg, J. A. Hodge, S. Ikarashi, R. J. Ivison, Y. Matsuda, S. M. Stach, F. Walter, W. -H Wang, P. van der Werf

Abstract: We present a survey of the molecular gas in 61 submillimetre galaxies (SMGs) selected from 870$μ$m continuum surveys of the COSMOS, UDS and ECDFS fields, using the Atacama Large Millimeter Array (ALMA) and the Northern Extended Millimeter Array (NOEMA). 46 $^{12}$CO ($J=$2-5) emission lines are detected in 45 of the targets at $z=$1.2-4.8, with redshifts indicating that those which are submillimet… ▽ More We present a survey of the molecular gas in 61 submillimetre galaxies (SMGs) selected from 870$μ$m continuum surveys of the COSMOS, UDS and ECDFS fields, using the Atacama Large Millimeter Array (ALMA) and the Northern Extended Millimeter Array (NOEMA). 46 $^{12}$CO ($J=$2-5) emission lines are detected in 45 of the targets at $z=$1.2-4.8, with redshifts indicating that those which are submillimetre bright and undetected/faint in the optical/near-infrared typically lie at higher redshifts, with a gradient of $Δz/ΔS_{870}=$0.11$\pm$0.04mJy$^{-1}$. We also supplement our data with literature sources to construct a statistical CO spectral line energy distribution and find the $^{12}$CO line luminosities in SMGs peak at $J_{\rm up}\sim$6, consistent with the Cosmic Eyelash, among similar studies. Our SMGs lie mostly on or just above the main sequence, displaying a decrease in their gas depletion timescales $t_{\rm dep} = M_{\rm gas}/{\rm SFR}$ with redshift in the range $z\sim$1-5 and a median of 200$\pm$50Myr at $z\sim$2.8. This coincides with an increase in molecular gas fraction $μ_{\rm gas} = M_{\rm gas}/M_\ast$ across the same redshift range. Finally we demonstrate that the $M_{\rm baryon}$-$σ$ distribution of our SMGs is consistent with that followed by early-type galaxies in the Coma cluster, providing strong support to the suggestion that SMGs are progenitors of massive local spheroidal galaxies. On the basis of this we suggest that the SMG populations above and below an 870-$μ$m flux limit of $S_{870}\sim$5mJy may correspond to the division between slow- and fast-rotators seen in local early-type galaxies. △ Less

Submitted 7 September, 2020; originally announced September 2020.

Comments: Submitted to MNRAS

A Comparison of the Stellar, CO and Dust-Continuum Emission from Three, Star-Forming HUDF Galaxies at $z\sim 2$

Authors: Melanie Kaasinen, Fabian Walter, Mladen Novak, Marcel Neeleman, Ian Smail, Leindert Boogaard, Elisabete da Cunha, Axel Weiss, Daizhong Liu, Roberto Decarli, Gergö Popping, Tanio Diaz-Santos, Paulo Cortés, Manuel Aravena, Paul van der Werf, Dominik Riechers, Hanae Inami, Jacqueline A. Hodge, Hans-Walter Rix, Pierre Cox

Abstract: We compare the extent of the dust, molecular gas and stars in three star-forming galaxies, at $z= 1.4, 1.6$ and $2.7$, selected from the Hubble Ultra Deep Field based on their bright CO and dust-continuum emission as well as their large rest-frame optical sizes. The galaxies have high stellar masses, $\mathrm{M}_*>10^{11}\mathrm{M}_\odot$, and reside on, or slightly below, the main sequence of sta… ▽ More We compare the extent of the dust, molecular gas and stars in three star-forming galaxies, at $z= 1.4, 1.6$ and $2.7$, selected from the Hubble Ultra Deep Field based on their bright CO and dust-continuum emission as well as their large rest-frame optical sizes. The galaxies have high stellar masses, $\mathrm{M}_*>10^{11}\mathrm{M}_\odot$, and reside on, or slightly below, the main sequence of star-forming galaxies at their respective redshifts. We probe the dust and molecular gas using subarcsecond Atacama Large Millimeter/submillimeter Array observations of the 1.3 mm continuum and CO line emission, respectively, and probe the stellar distribution using \emph{Hubble Space Telescope} observations at 1.6 \textmu m. We find that for all three galaxies the CO emission appears $\gtrsim 30\%$ more compact than the stellar emission. For the $z= 1.4$ and $2.7$ galaxies, the dust emission is also more compact, by $\gtrsim 50\%$, than the stellar emission, whereas for the $z=1.6$ galaxy, the dust and stellar emission have similar spatial extents. This similar spatial extent is consistent with observations of local disk galaxies. However, most high redshift observations show more compact dust emission, likely due to the ubiquity of central starbursts at high redshift and the limited sensitivity of many of these observations. Using the CO emission line, we also investigate the kinematics of the cold interstellar medium in the galaxies, and find that all three have kinematics consistent with a rotation-dominated disk. △ Less

Submitted 8 September, 2020; v1 submitted 7 July, 2020; originally announced July 2020.

VLA-ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (VLASPECS): Total Cold Gas Masses and CO Line Ratios for z=2-3 "Main Sequence" Galaxies

Authors: Dominik A. Riechers, Leindert A. Boogaard, Roberto Decarli, Jorge Gonzalez-Lopez, Ian Smail, Fabian Walter, Manuel Aravena, Christopher L. Carilli, Paulo C. Cortes, Pierre Cox, Tanio Diaz-Santos, Jacqueline A. Hodge, Hanae Inami, Rob J. Ivison, Melanie Kaasinen, Jeff Wagg, Axel Weiss, Paul van der Werf

Abstract: Using the NSF's Karl G. Jansky Very Large Array (VLA), we report six detections of CO(J=1-0) emission and one upper limit in z=2-3 galaxies originally detected in higher-J CO emission in the Atacama Large submillimeter/Millimeter Array (ALMA) Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). From the CO(J=1-0) line strengths, we measure total cold molecular gas masses of M_gas = 2.4-11… ▽ More Using the NSF's Karl G. Jansky Very Large Array (VLA), we report six detections of CO(J=1-0) emission and one upper limit in z=2-3 galaxies originally detected in higher-J CO emission in the Atacama Large submillimeter/Millimeter Array (ALMA) Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). From the CO(J=1-0) line strengths, we measure total cold molecular gas masses of M_gas = 2.4-11.6 x 10^10 (alpha_CO/3.6) Msun. We also measure a median CO(J=3-2) to CO(J=1-0) line brightness temperature ratio of r_31 = 0.84 +/- 0.26, and a CO(J=7-6) to CO(J=1-0) ratio range of r_71 <0.05 to 0.17. These results suggest that CO(J=3-2) selected galaxies may have a higher CO line excitation on average than CO(J=1-0) selected galaxies, based on the limited, currently available samples from the ASPECS and VLA CO Luminosity Density at High Redshift (COLDz) surveys. This implies that previous estimates of the cosmic density of cold gas in galaxies based on CO(J=3-2) measurements should be revised down by a factor of ~=2 on average based on assumptions regarding CO excitation alone. This correction further improves the agreement between the best currently existing constraints on the cold gas density evolution across cosmic history from line scan surveys, and the implied characteristic gas depletion times. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Comments: 9 pages, 4 figures, 1 table, to appear in ApJ Letters (accepted May 19, 2020)

COLDz: A High Space Density of Massive Dusty Starburst Galaxies ~1 Billion Years after the Big Bang

Authors: Dominik A. Riechers, Jacqueline A. Hodge, Riccardo Pavesi, Emanuele Daddi, Roberto Decarli, Rob J. Ivison, Chelsea E. Sharon, Ian Smail, Fabian Walter, Manuel Aravena, Peter L. Capak, Christopher L. Carilli, Pierre Cox, Elisabete da Cunha, Helmut Dannerbauer, Mark Dickinson, Roberto Neri, Jeff Wagg

Abstract: We report the detection of CO($J$=2$\to$1) emission from three massive dusty starburst galaxies at $z$$>$5 through molecular line scans in the NSF's Karl G. Jansky Very Large Array (VLA) CO Luminosity Density at High Redshift (COLDz) survey. Redshifts for two of the sources, HDF 850.1 ($z$=5.183) and AzTEC-3 ($z$=5.298), were previously known. We revise a previous redshift estimate for the third s… ▽ More We report the detection of CO($J$=2$\to$1) emission from three massive dusty starburst galaxies at $z$$>$5 through molecular line scans in the NSF's Karl G. Jansky Very Large Array (VLA) CO Luminosity Density at High Redshift (COLDz) survey. Redshifts for two of the sources, HDF 850.1 ($z$=5.183) and AzTEC-3 ($z$=5.298), were previously known. We revise a previous redshift estimate for the third source GN10 ($z$=5.303), which we have independently confirmed through detections of CO $J$=1$\to$0, 5$\to$4, 6$\to$5, and [CII] 158 $μ$m emission with the VLA and the NOrthern Extended Milllimeter Array (NOEMA). We find that two currently independently confirmed CO sources in COLDz are "optically dark", and that three of them are dust-obscured galaxies at $z$$>$5. Given our survey area of $\sim$60 arcmin$^2$, our results appear to imply a $\sim$6-55 times higher space density of such distant dusty systems within the first billion years after the Big Bang than previously thought. At least two of these $z$$>$5 galaxies show star-formation rate surface densities consistent with so-called "maximum" starbursts, but we find significant differences in CO excitation between them. This result may suggest that different fractions of the massive gas reservoirs are located in the dense, star-forming nuclear regions - consistent with the more extended sizes of the [CII] emission compared to the dust continuum and higher [CII]-to-far-infrared luminosity ratios in those galaxies with lower gas excitation. We thus find substantial variations in the conditions for star formation between $z$$>$5 dusty starbursts, which typically have dust temperatures $\sim$57%$\pm$25% warmer than starbursts at $z$=2-3 due to their enhanced star formation activity. △ Less

Submitted 5 May, 2020; v1 submitted 21 April, 2020; originally announced April 2020.

Comments: 30 pages, 15 figures, 8 tables, to appear in ApJ (accepted April 21, 2020) [v2: fixed one latex issue and a few typos, no content changes]

Journal ref: ApJ, 895, 81 (2020)

High-redshift star formation in the ALMA era

Authors: Jacqueline A. Hodge, Elisabete da Cunha

Abstract: The Atacama Large Millimetre/submillimetre Array (ALMA) is currently in the process of transforming our view of star-forming galaxies in the distant ($z\gtrsim1$) universe. Before ALMA, most of what we knew about dust-obscured star formation in distant galaxies was limited to the brightest submillimetre sources$-$the so-called submillimetre galaxies (SMGs)$-$and even the information on those sourc… ▽ More The Atacama Large Millimetre/submillimetre Array (ALMA) is currently in the process of transforming our view of star-forming galaxies in the distant ($z\gtrsim1$) universe. Before ALMA, most of what we knew about dust-obscured star formation in distant galaxies was limited to the brightest submillimetre sources$-$the so-called submillimetre galaxies (SMGs)$-$and even the information on those sources was sparse, with resolved (i.e., sub-galactic) observations of the obscured star formation and gas reservoirs typically restricted to the most extreme and/or strongly lensed sources. Starting with the beginning of early science operations in 2011, the last nine years of ALMA observations have ushered in a new era for studies of high-redshift star formation. With its long baselines, ALMA has allowed observations of distant dust-obscured star formation with angular resolutions comparable to$-$or even far surpassing$-$the best current optical telescopes. With its bandwidth and frequency coverage, it has provided an unprecedented look at the associated molecular and atomic gas in these distant galaxies through targeted follow-up and serendipitous detections/blind line scans. Finally, with its leap in sensitivity compared to previous (sub-)millimetre arrays, it has enabled the detection of these powerful dust/gas tracers much further down the luminosity function through both statistical studies of color/mass-selected galaxy populations and dedicated deep fields. We review the main advances ALMA has helped bring about in our understanding of the dust and gas properties of high-redshift ($z\gtrsim1$) star-forming galaxies during these first nine years of its science operations, and we highlight the interesting questions that may be answered by ALMA in the years to come. △ Less

Submitted 15 January, 2021; v1 submitted 2 April, 2020; originally announced April 2020.

Comments: Invited Review for Royal Society Open Science, Published version (+ table of contents)

Deceptively cold dust in the massive starburst galaxy GN20 at $z\sim4$

Authors: Isabella Cortzen, Georgios E. Magdis, Francesco Valentino, Emanuele Daddi, Daizhong Liu, Dimitra Rigopoulou, Mark Sargent, Dominik Riechers, Diane Cormier, Jacqueline A. Hodge, Fabian Walter, David Elbaz, Matthieu Béthermin, Thomas R. Greve, Vasily Kokorev, Sune Toft

Abstract: We present new observations, carried out with IRAM NOEMA, of the atomic neutral carbon transitions [CI](1-0) at 492 GHz and [CI](2-1) at 809 GHz of GN20, a well-studied star-bursting galaxy at $z=4.05$. The high luminosity line ratio [CI](2-1)/[CI](1-0) implies an excitation temperature of $48^{+14}_{-9}$ K, which is significantly higher than the apparent dust temperature of $T_{\rm d}=33\pm2$ K (… ▽ More We present new observations, carried out with IRAM NOEMA, of the atomic neutral carbon transitions [CI](1-0) at 492 GHz and [CI](2-1) at 809 GHz of GN20, a well-studied star-bursting galaxy at $z=4.05$. The high luminosity line ratio [CI](2-1)/[CI](1-0) implies an excitation temperature of $48^{+14}_{-9}$ K, which is significantly higher than the apparent dust temperature of $T_{\rm d}=33\pm2$ K ($β=1.9$) derived under the common assumption of an optically thin far-infrared dust emission, but fully consistent with $T_{\rm d}=52\pm5$ K of a general opacity model where the optical depth ($τ$) reaches unity at a wavelength of $λ_0=170\pm23$ $μ$m. Moreover, the general opacity solution returns a factor of $\sim 2\times$ lower dust mass and, hence, a lower molecular gas mass for a fixed gas-to-dust ratio, than with the optically thin dust model. The derived properties of GN20 thus provide an appealing solution to the puzzling discovery of starbursts appearing colder than main-sequence galaxies above $z>2.5$, in addition to a lower dust-to-stellar mass ratio that approaches the physical value predicted for starburst galaxies. △ Less

Submitted 7 February, 2020; originally announced February 2020.

Comments: Comments: 6 pages, 4 figures, accepted for publication in Astronomy and Astrophysics Letters on February 4th 2020

Full of Orions: a 200-pc mapping of the interstellar medium in the redshift-3 lensed dusty star-forming galaxy SDP.81

Authors: M. Rybak, J. A. Hodge, S. Vegetti, P. van der Werf, P. Andreani, L. Graziani, J. P. McKean

Abstract: We present a sub-kpc resolved study of the interstellar medium properties in SDP.81, a z=3.042 strongly gravitationally lensed dusty star-forming galaxy, based on high-resolution, multi-band ALMA observations of the FIR continuum, CO ladder and the [CII] line. Using a visibility-plane lens modelling code, we achieve a median source-plane resolution of ~200 pc. We use photon-dominated region (PDR)… ▽ More We present a sub-kpc resolved study of the interstellar medium properties in SDP.81, a z=3.042 strongly gravitationally lensed dusty star-forming galaxy, based on high-resolution, multi-band ALMA observations of the FIR continuum, CO ladder and the [CII] line. Using a visibility-plane lens modelling code, we achieve a median source-plane resolution of ~200 pc. We use photon-dominated region (PDR) models to infer the physical conditions - far-UV field strength, density, and PDR surface temperature - of the star-forming gas on 200-pc scales, finding a FUV field strength of ~10^3-10^4 G0, gas density of ~10^5 cm^-3 and cloud surface temperatures up to 1500 K, similar to those in the Orion Trapezium region. The [CII] emission is significantly more extended than that FIR continuum: ~50 per cent of [CII] emission arises outside the FIR-bright region. The resolved [CII]/FIR ratio varies by almost 2 dex across the source, down to ~2x10^-4 in the star-forming clumps. The observed [CII]/FIR deficit trend is consistent with thermal saturation of the C+ fine-structure level occupancy at high gas temperatures. We make the source-plane reconstructions of all emission lines and continuum data publicly available. △ Less

Submitted 6 April, 2020; v1 submitted 28 December, 2019; originally announced December 2019.

Comments: 27 pages, 22 figures. 24 MB in total. MNRAS accepted 2020 March 17. Data available on request