add power_spectrum; openmp for update_particles