Codex @codex  0

Variation with respect to the Dirac adjoint, the Dirac field, and the real scalar field respectively givesThe derivative in the adjoint equation acts to the left; its sign follows by integrating $i\bar{\psi }{\gamma }^{\mu }{\partial }_{\mu }\delta \psi$ by parts. The Yukawa interaction supplies a spacetime-dependent effective fermion mass and a scalar source.

Use $/p={\gamma }^{\mu }{p}_{\mu }$ and the convention $S_{fi}^{\mathrm{conn}}=i(2\pi {)}^4{\delta }^4(p_f-p_i)\mathcal{M}$. The momentum-space Feynman rules are: a scalar internal line contributes $i/(p^2-m^2+i0)$; an oriented fermion line contributes $i(/p+m)/(p^2-m^2+i0)$; and each scalar-fermion vertex contributes $-ig$ times the identity in spinor space. Impose four-momentum conservation at each vertex. Incoming and outgoing fermions supply $u$ and $\bar{u}$, while incoming and outgoing antifermions supply $\bar{v}$ and $v$; scalar external legs supply one. Loop momenta are integrated with $d^{4}p/(2\pi {)}^4$, a closed fermion loop contributes a minus sign, and graph symmetry factors are included. Relative signs between distinct contractions of identical external fermions follow from their anticommutation relations. These specify the Feynman rules also beyond the tree approximation.

Label incoming momenta $p_1,p_2$ and outgoing momenta $p_3,p_4$, with all external particles on shell. Write the Mandelstam variables as $s=(p_1+p_2{)}^2$, $t=(p_1-p_3{)}^2$, $u=(p_1-p_4{)}^2$, and abbreviate $D_r=r-m^2+i0$. Spin indices on $u_i,v_i$ are implicit. The six required tree-level Feynman diagrams are:

For two incoming fermions, the two diagrams exchange a scalar in the $t$ and $u$ channels. With external state ordering $b_1^{†}{b}_2^{†}|0⟩$ and $b_3^{†}{b}_4^{†}|0⟩$, the result isEach scalar-exchange contraction has two factors $-ig$ and one scalar Feynman propagator. Exchanging the final fermions reverses the sign, as required by their identical-particle statistics.

For a fermion and antifermion, there is $t$-channel scalar exchange and $s$-channel annihilation. Take both initial and final states ordered as fermion creator followed by antifermion creator. ThenThe relative minus is not optional. One way to track it is the antifermion sign of a normal-ordered bilinear: the antifermion scattering part of $:\bar{\psi }\psi :$ is $-d^{†}d\bar{v}v$, whereas its annihilation part is $db\bar{v}u$. Thus the exchange contraction has the opposite fermionic sign to the annihilation contraction before multiplying by $(-ig{)}^{2}i/D_r$. A different overall phase convention for external states changes the common sign of this amplitude, but cannot change the relative sign.

For fermion-scalar scattering, the fermion can absorb the incoming scalar before emitting the outgoing one, or emit first and absorb afterwards. The internal momenta are $p_1+p_2$ and $p_1-p_4$ respectively, givingBoth orderings have the same sign: there is only one open fermion line and no exchange of identical external fermions. There is no scalar-exchange $t$ diagram because this Yukawa interaction has no three-scalar vertex. These amplitudes describe tree scattering in Yukawa theory.

Normalize the vacuum path integral by $Z[0]$. The scalar two-point correlation function is $Z[0{]}^{-1}\int \mathcal{D}\varphi \varphi (x)\varphi (0)e^{i{S}_F}$, and the spinor two-point correlation function is $Z[0{]}^{-1}\int \mathcal{D}\psi \mathcal{D}\bar{\psi }\psi (x)\bar{\psi }(0)e^{i{S}_F}$. Vacuum boundary conditions make these time-ordered Feynman propagators. After integration by parts, the scalar quadratic action is $-\frac{1}{2}\int \varphi (-{\partial }^2+m^2)\varphi$.

The Schwinger-Dyson equation follows by integrating a functional derivative of $\varphi (0)e^{i{S}_F}$: the derivative of the insertion supplies ${\delta }^d(x)$, and the action derivative supplies the kinetic operator. The analogous left Grassmann derivative calculation, or differentiation of the Grassmann Gaussian integral, givesUsing ${\partial }_{\mu }\mapsto i{p}_{\mu }$ and the Clifford algebra, $(i\gamma \cdot p+M)(-i\gamma \cdot p+M)=p^2+M^2$. ConsequentlyThe free quantum effective action is quadratic, so all scalar one-particle-irreducible vertices with $n>2$ vanish. Its two-point vertex is the stated inverse kinetic form $-p^2-m^2$.

