The Quantum Vacuum Self-Consistency Principle: Emergent Dynamics of Spacetime and the Standard Model

Abstract

The principle of self-consistency of the quantum vacuum postulates that the classical backgrounds we observe—spacetime geometry, gauge fields, and the Higgs condensate—are macroscopic order parameters of a single quantum state whose existence is sustained by the vacuum expectation values of all quantum fluctuations living on it. Building on this postulate, a background-field, heat-kernel based derivation is developed that yields the coupled low-energy effective field equations for the metric, gauge fields, and the Higgs field as vacuum “equations of state.” The resulting framework rigorously recovers the Einstein, Yang-Mills, and Higgs equations, augmented by the higher-derivative operators required by quantum consistency, with renormalized couplings determined by the content of quantum fields. A Renormalization Group (RG) structure is obtained in which all couplings—including Newton’s constant, the cosmological constant, the coefficients of $R^{2}$ and $C_{\mu\nu\rho\sigma}^{2}$, gauge couplings, Yukawas, and the Higgs quartic—run coherently with a single effective action that respects background diffeomorphism and gauge invariance. A simplified solvable model is analyzed: an $O(N)$ scalar sector with nonminimal coupling on a constant-curvature background, for which the one-loop effective potential including the leading curvature corrections is computed and the vacuum gap equations are solved explicitly. Phenomenological consequences follow. First, a robust prediction: the anomaly- and loop-induced $R^{2}$ operator generically drives Starobinsky-type inflation, with $n_{s}\simeq 1-2/N_{e}$ and $r\simeq 12/N_{e}^{2}$, compatible with Planck data; the required coefficient corresponds to a scalaron mass $M\simeq(1.3\pm 0.1)\times 10^{-5}M_{\rm Pl}$. Second, the universally calculable quantum correction to Newton’s potential and the Yukawa tails from the massive spin-0 and spin-2 modes are quantified and shown to satisfy laboratory bounds. Third, constraints from GW170817 enforce luminal gravitational wave speed for the massless graviton in this framework, while higher-derivative effects remain suppressed at LIGO/Virgo/KAGRA frequencies. The theory is predictive in its inflationary sector and in its universal low-energy corrections to gravity, while it remains honest about open issues: the smallness of the observed cosmological constant, the nonperturbative completion of the higher-derivative sector, and the determination of threshold-matched couplings beyond one loop. The structure is sufficiently complete to be testable across cosmology, astrophysics, and precision gravity, and it reduces to General Relativity and the Standard Model at accessible scales with controlled corrections.

I Introduction

The Standard Model (SM) of particle physics and General Relativity (GR) are profoundly successful descriptions of nature [Glashow1961, Weinberg1967, Salam1968, Misner1973]. Yet their unification remains elusive. Quantizing GR perturbatively leads to nonrenormalizable divergences [tHooft1974], which signal the need for a deeper organizing principle. A natural path, inspired by Sakharov’s induced gravity [Sakharov1968], is that the classical backgrounds we observe are emergent, macroscopic order parameters of an underlying quantum vacuum. In this view, geometry, gauge fields, and the Higgs condensate are not separate axioms; rather, they are coherent manifestations of one vacuum that sustains itself through self-consistency.

The central postulate is that the vacuum state $|\Omega\rangle$ determined by classical backgrounds must, through quantum fluctuations, generate precisely the sources that maintain those same backgrounds. This self-consistency is enforced by the quantum effective action computed in the presence of background fields. The background field method and the heat-kernel expansion provide the technical backbone [DeWitt2003, Vassilevich2003, Barvinsky1985, BirrellDavies1982, ParkerToms2009, Buchbinder1992, Avramidi2000]. At one loop, the Seeley–DeWitt coefficient $a_{2}$ fixes the divergences and therefore the RG running of all couplings that can appear in a local effective action.

This work develops that program into a predictive framework. The self-consistency postulate is kept intact and is expressed via vacuum gap equations that are precisely the stationarity conditions of the renormalized 1PI effective action with respect to the background fields. We compute and reorganize the one-loop effective action to derive: (i) the renormalized Einstein equations plus controlled higher-derivative corrections, (ii) the renormalized Yang–Mills equations, and (iii) the renormalized Higgs equation and potential. We also construct the coupled RG equations for gravitational and SM couplings, which makes the framework calculable.

A simplified, analytically tractable model is solved explicitly to exhibit the mechanism: an $O(N)$ scalar with nonminimal coupling on a constant-curvature background. The one-loop effective potential including curvature-dependent terms is computed and the gap equations are solved for the vacuum expectation value and curvature. The analysis isolates which combinations of couplings are fixed by self-consistency and which must be matched to data.

Finally, the framework is confronted with data. The loop- and anomaly-induced $R^{2}$ term drives Starobinsky inflation [Starobinsky1980, Whitt1984], yielding benchmark predictions for the spectral tilt and tensor-to-scalar ratio consistent with Planck [Planck2018]. Quantum corrections to Newton’s potential are universal and match the effective field theory result [Donoghue1994]. Short-range Yukawa corrections from the massive scalar (the scalaron) and the heavy spin-2 mode are shown to obey submillimeter bounds [Kapner2007]. Gravitational wave propagation is luminal in the infrared, in agreement with GW170817 [GW170817].

The approach is conservative in spirit: it embraces effective field theory, honors all symmetries via the background field method, and leans on mathematically controlled tools. At the same time, it reframes the unification problem in terms of a single organizing idea: vacuum self-consistency.

