Dirac's factorization approach was aimed at obtaining a linear first-order differential equation suitable for spin-1/2 particles, contrasting the second-order Klein-Gordon equation, which couldn't adequately describe such particles. Dirac proposed an equation ((iγµ∂µ −m)ψ = 0) that when squared resulted in the Klein-Gordon equation ((iγµ∂µ −m)(iγν∂ν + m)ψ = 0). This required the introduction of gamma matrices, which satisfy the Clifford algebra {γµ, γν} = 2gµνI4, ensuring that when the Dirac equation is squared, it recovers the original Klein-Gordon form. This factorization thus bridges the gap between quantum mechanics and relativity while accommodating fermionic particles.

Dirac's motivation for introducing a first-order differential equation for spin-1/2 particles was to resolve the inadequacies of the Klein-Gordon equation, which was second-order and did not appropriately describe fermions like electrons. He sought an equation consistent both with quantum theory and special relativity, which required maintaining Lorentz invariance and accommodating the spin characteristics intrinsic to fermions. This necessitated a linear equation with respect to time and spatial derivatives, which the Dirac equation achieved, bridging a critical gap between quantum mechanics and relativistic physics.

Plane wave solutions of the Dirac equation, represented as ψ(x) = u(p)e−ip·x, provide crucial insights into the dual characteristics of particles and antiparticles. The term u(p) is a Dirac spinor that must satisfy the algebraic equation (γµpµ −m)u(p) = 0. For positive energy solutions, u(p) describes particles such as electrons, characterized by specific energy and momentum relations. Negative energy solutions, v(p), denote antiparticles like positrons. These solutions illustrate the intrinsic link between particles and their antiparticles within the Dirac framework, reinforced by the symmetry and algebraic properties of the equation that naturally predict antimatter.

The covariant form of the Dirac equation (iγµ∂µ −m)ψ = 0 is crucial for preserving Lorentz invariance, a key feature in relativity, allowing the equation to be applied consistently across different inertial frames. The Lagrangian representation, ¯ψ(iγµ∂µ −m)ψ = 0, formalizes it within the framework of field theory, enabling its use in constructing dynamical models. This representation is the starting point for deriving Feynman rules in quantum electrodynamics, thus making the equation integral to particle physics and theoretical models describing fundamental interactions. It provides a robust, systematic approach for analyzing relativistic particles and their interactions, broadening its application in theoretical physics.

Dirac spinors represent solutions to the Dirac equation for spin-1/2 particles. In a plane wave solution ψ(x) = u(p)e−ip·x, u(p) is the Dirac spinor. For positive energy solutions, it takes the form u(p) = (ϕ, (σ·p)/(E+m)ϕ)T, where ϕ is a 2-component spinor, and σ represents the Pauli matrices. These solutions describe electrons. Conversely, negative energy solutions v(p) correspond to antiparticles of electrons. Dirac spinors are thus essential in depicting the quantum states of particles and antiparticles within quantum electrodynamics, showing the duality of matter and antimatter as predicted by the Dirac equation.

The Dirac equation, formulated by Paul Dirac in 1928, is significant because it was the first equation to successfully unify quantum mechanics and special relativity for spin-1/2 particles, leading to a consistent framework for describing electron behavior at relativistic speeds. Its formulation as a first-order linear differential equation ensured compatibility with Lorentz invariance and quantum principles. This unification inherently predicted the existence of antimatter as solutions to the equation involve both positive and negative energy states, the latter corresponding to antiparticles, or antimatter.

The Dirac equation has several profound physical implications. It predicts the existence of antimatter by allowing for solutions with negative energy states, which correspond to antiparticles like positrons for electrons. It also naturally describes fermionic particles, specifically spin-1/2 particles like electrons, integrating their intrinsic spin properties into a relativistic framework. Furthermore, it ensures probability conservation through a conserved current jµ = ¯ψγµψ, and forms the theoretical foundation for quantum electrodynamics (QED), which describes electromagnetic interactions at the quantum level, reinforcing the broader implications of the Dirac equation in particle physics.

Gamma matrices are crucial in the Dirac equation as they allow for the linearization of the relativistic energy-momentum relation, resulting in a first-order differential equation. They are 4x4 matrices that satisfy the anticommutation relation {γµ, γν} = 2gµνI4, forming the Clifford algebra. This relation is fundamental for the equation's transformation properties under Lorentz transformations, ensuring its covariance. Each gamma matrix is associated with the transformation of spinor fields and, along with the Minkowski metric gµν = diag(1, −1, −1, −1), they maintain the relativistic invariance of the equation.

The gamma matrices used in the Dirac equation must satisfy the anticommutation relation {γµ, γν} = 2gµνI4, which defines the Clifford algebra. This mathematical condition ensures that when the Dirac equation is multiplied by its conjugate, the resulting expression is identical to the Klein-Gordon equation. Here, gµν is the Minkowski metric, diag(1, −1, −1, −1). This ensures that the product of the Dirac equation terms reflects the energy-momentum relation required for relativistic wave equations, maintaining their consistency and compatibility, thereby ensuring covariance under Lorentz transformations.

The Dirac equation ensures probability conservation through the existence of a conserved current, jµ = ¯ψγµψ. This current is constructed from the Dirac spinor ψ and its adjoint ¯ψ = ψ†γ0, ensuring that the probability density and current together obey the continuity equation. This reflects a conservation law for probability, which is a fundamental principle in quantum mechanics. The Lorentz invariance of the Dirac equation guarantees that this probability conservation holds under relativistic transformations, providing a consistent and unified framework for quantum mechanics and special relativity.