Showing 1–28 of 28 results for author: Hygate, A

Characterizing the contribution of dust-obscured star formation at $z \gtrsim$ 5 using 18 serendipitously identified [CII] emitters

Authors: I. F. van Leeuwen, R. J. Bouwens, P. P. van der Werf, J. A. Hodge, S. Schouws, M. Stefanon, H. S. B. Algera, M. Aravena, L. A. Boogaard, R. A . A. Bowler, E. da Cunha, P. Dayal, R. Decarli, V. Gonzalez, H. Inami, I. de Looze, L. Sommovigo, B. P. Venemans, F. Walter, L. Barrufet, A. Ferrara, L. Graziani, A. P. S. Hygate, P. Oesch, M. Palla , et al. (2 additional authors not shown)

Abstract: We present a new method to determine the star formation rate (SFR) density of the Universe at $z \gtrsim 5$ that includes the contribution of dust-obscured star formation. For this purpose, we use a [CII] (158 $μ$m) selected sample of galaxies serendipitously identified in the fields of known $z\gtrsim 4.5$ objects to characterize the fraction of obscured SFR. The advantage of a [CII] selection is… ▽ More We present a new method to determine the star formation rate (SFR) density of the Universe at $z \gtrsim 5$ that includes the contribution of dust-obscured star formation. For this purpose, we use a [CII] (158 $μ$m) selected sample of galaxies serendipitously identified in the fields of known $z\gtrsim 4.5$ objects to characterize the fraction of obscured SFR. The advantage of a [CII] selection is that our sample is SFR-selected, in contrast to a UV-selection that would be biased towards unobscured star formation. We obtain a sample of 23 [CII] emitters near star-forming (SF) galaxies and QSOs -- three of which we identify for the first time -- using previous literature and archival ALMA data. 18 of these serendipitously identified galaxies have sufficiently deep rest-UV data and are used to characterize the obscured fraction of the star formation in galaxies with SFRs $\gtrsim 30\ \text{M}_{\odot} \ \text{yr}^{-1}$. We find that [CII] emitters identified around SF galaxies have $\approx$63\% of their SFR obscured, while [CII] emitters around QSOs have $\approx$93\% of their SFR obscured. By forward modeling existing wide-area UV luminosity function (LF) determinations, we derive the intrinsic UV LF using our characterization of the obscured SFR. Integrating the intrinsic LF to $M_{UV}$ = $-$20 we find that the obscured SFRD contributes to $>3\%$ and $>10\%$ of the total SFRD at $z \sim 5$ and $z \sim 6$ based on our sample of companions galaxies near SFGs and QSOs, respectively. Our results suggest that dust obscuration is not negligible at $z\gtrsim 5$, further underlining the importance of far-IR observations of the $z\gtrsim 5$ Universe. △ Less

Submitted 17 September, 2024; originally announced September 2024.

Comments: Accepted for publication in MNRAS. 24 pages, 14 figures, 8 tables (including appendices)

Accurate Simultaneous Constraints on the Dust Mass, Temperature and Emissivity Index of a Galaxy at Redshift 7.31

Authors: Hiddo Algera, Hanae Inami, Ilse De Looze, Andrea Ferrara, Hiroyuki Hirashita, Manuel Aravena, Tom Bakx, Rychard Bouwens, Rebecca Bowler, Elisabete Da Cunha, Pratika Dayal, Yoshinobu Fudamoto, Jacqueline Hodge, Alexander Hygate, Ivana van Leeuwen, Themiya Nanayakkara, Marco Palla, Andrea Pallottini, Lucie Rowland, Renske Smit, Laura Sommovigo, Mauro Stefanon, Aswin Vijayan, Paul van der Werf

Abstract: We present new multi-frequency ALMA continuum observations of the massive [$\log_{10}(M_\star/M_\odot) = 10.3_{-0.2}^{+0.1}$], UV-luminous [$M_\mathrm{UV} = -21.7 \pm 0.2$] $z=7.31$ galaxy REBELS-25 in Bands 3, 4, 5, and 9. Combining the new observations with previously-taken data in Bands 6 and 8, we cover the dust continuum emission of the galaxy in six distinct bands -- spanning rest-frame… ▽ More We present new multi-frequency ALMA continuum observations of the massive [$\log_{10}(M_\star/M_\odot) = 10.3_{-0.2}^{+0.1}$], UV-luminous [$M_\mathrm{UV} = -21.7 \pm 0.2$] $z=7.31$ galaxy REBELS-25 in Bands 3, 4, 5, and 9. Combining the new observations with previously-taken data in Bands 6 and 8, we cover the dust continuum emission of the galaxy in six distinct bands -- spanning rest-frame $50-350\,μ$m -- enabling simultaneous constraints on its dust mass ($M_\mathrm{dust}$), temperature ($T_\mathrm{dust}$) and emissivity index ($β_\mathrm{IR}$) via modified blackbody fitting. Given a fiducial model of optically thin emission, we infer a cold dust temperature of $T_\mathrm{dust} = 32_{-6}^{+9}\,$K and a high dust mass of $\log_{10}(M_\mathrm{dust}/M_\odot) = 8.2_{-0.4}^{+0.6}$, and moderately optically thick dust does not significantly alter these estimates. If we assume dust production is solely through supernovae (SNe), the inferred dust yield would be high, $y = 0.7_{-0.4}^{+2.3}\,M_\odot$ per SN. Consequently, we argue grain growth in the interstellar medium of REBELS-25 also contributes to its dust build-up. This is supported by the steep dust emissivity index $β_\mathrm{IR} = 2.5 \pm 0.4$ we measure for REBELS-25, as well as by its high stellar mass, dense interstellar medium, and metal-rich nature. Our results suggest that constraining the dust emissivity indices of high-redshift galaxies is important not only to mitigate systematic uncertainties in their dust masses and obscured star formation rates, but also to assess if dust properties evolve across cosmic time. We present an efficient observing setup to do so with ALMA, combining observations of the peak and Rayleigh-Jeans tail of the dust emission. △ Less

Submitted 5 August, 2024; originally announced August 2024.

Comments: 14 pages main text, 6 figures + appendices; re-submitted to MNRAS after a positive referee report

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

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