II Unified framework and the self-consistency postulate

All fields are decomposed into classical backgrounds and quantum fluctuations,

$$\hat{g}_{\mu\nu}(x) =\bar{g}_{\mu\nu}(x)+\hat{h}_{\mu\nu}(x),$$ (1)

$$\hat{A}_{\mu}^{a}(x) =\bar{A}_{\mu}^{a}(x)+\hat{a}_{\mu}^{a}(x),$$ (2)

$$\hat{\Phi}(x) =\frac{1}{\sqrt{2}}\bigl(v_{H}(x)+\hat{\sigma}(x)\bigr),$$ (3)

$$\hat{\Psi}(x) =\hat{\psi}(x).$$ (4)

Fermions are purely quantum in the vacuum since they do not condense. Backgrounds $\{\bar{g},\bar{A},v_{H}\}$ define a vacuum $|\Omega\rangle=|0_{\bar{g},\bar{A},v_{H}}\rangle$.

The unified self-consistency postulate is that the backgrounds are sustained by the renormalized vacuum expectation values (VEVs) of their source operators in $|\Omega\rangle$. Equivalently, in the presence of external sources $J$, the renormalized 1PI effective action $\Gamma$ obeys the background-field stationarity conditions,

$$\frac{2}{\sqrt{-\bar{g}}}\frac{\delta\Gamma}{\delta\bar{g}^{\mu\nu}} =T_{\mu\nu}^{\rm ext},$$ (5)

$$\frac{1}{\sqrt{-\bar{g}}}\frac{\delta\Gamma}{\delta\bar{A}_{\mu}^{a}} =J^{a\,\mu}_{\rm ext},$$ (6)

$$\frac{1}{\sqrt{-\bar{g}}}\frac{\delta\Gamma}{\delta v_{H}} =J_{\sigma}^{\rm ext}.$$ (7)

In the absence of external sources, these become vacuum gap equations. Writing the bare (classical) actions with bare parameters and integrating out fluctuations generates the quantum terms that renormalize both parameters and the structure of the equations. The background field method ensures background diffeomorphism and gauge invariance, so that Ward identities are preserved [Abbott1981, Buchbinder1992].

To one loop, the quantum effective action is

$$\Gamma[\bar{g},\bar{A},v_{H}] =S_{\rm bare}[\bar{g},\bar{A},v_{H}]$$

$$\quad+\frac{i}{2}\mathrm{Tr}\ln\mathcal{D}_{\rm bos}-i\,\mathrm{Tr}\ln\mathcal{D}_{\rm ferm}$$

$$\quad-i\,\mathrm{Tr}\ln\mathcal{D}_{\rm ghost}+\cdots,$$ (8)

with minimal second-order operators of the form $\mathcal{D}=-\bar{\nabla}^{2}+\mathcal{P}$ acting on appropriate bundles and $\mathcal{P}$ containing masses, curvature couplings, and background field strengths. Divergences are governed by the heat-kernel coefficient $a_{2}$ [DeWitt2003, Vassilevich2003, Barvinsky1985, Avramidi2000]. For each species $j$, one finds a local divergent density

$$\mathcal{L}^{(1)}_{\rm div,j} =\frac{1}{16\pi^{2}\epsilon}\sqrt{-\bar{g}}\;\mathrm{tr}\,a_{2,j}(x),$$ (9)

with $\epsilon=4-d$ in dimensional regularization and the trace over internal indices. Summing over all fluctuating fields with appropriate statistics yields

$$\Gamma_{\rm div}^{(1)} =\frac{1}{16\pi^{2}\epsilon}\int d^{4}x\sqrt{-\bar{g}}\Bigl[C_{\Lambda}+C_{R}\bar{R}$$

$$\quad+C_{C^{2}}\,C_{\mu\nu\rho\sigma}C^{\mu\nu\rho\sigma}+C_{E}\,E_{4}$$

$$\quad+\sum_{i}C_{F_{i}}\,\mathrm{tr}\bigl(\bar{F}_{\mu\nu}^{(i)}\bar{F}^{(i)\mu\nu}\bigr)$$

$$\quad+C_{RH}\,\bar{R}\,v_{H}^{2}+C_{H}\,v_{H}^{4}+\cdots\Bigr],$$ (10)

where $E_{4}=R_{\mu\nu\rho\sigma}^{2}-4R_{\mu\nu}^{2}+R^{2}$ is the Euler density and $C_{\mu\nu\rho\sigma}$ is the Weyl tensor. The coefficients $C_{X}$ are calculable sums of species-dependent numbers and masses, and depend on the nonminimal couplings of scalars. The ellipsis denotes higher-derivative scalar terms and total derivatives such as $\Box R$ that can be discarded or absorbed into local counterterms.

Finite parts, including nonlocal form factors that resum infrared and threshold effects, are encoded in the full $\Gamma$ [BarvinskyVilkovisky1985b, BarvinskyVilkovisky1990] and are essential for precision predictions, but for many purposes the local part suffices.

III Stationarity conditions and emergent field equations

The gap equations (5)–(7) are the equations of state of the vacuum. They are obtained by varying the renormalized effective action,

