Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEnergy equipartition between stellar and dark matter particles in cosmological simulations results in spurious growth of galaxy sizes
144
0
0.0
(
0
)
Authors
Aaron D. Ludlow
Ask ChatGPT about the research
No Arabic abstract
The impact of 2-body scattering on the innermost density profiles of dark matter haloes is well established. We use a suite of cosmological simulations and idealised numerical experiments to show that 2-body scattering is exacerbated in situations where there are two species of unequal mass. This is a consequence of mass segregation and reflects a flow of kinetic energy from the more to less massive particles. This has important implications for the interpretation of galaxy sizes in cosmological hydrodynamic simulations, which nearly always model stars with less massive particles than are used for the dark matter. We compare idealised models as well as simulations from the EAGLE project that differ only in the mass resolution of the dark matter component, but keep sub-grid physics, baryonic mass resolution and gravitational force softening fixed. If the dark matter particle mass exceeds the mass of stellar particles, then galaxy sizes--quantified by their projected half-mass radii, ${rm R_{50}}$--increase systematically with time until ${rm R_{50}}$ exceeds a small fraction of the redshift-dependent mean inter-particle separation, $l$ (${rm R_{50}}geq 0.05times l$). Our conclusions should also apply to simulations that adopt different hydrodynamic solvers, subgrid physics or adaptive softening, but in that case may need quantitative revision. Any simulation employing a stellar-to-dark matter particle mass ratio greater than unity will escalate spurious energy transfer from dark matter to baryons on small scales.
rate research
Read More
We use idealized N-body simulations of equilibrium stellar disks embedded within course-grained dark matter haloes to study the effects of spurious collisional heating on disk structure and kinematics. Collisional heating drives a systematic increase in both the vertical and radial velocity dispersions of disk stars, and leads to an artificial increase in the thickness and size of disks; the effects are felt at all galacto-centric radii, and are not limited to the central regions of galaxies. We demonstrate that relaxation is driven primarily by the coarse-grained nature of simulated dark matter haloes, with bulges, stellar haloes and disk stars contributing little to disk heating. The integrated effects of collisional heating are determined primarily by the mass of dark matter particles (or equivalently by the number of particles at fixed halo mass), their local density and characteristic velocity, but are largely insensitive to the masses of stellar particles. This suggests that the effects of numerical relaxation on simulated galaxies can be reduced by increasing the mass resolution of the dark matter in cosmological simulations, with limited benefits from increasing the baryonic (or stellar) mass resolution. We provide a simple empirical model that accurately captures the effects of collisional heating on the vertical and radial velocity dispersions of disk stars, as well as on their scale heights. We use the model to assess the extent to which spurious collisional relaxation may have affected the structure of simulated galaxy disks. For example, we find that dark matter haloes resolved with fewer than $approx 10^6$ particles will collisionally heat stars near the stellar half-mass radius such that their vertical velocity dispersion increases by more than 10 per cent of the halos virial velocity in approximately one Hubble time.
We explain dark energy with equipartition theorem in string landscape.
Many phenomenologically successful cosmological galaxy formation simulations employ kinetic winds to model galactic outflows, a crucial ingredient in obtaining predictions that agree with various observations. Yet systematic studies of how variations in kinetic wind scalings might alter observable galaxy properties are rare. Here we employ GADGET-3 simulations to study how the baryon cycle, stellar mass function, and other galaxy and CGM predictions vary as a function of the assumed outflow speed $v_w$ and the scaling of the mass loading factor $eta$ with velocity dispersion $sigma$. We design our fiducial model to reproduce the measured wind properties at 25% of the virial radius from the Feedback In Realistic Environments (FIRE) simulations. We find that a strong dependence of $eta sim sigma^5$ in low mass haloes with $sigma < 106 mathrm{km s^{-1}}$ is required to match the faint end of the stellar mass functions at $z > 1$. The wind speed also has a major impact, with faster winds significantly reducing wind recycling and heating more halo gas. Both effects result in less stellar mass growth in massive haloes and impact high ionization absorption in halo gas. We cannot simultaneously match the stellar content at $z=2$ and $z=0$ within a single model, suggesting that an additional feedback source such as AGN might be required in massive galaxies at lower redshifts, but the amount needed depends strongly on assumptions regarding the outflow properties. We run a 50 $mathrm{Mpc/h}$, $2times576^3$ simulation with our fiducial parameters and show that it matches a range of star-forming galaxy properties at $zsim0-2$. In closing, the results from simulations of galaxy formation are much more sensitive to small changes in the feedback implementation than to the hydrodynamic technique.
Hydrodynamical simulations of galaxy formation and evolution attempt to fully model the physics that shapes galaxies. The agreement between the morphology of simulated and real galaxies, and the way the morphological types are distributed across galaxy scaling relations are important probes of our knowledge of galaxy formation physics. Here we propose an unsupervised deep learning approach to perform a stringent test of the fine morphological structure of galaxies coming from the Illustris and IllustrisTNG (TNG100 and TNG50) simulations against observations from a subsample of the Sloan Digital Sky Survey. Our framework is based on PixelCNN, an autoregressive model for image generation with an explicit likelihood. We adopt a strategy that combines the output of two PixelCNN networks in a metric that isolates the fine morphological details of galaxies from the sky background. We are able to emph{quantitatively} identify the improvements of IllustrisTNG, particularly in the high-resolution TNG50 run, over the original Illustris. However, we find that the fine details of galaxy structure are still different between observed and simulated galaxies. This difference is driven by small, more spheroidal, and quenched galaxies which are globally less accurate regardless of resolution and which have experienced little improvement between the three simulations explored. We speculate that this disagreement, that is less severe for quenched disky galaxies, may stem from a still too coarse numerical resolution, which struggles to properly capture the inner, dense regions of quenched spheroidal galaxies.
We present a method to flexibly and self-consistently determine individual galaxies star formation rates (SFRs) from their host haloes potential well depths, assembly histories, and redshifts. The method is constrained by galaxies observed stellar mass functions, SFRs (specific and cosmic), quenched fractions, UV luminosity functions, UV-SM relations, IRX-UV relations, auto- and cross-correlation functions (including quenched and star-forming subsamples), and quenching dependence on environment; each observable is reproduced over the full redshift range available, up to 0<z<10. Key findings include: galaxy assembly correlates strongly with halo assembly; quenching at z>1 correlates strongly with halo mass; quenched fractions at fixed halo mass decrease with increasing redshift; massive quenched galaxies reside in higher-mass haloes than star-forming galaxies at fixed galaxy mass; star-forming and quenched galaxies star formation histories at fixed mass differ most at z<0.5; satellites have large scatter in quenching timescales after infall, and have modestly higher quenched fractions than central galaxies; Planck cosmologies result in up to 0.3 dex lower stellar mass-halo mass ratios at early times; and, nonetheless, stellar mass-halo mass ratios rise at z>5. Also presented are revised stellar mass-halo mass relations for all, quenched, star-forming, central, and satellite galaxies; the dependence of star formation histories on halo mass, stellar mass, and galaxy SSFR; quenched fractions and quenching timescale distributions for satellites; and predictions for higher-redshift galaxy correlation functions and weak lensing surface densities. The public data release (DR1) includes the massively parallel (>10^5 cores) implementation (the UniverseMachine), the newly compiled and remeasured observational data, derived galaxy formation constraints, and mock catalogues including lightcones.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
