  Upcoming imaging surveys will allow for high signal-to-noise measurements of
galaxy clustering at small scales. In this work, we present the results of the
LSST bias challenge, the goal of which is to compare the performance of
different nonlinear galaxy bias models in the context of LSST Y10 data.
Specifically, we compare two perturbative approaches, Lagrangian perturbation
theory (LPT) and Eulerian PT (EPT) to two variants of Hybrid Effective Field
Theory (HEFT), with our fiducial implementation of these models including terms
up to second order in the bias expansion as well as nonlocal bias and
deviations from Poissonian stochasticity. We consider different simulated
galaxy samples and test the performance of the bias models in a tomographic
joint analysis of LSST-Y10-like galaxy clustering, galaxy-galaxy-lensing and
cosmic shear. We find both HEFT methods as well as LPT and EPT combined with
non-perturbative predictions for the matter power spectrum to yield unbiased
constraints on cosmological parameters up to at least a maximal scale of
$k_{\mathrm{max}}=0.4 \; \mathrm{Mpc}^{-1}$ for all samples considered, even in
the presence of assembly bias. While we find that we can reduce the complexity
of the bias model for HEFT without compromising fit accuracy, this is not
generally the case for the perturbative models. We find significant detections
of non-Poissonian stochasticity in all cases considered, and our analysis shows
evidence that small-scale galaxy clustering predominantly improves constraints
on galaxy bias rather than cosmological parameters. These results therefore
suggest that the systematic uncertainties associated with current nonlinear
bias models are likely to be subdominant compared to other sources of error for
tomographic analyses of upcoming photometric surveys, which bodes well for
future galaxy clustering analyses using these high signal-to-noise data.
[abridged]
