Star Formation Efficiency Research Articles

Molecular Gas and the Star-Formation Process on Cloud Scales in Nearby Galaxies

Observations that resolve nearby galaxies into individual regions across multiple phases of the gas–star formation–feedback “matter cycle” have provided a sharp new view of molecular clouds, star-formation efficiencies, timescales for region evolution, and stellar feedback. We synthesize these results, covering aspects relevant to the interpretation of observables, and conclude the following: ▪The observed cloud-scale molecular gas surface density, line width, and internal pressure all reflect the large-scale galactic environment while also appearing mostly consistent with properties of a turbulent medium strongly affected by self-gravity.▪Cloud-scale data allow for statistical inference of both evolutionary and physical timescales. These suggest a period of cloud collapse on the order of the free-fall or turbulent crossing time (∼10–30 Myr) followed by forming massive stars and subsequent rapid (≲5 Myr) gas clearing after the onset of star formation. The star-formation efficiency per free-fall time is well determined over thousands of individual regions at εff ≈ 0.5−0.3 +0.7%.▪The role of stellar feedback is now measured using multiple observational approaches. The net yield is constrained by the requirement to support the vertical weight of the galaxy disk. Meanwhile, the short gas-clearing timescales suggest a large role for presupernova feedback in cloud disruption. This leaves the supernovae free to exert a large influence on the larger galaxy, including stirring turbulence, launching galactic-scale winds, and carving superbubbles.

The Star Formation Histories, Star Formation Efficiencies and Ionizing Sources of ATLASGAL Clumps with H ii Regions

1226 ATLASGAL clumps with H ii regions (H ii-clumps) were matched with radio sources in the CORNISH-North/South surveys, and 392 of them have corresponding radio sources. We determined the stellar luminosity L *,T84 according to the Lyman continuum flux N Ly. When the bolometric luminosity of H ii-clumps is less than log10(L bol,obs/L ⊙) ≈3.7, corresponding to a clump mass log10(M cl/M ⊙) ≈ 2.55, the values of L *,T84 derived from N Ly overestimate the actual stellar luminosities, because the accretion onto the protostars contributes significantly to the radio emission. After subtracting the accretion luminosity from L *,T84, we obtained reasonable estimates of the stellar luminosity. Using the 0.5 Myr isochrone, we calculated the stellar masses according to the stellar luminosities, and found that they roughly follow the mmax−Mecl relation of embedded clusters, consistent with the ionizing sources representing the most massive stars in the embedded clusters of H ii-clumps. We also studied the contribution of the possible flaring activity to the observed stellar luminosity and found that they can be neglected. We further studied the change of SFE with the clump mass. According to the derived mass of the most massive star in each HII-clump, using the theoretical mmax−Mecl relation, we calculated the mass of the corresponding embedded cluster and then the SFE of the clump. The SFE decreases with increasing clump mass, with a median value of ≈0.3. We also independently derived the SFE for each H ii-clump based on the model developed in our previous work. The SFEs of H ii-clumps derived from the observation and the model are in good agreement. Concerning the star formation histories of the ATLASGAL clumps, low-mass clumps may reach the peak of star formation earlier than high-mass clumps, consistent with the shorter free-fall time of low-mass clumps.

FEASTS Combined with Interferometry. II. Significantly Changed H i Surface Densities and Even More Inefficient Star Formation in Galaxy Outer Disks