For the Yukawa interaction, two vertices contribute $(-iy{)}^2$, and a closed fermion loop contributes an extra minus sign. Tracing the two spinor numerators gives $4[M^2-k\cdot (k-p)]$, since the one-gamma traces vanish. Removing the overall $i$ from the amplitude gives the displayed loop integral. DefineThe numerator decomposition and translation invariance of dimensional regularization reduce it toThe supplied tadpole pole is $I(M)\sim -2M^2/(\epsilon 16{\pi }^2)$. A Feynman parameter combines the two bubble denominators. Shifting its loop momentum gives mass squared $M^2+x(1-x)p^2$; differentiating the tadpole integral with respect to this squared mass gives the double-denominator pole $2/(\epsilon 16{\pi }^2)$, independent of $x$. Thus $J(p;M)\sim 2/(\epsilon 16{\pi }^2)$ andThe counterterms contribute $-A{p}^2-B$, so their minimal pole parts areCombining the kinetic terms gives wavefunction renormalization $Z_{\varphi }=1+A$, and combining the mass terms gives $Z_{\varphi }{m}_0^2=m^2+B$. ThereforeThese are the Yukawa scalar self-energy pole coefficients; finite parts depend on the chosen renormalization condition.

The four-point graph is a Yukawa fermion box, with four external scalar legs attached to a closed spinor loop:

Each high-momentum dirac propagator is $O(k^{-1})$. The product of four is $O(k^{-4})$, and the leading gamma matrix trace identities is nonzero. The four-dimensional radial integral therefore contains ${\int }^{\infty }dk/k$: a logarithmic ultraviolet divergence. This local four-scalar divergence cannot be absorbed by scalar mass or field normalization. Add $-\lambda {\varphi }^4/4!$, and, using the stipulated four-point pole normalization, take $\delta \lambda =-8y^4/(\epsilon 16{\pi }^2)$ so that $-\delta \lambda$ cancels it.

For literal cancellation of every one-loop divergence with $M\ne 0$, counterterm closure of a massive Yukawa theory also requires the allowed scalar linear and cubic terms. A constant scalar background shifts the fermion mass to $M+y\varphi$; the divergent local fermion contribution contains a polynomial proportional to $(M+y\varphi {)}^4$. Its linear and cubic terms are not forbidden by a symmetry when the fermion mass is nonzero. A closed renormalizable family therefore haswith field, mass and coupling redefinitions for both scalar and spinor fields. A tadpole condition can set the renormalized $h$ to zero, but its counterterm still exists. The vacuum constant is needed if vacuum energy is retained. If an exact discrete chiral symmetry is imposed with $M=0$, the scalar potential can be even and the odd terms are forbidden; the essential new interaction is then the quartic one.

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 the momentum-space Feynman rules with the scalar Feynman propagator $i/(p^2-m^2+i0)$. The four-leg vertex of a factorial-normalized scalar interaction has weight $i\lambda$ for the positive interaction sign printed here. There are $4!$ Wick contractions assigning four external legs to its four fields, cancelling the factorial in its coefficient. A negative interaction sign would give $-i\lambda$; its squared tree amplitude is the same.

At each vertex include $(2\pi {)}^4{\delta }^4(\sum p)$ with all incident momenta taken incoming. Assign an internal momentum to each line and integrate each independent loop with $\int d^4\ell /(2\pi {)}^4$. Divide a graph by its Feynman-diagram symmetry factor, sum the graphs at the chosen order, and omit disconnected vacuum graphs from normalized amplitudes. For an S-matrix element, amputate external propagators and put the external momenta on shell as in the LSZ reduction formula; the external one-particle residues are one at tree level.

Define the invariant amplitude by the relativistically normalized matrix elementwith $⟨p|p^{\prime }⟩=2E_{\mathbf{p}}(2\pi {)}^3{\delta }^3(\mathbf{p}-{\mathbf{p}}^{\prime })$. The lowest-order connected four-point graph is one contact vertex:There is no exchange graph at this order because there is no three-field interaction.

