Galaxy Stellar Mass Research Articles

Dependence of peculiar velocity on the host properties of the gravitational wave sources and its impact on the measurement of Hubble constant

ABSTRACT Accurate measurement of the Hubble constant from standard sirens such as the gravitational wave (GW) sources with electromagnetic counterparts relies on the robust peculiar velocity correction of the redshift of the host galaxy. We show in this work that the peculiar velocity of the host galaxies exhibits a correlation with the properties of the host galaxy primarily such as its stellar mass and this correlation also evolves with redshift. As the galaxies of higher stellar mass tend to form in galaxies with higher halo masses which are located in spatial regions having a non-linear fluctuation in the density field of the matter distribution, the root mean square peculiar velocity of more massive galaxies is higher. As a result, depending on the formation channel of the binary compact objects, the peculiar velocity contamination to the galaxies will be different. The variation in the peculiar velocity of the host galaxies can lead to a significant variation in the estimation of the Hubble constant inferred using sources such as binary neutron stars. For the network of GW detectors such as LIGO-Virgo-KAGRA (LVK), LVK+LIGO-India, and Cosmic Explorer + Einstein Telescope, the variation in the precision of Hubble constant inferred from 10 bright siren events can vary from $\sim 5.4 - 6~{{\ \rm per \, cent}}$, $\sim 4.5 - 5.3~{{\ \rm per \, cent}}$, and $\sim 1.1 - 2.7~{{\ \rm per \, cent}}$, respectively. The impact of such a correlation between peculiar velocity and stellar mass on the inference of the Hubble constant is not only limited to GW sources but also applicable to type-Ia supernovae.

How many stars form in galaxy mergers?

ABSTRACT We forward model the difference in stellar age between post-coalescence mergers and a control sample with the same stellar mass, environmental density, and redshift. In particular, we use a pure sample of 445 post-coalescence mergers from the recent visually confirmed post-coalescence merger sample identified by Bickley et al. and find that post-coalescence mergers are on average younger than control galaxies for 10 < log (M⋆/M⊙) < 11. The difference in age from matched controls is up to 1.5 Gyr, highest for lower stellar mass galaxies. We forward model this difference using parametric star formation histories, accounting for the pre-coalescence inspiral phase of enhanced star formation using close pair data, and a final additive burst of star formation at coalescence. We find a best-fitting stellar mass burst fraction of fburst = ΔM⋆/M⋆, merger = 0.18 ± 0.02 for 10 < log (M⋆/M⊙) < 11 galaxies, with no evidence of a trend in stellar mass. The modelled burst fraction is robust to choice of parametric star formation history, as well as differences in burst duration. The result appears consistent with some prior observationally derived values, but is significantly higher than that found in hydrodynamical simulations. Using published luminous infrared galaxy star formation rates, we find a burst duration increasing with stellar mass, from 120 to 250 Myr. A comparison to published cold gas measurements indicates there is enough molecular gas available in very close pairs to fuel the burst. Additionally, given our stellar mass burst estimate, the predicted cold gas fraction remaining after the burst is consistent with observed post-coalescence mergers.

Transonic galactic wind model including stellar feedbacks and application to outflows in high/low-z galaxies

Abstract Galactic winds play a crucial role in the ejection of the interstellar medium (ISM) into intergalactic space. This study presents a model that classifies possible transonic solutions of galactic winds in the gravitational potential of the dark matter halo and stellar component under spherically symmetric and steady assumptions. Our model includes injections of mass and energy resulting from supernova feedback along a flow line. The mass flux in galactic winds is a critical factor in determining the acceleration process of the flow and revealing the impact of galactic winds on galaxy evolution. We apply the transonic galactic wind model to the observed outflow velocities of star-forming galaxies to estimate the mass flux. Dividing the mass flux by the star formation rate (SFR) yields the mass loading rate (and mass loading factor), which indicates the entrainment effect of the ISM by the hot gas flow. Our results demonstrate that the mass loading rate is inversely correlated with galaxy mass and SFR. In less massive galaxies (stellar mass ∼107−8 M⊙), the mass loading rate exceeds unity, indicating effective ejection of the ISM into intergalactic space. However, in massive galaxies (stellar mass ∼1010−11 M⊙), the mass loading rate falls below unity, meaning that the mass flux cannot exceed the injected mass by supernovae, thus resulting in ineffective ejection of the ISM.

