Galaxy bias renormalization: Two-loop Power Spectrum, One-loop Trispectrum and Bispectrum

We present a complete treatment of fifth-order renormalized galaxy bias at the one- and two-loop level in real space, including gradient corrections to deterministic bias operators at next-to-leading order. We then provide a complete computation of the two-loop power spectrum as well as the one-loop bispectrum and trispectrum of biased tracers, and demonstrate how to jointly model these statistics in a fully renormalized framework. These statistics also require stochastic renormalization of products of two, three or four operators at coincidence, which we include at leading order in gradients by means of an operator product expansion. We verify that all UV limits of loop integrals are absorbed by the counterterms we consider. Upon solving the resulting renormalization group equations, we find a pronounced scale-dependence of higher-gradient bias coefficients. Since our renormalization prescription is performed manifestly at the operator level, our results can also easily be extended to higher $N$-point functions, higher loop orders and field-level analyses.