Search

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

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

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., 17 pages, 3 figures. Authors' version. Published online by Nature on September 22, 2021

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., Published in ApJL, 8 pages, 5 figures

Matching the number counts of high-$z$ sub-millimetre-selected galaxies (SMGs) has been a long standing problem for galaxy formation models. In this paper, we use 3D dust radiative transfer to model the sub-mm emission from galaxies in the SIMBA cosmological hydrodynamic simulations, and compare predictions to the latest single-dish observational constraints on the abundance of 850$\mathrm{\mu m}$-selected sources. We find good agreement with the shape of the integrated 850$\mathrm{\mu m}$ luminosity function, and the normalisation is within 0.25 dex at $> 3 \; \mathrm{mJy}$, unprecedented for a fully cosmological hydrodynamic simulation, along with good agreement in the redshift distribution of bright SMGs. The agreement is driven primarily by SIMBA's good match to infrared measures of the star formation rate (SFR) function between $z = 2-4$ at high SFRs. Also important is the self-consistent on-the-fly dust model in SIMBA, which predicts, on average, higher dust masses (by up to a factor of 2.5) compared to using a fixed dust-to-metals ratio of 0.3. We construct a lightcone to investigate the effect of far-field blending, and find that 52% of sources are blends of multiple components, which makes a small contribution to the normalisation of the bright-end of the number counts. We provide new fits to the 850$\mathrm{\mu m}$ luminosity as a function of SFR and dust mass. Our results demonstrate that exotic solutions to the discrepancy between sub-mm counts in simulations and observations, such as a top-heavy IMF, are unnecessary, and that sub-millimetre-bright phases are a natural consequence of massive galaxy evolution., Comment: 23 pages, 21 figures, Accepted to MNRAS

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 nonuniform screen model parameterized by the fraction of unobscured stellar light. This nonuniform screen model allows for a nonzero 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 nonuniform 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 nonuniform screen model. With this increase in dust attenuation modeling accuracy, we also improve the star formation rates (SFRs) inferred with the nonuniform 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.

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., 33 pages, 24 figures, 4 tables, ApJ in press

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., 25 pages, 15 figures. Accepted for publication in ApJ

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 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 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., 20 pages, 13 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, Comment: Accepted for publication in ApJ Letters, 9 pages, 5 figures. Interactive versions of Fig. 2 and data products are available at https://github.com/spt-smg/publicdata

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

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$σ$ 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., Submitted to AAS Journals on December 15, 2021. Comments welcome

One of the most fundamental baryonic matter components of galaxies is the neutral atomic hydrogen (H i). At low redshifts, this component can be traced directly through the 21 cm transition, but to infer the H 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 (2 P 3/2 − 2 P 1/3) transition of singly ionized carbon C ii at 158 μm as a proxy for H i in a set simulated galaxies at z ≈ 6, following the work by Heintz et al. We select 11,125 star-forming galaxies from the simba simulations, with far-infrared line emissions postprocessed and modeled within the Sigame framework. We find a strong connection between C ii and H i, with the relation between this C ii-to-H i relation (β [C II]) being anticorrelated 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 ≈ 6, and specifically the H i gas mass fraction. We find mean values of M H I/M ⋆ = 1.4 and M H I/M bar,tot = 0.45. These results provide strong evidence for H i being the dominant baryonic matter component by mass in galaxies at z ≈ 6.

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

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, Comment: 23 pages, 16 figures, accepted to ApJ

