We address the problem of inferring the three-dimensional matter distribution from a sparse set of one-dimensional quasar absorption spectra of the Lyman-$\alpha$ forest. Using a Bayesian forward modelling approach, we focus on extending the dynamical model to a fully self-consistent hierarchical field-level prediction of redshift-space quasar absorption sightlines. Our field-level approach rests on a recently developed semiclassical analogue to Lagrangian perturbation theory (LPT), which improves over noise problems and interpolation requirements of LPT. It furthermore allows for a manifestly conservative mapping of the optical depth to redshift space. In addition, this new dynamical model naturally introduces a coarse-graining scale, which we exploit to accelerate the MCMC sampler using simulated annealing. By gradually reducing the effective temperature of the forward model, we can allow it to converge first on large spatial scales before the sampler becomes sensitive to the increasingly larger space of smaller scales. We demonstrate the advantages -- in terms of speed and noise properties -- of this field-level approach over using LPT as a forward model, and, using mock data, validate its performance to reconstruct three-dimensional primordial perturbations and matter distribution from sparse quasar sightlines.

At early times dark matter has a thermal velocity dispersion of unknown amplitude which, for warm dark matter models, can influence the formation of nonlinear structure on observable scales. We propose a new scheme to simulate cosmologies with a small-scale suppression of perturbations that combines two previous methods in a way that avoids the numerical artefacts which have so far prevented either from producing fully reliable results. At low densities and throughout most of the cosmological volume, we represent the dark matter phase-sheet directly using high-accuracy interpolation, thereby avoiding the artificial fragmentation which afflicts particle-based methods in this regime. Such phase-sheet methods are, however, unable to follow the rapidly increasing complexity of the denser regions of dark matter haloes, so for these we switch to an N-body scheme which uses the geodesic deviation equation to track phase-sheet properties local to each particle. In addition, we present a novel high-resolution force calculation scheme based on an oct-tree of cubic force resolution elements which is well suited to approximate the force-field of our combined sheet+particle distribution. Our hybrid simulation scheme enables the first reliable simulations of the internal structure of low-mass haloes in a warm dark matter cosmology.

In this paper, we present resolved observations of the Sunyaev-Zel'dovich (SZ) effect, obtained with the NIKA2 camera, towards the cluster of galaxies XLSSC102, a relatively low-mass system ($M_{500} \sim 2 \times 10^{14}$ M$_{\odot}$) at $z = 0.97$ detected from the XXL survey. We combine NIKA2 SZ data, XMM-Newton X-ray data, and Megacam optical data to explore, respectively, the spatial distribution of the gas electron pressure, the gas density, and the galaxies themselves. We find significant offsets between the X-ray peak, the SZ peak, the brightest cluster galaxy, and the peak of galaxy density. Additionally, the galaxy distribution and the gas present elongated morphologies. This is interpreted as the sign of a recent major merging event, which induced a local boost of the gas pressure towards the north of XLSSC102 and stripped the gas out of the galaxy group. The NIKA2 data are also combined with XXL data to construct the thermodynamic profiles of XLSSC102, obtaining relatively tight constraints up to about $\sim r_{500}$, and revealing properties that are typical of disturbed systems. We also explore the impact of the cluster centre definition and the implication of local pressure substructure on the recovered profiles. Finally, we derive the global properties of XLSSC102 and compare them to those of high-mass-and-low-redshift systems, finding no strong evidence for non-standard evolution. We also use scaling relations to obtain alternative mass estimates from our profiles. The variation between these different mass estimates reflects the difficulty to accurately measure the mass of low-mass clusters at z$\sim$1, especially with low signal-to-noise ratio (S/N) data and for a disturbed system. However, it also highlights the strength of resolved SZ observations alone and in combination with survey-like X-ray data.

Accurate cosmological simulations that include the effect of non-linear matter clustering as well as of massive neutrinos are essential for measuring the neutrino mass scale from upcoming galaxy surveys. Typically, Newtonian simulations are employed and the neutrino distribution is sampled with a large number of particles in order to beat down the shot noise. Here we perform very efficient simulations with light massive neutrinos which require virtually no extra cost over matter-only simulations, and furthermore do not require tracer particles for the neutrinos. Instead, we use a weak-field dictionary based on the recently developed Newtonion motion approach, where Newtonian simulations for matter are paired with a linear relativistic Boltzmann code to allow for an absorption of the neutrino evolution into a time-dependent coordinate transformation. For this, only minimal modifications in existing N-body codes are required, which we have explicitly implemented in $\textit{gevolution}$ and $\textit{gadget-2}$. Our fast method determines the non-linear matter power spectrum to permille-level precision when compared against state-of-the-art simulations that have been performed for $0.1\,{\rm eV} \leq \sum m_\nu \leq 0.3\,$eV.

We present anisotropic ``separate universe'' simulations which modify the N-body code Gadget-4 in order to represent a large-scale tidal field through an anisotropic expansion factor. These simulations are used to measure the linear, quasi-linear and nonlinear response of the matter power spectrum to a spatially uniform trace-free tidal field up to wavenumber $k = 10 h \text{Mpc}^{-1}$. Together with the response to a large-scale overdensity measured in previous work, this completely describes the nonlinear matter bispectrum in the squeezed limit. We find that the response amplitude does not approach zero on small scales in physical coordinates, but rather a constant value, $R_K\approx 0.5$ for $k \geq 3 h\text{Mpc}^{-1}$ at $z=0$. This shows that even the inner regions of haloes are affected by the large-scale tidal field. We also measure directly the alignment of halo shapes with the tidal field, finding a clear signal which increases with halo mass.

The coalescence of supermassive black holes (SMBHs) should generate the strongest sources of gravitational waves (GWs) in the Universe. However, the dynamics of their coalescence is the subject of much debate. In this study, we use a suite of $N$-body simulations to follow the merger of two nuclear star clusters (NSCs), each hosting a SMBH in their centre. We find that the presence of distinct star clusters around each SMBH has important consequences for the dynamical evolution of the SMBH binary: (i) The separation between the SMBHs decreases by a few orders of magnitude in the first few Myrs by the combined effects of dynamical friction and a drag force caused by tidally stripped stars. In fact, this is a significant speedup for equal mass ratio binaries, and becomes extreme for unequal mass ratios, e.g. 1:10 or 1:100, which traditional dynamical friction alone would not permit to bind. (ii) The subsequent binary hardening is driven by the gravitational slingshots between the SMBH binary and stars, and also depends on the mass ratio between the SMBHs. Thus, with this additional drag force, we find that all SMBHs in our suite coalesce within a Hubble time. Given that about 50% of Milky Way sized galaxies host NSCs, our results are encouraging for upcoming GW observations with the Laser Interferometer Space Antenna -- LISA -- which will detect SMBH coalescence in the $10^4-10^7$$M_{\rm \odot}$ mass range.

