-
Merger Response of Halo Anisotropy Properties
Authors:
Kuan Wang,
Philip Mansfield,
Dhayaa Anbajagane,
Camille Avestruz
Abstract:
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 haloes. These properties reflect the dynamical state of haloes, which in turn depends on the mass assembly history. In this work, we study the evolution of anisotropy pro…
▽ More
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 haloes. These properties reflect the dynamical state of haloes, 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 behaviour 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 behaviour of halo properties. This study paves the way towards a physical understanding of scaling relations, particularly to how systematics in their scatter are connected to the mass assembly histories of haloes.
△ Less
Submitted 14 November, 2023;
originally announced November 2023.
-
Symfind: Addressing the Fragility of Subhalo Finders and Revealing the Durability of Subhalos
Authors:
Philip Mansfield,
Elise Darragh-Ford,
Yunchong Wang,
Ethan O. Nadler,
Risa H. Wechsler
Abstract:
A major question in $Λ$CDM is what this theory actually predicts for the properties of subhalo populations. Subhalos are difficult to simulate and to find within simulations, and this propagates into uncertainty in theoretical predictions for satellite galaxies. We present Symfind, a new particle-tracking-based subhalo finder, and demonstrate that it can track subhalos to orders-of-magnitude lower…
▽ More
A major question in $Λ$CDM is what this theory actually predicts for the properties of subhalo populations. Subhalos are difficult to simulate and to find within simulations, and this propagates into uncertainty in theoretical predictions for satellite galaxies. We present Symfind, a new particle-tracking-based subhalo finder, and demonstrate that it can track subhalos to orders-of-magnitude lower masses than commonly used halo-finding tools, with a focus on Rockstar and consistent-trees. These longer survival mean that at a fixed peak subhalo mass, we find $\approx 15\%{-}40\%$ more subhalos within the virial radius, $R_\textrm{vir}$, and $\approx 35\%-120\%$ more subhalos within $R_\textrm{vir}/4$ in the Symphony dark-matter-only simulation suite. More subhalos are found as resolution is increased. We perform extensive numerical testing. In agreement with idealized simulations, we show that the $v_{\rm max}$ of subhalos is only resolved at high resolutions ($n_\textrm{peak}\gtrsim3\times 10^4$), but that mass loss itself can be resolved at much more modest particle counts ($n_\textrm{peak}\gtrsim4\times 10^3$). We show that Rockstar converges to false solutions for the mass function, radial distribution, and disruption masses of subhalos. We argue that our new method can trace resolved subhalos until the point of typical galaxy disruption without invoking ``orphan'' modeling. We outline a concrete set of steps for determining whether other subhalo finders meet the same criteria. We publicly release Symfind catalogs and particle data for the Symphony simulation suite at \url{http://web.stanford.edu/group/gfc/symphony}.
△ Less
Submitted 21 August, 2023;
originally announced August 2023.
-
Haunted haloes: tracking the ghosts of subhaloes lost by halo finders
Authors:
Benedikt Diemer,
Peter Behroozi,
Philip Mansfield
Abstract:
Dark matter subhaloes are key for the predictions of simulations of structure formation, but their existence frequently ends prematurely due to two technical issues, namely numerical disruption in N-body simulations and halo finders failing to identify them. Here we focus on the second issue, using the phase-space friends-of-friends halo finder ROCKSTAR as a benchmark (though we expect our results…
▽ More
Dark matter subhaloes are key for the predictions of simulations of structure formation, but their existence frequently ends prematurely due to two technical issues, namely numerical disruption in N-body simulations and halo finders failing to identify them. Here we focus on the second issue, using the phase-space friends-of-friends halo finder ROCKSTAR as a benchmark (though we expect our results to translate to comparable codes). We confirm that the most prominent cause for losing track of subhaloes is tidal distortion rather than a low number of particles. As a solution, we present a flexible post-processing algorithm that tracks all subhalo particles over time, computes subhalo positions and masses based on those particles, and progressively removes stripped matter. If a subhalo is lost by the halo finder, this algorithm keeps tracking its so-called ghost until it has almost no particles left or has truly merged with its host. We apply this technique to a large suite of N-body simulations and restore lost subhaloes to the halo catalogues, which has a dramatic effect on key summary statistics of large-scale structure. Specifically, the subhalo mass function increases by about 50% and the halo correlation function increases by a factor of two at small scales. While these quantitative results are somewhat specific to our algorithm, they demonstrate that particle tracking is a promising way to reliably follow haloes and reduce the need for orphan models. Our algorithm and augmented halo catalogues are publicly available.
△ Less
Submitted 1 May, 2023;
originally announced May 2023.
-
MultiCAM: A multivariable framework for connecting the mass accretion history of haloes with their properties
Authors:
Ismael Mendoza,
Philip Mansfield,
Kuan Wang,
Camille Avestruz
Abstract:
Models that connect galaxy and halo properties often summarize a halo's mass accretion history (MAH) with a single value, and use this value as the basis for predictions. However, a single-value summary fails to capture the complexity of MAHs and information can be lost in the process. We present MultiCAM, a generalization of traditional abundance matching frameworks, which can simultaneously conn…
▽ More
Models that connect galaxy and halo properties often summarize a halo's mass accretion history (MAH) with a single value, and use this value as the basis for predictions. However, a single-value summary fails to capture the complexity of MAHs and information can be lost in the process. We present MultiCAM, a generalization of traditional abundance matching frameworks, which can simultaneously connect the full MAH of a halo with multiple halo and/or galaxy properties. As a first case study, we apply MultiCAM to the problem of connecting dark matter halo properties to their MAHs in the context of a dark matter-only simulation. While some halo properties, such as concentration, are more strongly correlated to the early-time mass growth of a halo, others, like the virial ratio, have stronger correlations with late-time mass growth. This highlights the necessity of considering the impact of the entire MAH on halo properties. For most of the halo properties we consider, we find that MultiCAM models that use the full MAH achieve higher accuracy than conditional abundance matching models which use a single epoch. We also demonstrate an extension of MultiCAM that captures the covariance between predicted halo properties. This extension provides a baseline model for applications where the covariance between predicted properties is important.
△ Less
Submitted 10 July, 2023; v1 submitted 2 February, 2023;
originally announced February 2023.
-
Symphony: Cosmological Zoom-in Simulation Suites over Four Decades of Host Halo Mass
Authors:
Ethan O. Nadler,
Philip Mansfield,
Yunchong Wang,
Xiaolong Du,
Susmita Adhikari,
Arka Banerjee,
Andrew Benson,
Elise Darragh-Ford,
Yao-Yuan Mao,
Sebastian Wagner-Carena,
Risa H. Wechsler,
Hao-Yi Wu
Abstract:
We present Symphony, a compilation of $262$ cosmological, cold-dark-matter-only zoom-in simulations spanning four decades of host halo mass, from $10^{11}$$-$$10^{15}~M_{\mathrm{\odot}}$. This compilation includes three existing simulation suites at the cluster and Milky Way$-$mass scales, and two new suites: $39$ Large Magellanic Cloud-mass ($10^{11}~M_{\mathrm{\odot}}$) and $49$ strong-lens-anal…
▽ More
We present Symphony, a compilation of $262$ cosmological, cold-dark-matter-only zoom-in simulations spanning four decades of host halo mass, from $10^{11}$$-$$10^{15}~M_{\mathrm{\odot}}$. This compilation includes three existing simulation suites at the cluster and Milky Way$-$mass scales, and two new suites: $39$ Large Magellanic Cloud-mass ($10^{11}~M_{\mathrm{\odot}}$) and $49$ strong-lens-analog ($10^{13}~M_{\mathrm{\odot}}$) group-mass hosts. Across the entire host halo mass range, the highest-resolution regions in these simulations are resolved with a dark matter particle mass of $\approx 3\times 10^{-7}$ times the host virial mass and a Plummer-equivalent gravitational softening length of $\approx 9\times 10^{-4}$ times the host virial radius, on average. We measure correlations between subhalo abundance and host concentration, formation time, and maximum subhalo mass, all of which peak at the Milky Way host halo mass scale. Subhalo abundances are $\approx 50\%$ higher in clusters than in lower-mass hosts at fixed sub-to-host halo mass ratios. Subhalo radial distributions are approximately self-similar as a function of host mass and are less concentrated than hosts' underlying dark matter distributions. We compare our results to the semianalytic model $\mathrm{\texttt{Galacticus}}$, which predicts subhalo mass functions with a higher normalization at the low-mass end and radial distributions that are slightly more concentrated than Symphony. We use $\mathrm{\texttt{UniverseMachine}}$ to model halo and subhalo star formation histories in Symphony, and we demonstrate that these predictions resolve the formation histories of the halos that host nearly all currently observable satellite galaxies in the universe. To promote open use of Symphony, data products are publicly available at http://web.stanford.edu/group/gfc/symphony.
△ Less
Submitted 16 March, 2023; v1 submitted 6 September, 2022;
originally announced September 2022.
-
Approximating Density Probability Distribution Functions Across Cosmologies
Authors:
Huanqing Chen,
Nickolay Y. Gnedin,
Philip Mansfield
Abstract:
Using a suite of self-similar cosmological simulations, we measure the probability distribution functions (PDFs) of real-space density, redshift-space density, and their geometric mean. We find that the real-space density PDF is well-described by a function of two parameters: $n_s$, the spectral slope, and $σ_L$, the linear rms density fluctuation. For redshift-space density and the geometric mean…
▽ More
Using a suite of self-similar cosmological simulations, we measure the probability distribution functions (PDFs) of real-space density, redshift-space density, and their geometric mean. We find that the real-space density PDF is well-described by a function of two parameters: $n_s$, the spectral slope, and $σ_L$, the linear rms density fluctuation. For redshift-space density and the geometric mean of real- and redshift-space densities, we introduce a third parameter, $s_L={\sqrt{\langle(dv^L_{\rm pec}/dr)^2\rangle}}/{H}$. We find that density PDFs for the LCDM cosmology is also well-parameterized by these three parameters. As a result, we are able to use a suite of self-similar cosmological simulations to approximate density PDFs for a range of cosmologies. We make the density PDFs publicly available and provide an analytical fitting formula for them.
△ Less
Submitted 13 September, 2021;
originally announced September 2021.
-
How Biased Are Halo Properties in Cosmological Simulations?
Authors:
Philip Mansfield,
Camille Avestruz
Abstract:
Cosmological N-body simulations have been a major tool of theorists for decades, yet many of the numerical issues that these simulations face are still unexplored. This paper measures numerical biases in these large, dark matter-only simulations that affect the properties of their dark matter haloes. We compare many simulation suites in order to provide several tools for simulators and analysts wh…
▽ More
Cosmological N-body simulations have been a major tool of theorists for decades, yet many of the numerical issues that these simulations face are still unexplored. This paper measures numerical biases in these large, dark matter-only simulations that affect the properties of their dark matter haloes. We compare many simulation suites in order to provide several tools for simulators and analysts which help mitigate these biases. We summarise our comparisons with practical `convergence limits' that can be applied to a wide range of halo properties, including halo properties which are traditionally overlooked by the testing literature. We also find that the halo properties predicted by different simulations can diverge from one another at unexpectedly high resolutions. We demonstrate that many halo properties depend strongly on force softening scale and that this dependence leads to much of the measured divergence between simulations. We offer an empirical model to estimate the impact of such effects on the rotation curves of a halo population. This model can serve as a template for future empirical models of the biases in other halo properties.
△ Less
Submitted 19 August, 2020;
originally announced August 2020.
-
Debiased Galaxy Cluster Pressure Profiles from X-ray Observations and Simulations
Authors:
Yizhou He,
Philip Mansfield,
Markus M. Rau,
Hy Trac,
Nicholas Battaglia
Abstract:
We present an updated model for the average cluster pressure profile, adjusted for hydrostatic mass bias by combining results from X-ray observations with cosmological simulations. Our model estimates this bias by fitting a power-law to the relation between the "true" halo mass and X-ray cluster mass in hydrodynamic simulations (IllustrisTNG, BAHAMAS, and MACSIS). As an example application, we con…
▽ More
We present an updated model for the average cluster pressure profile, adjusted for hydrostatic mass bias by combining results from X-ray observations with cosmological simulations. Our model estimates this bias by fitting a power-law to the relation between the "true" halo mass and X-ray cluster mass in hydrodynamic simulations (IllustrisTNG, BAHAMAS, and MACSIS). As an example application, we consider the REXCESS X-ray cluster sample and the Universal Pressure Profile (UPP) derived from scaled and stacked pressure profiles. We find adjusted masses, $M_\mathrm{500c},$ that are $\lesssim$15% higher and scaled pressures $P/P_\mathrm{500c}$ that have $\lesssim$35% lower normalization than previously inferred. Our Debiased Pressure Profile (DPP) is well-fit by a Generalized Navarro-Frenk-White (GNFW) function, with parameters $[P_0,c_{500},α,β,γ]=[5.048,1.217,1.192,5.490,0.433]$ and does not require a mass-dependent correction term. When the DPP is used to model the Sunyaev-Zel'dovich (SZ) effect, we find that the integrated Compton $Y-M$ relation has only minor deviations from self-similar scaling. The thermal SZ angular power spectrum is lower in amplitude by approximately 30%, assuming nominal cosmological parameters (e.g. $Ω_\text{m}=0.3$, $σ_8 = 0.8$), and is broadly consistent with recent Planck results without requiring additional bias corrections.
△ Less
Submitted 7 July, 2021; v1 submitted 10 August, 2020;
originally announced August 2020.
-
The Sheet of Giants: Unusual Properties of the Milky Way's Immediate Neighbourhood
Authors:
Maria K. Neuzil,
Philip Mansfield,
Andrey V. Kravtsov
Abstract:
We quantify the shape and overdensity of the galaxy distribution in the `Local Sheet' within a sphere of $R=8$ Mpc, and compare these properties with the expectations of the $Λ$CDM model. We measure ellipsoidal axis ratios of $c/a\approx0.16$ and $b/a\approx0.79$, indicating that the distribution of galaxies in the Local Volume can be approximated by a flattened oblate ellipsoid, consistent with t…
▽ More
We quantify the shape and overdensity of the galaxy distribution in the `Local Sheet' within a sphere of $R=8$ Mpc, and compare these properties with the expectations of the $Λ$CDM model. We measure ellipsoidal axis ratios of $c/a\approx0.16$ and $b/a\approx0.79$, indicating that the distribution of galaxies in the Local Volume can be approximated by a flattened oblate ellipsoid, consistent with the `sheet'-like configuration noted in previous studies. In contrast with previous estimates that the Local Sheet has a density close to average, we find that the number density of faint and bright galaxies in the Local Volume is $\approx1.7$ and $\approx5.2$ times denser, respectively, than the mean number density of galaxies of the same luminosity. Comparison with simulations shows that the number density contrasts of bright and faint galaxies within $8$ Mpc alone make the Local Volume a $\approx 2.5σ$ outlier in the $Λ$CDM cosmology. Our results indicate that the cosmic neighbourhood of the Milky Way may be unusual for galaxies of similar luminosity. The impact of the peculiar properties of our neighbourhood on the properties of the Milky Way and other nearby galaxies is not yet understood and warrants further study.
△ Less
Submitted 27 March, 2020; v1 submitted 9 December, 2019;
originally announced December 2019.
-
The Three Causes of Low-Mass Assembly Bias
Authors:
Philip Mansfield,
Andrey V. Kravtsov
Abstract:
We present a detailed analysis of the physical processes that cause halo assembly bias -- the dependence of halo clustering on proxies of halo formation time. We focus on the origin of assembly bias in the mass range corresponding to the hosts of typical galaxies and use halo concentration as our chief proxy of halo formation time. We also repeat our key analyses across a broad range of halo masse…
▽ More
We present a detailed analysis of the physical processes that cause halo assembly bias -- the dependence of halo clustering on proxies of halo formation time. We focus on the origin of assembly bias in the mass range corresponding to the hosts of typical galaxies and use halo concentration as our chief proxy of halo formation time. We also repeat our key analyses across a broad range of halo masses and for alternative formation time definitions. We show that splashback subhaloes are responsible for two thirds of the assembly bias signal, but do not account for the entire effect. After splashback subhaloes have been removed, we find that the remaining assembly bias signal is due to a relatively small fraction ($\lesssim 10\%$) of haloes in dense regions. We test a number of additional physical processes thought to contribute to assembly bias and demonstrate that the two key processes are the slowing of mass growth by large-scale tidal fields and by the high velocities of ambient matter in sheets and filaments. We also rule out several other proposed physical causes of halo assembly bias. Based on our results, we argue that there are three processes that contribute to assembly bias of low-mass halos: large-scale tidal fields, gravitational heating due to the collapse of large-scale structures, and splashback subhaloes located outside the virial radius.
△ Less
Submitted 31 January, 2019;
originally announced February 2019.
-
The splashback radius of halos from particle dynamics. II. Dependence on mass, accretion rate, redshift, and cosmology
Authors:
Benedikt Diemer,
Philip Mansfield,
Andrey V. Kravtsov,
Surhud More
Abstract:
The splashback radius $R_{\rm sp}$, the apocentric radius of particles on their first orbit after falling into a dark matter halo, has recently been suggested as a physically motivated halo boundary that separates accreting from orbiting material. Using the SPARTA code presented in Paper I, we analyze the orbits of billions of particles in cosmological simulations of structure formation and measur…
▽ More
The splashback radius $R_{\rm sp}$, the apocentric radius of particles on their first orbit after falling into a dark matter halo, has recently been suggested as a physically motivated halo boundary that separates accreting from orbiting material. Using the SPARTA code presented in Paper I, we analyze the orbits of billions of particles in cosmological simulations of structure formation and measure $R_{\rm sp}$ for a large sample of halos that span a mass range from dwarf galaxy to massive cluster halos, reach redshift 8, and include WMAP, Planck, and self-similar cosmologies. We analyze the dependence of $R_{\rm sp}/R_{\rm 200m}$ and $M_{\rm sp}/M_{\rm 200m}$ on the mass accretion rate $Γ$, halo mass, redshift, and cosmology. The scatter in these relations varies between 0.02 and 0.1 dex. While we confirm the known trend that $R_{\rm sp}/R_{\rm 200m}$ decreases with $Γ$, the relationships turn out to be more complex than previously thought, demonstrating that $R_{\rm sp}$ is an independent definition of the halo boundary that cannot trivially be reconstructed from spherical overdensity definitions. We present fitting functions for $R_{\rm sp}/R_{\rm 200m}$ and $M_{\rm sp}/M_{\rm 200m}$ as a function of accretion rate, peak height, and redshift, achieving an accuracy of 5% or better everywhere in the parameter space explored. We discuss the physical meaning of the distribution of particle apocenters and show that the previously proposed definition of $R_{\rm sp}$ as the radius of the steepest logarithmic density slope encloses roughly three-quarters of the apocenters. Finally, we conclude that no analytical model presented thus far can fully explain our results.
△ Less
Submitted 14 July, 2017; v1 submitted 28 March, 2017;
originally announced March 2017.
-
Splashback Shells of Cold Dark Matter Halos
Authors:
Philip Mansfield,
Andrey V. Kravtsov,
Benedikt Diemer
Abstract:
The density field in the outskirts of dark matter halos is discontinuous due to a caustic formed by matter at its first apocenter after infall. In this paper, we present an algorithm to identify the "splashback shell" formed by these apocenters in individual simulated halos using only a single snapshot of the density field. We implement this algorithm in the code SHELLFISH (SHELL Finding In Sphero…
▽ More
The density field in the outskirts of dark matter halos is discontinuous due to a caustic formed by matter at its first apocenter after infall. In this paper, we present an algorithm to identify the "splashback shell" formed by these apocenters in individual simulated halos using only a single snapshot of the density field. We implement this algorithm in the code SHELLFISH (SHELL Finding In Spheroidal Halos) and demonstrate that the code identifies splashback shells correctly and measures their properties with an accuracy of $<5\%$ for halos with more than 50,000 particles and mass accretion rates of $Γ_\textrm{DK14}>0.5$. Using SHELLFISH, we present the first estimates for several basic properties of individual splashback shells, such as radius, $R_\textrm{sp}$, mass, and overdensity, and provide fits to the distribution of these quantities as functions of $Γ_\textrm{DK14}$, $ν_\textrm{200m}$, and $z.$ We confirm previous findings that $R_\textrm{sp}$ decreases with increasing $Γ_\textrm{DK14}$, but we show that independent of accretion rate, it also decreases with increasing $ν_\textrm{200m}$. We also study the 3D structures of these shells and find that they generally have non-ellipsoidal oval shapes. We find that splashback radii estimated by SHELLFISH are $20\%-30\%$ larger than those estimated in previous studies from stacked density profiles at high accretion rates. We demonstrate that the latter are biased low due to the contribution of high-mass subhalos to these profiles and show that using the median instead of mean density in each radial bin mitigates the effect of substructure on density profiles and removes the bias.
△ Less
Submitted 8 May, 2017; v1 submitted 5 December, 2016;
originally announced December 2016.
-
SCORCH I: The Galaxy-Halo Connection in the First Billion Years
Authors:
Hy Trac,
Renyue Cen,
Philip Mansfield
Abstract:
SCORCH (Simulations and Constructions of the Reionization of Cosmic Hydrogen) is a new project to study the Epoch of Reionization (EoR). In this first paper, we probe the connection between observed high-redshift galaxies and simulated dark matter halos to better understand the abundance and evolution of the primary source of ionizing radiation. High-resolution N-body simulations are run to quanti…
▽ More
SCORCH (Simulations and Constructions of the Reionization of Cosmic Hydrogen) is a new project to study the Epoch of Reionization (EoR). In this first paper, we probe the connection between observed high-redshift galaxies and simulated dark matter halos to better understand the abundance and evolution of the primary source of ionizing radiation. High-resolution N-body simulations are run to quantify the abundance of dark matter halos as a function of mass $M$, accretion rate $\dot{M}$, and redshift $z$. A new fit for the halo mass function $dn/dM$ is $\approx 20\%$ more accurate at the high-mass end where bright galaxies are expected to reside. A novel approach is used to fit the halo accretion rate function $dn/d\dot{M}$ in terms of the halo mass function. Abundance matching against the observed galaxy luminosity function is used to estimate the luminosity-mass relation and the luminosity-accretion-rate relation. The inferred star formation efficiency is not monotonic with $M$ nor $\dot{M}$, but reaches a maximum value at a characteristic mass $\sim 2 \times 10^{11}\ M_\odot$ and a characteristic accretion rate $\sim 6 \times 10^2\ M_\odot/{\rm yr}$ at $z \approx 6$. We find a universal EoR luminosity-accretion-rate relation and construct a fiducial model for the galaxy luminosity function. The Schechter parameters evolve such that $φ_\star$ decreases, $M_\star$ is more positive (fainter), and $α$ is more negative (steeper) at higher redshifts. We forecast for the upcoming James Webb Space Telescope and show that with apparent magnitude limit $m_{\rm AB} \approx 31\ (32)$, it can observe $\gtrsim 11\ (24)$ unlensed galaxies per square degree per unit redshift at least down to $M_\star$ at $z \lesssim 13\ (14)$.
△ Less
Submitted 29 October, 2015; v1 submitted 9 July, 2015;
originally announced July 2015.