Relativistic wave equations

Part of a series of articles about
Quantum mechanics
${\displaystyle i\hbar {\frac {d}{dt}}|\Psi \rangle ={\hat {H}}|\Psi \rangle }$ Schrödinger equation
Introduction Glossary History
Background Classical mechanics Old quantum theory Bra–ket notation Hamiltonian Interference
Fundamentals Complementarity Decoherence Entanglement Energy level Measurement Nonlocality Quantum number State Superposition Symmetry Tunnelling Uncertainty Wave function Collapse
Experiments Bell's inequality CHSH inequality Davisson–Germer Double-slit Elitzur–Vaidman Franck–Hertz Leggett inequality Leggett–Garg inequality Mach–Zehnder Popper Quantum eraser Delayed-choice Schrödinger's cat Stern–Gerlach Wheeler's delayed-choice
Formulations Overview Heisenberg Interaction Matrix Phase-space Schrödinger Sum-over-histories (path integral)
Equations Dirac Klein–Gordon Pauli Rydberg Schrödinger
Interpretations Bayesian Consistent histories Copenhagen de Broglie–Bohm Ensemble Hidden-variable Local Superdeterminism Many-worlds Objective-collapse Quantum logic Relational Transactional Von Neumann–Wigner
Advanced topics Relativistic quantum mechanics Quantum field theory Quantum information science Quantum computing Quantum chaos EPR paradox Density matrix Scattering theory Quantum statistical mechanics Quantum machine learning
Scientists Aharonov Bell Bethe Blackett Bloch Bohm Bohr Born Bose de Broglie Compton Dirac Davisson Debye Ehrenfest Einstein Everett Fock Fermi Feynman Glauber Gutzwiller Heisenberg Hilbert Jordan Kramers Lamb Landau Laue Moseley Millikan Onnes Pauli Planck Rabi Raman Rydberg Schrödinger Simmons Sommerfeld von Neumann Weyl Wien Wigner Zeeman Zeilinger
v t e

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 Wick's theorem 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
v t e

In physics, specifically relativistic quantum mechanics (RQM) and its applications to particle physics, relativistic wave equations predict the behavior of particles at high energies and velocities comparable to the speed of light. In the context of quantum field theory (QFT), the equations determine the dynamics of quantum fields. The solutions to the equations, universally denoted as ψ or Ψ (Greek psi), are referred to as "wave functions" in the context of RQM, and "fields" in the context of QFT. The equations themselves are called "wave equations" or "field equations", because they have the mathematical form of a wave equation or are generated from a Lagrangian density and the field-theoretic Euler–Lagrange equations (see classical field theory for background).

In the Schrödinger picture, the wave function or field is the solution to the Schrödinger equation; $${\displaystyle i\hbar {\frac {\partial }{\partial t}}\psi ={\hat {H}}\psi }$$ one of the postulates of quantum mechanics. All relativistic wave equations can be constructed by specifying various forms of the Hamiltonian operator Ĥ describing the quantum system. Alternatively, Feynman's path integral formulation uses a Lagrangian rather than a Hamiltonian operator.

More generally – the modern formalism behind relativistic wave equations is Lorentz group theory, wherein the spin of the particle has a correspondence with the representations of the Lorentz group.^[1]

The failure of classical mechanics applied to molecular, atomic, and nuclear systems and smaller induced the need for a new mechanics: quantum mechanics. The mathematical formulation was led by De Broglie, Bohr, Schrödinger, Pauli, and Heisenberg, and others, around the mid-1920s, and at that time was analogous to that of classical mechanics. The Schrödinger equation and the Heisenberg picture resemble the classical equations of motion in the limit of large quantum numbers and as the reduced Planck constant ħ, the quantum of action, tends to zero. This is the correspondence principle. At this point, special relativity was not fully combined with quantum mechanics, so the Schrödinger and Heisenberg formulations, as originally proposed, could not be used in situations where the particles travel near the speed of light, or when the number of each type of particle changes (this happens in real particle interactions; the numerous forms of particle decays, annihilation, matter creation, pair production, and so on).