Gravitational collapse of cold dark matter leads to infinite-density caustics that seed the primordial dark-matter halos in the large-scale structure. The development of these caustics begins, generically, as an almost one-dimensional phenomenon with the formation of pancakes. Focusing on the one-dimensional case, we identify a landscape of so far unknown singularities in the particle acceleration that emerge after the first crossing of particle trajectories. We complement our fully analytical studies by high-resolution N-body simulations and find outstanding agreement, particularly shortly after the first crossing. We develop the methods in 1D but outline briefly the necessary steps for the 3D case.

The internal properties of dark matter haloes correlate with the large-scale halo clustering strength at fixed halo mass $-$ an effect known as assembly bias $-$ and are also strongly affected by the local, non-linear cosmic web. Characterising a halo's local web environment by its tidal anisotropy $\alpha$ at scales $\sim4$ x the halo radius, we demonstrate that these multi-scale correlations represent two distinct statistical links: one between the internal property and $\alpha$, and the other between $\alpha$ and large-scale ( $>30h^{-1}$Mpc) halo bias $b_1$. We focus on scalar internal properties of haloes related to formation time (concentration $c_{\rm vir}$), shape (mass ellipsoid asphericity $c/a$), velocity dispersion structure (velocity ellipsoid asphericity $c_v/a_v$ and velocity anisotropy $\beta$) and angular momentum (dimensionless spin $\lambda$) in the mass range $8\times10^{11}< M_{\rm vir}/(h^{-1}M_\odot)<5\times10^{14}$. Using conditional correlation coefficients and other detailed tests, we show that the joint distribution of $\alpha$, $b_1$ and any of the internal properties $c\in\{\beta,c_v/a_v,c/a,c_{\rm vir},\lambda\}$ is consistent with $p(\alpha,b_1,c)\simeq p(\alpha)p(b_1|\alpha)p(c|\alpha)$, at all but the largest masses. $\textit{Thus, the assembly bias trends $c-b_1$ reflect the two fundamental correlations $c-\alpha$ and $b_1-\alpha$.}$ Our results are unaffected by the exclusion of haloes with recent major merger events or splashback objects, although the latter are distinguished by the fact that $\alpha$ does not explain their assembly bias trends. The overarching importance of $\alpha$ provides a new perspective on the nature of assembly bias of distinct haloes, with potential ramifications for incorporating realistic assembly bias effects into mock catalogs of future large-scale structure surveys and for detecting galaxy assembly bias.

Gravitational collapse in cosmological context produces an intricate cosmic web of voids, walls, filaments and nodes. The anisotropic nature of collisionless collapse leads to the emergence of an anisotropic velocity dispersion, or stress, that absorbs most of the kinetic energy after shell-crossing. In this paper, we measure this large-scale velocity dispersion tensor $\sigma^2_{ij}$ in $N$-body simulations using the phase-space interpolation technique. We study the environmental dependence of the amplitude and anisotropy of the velocity dispersion tensor field, and measure its spatial correlation and alignment. The anisotropy of $\sigma^2_{ij}$ naturally encodes the collapse history and thus leads to a parameter-free identification of the four dynamically distinct cosmic web components. We find this purely dynamical classification to be in good agreement with some of the existing classification methods. In particular, we demonstrate that $\sigma^2_{ij}$ is well aligned with the large-scale tidal field. We further investigate the influence of small scale density fluctuations on the large scale velocity dispersion, and find that the measured amplitude and alignments are dominated by the largest perturbations and thus remain largely unaffected. We anticipate that these results will give important new insight into the anisotropic nature of gravitational collapse on large scales, and the emergence of anisotropic stress in the cosmic web.

We chart a path toward solving for the nonlinear gravitational dynamics of cold dark matter by relying on a semiclassical description using the propagator. The evolution of the propagator is given by a Schr\"odinger equation, where the small parameter $\hbar$ acts as a softening scale that regulates singularities at shell-crossing. The leading-order propagator, called free propagator, is the semiclassical equivalent of the Zel'dovich approximation (ZA), that describes inertial particle motion along straight trajectories. At next-to-leading order, we solve for the propagator perturbatively and obtain, in the classical limit the displacement field from second-order Lagrangian perturbation theory (LPT). The associated velocity naturally includes an additional term that would be considered as third order in LPT. We show that this term is actually needed to preserve the underlying Hamiltonian structure, and ignoring it could lead to the spurious excitation of vorticity in certain implementations of second-order LPT. We show that for sufficiently small $\hbar$ the corresponding propagator solutions closely resemble LPT, with the additions that spurious vorticity is avoided and the dynamics at shell-crossing is regularised. Our analytical results possess a symplectic structure that allows us to advance numerical schemes for the large-scale structure. For times shortly after shell-crossing, we explore the generation of vorticity, which in our method does not involve any explicit multi-stream averaging, but instead arises naturally as a conserved topological charge.

The abundance and demographics of dark matter substructure is important for many areas in astrophysics and cosmological $N$-body simulations have been the primary tool used to investigate them. However, it has recently become clear that the simulations are subject to numerical artefacts, which hampers a proper treatment of the tidal evolution of subhaloes. Unfortunately, no analytical models that accurately describe subhalo evolution exist either. We therefore present a library of idealized, high resolution $N$-body simulations of the tidal evolution of individual subhaloes that can be used to calibrate semi-analytical models and to complement cosmological simulations. The simulations focus on minor mergers, i.e., the mass of the subhalo is much smaller than that of the host halo, such that the impact of dynamical friction is negligible. This setup allows the adoption of a fixed analytical potential for modelling the host halo. The dynamical evolution of subhaloes is followed with $N$-body computations. In the library, four parameters, two of which characterize the subhalo orbit with respect to the host halo, and the two concentrations of the host- and subhalo, are varied over the ranges encountered in cosmological simulations. We show several representative examples from the library that illustrate the evolution of the subhalo mass and velocity dispersion profiles. Additionally, we make publicly available a pre-trained non-parametric model of the subhalo mass evolution based on random forest regression. This model is able to interpolate the simulation data at the 0.1\,dex level and provides efficient access to the data for further use in modelling.

