The Adaptive TreePM: An Adaptive Resolution Code for Cosmological N-body Simulations

Abstract

Cosmological N-Body simulations are used for a variety of applications. Indeed progress in the study of large scale structures and galaxy formation would have been very limited without this tool. For nearly twenty years the limitations imposed by computing power forced simulators to ignore some of the basic requirements for modeling gravitational instability. One of the limitations of most cosmological codes has been the use of a force softening length that is much smaller than the typical inter-particle separation. This leads to departures from collisionless evolution that is desired in these simulations. We propose a particle based method with an adaptive resolution where the force softening length is reduced in high density regions while ensuring that it remains well above the local inter-particle separation. The method, called the Adaptive TreePM, is based on the TreePM code. We present the mathematical model and an implementation of this code, and demonstrate that the results converge over a range of options for parameters introduced in generalizing the code from the TreePM code. We explicitly demonstrate collisionless evolution in collapse of an oblique plane wave. We compare the code with the fixed resolution TreePM code and also an implementation that mimics adaptive mesh refinement methods and comment on the agreement, and disagreements in the results. We find that in most respects the ATreePM code performs at least as well as the fixed resolution TreePM in highly over-dense regions, from clustering and number density of haloes, to internal dynamics of haloes. We also show that the adaptive code is faster than the corresponding high resolution TreePM code.

18 pages, 11 figures. Accepted for publication in the MNRAS

References in corpus (11)

• Wilkinson Microwave Anisotropy Probe (WMAP) Three Year Results: Implications for Cosmology
• The shape of the SDSS DR5 galaxy power spectrum
• The Cosmic Code Comparison Project
• Simulation techniques for cosmological simulations
• Towards optimal parallel PM N-body codes: PMFAST
• Discreteness Effects in Lambda Cold Dark Matter Simulations: A Wavelet-Statistical View
• Towards quantitative control on discreteness error in the non-linear regime of cosmological N body simulations
• Comments on the size of the simulation box in cosmological N-Body simulations
• Effects of the size of cosmological N-Body simulations on physical quantities - III: Skewness
• On the Onset of Stochasticity in $\Lambda$CDM Cosmological Simulations
• Gravitational collapse in an expanding background and the role of substructure II: Excess power at small scales and its effect of collapse of structures at larger scales

Cited by in corpus (14)

• FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation
• Phantom: A smoothed particle hydrodynamics and magnetohydrodynamics code for astrophysics
• Haloes gone MAD: The Halo-Finder Comparison Project
• Smoothed Particle Hydrodynamics in Astrophysics
• Large-scale dark matter simulations
• A new approach to simulating collisionless dark matter fluids
• An adaptively refined phase-space element method for cosmological simulations and collisionless dynamics
• Tree-based solvers for adaptive mesh refinement code FLASH - I: gravity and optical depths
• Adaptive gravitational softening in GADGET
• Quasar Clustering in Cosmological Hydrodynamic Simulations: Evidence for mergers
• PHoToNs--A Parallel Heterogeneous & Threads oriented code for cosmological N-body simulation
• The effect of softening on dynamical simulations of galaxies
• Simulating the collapse of rotating primordial gas clouds to study the survival possibility of Pop III protostars
• Effect of the cubic torus topology on cosmological perturbations