Sampling Artifact in Volume Weighted Velocity Measurement.--- II. Detection in simulations and comparison with theoretical modelling

Measuring the volume weighted velocity power spectrum suffers from a severe systematic error, due to imperfect sampling of the velocity field from inhomogeneous distribution of dark matter particles/halos in simulations or galaxies with velocity measurement. This "sampling artifact" depends on both the mean particle number density $\bar{n}_P$ and the intrinsic large scale structure (LSS) fluctuation in the particle distribution. (1) We report robust detection of this sampling artifact in N-body simulations. It causes $\sim 12$% underestimation of the velocity power spectrum at $k=0.1$h/Mpc for samples with $\bar{n}_P=6\times10^{-3}$ (Mpc/h)$^{-3}$. This systematic underestimation increases with decreasing $\bar{n}_P$ and increasing $k$. Its dependence on the intrinsic LSS fluctuations is also robustly detected. (2) All these findings are expected by our theoretical modelling in paper I \cite{Zhang14}. In particular, the leading order theoretical approximation agrees quantitatively well with simulation result for $\bar{n}_P\gtrsim6\times 10^{-4}$(Mpc/h)$^{-3}$. Furthermore, we provide an ansatz to take high order terms into account. It improves the model accuracy to $\lesssim1$% at $k\lesssim0.1$h/Mpc over 3 orders of magnitude in $\bar{n}_P$ and over typical LSS clustering from $z=0$ to $z=2$. (3) The sampling artifact is determined by the deflection ${\bf D}$ field, which is straightforwardly available in both simulations and data of galaxy velocity. Hence the sampling artifact in the velocity power spectrum measurement can be self-calibrated within our framework. By applying such self-calibration in simulations, it becomes promising to determine the {\it real} large scale velocity bias of $10^{13}M_\odot$ halos with $\sim 1$% accuracy, and that of lower mass halos by better accuracy. ...[abridged]Comment: 11 pages, 11 figures. More arguments added, match the PRD accepted versio

Kriging Interpolating Cosmic Velocity Field

[abridged] Volume-weighted statistics of large scale peculiar velocity is preferred by peculiar velocity cosmology, since it is free of uncertainties of galaxy density bias entangled in mass-weighted statistics. However, measuring the volume-weighted velocity statistics from galaxy (halo/simulation particle) velocity data is challenging. For the first time, we apply the Kriging interpolation to obtain the volume-weighted velocity field. Kriging is a minimum variance estimator. It predicts the most likely velocity for each place based on the velocity at other places. We test the performance of Kriging quantified by the E-mode velocity power spectrum from simulations. Dependences on the variogram prior used in Kriging, the number $n_k$ of the nearby particles to interpolate and the density $n_P$ of the observed sample are investigated. First, we find that Kriging induces $1\%$ and $3\%$ systematics at $k\sim 0.1h{\rm Mpc}^{-1}$ when $n_P\sim 6\times 10^{-2} ({\rm Mpc}/h)^{-3}$ and $n_P\sim 6\times 10^{-3} ({\rm Mpc}/h)^{-3}$, respectively. The deviation increases for decreasing $n_P$ and increasing $k$. When $n_P\lesssim 6\times 10^{-4} ({\rm Mpc}/h)^{-3}$, a smoothing effect dominates small scales, causing significant underestimation of the velocity power spectrum. Second, increasing $n_k$ helps to recover small scale power. However, for $n_P\lesssim 6\times 10^{-4} ({\rm Mpc}/h)^{-3}$ cases, the recovery is limited. Finally, Kriging is more sensitive to the variogram prior for lower sample density. The most straightforward application of Kriging on the cosmic velocity field does not show obvious advantages over the nearest-particle method (Zheng et al. 2013) and could not be directly applied to cosmology so far. However, whether potential improvements may be achieved by more delicate versions of Kriging is worth further investigation.Comment: 11 pages, 5 figures, published in PR