A description of quantum mechanical systems which could account for relativistic effects was sought for by many theoretical physicists from the late 1920s to the mid-1940s.^[2] The first basis for relativistic quantum mechanics, i.e. special relativity applied with quantum mechanics together, was found by all those who discovered what is frequently called the Klein–Gordon equation:

$${\displaystyle -\hbar ^{2}{\frac {\partial ^{2}\psi }{\partial t^{2}}}+(\hbar c)^{2}\nabla ^{2}\psi =(mc^{2})^{2}\psi \,,}$$

(1)

by inserting the energy operator and momentum operator into the relativistic energy–momentum relation:

$${\displaystyle E^{2}-(pc)^{2}=(mc^{2})^{2}\,,}$$

(2)

The solutions to (1) are scalar fields. The KG equation is undesirable due to its prediction of negative energies and probabilities, as a result of the quadratic nature of (2) – inevitable in a relativistic theory. This equation was initially proposed by Schrödinger, and he discarded it for such reasons, only to realize a few months later that its non-relativistic limit (what is now called the Schrödinger equation) was still of importance. Nevertheless, (1) is applicable to spin-0 bosons.^[3]

Neither the non-relativistic nor relativistic equations found by Schrödinger could predict the fine structure in the Hydrogen spectral series. The mysterious underlying property was spin. The first two-dimensional spin matrices (better known as the Pauli matrices) were introduced by Pauli in the Pauli equation; the Schrödinger equation with a non-relativistic Hamiltonian including an extra term for particles in magnetic fields, but this was phenomenological. Weyl found a relativistic equation in terms of the Pauli matrices; the Weyl equation, for massless spin-⁠1/2⁠ fermions. The problem was resolved by Dirac in the late 1920s, when he furthered the application of equation (2) to the electron – by various manipulations he factorized the equation into the form:

$${\displaystyle \left({\frac {E}{c}}-{\boldsymbol {\alpha }}\cdot \mathbf {p} -\beta mc\right)\left({\frac {E}{c}}+{\boldsymbol {\alpha }}\cdot \mathbf {p} +\beta mc\right)\psi =0\,,}$$

(3A)

and one of these factors is the Dirac equation (see below), upon inserting the energy and momentum operators. For the first time, this introduced new four-dimensional spin matrices α and β in a relativistic wave equation, and explained the fine structure of hydrogen. The solutions to (3A) are multi-component spinor fields, and each component satisfies (1). A remarkable result of spinor solutions is that half of the components describe a particle while the other half describe an antiparticle; in this case the electron and positron. The Dirac equation is now known to apply for all massive spin-⁠1/2⁠ fermions. In the non-relativistic limit, the Pauli equation is recovered, while the massless case results in the Weyl equation.

Although a landmark in quantum theory, the Dirac equation is only true for spin-⁠1/2⁠ fermions, and still predicts negative energy solutions, which caused controversy at the time (in particular – not all physicists were comfortable with the "Dirac sea" of negative energy states).

The natural problem became clear: to generalize the Dirac equation to particles with any spin; both fermions and bosons, and in the same equations their antiparticles (possible because of the spinor formalism introduced by Dirac in his equation, and then-recent developments in spinor calculus by van der Waerden in 1929), and ideally with positive energy solutions.^[2]

This was introduced and solved by Majorana in 1932, by a deviated approach to Dirac. Majorana considered one "root" of (3A):

$${\displaystyle \left({\frac {E}{c}}+{\boldsymbol {\alpha }}\cdot \mathbf {p} -\beta mc\right)\psi =0\,,}$$

(3B)

where ψ is a spinor field now with infinitely many components, irreducible to a finite number of tensors or spinors, to remove the indeterminacy in sign. The matrices α and β are infinite-dimensional matrices, related to infinitesimal Lorentz transformations. He did not demand that each component of 3B satisfy equation (2); instead he regenerated the equation using a Lorentz-invariant action, via the principle of least action, and application of Lorentz group theory.^[4][5]