We update the H i surface density measurements for a subset of 17 THINGS galaxies by dealing with the short-spacing problem of the original Very Large Array (VLA) H i images. It is the same sample that Bigiel et al. used to study the relation between H i surface densities and star formation (SF) rate surface densities in galaxy outer disks, which are beyond the optical radius r 25. For 10 galaxies, the update is based on combining original THINGS VLA H i images with H i images taken by the single-dish FAST in the FEASTS program. The median increment of H i surface densities in outer disks is 0.15–0.4 dex at a given new H i surface density. Several galaxies change significantly in the shape of radial profiles H i surface densities, and seven galaxies are now more than 1σ below the H i size–mass relation. We update the H i star formation laws in outer disks. The median relation between the H i surface densities and SF rate surface densities based on pixel-wise measurements shifts downward by around 0.15 dex because the H i surface density values shift rightward, and the scatter increases significantly. The scatter of the relation, indicating the star-forming efficiency, exhibits a much stronger positive correlation with the stellar mass surface density than before. Thus, detecting the previously missed, diffuse H i due to the short-spacing problem of the Very Large Array observations is important in revealing the true condition and variation in SF possibly regulated by stellar feedbacks in the localized environment of outer disks.

Suppressed Cosmic-Ray Energy Densities in Molecular Clouds from Streaming Instability-regulated Transport

Cosmic rays (CRs) are the primary driver of ionization in star-forming molecular clouds (MCs). Despite their potential impacts on gas dynamics and chemistry, no simulations of star cluster formation following the creation of individual stars have included explicit cosmic-ray transport (CRT) to date. We conduct the first numerical simulations following the collapse of a 2000M ⊙ MC and the subsequent star formation including CRT using the STAR FORmation in Gaseous Environments framework implemented in the GIZMO code. We show that when CRT is streaming-dominated, the CR energy in the cloud is strongly attenuated due to energy losses from the streaming instability. Consequently, in a Milky Way–like environment the median CR ionization rate in the cloud is low (ζ ≲ 2 × 10−19 s−1) during the main star-forming epoch of the calculation and the impact of CRs on the star formation in the cloud is limited. However, in high-CR environments, the CR distribution in the cloud is elevated (ζ ≲ 6 × 10−18), and the relatively higher CR pressure outside the cloud causes slightly earlier cloud collapse and increases the star formation efficiency by 50% to ∼13%. The initial mass function is similar in all cases except with possible variations in a high-CR environment. Further studies are needed to explain the range of ionization rates observed in MCs and explore star formation in extreme CR environments.

Playing with FIRE: A Galactic Feedback-halting Experiment Challenges Star Formation Rate Theories

Stellar feedback influences the star formation rate (SFR) and the interstellar medium of galaxies in ways that are difficult to quantify numerically, because feedback is an essential ingredient of realistic simulations. To overcome this, we conduct a feedback-halting experiment starting with a Milky Way–mass galaxy in the second-generation Feedback In Realistic Environments (FIRE-2) simulation framework. By terminating feedback, and comparing to a simulation in which feedback is maintained, we monitor how the runs diverge. We find that without feedback, the interstellar turbulent velocities decay. There is a marked increase of dense material, while the SFR increases by over an order of magnitude. Importantly, this SFR boost is a factor of ∼15–20 larger than is accounted for by the increased freefall rate caused by higher densities. This implies that feedback moderates the star formation efficiency per freefall time more directly than simply through the density distribution. To probe changes at the scale of giant molecular clouds (GMCs), we identify GMCs using density and virial parameter thresholds, tracking clouds as the galaxy evolves. Halting feedback stimulates rapid changes, including a proliferation of new bound clouds, a decrease of turbulent support in loosely bound clouds, an overall increase in cloud densities, and a surge of internal star formation. Computing the cloud-integrated SFR using several theories of turbulence regulation, we show that these theories underpredict the surge in SFR by at least a factor of 3. We conclude that galactic star formation is essentially feedback regulated on scales that include GMCs, and that stellar feedback affects GMCs in multiple ways.

The properties and kinematics of HCN emission across the closest starburst galaxy NGC 253 observed with ALMA