We study the role of the local tidal environment in determining the assembly bias of dark matter haloes. Previous results suggest that the anisotropy of a halo's environment (i.e, whether it lies in a filament or in a more isotropic region) can play a significant role in determining the eventual mass and age of the halo. We statistically isolate this effect using correlations between the large-scale and small-scale environments of simulated haloes at $z=0$ with masses between $10^{11.6}\lesssim (m/h^{-1}M_{\odot})\lesssim10^{14.9}$. We probe the large-scale environment using a novel halo-by-halo estimator of linear bias. For the small-scale environment, we identify a variable $\alpha_R$ that captures the $\textit{tidal anisotropy}$ in a region of radius $R=4R_{\textrm{200b}}$ around the halo and correlates strongly with halo bias at fixed mass. Segregating haloes by $\alpha_R$ reveals two distinct populations. Haloes in highly isotropic local environments ($\alpha_R\lesssim0.2$) behave as expected from the simplest, spherically averaged analytical models of structure formation, showing a $\textit{negative}$ correlation between their concentration and large-scale bias at $\textit{all}$ masses. In contrast, haloes in anisotropic, filament-like environments ($\alpha_R\gtrsim0.5$) tend to show a $\textit{positive}$ correlation between bias and concentration at any mass. Our multi-scale analysis cleanly demonstrates how the overall assembly bias trend across halo mass emerges as an average over these different halo populations, and provides valuable insights towards building analytical models that correctly incorporate assembly bias. We also discuss potential implications for the nature and detectability of galaxy assembly bias.

While feedback from Active Galactic Nuclei (AGN) is an important heating source in the centre of galaxy clusters, it is still unclear how the feedback energy is injected into the intracluster medium (ICM) and what role different numerical approaches play. Here, we compare four hydrodynamical schemes in idealized simulations of a rising bubble inflated by AGN feedback in a hot stratified ICM: (traditional) smoothed particle hydrodynamics (TSPH), a pressure flavour of SPH (PSPH), a meshless finite mass (MFM) scheme, as well as an Eulerian code with adaptive mesh refinement. In the absence of magnetic fields, the bubble is Kelvin-Helmholtz unstable on short enough time scales to dissolve it fully in the ICM, which is captured by MFM and RAMSES simulations, while in the TSPH simulation the bubble survives. When the ICM is turbulent, mixing of the bubble with the ICM is accelerated. This occurs if the numerical scheme can capture the instabilities well. The differences in the evolution of the bubble has a surprisingly small influence on the thermal structure of the ICM. However, in the simulations with MFM and RAMSES the bubble disruption leads to turbulent stirring of the ICM which is suppressed in SPH. In the latter the thermal energy remains trapped in the bubble and is transported to large radii. We discuss if the choice of hydrodynamical schemes can lead to systematic differences in the outcomes of cosmological simulations.

The influence of the Cosmic Web on galaxy formation and evolution is of great observational and theoretical interest. We investigate whether the Cosmic Web leaves an imprint in the spatial clustering of galaxies in the Sloan Digital Sky Survey (SDSS), using the group catalog of Yang et al. and tidal field estimates at $\sim2h^{-1}$Mpc scales from the Mass-Tides-Velocity data set of Wang et al. We use the $\textit{tidal anisotropy}$ $\alpha$ (Paranjape et al.) to characterise the tidal environment of groups, and measure the redshift-space 2-point correlation function (2pcf) of group positions and the luminosity- and colour-dependent clustering of group galaxies using samples segregated by $\alpha$. We find that all the 2pcf measurements depend strongly on $\alpha$, with factors of $\sim20$ between the large-scale 2pcf of objects in the most and least isotropic environments. To test whether these strong trends imply `beyond halo mass' effects for galaxy evolution, we compare our results with corresponding 2pcf measurements in mock catalogs constructed using a halo occupation distribution that only uses halo mass as an input. We find that this prescription qualitatively reproduces $\textit{all}$ observed trends, and also quantitatively matches many of the observed results. Although there are some statistically significant differences between our `halo mass only' mocks and the data -- in the most and least isotropic environments -- which deserve further investigation, our results suggest that if the tidal environment induces additional effects on galaxy properties other than those inherited from their host halos, then these must be weak.

Sub-structures in the hot gas of galaxy clusters are related to their formation history and to the astrophysical processes at play in the intracluster medium (ICM). The thermal Sunyaev-Zel'dovich (tSZ) effect is directly sensitive to the line-of-sight integrated ICM pressure, and is thus particularly adapted to study ICM sub-structures. We apply structure-enhancement filtering algorithms to high resolution tSZ observations of distant clusters, in order to search for pressure discontinuities, compressions, as well as secondary peaks in the ICM. The same filters are applied to synthetic tSZ images extracted from RHAPSODY-G hydrodynamic simulations, in order to better interpret the extracted features. We also study the noise propagation trough the filters and quantify the impact of systematic effects, point source residuals being identified as the dominant potential contaminant. In 3 of our 6 NIKA clusters we identify features at high S/N that show clear evidence for merger events. In MACSJ0717 (z=0.55), three strong pressure gradients are observed on the E, SE and W sectors, and two main peaks in the pressure distribution are identified. We observe a lack of tSZ compact structure in the cool-core cluster PSZ1G045.85 (z=0.61), and a tSZ gradient ridge dominates in the SE. In the highest z cluster, CLJ1227 (z=0.89), we detect a ~45" (360 kpc) long ridge pressure gradient associated with a secondary pressure peak in the W region. Our results show that current tSZ facilities have now reached the angular resolution and sensitivity to allow an exploration of the details of pressure sub-structures in clusters, even at high z. This opens the possibility to quantify the impact of the dynamical state on the relation between the tSZ signal and the mass of clusters, which is important when using tSZ clusters to test cosmological models. This work also marks the first NIKA cluster sample data release.

Accurately predicting the demographics of dark matter (DM) substructure is of paramount importance for many fields of astrophysics, including gravitational lensing, galaxy evolution, halo occupation modeling, and constraining the nature of dark matter. Because of its strongly non-linear nature, DM substructure is typically modeled using N-body simulations, which reveal that large fractions of DM subhaloes undergo complete disruption. In this paper we use both analytical estimates and idealized numerical simulations to investigate whether this disruption is mainly physical, due to tidal heating and stripping, or numerical (i.e., artificial). We show that, contrary to naive expectation, subhaloes that experience a tidal shock $\Delta E$ that exceeds the subhalo's binding energy, $|E_{\rm b}|$, do not undergo disruption, even when $\Delta E/|E_{\rm b}|$ is as large as 100. Along the same line, and contrary to existing claims in the literature, instantaneously stripping matter from the outskirts of a DM subhalo also does not result in its complete disruption, even when the instantaneous remnant has positive binding energy. In addition, we show that tidal heating due to high-speed (impulsive) encounters with other subhaloes (`harassment'), is negligible compared to the tidal effects due to the host halo. Hence, we conclude that, in the absence of baryonic processes, the complete, physical disruption of CDM substructure is extremely rare, and that most disruption in numerical simulations therefore must be artificial. We discuss various processes that have been associated with numerical overmerging, and conclude that inadequate force-softening is the most likely culprit.

