Codex @codex  0

For an autonomous Lagrangian, the canonical momentum and energy areAlong a solution of the Euler-Lagrange equation,The absence of explicit time dependence is essential to this conservation law.

The classical action is $S_c(q_I,q_F,T)={\int }_0^{T}L(q_c,{\dot{q}}_c)dt$. If the final endpoint is moved along one fixed classical trajectory, while the initial endpoint is held fixed, the fundamental theorem of calculus givesThis total derivative is different from the partial derivative with fixed final position.

For a family of nearby classical trajectories, integration by parts gives the endpoint variationTo obtain the final term, the endpoint change in the variation at fixed parameter time is $\delta q(T)=\delta q_F-{\dot{q}}_F\delta T$, while the changing upper integration limit contributes $L_F\delta T$. The interior integral vanishes by the Euler-Lagrange equation. Therefore, on any smooth branch of classical boundary-value solutions,The chain rule along the endpoint trajectory now reads $d{S}_c/dT=p_F{\dot{q}}_F-E=L_F$, consistently with the direct calculation.

The canonical momentum is $p_i=\partial L/\partial {\dot{x}}_i=m{\dot{x}}_i+q{A}_i$. Its time derivative contains the total derivative along the trajectory: ${\dot{p}}_i=m{\ddot{x}}_i+q{\partial }_t{A}_i+q{\dot{x}}_j{\partial }_j{A}_i$. The Euler-Lagrange equation therefore givesThe magnetic field is $\mathbf{B}=\nabla \times \mathbf{A}$, so the last term is $(\dot{\mathbf{x}}\times \mathbf{B}{)}_i$. With the electric field $\mathbf{E}=-\nabla \varphi -{\partial }_t\mathbf{A}$, this is precisely the Lorentz force law:

For fields ${\Phi }_a$ with a first-derivative Lagrangian density, suppose an infinitesimal global transformation at fixed coordinates is $\delta {\Phi }_a=\alpha {\Delta }_a$ and $\delta \mathcal{L}=\alpha {\partial }_{\mu }{K}^{\mu }$, with constant infinitesimal $\alpha$. Noether theorem gives the on shell conserved currentThis follows by varying the Lagrangian, integrating the terms with ${\partial }_{\mu }{\Delta }_a$ by parts, and setting the Euler-Lagrange expressions to zero. If spatial boundary flux vanishes, $\int d^{3}x{J}^0$ is a conserved Noether charge. Spacetime symmetries are covered by including the induced field variation at fixed coordinates and the corresponding total derivative $K^{\mu }$.

An example of a theory with a Lorentz four-vector is the source-free Maxwell field with electromagnetic four-potential $A_{\mu }$. Under an Abelian gauge transformation,The change in $F_{\mu \nu }$ is ${\partial }_{\mu }{\partial }_{\nu }\chi -{\partial }_{\nu }{\partial }_{\mu }\chi =0$, so the Maxwell Lagrangian is gauge-invariant. For $A_{\mu }=(\Phi ,-\mathbf{A})$ define $E_i=F_{0i}$ and $F_{ij}=-{ϵ}_{ijk}{B}_k$. With signature $+---$,The minus sign in $-F^2/4$ gives the positive electric kinetic term and positive physical field energy; reversing it would reverse the Hamiltonian's sign for the physical modes. The normalization $1/4$ accounts for antisymmetry and yields the conventional equations. Varying $A_{\nu }$ gives ${\partial }_{\mu }{F}^{\mu \nu }=0$.

For time translation at fixed coordinates take ${\Delta }_{\nu }={\partial }_0{A}_{\nu }$. Since $\alpha$ is constant, $\delta F_{\mu \nu }=\alpha {\partial }_0{F}_{\mu \nu }$, and hence $\delta \mathcal{L}=\alpha {\partial }_0\mathcal{L}$. Thus $K^{\mu }={\delta }^{\mu }{}_0\mathcal{L}$. The momentum derivative isNoether theorem therefore gives the canonical stress-energy tensor's time-translation currentHere ${\dot{A}}_i=E_i+{\partial }_i{A}_0$ was used. The canonical momentum of $A_0$ vanishes, while that of the lower spatial component $A_i$ is $E_i$: $A_0$ acts as a multiplier for Gauss law, not a propagating degree of freedom. The density above is not manifestly gauge-invariant, but its charge agrees with the physical energy on the constraint surface and with appropriate boundary conditions.

