We introduce a new code for cosmological simulations, PHoToNs, which has features on performing massive cosmological simulations on heterogeneous high performance Computer (HPC) and threads oriented programming. PHoToNs adopts a hybrid scheme to compute gravity force, with the conventional PM to compute the long-range force, the Tree algorithm to compute the short range force, and the direct summation PP to compute the gravity from very close particles. A self-similar space filling Peano-Hilbert curve is used to decompose computing domain. Threads programming is highly used to more flexibly manage the domain communication, PM calculation and synchronization, as well as Dual Tree Traversal on the CPU+MIC platform. The scalability of the PHoToNs performs well and the efficiency of PP kernel achieves 68.6% of peak performance on MIC and 74.4% on CPU platforms. We also test the accuracy of the code against the much used Gadget-2 in the community and found excellent agreement.

Analysis of the Pangu N-body simulation validates that the bulk flow of halos follows a Maxwellian distribution which variance is consistent with the prediction of the linear theory of structure formation. We propose that the consistency between the observed bulk velocity and theories should be examined at the effective scale of the radius of a spherical top-hat window function yielding the same smoothed velocity variance in linear theory as the sample window function does. We compared some recently estimated bulk flows from observational samples with the prediction of the \Lambda CDM model we used; some results deviate from expectation at a level of ~ 3\sigma but the discrepancy is not as severe as previously claimed. We show that bulk flow is only weakly correlated with the dipole of the internal mass distribution, the alignment angle between the mass dipole and the bulk flow has a broad distribution peaked at ~ 30-50 deg., and also that the bulk flow shows little dependence on the mass of the halos used in the estimation. In a simulation of box size 1Gpc/h, for a cell of radius 100 Mpc/h the maximal bulk velocity is >500 km/s, dipoles of the environmental mass outside the cell are not tightly aligned with the bulk flow, but are rather located randomly around it with separation angles ~ 20-40 deg. In the fastest cell there is a slightly smaller number of low-mass halos; however halos inside are clustered more strongly at scales > ~ 20 Mpc/h, which might be a significant feature since the correlation between bulk flow and halo clustering actually increases in significance beyond such scales.