The abundance and demographics of dark matter substructure is important for many areas in astrophysics and cosmological N-body simulations have been the primary tool used to investigate them. However, it has recently become clear that the simulations are subject to numerical artefacts, which hampers a proper treatment of the tidal evolution of subhaloes. Unfortunately, no analytical models that accurately describe subhalo evolution exist either. We therefore present a library of idealized, high resolution N-body simulations of the tidal evolution of individual subhaloes that can be used to calibrate semi-analytical models and to complement cosmological simulations. The simulations focus on minor mergers, i.e., the mass of the subhalo is much smaller than that of the host halo, such that the impact of dynamical friction is negligible. This setup allows the adoption of a fixed analytical potential for modelling the host halo. The dynamical evolution of subhaloes is followed with N-body computations. In the library, four parameters, two of which characterize the subhalo orbit with respect to the host halo, and the two concentrations of the host-and subhalo, are varied over the ranges encountered in cosmological simulations. We show several representative examples from the library that illustrate the evolution of the subhalo mass and velocity dispersion profiles. Additionally, we make publicly available a pre-trained nonparametric model of the subhalo mass evolution based on random forest regression. This model is able to interpolate the simulation data at the 0.1 dex level and provides efficient access to the data for further use in modelling †. (GO) † DASH simulation data and machine learning model are available at https://cosmo.oca.eu/dash.

Several recent studies have indicated that artificial subhalo disruption (the spontaneous, non-physical disintegration of a subhalo) remains prevalent in state-of-the-art dark matter-only cosmological simulations. In order to quantify the impact of disruption on the inferred subhalo demographics, we augment the semi-analytical SatGen dynamical subhalo evolution model with an improved treatment of tidal stripping that is calibrated using the DASH database of idealized high-resolution simulations of subhalo evolution, which are free from artificial disruption. We also develop a model of artificial disruption that reproduces the statistical properties of disruption in the Bolshoi simulation. Using this framework, we predict subhalo mass functions (SHMFs), number density profiles, and substructure mass fractions and study how these quantities are impacted by artificial disruption and mass resolution limits. We find that artificial disruption affects these quantities at the $10-20\%$ level, ameliorating previous concerns that it may suppress the SHMF by as much as a factor of two. We demonstrate that semi-analytical substructure modeling must include orbit integration in order to properly account for splashback haloes, which make up roughly half of the subhalo population. We show that the resolution limit of N-body simulations, rather than artificial disruption, is the primary cause of the radial bias in subhalo number density found in dark matter-only simulations. Hence, we conclude that the mass resolution remains the primary limitation of using such simulations to study subhaloes. Our model provides a fast, flexible, and accurate alternative to studying substructure statistics in the absence of both numerical resolution limits and artificial disruption.

We present a semi-analytical model of satellite galaxies, SatGen, which can generate large statistical samples of satellite populations for a host halo of desired mass, redshift, and assembly history. The model combines dark matter (DM) halo merger trees, empirical relations for the galaxy–halo connection, and analytical prescriptions for tidal effects, dynamical friction, and ram-pressure stripping. SatGen emulates cosmological zoom-in hydrosimulations in certain aspects. Satellites can reside in cored or cuspy DM subhaloes, depending on the halo response to baryonic physics that can be formulated from hydrosimulations and physical modelling. The subhalo profile and the stellar mass and size of a satellite evolve depending on its tidal mass-loss and initial structure. The host galaxy can include a baryonic disc and a stellar bulge, each described by a density profile that allows analytic satellite orbit integration. SatGen complements simulations by propagating the effect of halo response found in simulated field galaxies to satellites (not properly resolved in simulations) and outperforms simulations by sampling the halo-to-halo variance of satellite statistics and overcoming artificial disruption due to insufficient resolution. As a first application, we use the model to study satellites of Milky Way (MW)- and M31-sized hosts, making it emulate simulations of bursty star formation and of smooth star formation, respectively, and to experiment with a disc potential in the host halo. We find that our model reproduces the observed satellite statistics reasonably well. Different physical recipes make a difference in satellite abundance and spatial distribution at the 25 per cent level, not large enough to be distinguished by current observations given the halo-to-halo variance. The MW/M31 disc depletes satellites by ${\sim } 20{{\ \rm per\ cent}}$ and has a subtle effect of diversifying the internal structure of satellites, which is important for alleviating certain small-scale problems. We discuss the conditions for a massive satellite to survive in MW/M31.

We present a method called Significant Cosmic Holes in Universe (SCHU) for identifying cosmic voids and loops of filaments in cosmological datasets and assigning their statistical significance using techniques from topological data analysis. In particular, persistent homology is used to find different dimensional holes. For dark matter halo catalogs and galaxy surveys, the 0-, 1-, and 2dimensional holes can be identified with connected components (i.e. clusters), loops of filaments, and voids. The procedure overlays dark matter halos/galaxies on a three-dimensional grid, and a distance-to-measure (DTM) function is calculated at each point of the grid. A nested set of simplicial complexes (a filtration) is generated over the lower-level sets of the DTM across increasing threshold values. The filtered simplicial complex can then be used to summarize the birth and death times of the different dimension homology group generators (i.e., the holes). Persistent homology summary diagrams, called persistence diagrams, are produced from the dimension, birth times, and death times of each homology group generator. Using the persistence diagrams and bootstrap sampling, we explain how p-values can be assigned to each homology group generator. The homology group generators on a persistence diagram are not, in general, uniquely located back in the original dataset volume so we propose a method for finding a representation of the homology group generators. This method provides a novel, statistically rigorous approach for locating informative generators in cosmological datasets, which may be useful for providing complementary cosmological constraints on the effects of, for example, the sum of the neutrino masses. The method is tested on a Voronoi foam simulation, and then subsequently applied to a subset of the SDSS galaxy survey and a cosmological simulation. Lastly, we calculate Betti functions for two of the MassiveNuS simulations and discuss implications for using the persistent homology of the density field to help break degeneracy in the cosmological parameters.

