Skip to content

Effective Average Actions and the Wetterich Equation

The effective average action Γk\Gamma_k is a scale-dependent 1PI generator constructed by adding an infrared quadratic kernel before taking a modified Legendre transform. For real bosons it obeys

tΓk[φ]=12Tr[(Γk(2)[φ]+Rk)1tRk],tlnkΛUV.\boxed{ \partial_t\Gamma_k[\varphi] = \frac12\operatorname{Tr} \left[ \left(\Gamma_k^{(2)}[\varphi]+R_k\right)^{-1} \partial_tR_k \right] }, \qquad t\equiv\ln\frac{k}{\Lambda_{\rm UV}}.

The equation is exact because the inverse contains the full running Hessian. Its one-trace appearance does not make a finite ansatz exact. This page derives the sign, the inverse Hessian, the bosonic factor 1/21/2, the fermionic supertrace qualification, both endpoints, and a scalar local-potential projection whose one-loop coefficient can be checked analytically.

Required background. Wilsonian Coarse Graining and Theory Space supplies the physical scale interpretation and the exact-versus-projected distinction. The 1PI Effective Action and Mean-Field Equations supplies the ordinary Legendre transform, connected Hessian, and 1PI inverse-propagator relation.

Helpful background. The Polchinski Exact RG Equation develops the complementary Wilson-action formulation and fixes this chapter’s tt convention.

For a real Euclidean bosonic field χ\chi, add

ΔSk[χ]=12χRkχ\Delta S_k[\chi] = \frac12\chi\cdot R_k\cdot\chi

to the microscopic action. With fAgf\cdot A\cdot g denoting integration over spacetime or momentum and contraction of internal indices, define

Zk[J]=DχeS[χ]ΔSk[χ]+Jχ,Wk[J]lnZk[J].Z_k[J] = \int\mathcal D\chi\, e^{-S[\chi]-\Delta S_k[\chi]+J\cdot\chi}, \qquad W_k[J]\equiv\ln Z_k[J].

The average field and connected two-point function are

φ=δWkδJ,Gk=δ2WkδJδJ=χχJφφ.\varphi = \frac{\delta W_k}{\delta J}, \qquad G_k = \frac{\delta^2W_k}{\delta J\,\delta J} = \langle\chi\chi\rangle_J-\varphi\varphi.

For a bosonic scalar, a useful regulator satisfies:

RequirementOperational meaning
Rk(q)>0R_k(q)>0 for q2k2q^2\ll k^2Slow modes acquire an effective gap and do not yet fluctuate freely.
Rk(q)0R_k(q)\to0 for q2/k2q^2/k^2\to\inftyModes well above kk are not distorted by the infrared regulator.
Rk(q)0R_k(q)\to0 as k0k\to0The ordinary unregulated generating functional is recovered at the infrared endpoint.
RkR_k large at the starting scaleFluctuations can be suppressed enough to match ΓkUV\Gamma_{k_{\rm UV}} to declared microscopic data.
(tRk)/(Γk(2)+Rk)(\partial_tR_k)/(\Gamma_k^{(2)}+R_k) ultraviolet integrableThe trace is finite or is accompanied by an explicit ultraviolet regulator and renormalization prescription.

Positivity is a bosonic stability condition, not a universal sign rule for every field species. Gauge fields, ghosts, and fermions require regulator matrices compatible with their quadratic forms and symmetry identities.

Let J=Jk[φ]J=J_k[\varphi] solve φ=δWk/δJ\varphi=\delta W_k/\delta J. The effective average action is

Γk[φ]=JφWk[J]12φRkφ.\Gamma_k[\varphi] = J\cdot\varphi-W_k[J] -\frac12\varphi\cdot R_k\cdot\varphi.

Subtracting the regulator term is essential. Variation at fixed kk gives

δΓkδφ=JRkφ.\frac{\delta\Gamma_k}{\delta\varphi} = J-R_k\varphi.

Differentiate once more and use δφ/δJ=Gk\delta\varphi/\delta J=G_k:

Γk(2)+Rk=δJδφ=Gk1,Gk=(Γk(2)+Rk)1.\Gamma_k^{(2)}+R_k = \frac{\delta J}{\delta\varphi} = G_k^{-1}, \qquad G_k = \left(\Gamma_k^{(2)}+R_k\right)^{-1}.

This is the exact inverse-Hessian identity. A zero eigenvalue of Γk(2)+Rk\Gamma_k^{(2)}+R_k makes the flow singular, so the domain on which this operator is invertible is part of any existence or numerical claim.

Wetterich defines the infrared-regulated Legendre functional, proves this Hessian relation, and derives the exact scale equation in Wetterich 1993, pp. 90–94.

Differentiate WkW_k at fixed source:

tWk[J]J=12χtRkχJ=12φtRkφ12Tr(tRkGk).\begin{aligned} \left.\partial_tW_k[J]\right|_J &= -\frac12 \left\langle \chi\cdot\partial_tR_k\cdot\chi \right\rangle_J\\ &= -\frac12\varphi\cdot\partial_tR_k\cdot\varphi -\frac12\operatorname{Tr}\left(\partial_tR_k\,G_k\right). \end{aligned}

At fixed φ\varphi, the terms containing tJ\partial_tJ cancel between JφJ\cdot\varphi and Wk[J]W_k[J]. Therefore

tΓk[φ]φ=tWk[J]J12φtRkφ=12Tr(tRkGk)=12Tr[(Γk(2)+Rk)1tRk].\begin{aligned} \left.\partial_t\Gamma_k[\varphi]\right|_\varphi &= -\left.\partial_tW_k[J]\right|_J -\frac12\varphi\cdot\partial_tR_k\cdot\varphi\\ &= \frac12\operatorname{Tr}\left(\partial_tR_k\,G_k\right)\\ &= \frac12\operatorname{Tr} \left[ \left(\Gamma_k^{(2)}+R_k\right)^{-1} \partial_tR_k \right]. \end{aligned}

No loop expansion was used. The full GkG_k depends on all vertices through Γk(2)\Gamma_k^{(2)}. Replacing that Hessian by the Hessian of a finite polynomial, derivative expansion, or vertex ansatz is where approximation enters.

For a field multiplet, a compact convention is

tΓk=12STr[(Γk(2)+Rk)1tRk].\partial_t\Gamma_k = \frac12\operatorname{STr} \left[ \left(\Gamma_k^{(2)}+\mathcal R_k\right)^{-1} \partial_t\mathcal R_k \right].

The supertrace includes momentum and internal labels and inserts a minus sign for Grassmann-odd blocks. In a doubled fermionic field basis the displayed factor 1/21/2 applies to the whole supermatrix. With an undoubled complex pair ψˉ,ψ\bar\psi,\psi and regulator ψˉRkψ\bar\psi R_k\psi, the same content is conventionally written as a fermionic term Tr[(Γk(2)+Rk)1tRk]-\operatorname{Tr}[(\Gamma_k^{(2)}+R_k)^{-1}\partial_tR_k] without the bosonic half. The field ordering and doubling convention must accompany the formula.

The figure reappears here to emphasize its right branch. Inspect the endpoint boxes: sending Rk0R_k\to0 removes the modification from the Legendre transform, whereas sending Ck0C_k\to0 in a Wilsonian representation does not turn its interaction action into Γk\Gamma_k.

Polchinski and Wetterich flows use different functionals and kernels; both exact identities sit above a separate dashed finite-projection layer, and only the Wetterich branch ends at the full 1PI action.

The solid branches are schematic exact identities for a regulated bosonic theory with t=ln(k/ΛUV)t=\ln(k/\Lambda_{\rm UV}). The effective average action uses an additive infrared kernel RkR_k and tends to the ordinary 1PI action only when Rk0R_k\to0 and the limit exists. A compatible modified Legendre map can relate it to a Wilsonian action, but the kernels, fields, boundary data, and endpoints remain distinct. The dashed PNP_N layer marks the optional finite truncation, not part of either exact derivation.

ItemWilson interaction SkIS_k^{\rm I}Effective average action Γk\Gamma_k
Independent fieldRetained integration field ϕ\phiMean field φ=δWk/δJ\varphi=\delta W_k/\delta J
Scale kernelCovariance CkC_kAdditive inverse-propagator term RkR_k
Exact flow structureFirst-derivative product minus second-derivative traceFull inverse-Hessian trace
Infrared statementWilson vertices encode eliminated fluctuations; no automatic 1PI endpointRk0R_k\to0 gives ΓkΓ\Gamma_k\to\Gamma under controlled limits
Typical projectionMomentum-dependent Wilson verticesEffective potential, derivative expansion, or 1PI vertices

For complementary cutoff data, a generalized Legendre transform relates Wilsonian and infrared-regulated 1PI functionals Morris 1994, § 3, eqs. (3.25)–(3.28). This relation justifies the dotted cross-arrow; it does not identify the two boxes.

Ultraviolet matching and infrared endpoint

Section titled “Ultraviolet matching and infrared endpoint”

If Rk(q)0R_k(q)\to0 for every fixed qq as k0k\to0, then

Γk0[φ]=Γ[φ],\Gamma_{k\to0}[\varphi] = \Gamma[\varphi],

provided the regulated functionals converge and the Legendre transform remains well defined. Under the usual positivity assumptions, the endpoint effective action has the convexity properties of a Legendre transform. A finite-kk potential need not yet be convex because Γk(2)\Gamma_k^{(2)} is bounded together with RkR_k rather than by itself.

