Renormalization Bias in Modern Theories

• Renormalization bias is the process by which nominal bias parameters are redefined or corrected via operator mixing and UV regularization in various scientific contexts.
• Methodologies like effective-theory closure, operator-level renormalization, and generalized eigenvalue problems are employed to absorb or mitigate small-scale sensitivities.
• Applications range from cosmological halo clustering and inflationary perturbation theory to bias correction in text embeddings, visual-inertial estimation, and deliberate bias potentials in Monte Carlo simulations.

Renormalization bias denotes several distinct but structurally related notions in contemporary research. In large-scale structure, it is the statement that the coefficients appearing in a halo or galaxy bias expansion are not directly observable because composite operators mix under renormalization, so physically meaningful bias parameters are renormalized couplings or finite-cutoff running coefficients rather than bare expansion coefficients. In inflation, it is the question of whether renormalization ambiguities can induce an observable distortion of the primordial scalar spectrum; the reported answer is that inflation exponentially suppresses the renormalization sector on observable scales. In representation learning and visual-inertial estimation, “renormalization” instead names explicit bias-reduction procedures that remove a shared mean direction in embeddings or the inherent statistical bias of linear least-squares initializers. A further usage appears in Monte Carlo renormalization-group methods, where a variational bias potential is introduced deliberately to reshape the coarse-grained ensemble rather than to denote an unwanted distortion (Nishizawa et al., 2012, Assassi et al., 2014, Romero et al., 18 Jun 2026, Ren et al., 14 Nov 2025, Micusik et al., 2020, Wu et al., 2017).

1. Renormalized bias in halo clustering

A foundational formulation appears in perturbative halo clustering, where the halo overdensity is expanded as a local functional of the matter overdensity and a stochastic term,

$$\delta_{\rm h}(\mathbf{k},z)= \sum_n \frac{b_n}{n!}\int d^3q_1\cdots d^3q_n\, \delta_D^3\!\left(\mathbf{k}-\sum_i\mathbf{q}_i\right) \delta_m(\mathbf{q}_1,z)\cdots\delta_m(\mathbf{q}_n,z)+\epsilon(\mathbf{k}).$$

At one loop, the formal halo-matter cross-spectrum contains a divergent contribution proportional to

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$

so the coefficient multiplying $P_m^{\rm L}(k)$ is not a physical large-scale bias. The renormalization step absorbs this small-scale sensitivity into an effective linear bias, yielding

$$P_{\rm hm}(k;M,z) \equiv b_1^{\rm eff} P_m^{\rm NL}(k;z) + b_2(M)\int\frac{d^3q}{(2\pi)^3} P_m^{\rm L}(q)P_m^{\rm L}(|\mathbf{k}-\mathbf{q}|) F_2(\mathbf{q},\mathbf{k}-\mathbf{q}),$$

while the halo auto-spectrum requires an additional constant residual shot-noise term $\delta N$. In this framework, $b_1^{\rm eff}$ is not identified with the theoretical halo bias $b_1(M)$; it is a renormalized effective bias fit from the large-scale cross-spectrum, and the explicit $b_2$-dependent loop terms generate scale-dependent bias in the weakly nonlinear regime (Nishizawa et al., 2012).

A complementary formulation emphasizes operator mixing directly. Composite operators such as $\delta^2$ and $s^2$ are built from fields evaluated at the same point, so their low-$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$0 contributions receive UV contamination from short-wavelength modes. For the simple expansion

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$1

one has $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$2, and the halo-matter cross-spectrum acquires a contribution

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$3

which mixes the quadratic operator into the apparent linear bias. This is removed by defining renormalized operators such as

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$4

so that the observable linear bias becomes

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$5

A non-perturbative implementation in $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$6-body simulations defines renormalized cross-spectra by subtracting the measured low-$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$7 projection onto the linear density field. With this procedure, the renormalized spectra become nearly independent of $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$8 on large scales. Bayesian model selection for $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$9 at $P_m^{\rm L}(k)$0 and $P_m^{\rm L}(k)$1 favors the operator set

$P_m^{\rm L}(k)$2

while next-to-leading-order perturbative solutions are reported to be inaccurate for $P_m^{\rm L}(k)$3 (Werner et al., 2019).

2. Effective-theory closure, regularization schemes, and RG flow

In the effective-theory treatment of halo bias, renormalization is not only a redefinition of coefficients but also a constraint on the operator basis. The naive Eulerian ansatz

$P_m^{\rm L}(k)$4

is not stable under coarse-graining, because loop contributions from short scales generate mixing with other operators allowed by the symmetries, including non-local tidal and velocity operators. The renormalized theory is instead written as

$P_m^{\rm L}(k)$5

with a complete operator basis through quartic order in the density contrast and at leading order in derivatives. At quadratic order, this requires $P_m^{\rm L}(k)$6 and