Journal ref: 2024, MNRAS, 528, 2407

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

The ALMA REBELS Survey: The First Infrared Luminosity Function Measurement at $\mathbf{z \sim 7}

Authors: L. Barrufet, P. A. Oesch, R. Bouwens, H. Inami, L. Sommovigo, H. Algera, E. da Cunha, M. Aravena, P. Dayal, A. Ferrara, Y. Fudamoto, V. Gonzalez, L. Graziani, A. Hygate, I. de Looze, T. Nanayakkara, A. Pallottini, R. Schneider, M. Stefanon, M. Topping, P. van Der Werf

Abstract: We present the first observational infrared luminosity function (IRLF) measurement in the Epoch of Reionization (EoR) based on a UV-selected galaxy sample with ALMA spectroscopic observations. Our analysis is based on the ALMA large program Reionization Era Bright Emission Line Survey (REBELS), which targets 42 galaxies at $\mathrm{z=6.4-7.7}$ with [CII] 158$\micron$ line scans. 16 sources exhibit… ▽ More We present the first observational infrared luminosity function (IRLF) measurement in the Epoch of Reionization (EoR) based on a UV-selected galaxy sample with ALMA spectroscopic observations. Our analysis is based on the ALMA large program Reionization Era Bright Emission Line Survey (REBELS), which targets 42 galaxies at $\mathrm{z=6.4-7.7}$ with [CII] 158$\micron$ line scans. 16 sources exhibit a dust detection, 15 of which are also spectroscopically confirmed through the [CII] line. The IR luminosities of the sample range from $\log L_{IR}/L_\odot=11.4$ to 12.2. Using the UVLF as a proxy to derive the effective volume for each of our target sources, we derive IRLF estimates, both for detections and for the full sample including IR luminosity upper limits. The resulting IRLFs are well reproduced by a Schechter function with the characteristic luminosity of $\log L_{*}/L_\odot=11.6^{+0.2}_{-0.1}$. Our observational results are in broad agreement with the average of predicted IRLFs from simulations at $z\sim7$. Conversely, our IRLFs lie significantly below lower redshift estimates, suggesting a rapid evolution from $z\sim4$ to $z\sim7$, into the reionization epoch. The inferred obscured contribution to the cosmic star-formation rate density at $z\sim7$ amounts to $\mathrm{log(SFRD/M_{\odot}/yr/Mpc^{3}) = -2.66^{+0.17}_{-0.14} }$ which is at least $\sim$10\% of UV-based estimates. We conclude that the presence of dust is already abundant in the EoR and discuss the possibility of unveiling larger samples of dusty galaxies with future ALMA and JWST observations. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: 9 pages, 5 figures

Cold Dust and Low [OIII]/[CII] Ratios: an Evolved Star-forming Population at Redshift 7

Authors: Hiddo Algera, Hanae Inami, Laura Sommovigo, Yoshinobu Fudamoto, Raffaella Schneider, Luca Graziani, Pratika Dayal, Rychard Bouwens, Manuel Aravena, Elisabete da Cunha, Andrea Ferrara, Alexander Hygate, Ivana van Leeuwen, Ilse De Looze, Marco Palla, Andrea Pallottini, Renske Smit, Mauro Stefanon, Michael Topping, Paul van der Werf

Abstract: We present new ALMA Band 8 (rest-frame $90\,μ$m) observations of three massive ($M_\star \approx 10^{10}\,M_\odot$) galaxies at $z\approx7$ previously detected in [CII]$158\,μ$m and underlying dust continuum emission in the Reionization Era Bright Emission Line Survey (REBELS). We detect the dust continuum emission of two of our targets in Band 8 (REBELS-25 and REBELS-38), while REBELS-12 remains… ▽ More We present new ALMA Band 8 (rest-frame $90\,μ$m) observations of three massive ($M_\star \approx 10^{10}\,M_\odot$) galaxies at $z\approx7$ previously detected in [CII]$158\,μ$m and underlying dust continuum emission in the Reionization Era Bright Emission Line Survey (REBELS). We detect the dust continuum emission of two of our targets in Band 8 (REBELS-25 and REBELS-38), while REBELS-12 remains undetected. Through modified blackbody fitting we determine cold dust temperatures ($T_\mathrm{dust} \approx 30 - 35\,$K) in both of the dual-band detected targets, given a fiducial model of optically thin emission with $β= 2.0$. Their dust temperatures are lower than most $z\sim7$ galaxies in the literature, and consequently their dust masses are higher ($M_\mathrm{dust} \approx 10^{8}\,M_\odot$). Nevertheless, these large dust masses are still consistent with predictions from models of dust production in the early Universe. In addition, we target and detect [OIII]$88\,μ$m emission in both REBELS-12 and REBELS-25, and find $L_\mathrm{[OIII]} / L_\mathrm{[CII]}$ ratios of approximately unity, low compared to the $L_\mathrm{[OIII]} / L_\mathrm{[CII]} \gtrsim 2 - 10$ observed in the known $z\gtrsim6$ population thus far. We argue the lower line ratios are due to a comparatively weaker ionizing radiation field resulting from the less starbursty nature of our targets. This low burstiness supports the cold dust temperatures and below average $\mathrm{[OIII]}λ\lambda4959,5007 + \mathrm{H}β$ equivalent widths of REBELS-25 and REBELS-38, compared to the known high-redshift population. Overall, this provides evidence for the existence of a massive, dust-rich galaxy population at $z\approx7$ which has previously experienced vigorous star formation, but is currently forming stars in a steady, as opposed to bursty, manner. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: 18 pages + appendices; 8 figures in main text; Submitted to MNRAS on 27 Dec 2022

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)

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

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)

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

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

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

PHANGS-ALMA Data Processing and Pipeline

Authors: Adam K. Leroy, Annie Hughes, Daizhong Liu, Jerome Pety, Erik Rosolowsky, Toshiki Saito, Eva Schinnerer, Andreas Schruba, Antonio Usero, Christopher M. Faesi, Cinthya N. Herrera, Melanie Chevance, Alexander P. S. Hygate, Amanda A. Kepley, Eric W. Koch, Miguel Querejeta, Kazimierz Sliwa, David Will, Christine D. Wilson, Gagandeep S. Anand, Ashley Barnes, Francesco Belfiore, Ivana Beslic, Frank Bigiel, Guillermo A. Blanc , et al. (43 additional authors not shown)