The gauge-covariant time translation instead has ${\Delta }_{\nu }={\partial }_0{A}_{\nu }-{\partial }_{\nu }{A}_0=F_{0\nu }$. It is time translation combined with a gauge transformation of parameter $\chi =-\alpha A_0$. Directly,using the definition of $F$, equivalently its Bianchi identity. The total derivative in $\delta \mathcal{L}$ is unchanged. The improved current isIt is the gauge-invariant Maxwell stress-energy tensor. In particular,The spatial current is the Poynting vector $\mathbf{E}\times \mathbf{B}$, so conservation gives ${\partial }_{t}u+\nabla \cdot (\mathbf{E}\times \mathbf{B})=0$. This is the familiar positive local electromagnetic field energy density, expressed entirely in measurable fields.

To compare the currents,The last term vanishes on the source-free Maxwell equations, leaving a stress-energy tensor improvement generated by a stress-energy superpotential. At the charge level,when the surface term vanishes and Gauss law holds. Thus the symmetries yield the same conserved energy under these conditions while the second supplies a gauge-invariant local density. Boundary terms would have to be retained if those boundary conditions were not imposed.

Take $\lambda \ne 0$ and use the complex sine-Gordon theory reduced density $\ell =\mathcal{L}/M^2$, with physical lengths and times $x/M,t/M$. Physical rest energy is $M\int dx{\mathcal{H}}_{\ell }$, and the dimensionless Noether charge is obtained from the reduced action. Put $c=\cos \alpha$, $s=\sin \alpha$, $f(x)=c\mathrm{sech}(cx)/\lambda$, so $\psi =f{e}^{ist}$. Initially take $0<|\alpha |<\pi /2$, avoiding the singular target-space boundary.

The canonical Hamiltonian density isHere $|\dot{\psi }{|}^2=s^2{f}^2$, $|{\psi }_x{|}^2=c^2{f}^2{\tanh }^2(cx)$, and $D=1-c^2{\mathrm{sech}}^2(cx)$. Since $s^2+c^2{\tanh }^2(cx)=D$, the kinetic contribution is simply $f^2$. ConsequentlyTo fix the charge sign, take the Noether current directly from $\delta \psi =i\epsilon \psi$, $\delta {\psi }^{\ast }=-i\epsilon {\psi }^{\ast }$:Thus positive internal rotation has positive charge in this convention. With $X=cx$ and the supplied integral,Reversing the definition of the generator would reverse all charge labels but not any mass. The limit $\alpha \to \pi /2$ is the vacuum with zero mass and charge. At $\alpha =0$ the field reaches $|\psi |=1/|\lambda |$ at its center and the current's coordinate expression becomes singular. The mass still has the finite limiting value $4M/{\lambda }^2$, butThe displayed local density does not assign a unique endpoint charge merely by substituting $\alpha =0$.

The four collective parameters. Translation, a Lorentz boost of rapidity $b$, and the global phase symmetry produceThe independent real parameters are $\alpha ,b,X_0,\chi$, with $\chi$ periodic modulo $2\pi$. A further time translation is absorbed into these parameters. The energy and momentum are $\mathcal{M}\cosh b$ and $\mathcal{M}\sinh b$, while $Q$ is invariant. The internal angle is the collective coordinate whose canonical momentum is this Noether charge.