Majorana produced other important contributions that were unpublished, including wave equations of various dimensions (5, 6, and 16). They were anticipated later (in a more involved way) by de Broglie (1934), and Duffin, Kemmer, and Petiau (around 1938–1939) see Duffin–Kemmer–Petiau algebra. The Dirac–Fierz–Pauli formalism was more sophisticated than Majorana's, as spinors were new mathematical tools in the early twentieth century, although Majorana's paper of 1932 was difficult to fully understand; it took Pauli and Wigner some time to understand it, around 1940.^[2]

Dirac in 1936, and Fierz and Pauli in 1939, built equations from irreducible spinors A and B, symmetric in all indices, for a massive particle of spin n + 1⁄2 for integer n (see Van der Waerden notation for the meaning of the dotted indices):

$${\displaystyle p_{\gamma {\dot {\alpha }}}A_{\epsilon _{1}\epsilon _{2}\cdots \epsilon _{n}}^{{\dot {\alpha }}{\dot {\beta }}_{1}{\dot {\beta }}_{2}\cdots {\dot {\beta }}_{n}}=mcB_{\gamma \epsilon _{1}\epsilon _{2}\cdots \epsilon _{n}}^{{\dot {\beta }}_{1}{\dot {\beta }}_{2}\cdots {\dot {\beta }}_{n}}}$$

(4A)

$${\displaystyle p^{\gamma {\dot {\alpha }}}B_{\gamma \epsilon _{1}\epsilon _{2}\cdots \epsilon _{n}}^{{\dot {\beta }}_{1}{\dot {\beta }}_{2}\cdots {\dot {\beta }}_{n}}=mcA_{\epsilon _{1}\epsilon _{2}\cdots \epsilon _{n}}^{{\dot {\alpha }}{\dot {\beta }}_{1}{\dot {\beta }}_{2}\cdots {\dot {\beta }}_{n}}}$$

(4B)

where p is the momentum as a covariant spinor operator. For n = 0, the equations reduce to the coupled Dirac equations and A and B together transform as the original Dirac spinor. Eliminating either A or B shows that A and B each fulfill (1).^[2] The direct derivation of the Dirac-Pauli-Fierz equations using the Bargmann-Wigner operators is given in.^[6]

In 1941, Rarita and Schwinger focussed on spin-3⁄2 particles and derived the Rarita–Schwinger equation, including a Lagrangian to generate it, and later generalized the equations analogous to spin n + 1⁄2 for integer n. In 1945, Pauli suggested Majorana's 1932 paper to Bhabha, who returned to the general ideas introduced by Majorana in 1932. Bhabha and Lubanski proposed a completely general set of equations by replacing the mass terms in (3A) and (3B) by an arbitrary constant, subject to a set of conditions which the wave functions must obey.^[7]

