Search

JWST has enabled the detection of the UV continuum of galaxies at z>10, evidencing a population of extremely blue, potentially dust-free galaxies. Interpreting the UV spectra of galaxies as they redden is complicated by the well-known degeneracy between stellar ages, dust, and nebular continuum. The main goal of this paper is to develop a theoretical model for the relationship between galaxy UV slopes, bursty star formation histories, dust evolution, and the contribution from nebular regions. We accomplish this via cosmological zoom-in simulations, and in specific, build a layered model where we simulate the UV slopes of galaxies with increasingly complex physics. Our main results follow. (i) Unattenuated stellar populations with no nebular emission exhibit a diverse range of intrinsic UV slopes, with values ranging from beta ~ -3 --> -2.2 due to long delays between bursts. This is manifested by an inverse correlation between the intrinsic UV slope and sSFR for early galaxies such that higher sSFR corresponds to bluer UV slopes. (ii) When including dust, our model galaxies demonstrate a rapid rise in dust obscuration between z ~ 8-10. This increase in dust mass is due to high grain-grain shattering rates, and enhanced growth per unit dust mass in very small grains, resulting in UV-detected galaxies at z ~ 12 descending into ALMA-detectable galaxies by z ~ 6. The rapid rise in dust content at z ~ 8-10 leads to a systematic reddening of the UV slopes during this redshift range. (iii) The inclusion of nebular continuum reddens the UV slope by a median factor Delta beta ~ 0.2-0.4. However, when including nebular continuum, our highest redshift galaxies (z~12) are insufficiently blue compared to observations; this may imply an evolving escape fraction from HII regions with redshift., Comment: Submitted to ApJ; comments welcome

We use JWST/NIRSpec observations from the Assembly of Ultradeep Rest-optical Observations Revealing Astrophysics (AURORA) survey to constrain the shape of the nebular attenuation curve of a star-forming galaxy at z=4.41, GOODSN-17940. We utilize 11 unblended HI recombination lines to derive the attenuation curve spanning optical to near-infrared wavelengths (3751-9550 \r{A}). We then leverage a high-S/N spectroscopic detection of the rest-frame ultraviolet continuum in combination with rest-UV photometric measurements to constrain the shape of the curve at ultraviolet wavelengths. While this UV constraint is predominantly based on stellar emission, the large measured equivalent widths of H$\alpha$ and H$\beta$ indicate that GOODSN-17940 is dominated by an extremely young stellar population <10 Myr in age such that the UV stellar continuum experiences the same attenuation as the nebular emission. The resulting combined nebular attenuation curve spans 1400-9550 \r{A} and has a shape that deviates significantly from commonly assumed dust curves in high-redshift studies. Relative to the Milky Way, SMC, and Calzetti curves, the new curve has a steeper slope at long wavelengths ($\lambda>5000$ \r{A}) while displaying a similar slope across blue-optical wavelengths ($\lambda=3750-5000$ \r{A}). In the ultraviolet, the new curve is shallower than the SMC and Calzetti curves and displays no significant 2175 \r{A} bump. This work demonstrates that the most commonly assumed dust curves are not appropriate for all high-redshift galaxies. These results highlight the ability to derive nebular attenuation curves for individual high-redshift sources with deep JWST/NIRSpec spectroscopy, thereby improving the accuracy of physical properties inferred from nebular emission lines., Comment: 21 pages, 6 figures, 1 table, submitted to ApJ

We present results on the emission-line properties of z=1.4-7.5 star-forming galaxies in the Assembly of Ultradeep Rest-optical Observations Revealing Astrophysics (AURORA) Cycle 1 JWST/NIRSpec program. Based on its depth, continuous wavelength coverage from 1--5 microns, and medium spectral resolution (R~1000), AURORA includes detections of a large suite of nebular emission lines spanning a broad range in rest wavelength. We investigate the locations of AURORA galaxies in multiple different emission-line diagrams, including traditional "BPT" diagrams of [OIII]/Hbeta vs. [NII]/Halpha, [SII]/Halpha, and [OI]/Halpha, and the "ionization-metallicity" diagram of [OIII]/[OII] (O32) vs. ([OIII]+[OII])/Hbeta (R23). We also consider a bluer rest-frame "ionization-metallicity" diagram introduced recently to characterize z>10 galaxies: [NeIII]/[OII] vs. ([NeIII]+[OII])/Hdelta; as well as longer-wavelength diagnostic diagrams extending into the rest-frame near-IR: [OIII]/Hbeta vs. [SIII]/[SII] (S32); and HeI/Pagamma and [SIII]/Pagamma vs. [FeII]/Pabeta. With a significant boost in signal-to-noise and large, representative samples of individual galaxy detections, the AURORA emission-line diagrams presented here definitively confirm a physical picture in which chemically-young, alpha-enhanced, massive stars photoionize the ISM in distant galaxies with a harder ionizing spectrum at fixed nebular metallicity than in their z~0 counterparts. We also uncover previously unseen evolution prior to z~2 in the [OIII]/Hbeta vs. [NII]/Halpha diagram, which motivates deep NIRSpec observations at even higher redshift. Finally, we present the first statistical sample of rest-frame near-IR emission-line diagnostics in star-forming galaxies at high redshift. In order to truly interpret rest-frame near-IR line ratios including [FeII], we must obtain better constraints on dust depletion in the high-redshift ISM., Comment: 17 pages, 8 figures, submitted to ApJ

This paper presents an overview of the Milky Way and nearby galaxies science case for the \textit{Single Aperture Large Telescope for Universe Studies} (SALTUS) far-infrared NASA probe-class mission concept. SALTUS offers enormous gains in spatial resolution and spectral sensitivity over previous far-IR missions, thanks to its cold ($<$40~K) 14-m primary mirror. Key Milky Way and nearby galaxies science goals for SALTUS focus on understanding the role of star formation in feedback in the Local Universe. In addition to this science case, SALTUS would open a new window to to of Galactic and extragalactic communities in the 2030s, enable fundamentally new questions to be answered, and be a far-IR analog to the near- and mid-IR capabilities of JWST. This paper summarizes the Milky Way and nearby galaxies science case and plans for notional observing programs in both guaranteed and guest (open) time., Comment: Resubmitted to JATIS as part of a special issue on SALTUS

This paper presents an overview of the high-redshift extragalactic science case for the Single Aperture Large Telescope for Universe Studies (SALTUS) far-infrared NASA probe-class mission concept. Enabled by its 14m primary reflector, SALTUS offers enormous gains in spatial resolution and spectral sensitivity over previous far-IR missions. SALTUS would be a versatile observatory capable of responding to the scientific needs of the extragalactic community in the 2030s, and a natural follow-on to the near- and mid-IR capabilities of JWST. Key early-universe science goals for SALTUS focus on understanding the role of galactic feedback processes in regulating galaxy growth across cosmic time, and charting the rise of metals and dust from the early universe to the present. This paper summarizes these science cases and the performance metrics most relevant for high-redshift observations., Comment: Submitted to JATIS as part of a special issue on SALTUS. Key performance metrics in Figs. 1 and 2

We update the dust model present within the Simba galaxy simulations with a self-consistent framework for the co-evolution of dust and molecular hydrogen populations in the interstellar medium, and use this to explore $z \geq 6$ galaxy evolution. In addition to tracking the evolution of dust and molecular hydrogen abundances, our model fully integrates these species into the Simba simulation, explicitly modelling their impact on physical processes such as star formation and cooling through the inclusion of a novel two-phase sub-grid model for interstellar gas. In running two high-resolution simulations down to $z \sim 6$ we find that our Simba-EoR model displays a generally tighter concordance with observational data than fiducial Simba. Additionally we observe that our Simba-EoR models increase star formation activity at early epochs, producing larger dust-to-gas ratios consequently. Finally, we discover a significant population of hot dust at $\sim 100$ K, aligning with contemporaneous observations of high-redshift dusty galaxies, alongside the large $\sim 20$ K population typically identified., Comment: 19 pages, 12 figures, submitted to MNRAS, comments welcome