FLAMINGO: calibrating large cosmological hydrodynamical simulations with machine learning.

To fully take advantage of the data provided by large-scale structure surveys, we need to quantify the potential impact of baryonic effects, such as feedback from active galactic nuclei (AGN) and star formation, on cosmological observables. In simulations, feedback processes originate on scales that remain unresolved. Therefore, they need to be sourced via subgrid models that contain free parameters. We use machine learning to calibrate the AGN and stellar feedback models for the FLAMINGO (Fullhydro Large-scale structure simulations with All-sky Mapping for the Interpretation of Next Generation Observations) cosmological hydrodynamical simulations. Using Gaussian process emulators trained on Latin hypercubes of 32 smaller volume simulations, we model how the galaxy stellar mass function (SMF) and cluster gas fractions change as a function of the subgrid parameters. The emulators are then fit to observational data, allowing for the inclusion of potential observational biases. We apply our method to the three different FLAMINGO resolutions, spanning a factor of 64 in particle mass, recovering the observed relations within the respective resolved mass ranges. We also use the emulators, which link changes in subgrid parameters to changes in observables, to find models that skirt or exceed the observationally allowed range for cluster gas fractions and the SMF. Our method enables us to define model variations in terms of the data that they are calibrated to rather than the values of specific subgrid parameters. This approach is useful, because subgrid parameters are typically not directly linked to particular observables, and predictions for a specific observable are influenced by multiple subgrid parameters.

Binaries drive high Type Ia supernova rates in dwarf galaxies

ABSTRACT The scaling of the specific Type Ia supernova (SN Ia) rate with host galaxy stellar mass $\dot{\text{N}}_\text{Ia} / \text{M}_\star \sim \text{M}_\star ^{-0.3}$ as measured in ASAS-SN and DES strongly suggests that the number of SNe Ia produced by a stellar population depends inversely on its metallicity. We estimate the strength of the required metallicity dependence by combining the average star formation histories (SFHs) of galaxies as a function of their stellar mass with the mass–metallicity relation (MZR) for galaxies and common parametrizations for the SN Ia delay-time distribution. The differences in SFHs can account for only ∼30 per cent of the increase in the specific SN Ia rate between stellar masses of M⋆ = 1010 and 107.2 M⊙. We find that an additional metallicity dependence of approximately ∼Z−0.5 is required to explain the observed scaling. This scaling matches the metallicity dependence of the close binary fraction observed in APOGEE, suggesting that the enhanced SN Ia rate in low-mass galaxies can be explained by a combination of their more extended SFHs and a higher binary fraction due to their lower metallicities. Due to the shape of the MZR, only galaxies below M⋆ ≈ 3 × 109 M⊙ are significantly affected by the metallicity-dependent SN Ia rates. The $\dot{\text{N}}_\text{Ia} / \text{M}_\star \sim \text{M}_\star ^{-0.3}$ scaling becomes shallower with increasing redshift, dropping by factor of ∼2 at 107.2 M⊙ between z = 0 and 1 with our ∼Z−0.5 scaling. With metallicity-independent rates, this decrease is a factor of ∼3. We discuss the implications of metallicity-dependent SN Ia rates for one-zone models of galactic chemical evolution.

IllustrisTNG in the HSC-SSP: image data release and the major role of mini mergers as drivers of asymmetry and star formation

