Solving Large Scale Structure in Ten Easy Steps with COLA

We present the COmoving Lagrangian Acceleration (COLA) method: an N-body method for solving for Large Scale Structure (LSS) in a frame that is comoving with observers following trajectories calculated in Lagrangian Perturbation Theory (LPT). Unlike standard N-body methods, the COLA method can straightforwardly trade accuracy at small-scales in order to gain computational speed without sacrificing accuracy at large scales. This is especially useful for cheaply generating large ensembles of accurate mock halo catalogs required to study galaxy clustering and weak lensing, as those catalogs are essential for performing detailed error analysis for ongoing and future surveys of LSS. As an illustration, we ran a COLA-based N-body code on a box of size 100Mpc/h with particles of mass ~5*10^9Msolar/h. Running the code with only 10 timesteps was sufficient to obtain an accurate description of halo statistics down to halo masses of at least 10^11Msolar/h. This is only at a modest speed penalty when compared to mocks obtained with LPT. A standard detailed N-body run is orders of magnitude slower than our COLA-based code. The speed-up we obtain with COLA is due to the fact that we calculate the large-scale dynamics exactly using LPT, while letting the N-body code solve for the small scales, without requiring it to capture exactly the internal dynamics of halos. Achieving a similar level of accuracy in halo statistics without the COLA method requires at least 3 times more timesteps than when COLA is employed.Comment: 18 pages, 7 figure

Cosmological Perturbations at Second Order and Recombination Perturbed

We derive the full set of second-order equations governing the evolution of cosmological perturbations, including the effects of the first-order electron number density perturbations, \delta_e. We provide a detailed analysis of the perturbations to the recombination history of the universe and show that a perturbed version of the Peebles effective 3-level atom is sufficient for obtaining the evolution of \delta_e for comoving wavenumbers smaller than 1Mpc^{-1}. We calculate rigorously the perturbations to the Ly\alpha escape probability and show that to a good approximation it is governed by the local baryon velocity divergence. For modes shorter than the photon diffusion scale, we find that \delta_e is enhanced during recombination by a factor of roughly 5 relative to other first-order quantities sourcing the CMB anisotropies at second order. Using these results, in a companion paper we calculate the CMB bispectrum generated during recombination.Comment: 47 pages, 6 figure

Estimating CDM Particle Trajectories in the Mildly Non-Linear Regime of Structure Formation. Implications for the Density Field in Real and Redshift Space

We obtain approximations for the CDM particle trajectories starting from Lagrangian Perturbation Theory. These estimates for the CDM trajectories result in approximations for the density in real and redshift space, as well as for the momentum density that are better than what standard Eulerian and Lagrangian perturbation theory give. For the real space density, we find that our proposed approximation gives a good cross-correlation (>95%) with the non-linear density down to scales almost twice smaller than the non-linear scale, and six times smaller than the corresponding scale obtained using linear theory. This allows for a speed-up of an order of magnitude or more in the scanning of the cosmological parameter space with N-body simulations for the scales relevant for the baryon acoustic oscillations. Possible future applications of our method include baryon acoustic peak reconstruction, building mock galaxy catalogs, momentum field reconstruction.Comment: 25 pages, 11 figures; reference adde

Towards an Optimal Reconstruction of Baryon Oscillations

The Baryon Acoustic Oscillations (BAO) in the large-scale structure of the universe leave a distinct peak in the two-point correlation function of the matter distribution. That acoustic peak is smeared and shifted by bulk flows and non-linear evolution. However, it has been shown that it is still possible to sharpen the peak and remove its shift by undoing the effects of the bulk flows. We propose an improvement to the standard acoustic peak reconstruction. Contrary to the standard approach, the new scheme has no free parameters, treats the large-scale modes consistently, and uses optimal filters to extract the BAO information. At redshift of zero, the reconstructed linear matter power spectrum leads to a markedly improved sharpening of the reconstructed acoustic peak compared to standard reconstruction.Comment: 20 pages, 5 figures; footnote adde

Kinematic Density Waves in Accretion Disks

When thin accretion disks around black holes are perturbed, the main restoring force is gravity. If gas pressure, magnetic stresses, and radiation pressure are neglected, the disk remains thin as long as orbits do not intersect. Intersections would result in pressure forces which limit the growth of perturbations. We find that a discrete set of perturbations is possible for which orbits remain non-intersecting for arbitrarily long times. These modes define a discrete set of frequencies. We classify all long-lived perturbations for arbitrary potentials and show how their mode frequencies are related to pattern speeds computed from the azimuthal and epicyclic frequencies. We show that modes are concentrated near radii where the pattern speed has vanishing radial derivative. We explore these modes around Kerr black holes as a possible explanation for the high-frequency quasi-periodic oscillations of black hole binaries such as GRO J1655-40. The long-lived modes are shown to coincide with diskoseismic waves in the limit of small sound speed. While the waves have long lifetime, they have the wrong frequencies to explain the pairs of high-frequency quasi-periodic oscillations observed in black hole binaries.Comment: 28 pages, 6 figures; extended comparison with diskoseismology; added reference to astro-ph/060368