Context. Investigating molecular gas tracers, such as hydrogen cyanide (HCN), to probe higher densities than CO emission across nearby galaxies remains challenging. This is due to the large observing times required to detect HCN at a high sensitivity and spatial resolution. Although approximate kiloparsec scales of HCN maps are available for tens of galaxies, higher-resolution maps still need to be available. Aims. We aim to study the properties of molecular gas, the contrast in intensity between two tracers that probe different density regimes (the HCN(1–0)/CO(2–1) ratio), and their kinematics across NGC 253, one of the closest starburst galaxies. With its advanced capabilities, the Atacama Large Millimeter/submillimeter Array (ALMA) can map these features at a high resolution across a large field of view and uncover the nature of such dense gas in extragalactic systems. Methods. We present new ALMA Atacama Compact Array and Total Power (ACA+TP) observations of the HCN emission across NGC 253. The observations cover the inner 8.6′ of the galaxy disk at a spatial resolution of 300 pc. Our study examines the distribution and kinematics of the HCN-traced gas and its relationship with the bulk molecular gas traced by CO(2–1). We analyze the integrated intensity and mean velocity of HCN and CO along each line of sight. We also used the SCOUSE software to perform spectral decomposition, which considers each velocity component separately. Results. We find that the denser molecular gas traced by HCN piles up in a ring-like structure at a radius of 2 kpc. The HCN emission is enhanced by two orders of magnitude in the central 2 kpc regions, beyond which its intensity decreases with increasing galactocentric distance. The number of components in the HCN spectra shows a robust environmental dependence, with multiple velocity features across the center and bar. The HCN spectra exhibit multiple velocity features across the center and bar, which shows a robust environmental dependence. We have identified an increase in the HCN/CO ratio in these regions, corresponding to a velocity component likely associated with a molecular outflow. We have also discovered that the ratio between the total infrared luminosity and dense gas mass, which is an indicator of the star formation efficiency of dense gas, is anticorrelated with the molecular gas surface density up to approximately 200 M⊙ pc−2. However, beyond this point, the ratio starts to increase. Conclusions. We argue that using information about spectroscopic features of molecular emission is an important aspect of understanding molecular properties in galaxies.

Constraints on Primordial Magnetic Fields from High Redshift Stellar Mass Density

Primordial magnetic fields (PMFs) play a pivotal role in influencing small-scale fluctuations within the primordial density field, thereby enhancing the matter power spectrum within the context of the ΛCDM model at small scales. These amplified fluctuations accelerate the early formation of galactic halos and stars, which can be observed through advanced high-redshift observational techniques. Therefore, stellar mass density (SMD) observations, which provide significant opportunities for detailed studies of galaxies at small scales and high redshifts, offer a novel perspective on small-scale cosmic phenomena and constrain the characteristics of PMFs. In this study, we compile 14 SMD data points at redshifts z > 6 and derive stringent constraints on the parameters of PMFs, which include the amplitude of the magnetic field at a characteristic scale of λ = 1 Mpc, denoted as B 0, and the spectral index of the magnetic field power spectrum, n B . At 95% confidence level, we establish upper limits of B 0 < 4.44 nG and n B < −2.24, along with a star formation efficiency of approximately f*0∼0.1 . If we fix n B at specific values, such as −2.85, −2.9, and −2.95, the 95% upper limits for the amplitude of the magnetic field can be constrained to 1.33, 2.21, and 3.90 nG, respectively. Finally, we attempt to interpret recent early observations provided by the James Webb Space Telescope using the theory of PMFs and find that by selecting appropriate PMF parameters, it is possible to explain these results without significantly increasing the star formation efficiency.

JWST PRIMER: a new multifield determination of the evolving galaxy UV luminosity function at redshifts z ≃ 9 – 15