Semiclassical quantization. For one rest-frame internal rotation, the period in dimensionless time is $T=2\pi /|s|$. The action around the orbit isEquivalently this is $S+\mathcal{E}T$, where $\mathcal{E}$ is the reduced rest energy. Bohr-Sommerfeld quantization with $\hslash =1$ gives $I=2\pi n$. Because the internal orbit is a cyclic angle, not a one-dimensional oscillator between turning points, there is no half-integer Maslov index shift. HenceOn the nonsingular branch $0<n<2\pi /{\lambda }^2$. The elementary charge-one mass approaches $M$ at weak coupling. For positive $a,b$ with $(a+b){\lambda }^2/4\le \pi /2$, $\sin a^{\prime }+\sin b^{\prime }>\sin (a^{\prime }+b^{\prime })$, where $a^{\prime }=a{\lambda }^2/4$, $b^{\prime }=b{\lambda }^2/4$. Thus ${\mathcal{M}}_{a+b}<{\mathcal{M}}_a+{\mathcal{M}}_b$: the charge-$n$ soliton cannot fragment into smaller like-charge particles. For a general charge-conserving partition, a constituent with $|q_i|\ge n$ already has mass at least ${\mathcal{M}}_n$; if every $|q_i|<n$, the decreasing ratio ${\mathcal{M}}_q/q$ gives ${\sum }_i{\mathcal{M}}_{|q_i|}\ge ({\mathcal{M}}_n/n){\sum }_i|q_i|\ge {\mathcal{M}}_n$, strictly for nontrivial fragmentation. Thus opposite-charge constituents cannot open a lower-energy decay either. This is charged rather than topological stability, and the semiclassical spectrum is reliable at weak coupling.

Counting states and the endpoint. Write $L=2\pi /{\lambda }^2$. If $L$ is not an integer, the regular branch givesIf $L$ is an integer, the local nonsingular calculation gives only $2(L-1)$ regular states. The formal closed-branch list has $2L$ charge labels, but its last two labels are limits of the same singular static solution. Counting them needs a global completion; it cannot be decided from a regular current integral at $\alpha =0$.

In the additional integer-level completion, ${\lambda }^2=4\pi /k$ at leading semiclassical order and charges are identified modulo integer $k>1$. There are then $k-1$ nonzero charge sectors. For odd $k$, positive and negative labels are distinct; for even $k$, the endpoint labels $\pm k/2$ describe one self-conjugate state. This supplies $2L-1$ states when $k=2L$, rather than $2L$ endpoint labels. This global interpretation, and the distinction between classical and renormalized couplings, are discussed in Dorey and Hollowood, section 2. A generic classical $\lambda$ does not automatically define an integer $k$ or specify this quantum completion.

Use signature $(-,+,+,+)$ and write $h_{ij}$ for the spatial induced metric, $D_i$ for its spatial covariant derivative, and $q^2=h^{ij}{D}_i\varphi D_j\varphi$. The negative shift vector convention gives $n^{\mu }=N^{-1}(1,N^i)$ and $\Pi =n^{\mu }{\partial }_{\mu }\varphi$. Decomposing the scalar field derivative givesVarying the scalar field matter action with respect to the inverse metric tensor, including the variation of $\sqrt{-g}$, gives the stress-energy tensorRaising both indices gives the requested $T^{\mu \nu }$. Contracting with the unit normal twice, or once with the spatial projection tensor, givesProjecting both indices of the stress-energy tensor onto the spatial hypersurface givesTwo printed identifications require correction. The canonical momentum of the displayed density is ${\pi }_{\varphi }=\sqrt{h}\Pi$, not $\Pi$: explicitly ${\mathcal{L}}_m=N\sqrt{h}[{\Pi }^2/2-q^2/2-V]$ and differentiation with respect to $\dot{\varphi }$ supplies $\sqrt{h}$. Thus $\Pi$ is the normal scalar-field momentum. Also, the printed trace-free stress has an extra factor $1/2$ and uses ${\delta }_{ij}$ where the general spatial induced metric requires $h_{ij}$. In an orthonormal spatial frame the latter becomes ${\delta }_{ij}$, but the coefficient is still one. For a gradient $(q,0,0)$ in that frame, direct projection gives ${\tilde{S}}_{11}=2q^2/3$, rather than the printed $q^2/3$. The PDF's energy-density gradient is spatial ${\partial }_i\varphi {\partial }^i\varphi$; the TeX aid's spacetime index there is an OCR defect.

