No Arabic abstract
We construct merger trees from the largest database of dark matter haloes to date provided by the Millennium simulation to quantify the merger rates of haloes over a broad range of descendant halo mass (1e12 < M0 < 1e15 Msun), progenitor mass ratio (1e-3 < xi < 1), and redshift (0 < z < 6). We find the mean merger rate per halo, B/n, to have very simple dependence on M0, xi, and z, and propose a universal fitting form for B/n that is accurate to 10-20%. Overall, B/n depends very weakly on the halo mass (proportional to M0^0.08) and scales as a power law in the progenitor mass ratio (proportional to xi^-2) for minor mergers (xi < 0.1) with a mild upturn for major mergers. As a function of time, we find the merger rate per Gyr to evolve as (1+z)^n with n=2-2.3, while the rate per unit redshift is nearly independent of z. Several tests are performed to assess how our merger rates are affected by, e.g. the time interval between Millennium outputs, binary vs multiple progenitor mergers, and mass conservation and diffuse accretion during mergers. In particular, we find halo fragmentations to be a general issue in merger tree construction from N-body simulations and compare two methods for handling these events. We compare our results with predictions of two analytical models for halo mergers based on the Extended Press-Schechter (EPS) model and the coagulation theory. We find the EPS model to overpredict the major merger rates and underpredict the minor merger rates by up to a factor of a few.
We analyze the cold dark matter density profiles of 54 galaxy halos simulated with FIRE-2 galaxy formation physics, each resolved within $0.5%$ of the halo virial radius. These halos contain galaxies with masses that range from ultra-faint dwarfs ($M_star simeq 10^{4.5} M_{odot}$) to the largest spirals ($M_star simeq 10^{11} M_{odot}$) and have density profiles that are both cored and cuspy. We characterize our results using a new analytic density profile that extends the standard Einasto form to allow for a pronounced constant-density core in the resolved innermost radius. With one additional core-radius parameter, $r_{c}$, this core-Einasto profile is able to characterize the shape and normalization of our feedback-impacted dark matter halos. In order to enable comparisons with observations, we provide fitting functions for $r_{c}$ and other profile parameters as a function of both $M_star$ and $M_{star}/M_{rm halo}$. In agreement with similar studies done in the literature, we find that dark matter core formation is most efficient at the characteristic stellar-mass to halo-mass ratio $M_star/M_{rm halo} simeq 5 times 10^{-3}$, or $M_{star} sim 10^9 , M_{odot}$, with cores that are roughly the size of the galaxy half-light radius, $r_{c} simeq 1-5$ kpc. Furthermore, we find no evidence for core formation at radii $gtrsim 100 rm pc$ in galaxies with $M_{star}/M_{rm halo} < 5times 10^{-4}$ or $M_star lesssim 10^6 , M_{odot}$. For Milky Way-size galaxies, baryonic contraction often makes halos significantly more concentrated and dense at the stellar half-light radius than dark matter only runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of $simeq 0.5-2$ kpc in size. Recent evidence for a ${sim} 2$ kpc core in the Milky Ways dark matter halo is consistent with this expectation.
We investigate the transverse modes of the gravitational and velocity fields in $Lambda$CDM, based on a high-resolution simulation performed using the adaptive-mesh refinement general-relativistic $N$-body code GRAMSES. We study the generation of vorticity in the dark matter velocity field at low redshift, providing fits to the shape and evolution of its power spectrum over a range of scales. By analysing the gravitomagnetic vector potential, which is absent in Newtonian simulations, in dark matter haloes with masses ranging from $sim10^{12.5}~h^{-1}{M}_{odot}$ to $sim10^{15}~h^{-1}{M}_{odot}$, we find that its magnitude correlates with the halo mass, peaking in the inner regions. Nevertheless, on average, its ratio against the scalar gravitational potential remains fairly constant, below percent level, decreasing roughly linearly with redshift and showing a weak dependence on halo mass. Furthermore, we show that the gravitomagnetic acceleration in haloes peaks towards the core and reaches almost $10^{-10}$ $h$ cm/s$^2$ in the most massive halo of the simulation. However, regardless of the halo mass, the ratio between the magnitudes of the gravitomagnetic force and the standard gravitational force is typically at around the $10^{-5}$ level inside the haloes, again without significant radius dependence. This result confirms that the gravitomagnetic effects have a negligible impact on structure formation, even for the most massive structures, although its behaviour in low density regions remains to be explored. Likewise, the impact on observations remains to be understood in the future.
We analyze constraints on parameters characterizing the pre-inflating universe in an open inflation model with a present slightly open $Lambda$CDM universe. We employ an analytic model to show that for a broad class of inflation-generating effective potentials, the simple requirement that some fraction of the observed dipole moment represents a pre-inflation isocurvature fluctuation allows one to set upper and lower limits on the magnitude and wavelength scale of pre-inflation fluctuations in the inflaton field, and the curvature of the pre-inflation universe, as a function of the fraction of the total initial energy density in the inflaton field as inflation begins. We estimate that if the pre-inflation contribution to the current CMB dipole is near the upper limit set by the {it Planck} Collaboration then the current constraints on $Lambda$CDM cosmological parameters allow for the possibility of a significantly open $Omega_{i} le 0.4$ pre-inflating universe for a broad range of the fraction of the total energy in the inflaton field at the onset of inflation. This limit to $Omega_{i}$ is even smaller if a larger dark-flow tilt is allowed.
The homogeneous, isotropic, and flat $Lambda$CDM universe favored by observations of the cosmic microwave background can be described using only Euclidean geometry, locally correct Newtonian mechanics, and the basic postulates of special and general relativity. We present simple derivations of the most useful equations connecting astronomical observables (redshift, flux density, angular diameter, brightness, local space density,...) with the corresponding intrinsic properties of distant sources (lookback time, distance, spectral luminosity, linear size, specific intensity, source counts,...). We also present an analytic equation for lookback time that is accurate within 0.1% for all redshifts $z$. The exact equation for comoving distance is an elliptic integral that must be evaluated numerically, but we found a simple approximation with errors $< 0.2$% for all redshifts up to $z approx 50$.
We use the Evolution and Assembly of GaLaxies and their Environments ( EAGLE ) suite of hydrodynamical cosmological simulations to measure offsets between the centres of stellar and dark matter components of galaxies. We find that the vast majority (>95%) of the simulated galaxies display an offset smaller than the gravitational softening length of the simulations (Plummer-equivalent $epsilon = 700$ pc), both for field galaxies and satellites in clusters and groups. We also find no systematic trailing or leading of the dark matter along a galaxys direction of motion. The offsets are consistent with being randomly drawn from a Maxwellian distribution with $sigma leq 196$ pc. Since astrophysical effects produce no feasible analogues for the $1.62^{+0.47}_{-0.49}$ kpc offset recently observed in Abell 3827, the observational result is in tension with the collisionless cold dark matter model assumed in our simulations.