Using integrated spectra for two gravitationally lensed galaxies from the JWST TEMPLATES Early Release Science program, we analyze faint auroral lines, which provide direct measurements of the gas-phase chemical abundance. For the brighter galaxy, SGAS1723$+$34 ($z = 1.3293$), we detect the [OIII]$\lambda4363$, [SIII]$\lambda6312$, and [OII]$\lambda\lambda$7320,7330 auroral emission lines, and set an upper limit for the [NII]$\lambda5755$ line. For the second galaxy, SGAS1226$+$21 ($z = 2.925$), we do not detect any auroral lines, and report upper limits. With these measurements and upper limits, we constrain the electron temperatures in different ionization zones within both of these galaxies. For SGAS1723$+$34, where auroral lines are detected, we calculate direct oxygen and nitrogen abundances, finding an N/O ratio consistent with observations of nearby ($z\sim 0$) galaxies. These observations highlight the potent combination of JWST and gravitational lensing to measure faint emission lines in individual distant galaxies and to directly study the chemical abundance patterns in those galaxies., Comment: Submitted to the Astrophysical Journal's Focus Issue on the TEMPLATES JWST Early Release Science Program. 15 pages, 6 figures. Comments welcome

We explore the cosmic evolution of the fraction of dust obscured star formation predicted by the \textsc{simba} cosmological hydrodynamic simulations featuring an on-the-fly model for dust formation, evolution, and destruction. We find that up to $z=2$, our results are broadly consistent with previous observational results of little to no evolution in obscured star formation. However, at $z>2$ we find strong evolution at fixed galaxy stellar mass towards greater amounts of obscured star formation. We explain the trend of increasing obscuration at higher redshifts by greater typical dust column densities along the line of sight to young stars. We additionally see that at a fixed redshift, more massive galaxies have a higher fraction of their star formation obscured, which is explained by increased dust mass fractions at higher stellar masses. Finally, we estimate the contribution of dust-obscured star formation to the total star formation rate budget and find that the dust obscured star formation history (SFH) peaks around $z\sim 2-3$, and becomes subdominant at $z\gtrsim 5$., Comment: 15 pages, 9 figures, submitted to ApJ; comments welcome

This paper gives an overview of TEMPLATES, a JWST Early Release Science program that targeted four extremely bright, gravitationally lensed galaxies: two extremely dusty, two with low attenuation, as templates for galaxy evolution studies with JWST. TEMPLATES obtains a common set of spectral diagnostics for these 1.3 < z < 4.2 galaxies, in particular H alpha, Paschen alpha, and the rest-frame optical and near-infrared continua. In addition, two of the four targets have JWST coverage of [O III] 5007 Angstrom and H beta; the other two targets have have JWST coverage of PAH 3.3 micron and complementary ALMA data covering the [C II] 158 micron emission line. The science goals of TEMPLATES are to demonstrate attenuation-robust diagnostics of star formation, map the distribution of star formation, compare the young and old stellar populations, and measure the physical conditions of star formation and their spatial variation across the galaxies. In addition, TEMPLATES has technical goals to establish best practices for the Integral Field Units (IFU) within the NIRSpec and MIRI instruments, both in terms of observing strategy and in terms of data reduction. The paper describes TEMPLATES's observing program, scientific and technical goals, data reduction methods, and deliverables, including high-level data products and data reduction cookbooks., Comment: 28 pages, including 8 figures and 3 tables. Submitted to The Astrophysical Journal's Focus Issue on the TEMPLATES JWST Early Release Science Program, https://iopscience.iop.org/collections/apj-231204-01-templates

We present {\sc slick} (the Scalable Line Intensity Computation Kit), a software package that calculates realistic CO, [\ion{C}{1}], and [\ion{C}{2}] luminosities for clouds and galaxies formed in hydrodynamic simulations. Built on the radiative transfer code {\sc despotic}, {\sc slick} computes the thermal, radiative, and statistical equilibrium in concentric zones of model clouds, based on their physical properties and individual environments. We validate our results applying {\sc slick} to the high-resolution run of the {\sc Simba} simulations, testing the derived luminosities against empirical and theoretical/analytic relations. To simulate the line emission from a universe of emitting clouds, we have incorporated random forest machine learning (ML) methods into our approach, allowing us to predict cosmologically evolving properties of CO, [\ion{C}{1}] and [\ion{C}{2}] emission from galaxies such as luminosity functions. We tested this model in 100,000 gas particles, and 2,500 galaxies, reaching an average accuracy of $\sim$99.8\% for all lines. Finally, we present the first model light cones created with realistic and ML-predicted CO, [\ion{C}{1}], and [\ion{C}{2}] luminosities in cosmological hydrodynamical simulations, from $z=0$ to $z=10$., Comment: 17 pages, 11 figures, comments are welcome

The ratios of strong rest-frame optical emission lines are the dominant indicator of metallicities in high-redshift galaxies. Since typical strong-line based metallicity indicators are calibrated on auroral lines at $z=0$, their applicability for galaxies in the distant Universe is unclear. In this paper, we make use of mock emission line data from cosmological simulations to investigate the calibration of rest-frame optical emission lines as metallicity indicators at high redshift. Our model, which couples the SIMBA cosmological galaxy formation simulation with cloudy photoionization calculations, includes contributions from HII regions, post-AGB stars and Diffuse Ionized Gas (DIG). We find mild redshift evolution in the 12 indicators that we study, which implies that the dominant physical properties that evolve in our simulations do have a discernible impact on the metallicity calibrations at high redshifts. When comparing our calibrations with high redshift auroral line observations from James Webb Space Telescope we find a slight offset between our model results and the observations and find that a higher ionization parameter at high redshifts can be one of the possible explanations. We explore the physics that drives the shapes of strong-line metallicity relationships and propose calibrations for hitherto unexplored low-metallicity regimes. Finally, we study the contribution of DIG to total line fluxes. We find that the contribution of DIG increases with metallicity at z $\sim$ 0 for singly ionized oxygen and sulfur lines and can be as high as 70% making it crucial to include their contribution when modeling nebular emission., Comment: 28 pages, 6 figures, typo corrected: in polynomial fits, x represents metallicity and y represents the line ratio

Environmental effects on the evolution of galaxies have been one of the leading questions in galaxy studies for decades. In this work, we investigate the relationship between the star formation activity of galaxies and their environmental matter density using the cosmological hydrodynamic simulation Simba. The star formation activity indicators we explore include the star formation efficiency (SFE), specific star formation rate (sSFR) and molecular hydrogen mass fraction ($f^*_{H_2}$) and the environment is considered as the large-scale environmental matter density, calculated based on the stellar mass of nearby galaxies on a 1 Mpc/h grid using the cloud in cell (CIC) method. Our sample includes galaxies with $9<\log(M_*/M_{\odot})$ at $0 1.5$). We present a scaling relation for the depletion time of molecular hydrogen (${t_{depl}}=1/SFE$) as a function of galaxy parameters including environmental density. Our findings provide a framework for quantifying the environmental effects on the star formation activities of galaxies as a function of stellar mass and redshift. The most significant environmental dependence is seen at later cosmic times ($z<1.5$) and towards lower stellar masses ($9<\log(M_*/M_{\odot})<10$). Future large galaxy surveys can use this framework to look for the environmental dependence of the star formation activity and examine our predictions., Comment: 17 pages, 6 figures

We present JWST and ALMA results for the lensing system SPT0418-47, which includes a strongly-lensed, dusty star-forming galaxy at redshift z=4.225 and an associated multiply-imaged companion. JWST NIRCam and MIRI imaging observations presented in this paper were acquired as part of the Early Release Science program Targeting Extremely Magnified Panchromatic Lensed Arcs and Their Extended Star Formation (TEMPLATES). This data set provides robust, mutiwavelength detection of stellar light in both the main (SPT0418A) and companion (SPT0418B) galaxies, while the ALMA detection of [C II] emission confirms that SPT0418B lies at the same redshift as SPT0418A. From a source plane reconstruction, we infer that the projected physical separation of the two galaxies is $4.42\pm 0.05$ kpc. We derive total magnifications of $\mu=29.5\pm1.2$ and $\mu=4.2\pm 0.9$ for SPT0418A and SPT0418B, respectively. We use both CIGALE and PROSPECTOR to derive stellar masses. The stellar mass ratio of SPT0418A and SPT0418B is approximately 4 to 1 ($4.5\pm 1.0$ for CIGALE and $4.2^{+1.9}_{-1.6}$ for PROSPECTOR). We also see evidence of extended structure associated with SPT0418A in the lensing reconstruction that is suggestive of a tidal feature. Interestingly, the star formation rates and stellar masses of both galaxies are consistent with the main sequence of star-forming galaxies at this epoch, indicating that this ongoing interaction has not noticeably elevated the star formation levels., Comment: Submitted to AAS Journals,14 pages, 7 figures