$P_m^{\rm L}(k)$7

At cubic order, the independent set includes

$P_m^{\rm L}(k)$8

where

$P_m^{\rm L}(k)$9

The paper further states that Galileon operators such as $$P_{\rm hm}(k;M,z) \equiv b_1^{\rm eff} P_m^{\rm NL}(k;z) + b_2(M)\int\frac{d^3q}{(2\pi)^3} P_m^{\rm L}(q)P_m^{\rm L}(|\mathbf{k}-\mathbf{q}|) F_2(\mathbf{q},\mathbf{k}-\mathbf{q}),$$0 and $$P_{\rm hm}(k;M,z) \equiv b_1^{\rm eff} P_m^{\rm NL}(k;z) + b_2(M)\int\frac{d^3q}{(2\pi)^3} P_m^{\rm L}(q)P_m^{\rm L}(|\mathbf{k}-\mathbf{q}|) F_2(\mathbf{q},\mathbf{k}-\mathbf{q}),$$1 are not renormalized at leading order in derivatives, so they organize the EFT basis naturally (Assassi et al., 2014).

Regularization-scheme dependence enters explicitly when renormalizing composite operators in the galaxy bias expansion with primordial local non-Gaussianity and nonlinear gravitational evolution. For the operator $$P_{\rm hm}(k;M,z) \equiv b_1^{\rm eff} P_m^{\rm NL}(k;z) + b_2(M)\int\frac{d^3q}{(2\pi)^3} P_m^{\rm L}(q)P_m^{\rm L}(|\mathbf{k}-\mathbf{q}|) F_2(\mathbf{q},\mathbf{k}-\mathbf{q}),$$2, three UV prescriptions are compared: a cutoff on linear modes,

$$P_{\rm hm}(k;M,z) \equiv b_1^{\rm eff} P_m^{\rm NL}(k;z) + b_2(M)\int\frac{d^3q}{(2\pi)^3} P_m^{\rm L}(q)P_m^{\rm L}(|\mathbf{k}-\mathbf{q}|) F_2(\mathbf{q},\mathbf{k}-\mathbf{q}),$$3

a cutoff on the full nonlinear field,

$$P_{\rm hm}(k;M,z) \equiv b_1^{\rm eff} P_m^{\rm NL}(k;z) + b_2(M)\int\frac{d^3q}{(2\pi)^3} P_m^{\rm L}(q)P_m^{\rm L}(|\mathbf{k}-\mathbf{q}|) F_2(\mathbf{q},\mathbf{k}-\mathbf{q}),$$4

and a symmetrized cutoff on the convolution. The reported result is that boost-invariant counterterms depend on the regularization scheme, whereas non-boost-invariant counterterms are scheme independent. The latter are tied to the appearance of

$$P_{\rm hm}(k;M,z) \equiv b_1^{\rm eff} P_m^{\rm NL}(k;z) + b_2(M)\int\frac{d^3q}{(2\pi)^3} P_m^{\rm L}(q)P_m^{\rm L}(|\mathbf{k}-\mathbf{q}|) F_2(\mathbf{q},\mathbf{k}-\mathbf{q}),$$5

which encodes that $$P_{\rm hm}(k;M,z) \equiv b_1^{\rm eff} P_m^{\rm NL}(k;z) + b_2(M)\int\frac{d^3q}{(2\pi)^3} P_m^{\rm L}(q)P_m^{\rm L}(|\mathbf{k}-\mathbf{q}|) F_2(\mathbf{q},\mathbf{k}-\mathbf{q}),$$6 appears at the Lagrangian position $$P_{\rm hm}(k;M,z) \equiv b_1^{\rm eff} P_m^{\rm NL}(k;z) + b_2(M)\int\frac{d^3q}{(2\pi)^3} P_m^{\rm L}(q)P_m^{\rm L}(|\mathbf{k}-\mathbf{q}|) F_2(\mathbf{q},\mathbf{k}-\mathbf{q}),$$7, not the Eulerian position $$P_{\rm hm}(k;M,z) \equiv b_1^{\rm eff} P_m^{\rm NL}(k;z) + b_2(M)\int\frac{d^3q}{(2\pi)^3} P_m^{\rm L}(q)P_m^{\rm L}(|\mathbf{k}-\mathbf{q}|) F_2(\mathbf{q},\mathbf{k}-\mathbf{q}),$$8. Differences between regularizations arise from surface terms such as $$P_{\rm hm}(k;M,z) \equiv b_1^{\rm eff} P_m^{\rm NL}(k;z) + b_2(M)\int\frac{d^3q}{(2\pi)^3} P_m^{\rm L}(q)P_m^{\rm L}(|\mathbf{k}-\mathbf{q}|) F_2(\mathbf{q},\mathbf{k}-\mathbf{q}),$$9 and $\delta N$0, whereas the non-boost-invariant structure associated with the Lagrangian shift is universal (Patrone et al., 2023).

