Capstones

The scalar capstone below connects a single action to a scattering rate, a one-loop running coupling, and a controlled low-energy expansion. It includes staged hints, a complete solution, numerical benchmarks, and criteria for checking your reasoning. Four further project blueprints offer other kinds of synthesis; they do not yet include complete solutions.

Choose one project, one result, and one stopping rule

Choose by the work you want to practice, not by the number of topics a project mentions.

Choose	Best fit and work product
Scalar action to running coupling — worked capstone	Connect canonical and functional formulations, scattering, one-loop renormalization, and heavy-field matching in one massive model.
Ward identity in Compton scattering	Make a diagram or sign error fail visibly; produce a two-diagram QED amplitude with polarization-replacement tests.
Heavy-field matching and running	Carry short-distance information across scales; produce a matched prediction with scale cancellation and a truncation check.
Infrared-safe observable and uncertainty	Combine real and virtual terms; produce a finite inclusive result with numerical and theoretical uncertainties separated.
Reproduce a lattice scalar fixture	Test a computational claim; produce a clean rerun, continuum fit, correlated uncertainty, and independent analytic check.

Before beginning, write three sentences:

1. Primary claim: the one result you will calculate, including its regime and requested accuracy.
2. Independent check: a comparison that does not reuse the decisive step of the main derivation or implementation.
3. Stopping rule: the first failed identity, uncontrolled limit, or missing input that will make you report a bounded partial result instead of extending the claim.

Adding another loop order, observable, field, or numerical method is not a substitute for finishing those three sentences.

Scalar action to running coupling

Question. What, precisely, can one normalized massive scalar action predict about two-particle scattering, scale dependence, and unresolved heavy physics? Produce a connected calculation in eight stages. The free-field identities are exact; the interacting amplitude is calculated through one loop; the heavy-field comparison is a separate tree-level expansion with an exact remainder. None of those results by itself constructs an interacting four-dimensional continuum theory.

Required background. Use classical actions to vary the model, canonical scalar quantization and Gaussian correlators to compare formulations, perturbation theory and LSZ to obtain a rate, and regularization, renormalization, and EFT matching to control the corrections. Helpful background. The compressed refresher helps identify a missing step without repeating the whole path.

Model, accuracy, and deliverable

In four-dimensional Minkowski vacuum, take a real scalar with

$$\mathcal L =\frac12(\partial\phi)^2-\frac12m^2\phi^2 -\frac{\lambda}{4!}\phi^4+\mathcal L_{\mathrm{ct}}.$$

Assume $m>0$, a weak positive coupling $\lambda$, and the unbroken $\phi\mapsto-\phi$ phase. The displayed classical potential is bounded below with its unique minimum at zero. The lightest odd particle is stable in this perturbative massive setting. Use compactly supported variations and smeared fields; introduce a finite regulator before evaluating coincident products. The global conventions apply. The interaction is not normal ordered, so its tadpole must be treated by mass renormalization.

For loops, use $d=4-2\epsilon$, subtract the quartic coupling in $\overline{\mathrm{MS}}$, and choose on-shell mass and field counterterms: $m$ is the pole mass and the renormalized field has unit pole residue. This mixed choice is intentional; an arbitrary interpolating field will be used separately to test LSZ. Require loop corrections and coupling-weighted logarithms to remain small compared with the preceding order.

Write your answer as eight short derivations with their checks. An initial working estimate is 6–10 hours after the listed lessons; this estimate has not been calibrated with learners. Work open-book, cite any borrowed result, and distinguish your calculation from a reproduced formula. Written equations with explanatory prose, or an oral derivation accompanied by an accessible transcript and equations, can demonstrate the same reasoning. Code is optional; the checks and benchmark inputs are specified below.

1. Action and states. Derive the field equation and canonical momentum; determine the oscillator coefficient from the equal-time commutator and normalize a one-particle state.
2. Two formulations. Complete a positive finite-dimensional Euclidean Gaussian, obtain its covariance and fourth moment, and recover the Lorentzian Feynman two-point function from both source differentiation and the canonical vacuum.
3. Interaction and LSZ. Count the quartic contractions, derive the leading connected four-point function, extract $\mathcal M$, and check a nonzero constant rescaling of the interpolating field.
4. A rate. Derive the identical-scalar differential and total cross sections from invariant flux and two-body phase space.
5. A loop and its cut. Evaluate one channel of the massive bubble with its symmetry factor. Determine its ultraviolet pole and absorptive part, and verify the optical theorem including identical intermediate particles.
6. Renormalization and running. Subtract all three channels, derive the beta function at fixed bare coupling, and check scale independence of the finite amplitude through the computed order.
7. Heavy-field matching. Add the heavy scalar specified in the solution, prove boundedness of its classical potential, match the quartic and first derivative operator, and bound the exact tree-level amplitude remainder.
8. The limit of the result. State what regulator removal would have to establish beyond this calculation, and the precise scope of the scalar triviality theorem used in the Research handoff.

Hints

Hint 1: choose the object at each stage

