Dynamical flows through Dark Matter Haloes II: one and two points statistics at the virial radius

In a serie of three papers, the dynamical interplay between environments and dark matter haloes is investigated, while focussing on the dynamical flows through their virial sphere. Our method relies on both cosmological simulations, to constrain the environments, and an extension to the classical matrix method to derive the response of the halo (see Pichon & Aubert (2006), paper I). The current paper focuses on the statistical characterisation of the environments surrounding haloes, using a set of large scale simulations. Our description relies on a `fluid' halocentric representation where the interactions between the halo and its environment are investigated in terms of a time dependent external tidal field and a source term characterizing the infall. The method is applied to 15000 haloes, with masses between 5 x 10^12 Ms and 10^14 Ms evolving between z = 1 and z = 0. The net accretion at the virial radius is found to decrease with time, resulting from both an absolute decrease of infall and from a growing contribution of outflows. Infall is found to be mainly radial and occurring at velocities ~ 0.75 V200. Outflows are also detected through the virial sphere and occur at lower velocities ~ 0.6 V200 on more circular orbits. The external tidal field is found to be strongly quadrupolar and mostly stationnary, possibly reflecting the distribution of matter in the halo's near environment. The coherence time of the small scale fluctuations of the potential hints a possible anisotropic distribution of accreted satellites. The flux density of mass on the virial sphere appears to be more clustered than the potential while the shape of its angular power spectrum seems stationnary.Comment: 34 pages, 29 figures, accepted for publication in MNRA

Statistics of cosmic density profiles from perturbation theory

The joint probability distribution function (PDF) of the density within multiple concentric spherical cells is considered. It is shown how its cumulant generating function can be obtained at tree order in perturbation theory as the Legendre transform of a function directly built in terms of the initial moments. In the context of the upcoming generation of large-scale structure surveys, it is conjectured that this result correctly models such a function for finite values of the variance. Detailed consequences of this assumption are explored. In particular the corresponding one-cell density probability distribution at finite variance is computed for realistic power spectra, taking into account its scale variation. It is found to be in agreement with $\Lambda$-CDM simulations at the few percent level for a wide range of density values and parameters. Related explicit analytic expansions at the low and high density tails are given. The conditional (at fixed density) and marginal probability of the slope -- the density difference between adjacent cells -- and its fluctuations is also computed from the two-cells joint PDF; it also compares very well to simulations, in particular in under-dense regions, with a significant reduced cosmic scatter compared to over-dense regions. It is emphasized that this could prove useful when studying the statistical properties of voids as it can serve as a statistical indicator to test gravity models and/or probe key cosmological parameters.Comment: 22 pages, 15 figures, submitted to PR

Non Gaussian Minkowski functionals and extrema counts for 2D sky maps

In the conference presentation we have reviewed the theory of non-Gaussian geometrical measures for the 3D Cosmic Web of the matter distribution in the Universe and 2D sky data, such as Cosmic Microwave Background (CMB) maps that was developed in a series of our papers. The theory leverages symmetry of isotropic statistics such as Minkowski functionals and extrema counts to develop post- Gaussian expansion of the statistics in orthogonal polynomials of invariant descriptors of the field, its first and second derivatives. The application of the approach to 2D fields defined on a spherical sky was suggested, but never rigorously developed. In this paper we present such development treating effects of the curvature and finiteness of the spherical space $S_2$ exactly, without relying on the flat-sky approximation. We present Minkowski functionals, including Euler characteristic and extrema counts to the first non-Gaussian correction, suitable for weakly non-Gaussian fields on a sphere, of which CMB is the prime example.Comment: 6 pages, to appear as proceedings of the IAU Symposium No. 308, 2014 The Zeldovich Universe, Genesis and Growth of the Cosmic Web Rien van de Weygaert, Sergei Shandarin, Enn Saar and Jaan Einast

Dynamical flows through Dark Matter Haloes: Inner perturbative dynamics, secular evolution, and applications

We investigate statistically the dynamical consequences of cosmological fluxes of matter and related moments on progenitors of today's dark matter haloes. Their dynamics is described via canonical perturbation theory which accounts for two types of perturbations: the tidal field corresponding to fly-bys and accretion of dark matter through the halo's outer boundary. he dynamical equations are solved linearly, order by order, projecting on a biorthogonal basis to consistently satisfy the field equation. Since our solution of the Boltzmann Poisson equations is explicit, it allows statistical predictions for the ensemble distribution of the inner dynamical features of haloes. The secular evolution of open galactic haloes is investigated: we derive the kinetic equation which governs the quasi-linear evolution of dark matter profile induced by infall and its corresponding gravitational correlations. This yields a Fokker Planck-like equation for the angle-averaged underlying distribution function. We show how these extensions to the classical theory could be used to (i) observationally constrain the statistical nature of the infall (ii) predict the observed distribution and correlations of substructures in upcoming surveys, (iii) predict the past evolution of the observed distribution of clumps, and finally (iv) weight the relative importance of the intrinsic (via the unperturbed distribution function) and external (tidal and/ or infall) influence of the environment in determining the fate of galaxies.Comment: 35 pages, 12 Postscript figures, accepted for publication by MNRA

The invariant joint distribution of a stationary random field and its derivatives: Euler characteristic and critical point counts in 2 and 3D

The full moments expansion of the joint probability distribution of an isotropic random field, its gradient and invariants of the Hessian is presented in 2 and 3D. It allows for explicit expression for the Euler characteristic in ND and computation of extrema counts as functions of the excursion set threshold and the spectral parameter, as illustrated on model examples.Comment: 4 pages, 2 figures. Corrected expansion coefficients for orders n>=5. Relation between Gram-Charlier and Edgeworth expansions is clarified

Secular resonant dressed orbital diffusion II : application to an isolated self similar tepid galactic disc

The main orbital signatures of the secular evolution of an isolated self-gravitating stellar Mestel disc are recovered using a dressed Fokker-Planck formalism in angle-action variables. The shot-noise-driven formation of narrow ridges of resonant orbits is recovered in the WKB limit of tightly wound transient spirals, for a tepid Toomre-stable tapered disc. The relative effect of the bulge, the halo, the disc temperature and the spectral properties of the shot noise are investigated in turn. For such galactic discs all elements seem to impact the locus and direction of the ridge. For instance, when the halo mass is decreased, we observe a transition between a regime of heating in the inner regions of the disc through the inner Lindblad resonance to a regime of radial migration of quasi-circular orbits via the corotation resonance in the outer part of the disc. The dressed secular formalism captures both the nature of collisionless systems (via their natural frequencies and susceptibility), and their nurture via the structure of the external perturbing power spectrum. Hence it provides the ideal framework in which to study their long term evolution.Comment: 15 pages, 11 figure

The cosmic evolution of massive black holes in the Horizon-AGN simulation

We analyse the demographics of black holes (BHs) in the large-volume cosmological hydrodynamical simulation Horizon-AGN. This simulation statistically models how much gas is accreted onto BHs, traces the energy deposited into their environment and, consequently, the back-reaction of the ambient medium on BH growth. The synthetic BHs reproduce a variety of observational constraints such as the redshift evolution of the BH mass density and the mass function. Strong self-regulation via AGN feedback, weak supernova feedback, and unresolved internal processes result in a tight BH-galaxy mass correlation. Starting at z~2, tidal stripping creates a small population of BHs over-massive with respect to the halo. The fraction of galaxies hosting a central BH or an AGN increases with stellar mass. The AGN fraction agrees better with multi-wavelength studies, than single-wavelength ones, unless obscuration is taken into account. The most massive halos present BH multiplicity, with additional BHs gained by ongoing or past mergers. In some cases, both a central and an off-centre AGN shine concurrently, producing a dual AGN. This dual AGN population dwindles with decreasing redshift, as found in observations. Specific accretion rate and Eddington ratio distributions are in good agreement with observational estimates. The BH population is dominated in turn by fast, slow, and very slow accretors, with transitions occurring at z=3 and z=2 respectively.Comment: Accepted for publication in MNRA

Peak exclusion, stochasticity and convergence of perturbative bias expansions in 1+1 gravity

The Lagrangian peaks of a 1D cosmological random field representing dark matter are used as a proxy for a catalogue of biased tracers in order to investigate the small-scale exclusion in the two-halo term. The two-point correlation function of peaks of a given height is numerically estimated and analytical approximations that are valid inside the exclusion zone are derived. The resulting power spectrum of these tracers is investigated and shows clear deviations from Poisson noise at low frequencies. On large scales, the convergence of a perturbative bias expansion is discussed. Finally, we go beyond Gaussian statistics for the initial conditions and investigate the subsequent evolution of the two-point clustering of peaks through their Zel'dovich ballistic displacement, to clarify how exclusion effects mix up with scale-dependencies induced by nonlinear gravitational evolution. While the expected large-scale separation limit is recovered, significant deviations are found in the exclusion zone that tends in particular to be reduced at later times. Even though these findings apply to the clustering of one-dimensional tracers, they provide useful insights into halo exclusion and its impact on the two-halo term.Comment: 16 pages, 9 figures, accepted for publication in MNRA

The persistent cosmic web and its filamentary structure II: Illustrations

The recently introduced discrete persistent structure extractor (DisPerSE, Soubie 2010, paper I) is implemented on realistic 3D cosmological simulations and observed redshift catalogues (SDSS); it is found that DisPerSE traces equally well the observed filaments, walls, and voids in both cases. In either setting, filaments are shown to connect onto halos, outskirt walls, which circumvent voids. Indeed this algorithm operates directly on the particles without assuming anything about the distribution, and yields a natural (topologically motivated) self-consistent criterion for selecting the significance level of the identified structures. It is shown that this extraction is possible even for very sparsely sampled point processes, as a function of the persistence ratio. Hence astrophysicists should be in a position to trace and measure precisely the filaments, walls and voids from such samples and assess the confidence of the post-processed sets as a function of this threshold, which can be expressed relative to the expected amplitude of shot noise. In a cosmic framework, this criterion is comparable to friend of friend for the identifications of peaks, while it also identifies the connected filaments and walls, and quantitatively recovers the full set of topological invariants (Betti numbers) {\sl directly from the particles} as a function of the persistence threshold. This criterion is found to be sufficient even if one particle out of two is noise, when the persistence ratio is set to 3-sigma or more. The algorithm is also implemented on the SDSS catalogue and used to locat interesting configurations of the filamentary structure. In this context we carried the identification of an ``optically faint'' cluster at the intersection of filaments through the recent observation of its X-ray counterpart by SUZAKU. The corresponding filament catalogue will be made available online.Comment: A higher resolution version is available at http://www.iap.fr/users/sousbie together with complementary material (movie and data). Submitted to MNRA

Functional integral derivation of the kinetic equation of two-dimensional point vortices

We present a brief derivation of the kinetic equation describing the secular evolution of point vortices in two-dimensional hydrodynamics, by relying on a functional integral formalism. We start from Liouville's equation which describes the exact dynamics of a two-dimensional system of point vortices. At the order ${1/N}$, the evolution of the system is characterised by the first two equations of the BBGKY hierarchy involving the system's 1-body distribution function and its 1-body correlation function. Thanks to the introduction of auxiliary fields, these two evolution constraints may be rewritten as a functional integral. When functionally integrated over the 2-body correlation function, this rewriting leads to a new constraint coupling the 1-body distribution function and the two auxiliary fields. Once inverted, this constraint provides, through a new route, the closed non-linear kinetic equation satisfied by the 1-body distribution function. Such a method sheds new lights on the origin of these kinetic equations complementing the traditional derivation methods.Comment: 7 pages, 1 figur