N-body simulations are essential for understanding the formation and evolution of structure in the Universe. However, the discrete nature of these simulations affects their accuracy when modelling collisionless systems. We introduce a new approach to
simulate the gravitational evolution of cold collisionless fluids by solving the Vlasov-Poisson equations in terms of adaptively refineable Lagrangian phase space elements. These geometrical elements are piecewise smooth maps between Lagrangian space and Eulerian phase space and approximate the continuum structure of the distribution function. They allow for dynamical adaptive splitting to accurately follow the evolution even in regions of very strong mixing. We discuss in detail various one-, two- and three-dimensional test problems to demonstrate the performance of our method. Its advantages compared to N-body algorithms are: i) explicit tracking of the fine-grained distribution function, ii) natural representation of caustics, iii) intrinsically smooth gravitational potential fields, thus iv) eliminating the need for any type of ad-hoc force softening. We show the potential of our method by simulating structure formation in a warm dark matter scenario. We discuss how spurious collisionality and large-scale discreteness noise of N-body methods are both strongly suppressed, which eliminates the artificial fragmentation of filaments. Therefore, we argue that our new approach improves on the N-body method when simulating self-gravitating cold and collisionless fluids, and is the first method that allows to explicitly follow the fine-grained evolution in six-dimensional phase space.
I review the nature of three-dimensional collapse in the Zeldovich approximation, how it relates to the underlying nature of the three-dimensional Lagrangian manifold and naturally gives rise to a hierarchical structure formation scenario that progre
sses through collapse from voids to pancakes, filaments and then halos. I then discuss how variations of the Zeldovich approximation (based on the gravitational or the velocity potential) have been used to define classifications of the cosmic large-scale structure into dynamically distinct parts. Finally, I turn to recent efforts to devise new approaches relying on tessellations of the Lagrangian manifold to follow the fine-grained dynamics of the dark matter fluid into the highly non-linear regime and both extract the maximum amount of information from existing simulations as well as devise new simulation techniques for cold collisionless dynamics.
Understanding the velocity field is very important for modern cosmology: it gives insights to structure formation in general, and also its properties are crucial ingredients in modelling redshift-space distortions and in interpreting measurements of
the kinetic Sunyaev-Zeldovich effect. Unfortunately, characterising the velocity field in cosmological N-body simulations is inherently complicated by two facts: i) The velocity field becomes manifestly multi-valued after shell-crossing and has discontinuities at caustics. This is due to the collisionless nature of dark matter. ii) N-body simulations sample the velocity field only at a set of discrete locations, with poor resolution in low-density regions. In this paper, we discuss how the associated problems can be circumvented by using a phase-space interpolation technique. This method provides extremely accurate estimates of the cosmic velocity fields and its derivatives, which can be properly defined without the need of the arbitrary coarse-graining procedure commonly used. We explore in detail the configuration-space properties of the cosmic velocity field on very large scales and in the highly nonlinear regime. In particular, we characterise the divergence and curl of the velocity field, present their one-point statistics, analyse the Fourier-space properties and provide fitting formulae for the velocity divergence bias relative to the non-linear matter power spectrum. We furthermore contrast some of the interesting differences in the velocity fields of warm and cold dark matter models. We anticipate that the high-precision measurements carried out here will help to understand in detail the dynamics of dark matter and the structures it forms.
We investigate the problem of predicting the halo mass function from the properties of the Lagrangian density field. We focus on a perturbation spectrum with a small-scale cut-off (as in warm dark matter cosmologies). This cut-off results in a strong
suppression of low mass objects, providing additional leverage to rigorously test which perturbations collapse and to what mass. We find that all haloes are consistent with forming near peaks of the initial density field, with a strong correlation between proto-halo density and ellipticity. We demonstrate that, while standard excursion set theory with correlated steps completely fails to reproduce the mass function, the inclusion of the peaks constraint leads to the correct number of haloes but significantly underpredicts the masses of low-mass objects (with the predicted halo mass function at low masses behaving like dn/dln m ~ m^{2/3}). This prediction is very robust and cannot be easily altered within the framework of a single collapse barrier. The nature of collapse in the presence of a small-scale cut-off thus reveals that excursion set calculations require a more detailed understanding of the collapse-time of a general ellipsoidal perturbation to predict the ultimate collapsed mass of a peak -- a problem that has been hidden in the large abundance of small-scale structure in CDM. We demonstrate how this problem can be resolved within the excursion set framework.
Recently, we have shown how current cosmological N-body codes already follow the fine grained phase-space information of the dark matter fluid. Using a tetrahedral tesselation of the three-dimensional manifold that describes perfectly cold fluids in
six-dimensional phase space, the phase-space distribution function can be followed throughout the simulation. This allows one to project the distribution function into configuration space to obtain highly accurate densities, velocities, and velocity dispersions. Here, we exploit this technique to show first steps on how to devise an improved particle-mesh technique. At its heart, the new method thus relies on a piecewise linear approximation of the phase space distribution function rather than the usual particle discretisation. We use pseudo-particles that approximate the masses of the tetrahedral cells up to quadrupolar order as the locations for cloud-in-cell (CIC) deposit instead of the particle locations themselves as in standard CIC deposit. We demonstrate that this modification already gives much improved stability and more accurate dynamics of the collisionless dark matter fluid at high force and low mass resolution. We demonstrate the validity and advantages of this method with various test problems as well as hot/warm-dark matter simulations which have been known to exhibit artificial fragmentation. This completely unphysical behaviour is much reduced in the new approach. The current limitations of our approach are discussed in detail and future improvements are outlined.
The primordial velocity dispersion of dark matter is small compared to the velocities attained during structure formation. The initial density distribution is close to uniform and it occupies an initial sheet in phase space that is single valued in v
elocity space. Because of gravitational forces this three dimensional manifold evolves in phase space without ever tearing, conserving phase-space volume and preserving the connectivity of nearby points. N-body simulations already follow the motion of this sheet in phase space. This fact can be used to extract full fine-grained phase-space-structure information from existing cosmological N-body simulations. Particles are considered as the vertices of an unstructured three dimensional mesh, moving in six dimensional phase-space. On this mesh, mass density and momentum are uniquely defined. We show how to obtain the space density of the fluid, detect caustics, and count the number of streams as well as their individual contributions to any point in configuration-space. We calculate the bulk velocity, local velocity dispersions, and densities from the sheet - all without averaging over control volumes. This gives a wealth of new information about dark matter fluid flow which had previously been thought of as inaccessible to N-body simulations. We outline how this mapping may be used to create new accurate collisionless fluid simulation codes that may be able to overcome the sparse sampling and unphysical two-body effects that plague current N-body techniques.
We use a 380 h-1 pc resolution hydrodynamic AMR simulation of a cosmic filament to investigate the orientations of a sample of ~100 well-resolved galactic disks spanning two orders of magnitude in both stellar and halo mass. We find: (i) At z=0, ther
e is an almost perfect alignment at a median angle of 18 deg, in the inner dark matter halo regions where the disks reside, between the spin vector of the gaseous and stellar galactic disks and that of their inner host haloes. The alignment between galaxy spin and spin of the entire host halo is however significantly weaker, ranging from a median of ~46 deg at z=1 to ~50 deg at z=0. (ii) The most massive galaxy disks have spins preferentially aligned so as to point along their host filaments. (iii) The spin of disks in lower-mass haloes shows, at redshifts above z~0.5 and in regions of low environmental density, a clear signature of alignment with the intermediate principal axis of the large-scale tidal field. This behavior is consistent with predictions of linear tidal torque theory. This alignment decreases with increasing environmental density, and vanishes in the highest density regions. Non-linear effects in the high density environments are plausibly responsible for establishing this density-alignment correlation. We expect that our numerical results provide important insights for both understanding intrinsic alignment in weak lensing from the astrophysical perspective and formation and evolution processes of galactic disks in a cosmological context.
New statistical properties of dark matter halos in Lagrangian space are presented. Tracing back the dark matter particles constituting bound halos resolved in a series of N-body simulations, we measure quantitatively the correlations of the proto-hal
os inertia tensors with the local tidal tensors and investigate how the correlation strength depends on the proto-halos sphericity, local density and filtering scale. It is shown that the majority of the proto-halos exhibit strong correlations between the two tensors provided that the tidal field is smoothed on the proto-halos mass scale. The correlation strength is found to increase as the proto-halos sphericity increases, as the proto-halos mass increases, and as the local density becomes close to the critical value, delta_{ec}. It is also found that those peculiar proto-halos which exhibit exceptionally weak correlations between the two tensors tend to acquire higher specific angular momentum in Eulerian space, which is consistent with the linear tidal torque theory. In the light of our results, it is intriguing to speculate a hypothesis that the low surface brightness galaxies observed at present epoch correspond to the peculiar proto-halos with extreme low-sphericity whose inertia tensors are weakly correlated with the local tidal tensors.
Galaxies on the largest scales of the Universe are observed to be embedded in the filamentary cosmic web which is shaped by the nonlinear tidal field. As an efficient tool to quantitatively describe the statistics of this cosmic web, we present the a
nisotropic two-point correlation functions of the nonlinear traceless tidal field in the principal-axis frame, which are measured using numerical data from an N-body simulation. We show that both of the nonlinear density and traceless tidal fields are more strongly correlated along the directions perpendicular to the eigenvectors associated with the largest eigenvalues of the local tidal field. The correlation length scale of the traceless tidal field is found to be ~20 Mpc/h, which is much larger than that of the density field ~5 Mpc/h. We also provide analytic fitting formulae for the anisotropic correlation functions of the traceless tidal field, which turn out to be in excellent agreement with the numerical results. We expect that our numerical results and analytic formula are useful to disentangle cosmological information from the filamentary network of the large-scale structures.
We explore a possible origin for the puzzling anti-correlation between the formation epoch of galactic dark-matter haloes and their environment density. This correlation has been revealed from cosmological N-body simulations and is in conflict with t
he Extended Press-Schechter model of halo clustering. Using similar simulations, we first quantify the straightforward association of an early formation epoch with a reduced mass growth rate at late times. We then find that a primary driver of suppressed growth, by accretion and mergers, is tidal effects dominated by a neighbouring massive halo. The tidal effects range from a slowdown of the assembly of haloes due to the shear along the large-scale filaments that feed the massive halo to actual mass loss in haloes that pass through the massive halo. Using the restricted three-body problem, we show that haloes are prone to tidal mass loss within 1.5 virial radii of a larger halo. Our results suggest that the dependence of formation epoch on environment density is a secondary effect induced by the enhanced density of haloes in filaments near massive haloes where the tides are strong. Our measures of assembly rate are particularly correlated with the tidal field at high redshifts z~1.