$$\Gamma =\int d^{4}x\sqrt{-\bar{g}}\Biggl[\frac{Z_{R}}{16\pi G}\bar{R}-\frac{Z_{\Lambda}}{8\pi G}\Lambda$$

$$\quad+\alpha_{C}C_{\mu\nu\rho\sigma}C^{\mu\nu\rho\sigma}+\alpha_{R}\bar{R}^{2}$$

$$\quad-\sum_{i}\frac{Z_{F_{i}}}{4g_{i}^{2}}\mathrm{tr}\bigl(\bar{F}_{\mu\nu}^{(i)}\bar{F}^{(i)\mu\nu}\bigr)$$

$$\quad+Z_{H}\frac{1}{2}\bar{g}^{\mu\nu}\partial_{\mu}v_{H}\partial_{\nu}v_{H}-V_{\rm eff}(v_{H},\bar{R})$$

$$\quad+\cdots\Biggr]+\Gamma_{\rm nl},$$ (11)

where $Z$-factors include finite renormalizations, $\alpha_{C}$ and $\alpha_{R}$ are the renormalized dimensionless couplings of the curvature-squared operators, $V_{\rm eff}$ is the renormalized effective potential for the Higgs condensate including curvature dependence, and $\Gamma_{\rm nl}$ collects nonlocal form factors such as $R\ln(-\bar{\Box}/\mu^{2})R$ that are known from the covariant perturbation theory of Barvinsky and Vilkovisky [BarvinskyVilkovisky1990]. Dots denote higher-order and higher-derivative terms suppressed at low energies.

Varying (11) with respect to $\bar{g}^{\mu\nu}$ yields

$$Z_{R}\,G_{\mu\nu}+Z_{\Lambda}\,\Lambda\,\bar{g}_{\mu\nu}$$

$$\quad+16\pi G\Bigl[\alpha_{C}\,B_{\mu\nu}+\alpha_{R}\,H^{(R^{2})}_{\mu\nu}\Bigr]$$

$$=8\pi G\Bigl(T_{\mu\nu}^{\rm ext}+T_{\mu\nu}^{(H)}+T_{\mu\nu}^{\rm nl}\Bigr),$$ (12)

where $B_{\mu\nu}$ is the Bach tensor,

$$B_{\mu\nu} =\nabla^{\rho}\nabla^{\sigma}C_{\mu\rho\nu\sigma}+\tfrac{1}{2}R^{\rho\sigma}C_{\mu\rho\nu\sigma},$$ (13)

and $H^{(R^{2})}_{\mu\nu}$ is the metric variation of $R^{2}$,

$$H^{(R^{2})}_{\mu\nu} =2RR_{\mu\nu}-\tfrac{1}{2}\bar{g}_{\mu\nu}R^{2}-2\nabla_{\mu}\nabla_{\nu}R+2\bar{g}_{\mu\nu}\Box R.$$ (14)

The tensor $T_{\mu\nu}^{(H)}$ contains the stress tensor of $v_{H}$,

$$T_{\mu\nu}^{(H)} =Z_{H}\partial_{\mu}v_{H}\partial_{\nu}v_{H}-\tfrac{1}{2}Z_{H}\bar{g}_{\mu\nu}(\partial v_{H})^{2}+\bar{g}_{\mu\nu}V_{\rm eff},$$ (15)

and $T_{\mu\nu}^{\rm nl}$ encodes finite nonlocal contributions. In the long-wavelength, weak-curvature regime, the higher-derivative and nonlocal terms are suppressed, and (12) reduces to the Einstein equations with renormalized constants.

Variation with respect to the gauge backgrounds yields

$$\Bigl(\frac{Z_{F_{i}}}{g_{i}^{2}}\Bigr)\,\bar{D}_{\nu}\bar{F}^{(i)\nu\mu} =J^{(i)\mu}_{\rm ext}+J^{(i)\mu}_{\rm nl},$$ (16)

with $J^{(i)\mu}_{\rm nl}$ the finite nonlocal current from integrating out matter and gravity.

Variation with respect to $v_{H}$ gives the renormalized Higgs equation,

$$Z_{H}\,\Box_{\bar{g}}v_{H}-\frac{\partial V_{\rm eff}(v_{H},\bar{R})}{\partial v_{H}} =J_{\sigma}^{\rm ext},$$ (17)

where $V_{\rm eff}(v_{H},\bar{R})$ includes quantum and curvature corrections. In the vacuum, $J_{\rm ext}=0$, and (12)–(17) become the coupled gap equations that determine the macroscopic order parameters.

IV Heat-kernel coefficients and renormalization

The coefficients $C_{X}$ in (10) are determined by the $a_{2}$ coefficient of the heat-kernel expansion. For a Laplace-type operator $\mathcal{D}=-\bar{\nabla}^{2}+\mathcal{P}$ acting on a vector bundle with curvature $\mathcal{F}_{\mu\nu}=[\bar{\nabla}_{\mu},\bar{\nabla}_{\nu}]$, one has [Vassilevich2003, Avramidi2000]

$$a_{2} =\tfrac{1}{180}\Bigl(R_{\mu\nu\rho\sigma}^{2}-R_{\mu\nu}^{2}+\Box R\Bigr)\mathbb{I}+\tfrac{1}{2}\mathcal{P}^{2}$$

$$\quad-\tfrac{1}{6}R\,\mathcal{P}+\tfrac{1}{12}\mathcal{F}_{\mu\nu}^{2}.$$ (18)

