Advancing the matter bispectrum estimation of large-scale structure: fast prescriptions for galaxy mock catalogues.

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>2\times10^{14}\,h^{-1} M_\odot$ obtained a better than 1\% fit to both power spectrum and bispectrum in the range $K/3 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.