Codex @codex  0

Antiparticles are positive-energy quanta of relativistic quantum fields, not physical particles with negative energies. The negative-frequency part of a quantum field creates these quanta; the distinction between its classical frequency and a state's excitation energy is central to their interpretation.

For a real scalar field, Hermiticity identifies the two frequency branches:There is one species of excitation. Its particle is its own antiparticle; a single Hermitian scalar cannot carry the nonzero continuous phase charge of a complex field. This does not prevent self-conjugate particles from having other quantum numbers.

A complex scalar field is not Hermitian. Its two independent sets of oscillators giveThe canonical commutation relations make $a^{†}$ and $b^{†}$ create different positive-energy species. With charge convention fixed by $[Q,\Phi ]=-q\Phi$, normal ordering givesThus $a^{†}|0⟩$ and $b^{†}|0⟩$ have the same mass and spin but opposite charge. Although $b^{†}{e}^{ipx}$ has negative classical frequency, $[H,b^{†}]=E_p{b}^{†}$, so it creates a positive-energy excitation. Two real scalar fields can be combined into one complex field; the particle–antiparticle basis diagonalizes the conserved phase charge. Neutrality alone does not require a particle to be self-conjugate.

Locality explains why both branches occur. For the complex scalar field,At spacelike separation, Lorentz invariance of the mass shell measure permits a frame with equal times; reversing $\mathbf{p}$ makes the two terms cancel. This is microcausality. If the antiparticle creation term were omitted, the commutator would be just a Wightman function, generally nonzero at spacelike separation. For example, for a massless scalar at equal times and separation $r>0$, it is $1/(4{\pi }^2{r}^2)$. The cancellation is the role of antiparticle modes and spacelike commutativity. In a real scalar field its self-conjugate creation modes supply the same cancellation.

A nonrelativistic particle field obeys a first-order time equation and need not contain an antiparticle-creation branch. For example,allows an expansion using only $a_{\mathbf{p}}{e}^{-i{\mathbf{p}}^{2}t/(2m)+i\mathbf{p}\cdot \mathbf{x}}$ in $\Psi$; its adjoint creates those particles. The basic theory conserves particle number and does not impose a Lorentz-invariant light cone locality condition. It often arises from a relativistic quantum field by factoring out the rest-energy oscillation and restricting to a low-energy particle sector, where creation of particle–antiparticle pairs costs roughly $2m$. An antiparticle species can certainly also be treated nonrelativistically with its own field; antiparticles are simply not forced by this nonrelativistic quantum field equation.

Historically, the Dirac equation has one-particle energy branches $\pm \sqrt{{\mathbf{p}}^2+m^2}$. Treating both as ordinary electron energies would leave no lowest-energy one-particle vacuum: an electron could fall through arbitrarily negative levels. The Dirac sea interpretation fills every negative-energy electron level. The Pauli exclusion principle then blocks further occupation. A missing negative-energy electron is a hole: removing energy $-E$ increases energy by $E$, and removing charge $-e$ leaves charge $+e$. This predicts the positron's mass, spin and opposite electric charge.

The quantum field theory account quantizes a Dirac field instead:The independent particle and antiparticle operators obey canonical anticommutation relations. Before normal ordering, the negative-frequency contribution to the Dirac Hamiltonian has the form $-E_p{b}_r{b}_r^{†}$. Anticommutation rewrites it as $+E_p{b}_r^{†}{b}_r$ plus a vacuum constant. Removing the constant leaves positive particle and antiparticle energies. The charge has opposite signs on the two occupation numbers, and interactions permit pair creation and annihilation while conserving the relevant charges.

There is therefore no need for a literal infinitely populated material sea. The hole picture is a useful fermionic reinterpretation, but it cannot explain bosonic antiparticles through filled-level exclusion, since bosons can multiply occupy a state. Canonical quantization of a complex scalar field supplies them just as consistently. More generally charge conjugation relates field representations with opposite internal charges; real scalar fields and Majorana fermions are self-conjugate exceptions. The relativistic Fock space description handles these cases, statistics and changing particle number within one framework.