For stages 1–2, compare the same vacuum two-point distribution, not just two inverse matrices. For stage 3, keep the four external insertions labeled. For stage 4, decide whether your angular domain counts each final pair once or twice. For stage 5, distinguish a local ultraviolet pole from a physical two-particle cut. For stage 6, differentiate bare data before removing $\epsilon$. For stage 7, expand the heavy propagator only inside its convergence domain. For stage 8, a finite perturbative coefficient is not a limit of positive interacting measures.

Hint 2: identify the decisive operation

The equal-time commutator gives $2E_{\mathbf p}c_{\mathbf p}^2=1$. Translate a Euclidean Gaussian by $K_E^{-1}J$ and differentiate the normalized result. The labeled quartic contact has $4!$ attachments, while one labeled bubble channel has residual factor $1/2$. An external LSZ factor divides by the square root of the pole residue. The full-sphere rate needs $1/2!$. Feynman parametrization gives $m^2-sx(1-x)-i0$; find where this is negative. In $d=4-2\epsilon$, the quartic bare coupling carries $\mu^{2\epsilon}$. Complete the heavy-field square before expanding its inverse differential operator. In the constructive question, ask whether the connected fourth correlation survives the controlled continuum limit.

Hint 3: use independent normalization checks

Check $(\Box+m^2)D_F=-i\delta$ and $Z_0[J]=\exp(-JD_FJ/2)$ together. For the fourth Gaussian moment there are three pairings; for the interacting contact the coefficient is $-i\lambda$. Relativistic normalization gives $\int\mathrm d\Phi_2=\sqrt{1-4m^2/s}/(8\pi)$. The negative-logarithm interval in the bubble has exactly that square-root length, so the cut and phase-space calculation must agree. If the quartic simple-pole residue is $a\lambda^2$, its four-dimensional beta coefficient is $2a\lambda^2$. The identity $1/(1-z)-(1+z)=z^2/(1-z)$ supplies an exact EFT remainder. Separate this identity and the finite Gaussian from claims about all orders or an interacting continuum limit.

Worked solution

Open the complete scalar solution

1. Action, modes, and particle normalization

The first variation, with its surface term zero, gives

$$(\Box+m^2)\phi+\frac{\lambda}{3!}\phi^3=0, \qquad \pi=\dot\phi.$$

In four dimensions $[\phi]=1$, $[m]=1$, and $[\lambda]=0$; in the continued dimension $[\phi]=1-\epsilon$ and the bare quartic has dimension $2\epsilon$. These statements follow directly from the kinetic and interaction terms, using the action-principle derivation. The free Hamiltonian density is $[\pi^2+(\boldsymbol\nabla\phi)^2+m^2\phi^2]/2$.

Let $[a(\mathbf p),a^\dagger(\mathbf q)]=(2\pi)^3\delta^3(\mathbf p-\mathbf q)$ and $a(\mathbf p)|0\rangle=0$. Insert a real even coefficient $c_{\mathbf p}$ in the mode expansion. The two mixed commutators yield

$$[\phi(t,\mathbf x),\pi(t,\mathbf y)] =i\int\frac{\mathrm d^3\mathbf p}{(2\pi)^3} 2E_{\mathbf p}c_{\mathbf p}^2e^{i\mathbf p\cdot(\mathbf x-\mathbf y)}.$$

Thus $c_{\mathbf p}=1/\sqrt{2E_{\mathbf p}}$, and

$$\begin{aligned} \phi(x)&=\int\frac{\mathrm d^3\mathbf p}{(2\pi)^3\sqrt{2E_{\mathbf p}}} \left[a(\mathbf p)e^{-ipx}+a^\dagger(\mathbf p)e^{ipx}\right],\\ |\mathbf p\rangle&=\sqrt{2E_{\mathbf p}}a^\dagger(\mathbf p)|0\rangle,\\ \langle\mathbf q|\mathbf p\rangle &=2E_{\mathbf p}(2\pi)^3\delta^3(\mathbf p-\mathbf q). \end{aligned}$$

The canonical Volume treatment uses covariantly normalized oscillators instead. Its operator is $a_{\mathrm{cov}}(\mathbf p)=\sqrt{2E_{\mathbf p}}a(\mathbf p)$ and its measure contains $1/(2E_{\mathbf p})$; substituting this map gives exactly the same field and states. A measure and an oscillator commutator must be translated together.

2. The Gaussian and canonical answers coincide

With dimensionless finite-regulator variables $\varphi\in\mathbb R^N$, let $K_E=K_E^{\mathsf T}>0$ and $C_E=K_E^{-1}$. Completing the square gives

$$\begin{aligned} -\frac12\varphi^{\mathsf T}K_E\varphi+J^{\mathsf T}\varphi &=-\frac12(\varphi-C_EJ)^{\mathsf T}K_E(\varphi-C_EJ) +\frac12J^{\mathsf T}C_EJ,\\ Z_E[J]&=\frac{\mathcal Z_E[J]}{\mathcal Z_E[0]} =e^{J^{\mathsf T}C_EJ/2}. \end{aligned}$$

The determinant cancels because the same kernel, real integration domain, and measure occur in numerator and denominator. Two and four derivatives give

$$\begin{aligned} \langle\varphi_a\varphi_b\rangle&=(C_E)_{ab},\\ \langle\varphi_a\varphi_b\varphi_c\varphi_d\rangle &=(C_E)_{ab}(C_E)_{cd}+(C_E)_{ac}(C_E)_{bd} +(C_E)_{ad}(C_E)_{bc}. \end{aligned}$$

