For with fixed noncommuting Hermitian operators, . A constant ratio of the envelopes is sufficient to remove time ordering. Independently changing that ratio generally produces commutator corrections in the Magnus expansion, and a product of successive pulses need not equal the exponential of their summed areas.
Impulsive quantum kick 2026-10-14
An impulsive quantum kick is the limit of a pulse with finite integrated Hamiltonian operator, rather than a bounded finite step. For a bounded self-adjoint operator with units of action, a rectangular pulse lasting has exact unitary time evolution for every . In the limiting notation , is the Dirac delta function and
The integral estimate for quantum state continuity at a finite Hamiltonian step no longer tends to zero: the pulse height diverges. A bounded background Hamiltonian operator contributes vanishing evolution during this limit. A more general pulse with noncommuting generators can require time ordering; integrated area alone then need not specify the result.
For example, , with a Pauli matrix, sends a qubit in to . At , the Born rule predicts with probability one, however short the regularized pulse is. Its finite rotation is precisely what a bounded finite step cannot accomplish at the switching instant.
For a Hermitian scalar field in a positive-norm Hilbert space with a Lorentz-invariant vacuum, the spectrum condition and completeness give the connected vacuum two-point function as a superposition of free mass-shell functions. Define the Källén–Lehmann spectral density by , where and its vacuum expectation value has been subtracted if nonzero. Squared intermediate-state matrix elements make a nonnegative measure supported on nonnegative mass squared. Applying time ordering to each free mass shell gives the displayed Feynman propagator representation, using and signature . Its positivity does not automatically extend to gauge-variant fields in an indefinite auxiliary space.
For the complex scalar field, the canonical momenta are and . Split . The interaction picture is obtained from the Schrödinger picture by and , choosing a common time origin. The fields therefore obey the free Klein-Gordon equation, while the states evolve through . The canonical quantization of a complex scalar field gives
The adjoint expansion contains and . Particle and antiparticle modes are independent; their nonzero commutators are
All cross-species commutators, and all pairs of creation operators or of annihilation operators, vanish. These formulas reproduce and the remaining vanishing equal-time commutators.
There are no derivative interactions, so the interaction Hamiltonian is . Differentiating the interaction picture state gives . Its time-evolution operator satisfies
Iteration gives ordered integration regions . Each is of the full region with a time-ordered product, yielding the Dyson series
With the usual asymptotic scattering prescription, taking and proves
Time ordering is essential because different interaction Hamiltonians need not commute.
At first order, use the connected tree scattering amplitude, equivalently the normal-ordered quartic vertex with the vacuum and one-particle mass contributions removed. This qualification matters if the local product is interpreted literally without normal ordering: single-vertex tadpole diagrams and vacuum diagrams are not the requested connected two-particle amplitude. The complex scalar quartic contact vertex is obtained from
In a fully connected matrix element, an incoming particle and outgoing antiparticle attach to the two factors; an outgoing particle and incoming antiparticle attach to the two factors. Each pair has two assignments, giving . The external relativistic normalization of a one-particle state cancels each mode measure: for example and . Consequently the connected integrand is , and
Thus the invariant connected amplitude is
The Dirac delta function expresses four-momentum conservation.
For the identity contribution, commute through and independently through ; all terms leaving a creation operator against the left Fock vacuum or an annihilation operator against the right Fock vacuum vanish. Hence
There is no exchange term between the two distinct particle and antiparticle species.
Use units and the Minkowski metric , so . The canonical momentum of the real scalar field is , and the Hamiltonian operator is
Canonical quantization of a real scalar field promotes and to operator-valued distributions satisfying the equal-time canonical commutation relation
Products and operator commutators can first be interpreted after spatial smearing. An additive vacuum constant in does not affect the following equations. In the Heisenberg picture, . Applying the canonical commutation relation to the quadratic Hamiltonian operator, and integrating the gradient term by parts, gives
Thus the Klein-Gordon equation is an operator equation. The spatial Fourier transform reduces it to , with . Its two frequency branches, together with the Hermitian adjoint condition , give the mode expansion of a free field
The measure is Lorentz invariant; its normalization fixes that of the creation and annihilation operators rather than imposing any extra physical assumption.
For invariant-normalized scalar mode extraction, take the Klein-Gordon inner product of with the field. The positive-frequency contribution has coefficient and spatial integral . The negative-frequency contribution has coefficient and spatial integral , so it vanishes because . Consequently
Taking the Hermitian adjoint, with the order of the scalar factors immaterial, yields
Equivalently and . At a common time their operator commutator therefore is
The same calculation gives : the surviving factor is an energy difference supported at opposite spatial momenta. These are the bosonic canonical commutation relations in relativistic normalization of a one-particle state.
Choose the Fock vacuum by and . Only the annihilator-creator product contributes to the Wightman function, giving
Using time ordering and the reverse product for , the Feynman propagator, with the convention in this calculation, is
In particular the denominator in the spatial measure is in both terms.
To verify the four-dimensional Fourier transform representation, let and integrate the energy variable first:
The Feynman i-epsilon prescription puts the positive-energy pole at and the negative-energy pole at . For , the contour integral closes in the lower half-plane with clockwise orientation. The residue theorem multiplies the residue by , giving . For , it closes counterclockwise above: multiplying the residue by gives . Thus . After inserting the spatial exponential, reversing in the negative-time term gives exactly the time ordering expression above. Hence
This equality is understood as a distribution limit; the same prescription gives .
All three fields have mass . In the interaction picture the quantum fields evolve with the free Hamiltonian operator, so their equations are
The interactions evolve the states instead. For the appropriate mode expansion of a free field is
Independent species have commuting creation and annihilation operators. The full set of canonical commutation relations is
These follow from the separate equal-time canonical commutation relations for the three real scalar fields. For comparison, in the Heisenberg picture the full Euler-Lagrange field equations contain interactions: , , and . These are distinct from the free equations used in the interaction picture.
There are no derivative couplings, so the interaction Hamiltonian is . The interaction picture evolution operator obeys
Iterating its integral equation gives the Dyson series
The ordered integration region is equivalent to the unrestricted region with time ordering and the factor . Taking the scattering limit , , with the usual adiabatic prescription for perturbation theory, therefore yields the S-matrix
All fields inside this expression are interaction picture fields.
Use relativistic normalization of a one-particle state and take , . Define the connected scattering amplitude by . The identity contribution describes no scattering; the term linear in cannot connect these four external particles. At quadratic order, the connected term in the Dyson series is
By Wick theorem, a connected Wick contraction must attach the two external legs to at one vertex, the two external legs to at the other, and contract the remaining two fields. The two assignments of species to cancel the Dyson factor . At either vertex, the assignments of its identical fields to the incoming and outgoing particle cancel the interaction factor . Equivalently, the coefficient is multiplied by , giving . Therefore, writing ,
Insert the Feynman propagator . The two vertex integrals impose and . With the Mandelstam invariant , this gives
This is distinguishable scalar scattering by a shared mediator with mediator mass . In physical elastic kinematics , so the denominator has no pole for . There is only the displayed -channel tree-level Feynman diagram: an - or -channel would require an absent mixed vertex. Each vertex contributes and the exchanged contributes , agreeing with the direct Wick contraction calculation.
Figure 1.
Tree-level elastic scattering of distinct scalar particles through a shared scalar mediator
. Details.
Use a real scalar field and Minkowski action , with the vacuum selected by a Feynman propagator pole prescription or adiabatic vacuum boundary conditions. The normalized correlation functions are
Here is the interacting vacuum and denotes time ordering. These are the canonical vacuum expectation values, rather than unordered products of noncommuting fields. The functional integral formula assumes a regulated measure and the same vacuum prescription on both sides.
Couple an external source linearly to the scalar field and choose the normalized vacuum generating functional
Then functional derivatives insert the fields:
In the convention used below, the connected generating functional is
The logarithm removes products of disconnected components. For example the connected correlation function ; at three points one subtracts all three products and adds . Normalization has already removed vacuum bubbles.
Split the free interaction picture field as , where the first part contains annihilation operators and the second creation operators. Time ordering places the operator at the later time on the left, while normal ordering places all creation operators on the left. For , commuting the annihilation part of the first field through the creation part of the second gives
The commutator is the c-number
For , the same calculation interchanges . Combining the two time orders proves the two-field Wick contraction identity
where
The equality with the four-dimensional integral follows by closing the energy contour at the positive-energy pole for and the negative-energy pole for the reverse order. This fixes the scalar Feynman propagator pole prescription; in this convention itself includes the factor .
Because these interactions have no derivatives, . Expand the Dyson series
Apply the Wick theorem to each term together with the external insertions. Every paired field contributes a propagator, and the remaining normal-ordered vacuum expectation vanishes unless no fields remain. Representing the propagators by Fourier integrals and integrating each vertex position gives the four-momentum conservation delta function. The and contractions into labelled vertex slots cancel the factorials in the interaction density. The resulting scalar-field cubic and quartic Feynman vertices are
Each vertex also has , each independent internal momentum is integrated with , and repeated equivalent contractions give the usual Feynman-diagram symmetry factor. For the unamputated correlation functions, retain the propagators joining external insertions to vertices. The displayed unnormalized numerator includes disconnected vacuum bubbles; dividing by cancels them. Amputated scattering amplitudes additionally use the LSZ reduction formula.
For stability, examine the classical scalar potential
The tree-level boundedness criterion is , with arbitrary real cubic coupling. If , boundedness requires and ; if , or with , the potential is unbounded below. Thus a positive quartic coupling permits a stable ground-state vacuum, although it need not be the perturbative vacuum at .
For and , the stability of a cubic-quartic scalar potential determines when zero is a global minimum:
The square bracket is nonnegative for every field value exactly when
For a nontrivial equality case there is a second degenerate minimum at . If and , zero remains a local minimum but is a false vacuum, with a lower minimum elsewhere. If and , the origin is the quartic minimum despite its vanishing quadratic curvature. These are classical or tree-level stability statements; renormalized couplings and the quantum effective potential refine the quantum analysis.
Figure 1.
Stable, degenerate and metastable scalar vacua
. Details.
Use signature and covariant one-particle state normalization . The future scattering condition uses the outgoing asymptotic scalar field:
As with the incoming condition, this is understood in suitably smeared scattering matrix elements, not as unrestricted pointwise strong convergence of interacting and free local operators.
The free incoming field has expansion
The creation operator defines , with , so . On the other hand, translation invariance gives
Matching the stable one-particle state component in the asymptotic condition fixes the constant overlap. Choose its state phase so that it is real and nonnegative. Then
The same stable particle appears in the outgoing description, with compatible one-particle phases. Thus is the squared overlap measured by wave-function renormalization.
Insert a complete set of four-momentum eigenstates between the fields. For , translation invariance and Hermiticity give
For , reversing the time-ordered product gives the same sum with . Therefore, with continuous-state integrals implicit in the completeness sum,
Here is the Heaviside step function.
For the usual spectral representation, take a Lorentz-invariant vacuum, a positive-norm Hilbert space and the future-cone spectrum condition. Also take , as for the unbroken real scalar field, or subtract this constant and use the connected correlation function. Otherwise the vacuum intermediate state adds a disconnected constant not represented by the ordinary mass-shell measure below.
Let for that centered field, and use the Fourier transform convention . Define the Källén–Lehmann spectral density by the distributional equality
Lorentz invariance makes the weight on each future mass shell a function or measure of ; the spectrum condition gives . Its nonnegativity follows from the displayed squared matrix elements. Thus
Applying time ordering to this expression mass shell by mass shell proves the Källén–Lehmann spectral representation:
The factor is outside under the convention in the question.
The stable one-particle state part of completeness is . Its overlap is , so its contribution to is times the free mass- two-point function. Hence
Other isolated particles as well as multiparticle states can belong to . The full scalar propagator accordingly has isolated pole , making its pole residue.
Finally use canonical field normalization: the canonical momentum is and the equal-time canonical commutation relation is . The vacuum commutator has the same spectral decomposition. For each free mass shell,
independently of . Comparing with the interacting equal-time commutator gives the canonical scalar spectral sum rule:
For the particle genuinely interpolated by the stated asymptotic field, ; the general nonnegative bound allows a field with zero overlap. The total weight and the upper bound require this canonical normalization and a well-defined distributional equal-time limit. They do not follow from the asymptotic condition alone: rescaling rescales both and the entire spectral measure.
An antilinear involution need not preserve norms, so the required physical assumption is that momentum reversal maps a complete normalized basis to a normalized basis. For a spinless sector take with and relativistic normalization . Any extra degeneracy labels must likewise be mapped by a norm-preserving permutation or unitary matrix. The condition requires in the scalar basis.
For arbitrary normalizable wave packets and , where , antilinearity conjugates their coefficients. Using the normalized reversed basis gives
Completeness proves for all states, so this quantum time-reversal operator is antiunitary. The same proof with sums over extra labels uses the unitarity of their reversal matrix.
For a Hermitian interaction-picture Hamiltonian, the condition of time-reversal symmetry is
Expand the scattering matrix as a Dyson series, using ordered integration rather than commuting operators:
Conjugation by conjugates to but preserves the product order. Substitute ; the integration domain becomes and
On the other hand, taking the adjoint of the original series conjugates its coefficient and reverses the operator order. Relabel in each adjoint integral; it is exactly the same series. Thus
Equivalently, time ordering turns into anti-time ordering after reversing the time arguments. One can first use a finite symmetric time interval and the same time-symmetric switching prescription before taking the scattering limit; the argument never assumes .
Antiunitarity relates amplitudes by
Using and taking a complex conjugate gives the time-reversal reciprocity of the scattering matrix
The reverse process starts with the reversed final state and ends with the reversed initial state. Intrinsic state phases can appear when these reversed states are rewritten in a separately chosen basis, but the corresponding squared amplitudes agree. Equal amplitudes do not mean identical numerical cross sections when the reversed process has different flux or phase-space factors.
For a spin-zero complex scalar field, choose the transformation convention . Its mode expansion of a free field is
Apply to the expansion. Its antiunitarity conjugates both numerical plane waves, giving coefficients multiplying and multiplying . In , replace by to obtain precisely these same plane waves. Comparing independent modes proves time reversal of a complex scalar field:
Taking adjoints gives . Positive energy is preserved, and an annihilation operator remains an annihilation operator. Time reversal does not exchange particles with antiparticles.
The coefficient is an intrinsic time-reversal phase. Applying twice gives , so ; it need not be restricted to . Write and redefine . Because the phase is conjugated by ,
Thus the scalar intrinsic phase can be absorbed into the annihilation-operator convention. Redefining simultaneously implements the corresponding whole-field rephasing and removes the phase from both mode transformations.
Write and . By time ordering for a bosonic field,
Reverse in the second integral. At fixed spatial momentum consider
The Feynman i-epsilon prescription places the poles at and . For , close in the lower half-plane, clockwise. The positive-energy residue and the orientation give . For , close counterclockwise in the upper half-plane; the negative pole has derivative denominator , and its residue gives . Thus , reproducing the time-ordered mode expression. We obtain
This is a distributional boundary value, not an absolutely convergent Fourier integral. As an independent normalization check, , so ; the numerator is with the question's definition of the propagator.
For one Grassmann variable, define the Berezin integral by and , extending linearly. Since , these rules determine it on every polynomial. With several Grassmann variables, integrate from the rightmost differential first. Equivalently, extracts the coefficient of in that order. Repeated indices give zero, while permuting distinct generators gives the fermionic sign of that permutation. Hence
where is the Levi-Civita symbol.
For the doubled variables, fix the paired measure orientation by
The exponential in the Grassmann Gaussian integral has a finite Taylor expansion. Only its term of degree survives integration. Each barred and unbarred index must occur exactly once; the orders of the even bilinears cancel the Taylor factorial. Reordering the unbarred indices gives the determinant permutation sign. Therefore
This identity holds even for a singular matrix. The normalized expectations below require invertible.
Put . Differentiating with respect to inserts . The derivative of a determinant then gives the ordered Grassmann Gaussian contractions
Differentiate once more, using . The result is
The minus sign is the interchange sign in the fermionic Wick theorem. In particular . Reordering each barred-unbarred pair in the product and contracting the indices gives
The covariance of Grassmann bilinears is the connected second term. Its Feynman diagram is one oriented fermion loop through both insertions; the first term has two separate fermion loops, one at each insertion. Each closed fermion loop contributes a minus sign, so two loops give a plus and one gives a minus.
Figure 1.
Connected and disconnected contractions of Grassmann bilinears
. Details.
For the Dirac field, write and use signature , consistent with the denominator in the question. A finite regulator turns the free fermionic path integral into the same Gaussian problem with , because . The time-ordered product covariance is consequently , where the subscript specifies the Feynman i-epsilon prescription.
With the Fourier transform convention , the derivative kernel is . The gamma matrices obey
since their anticommutator makes . Thus the mostly-plus Dirac propagator with a real derivative kernel is
For , the denominator is . Its positive-energy pole is below and negative-energy pole above the real axis, giving time ordering. Here is defined without the leading factor , as in the question; conventions that include that factor in the name Dirac propagator give the same covariance.
Split the free scalar field into its annihilation part and creation part . 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 . The oscillator commutator therefore gives the Wightman function
Put and . Combining the two time orders, and reversing in the second when necessary, gives
To 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 , close below clockwise: the residue at is , and the clockwise factor after dividing by gives the required positive coefficient. For , close above counterclockwise: the residue at is , and the factor gives the other time order. Equivalently the energy contour passes above the positive-energy pole and below the negative-energy pole. Hence
The 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 past the creation part of gives
The normal-ordered two-field product is symmetric in , 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
For each fixed phase , the single-transition Quantum Hamiltonian is , so its generators at different times commute. Its exact propagator depends only on . On the active two-level subspace, ; on the spectator state, . Thus it is an embedded two-level quantum rotation, not an involution on the whole three-dimensional space. Explicitly,
The spectator entries in both matrices are exactly one; the converted TeX damages the second matrix, so these entries use the PDF.
The proposed combined exponential is false in general for two independent envelopes. Direct multiplication gives
so
The constant pulse-axis condition for removing time ordering holds if both envelopes are proportional to the same function, but not for generic independent controls. For example, turn on only the first transition, then only the second. The result is , whose Baker--Campbell--Hausdorff formula includes a nonzero commutator correction, rather than simply .
There is also a literal source defect: the PDF's right-hand exponent omits the factor on its second term. With it would give , which has eigenvalues and is not unitary for real nonzero , whereas the true propagator is unitary. Even after repairing that missing factor, time ordering is still necessary unless the commutators vanish.
Use natural units and the Minkowski metric . The canonical momentum is , and the Legendre transform in mechanics gives
For canonical quantization of a real scalar field, impose the equal-time canonical commutation relations
The Heisenberg picture equation then gives . For the other equation, the momentum term commutes with ; commuting the gradient term past and integrating the derivative of the Dirac delta function by parts gives . Consequently
The Fourier transform reduces this Klein-Gordon equation to oscillators of frequency . Hermiticity of the real scalar field pairs their positive- and negative-frequency coefficients, so its general mode expansion can be normalized as
To verify its normalization rather than guess the oscillator algebra, use scalar field oscillator inversion at :
Its adjoint has the conjugate phase and the opposite sign in front of . The two mixed field-momentum commutators give
The Fourier representation of the Dirac delta function sets the energies equal in the last step. Similarly, has a factor multiplying and vanishes; the two-creator commutator vanishes as well.
Substituting the expansion in and integrating the spatial phases, the two-annihilator and two-creator coefficients cancel by . The free real scalar Hamiltonian in oscillator variables is
It differs from by its constant vacuum energy; normal ordering removes that constant. Using the oscillator commutators, for either Hamiltonian,
Choose the Fock vacuum annihilated by every . Each creation operator adds a particle of energy and momentum ; 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 and , vacuum contraction and time ordering give the Feynman propagator
In the scalar Feynman propagator pole prescription, the integral of has poles at and . Close clockwise below for and counterclockwise above for . The residue theorem gives exactly the corresponding term above, proving
As a direct check of the distributional source, is continuous, but . The derivative jump of a free scalar time-ordered two-point function therefore gives . Spatial Fourier inversion yields
The interaction contains no time derivatives, so . The Legendre transform in mechanics gives
Thus
Here is the free Hamiltonian operator of mass , and is the interaction Hamiltonian.
In the interaction picture, operators evolve with and the states with . With , the state evolution satisfies . Integrating once and iterating gives the Dyson series
Equivalently it is the time ordering of the exponential. Taking the infinite-time scattering limit, with adiabatic switching understood, gives the scattering matrix
The field in the interaction is the free interaction-picture field, which is essential for applying the free-field contractions.
Use relativistic normalization of a one-particle state, , so each external scalar contraction contributes its plane-wave phase without an extra factor. For the connected on-shell two-to-two scattering amplitude, the first Dyson series term is . The Wick theorem gives ways to attach the two incoming and two outgoing external particles to the four fields, canceling the denominator. Consequently
where
The delta function expresses four-momentum conservation. The connected scattering convention removes spectator, vacuum and external self-energy contributions; equivalently at this order one can normal-order the interaction and keep the physical external mass fixed. These contributions should not be mistaken for additional connected four-point diagrams.
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.
The sources are independent odd Grassmann variables and anticommute with the Dirac field. Use and the same mostly-plus gamma matrices and Dirac adjoint as in the interaction calculation. The inverse is selected with vacuum Feynman i-epsilon prescription boundary conditions. With integral kernels and spinor contractions understood, the exponent can be completed to a square:
Translations preserve the Berezin integral, so the Gaussian generating functional for a Dirac field is
The positions of the sources matter. To extract the Dirac propagator, use a left Grassmann derivative with respect to followed by a right Grassmann derivative with respect to :
The leading minus sign removes the two insertion factors . It can also be checked by differentiating the quadratic source exponential: its ordered second derivative is .
For the Fourier transform , and the Clifford algebra gives
Therefore
This equals . The fermionic time ordering is explicitly
with the minus sign supplied by exchanging odd fields. The poles put positive energy forward in time and negative energy backward, which is the antiparticle contribution. Finally, checks the numerator and overall sign. All formulas are distributional limits with the indicated boundary prescription.
Use units , the Minkowski metric , and the Fourier transform convention . Here gamma matrices satisfy , and 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, give
Before gauge fixing, the Maxwell Lagrangian has zero modes in field theory along , so its quadratic operator cannot be inverted on all potentials. In the linear Lorenz gauge, the Faddeev-Popov determinant is . It is independent of 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 is
The 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 is the photon quadratic operator and , completing the square gives
In the last expression and 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 , 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, , the momentum-space two-point functions are
The 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 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 ; integrating out the quadratic canonical momentum produces the bosonic action. The fermionic coherent state overlaps produce the first-order term 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 transformation
The gauge covariant derivative obeys . Replacing by in the Dirac action therefore gives the invariant matter density
The electromagnetic field tensor is unchanged by the local transformation. Thus the unfixed quantum electrodynamics action is gauge-invariant. The Dirac current 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 , expand in powers of . A term of order contains spacetime integrations and . Applying the Wick theorem pairs the free fields: an - Wick contraction supplies a photon propagator, and a - 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 , the same cancellation appears as
Vacuum 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 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 .
The resulting momentum-space QED Feynman rules for the bare theory can be stated in Feynman gauge as follows:
Figure 1.
An oriented electron line meets a photon at the QED interaction vertex
. Details.
The illustrated interaction vertex has incoming fermion momentum , incoming photon momentum , and outgoing fermion momentum . Its solid-line arrows indicate fermion flow. The photon wavy line has no fermion-flow arrow. The propagators, the vertex , momentum conservation, and the fermionic signs determine the perturbative amplitudes.
For ordinary bosonic operators at distinct times, time ordering places the later operator to the left:
Equivalently, multiply these two products by the corresponding Heaviside step functions. For two odd fermionic operators the second product receives a minus sign; more generally a graded interchange has sign . This sign is the fermionic sign convention for time-ordered field products, not a minus sign for arbitrary bosonic operators. No equal-time convention is needed for the two distinct times in this clause.