Statement

Let DM=γμμD_{\mathrm M}=\gamma^\mu\partial_\mu be the with

{γμ,γν}=2ημνI,η=ημνμν.\{\gamma^\mu,\gamma^\nu\}=-2\eta^{\mu\nu}I, \qquad \Box_\eta=-\eta^{\mu\nu}\partial_\mu\partial_\nu.

For every constant mm,

(iDMm)(iDM+m)=(iDM+m)(iDMm)=(η+m2)I.(iD_{\mathrm M}-m)(iD_{\mathrm M}+m) =(iD_{\mathrm M}+m)(iD_{\mathrm M}-m) =-(\Box_\eta+m^2)I.

Consequently every solution of either massive free satisfies the componentwise.

Proof

Commutativity of the constant mm with DMD_{\mathrm M} cancels the cross terms. The Clifford relation and commutativity of partial derivatives give

DM2=12{γμ,γν}μν=ημνμνI=ηI.D_{\mathrm M}^{\,2} =\frac12\{\gamma^\mu,\gamma^\nu\}\partial_\mu\partial_\nu =-\eta^{\mu\nu}\partial_\mu\partial_\nu I =\Box_\eta I.

Therefore

(iDM)2m2=DM2m2=(η+m2)I.(iD_{\mathrm M})^2-m^2=-D_{\mathrm M}^{\,2}-m^2 =-(\Box_\eta+m^2)I.
What the implication does not say

The converse is false without additional data: a tuple of Klein–Gordon solutions need not satisfy a first-order Dirac equation. The factorization selects spinor solutions obeying an extra linear constraint.

Curvature and gauge warning

On a curved spin manifold, the square of the contains a scalar-curvature term; coupling to a gauge field adds Clifford contraction of its curvature. Thus the free flat identity acquires lower-order corrections. The principal part remains a normally hyperbolic connection operator, but one cannot replace it by the scalar g\Box_g componentwise in an arbitrary frame.

References
  1. H. Blaine Lawson Jr. and Marie-Louise Michelsohn, Spin Geometry, Princeton University Press, 1989. Publisher record. Relevant: Chapter II, §8.
  2. Christian Bär, Nicolas Ginoux, and Frank Pfäffle, Wave Equations on Lorentzian Manifolds and Quantization, European Mathematical Society, 2007. Publisher record. Relevant: §§1.3 and 3.4.