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
Simulating the dust content of galaxies and their surrounding gas is challenging due to the wide range of physical processes affecting the dust evolution. Here we present cosmological ...hydrodynamical simulations of a cluster of galaxies, $M_\text{200,crit}=6 \times 10^{14}{\, \rm M_\odot }$, including a novel dust model for the moving mesh code arepo. This model includes dust production, growth, supernova-shock-driven destruction, ion-collision-driven thermal sputtering, and high-temperature dust cooling through far-infrared reradiation of collisionally deposited electron energies. Adopting a rather low thermal sputtering rate, we find, consistent with observations, a present-day overall dust-to-gas ratio of ∼2 × 10−5, a total dust mass of ${\sim } 2\times 10^9{\, \rm M_\odot }$, and a dust mass fraction of ∼3 × 10−6. The typical thermal sputtering time-scales within ${\sim } 100\, {\rm kpc}$ are around ${\sim } 10\, {\rm Myr}$, and increase towards the outer parts of the cluster to ${\sim } 10^3\, {\rm Myr}$ at a cluster-centric distance of $1\, {\rm Mpc}$. The condensation of gas-phase metals into dust grains reduces high-temperature metal-line cooling, but also leads to additional dust infrared cooling. The additional infrared cooling changes the overall cooling rate in the outer parts of the cluster, beyond ${\sim } 1\, {\rm Mpc}$, by factors of a few. This results in noticeable changes of the entropy, temperature, and density profiles of cluster gas once dust formation is included. The emitted dust infrared emission due to dust cooling is consistent with observational constraints.
ABSTRACT
We present a new set of cosmological zoom-in simulations of a Milky Way (MW)-like galaxy that for the first time include elastic velocity-dependent self-interacting dark matter (SIDM) and ...IllustrisTNG physics. With these simulations, we investigate the interaction between SIDM and baryons and its effects on the galaxy evolution process. We also introduce a novel set of modified dark matter-only simulations that can reasonably replicate the effects of fully realized hydrodynamics on the DM halo while simplifying the analysis and lowering the computational cost. We find that baryons change the thermal structure of the central region of the halo to a greater extent than the SIDM scatterings for MW-like galaxies. Additionally, we find that the new thermal structure of the MW-like halo causes SIDM to create cuspier central densities rather than cores because the SIDM scatterings remove the thermal support by transferring heat away from the centre of the galaxy. We find that this effect, caused by baryon contraction, begins to affect galaxies with a stellar mass of 108 M⊙ and increases in strength to the MW-mass scale.
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.
ABSTRACT Self-interacting dark matter (SIDM) is increasingly studied as a potential solution to small-scale discrepancies between simulations of cold dark matter (CDM) and observations. We examine a ...physically motivated two-state SIDM model with both elastic and inelastic scatterings. In particular, endothermic, exothermic, and elastic scattering have equal transfer cross-sections at high relative velocities ($v_{\rm rel}\gtrsim 400~{\rm km\, s}^{-1})$. In a suite of cosmological zoom-in simulation of Milky Way-size haloes, we vary the primordial state fractions to understand the impact of inelastic dark matter self-interactions on halo structure and evolution. In particular, we test how the initial conditions impact the present-day properties of dark matter haloes. Depending on the primordial state fraction, scattering reactions will be dominated by either exothermic or endothermic effects for high and low initial excited state fractions, respectively. We find that increasing the initial excited fraction reduces the mass of the main halo, as well as the number of subhaloes on all mass scales. The main haloes are cored, with lower inner densities and higher outer densities compared with CDM. Additionally, we find that the shape of the main halo becomes more spherical the higher the initial excited state fraction is. Finally, we show that the number of satellites steadily decreases with initial excited state fraction across all satellite masses.
ABSTRACT
We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an ...unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.
ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an ...unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.
Inferring warm dark matter masses with deep learning Rose, Jonah C; Torrey, Paul; Villaescusa-Navarro, Francisco ...
Monthly notices of the Royal Astronomical Society,
01/2024, Letnik:
527, Številka:
1
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
We present a new suite of over 1500 cosmological N-body simulations with varied warm dark matter (WDM) models ranging from 2.5 to 30 keV. We use these simulations to train Convolutional ...Neural Networks (CNNs) to infer WDM particle masses from images of DM field data. Our fiducial setup can make accurate predictions of the WDM particle mass up to 7.5 keV with an uncertainty of ±0.5 keV at a 95 per cent confidence level from (25 h−1Mpc)2 maps. We vary the image resolution, simulation resolution, redshift, and cosmology of our fiducial setup to better understand how our model is making predictions. Using these variations, we find that our models are most dependent on simulation resolution, minimally dependent on image resolution, not systematically dependent on redshift, and robust to varied cosmologies. We also find that an important feature to distinguish between WDM models is present with a linear size between 100 and 200 h−1 kpc. We compare our fiducial model to one trained on the power spectrum alone and find that our field-level model can make two times more precise predictions and can make accurate predictions to two times as massive WDM particle masses when used on the same data. Overall, we find that the field-level data can be used to accurately differentiate between WDM models and contain more information than is captured by the power spectrum. This technique can be extended to more complex DM models and opens up new opportunities to explore alternative DM models in a cosmological environment.
ABSTRACT
Self-interacting dark matter (SIDM) offers the potential to mitigate some of the discrepancies between simulated cold dark matter (CDM) and observed galactic properties. We introduce a ...physically motivated SIDM model to understand the effects of self interactions on the properties of Milky Way and dwarf galaxy sized haloes. This model consists of dark matter with a nearly degenerate excited state, which allows for both elastic and inelastic scattering. In particular, the model includes a significant probability for particles to up-scatter from the ground state to the excited state. We simulate a suite of zoom-in Milky Way-sized N-body haloes with six models with different scattering cross sections to study the effects of up-scattering in SIDM models. We find that the up-scattering reaction greatly increases the central densities of the main halo through the loss of kinetic energy. However, the physical model still results in significant coring due to the presence of elastic scattering and down-scattering. These effects are not as apparent in the subhalo population compared to the main halo, but the number of subhaloes is reduced compared to CDM.