A Wilsonian finite-cutoff formulation makes this dependence dynamical. In the finite-cutoff effective field theory of large-scale structure, the galaxy density is expanded in operators built from the smoothed field, and the bias coefficients run with the cutoff $\delta N$1. The shell integration $\delta N$2 yields

$\delta N$3

with

$\delta N$4

This formulation is used to argue that one should choose $\delta N$5 such that

$\delta N$6

and that the standard renormalized-$\delta N$7-point-function scheme is recovered in the infrared limit even though finite perturbative loop contributions are retained explicitly at finite cutoff (Rubira et al., 2023).

Once stochasticity is included, the effective action contains

$\delta N$8

where $\delta N$9, $b_1^{\rm eff}$0 gives Gaussian stochastic terms, and higher $b_1^{\rm eff}$1 encode higher stochastic moments. The corresponding RG flow closes on the full set of $b_1^{\rm eff}$2, and a single nonlinear bias operator, specifically $b_1^{\rm eff}$3, is stated to generate all stochastic moments through RG evolution. With primordial non-Gaussianity, the operator basis must be enlarged further to include $b_1^{\rm eff}$4-type operators, for example $b_1^{\rm eff}$5 for spin-0, or $b_1^{\rm eff}$6 for spin-2, and a new stochastic contribution appears from the annihilation part of the interaction shell integrals (Rubira et al., 2024, Nikolis et al., 2024).

3. Two-loop renormalization and multi-statistic galaxy bias

A two-loop extension systematizes this program for the full deterministic basis up to fifth order in the density field. The minimal complete leading-gradient basis contains 29 deterministic bias operators up to fifth order, organized through the tensors $b_1^{\rm eff}$7. The two-loop renormalization condition removes both single-hard and double-hard UV regions. For the fifth-order kernels relevant to the two-loop power spectrum,

$b_1^{\rm eff}$8

and the double-hard coefficient takes the universal form

$b_1^{\rm eff}$9

with

$b_1(M)$0

The one-loop RG equation is

$b_1(M)$1

and the two-loop correction is

$b_1(M)$2

Diagonalizing the one-loop RG matrix reveals a positive eigenvalue $b_1(M)$3 in the three-operator sector, identifying a linear combination with enhanced UV sensitivity. The same work explicitly develops the analogy with QFT running and large-log resummation (Bakx et al., 18 Jul 2025).

A subsequent operator-level treatment extends renormalization simultaneously to the two-loop power spectrum and the one-loop bispectrum and trispectrum, while also including gradient corrections to deterministic bias operators at next-to-leading order. The deterministic galaxy overdensity is expanded as

$b_1(M)$4

with leading-gradient operators $b_1(M)$5 and next-to-leading-gradient operators $b_1(M)$6. Renormalized operators are defined by

$b_1(M)$7

and the renormalization conditions are imposed directly on low-momentum correlators and on their quadratic-in-external-momentum pieces. The next-to-leading-gradient sector includes overall-gradient operators, non-overall-gradient operators, and additional building blocks such as

$b_1(M)$8

which generate operators like

$b_1(M)$9

The same framework uses an operator product expansion to renormalize products of two, three, and four operators at coincidence, thereby generating the stochastic counterterms needed for the power spectrum, bispectrum, and trispectrum. The final renormalized correlators take the schematic form

$b_2$0

$b_2$1

$b_2$2

The reported conclusion is that operator-level renormalization makes the power spectrum, bispectrum, and trispectrum mutually consistent and leads to a pronounced scale-dependence of higher-gradient bias coefficients (Bakx et al., 30 Jun 2026).

4. Inflationary renormalization bias

In inflationary perturbation theory, the relevant question is whether renormalization of the primordial power spectrum can induce a significant observable bias. The curvature perturbation is written as $b_2$3, with

$b_2$4

and for slow-roll inflation in the Bunch–Davies vacuum,

$b_2$5

The coincidence limit is UV divergent, so renormalization is required. Using asymptotic regularization rather than adiabatic subtraction, the large-$b_2$6 expansion

$b_2$7

isolates the quadratic and logarithmic UV divergences. In the super-Hubble regime $b_2$8, the renormalized spectrum becomes

$b_2$9

This explicitly separates the physical spectrum, a subtraction sector scaling like $\delta^2$0, and a finite renormalization ambiguity scaling like $\delta^2$1. During inflation,

$\delta^2$2

so for slow roll, $\delta^2$3, the leading renormalization correction is roughly $\delta^2$4. The relative deviation

$\delta^2$5

therefore satisfies

$\delta^2$6