ABSTRACT We present a new determination of the evolving galaxy ultraviolet (UV) luminosity function (LF) over the redshift range $8.5&amp;lt; z&amp;lt; 15.5$ using a combination of several major Cycle-1 JWST imaging programmes – Public Release IMaging for Extragalactic Research, JWST Advanced Deep Extragalactic Survey, and Next Generation Deep Extragalactic Exploratory Public Survey. This multifield approach yields a total of $\simeq 370$ arcmin2 of JWST/NIRCam imaging, reaching (5-$\sigma$) depths of $\simeq 30$ AB mag in the deepest regions. We select a sample of 2548 galaxies with a significant probability of lying at high redshift ($p(z&amp;gt; 8.5)&amp;gt; 0.05$) to undertake a statistical calculation of the UV LF. Our new measurements span $\simeq 4$ mag in UV luminosity at $z=9-12.5$, placing new constraints on both the shape and evolution of the LF at early times. Our measurements yield a new estimate of the early evolution of cosmic star-formation rate density ($\rho _{\rm {SFR}}$) confirming the gradual decline deduced from early JWST studies, at least out to $z \simeq 12$. Finally we show that the observed early evolution of the galaxy UV LF (and $\rho _{\rm {SFR}}$) can be reproduced in a ${\rm \Lambda }$cold dark matter Universe, with no change in dust properties or star-formation efficiency required out to $z \simeq 12$. Instead, a progressive trend towards younger stellar population ages can reproduce the observations, and the typical ages required at $z \simeq$ 8, 9, 10, and 11 all converge on $\simeq 380-330$ Myr after the big bang, indicative of a rapid emergence of early galaxies at $z \simeq 12 - 13$. This is consistent with the first indications of a steeper drop-off in $\rho _{\rm {SFR}}$ we find beyond $z \simeq 13$, possibly reflecting the rapid evolution of the halo mass function at earlier times.

Exploring the dark energy equation of state with JWST

Observations from the James Webb Space Telescope (JWST) have unveiled several galaxies with stellar masses M∗≳1010M⊙\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$M_*\\gtrsim 10^{10} M_\\odot$$\\end{document} at redshifts 7.4≲z≲9.1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$7.4\\lesssim z\\lesssim 9.1$$\\end{document}. These remarkable findings indicate an unexpectedly high stellar mass density, which contradicts the prediction of the ΛCDM\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Lambda \\rm{CDM}$$\\end{document} model. Our study utilizes the Chevallier–Polarski–Linder (CPL) parameterization, one of the dynamic dark energy models, to probe the role of dark energy on shaping galaxy formation. By considering varying star formation efficiencies within this framework, our analysis demonstrates that in a universe with a higher proportion of dark energy, more massive galaxies are formed at high redshifts, given a fixed perturbation amplitude observed today. These intriguing results highlight the promising prospect of revealing the nature of dark energy by analyzing the high-redshift massive galaxies.

The ratio of [Eu/α] differentiates accreted/in-situ Milky Way stars across metallicities, as indicated by both field stars and globular clusters

Abstract We combine stellar orbits with the abundances of the heavy, r-process element europium and the light, α-element, silicon to separate in-situ and accreted populations in the Milky Way across all metallicities. At high orbital energy, the accretion-dominated halo shows elevated values of [Eu/Si], while at lower energies, where many of the stars were born in-situ, the levels of [Eu/Si] are lower. These systematically different levels of [Eu/Si] in the MW and the accreted halo imply that the scatter in [Eu/α] within a single galaxy is smaller than previously thought. At the lowest metallicities, we find that both accreted and in-situ populations trend down in [Eu/Si], consistent with enrichment via neutron star mergers. Through compiling a large dataset of abundances for 54 globular clusters (GCs), we show that differences in [Eu/Si] extend to populations of in-situ/accreted GCs. We interpret this consistency as evidence that in r-process elements GCs trace the star formation history of their hosts, motivating their use as sub-Gyr timers of galactic evolution. Furthermore, fitting the trends in [Eu/Si] using a simple galactic chemical evolution model, we find that differences in [Eu/Si] between accreted and in-situ MW field stars cannot be explained through star formation efficiency alone. Finally, we show that the use of [Eu/Si] as a chemical tag between GCs and their host galaxies extends beyond the Local Group, to the halo of M31 - potentially offering the opportunity to do Galactic Archaeology in an external galaxy.

Understanding the star formation efficiency in dense gas: Initial results from the CAFFEINE survey with ArTéMiS