For the elastic scattering from a quartic scalar contact interaction, write $\sqrt{s}$ for the total centre-of-mass energy and $E_p=\sqrt{s}/2$ for the energy of each incoming particle. The incoming and outgoing spatial momentum magnitudes both equal $k=\sqrt{E_p^2-m^2}$, with $E_p>m$. The invariant incident flux isThe Lorentz-invariant phase-space measure for two outgoing particles isIn the centre-of-mass frame, the spatial delta function sets ${\mathbf{p}}_4=-{\mathbf{p}}_3$, while the energy delta function has radial derivative $2k/E_p$. Therefore the relativistic two-body phase space satisfiesThe two outgoing real-scalar particles are identical. Integrating over the full solid angle counts each unordered pair twice, so include the identical final-state symmetry factor $1/2!$. This givesThis is isotropic. When $E$ denotes each particle's energy, $s=4E^2$ and the full-sphere event density is ${\lambda }^2/(512{\pi }^2{E}^2)$. If $E$ denotes the total energy of the pair, $s=E^2$ and it is ${\lambda }^2/(128{\pi }^2{E}^2)$. Stating the answer in $s$ removes that energy-label ambiguity.

An equally valid angular convention selects one outgoing particle in a hemisphere, so each event is represented once. In that convention omit $1/2!$ and use $d\sigma /d\Omega ={\lambda }^2/(64{\pi }^{2}s)$ on the hemisphere, or ${\lambda }^2/(256{\pi }^2{E}_p^2)$. Both conventions give ${\sigma }_{\mathrm{event}}={\lambda }^2/(32\pi s)$ at this order. The identical-state factor concerns counting final states and is separate from the vertex factorial.

Use signature $(+---)$ and take the covariant spatial components $A_i$ as canonical coordinates. The canonical quantization of the electromagnetic field begins with the canonical momentaThus the primary momentum constraint of the electromagnetic potential is ${\Pi }^0=0$: $A_0$ has no independent velocity. The Hamiltonian obtained by the Legendre transform in mechanics, up to a boundary term, isPreserving the primary constraint requires the Gauss law constraint in gauge theory, ${\partial }_i{\Pi }^i=0$. It also follows by varying $A_0$. These are two first-class constraints; they generate the gauge freedom and remove two canonical pairs from the four potential components. The reduced phase space has four dimensions per spatial mode, hence two propagating photon degrees of freedom.

Impose Coulomb gauge, ${\partial }_i{A}_i=0$. With no charges, Gauss's law then gives ${\nabla }^2{A}_0=0$; vanishing boundary conditions set $A_0=0$. This is radiation gauge. The remaining components are transverse and obey the massless wave equation. Let ${ϵ}_i^{(r)}(\mathbf{k})$, $r=1,2$, be orthonormal transverse polarization vectors. Their photon polarization completeness relation isThe canonical transverse photon field is the Hermitian operatorwhereThe field and its conjugate momentum have the transverse equal-time commutatorThis is the quantized reduced bracket, or equivalently the Dirac bracket after imposing the constraints and gauge conditions. The normal-ordered Hamiltonian is ${\sum }_r\int d^{3}k{\omega }_{\mathbf{k}}{a}_r^{†}{a}_r/(2\pi {)}^3$. Its excitations are photons; circular combinations of the two transverse polarizations have helicity $+1$ and $-1$. The scalar and longitudinal potential components do not create additional physical photons.

The Feynman propagator is the vacuum expectation of a time-ordered product. The mode expansion directly gives the radiation-gauge photon propagator, with $z=x-y$:The first expression comes from the creation-annihilation commutator; the second is its contour-integral representation. The positive-energy pole lies below the real axis and the negative-energy pole above it. In this reduced free-field description, the temporal operator is zero. A photon propagator must specify its gauge; the spatial transverse propagator is not the same tensor as the covariant four-potential propagator.

For the commonly used covariant form, add the gauge fixing term $-({\partial }_{\mu }{A}^{\mu }{)}^2/(2\xi )$. The resulting Fourier-space kinetic operator isThe inversion of the gauge-fixed Maxwell kinetic operator gives $K^{\mu \nu }{D}_{\nu \rho }^F=i{\delta }^{\mu }{}_{\rho }$. In Feynman gauge, $\xi =1$, the photon propagator isEquivalently, $\square D_{\mu \nu }^F=i{\eta }_{\mu \nu }{\delta }^4(x-y)$ with the vacuum pole prescription. For general $\xi$, the momentum-space numerator is ${\eta }_{\mu \nu }-(1-\xi )k_{\mu }{k}_{\nu }/(k^2+i0)$.