Observing the interstellar medium (ISM) in $z \gtrsim 6$ quasars host galaxies is essential for understanding the co-evolution between the supermassive black holes and their hosts. To probe the gas physical conditions and search for imprints of Active Galactic Nuclei (AGN) on the ISM, we report ALMA observations of the $\rm [N\ II]_{122 \mu m}$ and $\rm [O\ I]_{146 \mu m}$ lines and the underlying continuum from the $z=6.003$ quasar SDSS J231038.88+185519.7. Together with previous $\rm [C\ II]_{158 \mu m}$ and $\rm [O\ III]_{88 \mu m}$ observations, we use the ratios of these fine-structure lines to probe the ISM properties. Similar to other high-$z$ systems, this object exhibits a $\rm [C\ II]_{158 \mu m}$/$\rm [O\ I]_{146 \mu m}$ ratio comparable to the lowest values found in local (Ultra) luminous infrared galaxies, suggesting a "warmer" and "denser" gas component compared to typical local systems. The $\rm [O\ III]_{88 \mu m}$/$\rm [O\ I]_{146 \mu m}$ ratio is lower than that of other local and high-$z$ systems, indicating a smaller ionized gas fraction in this quasar. The $\rm [O\ III]_{88 \mu m}$/$\rm [N\ II]_{122 \mu m}$ ratio is comparable to that of local systems, and suggests a metallicity of $Z/Z_{\odot}$=1.5$-$2.1. Based on the $\rm [N\ II]_{122 \mu m}$ detection, we estimate that $17\%$ of the $\rm [C\ II]_{158 \mu m}$ emission is associated with ionized gas. The $\rm [N\ II]_{122 \mu m}$ line shows a "flux deficit" comparable to local systems. The $\rm [O\ I]_{146 \mu m}$ line, with a $\rm [O\ I]_{146 \mu m}$/FIR ratio $\ge 2\times$ than expected from the local relation, indicates no $\rm [O\ I]_{\rm 146 \mu m}$ deficit. The low $\rm [C\ II]_{158 \mu m}$/$\rm [O\ I]_{146 \mu m}$ ratio, together with the high $\rm [O\ I]_{146 \mu m}$/FIR ratio in J2310+1855, reveals that the warm and dense gas is likely a result of AGN heating to the ISM., Comment: Accepted for publication in Apj. 15 pages, 4 figures

We develop a cosmological model for the evolution of dust grains in galaxies with a distribution of sizes in order to understand the origin of the Milky Way dust extinction curve. Our model considers the formation of active dust in evolved stars, growth by accretion and coagulation, and destruction processes via shattering, sputtering, and astration in the ISM of galaxies over cosmic time. Our main results follow. Galaxies in our cosmological model with masses comparable to the Milky Way's at z~0 exhibit a diverse range of extinction laws, though with slopes and bump strengths comparable to the range observed in the Galaxy. The progenitors of the Milky Way have steeper slopes, and only flatten to slopes comparable to the Galaxy at $z \approx 1$. This owes to increased grain growth rates at late times/in high-metallicity environments driving up the ratio of large to small grains, with a secondary dependence on the graphite to silicate ratio evolution. The UV bump strengths depend primarily on the graphite to silicate ratio, and remain broadly constant in MW-like galaxies between z=3 and z=0, though show slight variability. Our models span comparable regions of bump-slope space as sightlines in the Galaxy do, though there is a lack of clear relationship between the model slopes and bump strengths owing to small scale fluctuations in the bump strength. Our models naturally produce slopes for some non-Milky Way analogs as steep as those of the LMC and SMC in metal poor galaxies, though notably the bump strengths are, on average, too large when comparing to the Magellanic clouds. This owes to the fact that we evolve the grain size distributions of graphites and silicates simultaneously, which is an oversimplification. Our model provides a novel framework to study the origins and variations of dust extinction curves in galaxies over cosmic time., Comment: 12 pages, 6 figures, MNRAS submitted

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., 23 pages, 15 figures

