Source author record

Stéphane Colombi

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

astro-ph.CO astro-ph astro-ph.GA Computational Geometry physics.comp-ph physics.flu-dyn

Catalog footprint

Research overview

8works

6topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2016arXiv

Modeling Lyman-α Forest Cross-Correlations with LyMAS

We use the Ly- $α$ Mass Association Scheme (LyMAS; Peirani et al. 2014) to predict cross-correlations at $z=2.5$ between dark matter halos and transmitted flux in the Ly- $α$ forest, and compare to cross-correlations measured for quasars and damped Ly- $α$ systems (DLAs) from the Baryon Oscillation Spectroscopic Survey (BOSS) by Font-Ribera et al. (2012, 2013). We calibrate LyMAS using Horizon-AGN hydrodynamical cosmological simulations of a $(100\ h^{-1}\ \mathrm{Mpc})^3$ comoving volume. We apply this calibration to a $(1\ h^{-1}\ \mathrm{Gpc})^3$ simulation realized with $2048^3$ dark matter particles. In the 100 $h^{-1}$ Mpc box, LyMAS reproduces the halo-flux correlations computed from the full hydrodynamic gas distribution very well. In the 1 $h^{-1}$ Gpc box, the amplitude of the large scale cross-correlation tracks the halo bias $b_h$ as expected. We provide empirical fitting functions that describe our numerical results. In the transverse separation bins used for the BOSS analyses, LyMAS cross-correlation predictions follow linear theory accurately down to small scales. Fitting the BOSS measurements requires inclusion of random velocity errors; we find best-fit RMS velocity errors of 399 km s $^{-1}$ and 252 km s $^{-1}$ for quasars and DLAs, respectively. We infer bias-weighted mean halo masses of $M_h/10^{12}\ h^{-1}M_\odot=2.19^{+0.16}_{-0.15}$ and $0.69^{+0.16}_{-0.14}$ for the host halos of quasars and DLAs, with $\sim 0.2$ dex systematic uncertainty associated with redshift evolution, IGM parameters, and selection of data fitting range.

preprint2015arXiv

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

preprint2014arXiv

Dancing in the dark: galactic properties trace spin swings along the cosmic web

A large-scale hydrodynamical cosmological simulation, Horizon-AGN, is used to investigate the alignment between the spin of galaxies and the cosmic filaments above redshift 1.2. The analysis of more than 150 000 galaxies per time step in the redshift range 1.2<z<1.8 with morphological diversity shows that the spin of low-mass blue galaxies is preferentially aligned with their neighbouring filaments, while high-mass red galaxies tend to have a perpendicular spin. The reorientation of the spin of massive galaxies is provided by galaxy mergers, which are significant in their mass build-up. We find that the stellar mass transition from alignment to misalignment happens around 3.10^10 M_sun. Galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as active galactic nuclei is mandatory to quench in situ star formation in massive galaxies and promote various morphologies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping spins misaligned with cosmic filaments.

preprint2014arXiv

LyMAS: Predicting Large-Scale Lyman-alpha Forest Statistics from the Dark Matter Density Field

[abridged] We describe LyMAS (Ly-alpha Mass Association Scheme), a method of predicting clustering statistics in the Ly-alpha forest on large scales from moderate resolution simulations of the dark matter distribution, with calibration from high-resolution hydrodynamic simulations of smaller volumes. We use the "Horizon MareNostrum" simulation, a 50 Mpc/h comoving volume evolved with the adaptive mesh hydrodynamic code RAMSES, to compute the conditional probability distribution P(F_s|delta_s) of the transmitted flux F_s, smoothed (1-dimensionally) over the spectral resolution scale, on the dark matter density contrast delta_s, smoothed (3-dimensionally) over a similar scale. In this study we adopt the spectral resolution of the SDSS-III BOSS at z=2.5, and we find optimal results for a dark matter smoothing length sigma=0.3 Mpc/h (comoving). In extended form, LyMAS exactly reproduces both the 1-dimensional power spectrum and 1-point flux distribution of the hydro simulation spectra. Applied to the MareNostrum dark matter field, LyMAS accurately predicts the 2-point conditional flux distribution and flux correlation function of the full hydro simulation for transverse sightline separations as small as 1 Mpc/h, including redshift-space distortion effects. It is substantially more accurate than a deterministic density-flux mapping ("Fluctuating Gunn-Peterson Approximation"), often used for large volume simulations of the forest. With the MareNostrum calibration, we apply LyMAS to 1024^3 N-body simulations of a 300 Mpc/h and 1.0 Gpc/h cube to produce large, publicly available catalogs of mock BOSS spectra that probe a large comoving volume. LyMAS will be a powerful tool for interpreting 3-d Ly-alpha forest data, thereby transforming measurements from BOSS and other massive quasar absorption surveys into constraints on dark energy, dark matter, space geometry, and IGM physics.

preprint2014arXiv

Vlasov-Poisson in 1D: waterbags

