eoghan said:
But I could easily come up with a Lorentz invariant equation that is first order, e.g.
$$
(M^\mu\partial_\mu + m^2)\phi=0
$$
where M is a generic matrix.
Now, something should be wrong with this equation
How about, almost everything is wrong with that equation:
1) If, as you say, [itex]M^{\mu}[/itex] is not [itex]\partial^{\mu}[/itex] but a “generic matrix”, then the expression [itex]M^{\mu}\partial_{\mu} + m^{2}[/itex] is
meaningless because each term has different physical unit. In the natural units, the dimension of the first term is [itex]\mbox{cm}^{-1}[/itex] while the dimension of [itex]m^{2}[/itex] is [itex]\mbox{cm}^{-2}[/itex].
2) You said that [itex]\phi[/itex] is a scalar field. In 4-dimensional spacetime, a (real) scalar field can be described
either by a single function
or (equivalently) by a 5-component field treating [itex](\phi , \partial_{\mu}\phi )[/itex] as independent variables. In the first case (i.e., when [itex]\phi[/itex] is a 1-component field), your “matrices” [itex]M^{\mu}[/itex] must be [itex]1 \times 1[/itex] matrices. Thus, you must take [itex]M^{\mu} = \partial^{\mu}[/itex] so that the correct dispersion relation [itex]E^{2} = P^{2} + m^{2}[/itex] holds. In the second (5-component) case, the correct first-order equation (for a real scalar field) looks exactly like Dirac equation [tex]\left( i \Gamma^{\mu} \partial_{\mu} - m \right) \Psi (x) = 0 , \ \ \ \ \ \ \ \ \ \ (1)[/tex] with [tex]\Psi = \left( \varphi , \psi_{0} , \psi_{1}, \psi_{2} , \psi_{3} \right)^{T} ,[/tex] and the [itex]\Gamma[/itex]’s are a set of four [itex]5 \times 5[/itex] matrices satisfying the Duffin-Kemmer algebra [tex]\Gamma^{\mu}\Gamma^{\rho}\Gamma^{\nu} + \Gamma^{\nu}\Gamma^{\rho}\Gamma^{\mu} = \eta^{\mu \rho}\Gamma^{\nu} + \eta^{\nu\rho}\Gamma^{\mu} . \ \ \ \ \ \ (2)[/tex] Notice that (the Duffin-Kemmer equation) Eq(1) has no [itex]m^{2}[/itex] term in it. Within the representation theory, the Duffin-Kemmer equation has beautiful interpretation. However, since your knowledge in group theory is (unfortunately) poor, bellow I will only show you how to obtain the DK equation (1) from the following Klein-Gordon equation of real (1-component) scalar field [tex]\partial^{\mu}\partial_{\mu} \phi (x) + m^{2} \phi (x) = 0 . \ \ \ \ \ \ \ \ \ \ \ \ \ (3)[/tex] So we want to reduce (3) into a system of five (coupled) first-order equations. To do this we define a scalar field by [tex]\varphi (x) = \sqrt{m} \phi (x) ,[/tex] and a 4-vector field by [tex]\psi_{\mu}(x) = \frac{1}{\sqrt{m}} \partial_{\mu}\phi (x) .[/tex] Thus, the KG equation (3) can be replaced by the following equivalent system of first-order equations [tex]\partial^{\mu}\psi_{\mu} (x) + m \varphi (x) = 0 ,[/tex][tex]\partial_{\mu} \varphi (x) - m \psi_{\mu}(x) = 0 .[/tex] These can easily be rewritten as matrix equation
[tex]\begin{pmatrix} - m & - \partial_{0} & \partial_{1} & \partial_{2} & \partial_{3} \\ \partial_{0} & - m & 0 & 0 & 0 \\ \partial_{1} & 0 & - m & 0 & 0 \\ \partial_{2} & 0 & 0 & - m & 0 \\ \partial_{3} & 0 & 0 & 0 & - m \end{pmatrix} \begin{pmatrix} \varphi (x) \\ \psi_{0}(x) \\ \psi_{1}(x) \\ \psi_{2}(x) \\ \psi_{3}(x) \end{pmatrix} = 0 .[/tex] This is exactly the Duffin-Kemmer equation (1) with the [itex]\Gamma[/itex]’s given by [tex]\Gamma^{0} = \begin{pmatrix} 0 & i & 0 & 0 & 0 \\ - i & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \end{pmatrix} , \ \ \ \Gamma^{1} = \begin{pmatrix} 0 & 0 & - i & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ - i & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \end{pmatrix} ,[/tex] [tex]\Gamma^{2} = \begin{pmatrix} 0 & 0 & 0 & - i & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ - i & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \end{pmatrix} , \ \ \ \Gamma^{3} = \begin{pmatrix} 0 & 0 & 0 & 0 & - i \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ - i & 0 & 0 & 0 & 0 \end{pmatrix} ,[/tex] and the Duffin-Kemmer field [itex]\Psi (x)[/itex] is given by [tex]\Psi (x) = \begin{pmatrix} \varphi (x) \\ \psi_{0}(x) \\ \psi_{1}(x) \\ \psi_{2}(x) \\ \psi_{3}(x) \end{pmatrix} .[/tex]
Similarly, you can equivalently write the Proca equation [tex]\partial_{\mu}F^{\mu\nu} + m^{2}A^{\nu} = 0 ,[/tex] which describes a free massive spin-1 field, as a Duffin-Kemmer equation, with four [itex]10 \times 10[/itex] [itex]\Gamma[/itex]'s. I leave this as exercise for you.