Search

The James Webb Space Telescope (JWST) is capable of probing extremely early eras of our Universe, when the supersonic relative motions between dark matter and baryonic overdensities modulate structure formation ( z ≳ 10). We study low-mass galaxy formation, including this “stream velocity,” using high-resolution AREPO hydrodynamics simulations and present theoretical predictions of the UV luminosity function (UVLF) and galaxy stellar mass function down to extremely faint and low-mass galaxies ( M _UV ≳ −15, 10 ^4 M _⊙ ≤ M _* ≤ 10 ^8 M _⊙ ). We show that, although the stream velocity suppresses early star formation overall, it induces a short period of rapid star formation in some larger dwarfs, leading to an enhancement in the faint end of the UVLF at z = 12. We demonstrate that JWST observations are close to this enhanced regime and propose that the UVLF may constitute an important probe of the stream velocity at high redshift for JWST and future observatories.

We study the stellar properties of a sample of simulated ultradiffuse galaxies (UDGs) with stellar mass M _⋆ = 10 ^7.5 –10 ^9 M _⊙ , selected from the TNG50 simulation, where UDGs form mainly in high-spin dwarf-mass halos. We divide our sample into star-forming and quenched UDGs, finding good agreement with the stellar assembly history measured in observations. Star-forming UDGs and quenched UDGs with M _⋆ ≥ 10 ^8 M _⊙ in our sample are particularly inefficient at forming stars, having 2–10 times less stellar mass than non-UDGs for the same virial mass halo. These results are consistent with recent mass inferences in UDG samples and suggest that the most inefficient UDGs arise from a late assembly of the dark matter mass followed by a stellar growth that is comparatively slower (for star-forming UDGs) or that was interrupted due to environmental removal of the gas (for quenched UDGs). Regardless of efficiency, UDGs are 60% poorer in [Fe/H] than the population of non-UDGs at a fixed stellar mass, with the most extreme objects having metal content consistent with the simulated mass–metallicity relation at z ∼ 2. Quenched UDGs stop their star formation in shorter timescales than non-UDGs of similar mass and are, as a consequence, alpha enhanced with respect to non-UDGs. We identify metallicity profiles in UDGs as a potential avenue to distinguish between different formation paths for these galaxies, where gentle formation as a result of high-spin halos would present well-defined declining metallicity radial profiles while powerful-outflows or tidal stripping formation models would lead to flatter or constant metallicity as a function of radius due to the inherent mixing of stellar orbits.

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 SFH burstiness. We compute the mock SEDs for these model galaxies at z = 7 via calculations of 3D dust radiative transfer, and then backward 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 priors for SFH to model star formation in the early Universe.

Milky Way-type galaxies are surrounded by a warm-hot gaseous halo containing a considerable amount of baryons and metals. The kinematics and spatial distribution of highly ionized ion species such as O vi can be significantly affected by supernova (SN) explosions and early (pre-SN) stellar feedback (e.g., stellar winds, radiation pressure). Here we investigate effects of stellar feedback on O vi absorptions in Milky Way−like galaxies by analyzing the suites of high-resolution hydrodynamical simulations under the framework of SMUGGLE, a physically motivated subgrid interstellar medium and stellar feedback model for the moving-mesh code Arepo . We find that the fiducial run with the full suite of stellar feedback and moderate star formation activities can reasonably reproduce Galactic O vi absorptions observed by space telescopes such as the Far-Ultraviolet Spectroscopic Explorer, including the scale height of low-velocity (∣ v _LSR ∣ < 100 km s ^−1 ) O vi , the column density–line width relation for high-velocity (100 km s ^−1 ≤ ∣ v _LSR ∣ < 400 km s ^−1 ) O vi , and the cumulative O vi column densities. In contrast, model variations with more intense star formation activities deviate from observations further. Additionally, we find that the run considering only SN feedback is in broad agreement with the observations, whereas in runs without SN feedback this agreement is absent, which indicates a dominant role of SN feedback in heating and accelerating interstellar O vi . This is consistent with the current picture that interstellar O vi is predominantly produced by collisional ionization where mechanical feedback can play a central role. In contrast, photoionization is negligible for O vi production owing to the lack of high-energy (≳114 eV) photons required.