In this invited talk, we discuss the physics of the lifecycle of dust in the context of galaxy formation simulations. After outlining the basic physical processes, we apply algorithms for the formation, growth, and destruction of dust in the ISM to a state-of-the-art cosmological simulation to develop a model for the evolution of the dust to gas and dust to metals ratios in galaxies. We show that while modern simulations are able to match the observed dust mass function at redshift z = 0, most models underpredict the observed mass function at high-redshift (z = 2). We then show the power of these techniques by expanding our model to include a spectrum of dust sizes, and make initial predictions for extinction laws in local galaxies.

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., 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 Origins Space Telescope (Origins) concept is designed to investigate the creation and dispersal of elements essential to life, the formation of planetary systems, and the transport of water to habitable worlds and the atmospheres of exoplanets around nearby K- and M-dwarfs to identify potentially habitable—and even inhabited—worlds. These science priorities are aligned with NASA’s three major astrophysics science goals: How does the Universe work? How did we get here? and Are we alone? We briefly describe the science case that arose from the astronomical community and the science traceability matrix for Origins. The science traceability matrix prescribes the design of Origins and demonstrates that it will address the key science questions motivated by the science case.

We present the first detailed dissection of the circumgalactic medium (CGM) of massive starburst galaxies at z > 2. Our target is a submillimeter galaxy (SMG) at z = 2.674 that has a star formation rate of 1200 $M_\odot$/yr and a molecular gas reservoir of $1.3\times10^{11} M_\odot$. We characterize its CGM with two background QSOs at impact parameters of 93 kpc and 176 kpc. We detect strong HI and metal-line absorption near the redshift of the SMG towards both QSOs, each consisting of three main subsystems spanning over 1500 km/s. The absorbers show remarkable kinematic and metallicity coherence across a separation of 86 kpc. In particular, the cool gas in the CGM of the SMG exhibits high HI column densities ($\log N_{\rm HI}/{\rm cm}^{-2} = 20.2, 18.6$), low metallicities ([M/H] $\approx$ -2.0), and similar radial velocities ($\approx$ -300 km/s). While the HI column densities match previous results on the CGM around QSOs at z > 2, the metallicities are lower by more than an order of magnitude, making it an outlier in the line width$-$metallicity relation of damped Ly$��$ absorbers. The large physical extent, the velocity coherence, the high surface density, and the low metallicity are all consistent with the cool, inflowing, and near-pristine gas streams predicted to penetrate hot massive halos at z > 1.5. We estimate a total gas accretion rate of ~100 $M_\odot$/yr from three such streams, which falls short of the star formation rate but is consistent with simulations. At this rate, it takes about a gigayear to acquire the molecular gas reservoir of the central starburst., ApJ in press. Data files available at https://github.com/fuhaiastro/data

The Origins Space Telescope will trace the history of our origins from the time dust and heavy elements permanently altered the cosmic landscape to present-day life. How did galaxies evolve from the earliest galactic systems to those found in the Universe today? How do habitable planets form? How common are life-bearing worlds? To answer these alluring questions, Origins will operate at mid- and far-infrared (IR) wavelengths and offer powerful spectroscopic instruments and sensitivity three orders of magnitude better than that of the Herschel Space Observatory, the largest telescope flown in space to date. We describe the baseline concept for Origins recommended to the 2020 US Decadal Survey in Astronomy and Astrophysics. The baseline design includes a 5.9-m diameter telescope cryocooled to 4.5 K and equipped with three scientific instruments. A mid-infrared instrument (Mid-Infrared Spectrometer and Camera Transit spectrometer) will measure the spectra of transiting exoplanets in the 2.8 to 20 μm wavelength range and offer unprecedented spectrophotometric precision, enabling definitive exoplanet biosignature detections. The far-IR imager polarimeter will be able to survey thousands of square degrees with broadband imaging at 50 and 250 μm. The Origins Survey Spectrometer will cover wavelengths from 25 to 588 μm, making wide-area and deep spectroscopic surveys with spectral resolving power R ∼ 300, and pointed observations at R ∼ 40,000 and 300,000 with selectable instrument modes. Origins was designed to minimize complexity. The architecture is similar to that of the Spitzer Space Telescope and requires very few deployments after launch, while the cryothermal system design leverages James Webb Space Telescope technology and experience. A combination of current-state-of-the-art cryocoolers and next-generation detector technology will enable Origins’ natural background-limited sensitivity.