SEARCH · Engineering Papers

Results for “cosmological hydrodynamical simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Cosmological Hydrodynamics at Exascale: A Trillion-Particle Leap in Capability

Resolving the most fundamental questions in cosmology requires simulations that match the scale, fidelity, and physical complexity demanded by next-generation sky surveys. To achieve the realism needed for this critical scientific partnership, detailed gas dynamics must be treated self-consistently with gravity for end-to-end modeling of structure formation. Exascale computing enables simulations that span survey-scale volumes while incorporating key astrophysical processes that shape complex cosmic structures. We present results from CRK-HACC, a cosmological hydrodynamics code built for extreme scalability. Using separation-of-scale techniques, GPU-resident tree solvers, in situ analysis pipelines, and multi-tiered I/O, CRK-HACCexecuted Frontier-E: a four trillion particle full-sky simulation, over an order of magnitude larger than previous efforts. The run achieved 513.1 PFLOPs peak performance, processing 46.6 billion particles per second and writing more than 100 PB of data in just over one week of runtime. Frontier-E marks a significant advance in predictive modeling for next-generation cosmological science.

Frontiere, Nicholas [Argonne National Laboratory (↗

DOE OSTI2025

A Gigaparsec-scale Hydrodynamic Volume Reconstructed with Deep Learning

The next generation of spectroscopic surveys will map the large-scale structure of the Universe at high redshifts (2 ≤ z ≤ 5) using millions of quasar spectra, enabling major advances in constraining both the standard cosmological model and its extensions. Robust cosmological analyses of these data sets require numerical simulations that both cover gigaparsec volumes and resolve features on ∼10 kpc scales and smaller. However, running such large-volume, high-resolution hydrodynamic simulations is computationally prohibitive. We present a generative deep learning model that enhances a low-resolution, gigaparsec-scale (960 h −1 Mpc) hydrodynamic simulation using a smaller (80 h −1 Mpc) high-resolution input hydrodynamic simulation as training data. The resulting enhanced simulation reproduces the line-of-sight power spectrum to within ∼10% and the three-dimensional power spectrum at the ∼20% level at intermediate to small scales (k ≲ 2 h Mpc −1 ). Our method shows strong promise for producing realistic simulations for cosmological analyses with current surveys such as the Dark Energy Spectroscopic Instrument and upcoming next-generation experiments, but further improvements are needed to accurately recover the large-scale modes. We publicly release the enhanced hydrodynamic simulation, along with a halo catalog from a companion N-body dark matter simulation to support the calibration of data analysis pipelines for these large-scale surveys.

Convolutional neural networks↗

NASA NTRS2014

Kinematic Evolution of Simulated Star-Forming Galaxies

Recent observations have shown that star-forming galaxies like our own Milky Way evolve kinematically into ordered thin disks over the last approximately 8 billion years since z = 1.2, undergoing a process of "disk settling." For the first time, we study the kinematic evolution of a suite of four state of the art "zoom in" hydrodynamic simulations of galaxy formation and evolution in a fully cosmological context and compare with these observations. Until now, robust measurements of the internal kinematics of simulated galaxies were lacking as the simulations suffered from low resolution, overproduction of stars, and overly massive bulges. The current generation of simulations has made great progress in overcoming these difficulties and is ready for a kinematic analysis. We show that simulated galaxies follow the same kinematic trends as real galaxies: they progressively decrease in disordered motions (sigma(sub g)) and increase in ordered rotation (V(sub rot)) with time. The slopes of the relations between both sigma(sub g) and V(sub rot) with redshift are consistent between the simulations and the observations. In addition, the morphologies of the simulated galaxies become less disturbed with time, also consistent with observations. This match between the simulated and observed trends is a significant success for the current generation of simulations, and a first step in determining the physical processes behind disk settling.

Evolution↗

DOE OSTI2020

Generating synthetic cosmological data with GalSampler

As part of the effort to meet the needs of the Large Synoptic Survey Telescope Dark Energy Science Collaboration (LSST DESC) for accurate, realistically complex mock galaxy catalogues, in this work we have developed galsampler, an open-source python package that assists in generating large volumes of synthetic cosmological data. The key idea behind galsampler is to recast hydrodynamical simulations and semi-analytic models as physically motivated galaxy libraries. galsampler populates a new, larger volume halo catalogue with galaxies drawn from the baseline library; by using weighted sampling guided by empirical modelling techniques, galsampler inherits statistical accuracy from the empirical model and physically motivated complexity from the baseline library. We have recently used galsampler to produce the cosmoDC 2 extragalactic catalogue made for the LSST DESC Data Challenge 2 . Using cosmoDC 2 as a guiding example, we outline how galsampler can continue to support ongoing and near-future galaxy surveys such as the Dark Energy Survey, the Dark Energy Spectroscopic Instrument, WFIRST, and Euclid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

The mass and galaxy distribution around SZ-selected clusters

ABSTRACT We present measurements of the radial profiles of the mass and galaxy number density around Sunyaev–Zel’dovich (SZ)-selected clusters using both weak lensing and galaxy counts. The clusters are selected from the Atacama Cosmology Telescope Data Release 5 and the galaxies from the Dark Energy Survey Year 3 data set. With signal-to-noise ratio of 62 (45) for galaxy (weak lensing) profiles over scales of about 0.2–20 h−1 Mpc, these are the highest precision measurements for SZ-selected clusters to date. Because SZ selection closely approximates mass selection, these measurements enable several tests of theoretical models of the mass and light distribution around clusters. Our main findings are: (1) The splashback feature is detected at a consistent location in both the mass and galaxy profiles and its location is consistent with predictions of cold dark matter N-body simulations. (2) The full mass profile is also consistent with the simulations. (3) The shapes of the galaxy and lensing profiles are remarkably similar for our sample over the entire range of scales, from well inside the cluster halo to the quasilinear regime. We measure the dependence of the profile shapes on the galaxy sample, redshift, and cluster mass. We extend the Diemer & Kravtsov model for the cluster profiles to the linear regime using perturbation theory and show that it provides a good match to the measured profiles. We also compare the measured profiles to predictions of the standard halo model and simulations that include hydrodynamics. Applications of these results to cluster mass estimation, cosmology, and astrophysics are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Emulating the Lyman-Alpha forest 1D power spectrum from cosmological simulations: new models and constraints from the eBOSS measurement

We present the Lyssa suite of high-resolution cosmological simulations of the Lyman-α forest designed for cosmological analyses. These 18 simulations have been run using the Nyx code with 40963 hydrodynamical cells in a 120 Mpc (∼ 81 Mpc/h) comoving box and individually provide sub-percent level convergence of the Lyman-α forest 1d flux power spectrum. We build a Gaussian process emulator for the Lyssa simulations in the lym1d likelihood framework to interpolate the power spectrum at arbitrary parameter values. We validate this emulator based on leave-one-out tests and based on the parameter constraints for simulations outside of the training set. We also perform comparisons with a previous emulator, showing a percent level accuracy and a good recovery of the expected cosmological parameters. Using this emulator we derive constraints on the linear matter power spectrum amplitude and slope parameters A Lyα and n Lyα . While the best-fit Planck ΛCDM model has A Lyα = 8.79 and n Lyα = -2.363, from DR14 eBOSS data we find that A Lyα < 7.6 (95% CI) and n Lyα = -2.369 ± 0.008. The low value of A Lyα , in tension with Planck, is driven by the correlation of this parameter with the mean transmission of the Lyman-α forest. This tension disappears when imposing a well-motivated external prior on this mean transmission, in which case we find A Lyα = 9.8 ± 1.1 in accordance with Planck.

Walther, Michael↗

DOE OSTI2026

Modeling the Cosmological Lyman-𝛼 Forest at the Field Level

The distribution of absorption lines in the spectra of distant quasars, called the Lyman-𝛼 (Ly-𝛼) forest, is a unique probe of cosmology and the intergalactic medium at high redshifts and small scales. The statistical power of ongoing redshift surveys demands precise theoretical tools to model the Ly-𝛼 forest. We address this challenge by developing an analytic, perturbative forward model to predict the Ly-𝛼 forest at the field level for a given set of cosmological initial conditions. Our model shows a remarkable performance when compared with the Sherwood hydrodynamic simulations: it reproduces the Ly-𝛼 forest flux power spectrum, its cross-correlation with dark matter halos, and the one-point probability distribution function of both fields at the percent level down to scales of a few Mpc. Our work provides crucial tools that bridge analytic modeling on large scales with simulations on small scales, enabling field-level inference from Ly-𝛼 forest data and simulation-based priors for cosmological analyses. Furthermore, this is especially timely for realizing the full scientific potential of the Ly-𝛼 forest measurements by the dark energy spectroscopic instrument.

Cosmological parameters↗

DOE OSTI2023

Compressing the Cosmological Information in One-dimensional Correlations of the Lyman-$\alpha$ Forest

Observations of the Lyman-α forest from spectroscopic surveys such as the Baryon Oscillation Spectroscopic Survey or its extension, eBOSS, or the ongoing Dark Energy Spectroscopic Instrument (DESI) survey offer a unique window to study the growth of structure on megaparsec scales. Interpretation of these measurements is a complicated task, requiring hydrodynamical simulations to model and marginalize over the thermal and ionization state of the intergalactic medium. This complexity has limited the use of Lyα clustering measurements in joint cosmological analyses. In this work we show that the cosmological information content of the one-dimensional power spectrum (P 1D ) of the Lyman-α forest can be compressed into a simple two-parameter likelihood without any significant loss of constraining power. We simulate P 1D measurements from DESI using hydrodynamical simulations and show that the compressed likelihood is model independent and lossless, recovering unbiased results even in the presence of massive neutrinos or running of the primordial power spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Weak lensing combined with the kinetic Sunyaev–Zel’dovich effect: a study of baryonic feedback

ABSTRACT Extracting precise cosmology from weak lensing surveys requires modelling the non-linear matter power spectrum, which is suppressed at small scales due to baryonic feedback processes. However, hydrodynamical galaxy formation simulations make widely varying predictions for the amplitude and extent of this effect. We use measurements of Dark Energy Survey Year 3 weak lensing (WL) and Atacama Cosmology Telescope DR5 kinematic Sunyaev–Zel’dovich (kSZ) to jointly constrain cosmological and astrophysical baryonic feedback parameters using a flexible analytical model, ‘baryonification’. First, using WL only, we compare the $S_8$ constraints using baryonification to a simulation-calibrated halo model, a simulation-based emulator model, and the approach of discarding WL measurements on small angular scales. We find that model flexibility can shift the value of $S_8$ and degrade the uncertainty. The kSZ provides additional constraints on the astrophysical parameters, with the joint WL + kSZ analysis constraining $S_8=0.823^{+0.019}_{-0.020}$. We measure the suppression of the non-linear matter power spectrum using WL + kSZ and constrain a mean feedback scenario that is more extreme than the predictions from most hydrodynamical simulations. We constrain the baryon fractions and the gas mass fractions and find them to be generally lower than inferred from X-ray observations and simulation predictions. We conclude that the WL + kSZ measurements provide a new and complementary benchmark for building a coherent picture of the impact of gas around galaxies across observations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

The impact of varying inhomogeneous reionization histories on metrics of Ly α opacity

The epoch of hydrogen reionization is complete by z = 5⁠, but its progression at higher redshifts is uncertain. Measurements of Ly α forest opacity show large scatter at z < 6⁠, suggestive of spatial fluctuations in neutral fraction, temperature, or ionizing background, either individually or in combination. However, there are degeneracies in the impact of such fluctuations, necessitating careful modelling. We develop a framework for modelling the reionization history and associated temperature fluctuations, with the intention of incorporating ionizing background fluctuations at a later time. We generate several reionization histories using seminumerical code AMBER, and implement them in the Nyx cosmological hydrodynamics code to examine the impact on the evolution of gas within the simulation and the associated metrics of the Ly α forest opacity. We find that the pressure smoothing scale within the intergalactic medium is strongly correlated with the adiabatic index of the temperature–density relation. We find that while models with 20 000 K photoheating at reionization are better able to reproduce the shape of the observed z = 5 1D flux power spectrum than colder ones, they fail to match the highest wavenumbers. The simulated autocorrelation function and optical depth distributions are systematically low and narrow, respectively, compared to the observed values, but are in better agreement when the reionization history is longer in duration, more symmetric in its distribution of reionization redshifts, or if there are remaining neutral regions at z < 6.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

The AGORA High-Resolution Galaxy Simulations Comparison Project

Context. Satellite galaxies experience multiple physical processes when interacting with their host halos, often leading to the quenching of star formation. In the Local Group, satellite quenching has been shown to be highly efficient, affecting nearly all satellites except the most massive ones. While recent surveys study Milky Way-analogs to assess how representative our Local Group is, the dominant physical mechanisms behind satellite quenching in Milky Way-mass halos remain under debate. Aims. We analyze satellite quenching within the same Milky Way-mass halo simulated using various widely used astrophysical codes, each using different hydrodynamic methods and implementing different supernovae feedback recipes. The goal is to determine whether quenched fractions, quenching timescales, and the dominant quenching mechanisms are consistent across codes or if they show sensitivity to the specific hydrodynamic method and supernovae feedback physics employed. Methods. We used a subset of high-resolution cosmological zoom-in simulations of a Milky Way-mass halo from the multiple-code AGORA CosmoRun suite. Our analysis focuses on comparing satellite quenching across the different models and against observational data. We also analyzed the dominant mechanisms driving satellite quenching in each model. Results. We find that the quenched fraction is consistent with the latest SAGA Survey results within its 1σ host-to-host scatter across all the models. Regarding quenching timescales, all the models reproduce the trend observed in the ELVES survey, Local Group observations, and previous simulations: The less massive the satellite, the shorter its quenching timescale. All of our models converge on the dominant quenching mechanisms: Strangulation halts cold gas accretion in all satellites, while ram pressure stripping is the predominant mechanism for gas removal, and it is particularly effective in satellites with M * <10 8 M ⊙ . Nevertheless, the efficiency of the stripping mechanisms differs among the codes, showing a strong sensitivity to the different supernovae feedback implementations and/or hydrodynamic methods employed.

Local Group↗

DOE OSTI2025

Merger Response of Halo Anisotropy Properties

Anisotropy properties—halo spin, shape, position offset, velocity offset, and orientation—are an important family of dark matter halo properties that indicate the level of directional variation of the internal structures of halos. These properties reflect the dynamical state of halos, which in turn depends on the mass assembly history. In this work, we study the evolution of anisotropy properties in response to merger activity using the IllustrisTNG simulations. We find that the response trajectories of the anisotropy properties significantly deviate from secular evolution. These trajectories have the same qualitative features and timescales across a wide range of merger and host properties. We propose explanations for the behavior of these properties and connect their evolution to the relevant stages of merger dynamics. We measure the relevant dynamical timescales. We also explore the dependence of the strength of the response on time of merger, merger ratio, and mass of the main halo. These results provide insight into the physics of halo mergers and their effects on the statistical behavior of halo properties. This study paves the way toward a physical understanding of scaling relations, particularly to how systematics in their scatter are connected to the mass assembly histories of halos.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Demonstrating the viability of Lagrangian in situ reduction on supercomputers

Performing exploratory analysis and visualization of large-scale time-varying computational science applications is challenging due to inaccuracies that arise from under-resolved data. In recent years, Lagrangian representations of the vector field computed using in situ processing are being increasingly researched and have emerged as a potential solution to enable exploration. However, prior works have offered limited estimates of the encumbrance on the simulation code as they consider “theoretical” in situ environments. Further, the effectiveness of this approach varies based on the nature of the vector field, benefitting from an in-depth investigation for each application area. With this study, an extended version of Sane et al. (2021), we contribute an evaluation of Lagrangian analysis viability and efficacy for simulation codes executing at scale on a supercomputer. We investigated previously unexplored cosmology and seismology applications as well as conducted a performance benchmarking study by using a hydrodynamics mini-application targeting exascale computing. Here, to inform encumbrance, we integrated in situ infrastructure with simulation codes, and evaluated Lagrangian in situ reduction in representative homogeneous and heterogeneous HPC environments. To inform post hoc accuracy, we conducted a statistical analysis across a range of spatiotemporal configurations as well as a qualitative evaluation. Additionally, our study contributes cost estimates for distributed-memory post hoc reconstruction. In all, we demonstrate viability for each application — data reduction to less than 1% of the total data via Lagrangian representations, while maintaining accurate reconstruction and requiring under 10% of total execution time in over 90% of our experiments.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

The aluminium-26 distribution in a cosmological simulation of a Milky Way-type Galaxy

Context . The 1.8 MeV γ -rays corresponding to the decay of the radioactive isotope 26 Al (with a half-life of 0.72 Myr ) have been observed by the SPI detector on the INTEGRAL spacecraft and extensively used as a tracer of star formation and current nucleosynthetic activity in the Milky Way Galaxy. Further information is encoded in the observation related to the higher 26 Al content found in regions of the Galaxy with the highest line-of-sight (LoS) velocity relative to an observer located in the Solar System. However, this feature remains unexplained. Aims . We ran a cosmological “zoom-in” chemodynamical simulation of a Milky Way-type galaxy, including the production and decays of radioactive nuclei in a fully self-consistent way. We then analyzed the results to follow the evolution of 26 Al throughout the lifetime of the simulated galaxy to provide a new method for interpreting the 26 Al observations. Methods . We included the massive star sources of 26 Al in the Galaxy and its radioactive decay into a state-of-the-art galactic chemical evolution model, coupled with cosmological growth and hydrodynamics. This approach allowed us to follow the spatial and temporal evolution of the 26 Al content in the simulated galaxy. Results . Our results are in agreement with the observations with respect to the fact that gas particles in the simulation with relatively higher 26 Al content also have the highest LoS velocities. On the other hand, gas particles with relatively lower 26 Al content (i.e., not bright enough to be observed) generally display the lowest LoS velocities. However, this result is not conclusive because the overall rotational velocity of our simulated galaxy is higher than that observed for cold CO gas in the Milky Way Galaxy. Furthermore, we found no significant correlation between gas temperature, rotational velocity, and 26 Al content at any given radius. We also found the presence of transient 26 Al-rich spots at low LoS velocities and we show that one such spot had been captured by the INTEGRAL/SPI data. Based on our model, we present a prediction for the detection of 1.8 MeV γ -rays by the future COSI mission. We find that according to our model, the new instrument will be able to observe similar 26 Al-emission patterns to those seen by INTEGRAL/SPI.

Wehmeyer, B. (ORCID:0009000254149292)↗

DOE OSTI2024

The AGORA High-resolution Galaxy Simulations Comparison Project. IV. Halo and Galaxy Mass Assembly in a Cosmological Zoom-in Simulation at z ≤ 2

In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift z = 2 and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at z = 4 and 3, and before the last major merger, focusing on the formation of well-defined rotationally supported disks, the mass–metallicity relation, the specific star formation rate, the gas metallicity gradients, and the nonaxisymmetric structures in the stellar disks. Codes generally converge well to the stellar-to-halo mass ratios predicted by semianalytic models at z ~ 2. We see that almost all the hydro codes develop rotationally supported structures at low redshifts. Most agree within 0.5 dex with the observed mass–metallicity relation at high and intermediate redshifts, and reproduce the gas metallicity gradients obtained from analytical models and low-redshift observations. We confirm that the intercode differences in the halo assembly history reported in the first paper of the collaboration also exist in CosmoRun, making the code-to-code comparison more difficult. We show that such differences are mainly due to variations in code-dependent parameters that control the time stepping strategy of the gravity solver. We find that variations in the early stellar feedback can also result in differences in the timing of the low-redshift mergers. All the simulation data down to z = 2 and the auxiliary data will be made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1997

Gamma-Ray Bursts and Relativistic MHD around Coalescing Neutron Stars

In this proposal we performed the first relativistic hydrodynamic simulations of the binary neutron stars near coalescence. Calculations were performed in three spatial dimensions and were supplemented with magnetohydrodynamic simulations of the magnetic reconnection as the neutron star fluid responds to the relativistic forces of the binary. The ultimate goal of the proposed work was to develop a theoretical model for the temporal and spectral evolution of cosmological gamma-ray bursts produced during the magnetic reconnection. This is the first time such hydrodynamic calculations were performed and some surprising results were obtained. The neutrons stars were observed to compress and heat in the simulations. They could even collapse to black holes many orbits before coalescence. The question then arose as to whether this strange phenomenon could power a gamma-ray burst. The bulk of computing effort in this proposal was therefore devoted to an analysis of the implications of this possible new gamma-ray burst paradigm. The initial focus of the work performed under this proposal was therefore to utilize (3+1) dimensional and spherical numerical general relativistic hydrodynamics to study the origin, evolution, and parametric sensitivity of the collapse instability. We derived modified conditions of hydrostatic equilibrium for stars in the curved space of quasi-static orbits. We followed the magnetic field evolution for approximately one orbit (10 msec). The magnetic field energy was observed to exponentiate with an e-folding time of about 1 msec. Thus, the field should build up very quickly to a magnitude such that reconnection and back reaction of the fluid inhibits further growth.

Mathews, Grant J.↗

DOE OSTI2022

Combining Planck and SPT Cluster Catalogs: Cosmological Analysis and Impact on the Planck Scaling Relation Calibration

We provide the first combined cosmological analysis of the South Pole Telescope (SPT) and Planck cluster catalogs. The aim is to provide an independent calibration for Planck scaling relations, exploiting the cosmological constraining power of the SPT-SZ cluster catalog and its dedicated weak lensing (WL) and X-ray follow-up observations. We build a new version of the Planck cluster likelihood. In the νΛ CDM scenario, focusing on the mass slope and mass bias of Planck scaling relations, we find ${\alpha }_{\mathrm{SZ}}={1.49}_{-0.10}^{+0.07}$ and ${\left(1-b\right)}_{\mathrm{SZ}}={0.69}_{-0.14}^{+0.07}$, respectively. The results for the mass slope show a ~4 σ departure from the self-similar evolution, α SZ ~ 1.8. This shift is mainly driven by the matter density value preferred by SPT data, Ωm = 0.30 ± 0.03, lower than the one obtained by Planck data alone, ${{\rm{\Omega }}}_{m}={0.37}_{-0.06}^{+0.02}$. The mass bias constraints are consistent both with outcomes of hydrodynamical simulations and external WL calibrations, (1 – b) ~ 0.8, and with results required by the Planck cosmic microwave background cosmology, (1 – b) ~ 0.6. From this analysis, we obtain a new catalog of Planck cluster masses M 500 . We estimate the ratio between the published Planck M SZ masses and our derived masses M 500 , as a "measured mass bias," ${\left(1-b\right)}_{M}$. We analyze the mass, redshift, and detection noise dependence of ${\left(1-b\right)}_{M}$, finding an increasing trend toward high redshift and low mass. These results mimic the effect of departure from self-similarity in cluster evolution, showing different dependencies for the low-mass, high-mass, low-z, and high-z regimes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Dwarf Galaxies at Cosmic Noon: New JWST Constraints on Satellite Models and Subhalo Tidal Evolution

The advent of JWST has revolutionized the study of faint satellite galaxies at z ≳ 1, enabling statistical constraints on galaxy evolution and the galaxy–halo connection in a previously unexplored mass and redshift regime. We compare satellite abundances at 1 < z < 3.5 from recent JWST observations with predictions from cosmological dark-matter-only zoom-in simulations. We identify and quantify several sources of biases that can impact theoretical satellite counts, finding that assumptions about subhalo tidal evolution introduce the largest uncertainty in predictions for the satellite mass function. Using a flexible galaxy disruption model, we explore a range of disruption scenarios, spanning hydrodynamically motivated and idealized prescriptions, to bracket plausible physical outcomes. We show that varying galaxy durability can change the predicted satellite mass functions by a factor of ∼3.5. The JWST data and our fiducial model are consistent within 1σ–2σ across the full redshift (1 < z < 3.5) and stellar mass (M⋆ > 10 7 M ⊙ ) range probed. We find evidence that subhalos are at least as long-lived as predicted by hydrodynamic simulations. Our framework will enable robust constraints on the tidal evolution of subhalos with future observations. This work presents the first direct comparison between cosmological models and observations of the high-redshift satellite population in this low-mass regime. These results showcase JWST’s emerging power to test structure formation in the first half of the Universe in a new domain and to constrain the physical processes driving the evolution of low-mass galaxies across cosmic time.

79 ASTRONOMY AND ASTROPHYSICS↗