We present the first spatially resolved maps of gas-phase metallicity for two dust-obscured star-forming galaxies (DSFGs) at $z\sim$ 4, from the JWST TEMPLATES Early Release Science program, derived from NIRSpec integral field unit spectroscopy of the H$\alpha$ and [NII] emission lines. Empirical optical line calibrations are used to determine that the sources are globally enriched to near-solar levels. While one source shows elevated [NII]/H$\alpha$ ratios and broad H$\alpha$ emission consistent with the presence of an AGN in a $\gtrsim$1kpc region, we argue that both systems have already undergone significant metal enrichment as a result of their extremely high star formation rates. Utilizing ALMA rest-frame 380$\mu$m continuum and [CI]($^3$P$_2$-$^3$P$_1$) line maps we compare the spatial variation of the metallicity and gas-to-dust ratio in the two galaxies, finding the two properties to be anticorrelated on highly resolved spatial scales, consistent with various literature studies of $z\sim$ 0 galaxies. The data are indicative of the enormous potential of JWST to probe the enrichment of the interstellar medium on $\sim$kpc scales in extremely dust-obscured systems at $z\sim$ 4 and beyond., Comment: 14 pages, 6 figures, accepted for publication in ApJ

In this paper, we demonstrate that the inference of galaxy stellar masses via spectral energy distribution (SED) fitting techniques for galaxies formed in the first billion years after the Big Bang carries fundamental uncertainties owing to the loss of star formation history (SFH) information from the very first episodes of star formation in the integrated spectra of galaxies. While this early star formation can contribute substantially to the total stellar mass of high-redshift systems, ongoing star formation at the time of detection outshines the residual light from earlier bursts, hampering the determination of accurate stellar masses. As a result, order of magnitude uncertainties in stellar masses can be expected. We demonstrate this potential problem via direct numerical simulation of galaxy formation in a cosmological context. In detail, we carry out two cosmological simulations with significantly different stellar feedback models which span a significant range in star formation history burstiness. We compute the mock SEDs for these model galaxies at z=7 via 3D dust radiative transfer calculations, and then backwards fit these SEDs with Prospector SED fitting software. The uncertainties in derived stellar masses that we find for z>7 galaxies motivate the development of new techniques and/or star formation history priors to model early Universe star formation., Comment: ApJ Accepted

In the current era of high-z galaxy discovery with JWST and ALMA, our ability to study the stellar populations and ISM conditions in a diverse range of galaxies at Cosmic Dawn has rapidly improved. At the same time, the need to understand the current limitations in modeling galaxy formation processes and physical properties in order to interpret these observations is critical. Here, we study the challenges in modeling galaxy dust temperatures, both in the context of forward modeling galaxy spectral properties from a hydrodynamical simulation and via backwards modeling galaxy physical properties from mock observations of far-infrared dust emission. Using the Simba model for galaxy formation combined with Powderday radiative transfer, we can accurately predict the evolution of dust at high redshift, though several aspects of the model are essentially free parameters (dust composition, sub-resolution dust in star-forming regions) that dull the predictive power of the model dust temperature distributions. We also highlight the uncertainties in the backwards modeling methods, where we find the commonly used models and assumptions to fit FIR SEDs and infer dust temperatures (e.g., single temperature, optically thin modified blackbody) largely fail to capture the complexity of high-z dusty galaxies. We caution that conclusions inferred from both simulations -- limited by resolution and post-processing techniques -- and observations -- limited by sparse data and simplistic model parameterizations -- are susceptible to unique and nuanced uncertainties that can limit the usefulness of current high-z dust measurements., Comment: 20 pages main text, 12 figures. Accepted to ApJ, newest version

Dust grains absorb half of the radiation emitted by stars throughout the history of the universe, re-emitting this energy at infrared wavelengths. Polycyclic aromatic hydrocarbons (PAHs) are large organic molecules that trace millimeter-size dust grains and regulate the cooling of the interstellar gas within galaxies. Observations of PAH features in very distant galaxies have been difficult due to the limited sensitivity and wavelength coverage of previous infrared telescopes. Here we present JWST observations that detect the 3.3um PAH feature in a galaxy observed less than 1.5 billion years after the Big Bang. The high equivalent width of the PAH feature indicates that star formation, rather than black hole accretion, dominates the infrared emission throughout the galaxy. The light from PAH molecules, large dust grains, and stars and hot dust are spatially distinct from one another, leading to order-of-magnitude variations in the PAH equivalent width and the ratio of PAH to total infrared luminosity across the galaxy. The spatial variations we observe suggest either a physical offset between the PAHs and large dust grains or wide variations in the local ultraviolet radiation field. Our observations demonstrate that differences in the emission from PAH molecules and large dust grains are a complex result of localized processes within early galaxies., Comment: Published in Nature 5 June 2023 at https://www.nature.com/articles/s41586-023-05998-6. MIRI MRS reduction notebook is available at https://github.com/jwst-templates

Quantifying changes in galaxies' interstellar medium (ISM) abundance after quenching star formation is an important aspect of galaxy evolution, but it is poorly constrained beyond the local universe. We characterise the dust-related properties in 548 quiescent galaxies observed at $0.1

We present a new methodology for simulating mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs) in galaxy evolution simulations. To do this, we combine theoretical models of PAH emission features as they respond to varying interstellar radiation fields, grain size distributions, and ionization states with a new on-the-fly model for dust evolution in hydrodynamic galaxy simulations. We apply these models to 3 idealized arepo galaxy evolution simulations within the smuggle physics framework. We use these simulations to develop numerical experiments investigating the buildup of PAH masses and luminosities in galaxies in idealized analogs of the Milky Way, a dwarf galaxy, and starburst disk. Our main results follow. Galaxies with high specific star formation rates have increased feedback energy per unit mass, and are able to efficiently shatter dust grains, driving up the fraction of ultra small grains. At the same time, in our model large radiation fields per unit gas density convert aliphatic grains into aromatics. The fraction of dust grains in the form of PAHs (q_PAH) can be understood as a consequence of these processes, and in our model PAHs form primarily from interstellar processing (shattering) of larger grains rather than from the growth of smaller grains. We find that the hardness of the radiation field plays a larger role than variations in the grain size distribution in setting the total integrated PAH luminosities, though cosmological simulations are necessary to fully investigate the complex interplay of processes that drive PAH band luminosities in galaxies. Finally, we highlight feature PAH strength variations, cautioning against the usage of emission templates with constant feature strength ratios., Comment: Submitted to ApJ; comments welcome