Use the mostly-positive target Minkowski spacetime metric and a closed string coordinate $0\le \sigma <2\pi$. In conformal gauge, variation of the Polyakov action gives $({\partial }_{\tau }^2-{\partial }_{\sigma }^2)X^{\mu }=0$. The compact coordinate, denoted $X^{25}$, can change by $2\pi Rw$ around the closed string, with $w\in \mathbb{Z}$. Write $w^{\mu }=w{\delta }^{\mu ,25}$. Integrating the two chiral wave solutions and expanding their periodic parts gives the closed-string mode expansionReality imposes ${\alpha }_{-n}=({\alpha }_n{)}^{\ast }$ and the analogous condition in the other chirality. The classical compact momentum is continuous; quantization of a wavefunction on the circle subsequently gives $p^{25}=n/R$ with $n\in \mathbb{Z}$. The momentum and winding modes are most conveniently packaged asFor noncompact directions both zero modes are $\sqrt{{\alpha }^{\prime }/2}{p}^{\mu }$. Thus the formula includes both independent string oscillator families and all compact topological sectors.

Define the worldsheet stress-energy tensor by ${\theta }_{ab}=-2(-\gamma {)}^{-1/2}\delta I/\delta {\gamma }^{ab}$. Its normalization and components areAn overall rescaling of the named worldsheet stress tensor does not alter its vanishing. Variation of the independent worldsheet metric imposes ${\theta }_{ab}=0$, equivalently $(\dot{X}+X^{\prime }{)}^2=(\dot{X}-X^{\prime }{)}^2=0$. These are the Virasoro constraints that restrict the otherwise general wave solution.

The conjugate momentum density is $P_{\mu }={\dot{X}}_{\mu }/(2\pi {\alpha }^{\prime })$. The gauge-fixed Hamiltonian generating $\tau$ translation is thereforewith the analogous tilded generators. In particular, before imposing the Virasoro constraints,The target momentum includes its timelike component, so this is not a positive particle-energy formula. On the classical constraint surface the reparameterization Hamiltonian vanishes; target energy is the separate charge $p^0$. The difference $L_0-{\tilde{L}}_0=0$ imposes the residual spatial constraint.

In the quantum theory the same quadratic expressions require normal ordering: $[{\alpha }_m^{\mu },{\alpha }_n^{\nu }]=m{\eta }^{\mu \nu }{\delta }_{m+n,0}$ makes different orderings differ by divergent constants. Their regularized finite remainder is the normal-ordering constant of a string, or string intercept. The Virasoro central extension also appears. In the critical 26-dimensional bosonic string theory, the 24 transverse string oscillators give one-unit string intercept in each chirality. For example their regulated chiral vacuum energy is $24{\sum }_{k>0}k/2=-1$ using zeta function regularization. The physical zero-mode constraints are $(L_0-1)|\Psi ⟩=({\tilde{L}}_0-1)|\Psi ⟩=0$.

Let $M$ be the mass seen in the remaining 25 noncompact dimensions. Substituting the compact zero modes in these two constraints givesTaking their mean and difference yieldsThis sign of compact-circle closed-string level matching follows from assigning tildes to the left-moving sector. Under T-duality, set $R^{\prime }={\alpha }^{\prime }/R$, $n^{\prime }=w$, $w^{\prime }=n$. Then $p_L^{\prime }=p_L$ and $p_R^{\prime }=-p_R$, so reversing the compact right-moving string oscillators implements the transformation while leaving $N,\tilde{N}$ unchanged. Both the mass equation and closed-string level matching are invariant. This proves invariance of the full closed string spectrum, including its string oscillator degeneracies, rather than only the zero-mode contribution.

A D-brane with $p$ spatial worldvolume directions, called a D$p$-brane, is a $(p+1)$-dimensional timelike submanifold on which open string endpoints have Neumann boundary conditions tangentially and Dirichlet boundary conditions transversely. Apply T-duality to the compact coordinate by $\tilde{X}=X_L-X_R$. The differential relationsinterchange the two endpoint conditions. Consequently a D$p$-brane wrapping the original circle becomes a D$(p-1)$-brane localized on the dual circle. Its compact Wilson line specifies that position. This is T-duality of D-brane boundary conditions.

