In quantum field theory, a fermionic field is a quantum field whose quanta are fermions; that is, they obey Fermi–Dirac statistics. Fermionic fields obey canonical anticommutation relations rather than the canonical commutation relations of bosonic fields.

The most prominent example of a fermionic field is the Dirac field, which describes fermions with spin-1/2: electrons, protons, quarks, etc. The Dirac field can be described as either a 4-component spinor or as a pair of 2-component Weyl spinors. Spin-1/2 Majorana fermions, such as the hypothetical neutralino, can be described as either a dependent 4-component Majorana spinor or a single 2-component Weyl spinor. It is not known whether the neutrino is a Majorana fermion or a Dirac fermion; observing neutrinoless double-beta decay experimentally would settle this question.

Free (non-interacting) fermionic fields obey canonical anticommutation relations; i.e., involve the anticommutators {a, b} = ab + ba, rather than the commutators [a, b] = ab − ba of bosonic or standard quantum mechanics. Those relations also hold for interacting fermionic fields in the interaction picture, where the fields evolve in time as if free and the effects of the interaction are encoded in the evolution of the states.

It is these anticommutation relations that imply Fermi–Dirac statistics for the field quanta. They also result in the Pauli exclusion principle: two fermionic particles cannot occupy the same state at the same time.

The prominent example of a spin-1/2 fermion field is the Dirac field (named after Paul Dirac), and denoted by ${\displaystyle \psi (x)}$. The equation of motion for a free spin 1/2 particle is the Dirac equation,

${\displaystyle \left(i\gamma ^{\mu }\partial _{\mu }-m\right)\psi (x)=0.\,}$

where ${\displaystyle \gamma ^{\mu }}$ are gamma matrices and ${\displaystyle m}$ is the mass. The simplest possible solutions ${\displaystyle \psi (x)}$ to this equation are plane wave solutions, ${\displaystyle u(p)e^{-ip\cdot x}\,}$ and ${\displaystyle v(p)e^{ip\cdot x}\,}$. These plane wave solutions form a basis for the Fourier components of ${\displaystyle \psi (x)}$, allowing for the general expansion of the wave function as follows,

${\displaystyle \psi _{\alpha }(x)=\int {\frac {d^{3}p}{(2\pi )^{3}}}{\frac {1}{\sqrt {2E_{p}}}}\sum _{s}\left(a_{\mathbf {p} }^{s}u_{\alpha }^{s}(p)e^{-ip\cdot x}+b_{\mathbf {p} }^{s\dagger }v_{\alpha }^{s}(p)e^{ip\cdot x}\right).\,}$

u and v are spinors, labelled by spin, s and spinor indices ${\displaystyle \alpha \in \{0,1,2,3\}}$. For the electron, a spin 1/2 particle, s = +1/2 or s = −1/2. The energy factor is the result of having a Lorentz invariant integration measure. In second quantization, ${\displaystyle \psi (x)}$ is promoted to an operator, so the coefficients of its Fourier modes must be operators too. Hence, ${\displaystyle a_{\mathbf {p} }^{s}}$ and ${\displaystyle b_{\mathbf {p} }^{s\dagger }}$ are operators. The properties of these operators can be discerned from the properties of the field. ${\displaystyle \psi (x)}$ and ${\displaystyle \psi (y)^{\dagger }}$ obey the anticommutation relations:

${\displaystyle \left\{\psi _{\alpha }(\mathbf {x} ),\psi _{\beta }^{\dagger }(\mathbf {y} )\right\}=\delta ^{(3)}(\mathbf {x} -\mathbf {y} )\delta _{\alpha \beta }.}$

We impose an anticommutator relation (as opposed to a commutation relation as we do for the bosonic field) in order to make the operators compatible with Fermi–Dirac statistics. By putting in the expansions for ${\displaystyle \psi (x)}$ and ${\displaystyle \psi (y)}$, the anticommutation relations for the coefficients can be computed.

${\displaystyle \left\{a_{\mathbf {p} }^{r},a_{\mathbf {q} }^{s\dagger }\right\}=\left\{b_{\mathbf {p} }^{r},b_{\mathbf {q} }^{s\dagger }\right\}=(2\pi )^{3}\delta ^{3}(\mathbf {p} -\mathbf {q} )\delta ^{rs},\,}$

In a manner analogous to non-relativistic annihilation and creation operators and their commutators, these algebras lead to the physical interpretation that ${\displaystyle a_{\mathbf {p} }^{s\dagger }}$ creates a fermion of momentum p and spin s, and ${\displaystyle b_{\mathbf {q} }^{r\dagger }}$ creates an antifermion of momentum q and spin r. The general field ${\displaystyle \psi (x)}$ is now seen to be a weighted (by the energy factor) summation over all possible spins and momenta for creating fermions and antifermions. Its conjugate field, ${\displaystyle {\overline {\psi }}\ {\stackrel {\mathrm {def} }{=}}\ \psi ^{\dagger }\gamma ^{0}}$, is the opposite, a weighted summation over all possible spins and momenta for annihilating fermions and antifermions.