For a real scalar with nonminimal coupling $\xi$, $\mathcal{P}=m^{2}+\xi R+\cdots$. For a Dirac fermion, $\mathcal{P}=m^{2}+\tfrac{1}{4}R+\cdots$ and one must include a minus sign from statistics and an overall factor from spinor components. For gauge vectors in background gauge, the operator acts on vectors and ghosts contribute with their own operators. Combining all species, one obtains the standard trace-anomaly coefficients for the $C^{2}$ and $E_{4}$ structures [Duff1994, ChristensenDuff1978]. For $N_{s}$ real scalars, $N_{f}$ Dirac fermions, and $N_{v}$ gauge vectors,

$$C_{C^{2}} =\frac{1}{120}N_{s}+\frac{1}{20}N_{f}+\frac{1}{10}N_{v},$$ (19)

$$C_{E} =\frac{1}{360}N_{s}+\frac{11}{360}N_{f}+\frac{31}{180}N_{v}.$$ (20)

These appear as the coefficients of the $1/\epsilon$ divergences in front of $C^{2}$ and $E_{4}$ in (10). Mass-dependent terms contribute to $C_{\Lambda}$ and $C_{R}$,

$$C_{\Lambda} =\sum_{j}(-1)^{F_{j}}\,n_{j}\,m_{j}^{4},$$ (21)

$$C_{R} =\sum_{s}\Bigl(\xi_{s}-\tfrac{1}{6}\Bigr)m_{s}^{2}+\sum_{f}\Bigl(-\tfrac{1}{12}\Bigr)m_{f}^{2}+\sum_{v}\Bigl(\tfrac{1}{6}\Bigr)m_{v}^{2},$$ (22)

where $n_{j}$ counts internal degrees of freedom and $F_{j}$ is 0 (1) for bosons (fermions). In dimensional regularization, strictly massless fields do not contribute to $C_{R}$ at one loop, while $C_{C^{2}}$ and $C_{E}$ are mass-independent. The gauge-field divergences yield

$$C_{F_{i}} =\kappa_{i}-\sum_{j}\ell_{ij},$$ (23)

with $\kappa_{i}$ the adjoint contribution and $\ell_{ij}$ matter contributions in representation $j$ of gauge group $i$.

Renormalization proceeds by absorbing (10) into counterterms, defining the renormalized couplings at scale $\mu$ in minimal subtraction. The resulting RG equations take the schematic form

$$\mu\frac{d}{d\mu}\Bigl(\frac{Z_{R}}{16\pi G}\Bigr) =-\frac{C_{R}}{8\pi^{2}},$$ (24)

$$\mu\frac{d}{d\mu}\Bigl(\frac{Z_{\Lambda}\Lambda}{8\pi G}\Bigr) =\frac{C_{\Lambda}}{16\pi^{2}},$$ (25)

$$\mu\frac{d\alpha_{C}}{d\mu} =\frac{C_{C^{2}}}{8\pi^{2}},\qquad\mu\frac{d\alpha_{R}}{d\mu}=\frac{C_{R^{2}}}{8\pi^{2}},$$ (26)

$$\mu\frac{d}{d\mu}\Bigl(\frac{Z_{F_{i}}}{g_{i}^{2}}\Bigr) =-\frac{C_{F_{i}}}{8\pi^{2}},$$ (27)

with $C_{R^{2}}$ a linear combination of $C_{E}$ and species-dependent terms that accompany $R^{2}$. The precise mapping between $\{C_{E},C_{C^{2}}\}$ and $\{\beta_{\alpha_{R}},\beta_{\alpha_{C}}\}$ depends on the chosen basis, see [Buchbinder1992, Duff1994]. In the matter sector, the SM one-loop $\beta$ functions are standard [MachacekVaughnI, MachacekVaughnII, Buttazzo2013].

Nonlocal terms in $\Gamma$ carry logarithms of the form $\ln(-\bar{\Box}/\mu^{2})$ and govern both decoupling and threshold behavior [BarvinskyVilkovisky1990, Donoghue1994]. At scales $\mu$ below a particle’s mass, matching across thresholds must be performed to ensure the decoupling of heavy species in gauge and Yukawa sectors [AppelquistCarazzone1975]. In curved space this procedure is technically involved but conceptually the same, and the covariant perturbation theory provides the right tool.

V A solvable model: $O(N)$ scalar on constant curvature

To display the mechanism in a fully analytic setting, consider $N$ real scalars $\phi^{I}$ with $O(N)$ symmetry and classical potential $V(\phi)=\tfrac{\lambda}{4}(\phi^{2}-v_{0}^{2})^{2}$, nonminimally coupled with coupling $\xi$ to the curvature, on a background of constant scalar curvature $\bar{R}$. Decompose $\phi^{I}=(v+\sigma,\pi^{a})$ with $a=1,\dots,N-1$, so that $v=\langle\phi\rangle$ plays the role of $v_{H}$. The classical Euclidean action relevant for the effective potential is

$$S_{\rm cl} =\int d^{4}x\sqrt{\bar{g}}\Bigl[\frac{M_{\rm Pl}^{2}}{2}\bar{R}+\alpha_{R}\bar{R}^{2}-\frac{1}{2}\xi\bar{R}\,\phi^{2}$$

$$\quad+\frac{1}{2}(\bar{\nabla}\phi)^{2}+\frac{\lambda}{4}(\phi^{2}-v_{0}^{2})^{2}+\cdots\Bigr].$$ (28)

In the constant background, gradients vanish and the tree-level potential density reads

$$U_{\rm tree}(v,\bar{R}) =-\frac{M_{\rm Pl}^{2}}{2}\bar{R}-\alpha_{R}\bar{R}^{2}+\frac{1}{2}\xi\bar{R}\,v^{2}$$