Let $X={\partial }^{-2}\zeta$, with an inverse spatial Laplacian defined by a boundary or zero-mode prescription. It depends on $\zeta$, not on its time derivative. Differentiation gives the canonical momentumCount $\zeta ,\pi ,X$ as first order. Inverting perturbatively givesFor the Legendre transform in mechanics, write $v=\dot{\zeta }=\pi +\Delta$ with $\Delta$ second order. The quadratic kinetic contribution is $\pi v-v^2/2={\pi }^2/2-{\Delta }^2/2$, whose correction is fourth order; in the cubic terms one may consequently replace $v$ by $\pi$. The Hamiltonian density is thereforeThe cubic interaction Hamiltonian density isIn the free interaction picture, ${\pi }_I={\dot{\zeta }}_I$. Thus here ${\mathcal{H}}_{\mathrm{int}}=-{\mathcal{L}}_3$ with the time derivatives interpreted as free fields. This equality follows from the cubic Legendre transform for scalar derivative interactions; at higher orders additional terms can appear.

For the real scalar field, the mostly-plus Minkowski metric gives a positive kinetic energy. Its canonical momentum and Hamiltonian operator areThe Euler-Lagrange field equation is $(\square -m^2)\varphi =0$, or $({\partial }_t^2-{\nabla }^2+m^2)\varphi =0$. Canonical quantization of a real scalar field promotes the fields to Hermitian operators and replaces the equal-time Poisson brackets by canonical commutation relations:The Heisenberg equation of motion then gives $\dot{\varphi }=\pi$ and $\dot{\pi }={\nabla }^2\varphi -m^2\varphi$, so the classical field equation remains an operator identity. Products at coincident points need a regulator; equivalently, the field is an operator-valued distribution.

Resolve the free field into positive- and negative-frequency plane waves. With $E_{\mathbf{k}}=\sqrt{|\mathbf{k}{|}^2+m^2}$ and $k\cdot x=-E_{\mathbf{k}}t+\mathbf{k}\cdot \mathbf{x}$,Reality pairs the two terms by Hermitian conjugation. The equal-time canonical commutation relations are equivalent toFor example, the two mixed terms in $[\varphi ,\pi ]$ each supply half the Dirac delta function. Thus the normalization $1/\sqrt{2E_{\mathbf{k}}}$ is fixed, not optional. Each independent mode is a quantum harmonic oscillator, with $a$ and $a^{†}$ its annihilation operator and creation operator.

Choose the Fock vacuum by $a(\mathbf{k})|0⟩=0$. Repeated creation operators construct the bosonic Fock space. Substitution into the Hamiltonian operator yieldswhere a finite volume $V$ can be used as a regulator. Normal ordering sets the flat-space vacuum reference to zero. The spatial momentum operator isIn particular, $[H,a^{†}(\mathbf{k})]=E_{\mathbf{k}}{a}^{†}(\mathbf{k})$ and $[\mathbf{P},a^{†}(\mathbf{k})]=\mathbf{k}{a}^{†}(\mathbf{k})$. A one-particle state therefore has positive energy, momentum $\mathbf{k}$, and invariant mass $m$. A scalar field transforms in the trivial spin representation, so its particles have spin zero. Because the field is real there is only one set of oscillators: the particle is its own antiparticle, with no distinct conserved particle-minus-antiparticle charge. The free theory does conserve its occupation-number sum, though generic scalar interactions need not.

A momentum eigenstate extends throughout space. A localized particle is described by a superposition such asIts wave packet evolves through the phase $e^{-i{E}_{\mathbf{k}}t}$; a narrow packet has group velocity $\mathbf{v}=\mathbf{k}/E_{\mathbf{k}}$, whose magnitude is at most one. Localization and dispersion concern these superpositions, rather than classical trajectories attached to individual field modes.

Relativistic spacetime behavior is encoded by microcausality. Directly from the oscillator commutators,This is a Lorentz-invariant distribution and vanishes at equal time. Any spacelike separation can be transformed to an equal-time separation, so $[\varphi (x),\varphi (y)]=0$ for spacelike $x-y$. Thus local operations at spacelike-separated points commute. The vacuum Wightman function need not vanish outside the light cone, but that correlation does not transmit a controllable signal. The Feynman propagator iswhich propagates positive energy forward and negative energy backward in the time-ordered product. For causal response one uses the retarded Green function, whose support lies in or on the future light cone.