Dark matter (DM) haloes forming near the thermal cut-off scale of the density perturbations are unique, since they are the smallest objects and form through monolithic gravitational collapse, while larger haloes contrastingly have experienced mergers. While standard cold dark matter (CDM) simulations readily produce haloes that follow the universal Navarro-Frenk-White (NFW) density profile with an inner slope, $\rho \propto r^{-\alpha}$, with $\alpha=1$, recent simulations have found that when the free-streaming cut-off expected for the CDM model is resolved, the resulting haloes follow nearly power-law density profiles of $\alpha\sim1.5$. In this paper, we study the formation of density cusps in haloes using idealized $N$-body simulations of the collapse of proto-haloes. When the proto-halo profile is initially cored due to particle free-streaming at high redshift, we universally find $\sim r^{-1.5}$ profiles irrespective of the proto-halo profile slope outside the core and large-scale non-spherical perturbations. Quite in contrast, when the proto-halo has a power-law profile, then we obtain profiles compatible with the NFW shape when the density slope of the proto-halo patch is shallower than a critical value, $\alpha_{\rm ini} \sim 0.3$, while the final slope can be steeper for $\alpha_{\rm ini}\ga 0.3$. We further demonstrate that the $r^{-1.5}$ profiles are sensitive to small scale noise, which gradually drives them towards an inner slope of $-1$, where they become resilient to such perturbations. We demonstrate that the $r^{-1.5}$ solutions are in hydrostatic equilibrium, largely consistent with a simple analytic model, and provide arguments that angular momentum appears to determine the inner slope.

We simulate neutralino dark matter ($\chi$DM) haloes from their initial collapse, at $\sim$ earth mass, up to a few percent solar. Our results confirm that the density profiles of the first haloes are described by a $\sim r^{-1.5}$ power-law. As haloes grow in mass, their density profiles evolve significantly. In the central regions, they become shallower and reach on average $\sim r^{-1}$, the asymptotic form of an NFW profile. Using non-cosmological controlled simulations, we observe that temporal variations in the gravitational potential caused by major mergers lead to a shallowing of the inner profile. This transformation is more significant for shallower initial profiles and for a higher number of merging systems. Depending on the merger details, the resulting profiles can be shallower or steeper than NFW in their inner regions. Interestingly, mergers have a much weaker effect when the profile is given by a broken power-law with an inner slope of $-1$ (such as NFW or Hernquist profiles). This offers an explanation for the emergence of NFW-like profiles: after their initial collapse, $r^{-1.5}$ $\chi$DM haloes suffer copious major mergers, which progressively shallows the profile. Once an NFW-like profile is established, subsequent merging do not change the profile anymore. This suggests that halo profiles are not universal but rather a combination of (1) the physics of the formation of the microhaloes and (2) their early merger history -- both set by the properties of the dark matter particle -- as well as (3) the resilience of NFW-like profiles to perturbations.

The analytical formalism to obtain the probability distribution functions (PDFs) of spherically-averaged cosmic densities and velocity divergences in the mildly non-linear regime is presented. A large-deviation principle is applied to those cosmic fields assuming their most likely dynamics in spheres is set by the spherical collapse model. We validate our analytical results using state-of-the-art dark matter simulations with a phase-space resolved velocity field finding a 2% percent level agreement for a wide range of velocity divergences and densities in the mildly nonlinear regime (~10Mpc/h at redshift zero), usually inaccessible to perturbation theory. From the joint PDF of densities and velocity divergences measured in two concentric spheres, we extract with the same accuracy velocity profiles and conditional velocity PDF subject to a given over/under-density which are of interest to understand the non-linear evolution of velocity flows. Both PDFs are used to build a simple but accurate maximum likelihood estimators for the redshift evolution of the variance of both the density and velocity divergence fields, which have smaller relative errors than their sample variances when non-linearities appear. Given the dependence of the velocity divergence on the growth rate, there is a significant gain in using the full knowledge of both PDFs to derive constraints on the equation of state of dark energy. Thanks to the insensitivity of the velocity divergence to bias, its PDF can be used to obtain unbiased constraints on the growth of structures ($\sigma_8$,f) or it can be combined with the galaxy density PDF to extract bias parameters.

The cosmic web is one of the most striking features of the distribution of galaxies and dark matter on the largest scales in the Universe. It is composed of dense regions packed full of galaxies, long filamentary bridges, flattened sheets and vast low density voids. The study of the cosmic web has focused primarily on the identification of such features, and on understanding the environmental effects on galaxy formation and halo assembly. As such, a variety of different methods have been devised to classify the cosmic web -- depending on the data at hand, be it numerical simulations, large sky surveys or other. In this paper we bring twelve of these methods together and apply them to the same data set in order to understand how they compare. In general these cosmic web classifiers have been designed with different cosmological goals in mind, and to study different questions. Therefore one would not {\it a priori} expect agreement between different techniques however, many of these methods do converge on the identification of specific features. In this paper we study the agreements and disparities of the different methods. For example, each method finds that knots inhabit higher density regions than filaments, etc. and that voids have the lowest densities. For a given web environment, we find substantial overlap in the density range assigned by each web classification scheme. We also compare classifications on a halo-by-halo basis; for example, we find that 9 of 12 methods classify around a third of group-mass haloes (i.e. $M_{\rm halo}\sim10^{13.5}h^{-1}M_{\odot}$) as being in filaments. Lastly, so that any future cosmic web classification scheme can be compared to the 12 methods used here, we have made all the data used in this paper public.

Recent analytical work on the modelling of dark halo abundances and clustering has demonstrated the advantages of combining the excursion set approach with peaks theory. We extend these ideas and introduce a model of excursion set peaks that incorporates the role of initial tidal effects or shear in determining the gravitational collapse of dark haloes. The model -- in which the critical density threshold for collapse depends on the tidal influences acting on protohaloes -- is well motivated from ellipsoidal collapse arguments and is also simple enough to be analytically tractable. We show that the predictions of this model are in very good agreement with measurements of the halo mass function and traditional scale dependent halo bias in N-body simulations across a wide range of masses and redshift. The presence of shear in the collapse threshold means that halo bias is naturally predicted to be nonlocal, and that protohalo densities at fixed mass are naturally predicted to have Lognormal-like distributions. We present the first direct estimate of Lagrangian nonlocal bias in N-body simulations, finding broad agreement with the model prediction. Finally, the simplicity of the model (which has essentially a single free parameter) opens the door to building efficient and accurate non-universal fitting functions of halo abundances and bias for use in precision cosmology.