$$\quad+\frac{\lambda}{4}(v^{2}-v_{0}^{2})^{2}+\rho_{\Lambda},$$ (29)

where $\rho_{\Lambda}$ collects the cosmological constant term and constant counterterms.

The one-loop effective potential including the leading curvature dependence is obtained by summing the fluctuation determinants of $\sigma$ and $\pi^{a}$. Their field-dependent masses in the constant background are

$$m_{\pi}^{2}(v,\bar{R}) =\lambda(v^{2}-v_{0}^{2})+\xi\bar{R},$$ (30)

$$m_{\sigma}^{2}(v,\bar{R}) =\lambda(3v^{2}-v_{0}^{2})+\xi\bar{R}.$$ (31)

At one loop in $\overline{\rm MS}$,

$$U_{\rm 1\mbox{-}loop}(v,\bar{R}) =\sum_{i=\pi,\sigma}\frac{n_{i}}{64\pi^{2}}m_{i}^{4}(v,\bar{R})$$

$$\quad\times\Bigl[\ln\frac{m_{i}^{2}(v,\bar{R})}{\mu^{2}}-c_{i}\Bigr]$$

$$\quad+\mathcal{O}\bigl(\bar{R}^{2}\ln(\mu)\bigr),$$ (32)

with $n_{\pi}=N-1$, $n_{\sigma}=1$, and $c_{i}=\tfrac{3}{2}$ for scalars in the $\overline{\rm MS}$ scheme. The $\mathcal{O}(\bar{R}^{2})$ terms in $U_{\rm 1\mbox{-}loop}$ can be absorbed into a running of $\alpha_{R}(\mu)$ plus total derivatives. The full effective potential is

$$U_{\rm eff}(v,\bar{R}) =U_{\rm tree}(v,\bar{R})+U_{\rm 1\mbox{-}loop}(v,\bar{R}).$$ (33)

The vacuum gap equations are the stationarity conditions with respect to $v$ and $\bar{R}$,

$$\frac{\partial U_{\rm eff}}{\partial v} =0,$$ (34)

$$\frac{\partial U_{\rm eff}}{\partial\bar{R}} =0.$$ (35)

Explicitly,

$$0 =\lambda(v^{2}-v_{0}^{2})v+\xi\bar{R}\,v$$

$$\quad+\sum_{i}\frac{n_{i}}{32\pi^{2}}m_{i}^{2}(v,\bar{R})\,\frac{\partial m_{i}^{2}}{\partial v}\Bigl[\ln\frac{m_{i}^{2}}{\mu^{2}}-1\Bigr],$$ (36)

$$0 =-\frac{M_{\rm Pl}^{2}}{2}-2\alpha_{R}\bar{R}+\frac{\xi}{2}v^{2}$$

$$\quad+\sum_{i}\frac{n_{i}}{64\pi^{2}}\Biggl\{2m_{i}^{2}\,\frac{\partial m_{i}^{2}}{\partial\bar{R}}\Bigl[\ln\frac{m_{i}^{2}}{\mu^{2}}-1\Bigr]\Biggr\}.$$ (37)

To see the origin of the numerical factor in (36), note that for a single mode the one-loop term is $\Delta U=\frac{1}{64\pi^{2}}m^{4}[\ln(m^{2}/\mu^{2})-\tfrac{3}{2}]$, hence

$$\frac{\partial\Delta U}{\partial v} =\frac{1}{64\pi^{2}}\,2m^{2}\frac{\partial m^{2}}{\partial v}\Bigl[\ln\frac{m^{2}}{\mu^{2}}-1\Bigr],$$ (38)

which produces the overall $1/(32\pi^{2})$ factor after summing degeneracies $n_{i}$.

Using (30)–(31), the derivatives are $\partial m_{i}^{2}/\partial v=2\lambda c_{i}v$ with $c_{\pi}=1$, $c_{\sigma}=3$, and $\partial m_{i}^{2}/\partial\bar{R}=\xi$. Solving (36)–(37) for $(v,\bar{R})$ at given $(\lambda,\xi,\alpha_{R},\mu)$ demonstrates the self-consistent determination of the order parameters.

Two limiting cases are instructive.

First, if $\bar{R}=0$ and quantum corrections are small at $\mu\sim v$, then (36) gives $v\simeq v_{0}$ as usual, and (37) fixes the renormalized cosmological constant counterterm via $\rho_{\Lambda}$ so that Minkowski is a solution. This is the standard tuning familiar from semiclassical gravity.

Second, if $\alpha_{R}$ is large and positive, a de Sitter solution arises with $\bar{R}>0$ dominated by the $R^{2}$ term. At tree level, (37) yields $\bar{R}\simeq\tfrac{1}{4\alpha_{R}}\bigl(\xi v^{2}-M_{\rm Pl}^{2}\bigr)$, and (36) reduces to $(\lambda v^{2}+\xi\bar{R}-\lambda v_{0}^{2})v\simeq 0$. For $\xi>0$ and $v\neq 0$, these equations define a consistent de Sitter branch, which is lifted by the loop terms in the expected way. This branch underlies the $R^{2}$ inflation scenario discussed below.

This solvable model shows explicitly how the self-consistency postulate supplies coupled algebraic equations that fix the background curvature and the condensate in terms of renormalized couplings. The same logic extends to the full SM plus gravity, with the technical difference that gauge bosons and fermions contribute additional terms and thresholds.