ABSTRACT At fixed galaxy stellar mass, there is a clear observational connection between structural asymmetry and offset from the star-forming main sequence, ΔSFMS. Herein, we use the TNG50 simulation to investigate the relative roles of major mergers (stellar mass ratios μ ≥ 0.25), minor (0.1 ≤ μ < 0.25), and mini mergers (0.01 ≤ μ < 0.1) in driving this connection amongst star-forming galaxies (SFGs). We use dust radiative transfer post-processing with SKIRT to make a large, public collection of synthetic Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) images of simulated IllustrisTNG (TNG) galaxies over 0.1 ≤ z ≤ 0.7 with log (M⋆/M⊙) ≥ 9 (∼750 k images). Using their instantaneous star formation rates (SFRs), known merger histories/forecasts, and HSC-SSP asymmetries, we show (1) that TNG50 SFGs qualitatively reproduce the observed trend between ΔSFMS and asymmetry and (2) a strikingly similar trend emerges between ΔSFMS and the time-to-coalescence for mini mergers. Controlling for redshift, stellar mass, environment, and gas fraction, we show that individual mini merger events yield small enhancements in SFRs and asymmetries that are sustained on long time-scales (at least ∼3 Gyr after coalescence, on average) – in contrast to major/minor merger remnants which peak at much greater amplitudes but are consistent with controls only ∼1 Gyr after coalescence. Integrating the boosts in SFRs and asymmetries driven by μ ≥ 0.01 mergers since z = 0.7 in TNG50 SFGs, we show that mini mergers are responsible for (i) 55 per cent of all merger-driven star formation and (ii) 70 per cent of merger-driven asymmetric structure. Due to their relative frequency and prolonged boost time-scales, mini mergers dominate over their minor and major counterparts in driving star formation and asymmetry in SFGs.

The impact of stochastic modelling on the predictive power of galaxy formation simulations

ABSTRACT All modern galaxy formation models employ stochastic elements in their sub-grid prescriptions to discretize continuous equations across the time domain. In this paper, we investigate how the stochastic nature of these models, notably star formation, black hole accretion, and their associated feedback, that act on small (< kpc) scales, can back-react on macroscopic galaxy properties (e.g. stellar mass and size) across long (> Gyr) time-scales. We find that the scatter in scaling relations predicted by the EAGLE model implemented in the swift code can be significantly impacted by random variability between re-simulations of the same object, even when galaxies are resolved by tens of thousands of particles. We then illustrate how re-simulations of the same object can be used to better understand the underlying model, by showing how correlations between galaxy stellar mass and black hole mass disappear at the highest black hole masses (MBH > 108 M⊙), indicating that the feedback cycle may be interrupted by external processes. We find that although properties that are collected cumulatively over many objects are relatively robust against random variability (e.g. the median of a scaling relation), the properties of individual galaxies (such as galaxy stellar mass) can vary by up to 25 per cent, even far into the well-resolved regime, driven by bursty physics (black hole feedback), and mergers between galaxies. We suggest that studies of individual objects within cosmological simulations be treated with caution, and that any studies aiming to closely investigate such objects must account for random variability within their results.

Supermassive black holes in a mass-limited galaxy sample

ABSTRACT The observed scaling relations between supermassive black hole masses and their host galaxy properties indicate that supermassive black holes influence the evolution of galaxies. However, the scaling relations may be affected by selection biases. We propose to measure black hole masses in a mass-limited galaxy sample including all non-detections to inprove constraints on galaxy mass – black hole mass scaling relations and test for selection bias. We use high-spatial resolution spectroscopy from the Keck and Gemini telescopes, and the Jeans Anisotropic Modelling method to measure black hole masses in early-type galaxies from the Virgo Cluster. We present four new black hole masses and one upper limit in our mass-selected sample of galaxies of galaxy mass (1.0–3.2) $\times 10^{10} \, \mathrm{M}_\odot$. This brings the total measured to 11 galaxies out of a full sample of 18 galaxies, allowing us to constrain scaling relations. We calculate a lower limit for the average black hole mass in our sample of $3.7 \times 10^{7} \, \mathrm{M}_\odot$. This is at an average galaxy stellar mass of $(1.81 \pm 0.14)\times 10^{10} \, \mathrm{M}_\odot$ and an average bulge mass of $(1.31 \pm 0.15) \times 10^{10} \, \mathrm{M}_\odot$. This lower limit shows that black hole masses in early-type galaxies are not strongly affected by selection biases.