Finally, in the year 1948 (the same year as Feynman's path integral formulation was cast), Bargmann and Wigner formulated the general equation for massive particles which could have any spin, by considering the Dirac equation with a totally symmetric finite-component spinor, and using Lorentz group theory (as Majorana did): the Bargmann–Wigner equations.^[2][8] In the early 1960s, a reformulation of the Bargmann–Wigner equations was made by H. Joos and Steven Weinberg, the Joos–Weinberg equation. Various theorists at this time did further research in relativistic Hamiltonians for higher spin particles.^[1][9][10]

The relativistic description of spin particles has been a difficult problem in quantum theory. It is still an area of the present-day research because the problem is only partially solved; including interactions in the equations is problematic, and paradoxical predictions (even from the Dirac equation) are still present.^[5]

The following equations have solutions which satisfy the superposition principle, that is, the wave functions are additive.

Throughout, the standard conventions of tensor index notation and Feynman slash notation are used, including Greek indices which take the values 1, 2, 3 for the spatial components and 0 for the timelike component of the indexed quantities. The wave functions are denoted ψ, and ∂_μ are the components of the four-gradient operator.

In matrix equations, the Pauli matrices are denoted by σ^μ in which μ = 0, 1, 2, 3, where σ⁰ is the 2 × 2 identity matrix: $${\displaystyle \sigma ^{0}={\begin{pmatrix}1&0\\0&1\\\end{pmatrix}}}$$ and the other matrices have their usual representations. The expression $${\displaystyle \sigma ^{\mu }\partial _{\mu }\equiv \sigma ^{0}\partial _{0}+\sigma ^{1}\partial _{1}+\sigma ^{2}\partial _{2}+\sigma ^{3}\partial _{3}}$$ is a 2 × 2 matrix operator which acts on 2-component spinor fields.

The gamma matrices are denoted by γ^μ, in which again μ = 0, 1, 2, 3, and there are a number of representations to select from. The matrix γ⁰ is not necessarily the 4 × 4 identity matrix. The expression $${\displaystyle i\hbar \gamma ^{\mu }\partial _{\mu }+mc\equiv i\hbar (\gamma ^{0}\partial _{0}+\gamma ^{1}\partial _{1}+\gamma ^{2}\partial _{2}+\gamma ^{3}\partial _{3})+mc{\begin{pmatrix}1&0&0&0\\0&1&0&0\\0&0&1&0\\0&0&0&1\end{pmatrix}}}$$ is a 4 × 4 matrix operator which acts on 4-component spinor fields.

Note that terms such as "mc" scalar multiply an identity matrix of the relevant dimension, the common sizes are 2 × 2 or 4 × 4, and are conventionally not written for simplicity.

Particle spin quantum number s	Name	Equation	Typical particles the equation describes
0	Klein–Gordon equation	${\displaystyle (\hbar \partial _{\mu }+imc)(\hbar \partial ^{\mu }-imc)\psi =0}$	Massless or massive spin-0 particle (such as Higgs bosons).
1/2	Weyl equation	${\displaystyle \sigma ^{\mu }\partial _{\mu }\psi =0}$	Massless spin-1/2 particles.
	Dirac equation	${\displaystyle \left(i\hbar \partial \!\!\!/-mc\right)\psi =0}$	Massive spin-1/2 particles (such as electrons).
	Two-body Dirac equations	${\displaystyle [(\gamma _{1})_{\mu }(p_{1}-{\tilde {A}}_{1})^{\mu }+m_{1}+{\tilde {S}}_{1}]\Psi =0,}$ ${\displaystyle [(\gamma _{2})_{\mu }(p_{2}-{\tilde {A}}_{2})^{\mu }+m_{2}+{\tilde {S}}_{2}]\Psi =0.}$	Massive spin-1/2 particles (such as electrons).
	Majorana equation	${\displaystyle i\hbar \partial \!\!\!/\psi -mc\psi _{c}=0}$	Massive Majorana particles.
	Breit equation	${\displaystyle i\hbar {\frac {\partial \Psi }{\partial t}}=\left(\sum _{i}{\hat {H}}_{D}(i)+\sum _{i>j}{\frac {1}{r_{ij}}}-\sum _{i>j}{\hat {B}}_{ij}\right)\Psi }$	Two massive spin-1/2 particles (such as electrons) interacting electromagnetically to first order in perturbation theory.
1	Maxwell's equations (in QED using the Lorenz gauge)	${\displaystyle \partial _{\mu }\partial ^{\mu }A^{\nu }=e{\overline {\psi }}\gamma ^{\nu }\psi }$	Photons, massless spin-1 particles.
	Proca equation	${\displaystyle \partial _{\mu }(\partial ^{\mu }A^{\nu }-\partial ^{\nu }A^{\mu })+\left({\frac {mc}{\hbar }}\right)^{2}A^{\nu }=0}$	Massive spin-1 particle (such as W and Z bosons).
3/2	Rarita–Schwinger equation	${\displaystyle \epsilon ^{\mu \nu \rho \sigma }\gamma ^{5}\gamma _{\nu }\partial _{\rho }\psi _{\sigma }+m\psi ^{\mu }=0}$	Massive spin-3/2 particles.
s	Bargmann–Wigner equations	${\displaystyle {\begin{aligned}(-i\hbar \gamma ^{\mu }\partial _{\mu }+mc)_{\alpha _{1}\alpha _{1}'}\psi _{\alpha '_{1}\alpha _{2}\alpha _{3}\cdots \alpha _{2s}}&=0\\(-i\hbar \gamma ^{\mu }\partial _{\mu }+mc)_{\alpha _{2}\alpha _{2}'}\psi _{\alpha _{1}\alpha '_{2}\alpha _{3}\cdots \alpha _{2s}}&=0\\&\;\;\vdots \\(-i\hbar \gamma ^{\mu }\partial _{\mu }+mc)_{\alpha _{2s}\alpha '_{2s}}\psi _{\alpha _{1}\alpha _{2}\alpha _{3}\cdots \alpha '_{2s}}&=0\end{aligned}}}$ where ψ is a rank-2s 4-component spinor.	Free particles of arbitrary spin (bosons and fermions).[9][11]
	Joos–Weinberg equation	${\displaystyle [(i\hbar )^{2s}\gamma ^{\mu _{1}\mu _{2}\cdots \mu _{2s}}\partial _{\mu _{1}}\partial _{\mu _{2}}\cdots \partial _{\mu _{2s}}+(mc)^{2s}]\psi =0}$	Free particles of arbitrary spin (bosons and fermions).

The Duffin–Kemmer–Petiau equation is an alternative equation for spin-0 and spin-1 particles: $${\displaystyle (i\hbar \beta ^{a}\partial _{a}-mc)\psi =0}$$

Start with the standard special relativity (SR) 4-vectors

• 4-position ${\displaystyle X^{\mu }=\mathbf {X} =(ct,{\vec {\mathbf {x} }})}$
• 4-velocity ${\displaystyle U^{\mu }=\mathbf {U} =\gamma (c,{\vec {\mathbf {u} }})}$
• 4-momentum ${\displaystyle P^{\mu }=\mathbf {P} =\left({\frac {E}{c}},{\vec {\mathbf {p} }}\right)}$
• 4-wavevector ${\displaystyle K^{\mu }=\mathbf {K} =\left({\frac {\omega }{c}},{\vec {\mathbf {k} }}\right)}$
• 4-gradient ${\displaystyle \partial ^{\mu }=\mathbf {\partial } =\left({\frac {\partial _{t}}{c}},-{\vec {\mathbf {\nabla } }}\right)}$

Note that each 4-vector is related to another by a Lorentz scalar:

• ${\displaystyle \mathbf {U} ={\frac {d}{d\tau }}\mathbf {X} }$, where ${\displaystyle \tau }$ is the proper time
• ${\displaystyle \mathbf {P} =m_{0}\mathbf {U} }$, where ${\displaystyle m_{0}}$ is the rest mass
• ${\displaystyle \mathbf {K} =(1/\hbar )\mathbf {P} }$, which is the 4-vector version of the Planck–Einstein relation & the de Broglie matter wave relation
• ${\displaystyle \mathbf {\partial } =-i\mathbf {K} }$, which is the 4-gradient version of complex-valued plane waves

Now, just apply the standard Lorentz scalar product rule to each one:

• ${\displaystyle \mathbf {U} \cdot \mathbf {U} =(c)^{2}}$
• ${\displaystyle \mathbf {P} \cdot \mathbf {P} =(m_{0}c)^{2}}$
• ${\displaystyle \mathbf {K} \cdot \mathbf {K} =\left({\frac {m_{0}c}{\hbar }}\right)^{2}}$
• ${\displaystyle \mathbf {\partial } \cdot \mathbf {\partial } =\left({\frac {-im_{0}c}{\hbar }}\right)^{2}=-\left({\frac {m_{0}c}{\hbar }}\right)^{2}}$

The last equation is a fundamental quantum relation.

When applied to a Lorentz scalar field ${\displaystyle \psi }$, one gets the Klein–Gordon equation, the most basic of the quantum relativistic wave equations.

• ${\displaystyle \left[\mathbf {\partial } \cdot \mathbf {\partial } +\left({\frac {m_{0}c}{\hbar }}\right)^{2}\right]\psi =0}$: in 4-vector format
• ${\displaystyle \left[\partial _{\mu }\partial ^{\mu }+\left({\frac {m_{0}c}{\hbar }}\right)^{2}\right]\psi =0}$: in tensor format
• ${\displaystyle \left[(\hbar \partial _{\mu }+im_{0}c)(\hbar \partial ^{\mu }-im_{0}c)\right]\psi =0}$: in factored tensor format

The Schrödinger equation is the low-velocity limiting case (v ≪ c) of the Klein–Gordon equation.

When the relation is applied to a four-vector field ${\displaystyle A^{\mu }}$ instead of a Lorentz scalar field ${\displaystyle \psi }$, then one gets the Proca equation (in Lorenz gauge): $${\displaystyle \left[\mathbf {\partial } \cdot \mathbf {\partial } +\left({\frac {m_{0}c}{\hbar }}\right)^{2}\right]A^{\mu }=0}$$

If the rest mass term is set to zero (light-like particles), then this gives the free Maxwell equation (in Lorenz gauge) $${\displaystyle [\mathbf {\partial } \cdot \mathbf {\partial } ]A^{\mu }=0}$$

Under a proper orthochronous Lorentz transformation x → Λx in Minkowski space, all one-particle quantum states ψ^j_σ of spin j with spin z-component σ locally transform under some representation D of the Lorentz group:^[12][13] $${\displaystyle \psi (x)\rightarrow D(\Lambda )\psi (\Lambda ^{-1}x)}$$ where D(Λ) is some finite-dimensional representation, i.e. a matrix. Here ψ is thought of as a column vector containing components with the allowed values of σ. The quantum numbers j and σ as well as other labels, continuous or discrete, representing other quantum numbers are suppressed. One value of σ may occur more than once depending on the representation. Representations with several possible values for j are considered below.

The irreducible representations are labeled by a pair of half-integers or integers (A, B). From these all other representations can be built up using a variety of standard methods, like taking tensor products and direct sums. In particular, space-time itself constitutes a 4-vector representation (⁠1/2⁠, ⁠1/2⁠) so that Λ ∈ D^{(1/2, 1/2)}. To put this into context; Dirac spinors transform under the (⁠1/2⁠, 0) ⊕ (0, ⁠1/2⁠) representation. In general, the (A, B) representation space has subspaces that under the subgroup of spatial rotations, SO(3), transform irreducibly like objects of spin j, where each allowed value: $${\displaystyle j=A+B,A+B-1,\dots ,|A-B|,}$$ occurs exactly once.^[14] In general, tensor products of irreducible representations are reducible; they decompose as direct sums of irreducible representations.

The representations D^{(j, 0)} and D^{(0, j)} can each separately represent particles of spin j. A state or quantum field in such a representation would satisfy no field equation except the Klein–Gordon equation.

There are equations which have solutions that do not satisfy the superposition principle.

• Yang–Mills equation: describes a non-abelian gauge field
• Yang–Mills–Higgs equations: describes a non-abelian gauge field coupled with a massive spin-0 particle

• Einstein field equations: describe interaction of matter with the gravitational field (massless spin-2 field): $${\displaystyle R_{\mu \nu }-{\frac {1}{2}}g_{\mu \nu }\,R+g_{\mu \nu }\Lambda ={\frac {8\pi G}{c^{4}}}T_{\mu \nu }}$$ The solution is a metric tensor field, rather than a wave function.

• List of equations in nuclear and particle physics
• List of equations in quantum mechanics
• Lorentz transformation
• Mathematical descriptions of the electromagnetic field
  • Quantization of the electromagnetic field
• Minimal coupling
• Scalar field theory
• Status of special relativity