VI Phenomenology and predictions

VI.1 Inflation from the induced $R^{2}$ term

The universal generation of curvature-squared operators is unavoidable at one loop. The $R^{2}$ term is particularly important because it is ghost-free in the scalar sector and, when dominant, is dynamically equivalent to GR plus a scalar field (the scalaron) [Whitt1984, Stelle1977]. In the Einstein frame, the scalaron potential is

$$V(\varphi) =\frac{3}{4}M^{2}M_{\rm Pl}^{2}\Bigl(1-e^{-\sqrt{\frac{2}{3}}\frac{\varphi}{M_{\rm Pl}}}\Bigr)^{2},$$ (39)

and the mapping from the Jordan-frame action $\int\sqrt{-g}\,\bigl[\tfrac{M_{\rm Pl}^{2}}{2}R+\alpha_{R}R^{2}\bigr]$ to (39) gives

$$\alpha_{R} =\frac{M_{\rm Pl}^{2}}{12M^{2}}.$$ (40)

For slow roll with $N_{e}$ e-folds,

$$n_{s} \simeq 1-\frac{2}{N_{e}},\qquad r\simeq\frac{12}{N_{e}^{2}}.$$ (41)

Taking $N_{e}=50$–$60$ yields $n_{s}\simeq 0.962$–$0.967$ and $r\simeq 0.003$–$0.004$, consistent with Planck [Planck2018]. The scalar amplitude fixes $M\simeq(1.3\pm 0.1)\times 10^{-5}M_{\rm Pl}$, hence

$$\alpha_{R} \simeq\frac{M_{\rm Pl}^{2}}{12M^{2}}\simeq(4.5\text{--}5.5)\times 10^{8}.$$ (42)

This is a concrete, quantitative prediction: the self-consistent vacuum generates an $R^{2}$ term whose coefficient is large and positive in the inflationary regime. In the late universe, the $R^{2}$ term is inert and safely compatible with local tests.

VI.2 Short-distance gravity and laboratory bounds

Quadratic gravity around flat space contains, besides the massless graviton, a massive scalar (the scalaron) and a massive spin-2 mode from the $C^{2}$ term [Stelle1977]. The static potential between two point masses $m_{1}$ and $m_{2}$ at separation $r$ is

$$V(r) =-\frac{G_{N}m_{1}m_{2}}{r}\Biggl[1+\frac{1}{3}e^{-Mr}-\frac{4}{3}e^{-m_{2}r}$$

$$\quad+\frac{41}{10\pi}\frac{G_{N}\hbar}{r^{2}c^{3}}+\cdots\Biggr],$$ (43)

where $M^{2}=M_{\rm Pl}^{2}/(12\alpha_{R})$ and $m_{2}^{2}\simeq M_{\rm Pl}^{2}/(2\alpha_{C})$ are the scalar and spin-2 masses, respectively, and the last term is the universal quantum correction [Donoghue1994]. Submillimeter tests constrain extra Yukawa forces with strength $\mathcal{O}(1)$ to have range $\lambda\lesssim 0.1$ mm, implying $M,m_{2}\gtrsim 2\times 10^{-3}$ eV [Kapner2007]. The inflationary value (42) corresponds to $M\sim 10^{13}$ GeV, and a loop-induced $\alpha_{C}$ of order $10^{-2}$–$10^{0}$ yields $m_{2}\sim M_{\rm Pl}/\sqrt{\alpha_{C}}\gtrsim 10^{18}$ GeV, both well above laboratory reach. The sign of the $m_{2}$ residue is negative (a ghost) at the perturbative level [Stelle1977]. As discussed below, this is acceptable within an effective description valid below $m_{2}$, and several proposals exist for its nonperturbative treatment [AnselmiPiva2018, Percacci2017].

VI.3 Gravitational waves

In the infrared limit where higher derivatives are negligible, gravitational waves propagate at the speed of light, $c_{T}=1$, consistent with GW170817 and its electromagnetic counterpart [GW170817]. The massive scalar and spin-2 modes would be excited only at frequencies $\omega\gtrsim M,m_{2}$, far above LIGO/Virgo/KAGRA bands. Nonlocal form factors induce a mild running of $G_{N}$ with momentum but do not modify the luminal propagation of the massless graviton in the observed regime [Donoghue1994, BarvinskyVilkovisky1990].

VI.4 Running couplings and vacuum stability

The RG of the SM couplings is standard and can be combined with the gravitational running in (24)–(27). For the SM matter content at high scales, the trace-anomaly coefficients are

$$c_{\rm SM} =\frac{1}{120}N_{s}+\frac{1}{20}N_{f}+\frac{1}{10}N_{v}$$

$$=\frac{1}{120}\times 4+\frac{1}{20}\times 22.5+\frac{1}{10}\times 12$$

$$\simeq 2.36,$$ (44)

$$a_{\rm SM} =\frac{1}{360}N_{s}+\frac{11}{360}N_{f}+\frac{31}{180}N_{v}$$

$$=\frac{4}{360}+\frac{11}{360}\times 22.5+\frac{31}{180}\times 12$$

$$\simeq 2.77,$$ (45)

with $N_{s}=4$ real scalars, $N_{f}=22.5$ Dirac fermions (equivalently 45 Weyl), and $N_{v}=12$ vectors. These numbers govern the one-loop running of the gravitational higher-derivative couplings in minimal subtraction [Duff1994, Buchbinder1992].