We revisit the issue of interpreting the results of large volume cosmological simulations in the context of large scale general relativistic effects. We look for simple modifications to the nonlinear evolution of the gravitational potential $\psi$ that lead on large scales to the correct, fully relativistic description of density perturbations in the Newtonian gauge. We note that the relativistic constraint equation for $\psi$ can be cast as a diffusion equation, with a diffusion length scale determined by the expansion of the Universe. Exploiting the weak time evolution of $\psi$ in all regimes of interest, this equation can be further accurately approximated as a Helmholtz equation, with an effective relativistic 'screening' scale $\ell$ related to the Hubble radius. We demonstrate that it is thus possible to carry out N-body simulations in the Newtonian gauge by replacing Poisson's equation with this Helmholtz equation, involving a trivial change in the Green's function kernel. Our results also motivate a simple, approximate (but very accurate) gauge transformation - $\delta_{\rm N}(\mathbf{k}) \approx \delta_{\rm sim}(\mathbf{k})\times (k^2+\ell^{-2})/k^2$ - to convert the density field $\delta_{\rm sim}$ of standard collisionless N-body simulations (initialised in the comoving synchronous gauge) into the Newtonian gauge density $\delta_{\rm N}$ at arbitrary times. A similar conversion can also be written in terms of particle positions. Our results can be interpreted in terms of a Jeans stability criterion induced by the expansion of the Universe. The appearance of the screening scale $\ell$ in the evolution of $\psi$, in particular, leads to a natural resolution of the 'Jeans swindle' in the presence of super-horizon modes.

We study the evolution of the stellar component and the metallicity of both the intracluster medium and of stars in massive ($M_{\rm vir}\approx 6\times 10^{14}$ M$_{\odot}/h$) simulated galaxy clusters from the Rhapsody-G suite in detail and compare them to observational results. The simulations were performed with the AMR code RAMSES and include the effect of AGN feedback at the sub-grid level. AGN feedback is required to produce realistic galaxy and cluster properties and plays a role in mixing material in the central regions and regulating star formation in the central galaxy. In both our low and high resolution runs with fiducial stellar yields, we find that stellar and ICM metallicities are a factor of two lower than in observations. We find that cool core clusters exhibit steeper metallicity gradients than non-cool core clusters, in qualitative agreement with observations. We verify that the ICM metallicities measured in the simulation can be explained by a simple "regulator" model in which the metallicity is set by a balance of stellar yield and gas accretion. It is plausible that a combination of higher resolution and higher metal yield in AMR simulation would allow the metallicity of simulated clusters to match observed values; however this hypothesis needs to be tested with future simulations. Comparison to recent literature highlights that results concerning the metallicity of clusters and cluster galaxies might depend sensitively on the scheme chosen to solve the hydrodynamics.

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.

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 present the Rhapsody-G suite of cosmological hydrodynamic AMR zoom simulations of ten massive galaxy clusters at the $M_{\rm vir}\sim10^{15}\,{\rm M}_\odot$ scale. These simulations include cooling and sub-resolution models for star formation and stellar and supermassive black hole feedback. The sample is selected to capture the whole gamut of assembly histories that produce clusters of similar final mass. We present an overview of the successes and shortcomings of such simulations in reproducing both the stellar properties of galaxies as well as properties of the hot plasma in clusters. In our simulations, a long-lived cool-core/non-cool core dichotomy arises naturally, and the emergence of non-cool cores is related to low angular momentum major mergers. Nevertheless, the cool-core clusters exhibit a low central entropy compared to observations, which cannot be alleviated by thermal AGN feedback. For cluster scaling relations we find that the simulations match well the $M_{500}-Y_{500}$ scaling of Planck SZ clusters but deviate somewhat from the observed X-ray luminosity and temperature scaling relations in the sense of being slightly too bright and too cool at fixed mass, respectively. Stars are produced at an efficiency consistent with abundance matching constraints and central galaxies have star formation rates consistent with recent observations. While our simulations thus match various key properties remarkably well, we conclude that the shortcomings strongly suggest an important role for non-thermal processes (through feedback or otherwise) or thermal conduction in shaping the intra-cluster medium.