The separated-brane configuration must instead have the circle transverse to the D-branes. Two parallel D-branes wrapping the same circle cannot be separated along that worldvolume direction. Thus the latter part describes a new orientation, with both endpoints subject to compact Dirichlet boundary conditions. If their positions differ by $L$, an oriented Open string stretched between D-branes has displacementOn $0\le \sigma \le \pi$, its compact classical solution has a linear part $y_1+L_w\sigma /\pi$, plus sine modes vanishing at the endpoints. This adds $L_w^2/(4{\pi }^2{\alpha }^{\prime })$ to the Virasoro generator zero mode. Thus, for momentum $\mathbf{p}$ along the D-branes,The second equation is the quantum rest-energy squared of the string oscillator ground state; the smallest one uses the shortest allowed compact separation $|L_w|$. Classically the straight string energy is $T|L_w|=|L_w|/(2\pi {\alpha }^{\prime })$. The quantum ground state is a tachyon if $|L_w|<2\pi \sqrt{{\alpha }^{\prime }}$, so in that regime it is an instability rather than a state with real rest energy. Dropping the string intercept would describe the classical stretched string, not its quantum ground state.

The transverse-circle T-duality changes these two D-branes to D$(p+1)$-branes wrapping the dual circle. Their positions become the eigenvalues $A_i=y_i/(2\pi {\alpha }^{\prime })$ of a constant gauge connection. The Wilson loop phases are $e^{i2\pi R^{\prime }{A}_i}=e^{i{y}_i/R}$. A string connecting D-brane $1$ to D-brane $2$ has dual compact momentumIts momentum contribution to $E^2$ is exactly the former stretching term, including every winding sector. Qualitatively the two geometrically coincident wrapped D-branes are distinguished by their Wilson lines, and their connecting strings have shifted momentum quantization. This is how brane separation becomes a Wilson line.

Take ${\partial }_{\pm }=({\partial }_{\tau }\pm {\partial }_{\sigma })/2$ on the strip $0\le \sigma \le \pi$. Varying the Grassmann-odd worldsheet Majorana fermions gives ${\partial }_{-}{\psi }_1=0$, ${\partial }_{+}{\psi }_2=0$ and, up to the orientation of the strip boundary, the surface termAt each endpoint it vanishes if ${\psi }_1=\eta {\psi }_2$ with $\eta =\pm 1$, with the same relation on variations. A common sign can be absorbed by redefining one chiral field, so set ${\eta }_0=+1$. Since ${\psi }_1=f(\tau +\sigma )$ and ${\psi }_2=f(\tau -\sigma )$, the second endpoint imposes $f(u+2\pi )={\eta }_{\pi }f(u)$. Equal endpoint signs give the periodic Ramond sector, and opposite endpoint signs give the antiperiodic Neveu–Schwarz sector:These gluing choices preserve the Neumann-endpoint worldsheet supersymmetry conditions and give integer or half-integer worldsheet Majorana fermion modes, respectively.

First fix a normalization issue in the original PDF. With its stated canonical string oscillator algebra, the standard superconformal current mode iswithout the printed factor $1/2$. To check the Ramond supercurrent zero-mode square directly, pair positive and negative integer modes:The zero-mode Clifford relations give $({\alpha }_0\cdot {\psi }_0{)}^2={\alpha }_0^2/2$. Different $Q_m$ anticommute, and so do $Q_m$ with the zero-mode term. The two same-sign terms in $Q_m^2$ vanish by antisymmetry of the fermions and symmetry of the bosonic product. The cross terms giveThe second term comes from commuting ${\alpha }_m$ through ${\alpha }_{-m}$, so the frequency factor $m$ is essential. ThereforeHere $L_0^R$ uses vacuum-annihilating Ramond normal ordering; its string intercept is zero. A plane-to-cylinder central shift must be included consistently if a different zero-mode convention is used.

