Subhalo abundance matching in $f(R)$ gravity
  Using the liminality N-body simulations of Shi et. al., we present the first
predictions for galaxy clustering in $f(R)$ gravity using subhalo abundance
matching. We find that, for a given galaxy density, even for an $f(R)$ model
with $f_{R0}=-10^{-6}$, for which the cold dark matter clustering is very
similar to $\Lambda$CDM, the predicted clustering of galaxies in the $f(R)$
model is very different from $\Lambda$CDM. The deviation can be as large as
$40\%$ for samples with mean densities close to that of $L_*$ galaxies. This
large deviation is testable given the accuracy that future large-scale galaxy
surveys aim to achieve. Our result demonstrates that galaxy surveys can provide
a stringent test of General Relativity on cosmological scales, which is
comparable to the tests from local astrophysical observations.