At the opposite endpoint, a large positive RkR_k suppresses fluctuations and makes a saddle approximation accurate. In an idealized kk\to\infty limit this gives the field-dependent part of the microscopic action. In a theory initialized at a finite ultraviolet scale,

ΓkUV=Smicro+cutoff and matching corrections.\Gamma_{k_{\rm UV}} = S_{\rm micro} +\text{cutoff and matching corrections}.

Thus ΓkUV=Smicro\Gamma_{k_{\rm UV}}=S_{\rm micro} is a boundary approximation whose accuracy depends on the ultraviolet cutoff, regulator, masses, and normalization. Wetterich states both the Rk0R_k\to0 endpoint and the large-kk classical limit, including the finite-cutoff qualification, in Wetterich 1993, pp. 91–93.

For one Z2\mathbb Z_2-symmetric scalar, take the local-potential ansatz

Γk[φ]=ddx[12(μφ)2+Uk(φ)].\Gamma_k[\varphi] = \int d^dx \left[ \frac12(\partial_\mu\varphi)^2 +U_k(\varphi) \right].

This ansatz fixes the wave-function coefficient to one and omits higher derivatives. At a constant field,

Γk(2)(q,q;φ)=q2+Uk(φ),\Gamma_k^{(2)}(q,-q;\varphi) = q^2+U_k''(\varphi),

so the functional equation projects to

tUk(φ)=12qtRk(q)q2+Rk(q)+Uk(φ).\partial_tU_k(\varphi) = \frac12\int_q \frac{\partial_tR_k(q)} {q^2+R_k(q)+U_k''(\varphi)}.

For an analytic evaluation, choose

Rk(q)=(k2q2)θ(k2q2).R_k(q) = (k^2-q^2)\theta(k^2-q^2).

Then tRk=2k2θ(k2q2)\partial_tR_k=2k^2\theta(k^2-q^2); the distribution from differentiating the step function is multiplied by k2q2k^2-q^2 and vanishes at its support. Define

AdSd1d(2π)d=2d(4π)d/2Γ(d/2).A_d \equiv \frac{S_{d-1}}{d(2\pi)^d} = \frac{2}{d(4\pi)^{d/2}\Gamma(d/2)}.

The momentum integral is elementary:

tUk(φ)=Adkd+2k2+Uk(φ).\boxed{ \partial_tU_k(\varphi) = A_d\frac{k^{d+2}}{k^2+U_k''(\varphi)} }.

The parent equation is exact, but this local-potential equation inherits the ansatz and regulator choice. The compact-support profile is convenient at this order; derivative projections require additional care at its nonsmooth boundary.

Expand

Uk(φ)=Ek+12mk2φ2+λk4!φ4+g6,k6!φ6+,U_k(\varphi) = \mathcal E_k +\frac12m_k^2\varphi^2 +\frac{\lambda_k}{4!}\varphi^4 +\frac{g_{6,k}}{6!}\varphi^6 +\cdots,

and set Pk=k2+mk2P_k=k^2+m_k^2. Matching powers of φ\varphi gives

tEk=Adkd+2Pk,tmk2=Adkd+2λkPk2,tλk=Adkd+2(6λk2Pk3g6,kPk2).\begin{aligned} \partial_t\mathcal E_k &= A_d\frac{k^{d+2}}{P_k},\\ \partial_tm_k^2 &= -A_d\frac{k^{d+2}\lambda_k}{P_k^2},\\ \partial_t\lambda_k &= A_dk^{d+2} \left( \frac{6\lambda_k^2}{P_k^3} -\frac{g_{6,k}}{P_k^2} \right). \end{aligned}

The running vacuum term Ek\mathcal E_k is required for an absolute free energy. The quartic flow depends on g6,kg_{6,k}, so a quartic polynomial is not an invariant subspace. Nevertheless, g6g_6 first affects the weak-coupling quartic beta function beyond O(λ2)O(\lambda^2). At d=4d=4, mk2=0m_k^2=0, and this order,

A4=132π2,dλkdlnk=3λk216π2+O(λk3),A_4=\frac{1}{32\pi^2}, \qquad \frac{d\lambda_k}{d\ln k} = \frac{3\lambda_k^2}{16\pi^2} +O(\lambda_k^3),

reproducing the shell and Polchinski checks.

Zero dimensions remove momentum dependence without removing the Legendre transform or regulator logic. For

Zk(J)=dχexp ⁣[S(χ)12Rkχ2+Jχ],Z_k(J) = \int_{-\infty}^{\infty}d\chi\, \exp\!\left[ -S(\chi)-\frac12R_k\chi^2+J\chi \right],

the trace has one entry and the exact equation is