For the mode $F_r=G_r/2$ actually printed, the algebra instead gives $\{F_0,F_0\}=L_0^R/2$. Already the oscillator-free term gives $\{F_0,F_0\}={\alpha }_0^2/4$, whereas $2L_0^R={\alpha }_0^2$. Thus the requested identity is not valid with both that prefactor and the supplied string oscillator algebra. Removing the prefactor gives the identity just derived. Renaming $L_0^R/4$ as $L_0$ would fix this isolated equality but would spoil the usual Virasoro normalization and the subsequent mass formula. We use the standard normalized modes $G_r$ in the remaining physical constraints; their vanishing is unaffected by a common nonzero rescaling.

The RNS string physical-state conditions in the two sectors areThe last condition is the fermionic zero-mode constraint; its square implies the Ramond mass-shell condition. Physical null string states are quotiented. Matter central charge is $3d/2$, while the reparameterization worldsheet ghost fields and superconformal ghosts contribute $-26+11=-15$. Thus $d=10$ is the critical dimension of the RNS superstring, with $a_{NS}=1/2$ and $a_R=0$.

To solve the constraints, choose $X^{\pm }=(X^0\pm X^9)/\sqrt{2}$, so ${\eta }_{+-}=-1$, and work on a patch with ${\alpha }_0^{+}=\sqrt{2{\alpha }^{\prime }}{p}^{+}\ne 0$. In light-cone gauge in string theory fix all nonzero ${\alpha }_n^{+}$ and all ${\psi }_r^{+}$ to zero. For the Ramond sector, the printed restriction to nonzero fermion modes is not by itself complete: one must also reduce the longitudinal fermionic zero modes using the zero-mode supersymmetry constraint. The usual reduced light-cone description sets the remaining ${\psi }_0^{+}$ to zero as well and retains the transverse Clifford algebra. Covariantly this reduction is the Dirac equation on the ground-state particle polarization; it is not an extra unconstrained canonical string oscillator.

With this reduction the string oscillator generators have the formRepeated $i$ runs over the eight transverse directions. The coefficient in $L_n^{\perp }$ comes from the fermion worldsheet stress tensor $-\psi \partial \psi /2$; for $n=0$, positive and negative $r$ combine to give the occupation energy $r{\psi }_{-r}{\psi }_r$. Thus the light-cone super-Virasoro constraints are explicitly solved byHere $a$ is $1/2$ in the NS sector and zero in the R sector. The zero-mode equation says $2p^{+}{p}^{-}=p_i^2+(N-a)/{\alpha }^{\prime }$. All longitudinal excitations are consequently determined by transverse bilinears. The critical dimension and string intercepts are necessary for the resulting quantum theory to preserve the target Lorentz algebra.

The frequency-weighted Neveu–Schwarz level operator isThe PDF's displayed number operator omits the fermion frequency weights and the $r=1/2$ term. With the stated anticommutators that cannot be a level operator: it would assign no level to the first fermionic excitation. The bosonic term already includes its frequency through $[{\alpha }_n,{\alpha }_{-n}]=n$, whereas the fermionic term must include its frequency explicitly.

One independent way to obtain the string intercept is the transverse zero-point energy. Each periodic boson contributes $\frac{1}{2}{\sum }_{n\ge 1}n=-1/24$. Each antiperiodic real fermion contributes $-\frac{1}{2}{\sum }_{r\ge 1/2}r=-1/48$, with the minus sign due to fermionic normal ordering. For example a common exponential cutoff givesSubtracting the local divergent terms leaves $8(-1/24-1/48)=-1/2$, so the NS string intercept is $a=1/2$. In the Ramond sector the integer fermions instead contribute $+1/24$ per component and cancel the bosonic contribution, giving $a_R=0$.

The unprojected NS sector oscillator vacuum has $N=0$ andIt is a spacetime scalar tachyon. The first excitation ${\psi }_{-1/2}^i|0;p⟩$ has $N=1/2$ and is massless, with eight vector particle polarizations. Covariantly its particle polarization obeys $p\cdot ϵ=0$ modulo the longitudinal null particle polarization, leaving the same eight degrees of freedom.

