Advancing the matter bispectrum estimation of large-scale structure: fast prescriptions for galaxy mock catalogues. (arXiv:1909.03248v1 [astro-ph.CO])
<a href="http://arxiv.org/find/astro-ph/1/au:+Hung_J/0/1/0/all/0/1">Johnathan Hung</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Manera_M/0/1/0/all/0/1">Marc Manera</a>, <a href="http://arxiv.org/find/astro-ph/1/au:+Shellard_E/0/1/0/all/0/1">E.P.S. Shellard</a>
We investigate various phenomenological schemes for the rapid generation of
3D mock galaxy catalogues with a given power spectrum and bispectrum. We apply
the fast bispectrum estimator MODALLSS{} to these mock galaxy catalogues and
compare to $N$-body simulation data analysed with the halo-finder
texttt{ROCKSTAR} (our benchmark data). We propose an assembly bias model for
populating parent halos with subhalos by using a joint lognormal-Gaussian
probability distribution for the subhalo occupation number and the halo
concentration. This prescription enabled us to recover the benchmark power
spectrum from $N$-body simulations to within 1% and the bispectrum to within
4% across the entire range of scales of the simulation. A small further boost
adding an extra galaxy to all parent halos above the mass threshold
$M>2times10^{14},h^{-1} M_odot$ obtained a better than 1% fit to both power
spectrum and bispectrum in the range $K/3<1.1,h,text{Mpc}^{-1}$, where
$K=k_1+k_2+k_3$. This statistical model should be applicable to fast dark
matter codes, allowing rapid generation of mock catalogues which simultaneously
reproduce the halo power spectrum and bispectrum obtained from $N$-body
simulations. We also investigate alternative schemes using the Halo Occupation
Distribution (HOD) which depend only on halo mass, but these yield results
deficient in both the power spectrum (2%) and the bispectrum (>4%) at $k,K/3
approx 0.2,h,text{Mpc}^{-1}$, with poor scaling for the latter. Efforts to
match the power spectrum by modifying the standard four-parameter HOD model
result in overboosting the bispectrum (with a 10% excess). We also
characterise the effect of changing the halo profile on the power spectrum and
bispectrum.
We investigate various phenomenological schemes for the rapid generation of
3D mock galaxy catalogues with a given power spectrum and bispectrum. We apply
the fast bispectrum estimator MODALLSS{} to these mock galaxy catalogues and
compare to $N$-body simulation data analysed with the halo-finder
texttt{ROCKSTAR} (our benchmark data). We propose an assembly bias model for
populating parent halos with subhalos by using a joint lognormal-Gaussian
probability distribution for the subhalo occupation number and the halo
concentration. This prescription enabled us to recover the benchmark power
spectrum from $N$-body simulations to within 1% and the bispectrum to within
4% across the entire range of scales of the simulation. A small further boost
adding an extra galaxy to all parent halos above the mass threshold
$M>2times10^{14},h^{-1} M_odot$ obtained a better than 1% fit to both power
spectrum and bispectrum in the range $K/3<1.1,h,text{Mpc}^{-1}$, where
$K=k_1+k_2+k_3$. This statistical model should be applicable to fast dark
matter codes, allowing rapid generation of mock catalogues which simultaneously
reproduce the halo power spectrum and bispectrum obtained from $N$-body
simulations. We also investigate alternative schemes using the Halo Occupation
Distribution (HOD) which depend only on halo mass, but these yield results
deficient in both the power spectrum (2%) and the bispectrum (>4%) at $k,K/3
approx 0.2,h,text{Mpc}^{-1}$, with poor scaling for the latter. Efforts to
match the power spectrum by modifying the standard four-parameter HOD model
result in overboosting the bispectrum (with a 10% excess). We also
characterise the effect of changing the halo profile on the power spectrum and
bispectrum.
http://arxiv.org/icons/sfx.gif