Finally, commuting creation operators implyThe multiparticle state is symmetric under the particle exchange operator, which is Bose statistics. For a normalized single mode, the occupation states are $(a^{†}{)}^n|0⟩/\sqrt{n!}$ with $n=0,1,2,\dots$: any number of identical bosons can occupy it. There is no Pauli exclusion principle for these particles. This construction agrees with the Spin-statistics theorem relating integer spin to bosonic exchange symmetry under relativistic locality and positive-energy assumptions. It demonstrates the required scalar case without assuming that theorem as the quantization prescription.

Use natural units and the Minkowski metric $\eta =\mathrm{diag}(1,-1,-1,-1)$. A real scalar field assigns a real variable $\varphi (t,\mathbf{x})$ to each spatial point. Its Lagrangian density can be taken to beThe principle of stationary action gives the Euler-Lagrange equation $\square \varphi +V^{\prime }(\varphi )=0$. With $V(\varphi )=m^2{\varphi }^2/2$, this is the Klein-Gordon equation. An additional nonlinear part of $V$ describes interactions.

The canonical momentum is $\pi =\partial \mathcal{L}/\partial \dot{\varphi }=\dot{\varphi }$. The Legendre transform in mechanics gives the canonical Hamiltonian density of a real scalar fieldThe Hamiltonian equations $\dot{\varphi }=\pi$ and $\dot{\pi }={\nabla }^2\varphi -V^{\prime }(\varphi )$ recover the same field equation. In canonical quantization, the fields become operators satisfying the equal-time canonical commutation relationsA spatial lattice makes the analogy with many coupled quantum-mechanical coordinates precise. Each lattice field value is a coordinate, with its own conjugate momentum. The path integral is another representation of the same quantum evolution.

To see its origin, first consider one coordinate with $H=p^2/(2M)+V(q)$. Split a time interval into $N$ steps of length $\epsilon$ and insert position and momentum resolutions of the identity. The short-time kernel isMultiplying the kernels and integrating over intermediate positions gives the phase-space path integralThe endpoints of $q$ are fixed. The momentum integrals are Gaussian integrals; completing the square produces the configuration-space path integralAt finite slicing its normalization contains $(M/(2\pi i\epsilon ){)}^{N/2}{\prod }_{j=1}^{N-1}d{q}_j$. This fixes the composition law and the initial delta-function kernel. One sums over all paths, not merely solutions of the classical equation. Restoring $\hslash$ replaces the weight by $e^{iS/\hslash }$; stationary phase explains the emergence of classical trajectories.

For the field, use scalar field configuration eigenstates $|\phi ⟩$, satisfying $\hat{\varphi }(\mathbf{x})|\phi ⟩=\phi (\mathbf{x})|\phi ⟩$. Insert their completeness relations on every time slice. This givesThe endpoint field configurations are fixed. Integrating the Gaussian momentum variables leaves the scalar field path integral $\mathcal{N}\int \mathcal{D}\varphi e^{iS[\varphi ]}$. The functional measure means a regulated product over the field variables. A spacetime lattice or another ultraviolet cutoff makes this product finite before the continuum limit; interacting continuum calculations may require renormalization. The oscillatory Minkowski weight is an amplitude, not a positive probability density.

For vacuum expectation values, the boundaries must select the vacuum rather than arbitrary field configurations. Long imaginary-time evolution suppresses excited states: $e^{-TH}|n⟩=e^{-T{E}_n}|n⟩$, so after normalization only the lowest-energy component remains as $T\to \infty$. This is vacuum projection by imaginary time. The corresponding Feynman i-epsilon prescription in the real-time integral specifies the vacuum boundary conditions and the poles of the propagator. With $t=-i\tau$, the Euclidean path integral has the weight $e^{-S_E}$, whereIt is often a useful regulated starting point; analytic continuation returns the vacuum time-ordered quantities.

