Use and . 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 , with the positive-frequency part of a quantum field containing annihilation operators and containing creation operators. The mode normalization in the PDF corresponds to and .
The time-ordered product places fields with the latest time on the left. For two bosonic insertions,
The normal-ordered product 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 function
Define the Feynman propagator here to be the two-point function itself, rather than that function divided by . Then
Equivalently,
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 . In particular, : the difference between time ordering and normal ordering is the Wick contraction.
For free bosonic fields, the Wick theorem states
where runs over all sets of disjoint pairs of indices, including the empty set, and 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 is the latest. For any remaining set , commuting the annihilation part of through the creation parts in a normal-ordered product gives
The creation part of is already in the correct position. Every commutator on the right is a scalar, so no further operator rearrangement is needed for it. Because is latest, that commutator is . Assume the Wick theorem for the other fields and multiply its expansion by . 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 . 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 . These facts explain why propagators are the edges of Feynman diagrams.
For the interacting quartic scalar field theory, the interaction Hamiltonian density is . With a single time ordering acting on external insertions and interaction vertices, the unnormalized numerator expands as
Apply the proved Wick theorem term by term. Each complete pairing connects external insertions and four-legged vertices by Feynman propagators. There are assignments of four distinct incident legs at a vertex, canceling its factorial. Permuting identical integrated vertices cancels the corresponding part of . Permutations still leaving the labeled graph unchanged give its Feynman-diagram symmetry factor. For example,
whereas the two-point tadpole diagram has coefficient , because choosing the two vertex legs contracted with the labeled external fields gives possibilities and the other two contract with each other.
Thus the position-space rules attach to every edge, 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 cancels those disconnected vacuum factors. They should not be discarded without that normalization.
There is a time-ordering convention in the printed correlation-function notation worth making explicit. The usual interacting correlator uses the single global time ordering in above. If the displayed is read literally with outside that ordering, the external block and the vertex block are separately time ordered: their cross-contractions are because each external field stands to the left of every vertex field. The same pairing construction then uses for lines within either block and for lines between blocks; in momentum space the latter is . This literal product is different from the conventional globally time-ordered interacting correlator, and does not have the ordinary external Feynman propagators derived above.