Simulating the LOcal Web (SLOW)

Context. Several observations of the Local Universe point toward the existence of very prominent structures: massive galaxy clusters and local superclusters on the one hand, but also large local voids and underdensities on the other. However, it is highly nontrivial to connect such different observational selected tracers to the underlying dark matter (DM) distribution. Aims. Therefore, constructing mock catalogs of such observable tracers using cosmological hydrodynamics simulations is needed. These simulations have to follow galaxy formation physics and also have to be constrained to reproduce the Local Universe. Such constraints should be based on observables that directly probe the full underlying gravitational field, such as the observed peculiar velocity field, to provide an independent test on the robustness of these distinctive structures. Methods. We used a 500 h−1 Mpc constrained simulation of the Local Universe to investigate the anomalies in the local density field, as found in observations. Constructing the initial conditions based on peculiar velocities derived from the CosmicFlows-2 catalog makes the predictions of the simulations completely independent from the distribution of the observed tracer population, and following galaxy formation physics directly in the hydrodynamics simulations also allows the comparison to be based directly on the stellar masses of galaxies or X-ray luminosity of clusters. We also used the 2668 h−1 Mpc large cosmological box from the Magneticum simulations to evaluate the frequency of finding such anomalies in random patches within simulations. Results. We demonstrate that halos and galaxies in our constrained simulation trace the local dark matter density field very differently. Thus, this simulation reproduces the observed 50% underdensity of galaxy clusters and groups within the sphere of ≈100 Mpc when applying the same mass or X-ray luminosity limit used in the observed cluster sample (CLASSIX), which is consistent with a ≈1.5σ feature. At the same time, the simulation reproduces the observed overdensity of massive galaxy clusters within the same sphere, which on its own also corresponds to a ≈1.5σ feature. Interestingly, we find that only 44 out of 15 635 random realizations (i.e., 0.28%) match both anomalies, thus making the Local Universe a ≈3σ environment. We finally compared a mock galaxy catalog with the observed distribution of galaxies in the Local Universe, finding a match to the observed factor of 2 overdensity at ∼16 Mpc as well as the observed 15% underdensity at ∼40 Mpc. Conclusions. Constrained simulations of the Local Universe which reproduce the main features of the local density field open a new window for local field cosmology, where the imprint of the specific density field and the impact on the bias through the observational specific tracers can be investigated in detail.

The HERBAL Model: A Hierarchical Errors-in-variables Bayesian Lognormal Hurdle Model for Galactic Globular Cluster Populations

Galaxy stellar mass is known to be monotonically related to the size of the galaxy’s globular cluster (GC) population for Milky Way sized and larger galaxies. However, the relation becomes ambiguous for dwarf galaxies, where there is some evidence for a downturn in GC population size at low galaxy masses. Smaller dwarfs are increasingly likely to have no GCs, and these zeros cannot be easily incorporated into linear models. We introduce the Hierarchical Errors-in-variables ERrors-in-variables BAyesian Lognormal hurdle (HERBAL) model to represent the relationship between dwarf galaxies and their GC populations, and apply it to the sample of Local Group galaxies, where the luminosity range coverage is maximal. This bimodal model accurately represents the two populations of dwarf galaxies: those that have GCs and those that do not. Our model thoroughly accounts for all uncertainties, including measurement uncertainty, uncertainty in luminosity to stellar mass conversions, and intrinsic scatter. The hierarchical nature of our Bayesian model also allows us to estimate galaxy masses and individual mass-to-light ratios from luminosity data within the model. We find that 50% of galaxies are expected to host GC populations at a stellar mass of , and that the expected mass of GC populations remains linear down to the smallest galaxies. Our hierarchical model recovers an accurate estimate of the Milky Way stellar mass. Under our assumed error model, we find a nonzero intrinsic scatter of (95% credible interval) that should be accounted for in future models.