We explore the chemodynamical properties of the Galaxy in the azimuthal velocity V _ϕ and metallicity [Fe/H] space using red giant stars from Gaia Data Release 3. The row-normalized V _ϕ –[Fe/H] maps form a coherent sequence from the bulge to the outer disk, clearly revealing the thin/thick disk and the Splash. The metal-rich stars display bar-like kinematics, while the metal-poor stars show dispersion-dominated kinematics. The intermediate-metallicity population (−1 < [Fe/H]< − 0.4) can be separated into two populations, one that is bar-like, i.e., dynamically cold ( ${\sigma }_{{V}_{R}}\sim 80$ km s ^−1 ) and fast-rotating ( V _ϕ ≳ 100 km s ^−1 ), and the Splash, which is dynamically hot ( ${\sigma }_{{V}_{R}}\sim 110$ km s ^−1 ) and slow-rotating ( V _ϕ ≲ 100 km s ^−1 ). We compare the observations in the bulge region with an Auriga simulation where the last major merger event occurred ∼10 Gyr ago: only stars born around the time of the merger reveal a Splash-like feature in the V _ϕ –[Fe/H] space, suggesting that the Splash is likely merger-induced, predominantly made up of heated disk stars and the starburst associated with the last major merger. Since the Splash formed from the proto-disk, its lower metallicity limit coincides with that of the thick disk. The bar formed later from the dynamically hot disk with [Fe/H] > − 1 dex, with the Splash not participating in the bar formation and growth. Moreover, with a set of isolated evolving N -body disk simulations, we confirm that a nonrotating classical bulge can be spun up by the bar and develop cylindrical rotation, consistent with the observations for the metal-poor stars.

At high redshifts ( z ≳ 12), the relative velocity between baryons and dark matter (the so-called streaming velocity) significantly affects star formation in low-mass objects. Streaming substantially reduces the abundance of low-mass gas objects while simultaneously allowing for the formation of supersonically induced gas objects (SIGOs) and their associated star clusters outside of dark matter halos. Here, we present a study of the population-level effects of streaming on star formation within both halos and SIGOs in a set of simulations with and without streaming. Notably, we find that streaming actually enhances star formation within individual halos of all masses at redshifts between z = 12 and z = 20. This is demonstrated both as an increased star formation rate per object as well as an enhancement of the Kennicutt–Schmidt relation for objects with streaming. We find that our simulations are consistent with some observations at high redshift, but on a population level, they continue to underpredict star formation relative to the majority of observations. Notably, our simulations do not include feedback and so can be taken as an upper limit on the star formation rate, exacerbating these differences. However, simulations of overdense regions (both with and without streaming) agree with observations, suggesting a strategy for extracting information about the overdensity and streaming velocity in a given survey volume in future observations.

Simulations of galaxy formation are mostly unable to resolve the energy-conserving phase of individual supernova events, having to resort to subgrid models to distribute the energy and momentum resulting from stellar feedback. However, the properties of these simulated galaxies, including the morphology, stellar mass formed, and the burstiness of the star formation history, are highly sensitive to the numerical choices adopted in these subgrid models. Using the SMUGGLE stellar feedback model, we carry out idealized simulations of an M _vir ∼ 10 ^10 M _⊙ dwarf galaxy, a regime where most simulation codes predict significant burstiness in star formation, resulting in strong gas flows that lead to the formation of dark matter cores. We find that by varying only the directional distribution of momentum imparted from supernovae to the surrounding gas, while holding the total momentum per supernova constant, bursty star formation may be amplified or completely suppressed, and the total stellar mass formed can vary by as much as a factor of ∼3. In particular, when momentum is primarily directed perpendicular to the gas disk, less bursty and lower overall star formation rates result, yielding less gas turbulence, more disky morphologies, and a retention of cuspy dark matter density profiles. An improved understanding of the nonlinear coupling of stellar feedback into inhomogeneous gaseous media is thus needed to make robust predictions for stellar morphologies and dark matter core formation in dwarfs independent of uncertain numerical choices in the baryonic treatment.