The Higgs quartic coupling $\lambda$ runs to small or even negative values at $10^{9}$–$10^{11}$ GeV depending on inputs [Buttazzo2013]. In this framework, curvature and nonminimal coupling $\xi$ feed back into the effective potential and can stabilize the vacuum during the inflationary epoch. In the late universe, where $\bar{R}\to 0$, the standard SM running applies. A full two-loop study with gravitational corrections will sharpen these statements; the one-loop structure already indicates how the feedback operates.

VII Predictivity and renormalization conditions

Predictivity arises in several ways:

• The $R^{2}$-driven inflationary sector is sharply predictive. Once $\alpha_{R}$ is matched to $A_{s}$, $n_{s}$ and $r$ follow from (41) independently of microphysics [Starobinsky1980, Whitt1984, Planck2018].

• The universal quantum correction to Newton’s potential in (43) has a fixed coefficient [Donoghue1994].

• The strengths and ranges of Yukawa corrections to gravity are set by $\alpha_{R}$ and $\alpha_{C}$, which are constrained by cosmology and laboratory experiments [Kapner2007].

Within the full coupled system, renormalization conditions encode the choice of vacuum. A natural set is:

$$\text{(i)}\quad\frac{\partial U_{\rm eff}}{\partial v}\Big|_{v=v_{\star},\bar{R}=\bar{R}_{\star}}=0,\quad\text{(ii)}\quad\frac{\partial U_{\rm eff}}{\partial\bar{R}}\Big|_{v=v_{\star},\bar{R}=\bar{R}_{\star}}=0,$$

$$\text{(iii)}\quad U_{\rm eff}(v_{\star},\bar{R}_{\star})=\rho_{\rm vac}^{\rm obs}.$$ (46)

Condition (iii) fixes the finite part of the cosmological constant counterterm. In practice, $\rho_{\rm vac}^{\rm obs}$ is tiny compared to any particle physics scale, which requires a tuning familiar from semiclassical gravity. The first two conditions, however, are dynamical gap equations that constrain combinations of couplings and masses.

In applications, one matches the renormalized couplings to measured low-energy observables and then integrates the RG equations upward, solving the gap equations at the desired scale for the background. Conversely, given a high-scale boundary such as inflation, one runs down and predicts low-energy corrections. This is the standard effective field theory strategy implemented with the self-consistency postulate.

VIII Limitations and outlook

The framework developed here rests on a controlled one-loop calculation with the background field method and the heat-kernel expansion. Several limitations and tasks lie ahead.

First, while the one-loop structure is universal and robust, higher-loop and nonperturbative effects could modify the running of gravitational couplings. Functional RG studies suggest nontrivial fixed points for gravity [Weinberg1979, Reuter1998, Percacci2017, Codello2008], which, if realized, would render the entire system asymptotically safe and predictive at all scales. A careful match between the one-loop effective action and functional RG fixed points is an important next step.

Second, quadratic gravity is perturbatively renormalizable but contains a massive spin-2 ghost [Stelle1977]. As an effective theory valid below the ghost mass scale $m_{2}$, there is no inconsistency with observations. For UV completion, several avenues exist: the fakeon prescription [AnselmiPiva2018], asymptotic safety [Percacci2017], and nonlocal form factors that soften the propagator without introducing additional poles [BarvinskyVilkovisky1990]. Deciding among these requires further work.

Third, the smallness of the observed cosmological constant remains a deep puzzle. In this framework, it appears as a renormalization condition selecting the vacuum. Whether a symmetry or a dynamical mechanism enforces the tiny value is an open question not addressed here.

Fourth, threshold matching in curved space, nonlocal form factors, and the interplay with finite-temperature and nonequilibrium effects in the early universe are technically rich and phenomenologically relevant. The covariant perturbation theory provides the conceptual tools, and explicit computations in the full SM plus gravity are feasible.

Despite these caveats, the self-consistency principle, embedded in the effective action and realized with standard quantum field theory techniques, yields a cohesive and predictive picture that connects inflationary cosmology, precision gravity, and particle physics within a single calculational framework.

IX Conclusion

Starting from the postulate that the classical backgrounds of gravity and the Standard Model are macroscopic order parameters of a unified quantum vacuum sustained by its own fluctuations, a coherent, calculable framework is constructed. The renormalized Einstein, Yang–Mills, and Higgs equations emerge as vacuum equations of state derived from the stationarity of the one-particle-irreducible effective action. The heat-kernel expansion fixes the renormalization of all couplings and requires the presence of higher-derivative gravitational terms. A solvable $O(N)$ model illustrates how the vacuum gap equations determine both curvature and condensates. The theory makes concrete predictions, notably Starobinsky-type inflation and universal quantum corrections to gravity, and it is consistent with current astrophysical and laboratory constraints. A systematic program of higher-loop calculations, functional RG studies, and phenomenological applications can now test and refine the picture.

Appendix A Conventions and useful variations

We use signature $(-,+,+,+)$ and define $R^{\rho}{}_{\sigma\mu\nu}=\partial_{\mu}\Gamma^{\rho}_{\nu\sigma}-\partial_{\nu}\Gamma^{\rho}_{\mu\sigma}+\cdots$, $R_{\mu\nu}=R^{\rho}{}_{\mu\rho\nu}$, and $R=\bar{g}^{\mu\nu}R_{\mu\nu}$. The Weyl tensor is