Properties and merger signatures of galaxies hosting LISA coalescing massive black hole binaries

The gravitational wave (GW) antenna LISA will detect the signal from coalescing massive black hole binaries (MBHBs) of 104 − 107 M⊙, providing clues as to their formation and growth throughout cosmic history. Some of these events will be localised with a precision of several to less than a deg2, enabling the possible identification of their host galaxy. This work explores the properties of the host galaxies of LISA MBHBs below z ≲ 3. We generate a simulated lightcone using the semi-analytical model L-Galaxies applied to the merger trees of the high-resolution N-body cosmological simulation Millennium-II. The model shows that LISA MBHBs are expected to be found in optically dim (r > 20), star-forming (sSFR > 10−10 yr−1), gas-rich (fgas > 0.6), and disc-dominated (B/T < 0.7) low-mass galaxies of stellar masses 108 − 109 M⊙. However, these properties are indistinguishable from those of galaxies harbouring single massive black holes of comparable mass, making it difficult to select LISA hosts among the whole population of low-mass galaxies. Motivated by this, we explore the possibility of using merger signatures to select LISA hosts. We find that 40%−80% of the galaxies housing LISA MBHBs display merger features related to the interaction that brought the secondary MBH to the galaxy. Despite this, around 60% of dwarf galaxies placed in the surroundings of the LISA hosts will show these kinds of features as well, challenging the unequivocal detection of LISA hosts through the search for merger signatures. Consequently, the detection of an electromagnetic transient associated with the MBHB merger will be vital in order to pinpoint the star-forming dwarf galaxy where these binary systems evolve and coalesce.

COSMOS2020: The galaxy stellar mass function

Context. How galaxies form, assemble, and cease their star formation is a central question within the modern landscape of galaxy evolution studies. These processes are indelibly imprinted on the galaxy stellar mass function (SMF), and its measurement and understanding is key to uncovering a unified theory of galaxy evolution. Aims. We present constraints on the shape and evolution of the galaxy SMF, the quiescent galaxy fraction, and the cosmic stellar mass density across 90% of the history of the Universe from z = 7.5 → 0.2 as a means to study the physical processes that underpin galaxy evolution. Methods. The COSMOS survey is an ideal laboratory for studying representative galaxy samples. Now equipped with deeper and more homogeneous near-infrared coverage exploited by the COSMOS2020 catalog, we leverage the large 1.27 deg2 effective area to improve sample statistics and understand spatial variations (cosmic variance) – particularly for rare, massive galaxies – and push to higher redshifts with greater confidence and mass completeness than previous studies. We divide the total stellar mass function into star-forming and quiescent subsamples through NUVrJ color-color selection. The measurements are then fit with single- and double-component Schechter functions to infer the intrinsic galaxy stellar mass function, the evolution of its key parameters, and the cosmic stellar mass density out to z = 7.5. Finally, we compare our measurements to predictions from state-of-the-art cosmological simulations and theoretical dark matter halo mass functions. Results. We find a smooth, monotonic evolution in the galaxy stellar mass function since z = 7.5, in general agreement with previous studies. The number density of star-forming systems have undergone remarkably consistent growth spanning four decades in stellar mass from z = 7.5 → 2 whereupon high-mass systems become predominantly quiescent (“downsizing”). Meanwhile, the assembly and growth of low-mass quiescent systems only occurred recently, and rapidly. An excess of massive systems at z ≈ 2.5 − 5.5 with strikingly red colors, with some being newly identified, increase the observed number densities to the point where the SMF cannot be reconciled with a Schechter function. Conclusions. Systematics including cosmic variance and/or active galactic nuclei contamination are unlikely to fully explain this excess, and so we speculate that they may be dust-obscured populations similar to those found in far infrared surveys. Furthermore, we find a sustained agreement from z ≈ 3 − 6 between the stellar and dark matter halo mass functions for the most massive systems, suggesting that star formation in massive halos may be more efficient at early times.

