sCOLA: The N-body COLA Method Extended to the Spatial Domain

We present sCOLA -- an extension of the N-body COmoving Lagrangian Acceleration (COLA) method to the spatial domain. Similar to the original temporal-domain COLA, sCOLA is an N-body method for solving for large-scale structure in a frame that is comoving with observers following trajectories calculated in Lagrangian Perturbation Theory. Incorporating the sCOLA method in an N-body code allows one to gain computational speed by capturing the gravitational potential from the far field using perturbative techniques, while letting the N-body code solve only for the near field. The far and near fields are completely decoupled, effectively localizing gravity for the N-body side of the code. Thus, running an N-body code for a small simulation volume using sCOLA can reproduce the results of a standard N-body run for the same small volume embedded inside a much larger simulation. We demonstrate that sCOLA can be safely combined with the original temporal-domain COLA. sCOLA can be used as a method for performing zoom-in simulations. It also allows N-body codes to be made embarrassingly parallel, thus allowing for efficiently tiling a volume of interest using grid computing. Moreover, sCOLA can be useful for cheaply generating large ensembles of accurate mock halo catalogs required to study galaxy clustering. Surveys that will benefit the most are ones with large aspect ratios, such as pencil-beam surveys, where sCOLA can easily capture the effects of large-scale transverse modes without the need to substantially increase the simulated volume. As an illustration of the method, we present proof-of-concept zoom-in simulations using a freely available sCOLA-based N-body code.

[1]  B. Wandelt,et al.  Bayesian analysis of the dynamic cosmic web in the SDSS galaxy survey , 2015, 1502.02690.

[2]  P. Fosalba,et al.  nIFTy cosmology: galaxy/halo mock catalogue comparison project on clustering statistics , 2014, 1412.7729.

[3]  Chris Power,et al.  halogen: a tool for fast generation of mock halo catalogues , 2014, 1412.5228.

[4]  W. Percival,et al.  The clustering of the SDSS main galaxy sample – II. Mock galaxy catalogues and a measurement of the growth of structure from redshift space distortions at z = 0.15 , 2014, 1409.3238.

[5]  Ashley J. Ross,et al.  The clustering of the SDSS DR7 Main Galaxy Sample I: a 4 per cent distance measure at z=0.15 , 2014, 1409.3242.

[6]  Cheng Zhao,et al.  EZmocks: Extending the Zel'dovich approximation to generate mock galaxy catalogues with accurate clustering statistics , 2014, 1409.1124.

[7]  F. Kitaura,et al.  Constraining the halo bispectrum in real and redshift space from perturbation theory and non-linear stochastic bias , 2014, 1407.1236.

[8]  W. Percival,et al.  The clustering of the SDSS DR 7 main Galaxy sample – I . A 4 per cent distance measure at z = 0 , 2015 .

[9]  B. Wandelt,et al.  Dark matter voids in the SDSS galaxy survey , 2014, 1410.0355.

[10]  Scott Croom,et al.  The WiggleZ Dark Energy Survey: improved distance measurements to z = 1 with reconstruction of the baryonic acoustic feature , 2014, 1401.0358.

[11]  S. Tassev.  Lagrangian or Eulerian; real or Fourier? Not all approaches to large-scale structure are created equal , 2013, 1311.4884.

[12]  Martin White,et al.  Mock galaxy catalogues using the quick particle mesh method , 2013, 1309.5532.

[13]  M. Zaldarriaga,et al.  The Lagrangian-space Effective Field Theory of large scale structures , 2013, 1311.2168.

[14]  A. Jenkins.  A new way of setting the phases for cosmological multiscale Gaussian initial conditions. , 2013, 1306.5968.

[15]  S. Borgani,et al.  An accurate tool for the fast generation of dark matter halo catalogues. , 2013, 1305.1505.

[16]  Matias Zaldarriaga,et al.  Solving large scale structure in ten easy steps with COLA , 2013, 1301.0322.

[17]  B. Reid,et al.  Convolution Lagrangian perturbation theory for biased tracers , 2012, 1209.0780.

[18]  J. Brinkmann,et al.  The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey:a large sample of mock galaxy catalogues , 2012, 1203.6609.

[19]  Benjamin Dan Wandelt,et al.  One-point remapping of Lagrangian perturbation theory in the mildly non-linear regime of cosmic structure formation , 2013 .

[20]  L. Senatore,et al.  The effective field theory of cosmological large scale structures , 2012, 1206.2926.

[21]  M. Zaldarriaga,et al.  Estimating CDM particle trajectories in the mildly non-linear regime of structure formation. Implications for the density field in real and redshift space , 2012, 1203.5785.

[22]  Oliver Hahn,et al.  Multi-scale initial conditions for cosmological simulations , 2011, 1103.6031.

[23]  I. Szapudi,et al.  FAST GENERATION OF ENSEMBLES OF COSMOLOGICAL N-BODY SIMULATIONS VIA MODE RESAMPLING , 2011, 1103.2767.

[24]  M. White,et al.  EMBEDDING REALISTIC SURVEYS IN SIMULATIONS THROUGH VOLUME REMAPPING , 2010, 1003.3178.

[25]  Simon D. M. White,et al.  One simulation to fit them all - changing the background parameters of a cosmological N-body simulation , 2009, 0912.4277.

[26]  M. Crocce,et al.  Transients from initial conditions in cosmological simulations , 2006, astro-ph/0606505.

[27]  R. Sheth,et al.  PTHALOS: a fast method for generating mock galaxy distributions , 2001, astro-ph/0106120.

[28]  S. Colombi,et al.  Perturbative Lagrangian Approach to Gravitational Instability , 1994, astro-ph/9406013.