Context. Despite recent progress, the question of what regulates the star formation efficiency (SFE) in galaxies remains one of the most debated problems in astrophysics. According to the dominant picture, star formation (SF) is regulated by turbulence and feedback, and the SFE is ~1–2% or less per local free-fall time on all scales from Galactic clouds to high-redshift galaxies. In an alternate scenario, the star formation rate (SFR) in galactic disks is linearly proportional to the mass of dense gas above some critical density threshold ~104 cm–3. Aims. We aim to discriminate between these two pictures thanks to high-resolution submillimeter and mid-infrared imaging observations, which trace both dense gas and young stellar objects (YSOs) for a comprehensive sample of 49 nearby massive SF complexes out to a distance of d ~ 3 kpc in the Galactic disk. Methods. We used data from CAFFEINE, a complete 350/450 µm survey with APEX/ArTéMiS of the densest portions of all southern molecular clouds at d ≲ 3 kpc, in combination with Herschel data to produce column density maps at a factor of ~4 higher resolution (8") than standard Herschel column density maps (36″). Our maps are free of any saturation effect around luminous high-mass pro-tostellar objects and resolve the structure of dense gas and the typical ~0.1 pc width of molecular filaments out to 3 kpc, which is the most important asset of the present study and is impossible to achieve with Herschel data alone. Coupled with SFR estimates derived from Spitzer mid-infrared observations of the YSO content of the same clouds, this allowed us to study the dependence of the SFE on density in the CAFFEINE clouds. We also combine our findings with existing SF efficiency measurements in nearby clouds to extend our analysis down to lower column densities. Results. Our results suggest that the SFE does not increase with density above the critical threshold and support a scenario in which the SFE in dense gas is approximately constant (independent of free-fall time). However, the SF efficiency measurements traced by Class I YSOs in nearby clouds are more inconclusive, since they are consistent with both the presence of a density threshold and a dependence on density above the threshold. Overall, we suggest that the SF efficiency in dense gas is primarily governed by the physics of filament fragmentation into protostellar cores.

Star Formation by Supernova Implosion

Recent hydrodynamical simulations of the late stages of supernova remnant (SNR) evolution have revealed that as they merge with the ambient medium, SNRs implode, leading to the formation of dense clouds in their center. While being highly chemically enriched by their host SNR, these clouds appear to have similar properties as giant molecular clouds, which are believed to be the main site of star formation. Here, we develop a simple model in order to estimate the efficiency of the star formation that might be triggered by the implosion of SNRs. We separately consider two cases: cyclic star formation, maintained by the episodic driving of feedback from new generations of stars, and a single burst of star formation, triggered by a single explosion. We find that in the cyclic case, star formation is inefficient, with implosion-triggered star formation contributing a few percent of the observed star formation efficiency per freefall timescale. In the single-burst case, higher star formation efficiencies can be obtained. However, while the implosion-triggered process might not contribute much to the overall star formation, due to the high chemical enrichment of the birth clouds, it can explain the formation of a significant fraction of metal-rich stars.

Self-similar cluster structures in massive star-forming regions: Isolated evolution from clumps to embedded clusters

We used the dendrogram algorithm to decompose the surface density distributions of stars into hierarchical structures. These structures were tied to the multiscale structures of star clusters. A similar power-law for the mass-size relation of star clusters measured at different scales suggests a self-similar structure of star clusters. We used the minimum spanning tree method to measure the separations between clusters and gas clumps in each massive star-forming region. The separations between clusters, between clumps, and between clusters and clumps were comparable, which indicates that the evolution from clump to embedded cluster proceeds in isolation and locally, and does not affect the surrounding objects significantly. By comparing the mass functions of the ATLASGAL clumps and the identified embedded clusters, we confirm that a constant star formation efficiency of approx 0.33 can be a typical value for the ATLASGAL clumps.

Tests of subgrid models for star formation using simulations of isolated disc galaxies

ABSTRACT We use smoothed particle hydrodynamics simulations of isolated Milky Way-mass disc galaxies that include cold, interstellar gas to test subgrid prescriptions for star formation (SF). Our fiducial model combines a Schmidt law with a gravitational instability criterion, but we also test density thresholds and temperature ceilings. While SF histories are insensitive to the prescription for SF, the Kennicutt–Schmidt (KS) relations between SF rate and gas surface density can discriminate between models. We show that our fiducial model, with an SF efficiency per free-fall time of 1 per cent, agrees with spatially resolved and azimuthally averaged observed KS relations for neutral, atomic, and molecular gas. Density thresholds do not perform as well. While temperature ceilings selecting cold, molecular gas can match the data for galaxies with solar metallicity, they are unsuitable for very low-metallicity gas and hence for cosmological simulations. We argue that SF criteria should be applied at the resolution limit rather than at a fixed physical scale, which means that we should aim for numerical convergence of observables rather than of the properties of gas labelled as star-forming. Our fiducial model yields good convergence when the mass resolution is varied by nearly 4 orders of magnitude, with the exception of the spatially resolved molecular KS relation at low surface densities. For the gravitational instability criterion, we quantify the impact on the KS relations of gravitational softening, the SF efficiency, and the strength of supernova feedback, as well as of observable parameters such as the inclusion of ionized gas, the averaging scale, and the metallicity.

Climbing the Cliffs: Classifying Young Stellar Objects in the Cosmic Cliffs JWST Data Using a Probabilistic Random Forest

Among the first observations released to the public from the JWST was a section of the star-forming region NGC 3324 known colloquially as the “Cosmic Cliffs.” We build a photometric catalog of the region and test the ability of using the probabilistic random forest machine-learning method to identify its young stellar objects (YSOs). We find 450 candidate YSOs (cYSOs) out of 19,497 total objects within the field, 413 of which are cYSOs not found in previous works. These classifications are verified with several different metrics, including recall and precision. Using the obtained probabilities of objects being YSOs, we employ a Monte Carlo approach to determine the surface density of cYSOs in the Cosmic Cliffs, which we find to be largely coincident with column densities derived from Herschel data, up to a column density of 1.37 × 1022 cm−2. The newly determined number and spatial distribution of YSOs in the Cosmic Cliffs demonstrate that JWST is far more capable of detecting YSOs in dusty regions than Spitzer. Comparisons of the observed colors and brightness of faint cYSOs with those of pre-main-sequence models suggest JWST has detected a significant population of substellar YSOs in the Cosmic Cliffs. The size of this population further suggests previous estimates of star formation efficiencies in molecular clouds have been systematically low.

A Temperature or Far-ultraviolet Tracer? The HNC/HCN Ratio in M83 on the Scale of Giant Molecular Clouds

The HNC/HCN ratio is observationally known as a thermometer in Galactic interstellar molecular clouds. A recent study has alternatively suggested that the HNC/HCN ratio is affected by the ultraviolet (UV) field, not by the temperature. We aim to study this ratio on the scale of giant molecular clouds in the barred spiral galaxy M83 towards the southwestern bar end and the central region from Atacama Large Millimeter/submillimeter Array observations, and if possible, distinguish the above scenarios. We compare the high-resolution (40–50 pc) HNC/HCN ratios with the star formation rate from the 3 mm continuum intensity and the molecular mass inferred from the HCN intensities. Our results show that the HNC/HCN ratios do not vary with the star formation rates, star formation efficiencies, or column densities in the bar-end region. In the central region, the HNC/HCN ratios become higher with higher star formation rates, which tend to cause higher temperatures. This result is not consistent with the previously proposed scenario in which the HNC/HCN ratio decreases with increasing temperature. Spectral shapes suggest that this trend may be due to optically thick HCN and optically thin HNC. In addition, we compare the large-scale (∼200 pc) correlation between the dust temperature from the far-IR ratio and the HNC/HCN ratio for the southwestern bar-end region. The HNC/HCN ratio is lower when the dust temperatures are higher. We suggest from the above results that the HNC/HCN ratio depends on the UV radiation field that affects the interstellar medium on the ∼100 pc scale where the column densities are low.

Modelling the star-formation activity and ionizing properties of high-redshift galaxies

Early results from the JWST observations have reported a surprisingly high number of UV-bright galaxies at z ≥ 10, which appears to challenge the theoretical predictions from standard galaxy formation models in the ΛCDM framework at these redshifts. To alleviate this tension, several cosmological and astrophysical interpretations have been advanced. However, all of these proposed scenarios carry noteworthy consequences for other large-scale processes in the early Universe, particularly cosmic reionization, since high-redshift galaxies are believed to be the primary ionizing sources during the Epoch of Reionization (EoR). To investigate this, we introduce a semi-analytical model of galaxy formation and evolution that explains the evolving galaxy UV luminosity function (UVLF) over 6 ≲ z ≲ 15, and also jointly tracks the time evolution of the globally averaged neutral hydrogen fraction in the intergalactic medium. The model self-consistently accounts for the suppression of star formation in low-mass galaxies due to reionization feedback and is constrained by comparing the model predictions with various observational probes like the UVLF data from HST and JWST, recent measurements of the neutral hydrogen fraction, and the CMB scattering optical depth. Our analysis confirms that a rapid enhancement in the star-formation rate efficiency and/or UV luminosity per stellar mass formed is necessary for consistency with the JWST UVLF estimates at z ≥ 10. We further find that it is possible to jointly satisfy the current reionization constraints when the escape fraction is assumed to be halo-mass dependent, requiring higher Lyman-continuum leakage from low-mass galaxies. We also examine the relative contribution of galaxies with different UV luminosities towards the ionizing photon budget for the EoR and investigate the large-scale bias of high-z galaxies.

Emergence of high-mass stars in complex fiber networks (EMERGE)

Context. Recent molecular surveys have revealed the rich gas organization of sonic-like filaments at small scales (so-called fibers) in all types of environments prior to the formation of low- and high-mass stars. These fibers form at the end of the turbulent cascade and are identified as the fine substructure within the hierarchical nature of the gas in the interstellar medium (ISM). Aims. Isolated fibers provide the subsonic conditions for the formation of low-mass stars. This paper introduces the Emergence of high-mass stars in complex fiber networks (EMERGE) project, which investigates whether complex fiber arrangements (networks) can also explain the origin of high-mass stars and clusters. Methods. We analyzed the EMERGE Early ALMA Survey including seven star-forming regions in Orion (OMC-1,2,3, and 4 South, LDN 1641N, NGC 2023, and the Flame Nebula) that were homogeneously surveyed in three molecular lines (N2H+ J = 1–0, HNC J = 1–0, and HC3N J = 10–9) and in the 3 mm continuum using a combination of interferometric ALMA mosaics and IRAM-30 m single-dish (SD) maps, together with a series of Herschel, Spitzer, and WISE archival data. We also developed a systematic data reduction framework allowing the massive data processing of ALMA observations. Results. We obtained independent continuum maps and spectral cubes for all our targets and molecular lines at different (SD and interferometric) resolutions, and we explored multiple data combination techniques. Based on our low-resolution (SD) observations (30″ or ~12 000 au), we describe the global properties of our sample, which covers a wide range of physical conditions, including low-(OMC-4 South and NGC 2023), intermediate (OMC-2, OMC-3, and LDN 1641N), and high-mass (OMC-1 and Flame Nebula) star-forming regions in different evolutionary stages. The comparison between our single-dish maps and ancillary YSO catalogs denotes N2H+ (1–0) as the best proxy for the dense, star-forming gas in our targets, which show a constant star formation efficiency and a fast time evolution of ≲1 Myr. While apparently clumpy and filamentary in our SD data, all targets show a much more complex fibrous substructure at the enhanced resolution of our combined ALMA+IRAM-30 m maps (4″.5 or ~2000 au). A large number of filamentary features at subparsec scales are clearly recognized in the high-density gas (≳ 105 cm−3) that is traced by N2H+ (1–0) directly connected to the formation of individual protostars. Surprisingly, this complex gas organization appears to extend farther into the more diffuse gas (~103−104 cm−3) traced by HNC (1–0). Conclusions. This paper presents the EMERGE Early ALMA Survey, which includes a first data release of continuum maps and spectral products for this project that are to be analysed in future papers of this series. A first look at these results illustrates the need of advanced data combination techniques between high-resolution interferometric (ALMA) and high-sensitivity, low-resolution single-dish (IRAM-30 m) datasets to investigate the intrinsic multiscale, gas structure of the ISM.

A Comprehensive Metallicity Analysis of J0332−3557: Establishing a z ∼ 4 Anchor for Direct Gas Metallicity and C/O Abundance Investigations

We provide one of the most comprehensive metallicity studies at z ∼ 4 by analyzing the UV/optical Hubble Space Telescope photometry and rest-frame Very Large Telescope (VLT)/FORS2 UV and VLT/XSHOOTER optical spectra of J0332−3557, a gravitationally lensed galaxy magnified by a factor of 20. With a 5σ detection of the auroral O iii] λ1666 line, we are able to derive a direct gas metallicity estimate for our target. We find Z gas =12+log(O/H)=8.26±0.06 , which is compatible with an increase of both the gas fraction and the outflow metal loading factor from z ∼ 0 to z ∼ 4. J0332−3557 is the most metal-rich individual galaxy at z ∼ 4 for which the C/O ratio has been measured. We derive a low log(C/O) = −1.02 ± 0.2, which suggests that J0332−3557 is in the early stages of interstellar medium carbon enrichment driven mostly by massive stars. The low C/O abundance also indicates that J0332−3557 is characterized by a low star formation efficiency, higher yields of oxygen, and longer burst duration. We find that EWC III]1906,9 is as low as ∼3 Å, and the main drivers of the low EWC III]1906,9 are the higher gas metallicity and the low C/O abundance. J0332−3557 is characterized by one diffuse and two more compact regions ∼1 kpc in size. We find that the carbon emission mostly originates in the compact knots. Our study on J0332−3557 serves as an anchor for studies investigating the evolution of metallicity and C/O abundance across different redshifts.