For the Lorentzian vacuum, the inherited source is $+J\phi$. Write $P=\Box+m^2$, $S_0=-\phi P\phi/2$, and choose its Feynman inverse with $P-i0$. Completion of the oscillatory square gives

$$Z_0[J]=\exp\!\left(-\frac12JD_FJ\right), \qquad D_F(p)=\frac{i}{p^2-m^2+i0}.$$

Here $JD_FJ$ abbreviates the two spacetime integrals. The insertion rule is $G^{(n)}=(-i)^n\delta^n Z/\delta J^n|_0$, so $-\delta^2Z_0|_0=D_F$. The fourth derivative reproduces the three Wick pairings, and its connected part vanishes. The sign is also fixed by $PD_F=-i\delta^{(4)}$. These operations use the Gaussian construction, source convention, and free Wick identity.

Independently, the canonical vacuum contraction gives

$$D_F(t,\mathbf x)=\int\frac{\mathrm d^3\mathbf p}{(2\pi)^3\,2E_{\mathbf p}} \left[\theta(t)e^{-iE_{\mathbf p}t+i\mathbf p\cdot\mathbf x} +\theta(-t)e^{iE_{\mathbf p}t-i\mathbf p\cdot\mathbf x}\right].$$

Closing the energy contour around the positive-energy pole below the real axis for $t>0$, or the negative-energy pole above it for $t<0$, recovers this same expression from $D_F(p)$. This checks the residue, phase, and boundary prescription; it is more informative than merely comparing denominators. The propagator treatment distinguishes this distribution from retarded and Wightman functions. A positive finite Gaussian is exact, while continuing an interacting formal functional integral requires additional justification.

3. Wick contractions, the quartic vertex, and LSZ

Expand the normalized Dyson expectation to first order. Connecting each of the four labeled external insertions to a different field at the quartic vertex gives $4!$ contractions. They cancel its $1/4!$, leaving

$$G_c^{(4)}(x_1,x_2,x_3,x_4) =-i\lambda\int\mathrm d^4z\prod_{j=1}^4D_F(x_j-z)+O(\lambda^2).$$

Vacuum components cancel against the zero-source denominator. Products of two external two-point functions instead belong to the disconnected four-point function and are removed when selecting $G_c^{(4)}$. Interacting Wick expansion and momentum-space rules therefore give, with all four momenta incoming,

$$\widetilde G_c^{(4)} =(-i\lambda)(2\pi)^4\delta^{(4)}\!\left(\sum_jp_j\right) \prod_{j=1}^4\frac{i}{p_j^2-m^2+i0}+O(\lambda^2).$$

Strip the single conservation delta and multiply each external pole by $(p_j^2-m^2)/i$ before taking the on-shell limit. In the convention $S_{fi,c}=i(2\pi)^4\delta^{(4)}(P_f-P_i)\mathcal M_{fi}$,

$$i\mathcal M_{\mathrm{tree}}=-i\lambda, \qquad \mathcal M_{\mathrm{tree}}=-\lambda.$$

This is dimensionless, invariant under the identical-particle permutations, and independent of $s,t,u$: the checks of a scalar contact amplitude.

For an arbitrary Hermitian interpolating field $\chi$, let $\langle\Omega|\chi|p\rangle=\sqrt Z>0$. The two-point pole is $iZ/(p^2-m^2+i0)$, but each external factor in the multipoint pole factorization is $i\sqrt Z/(p^2-m^2+i0)$. LSZ divides by $\sqrt Z$ per leg. Under $\chi'=c\chi$ with real $c\ne0$, $G_c^{(4)}$ acquires $c^4$ and $Z'=c^2Z$; four LSZ factors acquire $|c|^{-4}$ and cancel the change exactly. For odd particle multiplicity and $c<0$, the corresponding state phase must also be translated. The LSZ owner requires a stable isolated pole and asymptotic wave packets; a resonance or an infraparticle branch point would stop this argument.

4. Identical two-particle phase space

Put $s=(p_1+p_2)^2>4m^2$ and $v_s=\sqrt{1-4m^2/s}$. With the state normalization above, the invariant flux and two-body measure are

$$\begin{aligned} \mathcal F&=4\sqrt{(p_1\cdot p_2)^2-m^4}=2s v_s,\\ \mathrm d\Phi_2 &=(2\pi)^4\delta^{(4)}(P-p_3-p_4) \prod_{j=3}^4\frac{\mathrm d^3\mathbf p_j}{(2\pi)^3\,2E_j} =\frac{v_s}{32\pi^2}\mathrm d\Omega. \end{aligned}$$

The spatial delta eliminates $\mathbf p_4$. The remaining energy delta $\delta(\sqrt s-2\sqrt{p^2+m^2})$ supplies the radial Jacobian $E/(2p)$, giving the displayed measure. Thus the two equal-mass velocities cancel between phase space and flux. Over the full labeled sphere,

$$\frac{\mathrm d\sigma_{\mathrm{tree}}}{\mathrm d\Omega} =\frac{1}{\mathcal F}\frac{1}{2!}\lambda^2\frac{v_s}{32\pi^2} =\frac{\lambda^2}{128\pi^2s}, \qquad \sigma_{\mathrm{tree}}=\frac{\lambda^2}{32\pi s}.$$