We revisit in one dimension the waterbag method to solve numerically Vlasov-Poisson equations. In this approach, the phase-space distribution function $f(x,v)$ is initially sampled by an ensemble of patches, the waterbags, where $f$ is assumed to be constant. As a consequence of Liouville theorem it is only needed to follow the evolution of the border of these waterbags, which can be done by employing an orientated, self-adaptive polygon tracing isocontours of $f$ . This method, which is entropy conserving in essence, is very accurate and can trace very well non linear instabilities as illustrated by specific examples. As an application of the method, we generate an ensemble of single waterbag simulations with decreasing thickness, to perform a convergence study to the cold case. Our measurements show that the system relaxes to a steady state where the gravitational potential profile is a power-law of slowly varying index $β$ , with $β$ close to $3/2$ as found in the literature. However, detailed analysis of the properties of the gravitational potential shows that at the center, $β> 1.54$ . Moreover, our measurements are consistent with the value $β=8/5=1.6$ that can be analytically derived by assuming that the average of the phase-space density per energy level obtained at crossing times is conserved during the mixing phase. These results are incompatible with the logarithmic slope of the projected density profile $β-2 \simeq -0.47$ obtained recently by Schulz et al. (2013) using a $N$ -body technique. This sheds again strong doubts on the capability of $N$ -body techniques to converge to the correct steady state expected in the continuous limit.

preprint2012arXiv

Non Gaussianity and Minkowski Functionals: forecasts for Planck

We study Minkowski Functionals as probes of primordial non-Gaussianity in the Cosmic Microwave Background, specifically for the estimate of the primordial `local' bi-spectrum parameter f_NL, with instrumental parameters which should be appropriate for the Planck experiment. We use a maximum likelihood approach, which we couple with various filtering methods and test thoroughly for convergence. We included the effect of inhomogeneous noise as well as astrophysical biases induced by point sources and by the contamination from the Galaxy. We find that, when Wiener filtered maps are used (rather than simply smoothed with Gaussian), the expected error on the measurement of f_NL should be as small as Δf_NL \simeq 10 when combining the 3 channels at 100, 143 and 217 GHz in the Planck extended mission setup. This result is fairly insensitive to the non homogeneous nature of the noise, at least for realistic hit-maps expected from Planck. We then estimate the bias induced on the measurement of f_NL by point sources in those 3 channels. With the appropriate masking of the bright sources, this bias can be reduced to a negligible level in the 100 and 143 GHz channels. It remains significant in the 217 GHz channel, but can be corrected for. The galactic foreground biases are quite important and present a complex dependence on sky coverage: making them negligible will depend strongly on the quality of the component separation methods.

preprint2011arXiv

Most massive halos with Gumbel Statistics

We present an analytical calculation of the extreme value statistics for dark matter halos - that is, the probability distribution of the most massive halo within some region of the universe of specified shape and size. Our calculation makes use of the counts-in-cells formalism for the correlation functions, and the halo bias derived from the Sheth-Tormen mass function. We demonstrate the power of the method on spherical regions, comparing the results to measurements in a large cosmological dark matter simulation and achieving good agreement. Particularly good fits are obtained for the most likely value of the maximum mass and for the high-mass tail of the distribution, relevant in constraining cosmologies by observations of most massive clusters.

preprint2008arXiv

Cosmic velocity--gravity relation in redshift space

We propose a simple way to estimate the parameter beta = Omega_m^(0.6)/b from three-dimensional galaxy surveys. Our method consists in measuring the relation between the cosmological velocity and gravity fields, and thus requires peculiar velocity measurements. The relation is measured *directly in redshift space*, so there is no need to reconstruct the density field in real space. In linear theory, the radial components of the gravity and velocity fields in redshift space are expected to be tightly correlated, with a slope given, in the distant observer approximation, by g / v = (1 + 6 beta / 5 + 3 beta^2 / 7)^(1/2) / beta. We test extensively this relation using controlled numerical experiments based on a cosmological N-body simulation. To perform the measurements, we propose a new and rather simple adaptive interpolation scheme to estimate the velocity and the gravity field on a grid. One of the most striking results is that nonlinear effects, including `fingers of God', affect mainly the tails of the joint probability distribution function (PDF) of the velocity and gravity field: the 1--1.5 sigma region around the maximum of the PDF is *dominated by the linear theory regime*, both in real and redshift space. This is understood explicitly by using the spherical collapse model as a proxy of nonlinear dynamics. Applications of the method to real galaxy catalogs are discussed, including a preliminary investigation on homogeneous (volume limited) `galaxy' samples extracted from the simulation with simple prescriptions based on halo and sub-structure identification, to quantify the effects of the bias between the galaxy and the total matter distibution, and of shot noise (ABRIDGED).

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Stéphane Colombi

Research overview

Is this your work?

8 publications

Modeling Lyman-α Forest Cross-Correlations with LyMAS

ColDICE: a parallel Vlasov-Poisson solver using moving adaptive simplicial tessellation

Dancing in the dark: galactic properties trace spin swings along the cosmic web

LyMAS: Predicting Large-Scale Lyman-alpha Forest Statistics from the Dark Matter Density Field

Vlasov-Poisson in 1D: waterbags

Non Gaussianity and Minkowski Functionals: forecasts for Planck

Most massive halos with Gumbel Statistics

Cosmic velocity--gravity relation in redshift space

See the researcher in context

Building this map preview