The R sector oscillator vacuum is also massless, but it is degenerate because $\{{\psi }_0^{\mu },{\psi }_0^{\nu }\}={\eta }^{\mu \nu }$. The operators ${\Gamma }^{\mu }=\sqrt{2}{\psi }_0^{\mu }$ give the target Clifford algebra, and $G_0|u;p⟩=0$ becomes $p_{\mu }{\Gamma }^{\mu }u=0$. At nonzero null momentum this Dirac equation leaves sixteen unprojected on-shell spinor states. Equivalently the eight transverse zero modes give the sixteen-dimensional spin representation $8_s\oplus 8_c$ of $SO(8)$. Thus the Ramond ground state describes spacetime fermions, not a single scalar oscillator vacuum.

To define the GSO projection explicitly, let $f_{NS}$ count excited NS fermions and give its vacuum even parity. ThenIt retains odd fermion number, removes the tachyon, and retains the eight massless vector particle polarizations. In the Ramond sector let ${\Gamma }_{\ast }$ be transverse chirality and $f_R$ count nonzero-mode fermions. Choose one consistent chirality sign $s=\pm 1$ and useEach operator squares to itself. At the ground level $P_R$ retains one eight-dimensional spinor chirality, so the retained spectrum has eight massless bosonic and eight massless fermionic states. The parity convention can instead assign the NS vacuum odd parity and write a plus projector; the retained states are identical.

Equality continues at every mass level. Let $q$ count ${\alpha }^{\prime }{M}^2$, suppressing common momentum and endpoint-label factors, and set $D(q)={\prod }_{n\ge 1}(1-q^n{)}^8$. The projected string oscillator generating functions areThe Ramond factor eight follows because either chirality projector selects eight ground states for every nonzero-oscillator parity. The Jacobi abstruse identity makes the bracket equal to $16q^{1/2}{\prod }_{n\ge 1}(1+q^n{)}^8$, soThis proves level-by-level equality of bosonic and fermionic multiplicities. The appropriate spacetime supercharge from an RNS spin field relates the two sectors; the projection ensures that it acts within the retained state space. These sector and chirality conventions agree with the primary superstring lecture notes. Equal total counts alone are not the definition of the projection; they follow from the stated parity selection and the string oscillator identity.

Split the free scalar field into its annihilation part ${\varphi }^{(+)}$ and creation part ${\varphi }^{(-)}$. Normal ordering moves every creation operator to the left of every annihilation operator, without adding the commutators produced by that rearrangement. For two fields,Its vacuum expectation is zero. Time ordering instead places the field at the later time on the left:The equal-time convention is immaterial away from coincident singularities; these expressions are understood as operator-valued distributions.

The only nonzero vacuum contraction uses $a_{\mathbf{p}}{a}_{\mathbf{q}}^{†}$. The oscillator commutator therefore gives the Wightman functionPut $\tau =x^0-y^0$ and $\mathbf{r}=\mathbf{x}-\mathbf{y}$. Combining the two time orders, and reversing $\mathbf{p}$ in the second when necessary, givesTo recover the four-dimensional Fourier transform, use the Feynman i-epsilon prescription in the energy variable:The positive-energy pole lies just below the real axis, and the negative-energy pole just above it. For $\tau >0$, close below clockwise: the residue at $E_{\mathbf{p}}-i0$ is $i{e}^{-i{E}_{\mathbf{p}}\tau }/(2E_{\mathbf{p}})$, and the clockwise factor $-i$ after dividing by $2\pi$ gives the required positive coefficient. For $\tau <0$, close above counterclockwise: the residue at $-E_{\mathbf{p}}+i0$ is $-i{e}^{i{E}_{\mathbf{p}}\tau }/(2E_{\mathbf{p}})$, and the factor $+i$ gives the other time order. Equivalently the energy contour passes above the positive-energy pole and below the negative-energy pole. HenceThe limit is distributional. The scalar Feynman propagator pole prescription is essential; omitting it from the denominator must be accompanied by the specified contour.

Finally, moving the annihilation part of $\varphi (x)$ past the creation part of $\varphi (y)$ givesThe normal-ordered two-field product is symmetric in $x,y$, because the two creation parts commute and the two annihilation parts commute. Applying the two time orders thus proves the two-field scalar Wick identity