The final $1/2!$ removes duplicate final configurations. There is no incoming $1/2!$ for two specified incident beams, and the vertex factorial was already used. Integrating over one permutation-ordered hemisphere with no final factor gives the same answer. The result is positive and has mass dimension $-2$. Its threshold value is a limit from $s>4m^2$; at exact threshold the incident flux itself vanishes.

5. The massive bubble and a physical cut

For one fixed $s$-channel partition there are $2(4\cdot3)^2\,2!=576$ contractions: choose which vertex receives the first pair, attach each labeled pair, and join the two remaining lines. Dividing by $2!(4!)^2=1152$ from the Dyson expansion gives $1/2$. The same counting and its renormalization interpretation appear in the massless illustrative calculation of Schwartz 2014, § 15.4, pp. 296–298. Here the mass and physical branch are retained throughout.

Factor one common $\mu^{2\epsilon}$ from the dimensionful amplitude; the following $\mathcal M$ has four-dimensional coupling normalization. For $p^2=s$,

$$\begin{aligned} i\mathcal M_s^{(1)} &=\frac{(-i\lambda)^2}{2}\mu^{2\epsilon} \int\frac{\mathrm d^d\ell}{(2\pi)^d} \frac{i}{\ell^2-m^2+i0} \frac{i}{(\ell+p)^2-m^2+i0}\\ &=\frac{i\lambda^2}{32\pi^2} \left[\frac1{\bar\epsilon}+F(s;\mu)\right]+O(\epsilon),\\ F(s;\mu)&=-\int_0^1\mathrm dx\, \ln\frac{m^2-sx(1-x)-i0}{\mu^2},\\ \frac1{\bar\epsilon}&=\frac1\epsilon-\gamma_E+\ln4\pi. \end{aligned}$$

To obtain the second line, combine denominators, shift $\ell\mapsto\ell+(1-x)p$, and Wick rotate the regulated squared denominator. Its integral without numerator factors is $i\Gamma(\epsilon)(4\pi)^{-2+\epsilon} [m^2-sx(1-x)-i0]^{-\epsilon}$. The two vertex phases and two propagator numerators multiply to $+\lambda^2$ before this final $i$. This supplies a direct sign check. The ultraviolet pole is independent of $s$ and is a local quartic counterterm in this amplitude. The mass prevents an infrared divergence. Dimensional regularization and subtraction explains the continued integral and the finite subtraction constants.

For $s>4m^2$, the logarithm’s argument is negative when $x_-<x<x_+$, with $x_\pm=(1\pm v_s)/2$. Since $\ln(-a-i0)=\ln a-i\pi$ for $a>0$,

$$\operatorname{Im}F(s;\mu)=\pi v_s, \qquad \operatorname{Im}\mathcal M_s^{(1)}=\frac{\lambda^2 v_s}{32\pi}.$$

In physical elastic scattering $t,u\le0$, so their bubbles have no absorptive part. Independently, unitarity in the forward direction gives

$$2\operatorname{Im}\mathcal M_{\mathrm{forward}}^{(1)} =\frac1{2!}\int\mathrm d\Phi_2\, |\mathcal M_{\mathrm{tree}}|^2 =\frac{\lambda^2 v_s}{16\pi} =\mathcal F\sigma_{\mathrm{tree}}.$$

Both calculations agree. A wrong logarithm branch reverses the sign; a missing intermediate-particle factorial doubles the right side. A real local counterterm cannot cancel this physical cut.

6. Subtraction, the beta function, and scale cancellation

The three channels have the same ultraviolet pole. With $\mathcal L_{\mathrm{ct}}\supset-\mu^{2\epsilon}\delta\lambda\phi^4/4!$, choose

$$\begin{aligned} \delta\lambda&=\frac{3\lambda^2}{32\pi^2\bar\epsilon},\\ \mathcal M(s,t,u;\mu) &=-\lambda(\mu)+\frac{\lambda(\mu)^2}{32\pi^2} \big[F(s;\mu)+F(t;\mu)+F(u;\mu)\big]+O(\lambda^3). \end{aligned}$$

The one-loop tadpole is independent of external momentum. The on-shell mass counterterm cancels that two-point correction at the pole, and no wave-function derivative counterterm is needed at order $\lambda$. Consequently it contributes no additional one-loop four-point amplitude after external pole and residue normalization. Vacuum energy is removed from the normalized correlators. These facts explain why the three bubbles and the quartic counterterm suffice here.

Through this order, $\lambda_0=\mu^{2\epsilon} [\lambda+a\lambda^2/\bar\epsilon+O(\lambda^3)]$, where $a=3/(32\pi^2)$. Let $\beta_d=-2\epsilon\lambda+b\lambda^2+O(\lambda^3)$. Differentiate at fixed bare data before taking $\epsilon\to0$:

$$\begin{aligned} 0&=2\epsilon\left(\lambda+\frac{a\lambda^2}{\bar\epsilon}\right) +\beta_d\left(1+\frac{2a\lambda}{\bar\epsilon}\right) +O(\lambda^3),\\ 0&=(b-2a)\lambda^2+O(\epsilon\lambda^2,\lambda^3),\\ \beta_\lambda&=\mu\frac{\mathrm d\lambda}{\mathrm d\mu}\bigg|_0 =\frac{3\lambda^2}{16\pi^2}+O(\lambda^3). \end{aligned}$$

