ABSTRACT
Smoothed particle hydrodynamics (SPH) is a Lagrangian method for solving the fluid equations that is commonplace in astrophysics, prized for its natural adaptivity and stability. The choice ...of variable to smooth in SPH has been the topic of contention, with smoothed pressure (P-SPH) being introduced to reduce errors at contact discontinuities relative to smoothed density schemes. Smoothed pressure schemes produce excellent results in isolated hydrodynamics tests; in more complex situations however, especially when coupling to the ‘sub-grid’ physics and multiple time-stepping used in many state-of-the-art astrophysics simulations, these schemes produce large force errors that can easily evade detection as they do not manifest as energy non-conservation. Here, two scenarios are evaluated: the injection of energy into the fluid (common for stellar feedback) and radiative cooling. In the former scenario, force and energy conservation errors manifest (of the same order as the injected energy), and in the latter large force errors that change rapidly over a few time-steps lead to instability in the fluid (of the same order as the energy lost to cooling). Potential ways to remedy these issues are explored with solutions generally leading to large increases in computational cost. Schemes using a density-based formulation do not create these instabilities and as such it is recommended that they are preferred over pressure-based solutions when combined with an energy diffusion term to reduce errors at contact discontinuities.
Abstract We examine the central stellar velocity dispersion of subhalos based on IllustrisTNG cosmological hydrodynamic simulations. The central velocity dispersion is a fundamental observable that ...links galaxies with their dark matter subhalos. We carefully explore simulated stellar velocity dispersions derived with different definitions to assess possible systematics. We explore the impact of variation in the identification of member stellar particles, the viewing axes, the velocity dispersion computation technique, and simulation resolution. None of these issues impact the velocity dispersion significantly; any systematic uncertainties are smaller than the random error. We examine the stellar mass–velocity dispersion relation as an observational test of the simulations. At fixed stellar mass, the observed velocity dispersions significantly exceed the simulation results. This discrepancy is an interesting benchmark for the IllustrisTNG simulations because the simulations are not explicitly tuned to match this relation. We demonstrate that the stellar velocity dispersion provides measures of the dark matter velocity dispersion and the dark matter subhalo mass.
Abstract
A fundamental requirement for reionizing the Universe is that a sufficient fraction of the ionizing photons emitted by galaxies successfully escapes into the intergalactic medium. However, ...due to the scarcity of high-redshift observational data, the sources driving reionization remain uncertain. In this work, we calculate the ionizing escape fractions (fesc) of reionization-era galaxies from the state-of-the-art thesan simulations, which combine an accurate radiation-hydrodynamic solver (arepo-rt) with the well-tested IllustrisTNG galaxy formation model to self-consistently simulate both small-scale galaxy physics and large-scale reionization throughout a large patch of the universe ($L_\text{box} = 95.5\, \text{cMpc}$). This allows the formation of numerous massive haloes ($M_\text{halo} \gtrsim 10^{10}\, {\text{M}_{\odot }}$), which are often statistically underrepresented in previous studies but are believed to be important to achieving rapid reionization. We find that low-mass galaxies ($M_\text{stars} \lesssim 10^7\, {\text{M}_{\odot }}$) are the main drivers of reionization above z ≳ 7, while high-mass galaxies ($M_\text{stars} \gtrsim 10^8\, {\text{M}_{\odot }}$) dominate the escaped ionizing photon budget at lower redshifts. We find a strong dependence of fesc on the effective star formation rate (SFR) surface density defined as the SFR per gas mass per escape area, i.e. $\bar{\Sigma }_\text{SFR} = \text{SFR}/M_\text{gas}/R_{200}^2$. The variation in halo escape fractions decreases for higher mass haloes, which can be understood from the more settled galactic structure, SFR stability, and fraction of sightlines within each halo significantly contributing to the escaped flux. Dust is capable of reducing the escape fractions of massive galaxies, but the impact on the global fesc depends on the dust model. Finally, active galactic nuclei are unimportant for reionization in thesan and their escape fractions are lower than stellar ones due to being located near the centres of galaxy gravitational potential wells.
ABSTRACT
Smoothed particle hydrodynamics (SPH) is a ubiquitous numerical method for solving the fluid equations, and is prized for its conservation properties, natural adaptivity, and simplicity. We ...introduce the Sphenix SPH scheme, which was designed with three key goals in mind: to work well with sub-grid physics modules that inject energy, be highly computationally efficient (both in terms of compute and memory), and to be Lagrangian. sphenix uses a Density-Energy equation of motion, along with a variable artificial viscosity and conduction, including limiters designed to work with common sub-grid models of galaxy formation. In particular, we present and test a novel limiter that prevents conduction across shocks, preventing spurious radiative losses in feedback events. Sphenix is shown to solve many difficult test problems for traditional SPH, including fluid mixing and vorticity conservation, and it is shown to produce convergent behaviour in all tests where this is appropriate. Crucially, we use the same parameters within sphenix for the various switches throughout, to demonstrate the performance of the scheme as it would be used in production simulations. sphenix is the new default scheme in the swift cosmological simulation code and is available open source.
ABSTRACT
Cloud–wind interactions are common in the interstellar and circumgalactic media. Many studies have used simulations of such interactions to investigate the effect of particular physical ...processes, but the impact of the choice of hydrodynamic solver has largely been overlooked. Here we study the cloud–wind interaction, also known as the ‘blob test’, using seven different hydrodynamic solvers: three flavours of SPH, a moving mesh, adaptive mesh refinement, and two meshless schemes. The evolution of masses in dense gas and intermediate-temperature gas, as well as the covering fraction of intermediate-temperature gas, are systematically compared for initial density contrasts of 10 and 100, and five numerical resolutions. To isolate the differences due to the hydrodynamic solvers, we use idealized non-radiative simulations without physical conduction. We find large differences between these methods. SPH methods show slower dispersal of the cloud, particularly for the higher density contrast, but faster convergence, especially for the lower density contrast. Predictions for the intermediate-temperature gas differ particularly strongly, also between non-SPH codes, and converge most slowly. We conclude that the hydrodynamical interaction between a dense cloud and a supersonic wind remains an unsolved problem. Studies aiming to understand the physics or observational signatures of cloud–wind interactions should test the robustness of their results by comparing different hydrodynamic solvers.
ABSTRACT
We explore how the splashback radius (Rsp) of galaxy clusters, measured using the number density of the subhalo population, changes based on various selection criteria using the IllustrisTNG ...cosmological galaxy formation simulation. We identify Rsp by extracting the steepest radial gradient in a stacked set of clusters in 0.5 dex wide mass bins, with our clusters having halo masses 1013 ≤ M200,mean/M⊙ ≤ 1015. We apply cuts in subhalo mass, galaxy stellar mass, i-band absolute magnitude, and specific star formation rate. We find that, generally, galaxies of increasing mass and luminosity trace smaller measured splashback radii relative to the intrinsic dark matter radius. We also show that quenched galaxies may be used to reliably reconstruct the dark matter splashback radius. This trend is likely due to changes in the galaxy population. Additionally, we are able to reconcile different observational predictions that Rsp based upon galaxy number counts and dark matter may either align or show significant offset (e.g. those using optically or SZ-selected clusters) through the selection functions that these studies employ. Finally, we demonstrate that changes in Rsp measured through number counts are not due to a simple change in galaxy abundance inside and outside of the cluster.
ABSTRACT
We present a framework for characterizing the large-scale movement of baryons relative to dark matter in cosmological simulations, requiring only the initial conditions and final state of ...the simulation. This is performed using the spread metric that quantifies the distance in the final conditions between initially neighbouring particles, and by analysing the baryonic content of final haloes relative to that of the initial Lagrangian regions (LRs) defined by their dark matter component. Applying this framework to the simba cosmological simulations, we show that 40 per cent (10 per cent) of cosmological baryons have moved $\gt 1\, h^{-1}\, {\rm Mpc}{}$ ($3\, h^{-1}\, {\rm Mpc}{}$) by z = 0, primarily due to entrainment of gas by jets powered by an active galactic nucleus, with baryons moving up to $12\, h^{-1}\, {\rm Mpc}{}$ away in extreme cases. Baryons decouple from the dynamics of the dark matter component due to hydrodynamic forces, radiative cooling, and feedback processes. As a result, only 60 per cent of the gas content in a given halo at z = 0 originates from its LR, roughly independent of halo mass. A typical halo in the mass range Mvir = 1012–1013 M⊙ only retains 20 per cent of the gas originally contained in its LR. We show that up to 20 per cent of the gas content in a typical Milky Way-mass halo may originate in the region defined by the dark matter of another halo. This inter-Lagrangian baryon transfer may have important implications for the origin of gas and metals in the circumgalactic medium of galaxies, as well as for semi-analytic models of galaxy formation and ‘zoom-in’ simulations.
ABSTRACT
Backsplash galaxies are galaxies that once resided inside a cluster, and have migrated back outside as they move towards the apocentre of their orbit. The kinematic properties of these ...galaxies are well understood, thanks to the significant study of backsplashers in dark matter-only simulations, but their intrinsic properties are not well-constrained due to modelling uncertainties in subgrid physics, ram pressure stripping, dynamical friction, and tidal forces. In this paper, we use the IllustrisTNG300-1 simulation, with a baryonic resolution of Mb ≈ 1.1 × 107 M⊙, to study backsplash galaxies around 1302 isolated galaxy clusters with mass 1013.0 < M200,mean/M⊙ < 1015.5. We employ a decision tree classifier to extract features of galaxies that make them likely to be backsplash galaxies, compared to nearby field galaxies, and find that backsplash galaxies have low gas fractions, high mass-to-light ratios, large stellar sizes, and low black hole occupation fractions. We investigate in detail the origins of these large sizes, and hypothesize their origins are linked to the tidal environments in the cluster. We show that the black hole recentring scheme employed in many cosmological simulations leads to the loss of black holes from galaxies accreted into clusters, and suggest improvements to these models. Generally, we find that backsplash galaxies are a useful population to test and understand numerical galaxy formation models due to their challenging environments and evolutionary pathways that interact with poorly constrained physics.
ABSTRACT We illuminate the altered evolution of galaxies in clusters compared to central galaxies by tracking galaxies in the IllustrisTNG300 simulation as they enter isolated clusters of mass 1013 < ...M200,mean/M⊙ < 1015 (at z = 0). We demonstrate significant trends in galaxy properties with residence time (time since first infall) and that there is a population of galaxies that remain star forming even many Gyr after their infall. By comparing the properties of galaxies at their infall time to their properties at z = 0, we show how scaling relations, like the stellar-to-halo mass ratio, shift as galaxies live in the cluster environment. Galaxies with a residence time of 10 Gyr increase their stellar-to-halo mass ratio, by around 1 dex. As measurements of the steepest slope of the galaxy cluster number density profile (Rst), frequently used as a proxy for the splashback radius, have been shown to depend strongly on galaxy selection, we show how Rst depends on galaxy residence time. Using galaxies with residence times less than one cluster crossing time (≈5 Gyr) to measure Rst leads to significant offsets relative to using the entire galaxy population. Galaxies must have had the opportunity to ‘splash back’ to the first caustic to trace out a representative value of Rst, potentially leading to issues for galaxy surveys using ultraviolet-selected galaxies. Our work demonstrates that the evolution of cluster galaxies continues well into their lifetime in the cluster and departs from a typical central galaxy evolutionary path.