With the field modes understood and the conjugate field defined, it is possible to construct Lorentz invariant quantities for fermionic fields. The simplest is the quantity ${\displaystyle {\overline {\psi }}\psi \,}$. This makes the reason for the choice of ${\displaystyle {\overline {\psi }}=\psi ^{\dagger }\gamma ^{0}}$clear. This is because the general Lorentz transform on ${\displaystyle \psi }$ is not unitary so the quantity ${\displaystyle \psi ^{\dagger }\psi }$ would not be invariant under such transforms, so the inclusion of ${\displaystyle \gamma ^{0}\,}$ is to correct for this. The other possible non-zero Lorentz invariant quantity, up to an overall conjugation, constructible from the fermionic fields is ${\displaystyle {\overline {\psi }}\gamma ^{\mu }\partial _{\mu }\psi }$.

Since linear combinations of these quantities are also Lorentz invariant, this leads naturally to the Lagrangian density for the Dirac field by the requirement that the Euler–Lagrange equation of the system recover the Dirac equation.

${\displaystyle {\mathcal {L}}_{D}={\overline {\psi }}\left(i\gamma ^{\mu }\partial _{\mu }-m\right)\psi \,}$

Such an expression has its indices suppressed. When reintroduced the full expression is

${\displaystyle {\mathcal {L}}_{D}={\overline {\psi }}_{a}\left(i\gamma _{ab}^{\mu }\partial _{\mu }-m\mathbb {I} _{ab}\right)\psi _{b}\,}$

The Hamiltonian (energy) density can also be constructed by first defining the momentum canonically conjugate to ${\displaystyle \psi (x)}$, called ${\displaystyle \Pi (x):}$

${\displaystyle \Pi \ {\overset {\mathrm {def} }{=}}\ {\frac {\partial {\mathcal {L}}_{D}}{\partial (\partial _{0}\psi )}}=i\psi ^{\dagger }\,.}$

With that definition of ${\displaystyle \Pi }$, the Hamiltonian density is:

${\displaystyle {\mathcal {H}}_{D}={\overline {\psi }}\left[-i{\vec {\gamma }}\cdot {\vec {\nabla }}+m\right]\psi \,,}$

where ${\displaystyle {\vec {\nabla }}}$ is the standard gradient of the space-like coordinates, and ${\displaystyle {\vec {\gamma }}}$ is a vector of the space-like ${\displaystyle \gamma }$ matrices. It is surprising that the Hamiltonian density doesn't depend on the time derivative of ${\displaystyle \psi }$, directly, but the expression is correct.

Given the expression for ${\displaystyle \psi (x)}$ we can construct the Feynman propagator for the fermion field:

${\displaystyle D_{F}(x-y)=\left\langle 0\left|T(\psi (x){\overline {\psi }}(y))\right|0\right\rangle }$

we define the time-ordered product for fermions with a minus sign due to their anticommuting nature

${\displaystyle T\left[\psi (x){\overline {\psi }}(y)\right]\ {\overset {\text{def}}{=}}\ \theta \left(x^{0}-y^{0}\right)\psi (x){\overline {\psi }}(y)-\theta \left(y^{0}-x^{0}\right){\overline {\psi }}(y)\psi (x).}$

Plugging our plane wave expansion for the fermion field into the above equation yields:

${\displaystyle D_{F}(x-y)=\int {\frac {d^{4}p}{(2\pi )^{4}}}{\frac {i({p\!\!\!/}+m)}{p^{2}-m^{2}+i\epsilon }}e^{-ip\cdot (x-y)}}$

where we have employed the Feynman slash notation. This result makes sense since the factor

${\displaystyle {\frac {i({p\!\!\!/}+m)}{p^{2}-m^{2}}}}$

is just the inverse of the operator acting on ${\displaystyle \psi (x)}$ in the Dirac equation. Note that the Feynman propagator for the Klein–Gordon field has this same property. Since all reasonable observables (such as energy, charge, particle number, etc.) are built out of an even number of fermion fields, the commutation relation vanishes between any two observables at spacetime points outside the light cone. As we know from elementary quantum mechanics two simultaneously commuting observables can be measured simultaneously. We have therefore correctly implemented Lorentz invariance for the Dirac field, and preserved causality.

More complicated field theories involving interactions (such as Yukawa theory, or quantum electrodynamics) can be analyzed too, by various perturbative and non-perturbative methods.

Dirac fields are an important ingredient of the Standard Model.

• Dirac equation
• Spin–statistics theorem
• Spinor
• Composite Field
• Auxiliary Field

• Edwards, D. (1981). "The Mathematical Foundations of Quantum Field Theory: Fermions, Gauge Fields, and Super-symmetry, Part I: Lattice Field Theories". Int. J. Theor. Phys. 20 (7): 503–517. Bibcode:1981IJTP...20..503E. doi:10.1007/BF00669437. S2CID 120108219.
• Peskin, M and Schroeder, D. (1995). An Introduction to Quantum Field Theory, Westview Press. (See pages 35–63.)
• Srednicki, Mark (2007). Quantum Field Theory Archived 2011-07-25 at the Wayback Machine, Cambridge University Press, ISBN 978-0-521-86449-7.
• Weinberg, Steven (1995). The Quantum Theory of Fields, (3 volumes) Cambridge University Press.