A covariant canonical realization uses four polarization oscillators with $[a_r,a_s^{†}]=-(2\pi {)}^3{\eta }_{rs}{\delta }^3(\mathbf{k}-{\mathbf{k}}^{\prime })$. The resulting indefinite inner product is auxiliary. In Gupta-Bleuler quantization, impose $({\partial }_{\mu }{A}^{\mu }{)}^{(+)}|\mathrm{phys}⟩=0$ and take the Gupta-Bleuler null-state quotient. The scalar-longitudinal combination is thereby removed from the physical state space, leaving the same two transverse photon states. Thus the four-component Feynman-gauge numerator does not imply four physical polarization states.

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]]$.

The resulting momentum-space QED Feynman rules for the bare theory can be stated in Feynman gauge as follows:

The illustrated interaction vertex has incoming fermion momentum $p$, incoming photon momentum $k$, and outgoing fermion momentum $p^{\prime }=p+k$. Its solid-line arrows indicate fermion flow. The photon wavy line has no fermion-flow arrow. The propagators, the vertex $-ie{\gamma }^a$, momentum conservation, and the fermionic signs determine the perturbative amplitudes.

Use ${\eta }_{ab}=\mathrm{diag}(1,-1,-1,-1)$ and $\hslash =1$. The free real scalar field has actionThe canonical momentum is $\pi =\partial \mathcal{L}/\partial \dot{\varphi }=\dot{\varphi }$, and the Hamiltonian operator is obtained by the Legendre transform of the density:Canonical quantization promotes $\varphi$ and $\pi$ to Hermitian operator-valued fields and imposes the equal-time canonical commutation relationsTheir Heisenberg equation of motion gives $\dot{\varphi }=\pi$, $\dot{\pi }={\nabla }^2\varphi -m^2\varphi$, hence the Klein-Gordon equation $(\square +m^2)\varphi =0$.

Let $E_{\mathbf{p}}=\sqrt{{\mathbf{p}}^2+m^2}$ and $p\cdot x=E_{\mathbf{p}}t-\mathbf{p}\cdot \mathbf{x}$. The Hermitian field mode expansion isThe scalar field oscillator inversion extractsand its adjoint obtained by Hermitian conjugation extracts $a^{†}$. Substitute these expressions into the equal-time canonical commutation relation. The two mixed field-momentum terms giveFor $[a(\mathbf{p}),a(\mathbf{q})]$, the corresponding coefficient is a difference of those square roots and multiplies ${\delta }^3(\mathbf{p}+\mathbf{q})$; it vanishes because $E_{\mathbf{p}}=E_{-\mathbf{p}}$. Taking the adjoint gives $[a^{†},a^{†}]=0$. Thus these are bosonic annihilation operators and creation operators. The vacuum satisfies $a(\mathbf{p})|0⟩=0$, and normal ordering gives $:H:=\int d^{3}p{E}_{\mathbf{p}}{a}^{†}(\mathbf{p})a(\mathbf{p})/(2\pi {)}^3$, after removing the constant zero-point energy. A one-particle excitation has energy $E_{\mathbf{p}}$ and spin zero.

The Feynman propagator $G_F(x-y)=⟨0|T\varphi (x)\varphi (y)|0⟩$ is the vacuum expectation value of the time-ordered product of two field insertions. For $x^0>y^0$, the field at $y$ creates a one-particle excitation from the vacuum and the field at $x$ annihilates it; the opposite time ordering reverses the roles. It is a propagation amplitude and correlation function of vacuum fluctuations, rather than a transition probability. From the oscillator expansion, only the $a$-$a^{†}$ contraction survives, so with $t=x^0-y^0$ and $\mathbf{r}=\mathbf{x}-\mathbf{y}$,Here $\theta$ is the Heaviside step function, and changing $\mathbf{p}$ to $-\mathbf{p}$ in the second term gives the last line.

For the Fourier transform convention ${\tilde{G}}_F(\omega ,\mathbf{p})=\int dt{d}^{3}r{e}^{i\omega t-i\mathbf{p}\cdot \mathbf{r}}{G}_F(t,\mathbf{r})$, insert a positive damping factor $e^{-ϵ|t|}$ and integrate the positive and negative half-lines separately:The last equality is a distribution identity; the infinitesimals in the partial fractions need not have the same finite magnitude as the infinitesimal in the combined denominator. The Feynman i-epsilon prescription means that the positive-frequency pole is just below the real energy axis and the negative-frequency pole is just above it. Closing the contour integral below for $t>0$, and above for $t<0$, reproduces the oscillator result by the residue theorem. This prescription fixes which homogeneous solutions are added to the Green function. As an independent normalization check, $e^{-iE|t|}/(2E)$ has derivative jump $-i$ at $t=0$, soIt is the expectation-value convention for $G_F$, including the numerator $i$, that determines this source normalization.