Introduce a classical source $J(x)$ and define the normalized vacuum generating functionalwith the same vacuum prescription in numerator and denominator. A functional derivative brings down $i\varphi (x)$. The order of the time slices makes the operator insertion time-ordered. Thus source differentiation inserts time-ordered field operators:The denominator removes vacuum diagrams and gives normalized expectation values. It is essential that these are time-ordered products; differentiating this vacuum functional does not directly give every possible operator ordering.

The free theory illustrates the method. Its quadratic kernel is $K=-\square -m^2$ with the vacuum pole prescription, and completing the square gives the Gaussian evaluation of a free scalar generating functionalTwo source derivatives give ${\Delta }_F=⟨\Omega |T\{\hat{\varphi }(x)\hat{\varphi }(y)\}|\Omega ⟩$. Higher derivatives give all pairings, the content of Wick theorem. For an interaction $S_{\mathrm{int}}$, one may use path-integral perturbation by source derivatives:Expanding this expression generates Feynman diagrams and their Wick contractions. The connected generating functional $W[J]=-i\log Z[J]$ retains connected contributions; in particular $\delta W/\delta J=⟨\varphi {⟩}_J$. These functionals turn the computation of field-operator expectations into source differentiation of an ordinary regulated integral.

The canonical momentum is $p=L_{\dot{q}}$, and a regular Legendre transform in mechanics defines $H(q,p,t)=p\dot{q}-L(q,\dot{q},t)$ after solving for $\dot{q}$. Its differential isbecause the $d\dot{q}$ terms cancel. With the Euler-Lagrange equation $\dot{p}=L_q$, this gives Hamilton's equationsIn particular $H_q=-L_q$ at corresponding velocities, so a cyclic coordinate for the Lagrangian is also cyclic for the Hamiltonian.

Use units $\hslash =c=1$, the Minkowski metric ${\eta }_{ab}=\mathrm{diag}(1,-1,-1,-1)$, and the Fourier transform convention $f(x)=\int d^{4}p{e}^{-ip\cdot x}\tilde{f}(p)/(2\pi {)}^4$. Here gamma matrices satisfy $\{{\gamma }^a,{\gamma }^b\}=2{\eta }^{ab}{I}_4$, and $\bar{\psi }={\psi }^{†}{\gamma }^0$ is the Dirac adjoint. A single Dirac field describes both Electrons and Positrons; these are the particle and antiparticle sectors of that field.

The free Maxwell Lagrangian and Dirac action, with a linear covariant gauge condition for the photon, giveBefore gauge fixing, the Maxwell Lagrangian has zero modes in field theory along $A_a\mapsto A_a+{\partial }_a\lambda$, so its quadratic operator cannot be inverted on all potentials. In the linear Lorenz gauge, the Faddeev-Popov determinant is $\det \square$. It is independent of $A$ and may be absorbed into the normalization; the corresponding Faddeev-Popov ghost fields have no interacting vertices in this Abelian linear gauge. A residual gauge symmetry is removed by the specified boundary conditions.

The free generating functional isThe Dirac field variables and their sources are independent Grassmann fields in this path integral. They anticommute; replacing them by ordinary commuting fields would give the wrong statistics and the wrong functional determinant. The bosonic Gaussian functional integral contributes an inverse square root of a determinant, and the Grassmann Gaussian integral contributes a determinant. If $K_A$ is the photon quadratic operator and $K_D=i{\gamma }^a{\partial }_a-m$, completing the square givesIn the last expression $D_F=i{K}_A^{-1}$ and $S_F=i{K}_D^{-1}$ are the photon propagator and Dirac propagator, with Feynman i-epsilon prescription. The products include the appropriate spacetime integrals and index contractions. Functional derivatives with respect to $J$, and consistently ordered left or right Grassmann derivatives with respect to the fermionic sources, generate the time ordering of the corresponding fields. In Feynman gauge, $\alpha =1$, the momentum-space two-point functions areThe Feynman i-epsilon prescription specifies vacuum boundary conditions rather than an arbitrary inverse of the differential operator.