This implements the simple-pole extraction rule. Dropping the canonical $-2\epsilon\lambda$ term too early loses the factor of two. A source using $d=4-\epsilon_{\mathrm{source}}$ must be translated with $\epsilon_{\mathrm{source}}=2\epsilon$.

Each $F$ has $\partial F/\partial\ln\mu=2$ at fixed pole mass and external momenta. Hence

$$\frac{\mathrm d\mathcal M}{\mathrm d\ln\mu} =-\frac{3\lambda^2}{16\pi^2} +\frac{\lambda^2}{32\pi^2}(2+2+2)+O(\lambda^3) =O(\lambda^3).$$

The one-loop running solution, with specified boundary data, is

$$\lambda(\mu) =\frac{\lambda(\mu_0)} {1-\dfrac{3\lambda(\mu_0)}{16\pi^2}\ln(\mu/\mu_0)}.$$

Use it only where the running coupling remains small. Its formal Landau singularity marks failure of this approximation and is not a proof of a continuum no-go theorem. Coupling values and finite subtraction constants depend on scheme; predictions expressed in the same physical inputs agree through the retained order. Analytic finite redefinitions $\lambda'=\lambda+c\lambda^2+\cdots$ preserve this leading beta coefficient, while changing higher coefficients and finite amplitude terms.

The rate through one-loop interference would retain $|\mathcal M|^2=\lambda^2-2\lambda\operatorname{Re}\mathcal M^{(1)} +O(\lambda^4)$. Squaring the truncated amplitude without re-expanding also includes selected higher-order terms; it is not a complete two-loop rate.

7. A stable heavy-field model and its exact remainder

Now ask whether the same light quartic can arise from unresolved heavy exchange. At tree level add a real scalar $H$ with

$$\mathcal L_{\mathrm{full}} =\frac12(\partial\phi)^2+\frac12(\partial H)^2-U(\phi,H),$$ $$\begin{aligned} U&=\frac12m^2\phi^2+\frac12M^2H^2 +\frac g2H\phi^2+\frac{\lambda_{\mathrm{full}}}{24}\phi^4\\ &=\frac12M^2\left(H+\frac{g\phi^2}{2M^2}\right)^2 +\frac12m^2\phi^2 +\frac{\lambda_{\mathrm{full}}-3g^2/M^2}{24}\phi^4. \end{aligned}$$

Take $M\gg m$, real $g$ with mass dimension one, and $\lambda_{\mathrm{full}}>3g^2/M^2$. This sufficient condition proves boundedness and a unique zero-field classical minimum; the negative quartic shift below does not imply an unstable model. Let $g^2/(16\pi^2M^2)$ and $\lambda_{\mathrm{full}}/(16\pi^2)$ be small. The heavy particle need not be stable when $M>2m$; it is an internal mediator, never an external LSZ state in this exercise.

Solving $(\Box+M^2)H=-g\phi^2/2$ with the scattering boundary prescription, or completing the heavy-field square, gives

$$\begin{aligned} \Delta\mathcal L_{\mathrm{tree}} &=\frac{g^2}{8}\phi^2\frac1{M^2+\Box}\phi^2\\ &=\frac{g^2}{8M^2}\phi^4 -\frac{g^2}{8M^4}\phi^2\Box(\phi^2)+\cdots,\\ \lambda(\mu_M)&=\lambda_{\mathrm{full}}(\mu_M)-\frac{3g^2}{M^2} \quad\text{at tree matching},\qquad \mu_M\simeq M. \end{aligned}$$

The first derivative term is kept in this unreduced operator basis. The corresponding full and truncated light-particle amplitudes are

$$\begin{aligned} \mathcal M_{\mathrm{full}}^{\mathrm{tree}} &=-\lambda_{\mathrm{full}} +g^2\sum_{x\in\{s,t,u\}}\frac1{M^2-x},\\ \mathcal M_{\mathrm{EFT}}^{\mathrm{tree}} &=-\lambda+\frac{g^2}{M^4}(s+t+u) =-\lambda+\frac{4g^2m^2}{M^4}. \end{aligned}$$

Each exchange follows from $(-ig)^2i/(x-M^2+i0)$, so the positive sign in $g^2/(M^2-x)$ is fixed. On shell, the derivative contribution is constant because $s+t+u=4m^2$; it vanishes in the massless limit. This kinematic fact does not remove the off-shell operator. An equation-of-motion basis change must also transform the quartic and higher interactions.

For real invariants with $\rho=\max(|s|,|t|,|u|)/M^2<1$, the geometric identity gives an exact tree-level absolute error:

$$\begin{aligned} R&=\mathcal M_{\mathrm{full}}^{\mathrm{tree}} -\mathcal M_{\mathrm{EFT}}^{\mathrm{tree}} =\frac{g^2}{M^4}\sum_{x\in\{s,t,u\}}\frac{x^2}{M^2-x},\\ |R|&\le\frac{g^2(s^2+t^2+u^2)}{M^6(1-\rho)}. \end{aligned}$$