$$C_{\mu\nu\rho\sigma} =R_{\mu\nu\rho\sigma}-\frac{1}{2}\bigl(\bar{g}_{\mu\rho}R_{\nu\sigma}-\bar{g}_{\mu\sigma}R_{\nu\rho}$$

$$\quad-\bar{g}_{\nu\rho}R_{\mu\sigma}+\bar{g}_{\nu\sigma}R_{\mu\rho}\bigr)+\frac{R}{6}\bigl(\bar{g}_{\mu\rho}\bar{g}_{\nu\sigma}-\bar{g}_{\mu\sigma}\bar{g}_{\nu\rho}\bigr).$$ (47)

The basic variations used in the main text are

$$\delta\bigl(\sqrt{-\bar{g}}\,R\bigr) =\sqrt{-\bar{g}}\Bigl[\bigl(G_{\mu\nu}+\bar{g}_{\mu\nu}\Box-\nabla_{\mu}\nabla_{\nu}\bigr)\delta\bar{g}^{\mu\nu}\Bigr],$$ (48)

$$\delta\bigl(\sqrt{-\bar{g}}\,R^{2}\bigr) =\sqrt{-\bar{g}}\,H^{(R^{2})}_{\mu\nu}\delta\bar{g}^{\mu\nu},$$ (49)

$$\delta\bigl(\sqrt{-\bar{g}}\,C^{2}\bigr) =2\sqrt{-\bar{g}}\,B_{\mu\nu}\delta\bar{g}^{\mu\nu}.$$ (50)

Appendix B Heat-kernel coefficients for standard fields

For a real scalar with nonminimal coupling $\xi$, Dirac fermion, and gauge vector (in background gauge), the $a_{2}$ densities contributing to $C^{2}$ and $E_{4}$ are, respectively,

$$a_{2}^{\rm scalar} \supset\frac{1}{120}C^{2}-\frac{1}{360}E_{4}+\frac{1}{2}\Bigl(\xi-\tfrac{1}{6}\Bigr)^{2}R^{2},$$ (51)

$$a_{2}^{\rm Dirac} \supset\frac{1}{20}C^{2}-\frac{11}{360}E_{4},$$ (52)

$$a_{2}^{\rm vector} \supset\frac{1}{10}C^{2}-\frac{31}{180}E_{4},$$ (53)

in agreement with [Duff1994, Buchbinder1992]. Gauge and ghost contributions combine to the vector result.

Appendix C Bach tensor and equations of motion

The Bach tensor defined in (13) is traceless, $B^{\mu}{}_{\mu}=0$, and covariantly conserved in conformally flat spacetimes. In four dimensions, the metric variation of $\int d^{4}x\sqrt{-g}\,C^{2}$ yields $2\int d^{4}x\sqrt{-g}\,B_{\mu\nu}\delta g^{\mu\nu}$. The variation of $R^{2}$ is given by (14).

Appendix D Standard Model one-loop beta functions

At one loop in the $\overline{\rm MS}$ scheme, the SM gauge couplings obey

$$\mu\frac{dg_{1}}{d\mu} =\frac{41}{6}\frac{g_{1}^{3}}{16\pi^{2}},$$ (54)

$$\mu\frac{dg_{2}}{d\mu} =-\frac{19}{6}\frac{g_{2}^{3}}{16\pi^{2}},$$ (55)

$$\mu\frac{dg_{3}}{d\mu} =-7\frac{g_{3}^{3}}{16\pi^{2}},$$ (56)

with GUT normalization for $g_{1}$. The Higgs quartic $\lambda$ and the top Yukawa $y_{t}$ run as [MachacekVaughnI, MachacekVaughnII, Buttazzo2013]

$$\mu\frac{d\lambda}{d\mu} =\frac{1}{16\pi^{2}}\Bigl[24\lambda^{2}-6y_{t}^{4}+\frac{3}{8}\bigl(2g_{2}^{4}+(g_{2}^{2}+g_{1}^{2})^{2}\bigr)$$

$$\quad-9\lambda g_{2}^{2}-3\lambda g_{1}^{2}+12\lambda y_{t}^{2}\Bigr],$$ (57)

$$\mu\frac{dy_{t}}{d\mu} =\frac{y_{t}}{16\pi^{2}}\Bigl[\frac{9}{2}y_{t}^{2}-\frac{17}{12}g_{1}^{2}-\frac{9}{4}g_{2}^{2}-8g_{3}^{2}\Bigr].$$ (58)

Appendix E One-loop effective potential with curvature dependence

For a scalar field with field-dependent mass $m^{2}(\phi,R)=m_{0}^{2}(\phi)+\xi R$, the one-loop contribution to the effective potential in the limit of slowly varying backgrounds is

$$\Delta U =\frac{1}{64\pi^{2}}m^{4}(\phi,R)\Bigl[\ln\frac{m^{2}(\phi,R)}{\mu^{2}}-\frac{3}{2}\Bigr]$$

$$\quad+\frac{1}{192\pi^{2}}\Bigl(\xi-\frac{1}{6}\Bigr)m^{2}(\phi,R)R$$

$$\quad\times\Bigl[\ln\frac{m^{2}(\phi,R)}{\mu^{2}}-1\Bigr]+\cdots,$$ (59)

where dots denote higher-order terms in curvature and total derivatives [ParkerToms2009, Buchbinder1992, BirrellDavies1982]. This formula underlies (32) and the gap equations (36)–(37).