The covariant photon propagator uses four potential components. Its operator formalism counterpart is Gupta-Bleuler quantization: impose $({\partial }_a{A}^a{)}^{(+)}|\mathrm{phys}⟩=0$ and take the Gupta-Bleuler null-state quotient. This leaves a positive physical state space with two transverse photon polarization vectors. The temporal and longitudinal oscillator components occur in intermediate covariant expressions; the Ward identity removes their dependence from physical amplitudes. The free Dirac field uses the canonical anticommutation relations, producing the same fermionic signs as its Grassmann fields in the path integral.

The operator-path-integral equivalence can be seen directly with a regulator. Divide time into small intervals and insert complete sets of field-coordinate states for bosons, and resolutions in fermionic coherent states for fermions. The bosonic matrix elements produce the phase-space factor $\exp i\sum [\pi \Delta \varphi -H\Delta t]$; integrating out the quadratic canonical momentum produces the bosonic action. The fermionic coherent state overlaps produce the first-order term $i{\psi }^{†}\dot{\psi }$ and the Berezin integral measure. Multiplying the short-time kernels recovers the path integral. Projecting the remote endpoints onto the Fock vacuum with an infinitesimal damping selects the same Feynman propagators as the operator formalism. Field insertions become time-ordered products under this construction. Conversely, their quadratic generating functional obeys the canonical free-field equations and has precisely the oscillator two-point functions, so its higher free correlators agree by the Wick theorem. This establishes the equivalence for the regulated free theory and order by order in the perturbation series.

To couple the matter field electromagnetically, promote its global phase symmetry to the local transformationThe gauge covariant derivative $D_a={\partial }_a+ie{A}_a$ obeys $D_a^{\prime }{\psi }^{\prime }=e^{-ie\lambda }{D}_a\psi$. Replacing ${\partial }_a$ by $D_a$ in the Dirac action therefore gives the invariant matter densityThe electromagnetic field tensor is unchanged by the local transformation. Thus the unfixed quantum electrodynamics action is gauge-invariant. The Dirac current $j^a=\bar{\psi }{\gamma }^a\psi$ is conserved by the matter equations, and the Electron and Positron excitations carry opposite charges. The added gauge fixing density selects a representative and is not itself invariant under arbitrary local transformations; it does not change gauge-invariant observables. At a free fermion vertex,For external on-shell Dirac spinors the corresponding current contraction vanishes. The quantum extension is the Ward identity, which makes physical amplitudes insensitive to adding a multiple of the photon momentum to its polarization vector. A gauge-compatible regularization preserves this vector-current identity.

With ${\mathcal{L}}_{\mathrm{int}}=-e\bar{\psi }{\gamma }^a\psi A_a$, expand $\exp (i\int {\mathcal{L}}_{\mathrm{int}})$ in powers of $e$. A term of order $n$ contains $n$ spacetime integrations and $1/n!$. Applying the Wick theorem pairs the free fields: an $A$-$A$ Wick contraction supplies a photon propagator, and a $\psi$-$\bar{\psi }$ Wick contraction supplies an oriented Dirac propagator. Each insertion supplies an interaction vertex. The permutations of contractions cancel the expansion factorials except for the Feynman-diagram symmetry factor. Interchanging Grassmann fields produces the fermionic sign, including a minus sign for every closed fermion loop. This is how Feynman diagrams arise from expectation values, rather than an extra dynamical assumption.

The normalized vacuum generating functional removes components with no external insertions. In the operator formalism, for an interacting-vacuum expectation value of an inserted product $\mathcal{O}$, the same cancellation appears asVacuum projection and the Feynman i-epsilon prescription are implicit. The linked-cluster theorem exponentiates all connected vacuum bubbles into the same factor in numerator and denominator, so it cancels. Normalize by $Z[0]$ to remove every vacuum component. This cancellation of vacuum bubbles still leaves products of disconnected diagrams that each contain external insertions. If only connected correlators are wanted, differentiate the connected generating functional $W=-i\log [Z[J,\bar{\eta },\eta ]/Z[0]]$.