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

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$\alpha$ 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$\alpha$ 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

The ALMA REBELS Survey: dust continuum detections at z > 6.5

We report 18 dust continuum detections (≥3.3σ) at ∼88 and 158 $\mu{\rm m}$ out of 49 ultraviolet (UV)-bright galaxies (MUV 6.5, observed by the Cycle-7 Atacama Large Millimeter/submillimeter Array (ALMA) Large Program, Reionization-Era Bright Emission Line Survey (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 [C ii]$_{\rm 158\,{\rm \mu 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. We estimate that all of the sources have an infrared (IR) luminosity (LIR) in a range of $3\!-\!8 \times 10^{11}\, {\rm L_\odot }$, except for one with $L_{\rm IR} = 1.5^{+0.8}_{-0.5} \times 10^{12}\, \, {\rm L_{\odot }}$. Their fraction of obscured star formation is significant at ${\gtrsim} 50{{\ \rm per\ cent}}$, despite being UV-selected galaxies. Some of the dust continuum detected galaxies show spatial offsets (∼0.5-1.5 arcsec) between the rest-UV and far-IR emission peaks. These separations could imply spatially decoupled phases of obscured and unobscured star formation, but a higher spatial resolution observation is required to confirm this. REBELS offers the best available statistical constraints on obscured star formation in UV-luminous galaxies at $z$ > 6.5

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

We present Atacama Large Millimeter/Submillimeter Array (ALMA) [C ii] and ∼158 continuum observations of REBELS-25, a massive, morphologically complex ultra-luminous infrared galaxy (ULIRG; LIR = L⊙) 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. From dust-continuum and ultraviolet observations, we determine a total obscured + unobscured star formation rate of SFR. This is about four times the SFR estimated from an extrapolated main sequence. We also infer a [C ii]-based molecular gas mass of, implying a molecular gas depletion time of Gyr. We observe a [C ii] 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 [C ii] emission at large positive velocities (∼500 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, which is consistent with expectations for a star-formation-driven outflow. Given its large stellar mass, SFR, and molecular gas reservoir ∼700 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 ∼4

The ALMA REBELS Survey. Epoch of Reionization giants: Properties of dusty galaxies at z ≈ 7

We analyse FIR dust continuum measurements for 14 galaxies (redshift z ≈ 7) in the ALMA REBELS Large Program to derive their physical properties. Our model uses three input data, i.e. (a) the UV spectral slope, β, (b) the observed UV continuum flux at 1500 Å, F1500, (c) the observed continuum flux at ≈158μm, F158, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have 28 − 90.5 per cent of their star formation obscured; the total (UV+IR) star formation rates are in the range 31.5 1 M⊙, which is likely inconsistent with pure SN production, and might require dust growth via accretion of heavy elements from the interstellar medium. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in [C II] nicely follow the local LCII −SFR relation, and are approximately located on the Kennicutt-Schmidt relation. The sample-averaged gas depletion time is of 0.11y−2P Gyr, where yP is the ratio of the gas-to-stellar distribution radius. For some systems a solution simultaneously matching the observed (β, F1500, F158) values cannot be found. This occurs when the index Im = (F158/F1500)/(β − βint), where βint is the intrinsic UV slope, exceeds I∗m≈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

The REBELS ALMA Survey: cosmic dust temperature evolution out to z ∼ 7

ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature Td 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 [C II] 158μm line and underlying dust continuum measurements, we derive Td in the continuum and [C II] detected z ≈ 7 galaxies in the ALMA Large Project REBELS sample. We find 39 K < Td < 58 K, and dust masses in the narrow range Md = (0.9 − 3.6) × 107M⊙. These results allow us to extend for the first time the reported Td(z) relation into the Epoch of Reionization. We produce a new physical model that explains the increasing Td(z) trend with the decrease of gas depletion time, tdep = Mg/SFR, induced by the higher cosmological accretion rate at early times; this hypothesis yields Td∝(1 + z)0.4. The model also explains the observed Td 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, Td only depends on the gas column density (metallicity), Td∝N1/6H (Td∝Z−1/6). REBELS galaxies are on average relatively transparent, with effective gas column densities around NH ≃ (0.03 − 1) × 1021cm−2. We predict that other high-z galaxies (e.g. MACS0416-Y1, A2744-YD4), with estimated Td ≫ 60 K, are significantly obscured, low-metallicity systems. In fact Td is higher in metal-poor systems due to their smaller dust content, which for fixed LIR results in warmer temperatures

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

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$\alpha$ 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 $\Sigma_{\rm H_2}\geqslant8$M$_{\odot}$pc$^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $\Sigma_{\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$\alpha$ (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

The ALMA REBELS Survey: Average [C ii] 158 μm Sizes of Star-forming Galaxies from z ∼ 7 to z ∼ 4

We present the average [C ii] 158 μm emission line sizes of UV-bright star-forming galaxies at z ∼ 7. Our results are derived from a stacking analysis of [C ii] 158 μm emission lines and dust continua observed by the Atacama Large Millimeter/submillimeter Array (ALMA), taking advantage of the large program Reionization Era Bright Emission Line Survey. We find that the average [C ii] emission at z ∼ 7 has an effective radius re of 2.2 ± 0.2 kpc. It is ≳2× larger than the dust continuum and the rest-frame UV emission, in agreement with recently reported measurements for z ≲ 6 galaxies. Additionally, we compared the average [C ii] size with 4 < z < 6 galaxies observed by the ALMA Large Program to INvestigate [C ii] at Early times (ALPINE). By analyzing [C ii] sizes of 4 < z < 6 galaxies in two redshift bins, we find an average [C ii] size of re = 2.2 ± 0.2 kpc and re = 2.5 ± 0.2 kpc for z ∼ 5.5 and z ∼ 4.5 galaxies, respectively. These measurements show that star-forming galaxies, on average, show no evolution in the size of the [C ii] 158 μm emitting regions at redshift between z ∼ 7 and z ∼ 4. This finding suggests that the star-forming galaxies could be morphologically dominated by gas over a wide redshift range