Accurately predicting the abundance and structural evolution of dark matter subhaloes is crucial for understanding galaxy formation, modelling galaxy clustering, and constraining the nature of dark matter. Due to the non-linear nature of subhalo evolution, cosmological N-body simulations remain its primary method of investigation. However, it has recently been demonstrated that such simulations are still heavily impacted by artificial disruption, diminishing the information content on small scales, and reducing the reliability of all simulation-calibrated semi-analytical models. In this paper, we utilize the recently released DASH library of high-resolution, idealized simulations of the tidal evolution of subhaloes, which are unhindered by numerical overmerging due to discreteness noise or force softening, to calibrate an improved, more accurate model of the evolution of the density profiles of subhaloes that undergo tidal heating and stripping within their host halo. By testing previous findings that the structural evolution of a tidally truncated subhalo depends solely on the fraction of mass stripped, independent of the details of the stripping, we identify an additional dependence on the initial subhalo concentration. We provide significantly improved fitting functions for the subhalo density profiles and structural parameters (Vmax and rmax) that are unimpeded by numerical systematics and applicable to a wide range of parameter space. This model will be an integral component of a future semi-analytical treatment of substructure evolution, which can be used to predict key quantities, such as the evolved subhalo mass function and annihilation boost factors, and validate such calculations performed with cosmological simulations.

We present a machine-learning approach for estimating galaxy cluster masses, trained using both Chandra and eROSITA mock X-ray observations of 2041 clusters from the Magneticum simulations. We train a random forest (RF) regressor, an ensemble learning method based on decision tree regression, to predict cluster masses using an input feature set. The feature set uses core-excised X-ray luminosity and a variety of morphological parameters, including surface brightness concentration, smoothness, asymmetry, power ratios, and ellipticity. The regressor is cross-validated and calibrated on a training sample of 1615 clusters (80% of sample), and then results are reported as applied to a test sample of 426 clusters (20% of sample). This procedure is performed for two different mock observation series in an effort to bracket the potential enhancement in mass predictions that can be made possible by including dynamical state information. The first series is computed from idealized Chandra-like mock cluster observations, with high spatial resolution, long exposure time (1 Ms), and the absence of background. The second series is computed from realistic-condition eROSITA mocks with lower spatial resolution, short exposures (2 ks), instrument effects, and background photons modeled. We report a 20% reduction in the mass estimation scatter when either series is used in our RF model compared to a standard regression model that only employs core-excised luminosity. The morphological parameters that hold the highest feature importance are smoothness, asymmetry, and surface brightness concentration. Hence these parameters, which encode the dynamical state of the cluster, can be used to make more accurate predictions of cluster masses in upcoming surveys, offering a crucial step forward for cosmological analyses.

We investigate dynamical self-friction, the process by which material that is stripped from a subhalo torques its remaining bound remnant, which causes it to lose orbital angular momentum. By running idealized simulations of a subhalo orbiting within an analytical host halo potential, we isolate the effect of self-friction from traditional dynamical friction due to the host halo. While at some points in a subhalo’s orbit the torque of the stripped material can boost the orbital angular momentum of the remnant, the net effect over the long term is orbital decay regardless of the initial orbital parameters or subhalo mass. In order to quantify the strength of self-friction, we run a suite of simulations spanning typical host-to-subhalo mass ratios and orbital parameters. We find that the time-scale for self-friction, defined as the exponential decay time of the subhalo’s orbital angular momentum, scales with mass ratio and orbital circularity similar to standard dynamical friction. The decay time due to self-friction is roughly an order of magnitude longer, suggesting that self-friction only contributes at the 10 per cent level. However, along more radial orbits, self-friction can occasionally dominate over dynamical friction close to pericentric passage, where mass stripping is intense. This is also the epoch at which the self-friction torque undergoes large and rapid changes in both magnitude and direction, indicating that self-friction is an important process to consider when modelling pericentric passages of subhaloes and their associated satellite galaxies.

ABSTRACT X-ray and microwave cluster scaling relations are immensely valuable for cosmological analysis. However, their power is limited by astrophysical systematics that bias mass estimates and introduce additional scatter. Turbulence injected into the intracluster medium via mass assembly contributes substantially to cluster non-thermal pressure support, a significant source of such uncertainties. We use an analytical model to compute the assembly-driven non-thermal pressure profiles of haloes based on Monte Carlo-generated accretion histories. We introduce a fitting function for the average non-thermal pressure fraction profile, which exhibits minimal dependence on redshift at fixed peak height. Using the model, we predict deviations from self-similarity and the intrinsic scatter in the Sunyaev–Zel’dovich effect observable-mass scaling relation (YSZ − M) due solely to inter-cluster variation in mass accretion histories. We study the dependence of YSZ − M on aperture radius, cosmology, redshift, and mass limit. The model predicts $5-9{{\ \rm per\ cent}}$ scatter in YSZ − M at z = 0, increasing as the aperture used to compute YSZ increases from R500c to 5R500c. The predicted scatter lies slightly below that of studies based on non-radiative hydro-simulations, illustrating that assembly history variance is likely responsible for a substantial fraction of scatter in YSZ − M. This should be regarded as a lower bound, which will likely increase with the use of an updated gas density model that incorporates a more realistic response to halo assembly. As redshift increases, YSZ − M deviates more from self-similarity and scatter increases. We show that the YSZ − M residuals correlate strongly with the recent halo mass accretion rate, potentially providing an opportunity to infer the latter.