-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