The formation mechanism of globular clusters (GCs) has long been debated by astronomers. It was recently proposed that supersonically induced gas objects (SIGOs)–which formed in the early Universe due to the supersonic relative motion of baryons and dark matter at recombination–could be the progenitors of early GCs. In order to become GCs, SIGOs must form stars relatively efficiently despite forming outside of dark matter halos. We investigate the potential for star formation in SIGOs using cosmological hydrodynamic simulations, including the aforementioned relative motions of baryons and dark matter, molecular hydrogen cooling in primordial gas clouds, and explicit star formation. We find that SIGOs do form stars and that the nascent star clusters formed through this process are accreted by dark matter halos on short timescales (∼a few hundred megayears). Thus, SIGOs may be found as intact substructures within these halos, analogous to many present-day GCs. From this result, we conclude that SIGOs are capable of forming star clusters with similar properties to globular clusters in the early Universe, and we discuss their detectability by upcoming JWST surveys.

Elongated bar-like features are ubiquitous in galaxies, occurring at the centers of approximately two-thirds of spiral disks in the nearby Universe. Due to gravitational interactions between the bar and the other components of galaxies, it is expected that angular momentum and matter will redistribute over long (Gyr) timescales in barred galaxies. Previous work ignoring the gas phase of galaxies has conclusively demonstrated that bars should slow their rotation over time due to their interaction with dark matter halos. We have performed a simulation of a Milky Way–like galactic disk hosting a strong bar, including a state-of-the-art model of the interstellar medium and a live dark matter halo. In this simulation, the bar pattern does not slow down over time, and instead it remains at a stable, constant rate of rotation. This behavior has been observed in previous simulations using more simplified models for the interstellar gas, but the apparent lack of secular evolution has remained unexplained. We find that the presence of the gas phase arrests the process by which the dark matter halo slows down a bar, a phenomenon we term bar locking. This locking is responsible for stabilizing the bar pattern speed. We find that, in a Milky Way–like disk, a gas fraction of only about 5% is necessary for this mechanism to operate. Our result naturally explains why nearly all observed bars rotate rapidly and is especially relevant for our understanding of how the Milky Way arrived at its present state.

Supersonically induced gas objects (SIGOs) are a class of early universe objects that have gained attention as a potential formation route for globular clusters. SIGOs have recently begun to be studied in the context of molecular hydrogen cooling, which is key to characterizing their structure and evolution. Studying the population-level properties of SIGOs with molecular cooling is important for understanding their potential for collapse and star formation, and for addressing whether SIGOs can survive to the present epoch. Here, we investigate the evolution of SIGOs before they form stars, using a combination of numerical and analytical analysis. We study timescales important to the evolution of SIGOs at a population level in the presence of molecular cooling. Revising the previous formulation for the critical density of collapse for SIGOs allows us to show that their prolateness tends to act as an inhibiting factor to collapse. We find that simulated SIGOs are limited by artificial two-body relaxation effects that tend to disperse them. We expect that SIGOs in nature will be longer lived compared to our simulations. Further, the fall-back timescale on which SIGOs fall into nearby dark matter halos, potentially producing a globular-cluster-like system, is frequently longer than their cooling timescale and the collapse timescale on which they shrink through gravity. Therefore, some SIGOs have time to cool and collapse outside of halos despite initially failing to exceed the critical density. From this analysis we conclude that SIGOs should form stars outside of halos in nonnegligible stream velocity patches in the universe.