yielding values like $\delta^2$7 for $\delta^2$8 and $\delta^2$9 for $s^2$0 when $s^2$1. Only for $s^2$2–$s^2$3 do observable effects become appreciable. The stated conclusion is that renormalization does not induce a significant, observable bias in the primordial inflationary power spectrum within physically admissible schemes, because the renormalization sector decays while the physical curvature perturbation freezes on super-Hubble scales (Romero et al., 18 Jun 2026).

5. Bias correction by renormalization in embeddings and visual-inertial estimation

In text embedding models, the term refers to a systematic mean bias in the representation space. If $s^2$4 is the embedding of text $s^2$5, the corpus mean vector

$s^2$6

is reported to be significantly nonzero, and each embedding is modeled as

$s^2$7

where $s^2$8 is the semantic signal and $s^2$9 is a nearly constant bias vector shared across sentences. The proposed post-processing step, called Renormalization, has two variants. Variant R1 subtracts the mean directly,

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$00

while Variant R2 subtracts only the projection onto the bias direction,

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$01

Under a first-order noise analysis with $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$02, the residual for R1 is

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$03

whereas for R2

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$04

The prediction is therefore that R2 should outperform R1. On MMTEB, using 100,000 randomly sampled sentences from 20220301.en to estimate $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$05, experiments across 38 publicly available embedding models report aggregate gains of $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$06 on retrieval tasks, $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$07 on classification tasks, and $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$08 on other tasks for R2, with R2 consistently outperforming R1. The reported correlation between improvement and $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$09 is positive, including a Pearson correlation around $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$10 in one figure (Ren et al., 14 Nov 2025).

In rolling-shutter visual-inertial odometry, renormalization denotes a statistical procedure derived from Kanatani’s scheme for reducing the inherent bias of a reduced linear least-squares initializer. After eliminating depth variables by Schur complement, the problem is reduced to

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$11

Ordinary least squares minimizes an algebraic residual and is described as statistically biased because the reduced rows depend nonlinearly on noisy image correspondences. Renormalization replaces the basic eigenproblem with the generalized eigenvalue problem

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$12

where $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$13 is built from first-order noise propagation through the image-to-row mapping. The iterative algorithm initializes $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$14, forms

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$15

and updates the weight matrices from the current solution. The paper reports improvements over least squares of roughly 9% to 35% in velocity error and 5% to 15% in gravity angle error depending on sequence, summarized as around 20% improvement for initial velocity and around 8% improvement for gravity. Renormalization converges in 2–5 iterations, solves only a $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$16 generalized eigenproblem at each iteration, and takes about 70% of BA runtime in the reported Matlab implementation (Micusik et al., 2020).

6. Deliberate bias potentials and voltage-bias-dependent renormalization

In variational Monte Carlo renormalization-group methods, “bias” refers neither to an observable tracer bias nor to a statistical estimation bias. It is an intentionally constructed auxiliary potential $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$17 acting on coarse-grained variables. The method is built on the convex functional

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$18

whose minimizer is unique up to a constant and satisfies

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$19

If the target distribution $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$20 is chosen to be constant, then

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$21

so the bias potential cancels the effective interactions among block spins and renders the coarse-grained variables effectively uncorrelated. The renormalized couplings are then read off directly as

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$22

from the expansion $$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$23. In this usage, the phrase “biased ensemble” means a controlled variational reweighting that overcomes critical slowing down, not an unwanted distortion (Wu et al., 2017).

A different contrast appears in nanoscale junctions, where “bias” denotes the applied voltage across the junction rather than an inferential or cosmological bias parameter. There, renormalization concerns the frequency, damping, and heating of vibrational modes coupled to charge transport. The dressed vibrational propagator is

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$24

and the renormalized vibrational frequency is determined by

$$\sigma^2\equiv \int \frac{d^3q}{(2\pi)^3}P_m^{\rm L}(q),$$25

The paper reports a strong bias dependence of frequency renormalization and vibrational damping in junctions with intermediate lead coupling, and argues that the Raman shifts and linewidths observed in an OPV3 junction may be explained by a combination of dynamic carrier screening and molecular charging. This usage is conceptually separate from “renormalization bias” in cosmology or statistics, but it shows that the juxtaposition of “bias” and “renormalization” can also arise because renormalized quantities depend on an external voltage bias (Kaasbjerg et al., 2013).

In this literature, “renormalization bias” is therefore not a single invariant concept. It can mean renormalized tracer couplings in cosmology, the absence of an observable renormalization-induced distortion in inflation, a training-free removal of geometric mean bias in embeddings, a statistical bias-reduction scheme in visual-inertial initialization, or a deliberate variational bias potential in Monte Carlo RG. The common element is that renormalization changes how a nominal bias parameter is defined, measured, or removed, but the technical content depends entirely on the underlying theory and observable.