The Complete CEERS Early Universe Galaxy Sample: A Surprisingly Slow Evolution of the Space Density of Bright Galaxies at z ∼ 8.5–14.5

We present a sample of 88 candidate z ∼ 8.5–14.5 galaxies selected from the completed NIRCam imaging from the Cosmic Evolution Early Release Science survey. These data cover ∼90 arcmin2 (10 NIRCam pointings) in six broadband imaging filters and one medium-band imaging filter. With this sample we confirm at higher confidence early JWST conclusions that bright galaxies in this epoch are more abundant than predicted by most theoretical models. We construct the rest-frame ultraviolet luminosity functions at z ∼ 9, 11, and 14 and show that the space density of bright (M UV = −20) galaxies changes only modestly from z ∼ 14 to z ∼ 9, compared to a steeper increase from z ∼ 8 to z ∼ 4. While our candidates are photometrically selected, spectroscopic follow-up has now confirmed 13 of them, with only one significant interloper, implying that the fidelity of this sample is high. Successfully explaining the evidence for a flatter evolution in the number densities of UV-bright z > 10 galaxies may thus require changes to the dominant physical processes regulating star formation. While our results indicate that significant variations of dust attenuation with redshift are unlikely to be the dominant factor at these high redshifts, they are consistent with predictions from models that naturally have enhanced star formation efficiency and/or stochasticity. An evolving stellar initial mass function could also bring model predictions into better agreement with our results. Deep spectroscopic follow-up of a large sample of early galaxies can distinguish between these competing scenarios.