Within a sufficiently large cosmic volume, conservation of baryons implies a simple `closed box' view in which the sum of the baryonic components must equal a constant fraction of the total enclosed mass. We present evidence from Rhapsody-G hydrodynamic simulations of massive galaxy clusters that the closed-box expectation may hold to a surprising degree within the interior, non-linear regions of haloes. At a fixed halo mass, we find a significant anti-correlation between hot gas mass fraction and galaxy mass fraction (cold gas + stars), with a rank correlation coefficient of -0.69 within $R_{500c}$. Because of this anti-correlation, the total baryon mass serves as a low-scatter proxy for total cluster mass. The fractional scatter of total baryon fraction scales approximately as $0.02 (\Delta_c/100)^{0.6}$, while the scatter of either gas mass or stellar mass is larger in magnitude and declines more slowly with increasing radius. We discuss potential observational tests using cluster samples selected by optical and hot gas properties; the simulations suggest that joint selection on stellar and hot gas has potential to achieve 5% scatter in total halo mass.

We study the angular-momentum (AM) buildup in high-$z$ massive galaxies using high-resolution cosmological simulations. The AM originates in co-planar streams of cold gas and merging galaxies tracing cosmic-web filaments, and it undergoes four phases of evolution. (I) Outside the halo virial radius ($R_{\rm v}\!\sim\!100\,{\rm kpc}$), the elongated streams gain AM by tidal torques with a specific AM (sAM) $\sim\!1.7$ times the dark-matter (DM) spin due to the gas' higher quadrupole moment. This AM is expressed as stream impact parameters, from $\sim\!0.3R_{\rm v}$ to counter rotation. (II) In the outer halo, while the incoming DM mixes with the existing halo of lower sAM to a spin $\lambda_{\rm dm}\!\sim\!0.04$, the cold streams transport the AM to the inner halo such that their spin in the halo is $\sim\!3\lambda_{\rm dm}$. (III) Near pericenter, the streams dissipate into an irregular rotating ring extending to $\sim\!0.3R_{\rm v}$ and tilted relative to the inner disc. Torques exerted partly by the disc make the ring gas lose AM, spiral in, and settle into the disc within one orbit. The ring is observable with 30\% probability as a damped Lyman-$\alpha$ absorber. (IV) Within the disc, $<\!0.1R_{\rm v}$, torques associated with violent disc instability drive AM out and baryons into a central bulge, while outflows remove low-spin gas, introducing certain sensitivity to feedback strength. Despite the different AM histories of gas and DM, the disc spin is comparable to the DM-halo spin. Counter rotation can strongly affect disc evolution.

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 progresses 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.

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

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.

We perform a suite of multimass cosmological zoom simulations of individual dark matter halos and explore how to best select Lagrangian regions for resimulation without contaminating the halo of interest with low-resolution particles. Such contamination can lead to significant errors in the gas distribution of hydrodynamical simulations, as we show. For a fixed Lagrange volume, we find that the chance of contamination increases systematically with the level of zoom. In order to avoid contamination, the Lagrangian volume selected for resimulation must increase monotonically with the resolution difference between parent box and the zoom region. We provide a simple formula for selecting Lagrangian regions (in units of the halo virial volume) as a function of the level of zoom required. We also explore the degree to which a halo's Lagrangian volume correlates with other halo properties (concentration, spin, formation time, shape, etc.) and find no significant correlation. There is a mild correlation between Lagrange volume and environment, such that halos living in the most clustered regions have larger Lagrangian volumes. Nevertheless, selecting halos to be isolated is not the best way to ensure inexpensive zoom simulations. We explain how one can safely choose halos with the smallest Lagrangian volumes, which are the least expensive to resimulate, without biasing one's sample.

We find that infalling dark matter halos (i.e., the progenitors of satellite halos) begin losing mass well outside the virial radius of their eventual host halos. The peak mass occurs at a range of clustercentric distances, with median and 68th percentile range of 1.8 +2.3/-1.0 R_(vir,host) for progenitors of z=0 satellites. The peak circular velocity for infalling halos occurs at significantly larger distances (3.7 +3.3/-2.2 R_(vir,host) at z=0). This difference arises because different physical processes set peak circular velocity (typically, ~1:5 and larger mergers which cause transient circular velocity spikes) and peak mass (typically, smooth accretion) for infalling halos. We find that infalling halos also stop having significant mergers well before they enter the virial radius of their eventual hosts. Mergers larger than a 1:40 ratio in halo mass end for infalling halos at similar clustercentric distances (~ 1.9 R_(vir,host)) as the end of overall mass accretion. However, mergers larger than 1:3 typically end for infalling halos at more than 4 virial radial away from their eventual hosts. This limits the ability of mergers to affect quenching and morphology changes in clusters. We also note that the transient spikes which set peak circular velocity may lead to issues with abundance matching on that parameter, including unphysical galaxy stellar mass growth profiles near clusters; we propose a simple observational test to check if a better halo proxy for galaxy stellar mass exists.

Warm Dark Matter (WDM) cosmologies are a viable alternative to the Cold Dark Matter (CDM) scenario. Unfortunately, an accurate scrutiny of the WDM predictions with N-body simulations has proven difficult due to numerical artefacts. Here, we report on cosmological simulations that, for the first time, are devoid of those problems, and thus, are able to accurately resolve the WDM halo mass function well below the cut-off. We discover a complex picture, with perturbations at different evolutionary stages populating different ranges in the halo mass function. On the smallest mass scales we can resolve, identified objects are typically centres of filaments that are starting to collapse. On intermediate mass scales, objects typically correspond to fluctuations that have collapsed and are in the process of relaxation, whereas the high mass end is dominated by objects similar to haloes identified in CDM simulations. We then explicitly show how the formation of low-mass haloes is suppressed, which translates into a strong cut-off in the halo mass function. This disfavours some analytic formulations that predict a halo mass function that would extend well below the free streaming mass. We argue for a more detailed exploration of the formation of the smallest structures expected to form in a given cosmology, which, we foresee, will advance our overall understanding of structure formation.

We investigate the large-scale clustering and gravitational interaction of baryons and dark matter (DM) over cosmic time using a set of collisionless N-body simulations. Both components, baryons and DM, are evolved from distinct primordial density and velocity power spectra as predicted by early-universe physics. We first demonstrate that such two-component simulations require an unconventional match between force and mass resolution (i.e. force softening on at least the mean particle separation scale). Otherwise, the growth on any scale is not correctly recovered because of a spurious coupling between the two species at the smallest scales. With these simulations, we then demonstrate how the primordial differences in the clustering of baryons and DM are progressively diminished over time. In particular, we explicitly show how the BAO signature is damped in the spatial distribution of baryons and imprinted in that of DM. This is a rapid process, yet it is still not fully completed at low redshifts. On large scales, the overall shape of the correlation function of baryons and DM differs by 2% at z = 9 and by 0.2% at z = 0. The differences in the amplitude of the BAO peak are approximately a factor of 5 larger: 10% at z = 9 and 1% at z = 0. These discrepancies are, however, smaller than effects expected to be introduced by galaxy formation physics in both the shape of the power spectrum and in the BAO peak, and are thus unlikely to be detected given the precision of the next generation of galaxy surveys. Hence, our results validate the standard practice of modelling the observed galaxy distribution using predictions for the total mass clustering in the Universe.

By combining galaxy tracers from high-resolution N-body and hydrodynamical simulations, we present a consistent picture of the behaviour of galaxy velocities in massive clusters. In haloes above ~ 10^14 Msun, the brightest satellite galaxies are slightly cooler compared to the dark matter, while fainter satellites are hotter. Within the virial radius of a cluster, the mean velocity dispersion based on the 100 brightest galaxies is a factor of 1.065 +/- 0.005 (stat) +/- 0.027 (sys) higher than that of the dark matter (corresponding to a ~10-15 per cent bias in the dynamical mass estimate) while that based on only the five brightest galaxies is 0.868 +/- 0.039 (stat) +/- 0.035 (sys). These trends are approximately independent of redshift. The velocity structure is sensitive to the modelling of galaxies in clusters, indicative of the complex interplay of tidal stripping, dynamical friction, and merging. Velocity dispersions derived from instantaneous subhalo properties are larger than those employing either peak subhalo properties or hydrodynamical galaxy tracers. The latter two methods are consistent, implying that stacked spectroscopic analysis of cluster samples should, after correction for projection, show a trend towards slightly higher velocities when fainter galaxies are included, with an unbiased measure of dark matter velocity dispersion coming from approximately 30 galaxies per cluster. We show evidence that the velocity distribution function of bright galaxies near the cluster centre has a low-velocity tail due to strong dynamical friction.

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

We present a semi-analytic, computationally inexpensive model to identify halos capable of forming a Population III star in cosmological simulations across a wide range of times and environments. This allows for a much more complete and representative set of Population III star forming halos to be constructed, which will lead to Population III star formation simulations that more accurately reflect the diversity of Population III stars, both in time and halo mass. This model shows that Population III and chemically enriched stars coexist beyond the formation of the first generation of stars in a cosmological simulation until at least z~10, and likely beyond, though Population III stars form at rates that are 4-6 orders of magnitude lower than chemically enriched stars by z=10. A catalog of more than 40,000 candidate Population III forming halos were identified, with formation times temporally ranging from z=30 to z=10, and ranging in mass from 2.3x10^5 M_sun to 1.2x10^10 M_sun. At early times, the environment that Population III stars form in is very similar to that of halos hosting chemically enriched star formation. At later times Population III stars are found to form in low-density regions that are not yet chemically polluted due to a lack of previous star formation in the area. Population III star forming halos become increasingly spatially isolated from one another at later times, and are generally closer to halos hosting chemically enriched star formation than to another halo hosting Population III star formation by z~10.

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.

We discuss the properties of subhalos in cluster-size halos, using a high-resolution statistical sample: the Rhapsody simulations introduced in Wu et al. (2012). We demonstrate that the criteria applied to select subhalos have significant impact on the inferred properties of the sample, including the scatter in the number of subhalos, the correlation between the subhalo number and formation time, and the shape of subhalos' spatial distribution and velocity structure. We find that the number of subhalos, when selected using the peak maximum circular velocity in their histories (a property expected to be closely related to the galaxy luminosity), is uncorrelated with the formation time of the main halo. This is in contrast to the previously reported correlation from studies where subhalos are selected by the current maximum circular velocity; we show that this difference is a result of the tidal stripping of the subhalos. We also find that the dominance of the main halo and the subhalo mass fraction are strongly correlated with halo concentration and formation history. These correlations are important to take into account when interpreting results from cluster samples selected with different criteria. Our sample also includes a fossil cluster, which is presented separately and placed in the context of the rest of the sample.

We present the first results from the Rhapsody cluster re-simulation project: a sample of 96 "zoom-in" simulations of dark matter halos of 10^14.8 +- 0.05 Msun/h, selected from a 1 (Gpc/h)^3 volume. This simulation suite is the first to resolve this many halos with ~5x10^6 particles per halo in the cluster-mass regime, allowing us to statistically characterize the distribution of and correlation between halo properties at fixed mass. We focus on the properties of the main halos and how they are affected by formation history, which we track back to z=12, over five decades in mass. We give particular attention to the impact of the formation history on the density profiles of the halos. We find that the deviations from the Navarro-Frenk-White (NFW) model and the Einasto model depend on formation time. Late-forming halos tend to have considerable deviations from both models, partly due to the presence of massive subhalos, while early-forming halos deviate less but still significantly from the NFW model and are better described by the Einasto model. We find that the halo shapes depend only moderately on formation time. Departure from spherical symmetry impacts the density profiles through the anisotropic distribution of massive subhalos. Further evidence of the impact of subhalos is provided by analyzing the phase-space structure. A detailed analysis of the properties of the subhalo population in Rhapsody is presented in a companion paper.

We examine the Velocity Distribution Function (VDF) in dark matter halos from Milky Way to cluster mass scales. We identify an empirical model for the VDF with a wider peak and a steeper tail than a Maxwell--Boltzmann distribution, and discuss physical explanations. We quantify sources of scatter in the VDF of cosmological halos and their implication for direct detection of dark matter. Given modern simulations and observations, we find that the most significant uncertainty in the VDF of the Milky Way arises from the unknown radial position of the solar system relative to the dark matter halo scale radius.

In the last decades cosmological N-body dark matter simulations have enabled ab initio studies of the formation of structure in the Universe. Gravity amplified small density fluctuations generated shortly after the Big Bang, leading to the formation of galaxies in the cosmic web. These calculations have led to a growing demand for methods to analyze time-dependent particle based simulations. Rendering methods for such N-body simulation data usually employ some kind of splatting approach via point based rendering primitives and approximate the spatial distributions of physical quantities using kernel interpolation techniques, common in SPH (Smoothed Particle Hydrodynamics)-codes. This paper proposes three GPU-assisted rendering approaches, based on a new, more accurate method to compute the physical densities of dark matter simulation data. It uses full phase-space information to generate a tetrahedral tessellation of the computational domain, with mesh vertices defined by the simulation's dark matter particle positions. Over time the mesh is deformed by gravitational forces, causing the tetrahedral cells to warp and overlap. The new methods are well suited to visualize the cosmic web. In particular they preserve caustics, regions of high density that emerge, when several streams of dark matter particles share the same location in space, indicating the formation of structures like sheets, filaments and halos. We demonstrate the superior image quality of the new approaches in a comparison with three standard rendering techniques for N-body simulation data.

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 velocity 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 study the feeding of massive galaxies at high redshift through streams from the cosmic web using the Mare Nostrum hydro-cosmological simulation. Our statistical sample consists of 350 dark-matter haloes of ~10^12 Msun at z = 2.5. We find that ~70% of the influx into the virial radius Rv is in narrow streams covering 10% of the virial shell. On average 64% of the stream influx is in one stream, and 95% is in three dominant streams. The streams that feed a massive halo tend to lie in a plane that extends from half to a few Rv, hereafter "the stream plane" (SP). The streams are typically embedded in a thin sheet of low-entropy gas, a Zel'dovich pancake, which carries ~20% of the influx into Rv. The filaments-in-a-plane configuration about the massive haloes at the nodes of the cosmic web differs from the large- scale structure of the web where the filaments mark the intersections of slanted sheets. The stream plane is only weakly aligned with the angular momentum (AM) near Rv, consistent with the fact that typically 80% of the AM is carried by one dominant stream. The galactic disc plane shows a weak tendency to be perpendicular to the large-scale SP, consistent with tidal-torque theory. Most interesting, the direction of the disc AM is only weakly correlated with the AM direction at Rv. This indicates a significant AM exchange at the interphase between streams and disc in the greater environment of the disc inside an "AM sphere" of radius ~0.3Rv . The required large torques are expected based on the perturbed morphology and kinematics within this interaction sphere. This AM exchange may or may not require a major modification of the standard disc modeling based on AM conservation, depending on the extent to which the amplitude of the disc AM is affected, which is yet to be studied.

We discuss a new algorithm to generate multi-scale initial conditions with multiple levels of refinements for cosmological "zoom-in" simulations. The method uses an adaptive convolution of Gaussian white noise with a real space transfer function kernel together with an adaptive multi-grid Poisson solver to generate displacements and velocities following first (1LPT) or second order Lagrangian perturbation theory (2LPT). The new algorithm achieves RMS relative errors of order 10^(-4) for displacements and velocities in the refinement region and thus improves in terms of errors by about two orders of magnitude over previous approaches. In addition, errors are localized at coarse-fine boundaries and do not suffer from Fourier-space induced interference ringing. An optional hybrid multi-grid and Fast Fourier Transform (FFT) based scheme is introduced which has identical Fourier space behaviour as traditional approaches. Using a suite of re-simulations of a galaxy cluster halo our real space based approach is found to reproduce correlation functions, density profiles, key halo properties and subhalo abundances with per cent level accuracy. Finally, we generalize our approach for two-component baryon and dark-matter simulations and demonstrate that the power spectrum evolution is in excellent agreement with linear perturbation theory. For initial baryon density fields, it is suggested to use the local Lagrangian approximation in order to generate a density field for mesh based codes that is consistent with Lagrangian perturbation theory instead of the current practice of using the Eulerian linearly scaled densities.

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, there 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.

We use ~8,600 >5e10 Msol COSMOS galaxies to study how the morphological mix of massive ellipticals, bulge-dominated disks, intermediate-bulge disks, bulge-less disks and irregular galaxies evolves from z=0.2 to z=1. The morphological evolution depends strongly on mass. At M>3e11 Msol, no evolution is detected in the morphological mix: ellipticals dominate since z=1, and the Hubble sequence has quantitatively settled down by this epoch. At the 1e11 Msol mass scale, little evolution is detected, which can be entirely explained with major mergers. Most of the morphological evolution from z=1 to z=0.2 takes place at masses 5e10 - 1e11 Msol, where: (i) The fraction of spirals substantially drops and the contribution of early-types increases. This increase is mostly produced by the growth of bulge-dominated disks, which vary their contribution from ~10% at z=1 to >30% at z=0.2 (cf. the elliptical fraction grows from ~15% to ~20%). Thus, at these masses, transformations from late- to early-types result in disk-less elliptical morphologies with a statistical frequency of only 30% - 40%. Otherwise, the processes which are responsible for the transformations either retain or produce a non-negligible disk component. (ii) The bulge-less disk galaxies, which contribute ~15% to the intermediate-mass galaxy population at z=1, virtually disappear by z=0.2. The merger rate since z=1 is too low to account for the disappearance of these massive bulge-less disks, which most likely grow a bulge via secular evolution.

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-halo's inertia tensors with the local tidal tensors and investigate how the correlation strength depends on the proto-halo's 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-halo's mass scale. The correlation strength is found to increase as the proto-halo's sphericity increases, as the proto-halo's 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 anisotropic 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 perform a series of high-resolution N-body simulations of cosmological structure formation starting from Gaussian and non-Gaussian initial conditions. We adopt the best-fitting cosmological parameters of WMAP (3rd- and 5th-year) and we consider non-Gaussianity of the local type parameterised by 8 different values of the non-linearity parameter F_NL. Building upon previous work based on the Gaussian case, we show that, expressed in terms of suitable variables, the mass function of friends-of-friends haloes is approximately universal (i.e. independent of redshift, cosmology, and matter transfer function) to good precision (nearly 10 per cent) also in non-Gaussian scenarios. We provide fitting formulae for the high-mass end (M>10^13 M_sol/h) of the universal mass function in terms of F_NL, and we also present a non-universal fit in terms of both F_NL and z to be used for applications requiring higher accuracy. In the Gaussian case, we extend our fit to a wider range of halo masses (M>2.4 x 10^10 M_sol/h) and we also provide a consistent fit of the linear halo bias. We show that, for realistic values of F_NL, the matter power-spectrum in non-Gaussian cosmologies departs from the Gaussian one by up to 2 per cent on the scales where the baryonic- oscillation features are imprinted on the 2-point statistics. We confirm the strong k-dependence of the halo bias on large scales (k<0.05 h Mpc^-1) which was already detected in previous studies. However, we find that commonly used parameterisations based on the peak-background split do not provide an accurate description of our simulations which present extra dependencies on the wavenumber, the non-linearity parameter and, possibly, the clustering strength. We provide an accurate fit of the simulation data that can be used as a benchmark for future determinations of F_NL with galaxy surveys.

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 the 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.

Despite a history that dates back at least a quarter of a century studies of voids in the large--scale structure of the Universe are bedevilled by a major problem: there exist a large number of quite different void--finding algorithms, a fact that has so far got in the way of groups comparing their results without worrying about whether such a comparison in fact makes sense. Because of the recent increased interest in voids, both in very large galaxy surveys and in detailed simulations of cosmic structure formation, this situation is very unfortunate. We here present the first systematic comparison study of thirteen different void finders constructed using particles, haloes, and semi--analytical model galaxies extracted from a subvolume of the Millennium simulation. The study includes many groups that have studied voids over the past decade. We show their results and discuss their differences and agreements. As it turns out, the basic results of the various methods agree very well with each other in that they all locate a major void near the centre of our volume. Voids have very underdense centres, reaching below 10 percent of the mean cosmic density. In addition, those void finders that allow for void galaxies show that those galaxies follow similar trends. For example, the overdensity of void galaxies brighter than $m_B = -20 $ is found to be smaller than about -0.8 by all our void finding algorithms.

We use a series of N-body simulations of the LCDM cosmology to investigate the redshift evolution since z=1 of the properties and alignment with the large-scale structure of haloes in clusters, filaments, sheets and voids. We find that: (i) Once a rescaling of the halo mass with M*, the mass scale collapsing at redshift z, is performed, there is no further significant redshift dependence in the halo properties; (ii) The environment influences halo shape and formation time at all investigated redshifts for haloes with masses M<M*; and (iii) There is a significant alignment of both spin and shape of haloes with filaments and sheets. In detail, at all redshifts up to z=1: a) Haloes with M<M* tend to be more oblate when located in clusters than in the other environments; this trend is reversed at higher masses: above about M*, halos in clusters are typically more prolate than similar massive haloes in sheets, filaments and voids. b) Haloes with M>M* in filaments spin more rapidly than similar mass haloes in clusters; haloes in voids have the lowest median spin parameters; c) Haloes with M<M* tend to be younger in voids and older in clusters; d) In sheets, halo spin vectors tend to lie within the sheet plane independent of mass; in filaments, instead, haloes with M<M* tend to spin parallel to the filament and haloes with M>M* perpendicular to it. For masses M>M*, the major axis of haloes in filaments and sheets is strongly aligned with the filament or the sheet. Such halo-LSS alignments may be of importance in weak lensing analyses of cosmic shear. A question that is opened by our study is why, in the 0 < z < 1 redshift regime that we have investigated, the mass scale M* sets roughly the threshold below which the LSS-environment either begins to affect, or reverses, the properties of dark matter haloes.

Using a series of high-resolution N-body simulations of the concordance cosmology we investigate how the formation histories, shapes and angular momenta of dark-matter haloes depend on environment. We first present a classification scheme that allows to distinguish between haloes in clusters, filaments, sheets and voids in the large-scale distribution of matter. This method is based on a local-stability criterion for the orbits of test particles and closely relates to the Zel'dovich approximation. Applying this scheme to our simulations we then find that: i) Mass assembly histories and formation redshifts strongly depend on environment for haloes of mass M<M* (haloes of a given mass tend to be older in clusters and younger in voids) and are independent of it for larger masses; ii) Low-mass haloes in clusters are generally less spherical and more oblate than in other regions; iii) Low-mass haloes in clusters have a higher median spin than in filaments and present a more prominent fraction of rapidly spinning objects; we identify recent major mergers as a likely source of this effect. For all these relations, we provide accurate functional fits as a function of halo mass and environment. We also look for correlations between halo-spin directions and the large-scale structures: the strongest effect is seen in sheets where halo spins tend to lie within the plane of symmetry of the mass distribution. Finally, we measure the spatial auto-correlation of spin directions and the cross-correlation between the directions of intrinsic and orbital angular momenta of neighbouring haloes. While the first quantity is always very small, we find that spin-orbit correlations are rather strong especially for low-mass haloes in clusters and high-mass haloes in filaments.