Abstract: We describe the processing of the PHANGS-ALMA survey and present the PHANGS-ALMA pipeline, a public software package that processes calibrated interferometric and total power data into science-ready data products. PHANGS-ALMA is a large, high-resolution survey of CO J=2-1 emission from nearby galaxies. The observations combine ALMA's main 12-m array, the 7-m array, and total power observations and… ▽ More We describe the processing of the PHANGS-ALMA survey and present the PHANGS-ALMA pipeline, a public software package that processes calibrated interferometric and total power data into science-ready data products. PHANGS-ALMA is a large, high-resolution survey of CO J=2-1 emission from nearby galaxies. The observations combine ALMA's main 12-m array, the 7-m array, and total power observations and use mosaics of dozens to hundreds of individual pointings. We describe the processing of the u-v data, imaging and deconvolution, linear mosaicking, combining interferometer and total power data, noise estimation, masking, data product creation, and quality assurance. Our pipeline has a general design and can also be applied to VLA and ALMA observations of other spectral lines and continuum emission. We highlight our recipe for deconvolution of complex spectral line observations, which combines multiscale clean, single scale clean, and automatic mask generation in a way that appears robust and effective. We also emphasize our two-track approach to masking and data product creation. We construct one set of "broadly masked" data products, which have high completeness but significant contamination by noise, and another set of "strictly masked" data products, which have high confidence but exclude faint, low signal-to-noise emission. Our quality assurance tests, supported by simulations, demonstrate that 12-m+7-m deconvolved data recover a total flux that is significantly closer to the total power flux than the 7-m deconvolved data alone. In the appendices, we measure the stability of the ALMA total power calibration in PHANGS--ALMA and test the performance of popular short-spacing correction algorithms. △ Less

Submitted 14 April, 2021; originally announced April 2021.

Comments: Accepted for publication in the Astrophysical Journal Supplement series. 65 pages, 33 figures. Software available at https://github.com/akleroy/phangs_imaging_scripts . For a full resolution version see https://sites.google.com/view/phangs/publications

Re-analysis of the 267-GHz ALMA observations of Venus: No statistically significant detection of phosphine

Authors: I. A. G. Snellen, L. Guzman-Ramirez, M. R. Hogerheijde, A. P. S. Hygate, F. F. S. van der Tak

Abstract: Context: ALMA observations of Venus at 267 GHz have been presented in the literature that show the apparent presence of phosphine (PH3) in its atmosphere. Phosphine has currently no evident production routes on the planet's surface or in its atmosphere. Aims: The aim of this work is to assess the statistical reliability of the line detection by independent re-analysis of the ALMA data. Methods… ▽ More Context: ALMA observations of Venus at 267 GHz have been presented in the literature that show the apparent presence of phosphine (PH3) in its atmosphere. Phosphine has currently no evident production routes on the planet's surface or in its atmosphere. Aims: The aim of this work is to assess the statistical reliability of the line detection by independent re-analysis of the ALMA data. Methods: The ALMA data were reduced as in the published study, following the provided scripts. First the spectral analysis presented in the study was reproduced and assessed. Subsequently, the spectrum was statistically evaluated, including its dependence on selected ALMA baselines. Results: We find that the 12th-order polynomial fit to the spectral passband utilised in the published study leads to spurious results. Following their recipe, five other >10 sigma lines can be produced in absorption or emission within 60 km/s from the PH3 1-0 transition frequency by suppressing the surrounding noise. Our independent analysis shows a feature near the PH3 frequency at a ~2 sigma level, below the common threshold for statistical significance. Since the spectral data have a non-Gaussian distribution, we consider a feature at such level as statistically unreliable that cannot be linked to a false positive probability. Conclusions: We find that the published 267-GHz ALMA data provide no statistical evidence for phosphine in the atmosphere of Venus. △ Less

Submitted 16 November, 2020; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: Accepted for publication in A&A after minor revisions

Towards a multi-tracer timeline of star formation in the LMC -- I.\ Deriving the lifetimes of H\,{\sc i} clouds

Authors: Jacob L. Ward, Mélanie Chevance, J. M. Diederik Kruijssen, Alexander P. S. Hygate, Andreas Schruba, Steven N. Longmore

Abstract: The time-scales associated with the various stages of the star formation process remain poorly constrained. This includes the earliest phases of star formation, during which molecular clouds condense out of the atomic interstellar medium. We present the first in a series of papers with the ultimate goal of compiling the first multi-tracer timeline of star formation, through a comprehensive set of… ▽ More The time-scales associated with the various stages of the star formation process remain poorly constrained. This includes the earliest phases of star formation, during which molecular clouds condense out of the atomic interstellar medium. We present the first in a series of papers with the ultimate goal of compiling the first multi-tracer timeline of star formation, through a comprehensive set of evolutionary phases from atomic gas clouds to unembedded young stellar populations. In this paper, we present an empirical determination of the lifetime of atomic clouds using the Uncertainty Principle for Star Formation formalism, based on the de-correlation of H$α$ and H\,{\sc i} emission as a function of spatial scale. We find an atomic gas cloud lifetime of 48$\substack{+13\\-8}$\,Myr. This timescale is consistent with the predicted average atomic cloud lifetime in the LMC (based on galactic dynamics) that is dominated by the gravitational collapse of the mid-plane ISM. We also determine the overlap time-scale for which both H\,{\sc i} and H$α$ emission are present to be very short ($t_{over}<1.7$\,Myr), consistent with zero, indicating that there is a near-to-complete phase change of the gas to a molecular form in an intermediary stage between H\,{\sc i} clouds and H\,{\sc ii} regions. We utilise the time-scales derived in this work to place empirically determined limits on the time-scale of molecular cloud formation. By performing the same analysis with and without the 30 Doradus region included, we find that the most extreme star forming environment in the LMC has little effect on the measured average atomic gas cloud lifetime. By measuring the lifetime of the atomic gas clouds, we place strong constraints on the physics that drives the formation of molecular clouds and establish a solid foundation for the development of a multi-tracer timeline of star formation in the LMC. △ Less

