Wick's theorem

Quantum field theory
Feynman diagram
History
Background Field theory Electromagnetism Weak force Strong force Quantum mechanics Special relativity General relativity Gauge theory Yang–Mills theory
Symmetries Symmetry in quantum mechanics C-symmetry P-symmetry T-symmetry Lorentz symmetry Poincaré symmetry Gauge symmetry Explicit symmetry breaking Spontaneous symmetry breaking Noether charge Topological charge
Tools Anomaly Background field method BRST quantization Correlation function Crossing Effective action Effective field theory Expectation value Feynman diagram Lattice field theory LSZ reduction formula Partition function Path Integral Formulation Propagator Quantization Regularization Renormalization Vacuum state Wightman axioms
Equations Dirac equation Klein–Gordon equation Proca equations Wheeler–DeWitt equation Bargmann–Wigner equations Schwinger-Dyson equation Renormalization group equation
Standard Model Quantum electrodynamics Electroweak interaction Quantum chromodynamics Higgs mechanism
Incomplete theories String theory Supersymmetry Technicolor Theory of everything Quantum gravity
Scientists Adler Anderson Anselm Bargmann Becchi Belavin Bell Berezin Bethe Bjorken Bleuer Bogoliubov Brodsky Brout Buchholz Cachazo Callan Cardy Coleman Connes Dashen DeWitt Dirac Doplicher Dyson Englert Faddeev Fadin Fayet Fermi Feynman Fierz Fock Frampton Fritzsch Fröhlich Fredenhagen Furry Glashow Gell-Mann Glimm Goldstone Gribov Gross Gupta Guralnik Haag Hagen Han Heisenberg Hepp Higgs 't Hooft Iliopoulos Ivanenko Jackiw Jaffe Jona-Lasinio Jordan Jost Källén Kendall Kinoshita Kim Klebanov Kontsevich Kreimer Kuraev Landau Lee Lee Lehmann Leutwyler Lipatov Łopuszański Low Lüders Maiani Majorana Maldacena Matsubara Migdal Mills Møller Naimark Nambu Neveu Nishijima Oehme Oppenheimer Osborn Osterwalder Parisi Pauli Peccei Peskin Plefka Polchinski Polyakov Pomeranchuk Popov Proca Quinn Rouet Rubakov Ruelle Sakurai Salam Schrader Schwarz Schwinger Segal Seiberg Semenoff Shifman Shirkov Skyrme Sommerfield Stora Stueckelberg Sudarshan Symanzik Takahashi Thirring Tomonaga Tyutin Vainshtein Veltman Veneziano Virasoro Ward Weinberg Weisskopf Wentzel Wess Wetterich Weyl Wick Wightman Wigner Wilczek Wilson Witten Yang Yukawa Zamolodchikov Zamolodchikov Zee Zimmermann Zinn-Justin Zuber Zumino

Wick's theorem is a method of reducing high-order derivatives to a combinatorics problem.^[1] It is named after Italian physicist Gian Carlo Wick.^[2] It is used extensively in quantum field theory to reduce arbitrary products of creation and annihilation operators to sums of products of pairs of these operators. This allows for the use of Green's function methods, and consequently the use of Feynman diagrams in the field under study. A more general idea in probability theory is Isserlis' theorem.

In perturbative quantum field theory, Wick's theorem is used to quickly rewrite each time ordered summand in the Dyson series as a sum of normal ordered terms. In the limit of asymptotically free ingoing and outgoing states, these terms correspond to Feynman diagrams.

For two operators ${\displaystyle {\hat {A}}}$ and ${\displaystyle {\hat {B}}}$ we define their contraction to be

${\displaystyle {\hat {A}}^{\bullet }\,{\hat {B}}^{\bullet }\equiv {\hat {A}}\,{\hat {B}}\,-{\mathopen {:}}{\hat {A}}\,{\hat {B}}{\mathclose {:}}}$

