Novel Adaptive softening for collisionless N-body simulations: Eliminating spurious halos

Abstract

We describe a NOVel form of Adaptive softening (NovA) for collisionless $N$-body simulations, implemented in the Ramses adaptive mesh refinement code. We introduce a refinement criterion that the particle distribution within each cell be sufficiently isotropic, as measured by its moment of inertia tensor. In this way, collapse is only refined if it occurs along all three axes, ensuring that the softening $\epsilon$ is always of order twice the largest inter-particle spacing in a cell. This more conservative force softening criterion is designed to minimise spurious two-body effects, while maintaining high force resolution in collapsed regions of the flow. We test NovA using an antisymmetric perturbed plane wave collapse (`Valinia' test) before applying it to warm dark matter (WDM) simulations. For the Valinia test, we show that -- unlike the standard $N$-body method -- NovA produces no numerical fragmentation while still being able to correctly capture fine caustics and shells around the collapsing regions. For the WDM simulations, we find that NovA converges significantly more rapidly than standard $N$-body, producing little or no spurious halos on small scales. We show, however, that determining whether or not halos exist below the free streaming mass $M_{\mathrm{fs}}$ is complicated by the fact that our halo finder (AHF) likely incorrectly labels some caustics and criss-crossing filaments as halos, while one or two particularly massive filaments appear to fragment in any version of NovA where refinement is allowed. Such massive filaments may be physically unstable to collapse, as is the case for infinite, static, self-gravitating cylinders. We will use NovA in forthcoming papers to study the issue of halo formation below $M_{\mathrm{fs}}$; filament stability; and to obtain new constraints on the temperature of dark matter.

12 pages, 7 figures. To be submitted to MNRAS

References in corpus (10)

• Stellar Feedback in Dwarf Galaxy Formation
• Discreteness effects in simulations of Hot/Warm dark matter
• A simple multigrid scheme for solving the Poisson equation with arbitrary domain boundaries
• N-body simulations of gravitational dynamics
• The Cosmic Code Comparison Project
• An adaptively refined phase-space element method for cosmological simulations and collisionless dynamics
• Discreteness Effects in Lambda Cold Dark Matter Simulations: A Wavelet-Statistical View
• Bar-Halo Friction in Galaxies III: Halo Density Changes
• Towards quantitative control on discreteness error in the non-linear regime of cosmological N body simulations
• An optimum time-stepping scheme for N-body simulations

Cited by in corpus (16)

• Dark Matter Substructure in Numerical Simulations: A Tale of Discreteness Noise, Runaway Instabilities, and Artificial Disruption
• Large-scale dark matter simulations
• Simulating nonlinear cosmological structure formation with massive neutrinos
• Numerical convergence of simulations of galaxy formation: the abundance and internal structure of cold dark matter haloes
• Substructure and galaxy formation in the Copernicus Complexio warm dark matter simulations
• Hints against the cold and collisionless nature of dark matter from the galaxy velocity function
• Towards a general parametrization of the warm dark matter halo mass function
• Warm Dark Matter and Cosmic Reionization
• The lensing properties of subhaloes in massive elliptical galaxies in sterile neutrino cosmologies
• Spurious Small-Scale Structure & Discreteness-Driven Relaxation in Cosmological Simulations
• Dark matter halos in the multicomponent model. I. Substructure
• Positive feedback at the disc-halo interface
• An MHD spectral theory approach to Jeans' magnetized gravitational instability
• Analytic growth rate of gravitational instability in self-gravitating planar polytropes
• Tidal adaptive softening and artificial fragmentation in cosmological simulations
• Numerical Study of Bar Suppression in Galaxy Models Due to Disc Heating