For these real invariants each denominator is positive, so $R\ge0$ as well. With $Q^2=\max(|s|,|t|,|u|)$, the remainder is of order $(g^2/M^2)(Q/M)^4$. A small absolute error need not be a small relative error near a cancellation of the full amplitude. This sharp bound belongs to the specified tree model, not to arbitrary unknown EFT coefficients; the controlled-expansion treatment explains that distinction.

Below the threshold, the leading quartic runs with the light-scalar beta function already derived. Tree matching plus this one-loop evolution resums the corresponding leading logarithms. It does not include finite one-loop heavy threshold terms, loops with derivative operators, or their operator mixing. Those would be needed for a complete matched prediction at that combined order. The exact error bound above compares the two trees only; it cannot bound those omitted quantum contributions.

8. What a constructive continuum limit would add

The finite Gaussian defines a probability measure. A stable finite lattice interaction also defines an ordinary finite-dimensional integral. Neither statement proves that the interacting Schwinger functions have a non-Gaussian limit when lattice spacing $a\to0$ and volume $L\to\infty$ with physical parameters held fixed. Such a construction must establish the limiting distributions and their covariance, locality after reconstruction, reflection positivity, regularity and growth conditions, and the required control of both limits. For the specified normalized scalar scaling field, a surviving connected correlation of order higher than two excludes a Gaussian limit for that field. This test cannot be applied indiscriminately to nonlinear composites: the normal-ordered square of a free scalar already has a nonzero connected three-point function. The definition of a QFT explains why the action is only part of these data.

There is a rigorous obstruction for a specified scalar class: Aizenman and Duminil-Copin 2021, arXiv v4, § 1.2, pp. 3–5, PDF consider hypercubic lattice fields with a constant nearest-neighbor ferromagnetic coupling and single-site weights proportional to $e^{-\lambda_{\mathrm{lat}}\varphi^4+b\varphi^2}$, including the Ising limiting case. In four dimensions their continuum scaling limits are Gaussian if the long-distance two-point correlation vanishes. The limit concerns joint distributions of smeared fields, normalized by the square root of the box-sum variance: first the system-size to observation-scale ratio $R/\ell\to\infty$, then $\ell\to\infty$. Model parameters may vary along this sequence. The theorem covers arbitrary positive lattice quartic strength in that class; it does not specify every possible limiting covariance, nor prove that all conceivable four-dimensional QFTs are Gaussian. This nonperturbative theorem is stronger than extrapolating our one-loop running formula.

The Aizenman and Duminil-Copin 2024 corrigendum, p. 479 corrects estimates in the proof; the theorem’s conclusions remain unchanged.

The conclusion of this capstone is therefore a convention-consistent perturbative prediction and a controlled tree matching example. The precise next question is what survives the continuum limit for a declared regulator family and observable class. Continue to the four-dimensional constructive-QFT question for its evidence and the distinct gauge-theory problem.

Benchmarks and self-review

For a reproducible massive fixture choose $m$ as the energy unit, $\lambda(m)=0.2$, $s=9m^2$, and $\cos\theta=1/3$. Then $t=-5m^2/3$, $u=-10m^2/3$, and $s+t+u=4m^2$. The tree rate and absorptive part must satisfy

$$\sigma_{\mathrm{tree}}m^2=\frac{0.04}{288\pi} \simeq4.420971\times10^{-5}, \qquad \operatorname{Im}\mathcal M^{(1)}=\frac{0.04\sqrt5}{96\pi} \simeq2.965677\times10^{-4}.$$

The finite one-loop amplitude at $\mu=m$ is $\mathcal M\simeq-0.200012997+0.000296568\,i$, and one-loop evolution gives $\lambda(2m)\simeq0.200528120$. These decimals describe the specified truncated formula; their printed precision is not an estimate of omitted loop effects.

For the separate heavy-tree comparison, set $M=10m$, $g=m$, and $\lambda_{\mathrm{full}}=0.23$, so the tree matched quartic is $0.20$. Treat these as tree parameters; this does not identify the loop coupling at $\mu=m$ with a matched boundary at $\mu=M$ without running. The truncated tree amplitude is $-0.1996$. The full tree amplitude is approximately $-0.199497504$, its remainder is $1.024959\times10^{-4}$, and the absolute bound is $1.042735\times10^{-4}$. Keep enough precision to resolve that difference before subtracting the two amplitudes.

For each row below, assign 0 if the result or its domain is missing, 1 if the result is obtained but the independent check or explanation is missing, and 2 if both are correct. This is a transparent self-review guide, not a calibrated examination or certification. The suggested stopping rule is to obtain 2 in every row before treating the chain as complete; an error in a propagator sign, symmetry factor, or limit cannot be compensated by points elsewhere.

