226 research outputs found
ColDICE: a parallel Vlasov-Poisson solver using moving adaptive simplicial tessellation
Resolving numerically Vlasov-Poisson equations for initially cold systems can
be reduced to following the evolution of a three-dimensional sheet evolving in
six-dimensional phase-space. We describe a public parallel numerical algorithm
consisting in representing the phase-space sheet with a conforming,
self-adaptive simplicial tessellation of which the vertices follow the
Lagrangian equations of motion. The algorithm is implemented both in six- and
four-dimensional phase-space. Refinement of the tessellation mesh is performed
using the bisection method and a local representation of the phase-space sheet
at second order relying on additional tracers created when needed at runtime.
In order to preserve in the best way the Hamiltonian nature of the system,
refinement is anisotropic and constrained by measurements of local Poincar\'e
invariants. Resolution of Poisson equation is performed using the fast Fourier
method on a regular rectangular grid, similarly to particle in cells codes. To
compute the density projected onto this grid, the intersection of the
tessellation and the grid is calculated using the method of Franklin and
Kankanhalli (1993) generalised to linear order. As preliminary tests of the
code, we study in four dimensional phase-space the evolution of an initially
small patch in a chaotic potential and the cosmological collapse of a
fluctuation composed of two sinusoidal waves. We also perform a "warm" dark
matter simulation in six-dimensional phase-space that we use to check the
parallel scaling of the code.Comment: Code and illustration movies available at:
http://www.vlasix.org/index.php?n=Main.ColDICE - Article submitted to Journal
of Computational Physic
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
The persistent cosmic web and its filamentary structure I: Theory and implementation
We present DisPerSE, a novel approach to the coherent multi-scale
identification of all types of astrophysical structures, and in particular the
filaments, in the large scale distribution of matter in the Universe. This
method and corresponding piece of software allows a genuinely scale free and
parameter free identification of the voids, walls, filaments, clusters and
their configuration within the cosmic web, directly from the discrete
distribution of particles in N-body simulations or galaxies in sparse
observational catalogues. To achieve that goal, the method works directly over
the Delaunay tessellation of the discrete sample and uses the DTFE density
computed at each tracer particle; no further sampling, smoothing or processing
of the density field is required.
The idea is based on recent advances in distinct sub-domains of computational
topology, which allows a rigorous application of topological principles to
astrophysical data sets, taking into account uncertainties and Poisson noise.
Practically, the user can define a given persistence level in terms of
robustness with respect to noise (defined as a "number of sigmas") and the
algorithm returns the structures with the corresponding significance as sets of
critical points, lines, surfaces and volumes corresponding to the clusters,
filaments, walls and voids; filaments, connected at cluster nodes, crawling
along the edges of walls bounding the voids. The method is also interesting as
it allows for a robust quantification of the topological properties of a
discrete distribution in terms of Betti numbers or Euler characteristics,
without having to resort to smoothing or having to define a particular scale.
In this paper, we introduce the necessary mathematical background and
describe the method and implementation, while we address the application to 3D
simulated and observed data sets to the companion paper.Comment: A higher resolution version is available at
http://www.iap.fr/users/sousbie together with complementary material.
Submitted to MNRA
Vlasov versus N-body: the H\'enon sphere
We perform a detailed comparison of the phase-space density traced by the
particle distribution in Gadget simulations to the result obtained with a
spherical Vlasov solver using the splitting algorithm. The systems considered
are apodized H\'enon spheres with two values of the virial ratio, R ~ 0.1 and
0.5. After checking that spherical symmetry is well preserved by the N-body
simulations, visual and quantitative comparisons are performed. In particular
we introduce new statistics, correlators and entropic estimators, based on the
likelihood of whether N-body simulations actually trace randomly the Vlasov
phase-space density. When taking into account the limits of both the N-body and
the Vlasov codes, namely collective effects due to the particle shot noise in
the first case and diffusion and possible nonlinear instabilities due to finite
resolution of the phase-space grid in the second case, we find a spectacular
agreement between both methods, even in regions of phase-space where nontrivial
physical instabilities develop. However, in the colder case, R=0.1, it was not
possible to prove actual numerical convergence of the N-body results after a
number of dynamical times, even with N=10 particles.Comment: 19 pages, 11 figures, MNRAS, in pres
The three dimensional skeleton: tracing the filamentary structure of the universe
The skeleton formalism aims at extracting and quantifying the filamentary
structure of the universe is generalized to 3D density fields; a numerical
method for computating a local approximation of the skeleton is presented and
validated here on Gaussian random fields. This method manages to trace well the
filamentary structure in 3D fields such as given by numerical simulations of
the dark matter distribution on large scales and is insensitive to monotonic
biasing. Two of its characteristics, namely its length and differential length,
are analyzed for Gaussian random fields. Its differential length per unit
normalized density contrast scales like the PDF of the underlying density
contrast times the total length times a quadratic Edgeworth correction
involving the square of the spectral parameter. The total length scales like
the inverse square smoothing length, with a scaling factor given by 0.21 (5.28+
n) where n is the power index of the underlying field. This dependency implies
that the total length can be used to constrain the shape of the underlying
power spectrum, hence the cosmology. Possible applications of the skeleton to
galaxy formation and cosmology are discussed. As an illustration, the
orientation of the spin of dark halos and the orientation of the flow near the
skeleton is computed for dark matter simulations. The flow is laminar along the
filaments, while spins of dark halos within 500 kpc of the skeleton are
preferentially orthogonal to the direction of the flow at a level of 25%.Comment: 17 pages, 11 figures, submitted to MNRA
Filaments in observed and mock galaxy catalogues
Context. The main feature of the spatial large-scale galaxy distribution is
an intricate network of galaxy filaments. Although many attempts have been made
to quantify this network, there is no unique and satisfactory recipe for that
yet. Aims. The present paper compares the filaments in the real data and in the
numerical models, to see if our best models reproduce statistically the
filamentary network of galaxies. Methods. We apply an object point process with
interactions (the Bisous process) to trace and describe the filamentary network
both in the observed samples (the 2dFGRS catalogue) and in the numerical models
that have been prepared to mimic the data.We compare the networks. Results. We
find that the properties of filaments in numerical models (mock samples) have a
large variance. A few mock samples display filaments that resemble the observed
filaments, but usually the model filaments are much shorter and do not form an
extended network. Conclusions. We conclude that although we can build numerical
models that are similar to observations in many respects, they may fail yet to
explain the filamentary structure seen in the data. The Bisous-built filaments
are a good test for such a structure.Comment: 13 pages, accepted for publication in Astronomy and Astrophysic
On the filamentary environment of galaxies
The correlation between the large-scale distribution of galaxies and their
spectroscopic properties at z=1.5 is investigated using the Horizon MareNostrum
cosmological run.
We have extracted a large sample of 10^5 galaxies from this large
hydrodynamical simulation featuring standard galaxy formation physics. Spectral
synthesis is applied to these single stellar populations to generate spectra
and colours for all galaxies. We use the skeleton as a tracer of the cosmic web
and study how our galaxy catalogue depends on the distance to the skeleton. We
show that galaxies closer to the skeleton tend to be redder, but that the
effect is mostly due to the proximity of large haloes at the nodes of the
skeleton, rather than the filaments themselves.
This effects translate into a bimodality in the colour distribution of our
sample. The origin of this bimodality is investigated and seems to follow from
the ram pressure stripping of satellite galaxies within the more massive
clusters of the simulation.
The virtual catalogues (spectroscopical properties of the MareNostrum
galaxies at various redshifts) are available online at
http://www.iap.fr/users/pichon/MareNostrum/cataloguesComment: 18 pages, 27 figures, accepted for publication in MNRA
- …