Are we surprised to find SMBHs with JWST at z ≥ 9?

ABSTRACT JWST is unveiling for the first time accreting black holes (BHs) with masses of 106 − 107 M⊙ at z > 4, with the most distant residing in GNz11 at z = 10.6. Are we really surprised to find them in the nuclei of z ≃ 5 − 11 galaxies? Here we predict the properties of 4 < z < 11 BHs and their host galaxies considering an Eddington-limited (EL) and a super-Eddington (SE) BH accretion scenario, using the Cosmic Archaeology Tool (CAT) semi-analytical model. We calculate the transmitted spectral energy distribution of CAT synthetic candidates, representative of the BH/galaxy properties of GNz11. We also examine the possibility that the z = 8.7 galaxy CEERS-1019 could host an active BH. We find that the luminosity of high-z JWST detected BHs are better reproduced by the SE model, where BHs descend from efficiently growing light and heavy seeds. Conversely, the host galaxy stellar masses are better matched in the EL model, in which all the systems detectable with JWST surveys JADES and CEERS descend from heavy BH seeds. We support the interpretation that the central point source of GNz11 could be powered by a SE (λEdd ≃ 2 − 3) accreting BH with mass 1.5 × 106 M⊙, while the emission from CEERS-1019 is dominated by the host galaxy; if it harbours an active BH, we find it to have a mass of MBH ≃ 107 M⊙, and to be accreting at sub-Eddington rates (λEdd ≃ 0.5).

Different higher order kinematics between star-forming and quiescent galaxies based on the SAMI, MAGPI, and LEGA-C surveys

ABSTRACT We present the first statistical study of spatially integrated non-Gaussian stellar kinematics spanning 7 Gyr in cosmic time. We use deep, rest-frame optical spectroscopy of massive galaxies (stellar mass $M_\star \gt 10^{10.5} \, \mathrm{M_\odot }$) at redshifts z = 0.05, 0.3, and 0.8 from the SAMI, MAGPI, and LEGA-C surveys, to measure the excess kurtosis h4 of the stellar velocity distribution, the latter parametrized as a Gauss–Hermite series. We find that at all redshifts where we have large enough samples, h4 anticorrelates with the ratio between rotation and dispersion, highlighting the physical connection between these two kinematic observables. In addition, and independently from the anticorrelation with rotation-to-dispersion ratio, we also find a correlation between h4 and M⋆, potentially connected to the assembly history of galaxies. In contrast, after controlling for mass, we find no evidence of independent correlation between h4 and aperture velocity dispersion or galaxy size. These results hold for both star-forming and quiescent galaxies. For quiescent galaxies, h4 also correlates with projected shape, even after controlling for the rotation-to-dispersion ratio. At any given redshift, star-forming galaxies have lower h4 compared to quiescent galaxies, highlighting the link between kinematic structure and star-forming activity.

Mass Dependence of Galaxy–Halo Alignment in LOWZ and CMASS

We measure the galaxy-ellipticity (GI) correlations for the Sloan Digital Sky Survey Data Release 12 LOWZ and CMASS samples with the shape measurements from the DESI Legacy Imaging Surveys. We model the GI correlations in an N-body simulation with our recent accurate stellar–halo mass relation from the Photometric object Around Cosmic webs (PAC) method. The large data set and our accurate modeling turns out an accurate measurement of the alignment angle between central galaxies and their host halos. We find that the alignment of central elliptical galaxies with their host halos increases monotonically with galaxy stellar mass or host halo mass, which can be well described by a power law for the massive galaxies. We also find that central elliptical galaxies are more aligned with their host halos in LOWZ than in CMASS, which might indicate an evolution of galaxy–halo alignment, though future studies are needed to verify this is not induced by the sample selections. In contrast, central disk galaxies are aligned with their host halos about 10 times more weakly in the GI correlation. These results have important implications for intrinsic alignment (IA) correction in weak lensing studies, IA cosmology, and theory of massive galaxy formation.

