No Arabic abstract
Lagrangian algorithms to simulate the evolution of cold dark matter (CDM) are invaluable tools to generate large suites of mock halo catalogues. In this paper, we first show that the main limitation of current semi-analytical schemes to simulate the displacement of CDM is their inability to model the evolution of overdensities in the initial density field, a limit that can be circumvented by detecting halo particles in the initial conditions. We thus propose `MUltiscale Spherical Collapse Lagrangian Evolution Using Press-Schechter (muscle-ups), a new scheme that reproduces the results from Lagrangian perturbation theory on large scales, while improving the modelling of overdensities on small scales. In muscle-ups, we adapt the extended Press and Schechter (EPS) formalism to Lagrangian algorithms of the displacement field. For regions exceeding a collapse threshold in the density smoothed at a radius $R$, we consider all particles within a radius $R$ collapsed. Exploiting a multi-scale smoothing of the initial density, we build a halo catalogue on the fly by optimizing the selection of halo candidates. This allows us to generate a density field with a halo mass function that matches one measured in $N$-body simulations. We further explicitly gather particles in each halo together in a profile, providing a numerical, Lagrangian-based implementation of the halo model. Compared to previous semi-analytical Lagrangian methods, we find that muscle-ups improves the recovery of the statistics of the density field at the level of the probability density function (PDF), the power spectrum, and the cross correlation with the $N$-body result.
We apply the model relating halo concentration to formation history proposed by Ludlow et al. to merger trees generated using an algorithm based on excursion set theory. We find that while the model correctly predicts the median relation between halo concentration and mass, it underpredicts the scatter in concentration at fixed mass. Since the same model applied to N-body merger trees predicts the correct scatter, we postulate that the missing scatter is due to the lack of any environmental dependence in merger trees derived from excursion set theory. We show that a simple modification to the merger tree construction algorithm, which makes merger rates dependent on environment, can increase the scatter by the required amount, and simultaneously provide a qualitatively correct correlation between environment and formation epoch in the excursion set merger trees.
Formalism of extended Lagrangian represent a systematic procedure to look for the local symmetries of a given Lagrangian action. In this work, the formalism is discussed and applied to a field theory. We describe it in detail for a field theory with first-class constraints present in the Hamiltonian formulation. The method is illustrated on examples of electrodynamics, Yang-Mills field and non-linear sigma model.
We present a modification of the Press-Schechter (PS) formalism to derive general mass functions for primordial black holes (PBHs), considering their formation as being associated to the amplitude of linear energy density fluctuations. To accommodate a wide range of physical relations between the linear and non-linear conditions for collapse, we introduce an additional parameter to the PS mechanism, and that the collapse occurs at either a given cosmic time, or as fluctuations enter the horizon. We study the case where fluctuations obey Gaussian statistics and follow a primordial power spectrum of broken power-law form with a blue spectral index for small scales. We use the observed abundance of super-massive black holes (SMBH) to constrain the extended mass functions taking into account dynamical friction. We further constrain the modified PS by developing a method for converting existing constraints on the PBH mass fraction, derived assuming monochromatic mass distributions for PBHs, into constraints applicable for extended PBH mass functions. We find that when considering well established monochromatic constraints there are regions in parameter space where all the dark matter can be made of PBHs. Of special interest is the region for the characteristic mass of the distribution ~10^2 M_Sun, for a wide range of blue spectral indices in the scenario where PBHs form as they enter the horizon, where the linear threshold for collapse is of the order of the typical overdensities, as this is close to the black hole masses detected by LIGO which are difficult to explain by stellar collapse.
We examine the power spectrum of clusters in the Press-Schechter (PS) theory and in N-body simulations to see how the power spectrum of clusters is related to the power spectrum of matter density fluctuations in the Universe. An analytic model for the power spectrum of clusters for their given number density is presented, both for real space and redshift space. We test this model against results from N-body simulations and find that the agreement between the analytic theory and the numerical results is good for wavelengths $lambda >60h^{-1}$ Mpc. On smaller scales non-linear processes that are not considered in the linear PS approximation influence the result. We also use our analytic model to study the redshift-space power spectrum of clusters in cold dark matter models with a cosmological constant ($Lambda$CDM) and with a scale-invariant Harrison-Zeldovich initial spectrum of density fluctuations. We find that power spectra of clusters in these models are not consistent with the observed power spectra of the APM and Abell-ACO clusters. One possible explanation for the observed power spectra of clusters is an inflationary scenario with a scalar field with the potential that has a localized steplike feature. We use the PS theory to examine the power spectrum of clusters in this model.
We consider the growth of primordial dark matter halos seeded by three crossed initial sine waves of various amplitudes. Using a Lagrangian treatment of cosmological gravitational dynamics, we examine the convergence properties of a high-order perturbative expansion in the vicinity of shell-crossing, by comparing the analytical results with state-of-the-art high resolution Vlasov-Poisson simulations. Based on a quantitative exploration of parameter space, we study explicitly for the first time the convergence speed of the perturbative series, and find, in agreement with intuition, that it slows down when going from quasi one-dimensional initial conditions (one sine wave dominating) to quasi triaxial symmetry (three sine waves with same amplitude). In most cases, the system structure at collapse time is, as expected, very similar to what is obtained with simple one-dimensional dynamics, except in the quasi-triaxial regime, where the phase-space sheet presents a velocity spike. In all cases, the perturbative series exhibits a generic convergence behavior as fast as an exponential of a power-law of the order of the expansion, allowing one to numerically extrapolate it to infinite order. The results of such an extrapolation agree remarkably well with the simulations, even at shell-crossing.