where ${\displaystyle {\mathopen {:}}{\hat {O}}{\mathclose {:}}}$ denotes the normal order of an operator ${\displaystyle {\hat {O}}}$. Alternatively, contractions can be denoted by a line joining ${\displaystyle {\hat {A}}}$ and ${\displaystyle {\hat {B}}}$, like ${\displaystyle {\overset {\sqcap }{{\hat {A}}{\hat {B}}}}}$.

We shall look in detail at four special cases where ${\displaystyle {\hat {A}}}$ and ${\displaystyle {\hat {B}}}$ are equal to creation and annihilation operators. For ${\displaystyle N}$ bosonic or fermionic modes we'll denote the creation operators by ${\displaystyle {\hat {a}}_{i}^{\dagger }}$ and the annihilation operators by ${\displaystyle {\hat {a}}_{i}}$ ${\displaystyle (i=1,2,3,\ldots ,N)}$. They satisfy the commutation relations for bosonic operators ${\displaystyle [{\hat {a}}_{i},{\hat {a}}_{j}^{\dagger }]=\delta _{ij}{\hat {\mathbf {1} }}}$, or the anti-commutation relations for fermionic operators ${\displaystyle \{{\hat {a}}_{i},{\hat {a}}_{j}^{\dagger }\}=\delta _{ij}{\hat {\mathbf {1} }}}$ where ${\displaystyle \delta _{ij}}$ denotes the Kronecker delta and ${\displaystyle {\hat {\mathbf {1} }}}$ denotes the identity operator.

We then have

${\displaystyle {\hat {a}}_{i}^{\bullet }\,{\hat {a}}_{j}^{\bullet }={\hat {a}}_{i}\,{\hat {a}}_{j}\,-{\mathopen {:}}\,{\hat {a}}_{i}\,{\hat {a}}_{j}\,{\mathclose {:}}\,=0}$

${\displaystyle {\hat {a}}_{i}^{\dagger \bullet }\,{\hat {a}}_{j}^{\dagger \bullet }={\hat {a}}_{i}^{\dagger }\,{\hat {a}}_{j}^{\dagger }\,-\,{\mathopen {:}}{\hat {a}}_{i}^{\dagger }\,{\hat {a}}_{j}^{\dagger }\,{\mathclose {:}}\,=0}$

${\displaystyle {\hat {a}}_{i}^{\dagger \bullet }\,{\hat {a}}_{j}^{\bullet }={\hat {a}}_{i}^{\dagger }\,{\hat {a}}_{j}\,-{\mathopen {:}}\,{\hat {a}}_{i}^{\dagger }\,{\hat {a}}_{j}\,{\mathclose {:}}\,=0}$

${\displaystyle {\hat {a}}_{i}^{\bullet }\,{\hat {a}}_{j}^{\dagger \bullet }={\hat {a}}_{i}\,{\hat {a}}_{j}^{\dagger }\,-{\mathopen {:}}\,{\hat {a}}_{i}\,{\hat {a}}_{j}^{\dagger }\,{\mathclose {:}}\,=\delta _{ij}{\hat {\mathbf {1} }}}$

where ${\displaystyle i,j=1,\ldots ,N}$.

These relationships hold true for bosonic operators or fermionic operators because of the way normal ordering is defined.

We can use contractions and normal ordering to express any product of creation and annihilation operators as a sum of normal ordered terms. This is the basis of Wick's theorem. Before stating the theorem fully we shall look at some examples.

Suppose ${\displaystyle {\hat {a}}_{i}}$ and ${\displaystyle {\hat {a}}_{i}^{\dagger }}$ are bosonic operators satisfying the commutation relations:

${\displaystyle \left[{\hat {a}}_{i}^{\dagger },{\hat {a}}_{j}^{\dagger }\right]=0}$

${\displaystyle \left[{\hat {a}}_{i},{\hat {a}}_{j}\right]=0}$

