Source author record

Stephane Colombi

Stephane Colombi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph astro-ph.GA physics.flu-dyn

Catalog footprint

What is connected

8works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2025arXiv

Multiscale Turbulence Synthesis: Validation in 2D Hydrodynamics

Numerical simulations can follow the evolution of fluid motions through the intricacies of developed turbulence. However, they are rather costly to run, especially in 3D. In the past two decades, generative models have emerged which produce synthetic random flows at a computational cost equivalent to no more than a few time-steps of a simulation. These simplified models qualitatively bear some characteristics of turbulent flows in specific contexts (incompressible 3D hydrodynamics or magnetohydrodynamics), but generally struggle with the synthesis of coherent structures. We aim at generating random fields (e.g. velocity, density, magnetic fields, etc.) with realistic physical properties for a large variety of governing partial differential equations and at a small cost relative to time-resolved simulations. We propose a set of approximations applied to given sets of partial differential equations, and test the validity of our method in the simplest framework: 2D decaying incompressible hydrodynamical turbulence. We compare results of 2D decaying simulations with snapshots of our synthetic turbulence. We assess quantitatively the difference first with standard statistical tools: power spectra, increments and structure functions. These indicators can be reproduced by our method during up to about a third of the turnover time scale. We also consider recently developed scattering transforms statistics, able to efficiently characterise non-Gaussian structures. This reveals more significant discrepancy, which can however be bridged by bootstrapping. Finally, the number of Fourier transforms necessary for one synthesis scales logarithmically in the resolution, compared to linearly for time-resolved simulations. We have designed a multiscale turbulence synthesis (MuScaTS) method to efficiently short-circuit costly numerical simulations to produce realistic instantaneous fields.

preprint2014arXiv

Vlasov-Poisson in 1D for initially cold systems: post-collapse Lagrangian perturbation theory

We study analytically the collapse of an initially smooth, cold, self-gravitating collisionless system in one dimension. The system is described as a central "S" shape in phase-space surrounded by a nearly stationary halo acting locally like a harmonic background on the S. To resolve the dynamics of the S under its self-gravity and under the influence of the halo, we introduce a novel approach using post-collapse Lagrangian perturbation theory. This approach allows us to follow the evolution of the system between successive crossing times and to describe in an iterative way the interplay between the central S and the halo. Our theoretical predictions are checked against measurements in entropy conserving numerical simulations based on the waterbag method. While our post-collapse Lagrangian approach does not allow us to compute rigorously the long term behavior of the system, i.e. after many crossing times, it explains the close to power-law behavior of the projected density observed in numerical simulations. Pushing the model at late time suggests that the system could build at some point a very small flat core, but this is very speculative. This analysis shows that understanding the dynamics of initially cold systems requires a fine grained approach for a correct description of their very central part. The analyses performed here can certainly be extended to spherical symmetry.

preprint2011arXiv

Haloes gone MAD: The Halo-Finder Comparison Project

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

preprint2010arXiv

Cell Count Moments in the Halo Model

We study cell count moments up to fifth order of the distributions of haloes, of halo substructures as a proxy for galaxies, and of mass in the context of the halo model and compare theoretical predictions to the results of numerical simulations. On scales larger than the size of the largest cluster, we present a simple point cluster model in which results depend only on cluster-cluster correlations and on the distribution of the number of objects within a cluster, or cluster occupancy. The point cluster model leads to expressions for moments of galaxy counts in which the volume-averaged moments on large scales approach those of the halo distribution and on smaller scales exhibit hierarchical clustering with amplitudes $S_k$ determined by moments of the occupancy distribution. In this limit, the halo model predictions are purely combinatoric, and have no dependence on halo profile, concentration parameter, or potential asphericity. The full halo model introduces only two additional effects: on large scales, haloes of different mass have different clustering strengths, introducing relative bias parameters; and on the smallest scales, halo structure is resolved and details of the halo profile become important, introducing shape-dependent form factors. Because of differences between discrete and continuous statistics, the hierarchical amplitudes for galaxies and for mass behave differently on small scales even if galaxy number is exactly proportional to mass, a difference that is not necessarily well described in terms of bias.

preprint2009arXiv

The Skeleton: Connecting Large Scale Structures to Galaxy Formation

We report on two quantitative, morphological estimators of the filamentary structure of the Cosmic Web, the so-called global and local skeletons. The first, based on a global study of the matter density gradient flow, allows us to study the connectivity between a density peak and its surroundings, with direct relevance to the anisotropic accretion via cold flows on galactic halos. From the second, based on a local constraint equation involving the derivatives of the field, we can derive predictions for powerful statistics, such as the differential length and the relative saddle to extrema counts of the Cosmic web as a function of density threshold (with application to percolation of structures and connectivity), as well as a theoretical framework to study their cosmic evolution through the onset of gravity-induced non-linearities.