Stage	Evidence for a complete answer	If the check fails, revisit
1. Action and states	Correct field equation, $2Ec^2=1$, and relativistic state norm	Scalar variation and Mode and state normalization
2. Formulations	Finite positive Gaussian, correct source sign, three fourth-moment pairings, and canonical/Feynman agreement	Finite Gaussian, Feynman source prescription, and Wick pairings
3. Interaction and LSZ	$4!$ attachments, $\mathcal M=-\lambda$, and invariant result under $\chi\mapsto c\chi$	Quartic contractions and LSZ residues
4. Rate	Correct flux, phase space, and exactly one final-state counting correction	Flux, phase space, and identical final states
5. Loop and cut	Bubble factor $1/2$, local UV pole, positive cut, and optical-theorem agreement	Bubble counting, Regulated bubble and cut, and Identical-scalar cut and symmetry factor (the two-body speed is denoted ρ there)
6. Running	Three-channel subtraction, fixed-bare beta derivation, and cancellation of the order-$\lambda^2$ scale derivative	Three-channel subtraction and Beta function and scale cancellation
7. Matching	Stable full potential, correct coefficient signs, remainder bound, and separate quantum omissions	Stable heavy-field matching and Remainder and omitted effects
8. Scope	Exact/free, perturbative, tree-expansion, and constructive claims distinguished; theorem hypotheses stated	Action versus QFT data, Constructive limit, and Scalar theorem scope

You have demonstrated a bounded scalar synthesis when every link is reproducible and every approximation has its stopping condition. This does not establish proficiency in the fermion, gauge, infrared, or other branches required by the complete graduate core.

Project blueprint 2: Ward identity in Compton scattering

Question. Can you derive the two tree diagrams for $e^-\gamma\to e^-\gamma$ and show that gauge-representative dependence cancels only when their signs, fermion ordering, and propagator momenta are consistent?

Concrete deliverable. Produce the tensor amplitude $\mathcal M^{\mu\nu}$, demonstrate both external-photon Ward checks, and compute either one fixed-polarization rate or the unpolarized differential cross section in a stated frame. Include a one-page table that maps each algebraic cancellation to the assumption it uses.

Preparation. Review fermions and spin, vector fields and gauge redundancy, symmetry, currents, and Ward identities, QED and Yang–Mills theory, and LSZ and tree amplitudes.

Minimum required steps.

1. Fix the covariant derivative, photon polarization, spinor, and gamma-matrix conventions.
2. Derive the two ordered fermion-line contributions with their internal momenta rather than importing an unlabelled formula.
3. Contract the incoming and outgoing photon indices with their momenta and use momentum conservation plus the on-shell Dirac equations.
4. Choose physical polarizations or a justified polarization sum and compute the selected rate.
5. Repeat the result after shifting either polarization by a multiple of its photon momentum.

Decisive checks. Replacing an external polarization by its momentum must give zero. Deliberately omit either diagram and confirm that the test fails. Check the electron mass dimension, the soft-photon behavior appropriate to the selected observable, and agreement between a trace calculation and at least one explicit polarization configuration.

Scope and failure boundary. This calculation tests a tree-level QED Ward identity with on-shell external states. It does not establish loop-level renormalization, anomaly cancellation, detector inclusiveness, or non-Abelian Slavnov–Taylor identities. Stop if a polarization shift changes the rate or if the kinematic limit leaves the assumptions used in the derivation.

Optional extension. Cross the amplitude into pair annihilation or add the leading soft-emission contribution and state the inclusive measurement needed for a finite prediction.

Project blueprint 3: heavy-field matching and running

Question. For a light scalar $\phi$ coupled to a heavy scalar $H$, can you replace heavy exchange by local operators, evolve the coefficients from the matching scale to a low scale, and predict an amplitude with a tested error in $Q/M$?

Concrete deliverable. Starting from the heavy-scalar model in the Core EFT lesson, derive the exact tree-level exchange amplitude, match the leading and next derivative operators near $\mu_M\simeq M$, run the retained coefficient or coefficient vector to $\mu_Q\simeq Q$, and compare the EFT prediction with the expanded full-theory result at several values of $Q/M$.

Preparation. Use perturbative expansion and Feynman rules, loops and regularization, renormalization and the renormalization group, and effective field theory and matching. For operator mixing or several scales, follow renormalization and EFT for working researchers.

Minimum required steps.

1. State the light degrees of freedom, symmetry, hierarchy, scheme, operator basis, and matching observable.
2. Expand the heavy propagator through a declared power in $Q/M$ and match all operators that contribute at that order.
3. Derive or independently verify the anomalous dimension needed for the retained light-theory operators.
4. Solve the RG equation between $\mu_M$ and $\mu_Q$ and combine the evolved coefficient with the low-energy matrix element at a consistent order.
5. Compare with the full-theory expansion and tabulate the observed error as $Q/M$ and the retained order change.

Decisive checks. Check coefficient dimensions and signs, reproduce the first omitted power in a full-theory comparison, and show cancellation of explicit and implicit renormalization-scale dependence through the calculated order. Translate the answer through one allowed finite scheme change and back.

Scope and failure boundary. The local expansion applies below the heavy threshold while $Q/M$ and the relevant couplings remain small. It cannot reproduce the heavy pole at finite order. Residual scale variation diagnoses omitted perturbative terms but is not automatically a probability interval.

Optional extension. Perform one-loop matching, include a second operator that mixes under running, or test whether a symmetry removes the nominally leading correction.

Project blueprint 4: infrared-safe observable and uncertainty

Question. Can you define one inclusive rate or event-shape bin so that its soft and collinear limits are physically meaningful, combine real and virtual terms into a finite result, and state what each uncertainty component tests?