${\displaystyle \left[{\hat {a}}_{i},{\hat {a}}_{j}^{\dagger }\right]=\delta _{ij}{\hat {\mathbf {1} }}}$

where ${\displaystyle i,j=1,\ldots ,N}$, ${\displaystyle \left[{\hat {A}},{\hat {B}}\right]\equiv {\hat {A}}{\hat {B}}-{\hat {B}}{\hat {A}}}$ denotes the commutator, and ${\displaystyle \delta _{ij}}$ is the Kronecker delta.

We can use these relations, and the above definition of contraction, to express products of ${\displaystyle {\hat {a}}_{i}}$ and ${\displaystyle {\hat {a}}_{i}^{\dagger }}$ in other ways.

${\displaystyle {\hat {a}}_{i}\,{\hat {a}}_{j}^{\dagger }={\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{i}+\delta _{ij}={\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{i}+{\hat {a}}_{i}^{\bullet }\,{\hat {a}}_{j}^{\dagger \bullet }=\,{\mathopen {:}}\,{\hat {a}}_{i}\,{\hat {a}}_{j}^{\dagger }\,{\mathclose {:}}+{\hat {a}}_{i}^{\bullet }\,{\hat {a}}_{j}^{\dagger \bullet }}$

Note that we have not changed ${\displaystyle {\hat {a}}_{i}\,{\hat {a}}_{j}^{\dagger }}$ but merely re-expressed it in another form as ${\displaystyle \,{\mathopen {:}}\,{\hat {a}}_{i}\,{\hat {a}}_{j}^{\dagger }\,{\mathclose {:}}+{\hat {a}}_{i}^{\bullet }\,{\hat {a}}_{j}^{\dagger \bullet }}$

${\displaystyle {\hat {a}}_{i}\,{\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{k}=({\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{i}+\delta _{ij}){\hat {a}}_{k}={\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{i}\,{\hat {a}}_{k}+\delta _{ij}{\hat {a}}_{k}={\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{i}\,{\hat {a}}_{k}+{\hat {a}}_{i}^{\bullet }\,{\hat {a}}_{j}^{\dagger \bullet }{\hat {a}}_{k}=\,{\mathopen {:}}\,{\hat {a}}_{i}\,{\hat {a}}_{j}^{\dagger }{\hat {a}}_{k}\,{\mathclose {:}}+{\mathclose {:}}\,{\hat {a}}_{i}^{\bullet }\,{\hat {a}}_{j}^{\dagger \bullet }\,{\hat {a}}_{k}{\mathclose {:}}}$

${\displaystyle {\hat {a}}_{i}\,{\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{k}\,{\hat {a}}_{l}^{\dagger }=({\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{i}+\delta _{ij})({\hat {a}}_{l}^{\dagger }\,{\hat {a}}_{k}+\delta _{kl})}$

${\displaystyle ={\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{i}\,{\hat {a}}_{l}^{\dagger }\,{\hat {a}}_{k}+\delta _{kl}{\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{i}+\delta _{ij}{\hat {a}}_{l}^{\dagger }{\hat {a}}_{k}+\delta _{ij}\delta _{kl}}$

${\displaystyle ={\hat {a}}_{j}^{\dagger }({\hat {a}}_{l}^{\dagger }\,{\hat {a}}_{i}+\delta _{il}){\hat {a}}_{k}+\delta _{kl}{\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{i}+\delta _{ij}{\hat {a}}_{l}^{\dagger }{\hat {a}}_{k}+\delta _{ij}\delta _{kl}}$

${\displaystyle ={\hat {a}}_{j}^{\dagger }{\hat {a}}_{l}^{\dagger }\,{\hat {a}}_{i}{\hat {a}}_{k}+\delta _{il}{\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{k}+\delta _{kl}{\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{i}+\delta _{ij}{\hat {a}}_{l}^{\dagger }{\hat {a}}_{k}+\delta _{ij}\delta _{kl}}$

${\displaystyle =\,{\mathopen {:}}{\hat {a}}_{i}\,{\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{k}\,{\hat {a}}_{l}^{\dagger }\,{\mathclose {:}}+{\mathopen {:}}\,{\hat {a}}_{i}^{\bullet }\,{\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{k}\,{\hat {a}}_{l}^{\dagger \bullet }\,{\mathclose {:}}+{\mathopen {:}}\,{\hat {a}}_{i}\,{\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{k}^{\bullet }\,{\hat {a}}_{l}^{\dagger \bullet }\,{\mathclose {:}}+{\mathopen {:}}\,{\hat {a}}_{i}^{\bullet }\,{\hat {a}}_{j}^{\dagger \bullet }\,{\hat {a}}_{k}\,{\hat {a}}_{l}^{\dagger }\,{\mathclose {:}}+\,{\mathopen {:}}{\hat {a}}_{i}^{\bullet }\,{\hat {a}}_{j}^{\dagger \bullet }\,{\hat {a}}_{k}^{\bullet \bullet }\,{\hat {a}}_{l}^{\dagger \bullet \bullet }{\mathclose {:}}}$

In the last line we have used different numbers of ${\displaystyle ^{\bullet }}$ symbols to denote different contractions. By repeatedly applying the commutation relations it takes a lot of work to express ${\displaystyle {\hat {a}}_{i}\,{\hat {a}}_{j}^{\dagger }\,{\hat {a}}_{k}\,{\hat {a}}_{l}^{\dagger }}$ in the form of a sum of normally ordered products. It is an even lengthier calculation for more complicated products.

Luckily Wick's theorem provides a shortcut.

A product of creation and annihilation operators ${\displaystyle {\hat {A}}{\hat {B}}{\hat {C}}{\hat {D}}{\hat {E}}{\hat {F}}\ldots }$ can be expressed as

${\displaystyle {\begin{aligned}{\hat {A}}{\hat {B}}{\hat {C}}{\hat {D}}{\hat {E}}{\hat {F}}\ldots &={\mathopen {:}}{\hat {A}}{\hat {B}}{\hat {C}}{\hat {D}}{\hat {E}}{\hat {F}}\ldots {\mathclose {:}}\\&\quad +\sum _{\text{singles}}{\mathopen {:}}{\hat {A}}^{\bullet }{\hat {B}}^{\bullet }{\hat {C}}{\hat {D}}{\hat {E}}{\hat {F}}\ldots {\mathclose {:}}\\&\quad +\sum _{\text{doubles}}{\mathopen {:}}{\hat {A}}^{\bullet }{\hat {B}}^{\bullet \bullet }{\hat {C}}^{\bullet \bullet }{\hat {D}}^{\bullet }{\hat {E}}{\hat {F}}\ldots {\mathclose {:}}\\&\quad +\ldots \end{aligned}}}$

In other words, a string of creation and annihilation operators can be rewritten as the normal-ordered product of the string, plus the normal-ordered product after all single contractions among operator pairs, plus all double contractions, etc., plus all full contractions.

Applying the theorem to the above examples provides a much quicker method to arrive at the final expressions.

A warning: In terms on the right hand side containing multiple contractions care must be taken when the operators are fermionic. In this case an appropriate minus sign must be introduced according to the following rule: rearrange the operators (introducing minus signs whenever the order of two fermionic operators is swapped) to ensure the contracted terms are adjacent in the string. The contraction can then be applied (See "Rule C" in Wick's paper).

Example:

If we have two fermions (${\displaystyle N=2}$) with creation and annihilation operators ${\displaystyle {\hat {f}}_{i}^{\dagger }}$ and ${\displaystyle {\hat {f}}_{i}}$ (${\displaystyle i=1,2}$) then

${\displaystyle {\begin{aligned}{\hat {f}}_{1}\,{\hat {f}}_{2}\,{\hat {f}}_{1}^{\dagger }\,{\hat {f}}_{2}^{\dagger }\,={}&\,{\mathopen {:}}{\hat {f}}_{1}\,{\hat {f}}_{2}\,{\hat {f}}_{1}^{\dagger }\,{\hat {f}}_{2}^{\dagger }\,{\mathclose {:}}\\[5pt]&-\,{\hat {f}}_{1}^{\bullet }\,{\hat {f}}_{1}^{\dagger \bullet }\,\,{\mathopen {:}}{\hat {f}}_{2}\,{\hat {f}}_{2}^{\dagger }\,{\mathclose {:}}+\,{\hat {f}}_{1}^{\bullet }\,{\hat {f}}_{2}^{\dagger \bullet }\,\,{\mathopen {:}}{\hat {f}}_{2}\,{\hat {f}}_{1}^{\dagger }\,{\mathclose {:}}+\,{\hat {f}}_{2}^{\bullet }\,{\hat {f}}_{1}^{\dagger \bullet }\,\,{\mathopen {:}}{\hat {f}}_{1}\,{\hat {f}}_{2}^{\dagger }\,{\mathclose {:}}-{\hat {f}}_{2}^{\bullet }\,{\hat {f}}_{2}^{\dagger \bullet }\,\,{\mathopen {:}}{\hat {f}}_{1}\,{\hat {f}}_{1}^{\dagger }\,{\mathclose {:}}\\[5pt]&-{\hat {f}}_{1}^{\bullet \bullet }\,{\hat {f}}_{1}^{\dagger \bullet \bullet }\,{\hat {f}}_{2}^{\bullet }\,{\hat {f}}_{2}^{\dagger \bullet }\,+{\hat {f}}_{1}^{\bullet \bullet }\,{\hat {f}}_{2}^{\dagger \bullet \bullet }\,{\hat {f}}_{2}^{\bullet }\,{\hat {f}}_{1}^{\dagger \bullet }\,\end{aligned}}}$

Note that the term with contractions of the two creation operators and of the two annihilation operators is not included because their contractions vanish.

We use induction to prove the theorem for bosonic creation and annihilation operators. The ${\displaystyle N=2}$ base case is trivial, because there is only one possible contraction:

${\displaystyle {\hat {A}}{\hat {B}}={\mathopen {:}}{\hat {A}}{\hat {B}}{\mathclose {:}}+({\hat {A}}\,{\hat {B}}\,-{\mathopen {:}}{\hat {A}}\,{\hat {B}}{\mathclose {:}})={\mathopen {:}}{\hat {A}}{\hat {B}}{\mathclose {:}}+{\hat {A}}^{\bullet }{\hat {B}}^{\bullet }}$

In general, the only non-zero contractions are between an annihilation operator on the left and a creation operator on the right. Suppose that Wick's theorem is true for ${\displaystyle N-1}$ operators ${\displaystyle {\hat {B}}{\hat {C}}{\hat {D}}{\hat {E}}{\hat {F}}\ldots }$, and consider the effect of adding an Nth operator ${\displaystyle {\hat {A}}}$ to the left of ${\displaystyle {\hat {B}}{\hat {C}}{\hat {D}}{\hat {E}}{\hat {F}}\ldots }$ to form ${\displaystyle {\hat {A}}{\hat {B}}{\hat {C}}{\hat {D}}{\hat {E}}{\hat {F}}\ldots }$. By Wick's theorem applied to ${\displaystyle N-1}$ operators, we have:

${\displaystyle {\begin{aligned}{\hat {A}}{\hat {B}}{\hat {C}}{\hat {D}}{\hat {E}}{\hat {F}}\ldots &={\hat {A}}{\mathopen {:}}{\hat {B}}{\hat {C}}{\hat {D}}{\hat {E}}{\hat {F}}\ldots {\mathclose {:}}\\&\quad +{\hat {A}}\sum _{\text{singles}}{\mathopen {:}}{\hat {B}}^{\bullet }{\hat {C}}^{\bullet }{\hat {D}}{\hat {E}}{\hat {F}}\ldots {\mathclose {:}}\\&\quad +{\hat {A}}\sum _{\text{doubles}}{\mathopen {:}}{\hat {B}}^{\bullet }{\hat {C}}^{\bullet \bullet }{\hat {D}}^{\bullet \bullet }{\hat {E}}^{\bullet }{\hat {F}}\ldots {\mathclose {:}}\\&\quad +{\hat {A}}\ldots \end{aligned}}}$

${\displaystyle {\hat {A}}}$ is either a creation operator or an annihilation operator. If ${\displaystyle {\hat {A}}}$ is a creation operator, all above products, such as ${\displaystyle {\hat {A}}{\mathopen {:}}{\hat {B}}{\hat {C}}{\hat {D}}{\hat {E}}{\hat {F}}\ldots {\mathclose {:}}}$, are already normal ordered and require no further manipulation. Because ${\displaystyle {\hat {A}}}$ is to the left of all annihilation operators in ${\displaystyle {\hat {A}}{\hat {B}}{\hat {C}}{\hat {D}}{\hat {E}}{\hat {F}}\ldots }$, any contraction involving it will be zero. Thus, we can add all contractions involving ${\displaystyle {\hat {A}}}$ to the sums without changing their value. Therefore, if ${\displaystyle {\hat {A}}}$ is a creation operator, Wick's theorem holds for ${\displaystyle {\hat {A}}{\hat {B}}{\hat {C}}{\hat {D}}{\hat {E}}{\hat {F}}\ldots }$.

Now, suppose that ${\displaystyle {\hat {A}}}$ is an annihilation operator. To move ${\displaystyle {\hat {A}}}$ from the left-hand side to the right-hand side of all the products, we repeatedly swap ${\displaystyle {\hat {A}}}$ with the operator immediately right of it (call it ${\displaystyle {\hat {X}}}$), each time applying ${\displaystyle {\hat {A}}{\hat {X}}={\mathopen {:}}{\hat {A}}{\hat {X}}{\mathclose {:}}+{\hat {A}}^{\bullet }{\hat {X}}^{\bullet }}$ to account for noncommutativity. Once we do this, all terms will be normal ordered. All terms added to the sums by pushing ${\displaystyle {\hat {A}}}$ through the products correspond to additional contractions involving ${\displaystyle {\hat {A}}}$. Therefore, if ${\displaystyle {\hat {A}}}$ is an annihilation operator, Wick's theorem holds for ${\displaystyle {\hat {A}}{\hat {B}}{\hat {C}}{\hat {D}}{\hat {E}}{\hat {F}}\ldots }$.

We have proved the base case and the induction step, so the theorem is true. By introducing the appropriate minus signs, the proof can be extended to fermionic creation and annihilation operators. The theorem applied to fields is proved in essentially the same way.^[3]

The correlation function that appears in quantum field theory can be expressed by a contraction on the field operators:

${\displaystyle {\mathcal {C}}(x_{1},x_{2})=\left\langle 0\mid {\mathcal {T}}\phi _{i}(x_{1})\phi _{i}(x_{2})\mid 0\right\rangle =\langle 0\mid {\overline {\phi _{i}(x_{1})\phi _{i}(x_{2})}}\mid 0\rangle =i\Delta _{F}(x_{1}-x_{2})=i\int {{\frac {d^{4}k}{(2\pi )^{4}}}{\frac {e^{-ik(x_{1}-x_{2})}}{(k^{2}-m^{2})+i\epsilon }}},}$

where the operator ${\displaystyle {\overline {\phi _{i}(x_{1})\phi _{i}(x_{2})}}}$ are the amount that do not annihilate the vacuum state ${\displaystyle |0\rangle }$. Which means that ${\displaystyle {\overline {AB}}={\mathcal {T}}AB-{\mathopen {:}}{\mathcal {T}}AB{\mathclose {:}}}$. This means that ${\displaystyle {\overline {AB}}}$ is a contraction over ${\displaystyle {\mathcal {T}}AB}$. Note that the contraction of a time-ordered string of two field operators is a C-number.

In the end, we arrive at Wick's theorem:

The T-product of a time-ordered free fields string can be expressed in the following manner:

${\displaystyle {\mathcal {T}}\prod _{k=1}^{m}\phi (x_{k})={\mathopen {:}}{\mathcal {T}}\prod \phi _{i}(x_{k}){\mathclose {:}}+\sum _{\alpha ,\beta }{\overline {\phi (x_{\alpha })\phi (x_{\beta })}}{\mathopen {:}}{\mathcal {T}}\prod _{k\not =\alpha ,\beta }\phi _{i}(x_{k}){\mathclose {:}}+{}}$

${\displaystyle {\mathcal {}}+\sum _{(\alpha ,\beta ),(\gamma ,\delta )}{\overline {\phi (x_{\alpha })\phi (x_{\beta })}}\;{\overline {\phi (x_{\gamma })\phi (x_{\delta })}}{\mathopen {:}}{\mathcal {T}}\prod _{k\not =\alpha ,\beta ,\gamma ,\delta }\phi _{i}(x_{k}){\mathclose {:}}+\cdots .}$

Applying this theorem to S-matrix elements, we discover that normal-ordered terms acting on vacuum state give a null contribution to the sum. We conclude that m is even and only completely contracted terms remain.

${\displaystyle F_{m}^{i}(x)=\left\langle 0\mid {\mathcal {T}}\phi _{i}(x_{1})\phi _{i}(x_{2})\mid 0\right\rangle =\sum _{\mathrm {pairs} }{\overline {\phi (x_{1})\phi (x_{2})}}\cdots {\overline {\phi (x_{m-1})\phi (x_{m}}})}$

${\displaystyle G_{p}^{(n)}=\left\langle 0\mid {\mathcal {T}}{\mathopen {:}}v_{i}(y_{1}){\mathclose {:}}\dots {\mathopen {:}}v_{i}(y_{n}){\mathclose {:}}\phi _{i}(x_{1})\cdots \phi _{i}(x_{p})\mid 0\right\rangle }$

where p is the number of interaction fields (or, equivalently, the number of interacting particles) and n is the development order (or the number of vertices of interaction). For example, if ${\displaystyle v=gy^{4}\Rightarrow {\mathopen {:}}v_{i}(y_{1}){\mathclose {:}}={\mathopen {:}}\phi _{i}(y_{1})\phi _{i}(y_{1})\phi _{i}(y_{1})\phi _{i}(y_{1}){\mathclose {:}}}$

This is analogous to the corresponding Isserlis' theorem in statistics for the moments of a Gaussian distribution.

Note that this discussion is in terms of the usual definition of normal ordering which is appropriate for the vacuum expectation values (VEV's) of fields. (Wick's theorem provides as a way of expressing VEV's of n fields in terms of VEV's of two fields.^[4]) There are any other possible definitions of normal ordering, and Wick's theorem is valid irrespective. However Wick's theorem only simplifies computations if the definition of normal ordering used is changed to match the type of expectation value wanted. That is we always want the expectation value of the normal ordered product to be zero. For instance in thermal field theory a different type of expectation value, a thermal trace over the density matrix, requires a different definition of normal ordering.^[5]

• Isserlis' theorem

• Peskin, M. E.; Schroeder, D. V. (1995). An Introduction to Quantum Field Theory. Perseus Books. (§4.3)
• Schweber, Silvan S. (1962). An Introduction to Relativistic Quantum Field Theory. New York: Harper and Row. (Chapter 13, Sec c)