preprint2006arXiv

Baryon Dynamics, Dark Matter Substructure, and Galaxies

By comparing a collisionless cosmological N-body simulation (DM) to an SPH simulation with the same initial conditions, we investigate the correspondence between the dark matter subhalos produced by collisionless dynamics and the galaxies produced by dissipative gas dynamics in a dark matter background. When galaxies in the SPH simulation become satellites in larger groups, they retain local dark matter concentrations (SPH subhalos) whose mass is typically five times their baryonic mass. The more massive subhalos of the SPH simulation have corresponding subhalos of similar mass and position in the DM simulation; at lower masses, there is fairly good correspondence, but some DM subhalos are in different spatial positions and some suffer tidal stripping or disruption. The halo occupation statistics of DM subhalos -- the mean number of subhalos, pairs, and triples as a function of host halo mass -- are very similar to those of SPH subhalos and SPH galaxies. Gravity of the dissipative baryon component amplifies the density contrast of subhalos in the SPH simulation, making them more resistant to tidal disruption. Relative to SPH galaxies and SPH subhalos, the DM subhalo population is depleted in the densest regions of the most massive halos. The good agreement of halo occupation statistics between the DM subhalo and SPH galaxy populations leads to good agreement of their two-point correlation functions and higher order moments on large scales. The depletion of DM subhalos in dense regions depresses their clustering at R<1 Mpc/h. In these simulations, the "conversation" between dark matter and baryons is mostly one-way, with dark matter dynamics telling galaxies where to form and how to cluster, but the "back talk" of the baryons influences small scale clustering by enhancing the survival of substructure in the densest environments.

preprint1997arXiv

Non-Linear Evolution of the Bispectrum of Cosmological Perturbations

The bispectrum, the three-point function of density fluctuations in Fourier space, is the lowest order statistic that carries information about the spatial coherence of large-scale structures. For Gaussian initial conditions, when the density fluctuation amplitude is small, tree-level (leading order) perturbation theory predicts a characteristic dependence of the bispectrum on the shape of the triangle formed by the three wave vectors. This configuration dependence provides a signature of gravitational instability, and departures from it in galaxy catalogs can be interpreted as due to bias, that is, non-gravitational effects. On the other hand, N-body simulations indicate that the reduced three-point function becomes relatively shape-independent in the strongly non-linear regime. In order to understand this non-linear transition and assess the domain of reliability of shape-dependence as a probe of bias, we calculate the one-loop (next-to-leading order) corrections to the bispectrum in perturbation theory. We compare these results with measurements in numerical simulations with scale-free and CDM initial power spectra. We find that the one-loop corrections account very well for the departures from the tree-level results measured in numerical simulations on weakly non-linear scales. In this regime, the reduced bispectrum qualitatively retains its tree-level shape, but the amplitude can change significantly. At smaller scales, the reduced bispectrum in the simulations starts to flatten, an effect which can be partially understood from the one-loop results. In the strong clustering regime, where perturbation theory breaks down entirely, the simulation results confirm that the reduced bispectrum has almost no dependence on triangle shape, in rough agreement with the hierarchical ansatz.

preprint1995arXiv

Large Scale Structure Tests of Warm Dark Matter

Warm dark matter (WDM) is an intriguing model of structure formation from the point of view of both cosmology and particle physics. We consider a one-parameter family of WDM models. The linear power spectra for these models is calculated and compared with the corresponding spectra for cold dark matter (CDM), hot dark matter (HDM) and mixed dark matter (MDM) as well as the power spectrum derived from observations. Our linear analyses suggest that a model universe dominated by a particle whose mass to temperature ratio $m_x/T_x$ is increased by a factor of two as compared with the standard HDM neutrino gives a reasonable fit to the data on large $(>8h^{-1} {\rm ~Mpc})$ scales. $N$-body simulations for this particular WDM model show features of both HDM and CDM. As in HDM, the first objects to collapse are large pancake-like structures. The final matter distribution is rather smooth and structures as small as galaxy halos are excluded. However, there appear to be virialized rich clusters evident in the CDM but not the HDM simulations. Unfortunately, a simple comparison of the matter distribution and its statistical properties with observations indicates that WDM, like CDM, has too much power at small scales. This is particularly evident in the small-scale pairwize velocity dispersion. The cluster multiplicity function has the wrong shape with too many rich clusters being produced, though this conclusion is based on the simple assumption that light traces mass in groups of galaxies.