Submitted 7 July, 2020; originally announced July 2020.

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

When Gas Dynamics Decouples from Galactic Rotation: Characterizing ISM Circulation in Disk Galaxies

Authors: Jose Utreras, Guillermo A. Blanc, Andres Escala, Sharon Meidt, Eric Emsellem, Frank Bigiel, Simon C. O. Glover, Jonathan Henshaw, Alex Hygate, J. M. Diederik Kruijssen, Erik Rosolowsky, Eva Schinnerer, Andreas Schruba

Abstract: In galactic disks, galactic rotation sets the bulk motion of gas, and its energy and momentum can be transferred toward small scales. Additionally, in the interstellar medium, random and noncircular motions arise from stellar feedback, cloud-cloud interactions, and instabilities, among other processes. Our aim is to comprehend to which extent small-scale gas dynamics is decoupled from galactic rot… ▽ More In galactic disks, galactic rotation sets the bulk motion of gas, and its energy and momentum can be transferred toward small scales. Additionally, in the interstellar medium, random and noncircular motions arise from stellar feedback, cloud-cloud interactions, and instabilities, among other processes. Our aim is to comprehend to which extent small-scale gas dynamics is decoupled from galactic rotation. We study the relative contributions of galactic rotation and local noncircular motions to the circulation of gas, $Γ$, a macroscopic measure of local rotation, defined as the line integral of the velocity field around a closed path. We measure the circulation distribution as a function of spatial scale in a set of simulated disk galaxies and we model the velocity field as the sum of galactic rotation and a Gaussian random field. The random field is parameterized by a broken power law in Fourier space, with a break at the scale $λ_c$. We define the spatial scale $λ_{\rm eq}$ at which galactic rotation and non-circular motions contribute equally to $Γ$. For our simulated galaxies, the gas dynamics at the scale of molecular clouds is usually dominated by noncircular motions, but in the center of galactic disks galactic rotation is still relevant. Our model shows that the transfer of rotation from large scales breaks at the scale $λ_c$ and this transition is necessary to reproduce the circulation distribution. We find that $λ_{\rm eq}$, and therefore the structure of the gas velocity field, is set by the local conditions of gravitational stability and stellar feedback. △ Less

Submitted 10 July, 2020; v1 submitted 2 March, 2020; originally announced March 2020.

Comments: The Astrophysical Journal, 20 pages, 16 figures

Journal ref: The Astrophysical Journal, Volume 892, 2020, Issue 2, id.94

The lifecycle of molecular clouds in nearby star-forming disc galaxies

Authors: Mélanie Chevance, J. M. Diederik Kruijssen, Alexander P. S. Hygate, Andreas Schruba, Steven N. Longmore, Brent Groves, Jonathan D. Henshaw, Cinthya N. Herrera, Annie Hughes, Sarah M. R. Jeffreson, Philipp Lang, Adam K. Leroy, Sharon E. Meidt, Jérôme Pety, Alessandro Razza, Erik Rosolowsky, Eva Schinnerer, Frank Bigiel, Guillermo A. Blanc, Eric Emsellem, Christopher M. Faesi, Simon C. O. Glover, Daniel T. Haydon, I-Ting Ho, Kathryn Kreckel , et al. (7 additional authors not shown)

Abstract: It remains a major challenge to derive a theory of cloud-scale ($\lesssim100$ pc) star formation and feedback, describing how galaxies convert gas into stars as a function of the galactic environment. Progress has been hampered by a lack of robust empirical constraints on the giant molecular cloud (GMC) lifecycle. We address this problem by systematically applying a new statistical method for meas… ▽ More It remains a major challenge to derive a theory of cloud-scale ($\lesssim100$ pc) star formation and feedback, describing how galaxies convert gas into stars as a function of the galactic environment. Progress has been hampered by a lack of robust empirical constraints on the giant molecular cloud (GMC) lifecycle. We address this problem by systematically applying a new statistical method for measuring the evolutionary timeline of the GMC lifecycle, star formation, and feedback to a sample of nine nearby disc galaxies, observed as part of the PHANGS-ALMA survey. We measure the spatially-resolved ($\sim100$ pc) CO-to-H$α$ flux ratio and find a universal de-correlation between molecular gas and young stars on GMC scales, allowing us to quantify the underlying evolutionary timeline. GMC lifetimes are short, typically 10-30 Myr, and exhibit environmental variation, between and within galaxies. At kpc-scale molecular gas surface densities $Σ_{\rm H_2}\geqslant8$M$_{\odot}$pc$^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $Σ_{\rm H_2}\leqslant8$M$_{\odot}$pc$^{-2}$ GMCs decouple from galactic dynamics and live for an internal dynamical time-scale. After a long inert phase without massive star formation traced by H$α$ (75-90% of the cloud lifetime), GMCs disperse within just 1-5 Myr once massive stars emerge. The dispersal is most likely due to early stellar feedback, causing GMCs to achieve integrated star formation efficiencies of 4-10% These results show that galactic star formation is governed by cloud-scale, environmentally-dependent, dynamical processes driving rapid evolutionary cycling. GMCs and HII regions are the fundamental units undergoing these lifecycles, with mean separations of 100-300 pc in star-forming discs. Future work should characterise the multi-scale physics and mass flows driving these lifecycles. △ Less

Submitted 8 November, 2019; originally announced November 2019.

Comments: 39 pages, 14 figures; resubmitted to MNRAS after a favourable referee report (November 6, 2019)

The headlight cloud in NGC 628: An extreme giant molecular cloud in a typical galaxy disk

Authors: Cinthya N. Herrera, Jérôme Pety, Annie Hughes, Sharon E. Meidt, Kathryn Kreckel, Miguel Querejeta, Toshiki Saito, Philipp Lang, María Jesús Jiménez-Donaire, Ismael Pessa, Diane Cormier, Antonio Usero, Kazimierz Sliwa, Christopher Faesi, Guillermo A. Blanc, Frank Bigiel, Mélanie Chevance, Daniel A. Dale, Kathryn Grasha, Simon C. O. Glover, Alexander P. S. Hygate, J. M. Diederik Kruijssen, Adam K. Leroy, Erik Rosolowsky, Eva Schinnerer , et al. (3 additional authors not shown)

