No Arabic abstract
We investigate the correlation between nine different dark matter halo properties using a rank correlation analysis and a Principal Component Analysis for a sample of haloes spanning five orders of magnitude in mass. We consider mass and dimensionless measures of concentration, age, relaxedness, sphericity, triaxiality, substructure, spin, and environment, where the latter is defined in a way that makes it insensitive to mass. We find that concentration is the most fundamental property. Except for environment, all parameters are strongly correlated with concentration. Concentration, age, substructure, mass, sphericity and relaxedness can be considered a single family of parameters, albeit with substantial scatter. In contrast, spin, environment, and triaxiality are more independent, although spin does correlate strongly with substructure and both spin and triaxiality correlate substantially with concentration. Although mass sets the scale of a halo, all other properties are more sensitive to concentration.
Generating pre-initial conditions (or particle loads) is the very first step to set up a cosmological N-body simulation. In this work, we revisit the numerical convergence of pre-initial conditions on dark matter halo properties using a set of simulations which only differs in initial particle loads, i.e. grid, glass, and the newly introduced capacity constrained Voronoi tessellation (CCVT). We find that the median halo properties agree fairly well (i.e. within a convergence level of a few per cent) among simulations running from different initial loads. We also notice that for some individual haloes cross-matched among different simulations, the relative difference of their properties sometimes can be several tens of per cent. By looking at the evolution history of these poorly converged haloes, we find that they are usually merging haloes or haloes have experienced recent merger events, and their merging processes in different simulations are out-of-sync, making the convergence of halo properties become poor temporarily. We show that, comparing to the simulation starting with an anisotropic grid load, the simulation with an isotropic CCVT load converges slightly better to the simulation with a glass load, which is also isotropic. Among simulations with different pre-initial conditions, haloes in higher density environments tend to have their properties converged slightly better. Our results confirm that CCVT loads behave as well as the widely used grid and glass loads at small scales, and for the first time we quantify the convergence of two independent isotropic particle loads (i.e. glass and CCVT) on halo properties.
We present a study of unprecedented statistical power regarding the halo-to-halo variance of dark matter substructure. Using a combination of N-body simulations and a semi-analytical model, we investigate the variance in subhalo mass fractions and subhalo occupation numbers, with an emphasis on how these statistics scale with halo formation time. We demonstrate that the subhalo mass fraction, f_sub, is mainly a function of halo formation time, with earlier forming haloes having less substructure. At fixed formation redshift, the average f_sub is virtually independent of halo mass, and the mass dependence of f_sub is therefore mainly a manifestation of more massive haloes assembling later. We compare observational constraints on f_sub from gravitational lensing to our model predictions and simulation results. Although the inferred f_sub are substantially higher than the median LCDM predictions, they fall within the 95th percentile due to halo-to-halo variance. We show that while the halo occupation distribution of subhaloes, P(N|M), is super-Poissonian for large <N>, a well established result, it becomes sub-Poissonian for <N> < 2. Ignoring the non-Poissonity results in systematic errors of the clustering of galaxies of a few percent, and with a complicated scale- and luminosity-dependence. Earlier-formed haloes have P(N|M) closer to a Poisson distribution, suggesting that the dynamical evolution of subhaloes drives the statistics towards Poissonian. Contrary to a recent claim, the non-Poissonity of subhalo occupation statistics does not vanish by selecting haloes with fixed mass and fixed formation redshift. Finally, we use subhalo occupation statistics to put loose constraints on the mass and formation redshift of the Milky Way halo. Using observational constraints on the V_max of the most massive satellites, we infer that 0.25<M_vir/10^12M_sun/h<1.4 and 0.1<z_f<1.4 at 90% confidence.
Dissipative dark matter self-interactions can affect halo evolution and change its structure. We perform a series of controlled N-body simulations to study impacts of the dissipative interactions on halo properties. The interplay between gravitational contraction and collisional dissipation can significantly speed up the onset of gravothermal collapse, resulting in a steep inner density profile. For reasonable choices of model parameters controlling the dissipation, the collapse timescale can be a factor of 10-100 shorter than that predicted in purely elastic self-interacting dark matter. The effect is maximized when energy loss per collision is comparable to characteristic kinetic energy of dark matter particles in the halo. Our simulations provide guidance for testing the dissipative nature of dark matter with astrophysical observations.
Spectroscopic observations at the Russian 6-m telescope are used to study the two polar ring galaxies (PRGs) from the catalogue by Moiseev et al.: SPRC-7 and SPRC-260. We have analyzed the kinematics of the stellar component of the central galaxies as well as the ionized gas kinematics in the external ring structures. The disc-halo decomposition of rotation curves in two perpendicular directions are considered. The observed 2D velocity fields are compared with the model predictions for different dark halo shapes. Based on these data, we constrain that for potential of DM halo semiaxis ratios is $s=0.8$, $q=1$ for SPRC-7 and $s=0.95$, $q=1.1$ for SPRC-260. Using 3D hydrodynamic simulations we also study the dynamics and evolution of the polar component in the potential of the galactic disc and dark halo for these two galaxies. We show that the polar component is dynamically quasi-stable on the scale of $sim10$ dynamical times (about a few Gyr). This is demonstrate the possibility for the growth of a spiral structure, which then steadily transforms to a lopsided gaseous system in the polar pane.
Studies of flux anomalies statistics and perturbations in stellar streams have the potential to constrain models of warm dark matter (WDM), including sterile neutrinos. Producing these constraints requires a parametrization of the WDM mass function relative to that of the cold dark matter (CDM) equivalent. We use five WDM models with half-mode masses, $M_mathrm{hm}=[1.3,35]times10^{8}$~$M_{odot}$, spread across simulations of the Local Group, lensing ellipticals and the $z=2$ universe, to generate such a parametrization: we fit parameters to a functional form for the WDM-to-CDM halo mass function ratio, $n_mathrm{WDM}(M_{X})/n_mathrm{CDM}(M_{X})$, of ($1+(alpha M_mathrm{hm}/M_{X})^{beta})^{gamma}$. For $M_{X}equiv$ virial mass of central halos we obtain $alpha=2.3$, $beta=0.8$, and $gamma=-1.0$, and this fit is steeper than the extended Press-Schechter formalism predicts. For $M_{X}equiv$ mass of subhalos we instead obtain $alpha=4.2$, $beta=2.5$ and $gamma=-0.2$; in both mass definitions the scatter is $sim20$~per~cent. The second fit typically underestimates the relative abundance of $z=2$ WDM subhaloes at the tens of per cent level. We caution that robust constraints will require bespoke simulations and a careful definition of halo mass, particularly for subhalos of mass $<10^{8}M_{odot}$.