We compute the cosmological boost factor at high redshifts of $z=10$-$100$ by integrating non-linear matter power spectrum measured from high-resolution $N$-body simulations. An accurate is required to estimate energy injection dark annihilation, which may affect re-ionization process. combined various box-size simulations (side lengths $1 \, {\rm kpc}$-$10 Mpc}$) cover a wide range scales, i.e...