Modeling the BAO feature in Bispectrum

We investigate how well a simple leading order perturbation theory model of the bispectrum can fit the BAO feature in the measured bispectrum monopole of galaxies. Previous works showed that perturbative models of galaxy bispectrum start failing at the wavenumbers of k ~ 0.1 Mpc/h. We show that when the BAO feature in the bispectrum is separated it can be successfully modeled up to much higher wavenumbers. We validate our modeling on GLAM simulations that were run with and without the BAO feature in the initial conditions. We also quantify the amount of systematic error due to BAO template being offset from the true cosmology. We find that the systematic errors do not exceed 0.3 per cent for reasonable deviations from the true cosmology.