Showing papers on "Star formation published in 2022"

The First Stars

Abstract: ▪ Abstract We review recent theoretical results on the formation of the first stars in the universe, and emphasize related open questions. In particular, we discuss the initial conditions for Popul...

The PHANGS-MUSE survey

Abstract: We present the PHANGS-MUSE survey, a programme using the MUSE IFS at the ESO VLT to map 19 massive $(9.4 < \log(M_{*}/M_\odot) < 11.0)$ nearby (D < 20 Mpc) star-forming disc galaxies. The survey consists of 168 MUSE pointings (1'x1' each), a total of nearly 15 Million spectra, covering ~1.5 Million independent spectra. PHANGS-MUSE provides the first IFS view of star formation across different local environments (including galaxy centres, bars, spiral arms) in external galaxies at a median resolution of 50~pc, better than the mean inter-cloud distance in the ionised interstellar medium. This `cloud-scale' resolution allows detailed demographics and characterisations of HII regions and other ionised nebulae. PHANGS-MUSE further delivers a unique view on the associated gas and stellar kinematics, and provides constraints on the star formation history. The PHANGS-MUSE survey is complemented by dedicated ALMA CO(2-1) and multi-band HST observations, therefore allowing us to probe the key stages of the star formation process from molecular clouds to HII regions and star clusters. This paper describes the scientific motivation, sample selection, observational strategy, data reduction and analysis process of the PHANGS-MUSE survey. We present our bespoke automated data-reduction framework, which is built on the reduction recipes provided by ESO, but additionally allows for mosaicking and homogenisation of the point spread function. We further present a detailed quality assessment and a brief illustration of the potential scientific applications of the large set of PHANGS-MUSE data products generated by our data analysis framework. The data cubes and analysis data products described in this paper represent the basis for the first PHANGS-MUSE public data release and are available in the ESO archive and via the Canadian Astronomy Data Centre.

A Comprehensive Study of Galaxies at z ∼ 9–16 Found in the Early JWST Data: Ultraviolet Luminosity Functions and Cosmic Star Formation History at the Pre-reionization Epoch

Abstract: We conduct a comprehensive study on dropout galaxy candidates at z ∼ 9–16 using the first 90 arcmin2 James Webb Space Telescope (JWST) Near Infrared Camera images taken by the early release observations (ERO) and early release science programs. With the JWST simulation images, we find that a number of foreground interlopers are selected with a weak photo-z determination (Δχ 2 > 4). We thus carefully apply a secure photo-z selection criterion (Δχ 2 > 9) and conventional color criteria with confirmations of the ERO Near Infrared Spectrograph spectroscopic redshifts, and obtain a total of 23 dropout galaxies at z ∼ 9–16, including two candidates at zphot=16.25−0.46+0.24 and 16.41−0.55+0.66 . We perform thorough comparisons of dropout galaxies found in our work with recent JWST studies, and conclude that our galaxy sample is reliable enough for statistical analyses. We derive the UV luminosity functions at z ∼ 9–16, and confirm that our UV luminosity functions at z ∼ 9 and 12 agree with those determined by other Hubble Space Telescope and JWST studies. The cosmic star formation rate (SFR) density decreases from z ∼ 9 to 12, and perhaps to 16, but the densities at z ∼ 12–16 are higher than the constant star formation efficiency model. Interestingly, there are six bright galaxy candidates at z ∼ 10–16 with M UV < −19.5 mag and M * ∼ 108−9 M ⊙. Because a majority (∼80%) of these galaxies show no signatures of active galactic nuclei in their morphologies, the high cosmic SFR densities and the existence of these UV-luminous galaxies are explained by the lack of suppression of star formation by the UV background radiation at the pre-reionization epoch and/or an efficient UV radiation production by a top-heavy initial mass function with Population III–like star formation.

Star formation near the Sun is driven by expansion of the Local Bubble

Abstract: For decades we have known that the Sun lies within the Local Bubble, a cavity of low-density, high-temperature plasma surrounded by a shell of cold, neutral gas and dust1-3. However, the precise shape and extent of this shell4,5, the impetus and timescale for its formation6,7, and its relationship to nearby star formation8 have remained uncertain, largely due to low-resolution models of the local interstellar medium. Here we report an analysis of the three-dimensional positions, shapes and motions of dense gas and young stars within 200 pc of the Sun, using new spatial9-11 and dynamical constraints12. We find that nearly all of the star-forming complexes in the solar vicinity lie on the surface of the Local Bubble and that their young stars show outward expansion mainly perpendicular to the bubble's surface. Tracebacks of these young stars' motions support a picture in which the origin of the Local Bubble was a burst of stellar birth and then death (supernovae) taking place near the bubble's centre beginning approximately 14 Myr ago. The expansion of the Local Bubble created by the supernovae swept up the ambient interstellar medium into an extended shell that has now fragmented and collapsed into the most prominent nearby molecular clouds, in turn providing robust observational support for the theory of supernova-driven star formation.

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

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

Star formation near the Sun is driven by expansion of the Local Bubble

Abstract: For decades we have known that the Sun lies within the Local Bubble, a cavity of low-density, high-temperature plasma surrounded by a shell of cold, neutral gas and dust1-3. However, the precise shape and extent of this shell4,5, the impetus and timescale for its formation6,7, and its relationship to nearby star formation8 have remained uncertain, largely due to low-resolution models of the local interstellar medium. Here we report an analysis of the three-dimensional positions, shapes and motions of dense gas and young stars within 200 pc of the Sun, using new spatial9-11 and dynamical constraints12. We find that nearly all of the star-forming complexes in the solar vicinity lie on the surface of the Local Bubble and that their young stars show outward expansion mainly perpendicular to the bubble's surface. Tracebacks of these young stars' motions support a picture in which the origin of the Local Bubble was a burst of stellar birth and then death (supernovae) taking place near the bubble's centre beginning approximately 14 Myr ago. The expansion of the Local Bubble created by the supernovae swept up the ambient interstellar medium into an extended shell that has now fragmented and collapsed into the most prominent nearby molecular clouds, in turn providing robust observational support for the theory of supernova-driven star formation.

Fast, Slow, Early, Late: Quenching Massive Galaxies at z ∼ 0.8

Abstract: We investigate the stellar populations for a sample of 161 massive, mainly quiescent galaxies at $\langle z_{\rm obs} \rangle=0.8$ with deep Keck/DEIMOS rest-frame optical spectroscopy (HALO7D survey). With the fully Bayesian framework Prospector, we simultaneously fit the spectroscopic and photometric data with an advanced physical model (including non-parametric star-formation histories, emission lines, variable dust attenuation law, and dust and AGN emission) together with an uncertainty and outlier model. We show that both spectroscopy and photometry are needed to break the dust-age-metallicity degeneracy. We find a large diversity of star-formation histories: although the most massive ($M_{\star}>2\times10^{11}~M_{\odot}$) galaxies formed the earliest (formation redshift of $z_{\rm f}\approx5-10$ with a short star-formation timescale of $\tau_{\rm SF}\lesssim1~\mathrm{Gyr}$), lower-mass galaxies have a wide range of formation redshifts, leading to only a weak trend of $z_{\rm f}$ with $M_{\star}$. Interestingly, several low-mass galaxies with have formation redshifts of $z_{\rm f}\approx5-8$. Star-forming galaxies evolve about the star-forming main sequence, crossing the ridgeline several times in their past. Quiescent galaxies show a wide range and continuous distribution of quenching timescales ($\tau_{\rm quench}\approx0-5~\mathrm{Gyr}$) with a median of $\langle\tau_{\rm quench}\rangle=1.0_{-0.9}^{+0.8}~\mathrm{Gyr}$ and of quenching epochs of $z_{\rm quench}\approx0.8-5.0$ ($\langle z_{\rm quench}\rangle=1.3_{-0.4}^{+0.7}$). This large diversity of quenching timescales and epochs points toward a combination of internal and external quenching mechanisms. In our sample, rejuvenation and "late bloomers" are uncommon. In summary, our analysis supports the "grow & quench" framework and is consistent with a wide and continuously-populated diversity of quenching timescales.

Ram pressure stripping in high-density environments

Abstract: Galaxies living in rich environments are suffering different perturbations able to drastically affect their evolution. Among these, ram pressure stripping, i.e. the pressure exerted by the hot and dense intracluster medium (ICM) on galaxies moving at high velocity within the cluster gravitational potential well, is a key process able to remove their interstellar medium (ISM) and quench their activity of star formation. This review is aimed at describing this physical mechanism in different environments, from rich clusters of galaxies to loose and compact groups. We summarise the effects of this perturbing process on the baryonic components of galaxies, from the different gas phases (cold atomic and molecular, ionised, hot) to magnetic fields and cosmic rays, and describe their induced effects on the different stellar populations, with a particular attention to its role in the quenching episode generally observed in high-density environments. We also discuss on the possible fate of the stripped material once removed from the perturbed galaxies and mixed with the ICM, and we try to estimate its contribution to the pollution of the surrounding environment. Finally, combining the results of local and high-redshift observations with the prediction of tuned models and simulations, we try to quantify the importance of this process on the evolution of galaxies of different mass, from dwarfs to giants, in various environments and at different epochs.

A time-resolved picture of our Milky Way’s early formation history

Abstract: The formation of our Milky Way can be split up qualitatively into different phases that resulted in its structurally different stellar populations: the halo and the disk components1-3. Revealing a quantitative overall picture of our Galaxy's assembly requires a large sample of stars with very precise ages. Here we report an analysis of such a sample using subgiant stars. We find that the stellar age-metallicity distribution p(τ, [Fe/H]) splits into two almost disjoint parts, separated at age τ ≃ 8 Gyr. The younger part reflects a late phase of dynamically quiescent Galactic disk formation with manifest evidence for stellar radial orbit migration4-6; the other part reflects the earlier phase, when the stellar halo7 and the old α-process-enhanced (thick) disk8,9 formed. Our results indicate that the formation of the Galaxy's old (thick) disk started approximately 13 Gyr ago, only 0.8 Gyr after the Big Bang, and 2 Gyr earlier than the final assembly of the inner Galactic halo. Most of these stars formed around 11 Gyr ago, when the Gaia-Sausage-Enceladus satellite merged with our Galaxy10,11. Over the next 5-6 Gyr, the Galaxy experienced continuous chemical element enrichment, ultimately by a factor of 10, while the star-forming gas managed to stay well mixed.

GOLDRUSH. IV. Luminosity Functions and Clustering Revealed with ∼4,000,000 Galaxies at z ∼ 2–7: Galaxy–AGN Transition, Star Formation Efficiency, and Implication for Evolution at z &gt; 10

Abstract: We present new measurements of rest-UV luminosity functions and angular correlation functions from 4,100,221 galaxies at z~2-7 identified in the Subaru/Hyper Suprime-Cam survey and CFHT Large-Area U-band Survey. The obtained luminosity functions at z~4-7 cover a very wide UV luminosity range of ~0.002-2000L*uv combined with previous studies, revealing that the dropout luminosity function is a superposition of the AGN luminosity function dominant at Muv<-24 mag and the galaxy luminosity function dominant at Muv>-22 mag, consistent with galaxy fractions based on 1037 spectroscopically-identified sources. Galaxy luminosity functions estimated from the spectroscopic galaxy fractions show the bright end excess beyond the Schechter function at >2sigma levels, which is possibly made by inefficient mass quenching, low dust obscuration, and/or hidden AGN activity. By analyzing the correlation functions at z~2-6 with halo occupation distribution models, we find a weak redshift evolution (within 0.3 dex) of the ratio of the star formation rate (SFR) to the dark matter accretion rate, SFR/(dMh/dt), indicating the almost constant star formation efficiency at z~2-6, as suggested by our earlier work at z~4-7. Meanwhile, the ratio gradually increases with decreasing redshift at z<5 within 0.3 dex, which quantitatively reproduces the redshift evolution of the cosmic SFR density, suggesting that the evolution is primarily driven by the increase of the halo number density due to the structure formation, and the decrease of the accretion rate due to the cosmic expansion. Extrapolating this calculation to higher redshifts assuming the constant efficiency suggests a rapid decrease of the SFR density at z>10 with $\propto10^{-0.5(1+z)}$, which will be directly tested with JWST.

On the Stellar Populations of Galaxies at z = 9–11: The Growth of Metals and Stellar Mass at Early Times

Abstract: We present a detailed stellar population analysis of 11 bright ($H<26.6$) galaxies at $z=9-11$ (three spectroscopically confirmed) to constrain the chemical enrichment and growth of stellar mass of early galaxies. We use the flexible Bayesian spectral energy distribution (SED) fitting code Prospector with a range of star-formation histories (SFHs), a flexible dust attenuation law and a self-consistent modeling of emission lines. This approach allows us to assess how different priors affect our results, and how well we can break degeneracies between dust attenuation, stellar ages, metallicity and emission lines using data which probe only the rest-frame ultraviolet to optical wavelengths. We measure a median observed ultraviolet spectral slope $\beta=-1.87_{-0.43}^{+0.35}$ for relatively massive star-forming galaxies ($9<\log(M_{\star}/M_{\odot})<10$), consistent with no change from $z=4$ to $z=9-10$ at these stellar masses, implying rapid enrichment. Our SED-fitting results are consistent with a star-forming main sequence with sub-linear slope ($0.7\pm0.2$) and specific star-formation rates of $3-10~\mathrm{Gyr}^{-1}$. However, the stellar ages and SFHs are less well constrained. Using different SFH priors, we cannot distinguish between median mass-weighted ages of $\sim50-150$ Myr, which corresponds to 50\% formation redshifts of $z_{50}\sim10-12$ at $z\sim9$ and is of the order of the dynamical timescales of these systems. Importantly, the models with different SFH priors are able to fit the data equally well. We conclude that the current observational data cannot tightly constrain the mass-buildup timescales of these $z=9-11$ galaxies, with our results consistent with SFHs implying both a shallow and steep increase of the cosmic SFR density with time at $z>10$.

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

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

Galaxy Formation and Reionization: Key Unknowns and Expected Breakthroughs by the <i>James Webb Space Telescope</i>

Abstract: The launch of the James Webb Space Telescope (JWST) in late 2021 marks a new start for studies of galaxy formation at high redshift ( z ≳ 6) during the era of cosmic reionization. JWST can capture sensitive, high-resolution images and multiobject spectroscopy in the IR that will transform our view of galaxy formation during the first billion years of cosmic history. This review summarizes our current knowledge of the role of galaxies in reionizing intergalactic hydrogen ahead of JWST, achieved through observations with the Hubble Space Telescope and ground-based facilities including Keck, the Very Large Telescope, Subaru, and the Atacama Large Millimeter/Submillimeter Array. We identify outstanding questions in the field that JWST can address during its mission lifetime, including with the planned JWST Cycle 1 programs. These findings include the following: ▪ Surveys with JWST have sufficient sensitivity and area to complete the census of galaxy formation at the current redshift frontier ( z ∼ 8–10). ▪ Rest-frame optical spectroscopy with JWST of galaxies will newly enable measures of star-formation rate, metallicity, and ionization at z ∼ 8–9, allowing for the astrophysics of early galaxies to be constrained. ▪ The presence of evolved stellar populations at z ∼ 8–10 can be definitively tested by JWST, which would provide evidence of star formation out to z ∼ 15.

The PHANGS-HST Survey: Physics at High Angular Resolution in Nearby Galaxies with the Hubble Space Telescope

Abstract: The PHANGS program is building the first dataset to enable the multi-phase, multi-scale study of star formation across the nearby spiral galaxy population. This effort is enabled by large survey programs with ALMA, VLT/MUSE, and HST, with which we have obtained CO(2-1) imaging, optical spectroscopic mapping, and high resolution UV-optical imaging, respectively. Here, we present PHANGS-HST, which has obtained five band NUV-U-B-V-I imaging of the disks of 38 spiral galaxies at distances of 4-23 Mpc, and parallel V and I band imaging of their halos, to provide a census of tens of thousands of compact star clusters, and multi-scale stellar associations. The combination of HST, ALMA, and VLT/MUSE observations will yield an unprecedented joint catalog of the observed and physical properties of ~100,000 star clusters, associations, HII regions, and molecular clouds. With these basic units of star formation, PHANGS will systematically chart the evolutionary cycling between gas and stars, across a diversity of galactic environments found in nearby galaxies. We discuss the design of the PHANGS-HST survey, and provide an overview of the HST data processing pipeline and first results, highlighting new methods for selecting star cluster candidates, morphological classification of candidates with convolutional neural networks, and identification of stellar associations over a range of physical scales with a watershed algorithm. We describe the cross-observatory imaging, catalogs, and software products to be released. These high-level science products will seed a broad range of investigations, in particular, the study of embedded stellar populations and dust with JWST, for which a PHANGS Cycle 1 Treasury program to obtain eight band 2-21 $\mu$m imaging has been approved.

The Cold Interstellar Medium of Galaxies in the Local Universe

Abstract: The cold interstellar medium (ISM) plays a central role in the galaxy evolution process. It is the reservoir that fuels galaxy growth via star formation, the repository of material formed by these stars, and a sensitive tracer of internal and external processes that affect entire galaxies. Consequently, significant efforts have gone into systematic surveys of the cold ISM of the galaxies in the local Universe. This review discusses the resulting network of scaling relations connecting the atomic and molecular gas masses of galaxies with their other global properties (stellar masses, morphologies, metallicities, star-formation activity…) and their implications for our understanding of galaxy evolution. Key take-home messages are as follows: ▪ From a gas perspective, there are three main factors that determine the star-formation rate of a galaxy: the total mass of its cold ISM, how much of that gas is molecular, and the rate at which any molecular gas is converted into stars. All three of these factors vary systematically across the local galaxy population. ▪ The shape and scatter of both the star-formation main sequence and the mass–metallicity relation are deeply linked to the availability of atomic and molecular gas. ▪ Future progress will come from expanding our exploration of scaling relations into new parameter space (in particular, the regime of dwarf galaxies), better connecting the cold ISM of large samples of galaxies with the environment that feeds them (the circumgalactic medium, in particular), and understanding the impact of these large scales on the efficiency of the star-formation process on molecular cloud scales. Expected final online publication date for the Annual Review of Astronomy and Astrophysics Volume 60 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

A first look at the SMACS0723 JWST ERO: spectroscopic redshifts, stellar masses and star-formation histories

Abstract: We present a first-look analysis of the JWST ERO data in the SMACS J0723.3-7327 cluster field. We begin by reporting 10 new spectroscopic redshifts from 𝜆 obs = 1 . 8 − 5 . 2 𝜇 m NIRSpec medium-resolution ( 𝑅 = 𝜆 / Δ 𝜆 = 1000) data. These are all determined via multiple high-SNR emission line detections, with 5 objects at 1 < 𝑧 < 3 displaying multiple rest-frame near-infrared Hydrogen Paschen lines, and 5 objects at 5 < 𝑧 < 9 displaying rest-frame optical Oxygen and Hydrogen Balmer lines. For the 5 higher-redshift galaxies we extract fluxes in 6 NIRCam bands spanning 𝜆 obs = 0 . 8 − 5 𝜇 m and perform spectral energy distribution fitting using these data in combination with existing HST photometry. The 7 < 𝑧 < 9 objects exhibit a U-shaped pattern across the F277W, F356W and F444W bands, indicating the presence of a Balmer break seen in emission (Balmer jump) and high-equivalent-width [O iii] emission. This is indicative of an extremely young stellar population, with the bulk of the current mass having formed within the past 10 Myr. We report robust stellar masses and mean stellar ages from our spectral fitting, with the four 𝑧 > 6 galaxies exhibiting low stellar masses from log 10 ( 𝑀 ∗ / M (cid:12) ) = 7 . 4 − 8 . 6 and correspondingly young mean stellar ages of only a few Myr. This work highlights the critical importance of combining large upcoming NIRCam surveys with NIRSpec follow-up to measure the spectroscopic redshifts necessary to robustly constrain physical parameters.

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

Abstract: 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 ∝ N H (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.

CEERS Key Paper. I. An Early Look into the First 500 Myr of Galaxy Formation with JWST

Abstract: We present an investigation into the first 500 Myr of galaxy evolution from the Cosmic Evolution Early Release Science (CEERS) survey. CEERS, one of 13 JWST ERS programs, targets galaxy formation from z ∼ 0.5 to >10 using several imaging and spectroscopic modes. We make use of the first epoch of CEERS NIRCam imaging, spanning 35.5 arcmin2, to search for candidate galaxies at z > 9. Following a detailed data reduction process implementing several custom steps to produce high-quality reduced images, we perform multiband photometry across seven NIRCam broad- and medium-band (and six Hubble broadband) filters focusing on robust colors and accurate total fluxes. We measure photometric redshifts and devise a robust set of selection criteria to identify a sample of 26 galaxy candidates at z ∼ 9–16. These objects are compact with a median half-light radius of ∼0.5 kpc. We present an early estimate of the z ∼ 11 rest-frame ultraviolet (UV) luminosity function, finding that the number density of galaxies at M UV ∼ −20 appears to evolve very little from z ∼ 9 to 11. We also find that the abundance (surface density [arcmin−2]) of our candidates exceeds nearly all theoretical predictions. We explore potential implications, including that at z > 10, star formation may be dominated by top-heavy initial mass functions, which would result in an increased ratio of UV light per unit halo mass, though a complete lack of dust attenuation and/or changing star formation physics may also play a role. While spectroscopic confirmation of these sources is urgently required, our results suggest that the deeper views to come with JWST should yield prolific samples of ultrahigh-redshift galaxies with which to further explore these conclusions.

OUP accepted manuscript

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 formation. Using a new method based on simultaneous [CII] 158$\mu$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.

Photodissociation and X-Ray-Dominated Regions

Abstract: The radiation from stars and active galactic nuclei (AGNs) creates photodissociation regions (PDRs) and X-ray-dominated regions (XDRs), where the chemistry or heating are dominated by far-ultraviolet (FUV) radiation or X-ray radiation, respectively. PDRs include a wide range of environments, from the diffuse interstellar medium to dense star-forming regions. XDRs are found in the center of galaxies hosting AGNs, in protostellar disks, and in the vicinity of X-ray binaries. In this review, we describe the dominant thermal, chemical, and radiation transfer processes in PDRs and XDRs, as well as give a brief description of models and their use for analyzing observations. We then present recent results from Milky Way, nearby extragalactic, and high-redshift observations. Several important results include the following: ▪ Velocity-resolved PDR lines reveal the kinematics of the neutral atomic gas and provide constraints on the stellar feedback process. Their interpretation is, however, in dispute, as observations suggest a prominent role for stellar winds, whereas they are much less important in theoretical models. ▪ A significant fraction of molecular mass resides in CO-dark gas especially in low-metallicity and/or highly irradiated environments. ▪ The CO ladder and [Ci][Formula: see text][Cii] ratios can determine if FUV or X rays dominate the ISM heating of extragalactic sources. ▪ With Atacama Large Millimeter/submillimeter Array, PDR and XDR tracers are now routinely detected on galactic scales over cosmic time. This makes it possible to link the star-formation history of the Universe to the evolution of the physical and chemical properties of the gas. Expected final online publication date for the Annual Review of Astronomy and Astrophysics Volume 60 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

The Low-redshift Lyman Continuum Survey. II. New Insights into LyC Diagnostics

Abstract: The Lyman continuum (LyC) cannot be observed at the epoch of reionization (z ≳ 6) owing to intergalactic H i absorption. To identify LyC emitters (LCEs) and infer the fraction of escaping LyC, astronomers have developed various indirect diagnostics of LyC escape. Using measurements of the LyC from the Low-redshift Lyman Continuum Survey (LzLCS), we present the first statistical test of these diagnostics. While optical depth indicators based on Lyα, such as peak velocity separation and equivalent width, perform well, we also find that other diagnostics, such as the [O iii]/[O ii] flux ratio and star formation rate surface density, predict whether a galaxy is an LCE. The relationship between these galaxy properties and the fraction of escaping LyC flux suggests that LyC escape depends strongly on H i column density, ionization parameter, and stellar feedback. We find that LCEs occupy a range of stellar masses, metallicities, star formation histories, and ionization parameters, which may indicate episodic and/or different physical causes of LyC escape.

The dynamics and outcome of star formation with jets, radiation, winds, and supernovae in concert

Abstract: We analyze the first giant molecular cloud (GMC) simulation to follow the formation of individual stars and their feedback from jets, radiation, winds, and supernovae, using the STARFORGE framework in the GIZMO code. We evolve the GMC for ∼ 9Myr, from initial turbulent collapse to dispersal by feedback. Protostellar jets dominate feedback momentum initially, but radiation and winds cause cloud disruption at ∼ 8% star formation efficiency (SFE), and the first supernova at 8 . 3Myr comes too late to influence star formation significantly. The per-freefall SFE is dynamic, accelerating from 0 to ∼ 18% before dropping quickly to <1%, but the estimate from YSO counts compresses it to a narrower range. The primary cluster forms hierarchically and condenses to a brief ( ∼ 1 Myr) compact ( ∼ 1pc) phase, but does not virialize before the cloud disperses, and the stars end as an unbound expanding association. The initial mass function resembles the Chabrier (2005) form with a high-mass slope 𝛼 = − 2 and a maximum mass of 55 𝑀 (cid:12) . Stellar accretion takes ∼ 400kyr on average, but (cid:38) 1Myr for > 10 𝑀 (cid:12) stars, so massive stars finish growing latest. The fraction of stars in multiples increases as a function of primary mass, as observed. Overall, the simulation much more closely resembles reality, compared to previous versions that neglected different feedback physics entirely. But more detailed comparison with synthetic observations will be needed to constrain the theoretical uncertainties.

OUP accepted manuscript

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, $\beta$, (b) the observed UV continuum flux at $1500$A, $F_{\rm UV}$, (c) the observed continuum flux at $\approx 158\mu$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 ($\beta, F_{\rm UV}, F_{158}$) values cannot be found. This occurs when the index $I_m = (F_{158}/F_{\rm UV})/(\beta-\beta_{\rm int})$, where $\beta_{\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-$\beta$ relation.

A New Census of the 0.2 &lt; z &lt; 3.0 Universe. II. The Star-forming Sequence

Abstract: Abstract We use the panchromatic spectral energy distribution (SED)-fitting code Prospector to measure the galaxy log M *–logSFR relationship (the star-forming sequence ) across 0.2 < z < 3.0 using the COSMOS-2015 and 3D-HST UV-IR photometric catalogs. We demonstrate that the chosen method of identifying star-forming galaxies introduces a systematic uncertainty in the inferred normalization and width of the star-forming sequence, peaking for massive galaxies at ∼0.5 and ∼0.2 dex, respectively. To avoid this systematic, we instead parameterize the density of the full galaxy population in the log M *–logSFR–redshift plane using a flexible neural network known as a normalizing flow. The resulting star-forming sequence has a low-mass slope near unity and a much flatter slope at higher masses, with a normalization 0.2–0.5 dex lower than typical inferences in the literature. We show this difference is due to the sophistication of the Prospector stellar populations modeling: the nonparametric star formation histories naturally produce higher masses while the combination of individualized metallicity, dust, and star formation history constraints produce lower star formation rates (SFRs) than typical UV+IR formulae. We introduce a simple formalism to understand the difference between SFRs inferred from SED fitting and standard template-based approaches such as UV+IR SFRs. Finally, we demonstrate the inferred star-forming sequence is consistent with predictions from theoretical models of galaxy formation, resolving a long-standing ∼ 0.2–0.5 dex offset with observations at 0.5 < z < 3. The fully trained normalizing flow including a nonparametric description of ρ ( log M * , logSFR , z ) is available online 20 20 https://github.com/jrleja/sfs_leja_trained_flow to facilitate straightforward comparisons with future work.

OUP accepted manuscript

Abstract: ABSTRACT We analyse the first giant molecular cloud (GMC) simulation to follow the formation of individual stars and their feedback from jets, radiation, winds, and supernovae, using the STARFORGE framework in the GIZMO code. We evolve the GMC for $\sim 9 \rm Myr$, from initial turbulent collapse to dispersal by feedback. Protostellar jets dominate feedback momentum initially, but radiation and winds cause cloud disruption at $\sim 8{{\ \rm per\ cent}}$ star formation efficiency (SFE), and the first supernova at $8.3\, \rm Myr$ comes too late to influence star formation significantly. The per-free-fall SFE is dynamic, accelerating from 0 per cent to $\sim 18{{\ \rm per\ cent}}$ before dropping quickly to &lt;1 per cent, but the estimate from YSO counts compresses it to a narrower range. The primary cluster forms hierarchically and condenses to a brief ($\sim 1\, \mathrm{Myr}$) compact ($\sim 1\, \rm pc$) phase, but does not virialize before the cloud disperses, and the stars end as an unbound expanding association. The initial mass function resembles the Chabrier (2005) form with a high-mass slope α = −2 and a maximum mass of 55 M⊙. Stellar accretion takes $\sim 400\, \rm kyr$ on average, but $\gtrsim 1\,\rm Myr$ for &gt;10 M⊙ stars, so massive stars finish growing latest. The fraction of stars in multiples increase as a function of primary mass, as observed. Overall, the simulation much more closely resembles reality, compared to previous versions that neglected different feedback physics entirely. But more detailed comparison with synthetic observations will be needed to constrain the theoretical uncertainties.

Toward a Better Understanding of Cosmic Chronometers: A New Measurement of H(z) at z ∼ 0.7

Abstract: We analyze the stellar ages obtained from a combination of Lick indices in Borghi et al. for 140 massive and passive galaxies selected in the LEGA-C survey at $0.6

The ALMA REBELS Survey: Specific Star-Formation Rates in the Reionization Era

Abstract: We present specific star-formation rates (sSFRs) for 40 UV-bright galaxies at z ∼ 7 − 8 observed as part of the Reionization Era Bright Emission Line Survey (REBELS) ALMA large program. The sSFRs are derived using improved SFR calibrations and SED-based stellar masses, made possible by measurements of far-infrared (FIR) continuum emission and [CII]-based spectroscopic redshifts. The median sSFR of the sample is $18_{-5}^{+7}~$Gyr−1, significantly larger than literature measurements lacking constraints in the FIR, reflecting the larger obscured SFRs derived from the dust continuum relative to that implied by the UV+optical SED. We suggest that such differences may reflect spatial variations in dust across these luminous galaxies, with the component dominating the FIR distinct from that dominating the UV. We demonstrate that the inferred stellar masses (and hence sSFRs) are strongly-dependent on the assumed star-formation history in reionization-era galaxies. When large sSFR galaxies (a population which is common at z > 6) are modeled with non-parametric star-formation histories, the derived stellar masses can increase by an order-of-magnitude relative to constant star-formation models, owing to the presence of a significant old stellar population that is outshined by the recent burst. The [CII] line widths in the largest sSFR systems are often very broad, suggesting dynamical masses capable of accommodating an old stellar population suggested by non-parametric models. Regardless of these systematic uncertainties among derived parameters, we find that sSFRs increase rapidly toward higher redshifts for massive galaxies (9.6 < log (M*/M⊙) < 9.8), evolving as (1 + z)1.7 ± 0.3, broadly consistent with expectations from the evolving baryon accretion rates.

The Structure of Multiphase Galactic Winds

Abstract: Abstract We present a novel analytic framework to model the steady-state structure of multiphase galactic winds comprised of a hot, volume-filling component and a cold, clumpy component. We first derive general expressions for the structure of the hot phase for arbitrary mass, momentum, and energy source terms. Next, informed by recent simulations, we parameterize the cloud–wind mass transfer rates, which are set by the competition between turbulent mixing and radiative cooling. This enables us to cast the cloud–wind interaction as a source term for the hot phase and thereby simultaneously solve for the evolution of both phases, fully accounting for their bidirectional influence. With this model, we explore the nature of galactic winds over a broad range of conditions. We find that (i) with realistic parameter choices, we naturally produce a hot, low-density wind that transports energy while entraining a significant flux of cold clouds, (ii) mixing dominates the cold cloud acceleration and decelerates the hot wind, (iii) during mixing thermalization of relative kinetic energy provides significant heating, (iv) systems with low hot phase mass loading factors and/or star formation rates can sustain higher initial cold phase mass loading factors, but the clouds are quickly shredded, and (v) systems with large hot phase mass loading factors and/or high star formation rates cannot sustain large initial cold phase mass loading factors, but the clouds tend to grow with distance from the galaxy. Our results highlight the necessity of accounting for the multiphase structure of galactic winds, both physically and observationally, and have important implications for feedback in galactic systems.

A survey of high-$z$ galaxies: SERRA simulations

Abstract: We introduce serra , a suite of zoom-in high-resolution (1 . 2 × 10 4 M (cid:12) , (cid:39) 25 pc at z = 7 . 7) cosmological simulations including non-equilibrium chemistry and on-the-fly radiative transfer. The outputs are post-processed to derive galaxy UV+FIR continuum and emission line properties. Results are compared with available multi-wavelength data to constrain the physical properties (e.g., star formation rates, stellar/gas/dust mass, metallicity) of high-redshift 6 < ∼ z < ∼ 15 galaxies. This flagship paper focuses on the z = 7 . 7 sub-sample, including 202 galaxies with stellar mass 10 7 M (cid:12) < ∼ M (cid:63) < ∼ 5 × 10 10 M (cid:12) , and specific star formation ranging from sSFR ∼ 100 Gyr − 1 in young, low-mass galaxies to ∼ 10 Gyr − 1 for older, massive ones. At this redshift, serra galaxies are typically bursty, i.e. they are located above the Schmidt-Kennicutt relation by a factor κ s = 3 . 03 +4 . 9 − 1 . 8 , consistent with recent findings for [O III ] and [C II ] emitters at high- z . They also show relatively large IRX = L FIR /L UV values as a result of their compact/clumpy morphology effectively blocking the stellar UV luminosity. Note that this conclusion might be affected by insufficient spatial resolution at the molecular cloud level. We confirm that early galaxies lie on the standard [C II ] − SFR relation; their observed L [OIII] /L [CII] (cid:39) 1 − 10 ratios can be reproduced by a part of the serra galaxies without the need of a top-heavy IMF and/or anomalous C/O abundances. [O I ] line intensities are similar to local ones, making ALMA high- z detections challenging but feasible ( ∼ 6 hr for an SFR of 50 M (cid:12) yr − 1 ).

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

Abstract: 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% of their star formation obscured; the total (UV+IR) star formation rates are in the range 31.5 < SFR/(M yr −1) < 129.5. The sampleaveraged dust mass and temperature are (1.3± 1.1)× 10M and 52± 11 K, respectively. However, in some galaxies dust is particularly abundant (REBELS-14, M ′ d ≈ 3.4 × 10M ), or hot (REBELS-18, Corresponding author: Andrea Ferrara andrea.ferrara@sns.it ar X iv :2 20 2. 07 66 6v 1 [ as tr oph .G A ] 1 5 Fe b 20 22