Observations of local star-forming galaxies (SFGs) show a tight correlation between their singly ionized carbon line luminosity ($L_{\rm [C_{II}]}$) and star formation rate (SFR), suggesting that $L_{\rm [C_{II}]}$ may be a useful SFR tracer for galaxies. Some other galaxy populations, however, are found to have lower $L_{\rm [C_{II}]}{}/{}\rm SFR$ than the local SFGs, including the infrared-luminous, starburst galaxies at low and high redshifts, as well as some moderately star-forming galaxies at the epoch of re-ionization (EoR). The origin of this `$\rm [C_{II}]$ deficit' is unclear. In this work, we study the $L_{\rm [C_{II}]}$-SFR relation of galaxies using a sample of $z=0-8$ galaxies with $M_*\approx10^7-5\times10^{11}\,M_\odot$ extracted from cosmological volume and zoom-in simulations from the Feedback in Realistic Environments (FIRE) project. We find a simple analytic expression for $L_{\rm [C_{II}]}$/SFR of galaxies in terms of the following parameters: mass fraction of $\rm [C_{II}]$-emitting gas ($f_{\rm [C_{II}]}$), gas metallicity ($Z_{\rm gas}$), gas density ($n_{\rm gas}$) and gas depletion time ($t_{\rm dep}{}={}M_{\rm gas}{}/{}\rm SFR$). We find two distinct physical regimes, where $t_{\rm dep}$ ($Z_{\rm gas}$) is the main driver of the $\rm [C_{II}]$ deficit in $\rm H_2$-rich ($\rm H_2$-poor) galaxies. The observed $\rm [C_{II}]$ deficit of IR-luminous galaxies and early EoR galaxies, corresponding to the two different regimes, is due to short gas depletion time and low gas metallicity, respectively. Our result indicates that $\rm [C_{II}]$ deficit is a common phenomenon of galaxies, and caution needs to be taken when applying a constant $L_{\rm [C_{II}]}$-to-SFR conversion factor derived from local SFGs to estimate cosmic SFR density at high redshifts and interpret data from upcoming $\rm [C_{II}]$ line intensity mapping experiments., Comment: 44 page, 24 figures, 8 tables and 8 appendices, accepted for publication in MNRAS. Key figures: fig. 16 & 17. Comments are welcome

Context. Deep spectroscopic surveys with the Atacama Large Millimeter/submillimeter Array (ALMA) have revealed that some of the brightest infrared sources in the sky correspond to concentrations of dusty star-forming galaxies (DSFG) at high redshift. Among these, the SPT2349-56 protocluster system at z = 4.304 is amongst the most extreme examples due to its high source density and integrated star formation rate. Aims. We conducted a deep Lyman-$\alpha$ line emission survey around SPT2349-56 using the Multi-Unit Spectroscopic Explorer (MUSE) at Very Large Telescope (VLT) in order to characterize this uniquely dense environment. Methods. Taking advantage of the deep three-dimensional nature of this survey, we performed a sensitive search for Lyman-$\alpha$ emitters (LAEs) toward the core and northern extension of the protocluster, which correspond to the brightest infrared regions in this field. Using a smoothed narrowband image extracted from the MUSE datacube around the protocluster redshift, we searched for possible extended structures. Results. We identify only three LAEs at z = 4.3 in this field, in concordance with expectations for blank-fields, and an extended Lyman-$\alpha$ structure spatially associated with core of the protocluster. All the previously-identified DSFGs in this field are undetected in Lyman-$\alpha$ emission, consistent with the conspicuous dust obscuration in these systems. We find an extended Lyman-$\alpha$ structure, about $60 \times 60$ kpc$^{2}$ in size, and located 56 kpc west of the protocluster core. Three DSFGs coincide spatially with the location of this structure. We conclude that either the three co-spatial DSFGs or the protocluster core itself are feeding ionizing photons to the Lyman-$\alpha$ structure., Comment: 14 pages, 9 figures, 2 tables, submitted to A&A

We present the Cosmic Sands suite of cosmological zoom-in simulations based on the Simba galaxy formation model in order to study the build up of the first massive and dusty galaxies in the early Universe. Residing in the most massive halos, we find that the compact proto-massive galaxies undergo nearly continuous mergers with smaller subhalos, boosting star formation rates (SFRs) and the build up of stellar mass. The galaxies are already appreciably chemically evolved by z=10, with modeled dust masses comparable to those inferred from observations in the same epoch. We track gas accretion onto the galaxies to understand how extreme SFRs can be sustained by these early systems. We find that smooth gas accretion can maintain SFRs above 250 M$_{\odot}$ / yr but to achieve SFRs that boost galaxies well above the main sequence, a larger perturbation like a gas-rich major merger is necessary to trigger a starburst episode. Post-processing the Cosmic Sands simulations with dust radiative transfer, we find that while the infrared luminosities of the most dust rich galaxies are comparable to local ULIRGs, they are substantially dimmer than classical z=2 sub-millimeter galaxies. We end with a discussion on the possible reasons for this discrepancy at the highest masses and the future work we intend to carry out to study the chemical enrichment of the earliest dusty galaxies., Comment: Submitted to ApJ, main text 17 pages, 12 figures. Comments welcome!

We present visual classifications of merger-induced tidal disturbances in 143 $\rm{M}_* \sim 10^{11}\rm{M}_\odot$ post-starburst galaxies at z$\sim$0.7 identified in the SQuIGG$\vec{L}$E Sample. This sample spectroscopically selects galaxies from the Sloan Digital Sky Survey that have stopped their primary epoch of star formation within the past $\sim$500 Myrs. Visual classifications are performed on Hyper Suprime Cam (HSC) i-band imaging. We compare to a control sample of mass- and redshift-matched star-forming and quiescent galaxies from the Large Early Galaxy Census and find that post-starburst galaxies are more likely to be classified as disturbed than either category. This corresponds to a factor of $3.6^{+2.9}_{-1.3}$ times the disturbance rate of older quiescent galaxies and $2.1^{+1.9}_{-.73}$ times the disturbance rate of star-forming galaxies. Assuming tidal features persist for $\lesssim500$ Myr, this suggests merging is coincident with quenching in a significant fraction of these post-starbursts. Galaxies with tidal disturbances are younger on average than undisturbed post-starburst galaxies in our sample, suggesting tidal features from a major merger may have faded over time. This may be exacerbated by the fact that, on average, the undisturbed subset is fainter, rendering low surface brightness tidal features harder to identify. However, the presence of ten young ($\lesssim150$ Myr since quenching) undisturbed galaxies suggests that major mergers are not the only fast physical mechanism that shut down the primary epoch of star formation in massive galaxies at intermediate redshift., Comment: 13 pages, 8 figures, accepted to ApJ

We present ALMA sub-kpc- to kpc-scale resolution observations of the [CII], CO(9-8), and OH$^{+}$\,($1_{1}$--$0_{1}$) lines along with their dust continuum emission toward the FIR luminous quasar SDSS J231038.88+185519.7 at $z = 6.0031$. The [CII] brightness follows a flat distribution with a Sersic index of 0.59. The CO(9-8) line and the dust continuum can be fit with an unresolved nuclear component and an extended Sersic component with a Sersic index of ~1. The dust temperature drops with distance from the center. The effective radius of the dust continuum is smaller than that of the line emission and the dust mass surface density, but is consistent with that of the star formation rate surface density. The OH$^{+}$\,($1_{1}$--$0_{1}$) line shows a P-Cygni profile with an absorption, which may indicate an outflow with a neutral gas mass of $(6.2\pm1.2)\times10^{8} M_{\odot}$ along the line of sight. We employed a 3D tilted ring model to fit the [CII] and CO(9-8) data cubes. The two lines are both rotation dominated and trace identical disk geometries and gas motions. We decompose the circular rotation curve measured from the kinematic model fit to the [CII] line into four matter components (black hole, stars, gas, and dark matter). The quasar-starburst system is dominated by baryonic matter inside the central few kiloparsecs. We constrain the black hole mass to be $2.97^{+0.51}_{-0.77}\times 10^{9}\,M_{\odot}$; this is the first time that the dynamical mass of a black hole has been measured at $z\sim6$. A massive stellar component (on the order of $10^{9}\,M_{\odot}$) may have already existed when the Universe was only ~0.93 Gyr old. The relations between the black hole mass and the baryonic mass of this quasar indicate that the central supermassive black hole may have formed before its host galaxy. [Abridged version. Please see the full abstract in the manuscript.], Comment: Accepted for publication in A&A

One of the most fundamental baryonic matter components of galaxies is the neutral atomic hydrogen (H$\,{\rm \scriptsize I}$). At low redshifts, this component can be traced directly through the 21-cm transition, but to infer H$\,{\rm \scriptsize I}$ gas content of the most distant galaxies, a viable tracer is needed. We here investigate the fidelity of the fine structure transition of the ($^2P_{3/2} - ^2P_{1/3}$) transition of singly-ionized carbon [C$\,{\rm \scriptsize II}$] at $158\,\mu$m as a proxy for H$\,{\rm \scriptsize I}$ in a set simulated galaxies at $z\approx 6$, following the work by Heintz et al. (2021). We select 11,125 star-forming galaxies from the SIMBA simulations, with far-infrared line emissions post-processed and modeled within the SIGAME framework. We find a strong connection between [C$\,{\rm \scriptsize II}$] and H$\,{\rm \scriptsize I}$, with the relation between this [C$\,{\rm \scriptsize II}$]-to-H$\,{\rm \scriptsize I}$ relation ($\beta_{\rm [C\,{\rm \scriptsize II}]}$) being anti-correlated with the gas-phase metallicity of the simulated galaxies. We further use these simulations to make predictions for the total baryonic matter content of galaxies at $z\approx 6$, and specifically the HI gas mass fraction. We find mean values of $M_{\rm HI}/M_\star = 1.4$, and $M_{\rm HI}/M_{\rm bar,tot} = 0.45$. These results provide strong evidence for H$\,{\rm \scriptsize I}$ being the dominant baryonic matter component by mass in galaxies at $z\approx 6$., Comment: 7 pages, 3 figures. Accepted for publication by ApJL

Observations and simulations have demonstrated that star formation in galaxies must be actively suppressed to prevent the formation of over-massive galaxies. Galactic outflows driven by stellar feedback or supermassive black hole accretion are often invoked to regulate the amount of cold molecular gas available for future star formation, but may not be the only relevant quenching processes in all galaxies. We present the discovery of vast molecular tidal features extending up to 64 kpc outside of a massive z=0.646 post-starburst galaxy that recently concluded its primary star-forming episode. The tidal tails contain (1.2 +/- 0.1)x10^10 Msun of molecular gas, 47 +/- 5 % of the total cold gas reservoir of the system. Both the scale and magnitude of the molecular tidal features are unprecedented compared to all known nearby or high-redshift merging systems. We infer that the cold gas was stripped from the host galaxies during the merger, which is most likely responsible for triggering the initial burst phase and the subsequent suppression of star formation. While only a single example, this result shows that galaxy mergers can regulate the cold gas contents in distant galaxies by directly removing a large fraction of the molecular gas fuel, and plausibly suppress star formation directly, a qualitatively different physical mechanism than feedback-driven outflows., Comment: Published in ApJL, 8 pages, 5 figures

Accurate models of the star formation histories (SFHs) of recently-quenched galaxies can provide constraints on when and how galaxies shut down their star formation. The recent development of "non-parametric" SFH models promises the flexibility required to make these measurements. However, model and prior choices significantly affect derived SFHs, particularly for post-starburst galaxies (PSBs) which have sharp changes in their recent SFH. In this paper, we create mock PSBs, then use the Prospector SED fitting software to test how well four different SFH models recover key properties. We find that a two-component parametric model performs well for our simple mock galaxies, but is sensitive to model mismatches. The fixed- and flexible-bin non-parametric models included in Prospector are able to rapidly quench a major burst of star formation, but systematically underestimate the post-burst age by up to 200 Myr. We develop a custom SFH model that allows for additional flexibility in the recent SFH. Our flexible non-parametric model is able to constrain post-burst ages with no significant offset and just ~90 Myr of scatter. Our results suggest that while standard non-parametric models are able to recover first-order quantities of the SFH (mass, SFR, average age), accurately recovering higher-order quantities (burst fraction, quenching time) requires careful consideration of model flexibility. These mock recovery tests are a critical part of future SFH studies. Finally, we show that our new, public SFH model is able to accurately recover the properties of mock star-forming and quiescent galaxies and is suitable for broader use in the SED fitting community., Comment: 24 pages, 10 figures. Posted after referee report; comments welcome. Key results in Figures 8 & 9. New SFH model publicly available as "continuity_psb_sfh" in Prospector: https://github.com/bd-j/prospector

The molecular gas in galaxies traces both the fuel for star formation and the processes that can enhance or suppress star formation. Observations of the molecular gas state can thus point to when and why galaxies stop forming stars. In this study, we present ALMA observations of the molecular gas in galaxies evolving through the post-starburst phase. These galaxies have low current star formation rates, regardless of the SFR tracer used, with recent starbursts ending within the last 600 Myr. We present CO (3-2) observations for three post-starburst galaxies, and dense gas HCN/HCO+/HNC (1-0) observations for six (four new) post-starburst galaxies. The post-starbursts have low excitation traced by the CO spectral line energy distribution (SLED) up to CO (3-2), more similar to early-type than starburst galaxies. The low excitation indicates that lower density rather than high temperatures may suppress star formation during the post-starburst phase. One galaxy displays a blueshifted outflow traced by CO (3-2). MaNGA observations show that the ionized gas velocity is disturbed relative to the stellar velocity field, with a blueshifted component aligned with the molecular gas outflow, suggestive of a multiphase outflow. Low ratios of HCO+/CO, indicating low fractions of dense molecular gas relative to the total molecular gas, are seen throughout post-starburst phase, except for the youngest post-starburst galaxy considered here. These observations indicate that the impact of any feedback or quenching processes may be limited to low excitation and weak outflows in the cold molecular gas during the post-starburst phase., Comment: 33 pages, 24 figures, 4 tables, ApJ in press

We present a high-resolution study of the cold molecular gas as traced by CO(1-0) in the unlensed z$\sim$3.4 submillimeter galaxy SMM J13120+4242, using multi-configuration observations with the Karl G. Jansky Very Large Array (JVLA). The gas reservoir, imaged on 0.39" ($\sim$3 kpc) scales, is resolved into two components separated by $\sim$11 kpc with a total extent of 16 $\sim$3 kpc. Despite the large spatial extent of the reservoir, the observations show a CO(1-0) FWHM linewidth of only 267 $\pm$ 64 km s$^{-1}$. We derive a revised line luminosity of L'$_\mathrm{CO(1-0)}$ = (10 $\pm$ 3) $\times$ 10$^{10}$ K km s$^{-1}$ pc$^2$ and a molecular gas mass of M$_\mathrm{gas}$ = (13 $\pm$ 3) $\times$ 10$^{10}$ ($\alpha_\mathrm{CO}$/1) M$_{\odot}$. Despite the presence of a velocity gradient (consistent with previous resolved CO(6-5) imaging), the CO(1-0) imaging shows evidence for significant turbulent motions which are preventing the gas from fully settling into a disk. The system likely represents a merger in an advanced stage. Although the dynamical mass is highly uncertain, we use it to place an upper limit on the CO-to-H$_2$ mass conversion factor $\alpha_\mathrm{CO}$ of 1.4. We revisit the SED fitting, finding that this galaxy lies on the very massive end of the main sequence at z = 3.4. Based on the low gas fraction, short gas depletion time and evidence for a central AGN, we propose that SMM J13120 is in a rapid transitional phase between a merger-driven starburst and an unobscured quasar. The case of SMM J13120 highlights the how mergers may drive important physical changes in galaxies without pushing them off the main sequence., Comment: 19 pages, 12 figures, accepted for publication in ApJ

The SPT0311-58 system at z=6.900 is an extremely massive structure within the reionization epoch, and offers a chance to understand the formation of galaxies in an extreme peak in the primordial density field. We present 70mas Atacama Large Millimeter/submillimeter Array observations of the dust continuum and CII 158um emission in the central pair of galaxies and reach physical resolution ~100-350pc, among the most detailed views of any reionization-era system to date. The observations resolve the source into at least a dozen kiloparsec-size clumps. The global kinematics and high turbulent velocity dispersion within the galaxies present a striking contrast to recent claims of dynamically cold thin-disk kinematics in some dusty galaxies just 800Myr later at z~4. We speculate that both gravitational interactions and fragmentation from massive parent disks have likely played a role in the overall dynamics and formation of clumps in the system. Each clump individually is comparable in mass to other 6

We present structural measurements of 145 spectroscopically selected intermediate-redshift (z$\sim$0.7), massive ($M_\star \sim 10^{11} \ M_\odot$) post-starburst galaxies from the SQuIGG$\vec{L}$E Sample measured using wide-depth Hyper Suprime-Cam i-band imaging. This deep imaging allows us to probe the sizes and structures of these galaxies, which we compare to a control sample of star forming and quiescent galaxies drawn from the LEGA-C Survey. We find that post-starburst galaxies systematically lie $\sim0.1$ dex below the quiescent mass-size (half-light radius) relation, with a scatter of $\sim0.2$ dex. This finding is bolstered by non-parametric measures, such as the Gini coefficient and the concentration, which also reveal these galaxies to have more compact light profiles than both quiescent and star-forming populations at similar mass and redshift. The sizes of post-starburst galaxies show either negative or no correlation with the time since quenching, such that more recently quenched galaxies are larger or similarly sized. This empirical finding disfavors the formation of post-starburst galaxies via a purely central burst of star formation that simultaneously shrinks the galaxy and shuts off star formation. We show that the central densities of post-starburst and quiescent galaxies at this epoch are very similar, in contrast with their effective radii. The structural properties of z$\sim$0.7 post-starburst galaxies match those of quiescent galaxies that formed in the early universe, suggesting that rapid quenching in the present epoch is driven by a similar mechanism to the one at high redshift., Comment: 28 pages, 19 figures, 2 appendices, updated to accepted ApJ version

We investigate the fine-structure [C${\rm \scriptsize II}$] line at $158\,\mu$m as a molecular gas tracer by analyzing the relationship between molecular gas mass ($M_{\rm mol}$) and [C${\rm \scriptsize II}$] line luminosity ($L_{\rm [CII]}$) in 11,125 $z\simeq 6$ star-forming, main sequence galaxies from the SIMBA simulations, with line emission modeled by S\'IGAME. Though most ($\sim 50-100\,\%$) of the gas mass in our simulations is ionized, the bulk ($> 50\,\%$) of the [C${\rm \scriptsize II}$] emission comes from the molecular phase. We find a sub-linear (slope $0.78\pm 0.01$) $\log L_{\rm [CII]}-\log M_{\rm mol}$ relation, in contrast with the linear relation derived from observational samples of more massive, metal-rich galaxies at $z \lesssim 6$. We derive a median [C${\rm \scriptsize II}$]-to-$M_{\rm mol}$ conversion factor of $\alpha_{\rm [CII]} \simeq 18\,{\rm M_{\rm \odot}/L_{\rm \odot}}$. This is lower than the average value of $\simeq 30\,{\rm M_{\rm \odot}/L_{\rm \odot}}$ derived from observations, which we attribute to lower gas-phase metallicities in our simulations. Thus, a lower, luminosity-dependent, conversion factor must be applied when inferring molecular gas masses from [C${\rm \scriptsize II}$] observations of low-mass galaxies. For our simulations, [C${\rm \scriptsize II}$] is a better tracer of the molecular gas than CO $J=1-0$, especially at the lowest metallicities, where much of the gas is 'CO-dark'. We find that $L_{\rm [CII]}$ is more tightly correlated with $M_{\rm mol}$ than with star-formation rate (${\rm SFR}$), and both the $\log L_{\rm [CII]}-\log M_{\rm mol}$ and $\log L_{\rm [CII]}-\log {\rm SFR}$ relations arise from the Kennicutt-Schmidt relation. Our findings suggest that $L_{\rm [CII]}$ is a promising tracer of the molecular gas at the earliest cosmic epochs., Comment: 13 pages, 9 figures. Accepted for publication in ApJ

One of the most common methods for inferring galaxy attenuation curves is via spectral energy distribution (SED) modeling, where the dust attenuation properties are modeled simultaneously with other galaxy physical properties. In this paper, we assess the ability of SED modeling to infer these dust attenuation curves from broadband photometry, and suggest a new flexible model that greatly improves the accuracy of attenuation curve derivations. To do this, we fit mock SEDs generated from the Simba cosmological simulation with the Prospector SED fitting code. We consider the impact of the commonly-assumed uniform screen model and introduce a new non-uniform screen model parameterized by the fraction of unobscured stellar light. This non-uniform screen model allows for a non-zero fraction of stellar light to remain unattenuated, resulting in a more flexible attenuation curve shape by decoupling the shape of the UV attenuation curve from the optical attenuation curve. The ability to constrain the dust attenuation curve is significantly improved with the use of a non-uniform screen model, with the median offset in UV attenuation decreasing from $-0.30$ dex with a uniform screen model to $-0.17$ dex with the non-uniform screen model. With this increase in dust attenuation modeling accuracy, we also improve the star formation rates (SFRs) inferred with the non-uniform screen model, decreasing the SFR offset on average by $0.12$ dex. We discuss the efficacy of this new model, focusing on caveats with modeling star-dust geometries and the constraining power of available SED observations., Comment: Submitted to ApJ. 15 pages, 10 figures

SPT0346-52 (z=5.7) is the most intensely star-forming galaxy discovered by the South Pole Telescope, with Sigma_SFR ~ 4200 Msol yr^-1 kpc^-2. In this paper, we expand on previous spatially-resolved studies, using ALMA observations of dust continuum, [NII]205 micron, [CII]158 micron, [OI]146 micron, and undetected [NII]122 micron and [OI]63 micron emission to study the multi-phase interstellar medium (ISM) in SPT0346-52. We use pixelated, visibility-based lens modeling to reconstruct the source-plane emission. We also model the source-plane emission using the photoionization code CLOUDY and find a supersolar metallicity system. We calculate T_dust = 48.3 K and lambda_peak = 80 micron, and see line deficits in all five lines. The ionized gas is less dense than comparable galaxies, with n_e < 32 cm^-3, while ~20% of the [CII]158 emission originates from the ionized phase of the ISM. We also calculate the masses of several phases of the ISM. We find that molecular gas dominates the mass of the ISM in SPT0346-52, with the molecular gas mass ~4x higher than the neutral atomic gas mass and ~100x higher than the ionized gas mass., Comment: 25 pages, 15 figures. Accepted for publication in ApJ

There is now ample evidence that dust is already present in abundance at high z. However, given the faintness of distant galaxies in the optical and the NIR, datasets are still limited and how the dust affects the emerging radiation of galaxies at very high redshift is not yet fully understood. Using the ALPINE survey, our objective is to quantify the dust attenuation properties in galaxies at z=4.4-5.5, and in particular the shape of their attenuation curve. Using the CIGALE code, we model the stellar populations and their interaction with the dust in order to measure some of the physical properties of a subsample of 23 main-sequence ALPINE galaxies. We find that the attenuation curves span a broad range of properties, from curves that are much steeper than the SMC extinction curve, to shallower than the starburst attenuation curve. The shape of the attenuation curves strongly depends on the V-band attenuation. Galaxies with the lowest attenuation also present the steepest curves. The steepness of such curves is probably the consequence of the combination of the intrinsic physical properties of the dust, the relative distribution of stars and dust in the interstellar medium, and the differential reddening. The broad range of attenuation curves found at z~5 shows that no single attenuation curve is appropriate for main sequence galaxies and that assuming a fixed curve can lead to large errors, for instance in the interpretation and use of the IRX-beta diagram, if SED modeling is not feasible. Great caution should be exercised when correcting high redshift galaxies for the presence of dust using the UV slope beta as it can affect the estimation of both SFR and stellar mass even at low V-band attenuation due to the steepness of the attenuation curve. However, when SED modeling can be used, the impact of the choice of the attenuation curve on the SFR and the stellar mass is limited., Comment: 20 pages, 14 figures. Accepted for publication in A&A

The Baldwin, Philips, & Terlevich diagram of [O III]/H$\beta$ vs. [N II]/H$\alpha$ (hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies at $z\sim2$ reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understand the potential physical processes driving this offset. We model nebular line emission in a large sample of galaxies, taken from the SIMBA cosmological hydrodynamic galaxy formation simulation, using the CLOUDY photoionization code to compute the nebular line luminosities from H II regions. We find that the observed shift toward higher [O III]/H$\beta$ and [N II]/H$\alpha$ values at high redshift arises from sample selection: when we consider only the most massive galaxies $M_* \sim 10^{10-11} M_\odot$, the offset naturally appears, due to their high metallicities. We predict that deeper observations that probe lower-mass galaxies will reveal galaxies that lie on a locus comparable to $z\sim 0$ observations. Even when accounting for sample selection effects, we find that there is a subtle mismatch between simulations and observations. To resolve this discrepancy, we investigate the impact of varying ionization parameters, H II region densities, gas-phase abundance patterns, and increasing radiation field hardness on N2-BPT diagrams. We find that either decreasing the ionization parameter or increasing the N/O ratio of galaxies at fixed O/H can move galaxies along a self-similar arc in N2-BPT space that is occupied by high-redshift galaxies., Comment: 25 pages, 15 figures. Accepted for publication in ApJ

We report the serendipitous discovery of an overdensity of CO emitters in an X-ray-identified cluster (Log$_{10}$M$_{\rm halo}/M_{\odot}\sim13.6$ at z=1.3188) using ALMA. We present spectroscopic confirmation of 6 new cluster members exhibiting CO(2-1) emission, adding to 2 existing optical/IR spectroscopic members undetected in CO. This is the lowest mass cluster to date at z>1 with molecular gas measurements, bridging the observational gap between galaxies in the more extreme, well-studied clusters (Log$_{10}$~M$_{\rm halo}/M_{\odot}\gtrsim14$) and those in group or field environments at cosmic noon. The CO sources are concentrated on the sky (within ~1-arcmin diameter) and phase space analysis indicates the gas resides in galaxies already within the cluster environment. We find that CO sources sit in similar phase space as CO-rich galaxies in more massive clusters at similar redshifts (have similar accretion histories) while maintaining field-like molecular gas reservoirs, compared to scaling relations. This work presents the deepest CO survey to date in a galaxy cluster at z>1, uncovering gas reservoirs down to M$_{\rm H_{2}}>1.6\times10^{10}$M$_{\odot}$ (5$\sigma$ at 50% primary beam). Our deep limits rule out the presence of gas content in excess of the field scaling relations; however, combined with literature CO detections, cluster gas fractions in general appear systematically high, on the upper envelope or above the field. This study is the first demonstration that low mass clusters at z~1-2 can host overdensities of CO emitters with surviving gas reservoirs, in line with the prediction that quenching is delayed after first infall while galaxies consume the gas bound to the disk., Comment: Submitted to AAS Journals on December 15, 2021. Comments welcome

The Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4,233 cosmological simulations, 2,049 N-body and 2,184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper we present the CAMELS public data release, describing the characteristics of the CAMELS simulations and a variety of data products generated from them, including halo, subhalo, galaxy, and void catalogues, power spectra, bispectra, Lyman-$\alpha$ spectra, probability distribution functions, halo radial profiles, and X-rays photon lists. We also release over one thousand catalogues that contain billions of galaxies from CAMELS-SAM: a large collection of N-body simulations that have been combined with the Santa Cruz Semi-Analytic Model. We release all the data, comprising more than 350 terabytes and containing 143,922 snapshots, millions of halos, galaxies and summary statistics. We provide further technical details on how to access, download, read, and process the data at \url{https://camels.readthedocs.io}., Comment: 18 pages, 3 figures. More than 350 Tb of data from thousands of simulations publicly available at https://www.camel-simulations.org

We describe the SQuIGGLE survey of intermediate-redshift post-starburst galaxies. We leverage the large sky coverage of the SDSS to select ~1300 recently-quenched galaxies at 0.575% of SQuIGGLE galaxies formed at least a quarter of their total stellar mass in the recent burst, which ended just ~200Myr before observation. We find that SQuIGGLE galaxies are on average younger and more burst-dominated than most other z<~1 post-starburst samples. This large sample of bright post-starburst galaxies at intermediate redshift opens a wide range of studies into the quenching process. In particular, the full SQuIGGLE survey will investigate the molecular gas reservoirs, morphologies, kinematics, resolved stellar populations, AGN incidence, and infrared properties of this unique sample of galaxies in order to place definitive constraints on the quenching process., Comment: 23 pages, 16 figures, accepted to ApJ

We use ALMA observations of CO(2-1) in 13 massive ($M_{\star}\gtrsim 10^{11} M_{\odot}$) post-starburst galaxies at $z\sim0.6$ to constrain the molecular gas content in galaxies shortly after they quench their major star-forming episode. The post-starburst galaxies in this study are selected from the Sloan Digital Sky Survey spectroscopic samples (DR14) based on their spectral shapes, as part of the SQuIGGLE program. Early results showed that two post-starburst galaxies host large H$_2$ reservoirs despite their low inferred star formation rates. Here we expand this analysis to a larger statistical sample of 13 galaxies. Six of the primary targets (45%) are detected, with $M_{H_2}\gtrsim10^9 M_{\odot}$. Given their high stellar masses, this mass limit corresponds to an average gas fraction of $\langle f_{H_2} \equiv M_{H_2}/M_{\star} \rangle \sim7\%$, or ${\sim}14\%$ using lower stellar masses estimates derived from analytic, exponentially declining star formation histories. The gas fraction correlates with the $D_n4000$ spectral index, suggesting that the cold gas reservoirs decrease with time since burst, as found in local K+A galaxies. Star formation histories derived from flexible stellar population synthesis modeling support this empirical finding: galaxies that quenched $\lesssim 150$ Myr prior to observation host detectable CO(2-1) emission, while older post-starburst galaxies are undetected. The large $\mathrm{H_2}$ reservoirs and low star formation rates in the sample imply that the quenching of star formation precedes the disappearance of the cold gas reservoirs. However, within the following 100-200 Myrs, the SQuIGGLE galaxies require the additional and efficient heating or removal of cold gas to bring their low star formation rates in line with standard $\mathrm{H_2}$ scaling relations., Comment: 20 pages, 13 figures, accepted for publication in ApJ

We present new constraints on the masses of the halos hosting the Milky Way and Andromeda galaxies derived using graph neural networks. Our models, trained on thousands of state-of-the-art hydrodynamic simulations of the CAMELS project, only make use of the positions, velocities and stellar masses of the galaxies belonging to the halos, and are able to perform likelihood-free inference on halo masses while accounting for both cosmological and astrophysical uncertainties. Our constraints are in agreement with estimates from other traditional methods., Comment: 2 figures, 2 tables, 7 pages. Code publicly available at https://github.com/PabloVD/HaloGraphNet

Understanding the halo-galaxy connection is fundamental in order to improve our knowledge on the nature and properties of dark matter. In this work we build a model that infers the mass of a halo given the positions, velocities, stellar masses, and radii of the galaxies it hosts. In order to capture information from correlations among galaxy properties and their phase-space, we use Graph Neural Networks (GNNs), that are designed to work with irregular and sparse data. We train our models on galaxies from more than 2,000 state-of-the-art simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project. Our model, that accounts for cosmological and astrophysical uncertainties, is able to constrain the masses of the halos with a $\sim$0.2 dex accuracy. Furthermore, a GNN trained on a suite of simulations is able to preserve part of its accuracy when tested on simulations run with a different code that utilizes a distinct subgrid physics model, showing the robustness of our method. The PyTorch Geometric implementation of the GNN is publicly available on Github at https://github.com/PabloVD/HaloGraphNet, Comment: 20 pages, 8 figures, code publicly available at https://github.com/PabloVD/HaloGraphNet

Observations of cold molecular gas reservoirs are critical for understanding the shutdown of star formation in massive galaxies. While dust continuum is an efficient and affordable tracer, this method relies upon the assumption of a "normal" molecular-gas to dust mass ratio, $\delta_{\mathrm{GDR}}$, typically of order one hundred. Recent null detections of quiescent galaxies in deep dust continuum observations support a picture where the cold gas and dust has been rapidly depleted or expelled. In this work, we present another viable explanation: a significant fraction of galaxies with low star formation per unit stellar mass are predicted to have extreme $\delta_{\mathrm{GDR}}$ ratios. We show that simulated massive quiescent galaxies at $0 < z < 3$ in the \textsc{simba} cosmological simulations have $\delta_{\mathrm{GDR}}$ values that extend $>$4 orders of magnitude. The dust in most simulated quiescent galaxies is destroyed significantly more rapidly than the molecular gas depletes, and cannot be replenished. The transition from star-forming to quiescent halts dust formation via star formation processes, with dust subsequently destroyed by supernova shocks and thermal sputtering of dust grains embedded in hot plasma. After this point, the dust growth rate in the models is not sufficient to overcome the loss of $>$3 orders of magnitude in dust mass to return to normal values of $\delta_{\mathrm{GDR}}$ despite having high metallicity. Our results indicate that it is not straight forward to use a single observational indicator to robustly pre-select exotic versus normal ratios. These simulations make strong predictions that can be tested with millimeter facilities., Comment: 9 pages, 5 figures, version accepted for publication in the Astrophysical Journal Letters

Detecting the line-intensity mapping (LIM) signal from the galaxies of the Epoch of Reionization is an emerging tool to constrain their role in reionization. Ongoing and upcoming experiments target the signal fluctuations across the sky to reveal statistical and astrophysical properties of these galaxies via signal statistics, e.g., the power spectrum. Here, we revisit the [C II]$_{158 \mu \text{m}}$ LIM power spectrum under non-uniform line-luminosity scatter, which has a halo-mass variation of statistical properties. Line-luminosity scatter from a cosmological hydrodynamic and radiative transfer simulation of galaxies at $z=6$ is considered in this study. We test the robustness of different model frameworks that interpret the impact of the line-luminosity scatter on the signal statistics. We use a simple power-law model to fit the scatter and demonstrate that the mean luminosity-halo mass correlation fit cannot preserve the mean intensity of the LIM signal (hence the clustering power spectrum) under non-uniform scatter. In our case, the mean intensity changes by $\sim 48$ per cent compared to the mean correlation fit in contrast to the general case with semi-analytic scatter. However, we find that the prediction for the mean intensity from the most-probable fit can be modelled robustly, considering the generalized and more realistic non-uniform scatter. We also explore the possibility of diminishing luminosity bias under non-uniform scatter, affecting the clustering power spectrum, although this phenomenon might not be statistically significant. Therefore, we should adopt appropriate approaches that can consistently interpret the LIM power spectrum from observations., Comment: 9 pages, 8 figures, comments are welcome, accepted in MNRAS

A wealth of cosmological and astrophysical information is expected from many ongoing and upcoming large-scale surveys. It is crucial to prepare for these surveys now and develop tools that can efficiently extract most information. We present HIFlow: a fast generative model of the neutral hydrogen (HI) maps that is conditioned only on cosmology ($\Omega_{m}$ and $\sigma_{8}$) and designed using a class of normalizing flow models, the Masked Autoregressive Flow (MAF). HIFlow is trained on the state-of-the-art simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project. HIFlow has the ability to generate realistic diverse maps without explicitly incorporating the expected 2D maps structure into the flow as an inductive bias. We find that HIFlow is able to reproduce the CAMELS average and standard deviation HI power spectrum (Pk) within a factor of $\lesssim$ 2, scoring a very high $R^{2} > 90\%$. By inverting the flow, HIFlow provides a tractable high-dimensional likelihood for efficient parameter inference. We show that the conditional HIFlow on cosmology is successfully able to marginalize over astrophysics at the field level, regardless of the stellar and AGN feedback strengths. This new tool represents a first step toward a more powerful parameter inference, maximizing the scientific return of future HI surveys, and opening a new avenue to minimize the loss of complex information due to data compression down to summary statistics., Comment: updated: 14 pages, 10 figures, a new section on inference has been added during revision. Accepted for publication in ApJ

We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) Multifield Dataset, CMD, a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from 2,000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span $\sim$100 million light years and have been generated from thousands of state-of-the-art hydrodynamic and gravity-only N-body simulations from the CAMELS project. Designed to train machine learning models, CMD is the largest dataset of its kind containing more than 70 Terabytes of data. In this paper we describe CMD in detail and outline a few of its applications. We focus our attention on one such task, parameter inference, formulating the problems we face as a challenge to the community. We release all data and provide further technical details at https://camels-multifield-dataset.readthedocs.io., Comment: 17 pages, 1 figure. Third paper of a series of four. Hundreds of thousands of labeled 2D maps and 3D grids from thousands of simulated universes publicly available at https://camels-multifield-dataset.readthedocs.io

Star formation in half of massive galaxies was quenched by the time the Universe was three billion years old. Very low amounts of molecular gas appear responsible for this, at least in some cases, though morphological gas stabilization, shock heating, or activity associated with accretion onto a central supermassive black hole is invoked in other cases. Recent studies of quenching by gas depletion have been based upon upper limits that are insufficiently sensitive to determine this robustly, or stacked emission with its problems of averaging. Here we report 1.3mm observations of dust emission from six strongly lensed galaxies where star formation has been quenched, with magnifications of up to a factor of 30. Four of the six galaxies are undetected in dust emission, with an estimated upper limit on the dust mass of 0.0001 times the stellar mass, and by proxy (assuming a Milky Way molecular gas-to-dust ratio) 0.01 times the stellar mass in molecular gas. This is two orders of magnitude less molecular gas per unit stellar mass than seen in star forming galaxies at similar redshifts. It remains difficult to extrapolate from these small samples, but these observations establish that gas depletion is responsible for a cessation of star formation in some fraction of high-redshift galaxies., Comment: 17 pages, 3 figures. Authors' version. Published online by Nature on September 22, 2021

We train neural networks to perform likelihood-free inference from $(25\,h^{-1}{\rm Mpc})^2$ 2D maps containing the total mass surface density from thousands of hydrodynamic simulations of the CAMELS project. We show that the networks can extract information beyond one-point functions and power spectra from all resolved scales ($\gtrsim 100\,h^{-1}{\rm kpc}$) while performing a robust marginalization over baryonic physics at the field level: the model can infer the value of $\Omega_{\rm m} (\pm 4\%)$ and $\sigma_8 (\pm 2.5\%)$ from simulations completely different to the ones used to train it., Comment: 7 pages, 4 figures. Second paper of a series of four. The 2D maps, codes, and network weights used in this paper are publicly available at https://camels-multifield-dataset.readthedocs.io

Astrophysical processes such as feedback from supernovae and active galactic nuclei modify the properties and spatial distribution of dark matter, gas, and galaxies in a poorly understood way. This uncertainty is one of the main theoretical obstacles to extract information from cosmological surveys. We use 2,000 state-of-the-art hydrodynamic simulations from the CAMELS project spanning a wide variety of cosmological and astrophysical models and generate hundreds of thousands of 2-dimensional maps for 13 different fields: from dark matter to gas and stellar properties. We use these maps to train convolutional neural networks to extract the maximum amount of cosmological information while marginalizing over astrophysical effects at the field level. Although our maps only cover a small area of $(25~h^{-1}{\rm Mpc})^2$, and the different fields are contaminated by astrophysical effects in very different ways, our networks can infer the values of $\Omega_{\rm m}$ and $\sigma_8$ with a few percent level precision for most of the fields. We find that the marginalization performed by the network retains a wealth of cosmological information compared to a model trained on maps from gravity-only N-body simulations that are not contaminated by astrophysical effects. Finally, we train our networks on multifields -- 2D maps that contain several fields as different colors or channels -- and find that not only they can infer the value of all parameters with higher accuracy than networks trained on individual fields, but they can constrain the value of $\Omega_{\rm m}$ with higher accuracy than the maps from the N-body simulations., Comment: 11 pages, 7 figures. First paper of a series of four. All 2D maps, codes, and networks weights publicly available at https://camels-multifield-dataset.readthedocs.io

We use a generic formalism designed to search for relations in high-dimensional spaces to determine if the total mass of a subhalo can be predicted from other internal properties such as velocity dispersion, radius, or star-formation rate. We train neural networks using data from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project and show that the model can predict the total mass of a subhalo with high accuracy: more than 99% of the subhalos have a predicted mass within 0.2 dex of their true value. The networks exhibit surprising extrapolation properties, being able to accurately predict the total mass of any type of subhalo containing any kind of galaxy at any redshift from simulations with different cosmologies, astrophysics models, subgrid physics, volumes, and resolutions, indicating that the network may have found a universal relation. We then use different methods to find equations that approximate the relation found by the networks and derive new analytic expressions that predict the total mass of a subhalo from its radius, velocity dispersion, and maximum circular velocity. We show that in some regimes, the analytic expressions are more accurate than the neural networks. We interpret the relation found by the neural network and approximated by the analytic equation as being connected to the virial theorem., Comment: 23 pages, 15 figures