The FLAMINGO project: cosmological hydrodynamical simulations for large-scale structure and galaxy cluster surveys

ABSTRACT We introduce the Virgo Consortium’s FLAMINGO suite of hydrodynamical simulations for cosmology and galaxy cluster physics. To ensure the simulations are sufficiently realistic for studies of large-scale structure, the subgrid prescriptions for stellar and AGN feedback are calibrated to the observed low-redshift galaxy stellar mass function and cluster gas fractions. The calibration is performed using machine learning, separately for each of FLAMINGO’s three resolutions. This approach enables specification of the model by the observables to which they are calibrated. The calibration accounts for a number of potential observational biases and for random errors in the observed stellar masses. The two most demanding simulations have box sizes of 1.0 and 2.8 Gpc on a side and baryonic particle masses of 1 × 108 and $1\times 10^9\, \text{M}_\odot$, respectively. For the latter resolution, the suite includes 12 model variations in a 1 Gpc box. There are 8 variations at fixed cosmology, including shifts in the stellar mass function and/or the cluster gas fractions to which we calibrate, and two alternative implementations of AGN feedback (thermal or jets). The remaining 4 variations use the unmodified calibration data but different cosmologies, including different neutrino masses. The 2.8 Gpc simulation follows 3 × 1011 particles, making it the largest ever hydrodynamical simulation run to z = 0. Light-cone output is produced on-the-fly for up to 8 different observers. We investigate numerical convergence, show that the simulations reproduce the calibration data, and compare with a number of galaxy, cluster, and large-scale structure observations, finding very good agreement with the data for converged predictions. Finally, by comparing hydrodynamical and ‘dark-matter-only’ simulations, we confirm that baryonic effects can suppress the halo mass function and the matter power spectrum by up to ≈20 per cent.

The JWST FRESCO survey: legacy NIRCam/grism spectroscopy and imaging in the two GOODS fields

ABSTRACT We present the JWST cycle 1 53.8 h medium program FRESCO, short for ‘First Reionization Epoch Spectroscopically Complete Observations’. FRESCO covers 62 arcmin2 in each of the two GOODS/CANDELS fields for a total area of 124 arcmin2 exploiting JWST’s powerful new grism spectroscopic capabilities at near-infrared wavelengths. By obtaining ∼2 h deep NIRCam/grism observations with the F444W filter, FRESCO yields unprecedented spectra at R ∼ 1600 covering 3.8–5.0 µm for most galaxies in the NIRCam field of view. This setup enables emission line measurements over most of cosmic history, from strong PAH lines at z ∼ 0.2–0.5, to Pa α and Pa β at z ∼ 1–3, He i and [S iii] at z ∼ 2.5–4.5, H α and [N ii] at z ∼ 5–6.5, up to [O iii] and H β for z ∼ 7–9 galaxies. FRESCO’s grism observations provide total line fluxes for accurately estimating galaxy stellar masses and calibrating slit-loss corrections of NIRSpec/MSA spectra in the same field. Additionally, FRESCO results in a mosaic of F182M, F210M, and F444W imaging in the same fields to a depth of ∼28.2 mag (5σ in 0${_{.}^{\prime\prime}}$32 diameter apertures). Here, we describe the overall survey design and the key science goals that can be addressed with FRESCO. We also highlight several, early science results, including: spectroscopic redshifts of Lyman break galaxies that were identified almost 20 yr ago, the discovery of broad-line active galactic nuclei at z > 4, and resolved Pa α maps of galaxies at z ∼ 1.4. These results demonstrate the enormous power for serendipitous discovery of NIRCam/grism observations.

SIMBA-C: an updated chemical enrichment model for galactic chemical evolution in the SIMBA simulation