Concrete deliverable. Reproduce one benchmark fixed-order inclusive rate or normalized event-shape bin. Provide the measurement function or cut definition, cancellation or subtraction calculation, phase-space integration, normalization, and a table separating Monte Carlo uncertainty, input covariance, residual scale dependence, and any nonperturbative or detector correction actually used.

Preparation. Begin with infrared-safe observables and synthesis and the scattering phenomenology pathway. For the calculation itself, use Real–Virtual Cancellation and Local Subtraction and Phase-Space Integration and Monte Carlo Estimators.

Minimum required steps.

1. Define the observable and show how its measurement function behaves when a particle becomes soft or a pair becomes collinear.
2. Display the virtual poles and the integrated or local counterterms that cancel them, using one regulator and convention consistently.
3. Integrate the finite remainder with recorded cuts, maps, tolerances, and random-stream construction.
4. Normalize the result and vary renormalization and factorization scales coherently where both occur.
5. Report numerical, parametric, perturbative, matching, and model effects separately before choosing any combination rule.

Decisive checks. Verify independence from unphysical subtraction or slicing parameters over a resolved range, reproduce a known inclusive or soft/collinear limit, repeat the integration with a different phase-space map or independent stream, and show that real–virtual poles cancel before numerical rounding obscures them.

Scope and failure boundary. The conclusion applies only to the declared observable, cuts, fixed order, scale range, and treatment of long-distance physics. Scale variation is a sensitivity test, not a universal confidence level. Do not compare directly with detector-level data unless the transfer, hadronization, and calibration assumptions are included.

Optional extension. Add logarithmic resummation and verify fixed-order expansion plus matching, or compare two infrared-safe measurement definitions that emphasize different physical scales.

Project blueprint 5: reproduce a lattice scalar fixture

Question. Can you reproduce a finite-lattice scalar pole mass in a clean environment, recover its continuum order, propagate correlated uncertainty, and distinguish successful execution from evidence for the stated numerical claim?

Concrete deliverable. Reproduce the free-scalar fixture in computational field theory onboarding at several values of $am_0$. Extract $E_a$ once from the lattice propagator and once from a periodic correlator fit, extrapolate $E_a/m_0$ to $a=0$, and provide a complete rerun packet plus a separately derived analytic reference.

Preparation. Use the numerical and reproducibility review, Lattice Regulators and Continuum Targets, Statistical Inference and Error Budgets, and Reproduce and Validate a Result.

Minimum required steps.

1. Write the finite mathematical specification, units, boundary conditions, expected $O(a^2)$ artifact, and acceptance tolerances before writing code.
2. Implement direct momentum-space and correlator-based extractions without sharing the decisive mass-extraction routine.
3. Sweep lattice spacing, temporal extent, numeric precision, fit window, and statistics; use the full stable covariance for correlated fits.
4. Fit the predicted continuum order, test plausible higher-order terms, and compare with the exact pole relation.
5. Rerun from recorded source, dependencies, command, inputs, random streams, raw outputs, and analysis choices in a clean environment.

Decisive checks. Recover the sign and coefficient of the leading free cutoff term, observe second-order convergence in its asymptotic range, inject an incorrect momentum stencil or periodic-image omission and confirm that a test fails, and verify that a diagonal-covariance fit does not silently replace the stated analysis.

Scope and failure boundary. This project verifies a free scalar regulator, implementation, and continuum-inference procedure. It is not evidence that an interacting sampler equilibrates, that a sign problem is controlled, or that the same fit model resolves excited states. Report those as new questions.

Optional extension. Add a weak interaction, Hamiltonian truncation, or a tensor-network calculation while retaining the exact free fixture and adding method-specific truncation checks.

Copyable final-report template

Use this template for any blueprint. Keep the strongest supported conclusion near the top so a reader does not have to infer it from the final plot.

project title:

question:
strongest supported conclusion:
scope and failure boundary:

theory, state, conventions, and boundary conditions:
observable, normalization, units, and kinematic regime:
preparation pages actually used:

derivation or computational method:
approximations, truncations, and first omitted terms:
inputs and their covariance:

primary result:
uncertainty components and their meanings:
decisive checks and whether each passed:
independent cross-check and shared assumptions:

source, environment, data, and rerun instructions:
unexpected failure or negative result:
stopping rule reached:

optional extension not attempted:
next exact page or research question:

A useful report remains valuable when the target calculation fails: preserve the failed check, reduce the claim, and explain which missing derivation, input, or limiting sequence would be needed next. Use Scope a First Research Project to turn that bounded gap into a new question, or return to Learning pathways to choose a specialist route for the completed result.

References

• Aizenman, Michael, and Hugo Duminil-Copin. “Marginal Triviality of the Scaling Limits of Critical 4D Ising and $\phi^4_4$ Models.” Annals of Mathematics 194, no. 1 (2021): 163–235. DOI. Open PDF, arXiv v4. Corrigendum, Annals of Mathematics 199, no. 1 (2024): 479. Corrigendum DOI; Open PDF, complete one-page correction in publisher preview.
• Schwartz, Matthew D. Quantum Field Theory and the Standard Model. Cambridge: Cambridge University Press, 2014. DOI.