A supersonic relative velocity between dark matter (DM) and baryons (the stream velocity) at the time of recombination induces the formation of low-mass objects with anomalous properties in the early universe. We widen the scope of the “Supersonic Project” paper series to include objects we term Dark Matter + Gas Halos Offset by Streaming (DM GHOSts)—diffuse, DM-enriched structures formed because of a physical offset between the centers of mass of DM and baryonic overdensities. We present an updated numerical investigation of DM GHOSts and Supersonically Induced Gas Objects (SIGOs), including the effects of molecular cooling, in high-resolution hydrodynamic simulations using the AREPO code. Supplemented by an analytical understanding of their ellipsoidal gravitational potentials, we study the population-level properties of these objects, characterizing their morphology, spin, radial mass, and velocity distributions in comparison to classical structures in non-streaming regions. The stream velocity causes deviations from sphericity in both the gas and DM components and lends greater rotational support to the gas. Low-mass (≲10 ^5.5 M _⊙ ) objects in regions of streaming demonstrate core-like rotation and mass profiles. Anomalies in the rotation and morphology of DM GHOSts could represent an early universe analog to observed ultra-faint dwarf galaxies with variations in DM content and unusual rotation curves.

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 4233 cosmological simulations, 2049 N -body simulations, and 2184 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 catalogs, power spectra, bispectra, Ly α spectra, probability distribution functions, halo radial profiles, and X-rays photon lists. We also release over 1000 catalogs that contain billions of galaxies from CAMELS-SAM: a large collection of N -body simulations that have been combined with the Santa Cruz semianalytic 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 https://camels.readthedocs.io .

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 model for dust evolution in galaxy simulations. We apply these models to three 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 a starburst disk. Our main results are as follows. Galaxies with high specific star formation rates have increased feedback energy per unit mass, and are able to shatter grains efficiently, driving up the fraction of ultrasmall 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 investigate fully the complex interplay of processes that drive PAH band luminosities in galaxies.

We explore the role of galactic feedback on the low-redshift Ly α (Ly α ) forest ( z ≲ 2) statistics and its potential to alter the thermal state of the intergalactic medium. Using the Cosmology and Astrophysics with Machine Learning Simulations (CAMELS) suite, we explore variations of the AGN and stellar feedback models in the IllustrisTNG and Simba subgrid models. We find that both AGN and stellar feedback in Simba play a role in setting the Ly α forest column density distribution function (CDD) and the Doppler width ( b -value) distribution. The Simba AGN jet feedback mode is able to efficiently transport energy out to the diffuse IGM, causing changes in the shape and normalization of the CDD and a broadening of the b -value distribution. We find that stellar feedback plays a prominent role in regulating supermassive black hole growth and feedback, highlighting the importance of constraining stellar and AGN feedback simultaneously. In IllustrisTNG, the AGN feedback variations explored in CAMELS do not affect the Ly α forest, but varying the stellar feedback model does produce subtle changes. Our results imply that the low- z Ly α forest can be sensitive to changes in the ultraviolet background, stellar and black hole feedback, and that AGN jet feedback in particular can have a strong effect on the thermal state of the IGM.

Gas flows in the presence of two independently rotating nested bars remain not fully understood but are likely to play an important role in fueling the central black hole. We use high-resolution hydrodynamical simulations with detailed models of subgrid physics to study this problem. Our results show that the inner bar in double-barred galaxies can help drive gas flow from the nuclear ring to the center. In contrast, gas inflow usually stalls at the nuclear ring in single-barred galaxies. The inner bar causes a quasiperiodic inflow with a frequency determined by the difference between the two bar pattern speeds. We find that the star formation rate is higher in the model with two bars than in that with one bar. The inner bar in our model gradually weakens and dissolves due to gas inflow over a few billion years. Star formation produces metal-rich/ α -poor stars, which slows the weakening of the inner bar but does not halt its eventual decay. We also present a qualitative comparison of the gas morphology and kinematics in our simulations with those of observed double-barred galaxies.