Abstract: Cloud-scale surveys of molecular gas reveal the link between molecular clouds properties and star formation (SF) across a range of galactic environments. Cloud populations in galaxy disks are considered to be representative of the `normal' SF. At high resolution, however, clouds with exceptional gas properties and SF activity may also be observed in normal disk environments. In this paper, we stud… ▽ More Cloud-scale surveys of molecular gas reveal the link between molecular clouds properties and star formation (SF) across a range of galactic environments. Cloud populations in galaxy disks are considered to be representative of the `normal' SF. At high resolution, however, clouds with exceptional gas properties and SF activity may also be observed in normal disk environments. In this paper, we study the brightest cloud traced in CO emission in the disk of NGC628. The cloud is spatially coincident with an extremely bright HII region. We characterize its molecular gas properties and investigate how feedback and large-scale processes influence the properties of the molecular gas. High resolution CO ALMA observations are used to characterize its mass and dynamical state, which are compared to other clouds in NGC628. A LVG analysis is used to constrain the beam-diluted density and temperature of the molecular gas. We analyze the MUSE spectrum using Starburst99 to characterize the young stellar population associated with the HII region. The cloud is massive ($1-2\times10^7$M$_{\odot}$), with a beam-diluted density of $n_{\rm H_2}=5\times10^4$ cm$^{-3}$. It has a low virial parameter, suggesting that its CO emission may be overluminous due to heating by the HII region. A young ($2-4$ Myr), massive $3\times10^{5}$ M$_{\odot}$ stellar population is associated. We argue that the cloud is currently being destroyed by feedback from young massive stars. Due to the cloud's large mass, this phase of the cloud's evolution is long enough for the impact of feedback on the excitation of the gas to be observed. Its high mass may be related to its location at a spiral co-rotation radius, where gas experiences reduced galactic shear compared to other regions of the disk, and receives a sustained inflow of gas that can promote the cloud's mass growth. △ Less

Submitted 4 November, 2019; v1 submitted 31 October, 2019; originally announced October 2019.

Comments: A&A accepted

Journal ref: A&A 634, A121 (2020)

The Gas-Star Formation Cycle in Nearby Star-Forming Galaxies I. Assessment of Multi-scale Variations

Authors: Eva Schinnerer, Annie Hughes, Adam Leroy, Brent Groves, Guillermo A. Blanc, Kathryn Kreckel, Frank Bigiel, Melanie Chevance, Daniel Dale, Eric Emsellem, Christopher Faesi, Simon Glover, Kathryn Grasha, Jonathan Henshaw, Alexander Hygate, J. M. Diederik Kruijssen, Sharon Meidt, Jerome Pety, Miguel Querejeta, Erik Rosolowsky, Toshiki Saito, Andreas Schruba, Jiayi Sun, Dyas Utomo

Abstract: The processes regulating star formation in galaxies are thought to act across a hierarchy of spatial scales. To connect extragalactic star formation relations from global and kpc-scale measurements to recent cloud-scale resolution studies, we have developed a simple, robust method that quantifies the scale dependence of the relative spatial distributions of molecular gas and recent star formation.… ▽ More The processes regulating star formation in galaxies are thought to act across a hierarchy of spatial scales. To connect extragalactic star formation relations from global and kpc-scale measurements to recent cloud-scale resolution studies, we have developed a simple, robust method that quantifies the scale dependence of the relative spatial distributions of molecular gas and recent star formation. In this paper, we apply this method to eight galaxies with roughly 1 arcsec resolution molecular gas imaging from the PHANGS-ALMA and PAWS surveys that have matched resolution, high quality narrowband Halpha imaging. At a common scale of 140pc, our massive (log(Mstar/Msun)=9.3-10.7), normally star-forming (SFR/Msun/yr=0.3-5.9) galaxies exhibit a significant reservoir of quiescent molecular gas not associated with star formation as traced by Halpha emission. Galactic structures act as backbones for both molecular and HII region distributions. As we degrade the spatial resolution, the quiescent molecular gas disappears, with the most rapid changes occurring for resolutions up to about 0.5kpc. As the resolution becomes poorer, the morphological features become indistinct for spatial scales larger than about 1kpc. The method is a promising tool to search for relationships between the quiescent or star-forming molecular reservoir and galaxy properties, but requires a larger sample size to identify robust correlations between the star-forming molecular gas fraction and global galaxy parameters. △ Less

Submitted 23 October, 2019; originally announced October 2019.

Comments: Accepted to ApJ, 29 pages with 9 figures and 7 tables

Fast and inefficient star formation due to short-lived molecular clouds and rapid feedback

Authors: J. M. Diederik Kruijssen, Andreas Schruba, Mélanie Chevance, Steven N. Longmore, Alexander P. S. Hygate, Daniel T. Haydon, Anna F. McLeod, Julianne J. Dalcanton, Linda J. Tacconi, Ewine F. van Dishoeck