ABSTRACT We introduce a new chemical enrichment and stellar feedback model into GIZMO, using the SIMBA sub-grid models as a base. Based on the state-of-the-art chemical evolution model of Kobayashi et al., SIMBA-C tracks 34 elements from H→Ge and removes SIMBA’s instantaneous recycling approximation. Furthermore, we make some minor improvements to SIMBA’s base feedback models. SIMBA-C provides significant improvements on key diagnostics such as the knee of the z = 0 galaxy stellar mass function, the faint end of the main sequence, and the ability to track black holes in dwarf galaxies. SIMBA-C also matches better with recent observations of the mass–metallicity relation at z = 0, 2. By not assuming instantaneous recycling, SIMBA-C provides a much better match to galactic abundance ratio measures such as [O/Fe] and [N/O]. SIMBA-C thus opens up new avenues to constrain feedback models using detailed chemical abundance measures across cosmic time.

Missing Giants: Predictions on Dust-obscured Galaxy Stellar Mass Assembly Throughout Cosmic Time

Due to their extremely dust-obscured nature, much uncertainty still exists surrounding the stellar mass growth and content in dusty, star-forming galaxies (DSFGs) at z > 1. In this work, we present a numerical model built using empirical data on DSFGs to estimate their stellar mass contributions across the first ∼10 Gyr of cosmic time. We generate a dust-obscured stellar mass function that extends beyond the mass limit of star-forming stellar mass functions in the literature, and predict that massive DSFGs constitute as much as 50%–100% of all star-forming galaxies with M ≥1011 M ⊙ at z > 1. We predict the number density of massive DSFGs and find general agreement with observations, although more data is needed to narrow wide observational uncertainties. We forward-model mock massive DSFGs to their quiescent descendants and find remarkable agreement with observations from the literature demonstrating that, to first order, massive DSFGs are a sufficient ancestral population to describe the prevalence of massive quiescent galaxies at z > 1. We predict that massive DSFGs and their descendants contribute as much as 25%–60% to the cosmic stellar mass density during the peak of cosmic star formation, and predict an intense epoch of population growth during the ∼1 Gyr from z = 6 to 3 during which the majority of the most massive galaxies at high-z grow and then quench. Future studies seeking to understand massive galaxy growth and evolution in the early universe should strategize synergies with data from the latest observatories (e.g., JWST and the Atacama Large Millimeter/submillimeter Array) to better include the heavily dust-obscured galaxy population.

A Spatially Resolved Analysis of Star Formation Burstiness by Comparing UV and Hα in Galaxies at z ∼ 1 with UVCANDELS

The UltraViolet imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) program provides Hubble Space Telescope (HST)/UVIS F275W imaging for four CANDELS fields. We combine this UV imaging with existing HST/near-IR grism spectroscopy from 3D-HST+AGHAST to directly compare the resolved rest-frame UV and Hα emission for a sample of 979 galaxies at 0.7 < z < 1.5, spanning a range in stellar mass of 108−11.5 M ⊙. Using a stacking analysis, we perform a resolved comparison between homogenized maps of rest-UV and Hα to compute the average UV-to-Hα luminosity ratio (an indicator of burstiness in star formation) as a function of galactocentric radius. We find that galaxies below stellar mass of ∼109.5 M ⊙, at all radii, have a UV-to-Hα ratio higher than the equilibrium value expected from constant star formation, indicating a significant contribution from bursty star formation. Even for galaxies with stellar mass ≳109.5 M ⊙, the UV-to-Hα ratio is elevated toward their outskirts (R/R eff > 1.5), suggesting that bursty star formation is likely prevalent in the outskirts of even the most massive galaxies, but is likely overshadowed by their brighter cores. Furthermore, we present the UV-to-Hα ratio as a function of galaxy surface brightness, a proxy for stellar mass surface density, and find that regions below ∼107.5 M ⊙ kpc−2 are consistent with bursty star formation, regardless of their galaxy stellar mass, potentially suggesting that local star formation is independent of global galaxy properties at the smallest scales. Last, we find galaxies at z > 1.1 to have bursty star formation, regardless of radius or surface brightness.