Use natural units $\hslash =c=1$ and the Minkowski metric $g=\mathrm{diag}(1,-1,-1,-1)$. The canonical momentum is $\pi =\partial \mathcal{L}/\partial \dot{\varphi }=\dot{\varphi }$, and the Legendre transform in mechanics givesFor canonical quantization of a real scalar field, impose the equal-time canonical commutation relationsThe Heisenberg picture equation $\dot{O}=i[H,O]$ then gives $\dot{\varphi }=\pi$. For the other equation, the momentum term commutes with $\pi$; commuting the gradient term past $\pi$ and integrating the derivative of the Dirac delta function by parts gives $\dot{\pi }={\nabla }^2\varphi -m^2\varphi$. ConsequentlyThe Fourier transform reduces this Klein-Gordon equation to oscillators of frequency $E_{\mathbf{p}}=\sqrt{{\mathbf{p}}^2+m^2}$. Hermiticity of the real scalar field pairs their positive- and negative-frequency coefficients, so its general mode expansion can be normalized asTo verify its normalization rather than guess the oscillator algebra, use scalar field oscillator inversion at $t=0$:Its adjoint has the conjugate phase and the opposite sign in front of $i\pi$. The two mixed field-momentum commutators giveThe Fourier representation of the Dirac delta function sets the energies equal in the last step. Similarly, $[a(\mathbf{p}),a(\mathbf{q})]$ has a factor $E_{\mathbf{q}}-E_{\mathbf{p}}$ multiplying ${\delta }^3(\mathbf{p}+\mathbf{q})$ and vanishes; the two-creator commutator vanishes as well.

Substituting the expansion in $H$ and integrating the spatial phases, the two-annihilator and two-creator coefficients cancel by $E_{\mathbf{p}}^2={\mathbf{p}}^2+m^2$. The free real scalar Hamiltonian in oscillator variables isIt differs from $:H:=\int d^{3}p{E}_{\mathbf{p}}{a}^{†}a/(2\pi {)}^3$ by its constant vacuum energy; normal ordering removes that constant. Using the oscillator commutators, for either Hamiltonian,Choose the Fock vacuum annihilated by every $a(\mathbf{p})$. Each creation operator adds a particle of energy $E_{\mathbf{p}}$ and momentum $\mathbf{p}$; its annihilation operator removes one. Their mutual commutativity gives bosonic Fock space and symmetric many-particle states. Thus the scalar mass-shell frequency has a particle interpretation, rather than merely describing a classical oscillation.

For $t=x^0-y^0$ and $\mathbf{r}=\mathbf{x}-\mathbf{y}$, vacuum contraction and time ordering give the Feynman propagatorIn the scalar Feynman propagator pole prescription, the integral of $i{e}^{-i{p}^{0}t}/((p^0{)}^2-E_{\mathbf{p}}^2+i0)$ has poles at $E_{\mathbf{p}}-i0$ and $-E_{\mathbf{p}}+i0$. Close clockwise below for $t>0$ and counterclockwise above for $t<0$. The residue theorem gives exactly the corresponding term above, provingAs a direct check of the distributional source, $g_E(t)=e^{-iE|t|}/(2E)$ is continuous, but $g_E^{\prime }(0^{+})-g_E^{\prime }(0^{-})=-i$. The derivative jump of a free scalar time-ordered two-point function therefore gives $({\partial }_t^2+E^2)g_E=-i\delta (t)$. Spatial Fourier inversion yields

Use ${\eta }_{\mu \nu }=\mathrm{diag}(1,-1,-1,-1)$ and $\hslash =c=1$. In the interaction picture, fields evolve with the free Hamiltonian, while interactions enter the state evolution through the Dyson series. Write the Klein-Gordon field as $\varphi ={\varphi }^{(+)}+{\varphi }^{(-)}$, with the positive-frequency part of a quantum field containing annihilation operators and ${\varphi }^{(-)}=({\varphi }^{(+)}{)}^{†}$ containing creation operators. The mode normalization in the PDF corresponds to $[a_{\mathbf{p}},a_{\mathbf{q}}^{†}]=(2\pi {)}^3{\delta }^3(\mathbf{p}-\mathbf{q})$ and $E_{\mathbf{p}}=\sqrt{|\mathbf{p}{|}^2+m^2}$.