Abstract: The physics of star formation and the deposition of mass, momentum, and energy into the interstellar medium by massive stars (`feedback') are the main uncertainties in modern cosmological simulations of galaxy formation and evolution. These processes determine the properties of galaxies, but are poorly understood on the $\lesssim$100 pc scale of individual giant molecular clouds (GMCs) resolved in… ▽ More The physics of star formation and the deposition of mass, momentum, and energy into the interstellar medium by massive stars (`feedback') are the main uncertainties in modern cosmological simulations of galaxy formation and evolution. These processes determine the properties of galaxies, but are poorly understood on the $\lesssim$100 pc scale of individual giant molecular clouds (GMCs) resolved in modern galaxy formation simulations. The key question is why the timescale for depleting molecular gas through star formation in galaxies ($t_{\rm dep}\approx 2$ Gyr) exceeds the dynamical timescale of GMCs by two orders of magnitude. Either most of a GMC's mass is converted into stars over many dynamical times, or only a small fraction turns into stars before the GMC is dispersed on a dynamical timescale. Here we report our observation that molecular gas and star formation are spatially de-correlated on GMC scales in the nearby flocculent spiral galaxy NGC300, contrary to their tight correlation on galactic scales. We demonstrate that this de-correlation implies rapid evolutionary cycling between GMCs, star formation, and feedback. We apply a novel statistical method to quantify the evolutionary timeline and find that star formation is regulated by efficient stellar feedback, driving GMC dispersal on short timescales (~1.5 Myr) due to radiation and stellar winds, prior to supernova explosions. This feedback limits GMC lifetimes to about one dynamical timescale (~10 Myr), with integrated star formation efficiencies of only 2-3%. Our findings reveal that galaxies consist of building blocks undergoing vigorous, feedback-driven lifecycles, that vary with the galactic environment and collectively define how galaxies form stars. Systematic applications of this multi-scale analysis to large galaxy samples will provide key input for a predictive, bottom-up theory of galaxy formation and evolution. △ Less

Submitted 23 May, 2019; v1 submitted 21 May, 2019; originally announced May 2019.

Comments: Published in the May 23 issue of Nature; this is the authors' version before final edits; low-resolution versions of the Supplementary Videos are available as ancillary files in the arXiv source (see the link in the top right of this page), whereas the full-resolution versions are available with the Nature publication

Journal ref: Nature 569, 519 (2019)

Dense gas is not enough: environmental variations in the star formation efficiency of dense molecular gas at 100pc scales in M51

Authors: M. Querejeta, E. Schinnerer, A. Schruba, E. Murphy, S. Meidt, A. Usero, A. K. Leroy, J. Pety, F. Bigiel, M. Chevance, C. M. Faesi, M. Gallagher, S. García-Burillo, S. C. O. Glover, A. P. S. Hygate, M. J. Jiménez-Donaire, J. M. D. Kruijssen, E. Momjian, E. Rosolowsky, D. Utomo

Abstract: It remains unclear what sets the efficiency with which molecular gas transforms into stars. Here we present a new VLA map of the spiral galaxy M51 in 33GHz radio continuum, an extinction-free tracer of star formation, at 3" scales (~100pc). We combined this map with interferometric PdBI/NOEMA observations of CO(1-0) and HCN(1-0) at matched resolution for three regions in M51 (central molecular rin… ▽ More It remains unclear what sets the efficiency with which molecular gas transforms into stars. Here we present a new VLA map of the spiral galaxy M51 in 33GHz radio continuum, an extinction-free tracer of star formation, at 3" scales (~100pc). We combined this map with interferometric PdBI/NOEMA observations of CO(1-0) and HCN(1-0) at matched resolution for three regions in M51 (central molecular ring, northern and southern spiral arm segments). While our measurements roughly fall on the well-known correlation between total infrared and HCN luminosity, bridging the gap between Galactic and extragalactic observations, we find systematic offsets from that relation for different dynamical environments probed in M51, e.g. the southern arm segment is more quiescent due to low star formation efficiency (SFE) of the dense gas, despite having a high dense gas fraction. Combining our results with measurements from the literature at 100pc scales, we find that the SFE of the dense gas and the dense gas fraction anti-correlate and correlate, respectively, with the local stellar mass surface density. This is consistent with previous kpc-scale studies. In addition, we find a significant anti-correlation between the SFE and velocity dispersion of the dense gas. Finally, we confirm that a correlation also holds between star formation rate surface density and the dense gas fraction, but it is not stronger than the correlation with dense gas surface density. Our results are hard to reconcile with models relying on a universal gas density threshold for star formation and suggest that turbulence and galactic dynamics play a major role in setting how efficiently dense gas converts into stars. △ Less

Submitted 9 April, 2019; v1 submitted 27 February, 2019; originally announced February 2019.

Comments: 23 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 625, A19 (2019)

An uncertainty principle for star formation - IV. On the nature and filtering of diffuse emission

Authors: Alexander P. S. Hygate, J. M. Diederik Kruijssen, Mélanie Chevance, Andreas Schruba, Daniel T. Haydon, Steven N. Longmore

Abstract: Diffuse emission is observed in galaxies in many tracers across the electromagnetic spectrum, including tracers of star formation, such as H$α$ and ultraviolet (UV), and tracers of gas mass, such as carbon monoxide (CO) transition lines and the 21-cm line of atomic hydrogen (HI). Its treatment is key to extracting meaningful information from observations such as cloud-scale star formation rates. F… ▽ More Diffuse emission is observed in galaxies in many tracers across the electromagnetic spectrum, including tracers of star formation, such as H$α$ and ultraviolet (UV), and tracers of gas mass, such as carbon monoxide (CO) transition lines and the 21-cm line of atomic hydrogen (HI). Its treatment is key to extracting meaningful information from observations such as cloud-scale star formation rates. Finally, studying diffuse emission can reveal information about the physical processes taking place in the ISM, such as chemical transitions and the nature of stellar feedback (through the photon escape fraction). We present a physically-motivated method for decomposing astronomical images containing both diffuse emission and compact regions of interest, such as HII regions or molecular clouds, into diffuse and compact component images through filtering in Fourier space. We have previously presented a statistical method for constraining the evolutionary timeline of star formation and mean separation length between compact star forming regions with galaxy-scale observations. We demonstrate how these measurements are biased by the presence of diffuse emission in tracer maps and that by using the mean separation length as a critical length scale to separate diffuse emission from compact emission, we are able to filter out this diffuse emission, thus removing its biasing effect. Furthermore, this method provides, without the need for interferometry or ancillary spectral data, a measurement of the diffuse emission fraction in input tracer maps and decomposed diffuse and compact emission maps for further analysis. △ Less

Submitted 26 October, 2018; originally announced October 2018.

Comments: 23 pages, 20 figures and 7 tables. Submitted to MNRAS (14th of September 2018)

An uncertainty principle for star formation -- III. The characteristic emission time-scales of star formation rate tracers

Authors: Daniel T. Haydon, J. M. Diederik Kruijssen, Mélanie Chevance, Alexander P. S. Hygate, Mark R. Krumholz, Andreas Schruba, Steven N. Longmore

Abstract: We recently presented a new statistical method to constrain the physics of star formation and feedback on the cloud scale by reconstructing the underlying evolutionary timeline. However, by itself this new method only recovers the relative durations of different evolutionary phases. To enable observational applications, it therefore requires knowledge of an absolute 'reference time-scale' to conve… ▽ More We recently presented a new statistical method to constrain the physics of star formation and feedback on the cloud scale by reconstructing the underlying evolutionary timeline. However, by itself this new method only recovers the relative durations of different evolutionary phases. To enable observational applications, it therefore requires knowledge of an absolute 'reference time-scale' to convert relative time-scales into absolute values. The logical choice for this reference time-scale is the duration over which the star formation rate (SFR) tracer is visible because it can be characterised using stellar population synthesis (SPS) models. In this paper, we calibrate this reference time-scale using synthetic emission maps of several SFR tracers, generated by combining the output from a hydrodynamical disc galaxy simulation with the SPS model SLUG2. We apply our statistical method to obtain self-consistent measurements of each tracer's reference time-scale. These include H$α$ and 12 ultraviolet (UV) filters (from GALEX, Swift, and HST), which cover a wavelength range 150-350 nm. At solar metallicity, the measured reference time-scales of H$α$ are ${4.32^{+0.09}_{-0.23}}$ Myr with continuum subtraction, and 6-16 Myr without, where the time-scale increases with filter width. For the UV filters we find 17-33 Myr, nearly monotonically increasing with wavelength. The characteristic time-scale decreases towards higher metallicities, as well as to lower star formation rate surface densities, owing to stellar initial mass function sampling effects. We provide fitting functions for the reference time-scale as a function of metallicity, filter width, or wavelength, to enable observational applications of our statistical method across a wide variety of galaxies. △ Less

Submitted 24 September, 2020; v1 submitted 25 October, 2018; originally announced October 2018.

Comments: 24 pages, 18 figures, 7 tables (including Appendices); published in MNRAS

Journal ref: 2020, MNRAS, 498, 235-257

The Star Formation Efficiency per Free Fall Time in Nearby Galaxies

Authors: Dyas Utomo, Jiayi Sun, Adam K. Leroy, J. M. Diederik Kruijssen, Eva Schinnerer, Andreas Schruba, Frank Bigiel, Guillermo A. Blanc, Melanie Chevance, Eric Emsellem, Cinthya Herrera, Alexander P. S. Hygate, Kathryn Kreckel, Eve C. Ostriker, Jerome Pety, Miguel Querejeta, Erik Rosolowsky, Karin M. Sandstrom, Antonio Usero

Abstract: We estimate the star formation efficiency per gravitational free fall time, $ε_{\rm ff}$, from observations of nearby galaxies with resolution matched to the typical size of a Giant Molecular Cloud. This quantity, $ε_{\rm ff}$, is theoretically important but so far has only been measured for Milky Way clouds or inferred indirectly in a few other galaxies. Using new, high resolution CO imaging from… ▽ More We estimate the star formation efficiency per gravitational free fall time, $ε_{\rm ff}$, from observations of nearby galaxies with resolution matched to the typical size of a Giant Molecular Cloud. This quantity, $ε_{\rm ff}$, is theoretically important but so far has only been measured for Milky Way clouds or inferred indirectly in a few other galaxies. Using new, high resolution CO imaging from the PHANGS-ALMA survey, we estimate the gravitational free-fall time at 60 to 120 pc resolution, and contrast this with the local molecular gas depletion time to estimate $ε_{\rm ff}$. Assuming a constant thickness of the molecular gas layer ($H = 100$ pc) across the whole sample, the median value of $ε_{\rm ff}$ in our sample is $0.7\%$. We find a mild scale-dependence, with higher $ε_{\rm ff}$ measured at coarser resolution. Individual galaxies show different values of $ε_{\rm ff}$, with the median $ε_{\rm ff}$ ranging from $0.3\%$ to $2.6\%$. We find the highest $ε_{\rm ff}$ in our lowest mass targets, reflecting both long free-fall times and short depletion times, though we caution that both measurements are subject to biases in low mass galaxies. We estimate the key systematic uncertainties, and show the dominant uncertainty to be the estimated line-of-sight depth through the molecular gas layer and the choice of star formation tracers. △ Less

Submitted 28 June, 2018; originally announced June 2018.

Comments: 10 pages, 3 figures, and 2 tables. Accepted for publication in ApJ Letters

Journal ref: ApJL Volume 861 (2018) Number 2

Cloud-Scale Molecular Gas Properties in 15 Nearby Galaxies

Authors: Jiayi Sun, Adam K. Leroy, Andreas Schruba, Erik Rosolowsky, Annie Hughes, J. M. Diederik Kruijssen, Sharon Meidt, Eva Schinnerer, Guillermo A. Blanc, Frank Bigiel, Alberto D. Bolatto, Mélanie Chevance, Brent Groves, Cinthya N. Herrera, Alexander P. S. Hygate, Jérôme Pety, Miguel Querejeta, Antonio Usero, Dyas Utomo

Abstract: We measure the velocity dispersion, $σ$, and surface density, $Σ$, of the molecular gas in nearby galaxies from CO spectral line cubes with spatial resolution $45$-$120$ pc, matched to the size of individual giant molecular clouds. Combining $11$ galaxies from the PHANGS-ALMA survey with $4$ targets from the literature, we characterize ${\sim}30,000$ independent sightlines where CO is detected at… ▽ More We measure the velocity dispersion, $σ$, and surface density, $Σ$, of the molecular gas in nearby galaxies from CO spectral line cubes with spatial resolution $45$-$120$ pc, matched to the size of individual giant molecular clouds. Combining $11$ galaxies from the PHANGS-ALMA survey with $4$ targets from the literature, we characterize ${\sim}30,000$ independent sightlines where CO is detected at good significance. $Σ$ and $σ$ show a strong positive correlation, with the best-fit power law slope close to the expected value for resolved, self-gravitating clouds. This indicates only weak variation in the virial parameter $α_\mathrm{vir}\proptoσ^2/Σ$, which is ${\sim}1.5$-$3.0$ for most galaxies. We do, however, observe enormous variation in the internal turbulent pressure $P_\mathrm{turb}\proptoΣ\,σ^2$, which spans ${\sim}5\rm\;dex$ across our sample. We find $Σ$, $σ$, and $P_\mathrm{turb}$ to be systematically larger in more massive galaxies. The same quantities appear enhanced in the central kpc of strongly barred galaxies relative to their disks. Based on sensitive maps of M31 and M33, the slope of the $σ$-$Σ$ relation flattens at $Σ\lesssim10\rm\;M_\odot\,pc^{-2}$, leading to high $σ$ for a given $Σ$ and high apparent $α_\mathrm{vir}$. This echoes results found in the Milky Way, and likely originates from a combination of lower beam filling factors and a stronger influence of local environment on the dynamical state of molecular gas in the low density regime. △ Less

Submitted 2 May, 2018; originally announced May 2018.

Comments: Accepted for publication in ApJ. 45 pages, 11 figures, 8 tables, 4 Appendices; key results summarized in Figure 10. Machine-readable table can be downloaded at http://www.astronomy.ohio-state.edu/~sun.1608/datafile3.txt prior to publication. For a brief video describing the main results of this paper, please see https://www.youtube.com/watch?v=-_eL7t1PVq8&t

Journal ref: ApJ, 860, 172 (2018)

An uncertainty principle for star formation - II. A new method for characterising the cloud-scale physics of star formation and feedback across cosmic history

Authors: J. M. Diederik Kruijssen, Andreas Schruba, Alexander P. S. Hygate, Chia-Yu Hu, Daniel T. Haydon, Steven N. Longmore

Abstract: The cloud-scale physics of star formation and feedback represent the main uncertainty in galaxy formation studies. Progress is hampered by the limited empirical constraints outside the restricted environment of the Local Group. In particular, the poorly-quantified time evolution of the molecular cloud lifecycle, star formation, and feedback obstructs robust predictions on the scales smaller than t… ▽ More The cloud-scale physics of star formation and feedback represent the main uncertainty in galaxy formation studies. Progress is hampered by the limited empirical constraints outside the restricted environment of the Local Group. In particular, the poorly-quantified time evolution of the molecular cloud lifecycle, star formation, and feedback obstructs robust predictions on the scales smaller than the disc scale height that are resolved in modern galaxy formation simulations. We present a new statistical method to derive the evolutionary timeline of molecular clouds and star-forming regions. By quantifying the excess or deficit of the gas-to-stellar flux ratio around peaks of gas or star formation tracer emission, we directly measure the relative rarity of these peaks, which allows us to derive their lifetimes. We present a step-by-step, quantitative description of the method and demonstrate its practical application. The method's accuracy is tested in nearly 300 experiments using simulated galaxy maps, showing that it is capable of constraining the molecular cloud lifetime and feedback time-scale to $<0.1$ dex precision. Access to the evolutionary timeline provides a variety of additional physical quantities, such as the cloud-scale star formation efficiency, the feedback outflow velocity, the mass loading factor, and the feedback energy or momentum coupling efficiencies to the ambient medium. We show that the results are robust for a wide variety of gas and star formation tracers, spatial resolutions, galaxy inclinations, and galaxy sizes. Finally, we demonstrate that our method can be applied out to high redshift ($z\lesssim4$) with a feasible time investment on current large-scale observatories. This is a major shift from previous studies that constrained the physics of star formation and feedback in the immediate vicinity of the Sun. △ Less

Submitted 4 September, 2018; v1 submitted 30 April, 2018; originally announced May 2018.

Comments: 86 pages (including appendices), 30 figures, 16 tables; accepted by MNRAS (April 27, 2018). The method is summarised in Figure 2 and 3, with the main results shown in Figures 18, 26, and 27. For a cursory overview of the paper (of about 10 pages), please consult Sections 1, 2, 3.1, 4.4, 6, and 8. Figure 27 replaced with authors' original version

Journal ref: MNRAS, 479, 1866 (2018)

New and updated convex shape models of asteroids based on optical data from a large collaboration network

Authors: J. Hanuš, J. Ďurech, D. A. Oszkiewicz, R. Behrend, B. Carry, M. Delbo', O. Adam, V. Afonina, R. Anquetin, P. Antonini, L. Arnold, M. Audejean, P. Aurard, M. Bachschmidt, B. Badue, E. Barbotin, P. Barroy, P. Baudouin, L. Berard, N. Berger, L. Bernasconi, J-G. Bosch, S. Bouley, I. Bozhinova, J. Brinsfield , et al. (144 additional authors not shown)

Abstract: Asteroid modeling efforts in the last decade resulted in a comprehensive dataset of almost 400 convex shape models and their rotation states. This amount already provided a deep insight into physical properties of main-belt asteroids or large collisional families. We aim to increase the number of asteroid shape models and rotation states. Such results are an important input for various further stu… ▽ More Asteroid modeling efforts in the last decade resulted in a comprehensive dataset of almost 400 convex shape models and their rotation states. This amount already provided a deep insight into physical properties of main-belt asteroids or large collisional families. We aim to increase the number of asteroid shape models and rotation states. Such results are an important input for various further studies such as analysis of asteroid physical properties in different populations, including smaller collisional families, thermophysical modeling, and scaling shape models by disk-resolved images, or stellar occultation data. This provides, in combination with known masses, bulk density estimates, but constrains also theoretical collisional and evolutional models of the Solar System. We use all available disk-integrated optical data (i.e., classical dense-in-time photometry obtained from public databases and through a large collaboration network as well as sparse-in-time individual measurements from a few sky surveys) as an input for the convex inversion method, and derive 3D shape models of asteroids, together with their rotation periods and orientations of rotation axes. The key ingredient is the support of more that one hundred observers who submit their optical data to publicly available databases. We present updated shape models for 36 asteroids, for which mass estimates are currently available in the literature or their masses will be most likely determined from their gravitational influence on smaller bodies, which orbital deflection will be observed by the ESA Gaia astrometric mission. This was achieved by using additional optical data from recent apparitions for the shape optimization. Moreover, we also present new shape model determinations for 250 asteroids, including 13 Hungarias and 3 near-Earth asteroids. △ Less

Submitted 26 October, 2015; originally announced October 2015.

Comments: Accepted for publication in A&A

Journal ref: A&A 586, A108 (2016)