tΓk(φ)=12tRkΓk(φ)+Rk.\partial_t\Gamma_k(\varphi) = \frac12 \frac{\partial_tR_k} {\Gamma_k''(\varphi)+R_k}.

For the Gaussian action S(χ)=m2χ2/2S(\chi)=m^2\chi^2/2, direct integration gives

Γk(φ)=12m2φ2+12ln(m2+Rk)+constant independent of k.\Gamma_k(\varphi) = \frac12m^2\varphi^2 +\frac12\ln(m^2+R_k) +\text{constant independent of }k.

Its derivative is exactly (tRk)/[2(m2+Rk)](\partial_tR_k)/[2(m^2+R_k)], matching the flow because Γk=m2\Gamma_k''=m^2. This checks the sign, factor 1/21/2, inverse Hessian, endpoint, and running vacuum contribution independently of a differential solver.

A reproducible calculation extends this benchmark to a quartic integral: exact quadrature is compared with matched polynomial flows, a running vacuum term, nested truncations, and more than one regulator. Zero dimensions do not test momentum or derivative expansions, but they sharply expose unmatched boundary data, Hessian singularities, and the mistake of calling a polynomial solution exact.

  • Regulator removal. Setting Rk=0R_k=0 before differentiating makes tRk=0\partial_tR_k=0 and removes the flow; the endpoint is reached by integration followed by a controlled limit.
  • Gaussian check. The zero-dimensional result verifies the complete normalization-sensitive equation, not only its field derivatives.
  • One-loop check. The scalar local-potential projection reproduces 3/(16π2)3/(16\pi^2) in four dimensions with the declared coupling normalization.
  • Exactness ceiling. The functional trace is exact, while the local-potential ansatz, polynomial expansion, regulator profile, and numerical solution are separate choices.
  • Scope ceiling. The derivation is Euclidean and formal at the functional-integral level. Gauge identities, real-time contours, global existence, and subject-specific phenomenology require additional structures.

Using the ordinary Legendre transform at finite kk. Omitting the subtraction ΔSk[φ]\Delta S_k[\varphi] changes the stationarity and Hessian identities. The flow then contains extra field-quadratic terms and is not the displayed Wetterich equation.

Reading “one trace” as “one loop only.” The line in the trace is the full field- and scale-dependent inverse Hessian. Expanding that inverse in microscopic couplings generates arbitrarily high perturbative loop orders.

Equating finite ultraviolet data with the classical action automatically. At a finite starting scale, regulator and cutoff corrections can remain. State how ΓkUV\Gamma_{k_{\rm UV}} is matched and test the sensitivity to moving that scale.

Derive Gk=(Γk(2)+Rk)1G_k=(\Gamma_k^{(2)}+R_k)^{-1} from the modified Legendre transform.

Solution

Stationarity gives δΓk/δφ=JRkφ\delta\Gamma_k/\delta\varphi=J-R_k\varphi. Differentiating with respect to φ\varphi gives Γk(2)+Rk=δJ/δφ\Gamma_k^{(2)}+R_k=\delta J/\delta\varphi. Since δφ/δJ=Wk(2)=Gk\delta\varphi/\delta J=W_k^{(2)}=G_k, the two Jacobians are inverse operators, proving the identity.

For the compact-support regulator, derive the local-potential flow and the coefficient AdA_d.

Solution

Inside q<k|q|<k, q2+Rk(q)=k2q^2+R_k(q)=k^2 and tRk=2k2\partial_tR_k=2k^2. Outside, the numerator vanishes. Hence

tUk=k2k2+Ukq<kddq(2π)d.\partial_tU_k = \frac{k^2}{k^2+U_k''} \int_{|q|<k}\frac{d^dq}{(2\pi)^d}.

The ball volume is Sd1kd/dS_{d-1}k^d/d, so the remaining factor is AdkdA_dk^d. Multiplying by k2/(k2+Uk)k^2/(k^2+U_k'') gives the displayed result.

Show why omitting the vacuum term fails the zero-dimensional Gaussian benchmark.

Solution

The field-dependent Gaussian action has Γk=m2\Gamma_k''=m^2, so the exact right-hand side is the nonzero field-independent quantity (tRk)/[2(m2+Rk)](\partial_tR_k)/[2(m^2+R_k)]. A truncation containing only m2φ2/2m^2\varphi^2/2 has zero left-hand side at fixed m2m^2 and cannot satisfy the equation. Adding Ek=12ln(m2+Rk)\mathcal E_k=\frac12\ln(m^2+R_k) restores equality.

  • Morris, Tim R. “The Exact Renormalization Group and Approximate Solutions.” International Journal of Modern Physics A 9 (1994): 2411–2450. DOI. Open PDF.
  • Wetterich, Christof. “Exact Evolution Equation for the Effective Potential.” Physics Letters B 301 (1993): 90–94. DOI. Open PDF.