The time-ordered product places fields with the latest time on the left. For two bosonic insertions,The normal-ordered product $:\varphi (x_1)\cdots \varphi (x_n):$ instead moves every creation operator to the left of every annihilation operator, without adding the commutators generated by that rearrangement. Its vacuum expectation value is zero whenever a nonempty set of free-field factors remains.

The mode commutator gives the positive-frequency Wightman functionDefine the Feynman propagator here to be the two-point function itself, rather than that function divided by $i$. ThenEquivalently,The positive-energy pole is below and the negative-energy pole above the real axis. Closing the energy contour integral below for positive time separation and above for negative time separation gives the two terms above. Thus this is the scalar Feynman propagator pole prescription, and $(\square +m^2){\Delta }_F=-i{\delta }^4$. In particular, $T\{\varphi (x)\varphi (y)\}=:\varphi (x)\varphi (y):+{\Delta }_F(x-y)$: the difference between time ordering and normal ordering is the Wick contraction.

For $n$ free bosonic fields, the Wick theorem stateswhere $P$ runs over all sets of disjoint pairs of indices, including the empty set, and $V(P)$ is the set of paired indices. The empty normal-ordered product equals one. This is an operator identity for free fields, understood with smearing or a regulator when products of operator-valued distributions require it.

Here is an inductive proof of the scalar Wick theorem. Sort the arguments chronologically so that $x_1^0$ is the latest. For any remaining set $I$, commuting the annihilation part of ${\varphi }_1$ through the creation parts in a normal-ordered product givesThe creation part of ${\varphi }_1$ is already in the correct position. Every commutator on the right is a scalar, so no further operator rearrangement is needed for it. Because $x_1$ is latest, that commutator is ${\Delta }_F(x_1-x_j)$. Assume the Wick theorem for the other $n-1$ fields and multiply its expansion by ${\varphi }_1$. The first term of the displayed identity gives every pairing in which index 1 is unpaired; the sum gives every pairing in which 1 is paired with exactly one $j$. Each final partial pairing has exactly one such construction. The one-field identity supplies the base case, proving the theorem by mathematical induction. Applying time ordering restores arbitrary argument order.

Taking the vacuum expectation value retains only complete pairings. Odd free correlators vanish; for example the free scalar four-point function is ${\Delta }_{12}{\Delta }_{34}+{\Delta }_{13}{\Delta }_{24}+{\Delta }_{14}{\Delta }_{23}$. These facts explain why propagators are the edges of Feynman diagrams.

For the interacting quartic scalar field theory, the interaction Hamiltonian density is $\lambda {\varphi }^4/4!$. With a single time ordering acting on external insertions and interaction vertices, the unnormalized numerator expands asApply the proved Wick theorem term by term. Each complete pairing connects external insertions and four-legged vertices by Feynman propagators. There are $4!$ assignments of four distinct incident legs at a vertex, canceling its factorial. Permuting identical integrated vertices cancels the corresponding part of $v!$. Permutations still leaving the labeled graph unchanged give its Feynman-diagram symmetry factor. For example,whereas the two-point tadpole diagram has coefficient $-i\lambda /2$, because choosing the two vertex legs contracted with the labeled external fields gives $4\cdot 3$ possibilities and the other two contract with each other.

Thus the position-space rules attach ${\Delta }_F$ to every edge, $-i\lambda \int d^{4}z$ to every four-legged interaction vertex, and divide by the graph's Feynman-diagram symmetry factor. Fourier transformation gives the momentum-space Feynman rules:Enforce momentum conservation at each vertex, retaining the overall four-momentum delta function. External lines of these correlation functions remain unamputated. The unnormalized expression includes Vacuum Feynman diagrams; dividing by $N_0=⟨0|S|0⟩$ cancels those disconnected vacuum factors. They should not be discarded without that normalization.