# Published in Foundations of Physics Letters, vol. 12, pp. 291-8 (1998) Maxwell equations as the one-photon quantum equation

## Abstract

Maxwell equations (Faraday and Ampere-Maxwell laws) can be presented as a three component equation in a way similar to the two component neutrino equation. However, in this case, the electric and magnetic Gauss’s laws can not be derived from first principles. We have shown how all Maxwell equations can be derived simultaneously from first principles, similar to those which have been used to derive the Dirac relativistic electron equation. We have also shown that equations for massless particles, derived by Dirac in 1936, lead to the same result. The complex wave function, being a linear combination of the electric and magnetic fields, is a locally measurable and well understood quantity. Therefore Maxwell equations should be used as a guideline for proper interpretations of quantum theories.

The Maxwell equations (except for the electric and magnetic Gauss’s law) can be presented by a three component equation in a way similar to the two component neutrino equation. This was already known to Oppenheimer  and to Majorana , . Also this type of equation is a particular case of a more general equation for any spin derived by Weinberg . There is a continuous interest in this equation even to this day , , , . However one of the drawbacks of the above derivations is that the electric and magnetic Gauss’s laws are not derived from first principles.

The aim of the present latter is to complement the above mentioned works, and to derive all Maxwell equations directly from a decomposition similar to that which was used to derive the Dirac relativistic electron equation.

The Dirac equation is derived from the relativistic condition on the Energy $`E,`$ mass $`m,`$ and momentum $`\stackrel{}{𝐩}`$:

$$\left(E^2c^2\stackrel{}{𝐩}^2m^2c^4\right)I^{(4)}\mathrm{\Psi }=0,$$
(1)   
where $`I^{(4)}`$ is the $`4\times 4`$ unit matrix and $`\mathrm{\Psi }`$ is a four component column (bispinor) wave function. Eq. (1) is decomposed into

$$\left[EI^{(4)}+\left(\begin{array}{cc}mc^2I^{(2)}& c\stackrel{}{𝐩}\stackrel{}{\sigma }\\ c\stackrel{}{𝐩}\stackrel{}{\sigma }& mc^2I^{(2)}\end{array}\right)\right]\left[EI^{(4)}\left(\begin{array}{cc}mc^2I^{\left(2\right)}& c\stackrel{}{𝐩}\stackrel{}{\sigma }\\ c\stackrel{}{𝐩}\stackrel{}{\sigma }& mc^2I^{\left(2\right)}\end{array}\right)\right]\mathrm{\Psi }=0,$$
(2)   
where $`I^{\left(2\right)}`$ is the $`2\times 2`$ unit matrix and $`\stackrel{}{\sigma }`$ is the Pauli spin one-half vector matrix with the components

$$\sigma _x=\left(\begin{array}{cc}0& 1\\ 1& 0\end{array}\right),\sigma _y=\left(\begin{array}{cc}0& i\\ i& 0\end{array}\right),\sigma _z=\left(\begin{array}{cc}1& 0\\ 0& 1\end{array}\right),I^{\left(2\right)}=\left(\begin{array}{cc}1& 0\\ 0& 1\end{array}\right).$$
(3)   
The two component neutrino equation can be derived from the decomposition

$$\left(E^2c^2\stackrel{}{𝐩}^2\right)I^{\left(2\right)}\psi =\left[EI^{\left(2\right)}c\stackrel{}{𝐩}\stackrel{}{\sigma }\right]\left[EI^{\left(2\right)}+c\stackrel{}{𝐩}\stackrel{}{\sigma }\right]\psi =0,$$
(4)   
where $`\psi `$ is a two component spinor wavefunction.

We shall derive the photon equation from the following decomposition

$$\left(\frac{E^2}{c^2}\stackrel{}{𝐩}^2\right)I^{\left(3\right)}=\left(\frac{E}{c}I^{\left(3\right)}\stackrel{}{𝐩}\stackrel{}{𝐒}\right)\left(\frac{E}{c}I^{\left(3\right)}+\stackrel{}{𝐩}\stackrel{}{𝐒}\right)\left(\begin{array}{ccc}p_x^2& p_xp_y& p_xp_z\\ p_yp_x& p_y^2& p_yp_z\\ p_zp_x& p_zp_y& p_z^2\end{array}\right)=0,$$
(5)   
where $`I^{\left(3\right)}`$ is a $`3\times 3`$ unit matrix, and $`\stackrel{}{𝐒}`$ is a spin one vector matrix with components

$$S_x=\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& i\\ 0& i& 0\end{array}\right),S_y=\left(\begin{array}{ccc}0& 0& i\\ 0& 0& 0\\ i& 0& 0\end{array}\right),S_z=\left(\begin{array}{ccc}0& i& 0\\ i& 0& 0\\ 0& 0& 0\end{array}\right),I^{\left(3\right)}=\left(\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\end{array}\right),$$
(6)   
and with the properties

$$[S_x,S_y]=iS_z,[S_z,S_x]=iS_y,[S_y,S_z]=iS_x,\stackrel{}{𝐒}^2=2I^{\left(3\right)}.$$
(7)   
The decomposition (5) can be verified directly by substitution.

It will be crucial to note that the matrix on the right hand side of Eq. (5) can be rewritten as:

$$\left(\begin{array}{ccc}p_x^2& p_xp_y& p_xp_z\\ p_yp_x& p_y^2& p_yp_z\\ p_zp_x& p_zp_y& p_z^2\end{array}\right)=\left(\begin{array}{c}p_x\\ p_y\\ p_z\end{array}\right)\left(\begin{array}{ccc}p_x& p_y& p_z\end{array}\right).$$
(8)   
From Eqs. $`\left(\text{5}\text{6}\right)`$ and $`\left(\text{8}\right)`$, the photon equation can be obtained form

$$\left(\frac{E^2}{c^2}\stackrel{}{𝐩}^2\right)\stackrel{}{𝚿}=\left(\frac{E}{c}I^{\left(3\right)}\stackrel{}{𝐩}\stackrel{}{𝐒}\right)\left(\frac{E}{c}I^{\left(3\right)}+\stackrel{}{𝐩}\stackrel{}{𝐒}\right)\stackrel{}{𝚿}\left(\begin{array}{c}p_x\\ p_y\\ p_z\end{array}\right)\left(\stackrel{}{𝐩}\stackrel{}{𝚿}\right)=0,$$
(9)   
where $`\stackrel{}{𝚿}`$ is a 3 component (column) wave function. Eq. (9) will be satisfied if the two equations

$`\left({\displaystyle \frac{E}{c}}I^{\left(3\right)}+\stackrel{}{𝐩}\stackrel{}{𝐒}\right)\stackrel{}{𝚿}`$ $`=`$ $`0,`$ (10)   
$`\stackrel{}{𝐩}\stackrel{}{𝚿}`$ $`=`$ $`0,`$ (11)   
will be simultaneously satisfied. For real energies and momenta complex conjugation of Eqs. (9) and (6) leads to

$$\left(\frac{E^2}{c^2}\stackrel{}{𝐩}^2\right)\stackrel{}{𝚿}^{}=\left(\frac{E}{c}I^{\left(3\right)}+\stackrel{}{𝐩}\stackrel{}{𝐒}\right)\left(\frac{E}{c}I^{\left(3\right)}\stackrel{}{𝐩}\stackrel{}{𝐒}\right)\stackrel{}{𝚿}^{}\left(\begin{array}{c}p_x\\ p_y\\ p_z\end{array}\right)\left(\stackrel{}{𝐩}\stackrel{}{𝚿}^{}\right)=0,$$
(12)   
where $`\stackrel{}{𝚿}^{}`$ is the complex conjugate of $`\stackrel{}{𝚿}.`$ Eq. (12) will be satisfied if the two equations

$`\left({\displaystyle \frac{E}{c}}I^{\left(3\right)}\stackrel{}{𝐩}\stackrel{}{𝐒}\right)\stackrel{}{𝚿}^{}`$ $`=`$ $`0,`$ (13)   
$`\stackrel{}{𝐩}\stackrel{}{𝚿}^{}`$ $`=`$ $`0,`$ (14)   
will be simultaneously satisfied. Eqs. (9) and (12) are the two different possible decompositions of their left hand side. Eqs. (12-14) do not contain new information as they are only the complex conjugates of Eqs. (9-11). On the other hand the physical interpretation is different, namely Eq. (10) is the negative helicity equation, while Eq. (13) is the positive helicity equation. It will be interesting to note that also other set of equivalent equations is possible. Eqs. (9) and (11) can be rewritten as

$$\left(\begin{array}{c}p_x\\ p_y\\ p_z\end{array}\right)\left(\stackrel{}{𝐩}\stackrel{}{𝚿}\right)=\left(\frac{E}{c}I^{\left(3\right)}\stackrel{}{𝐩}\stackrel{}{𝐒}\right)\left(\frac{E}{c}I^{\left(3\right)}+\stackrel{}{𝐩}\stackrel{}{𝐒}\right)\stackrel{}{𝚿}\left(\frac{E^2}{c^2}\stackrel{}{𝐩}^2\right)\stackrel{}{𝚿}=\mathrm{𝟎},$$
(15)   
which will be satisfied if the two equations

$$\left(\frac{E}{c}I^{\left(3\right)}+\stackrel{}{𝐩}\stackrel{}{𝐒}\right)\stackrel{}{𝚿}=0,\left(\frac{E^2}{c^2}\stackrel{}{𝐩}^2\right)\stackrel{}{𝚿}=\mathrm{𝟎},$$
(16)   
or their equivalents

$$\left(\frac{E}{c}I^{\left(3\right)}\stackrel{}{𝐩}\stackrel{}{𝐒}\right)\stackrel{}{𝚿}^{}=0,\left(\frac{E^2}{c^2}\stackrel{}{𝐩}^2\right)\stackrel{}{𝚿}^{}=\mathrm{𝟎},$$
(17)   
will be simultaneously satisfied. Maxwell equations will be derived from Eqs. (10-11).

We will show below that if in Eqs. (10) and (11) the quantum operator substitutions

$$Ei\mathrm{}\frac{}{t},\stackrel{}{𝐩}i\mathrm{},$$
(18)   
and the wavefunction substitution

$$\stackrel{}{𝚿}=\stackrel{}{𝐄}i\stackrel{}{𝐁},$$
(19)   
are made, as a result the Maxwell equations will be obtained. In Eq. (19) $`\stackrel{}{𝐄}`$ and $`\stackrel{}{𝐁}`$ are the electric and magnetic fields respectively. Indeed, one can easily check from Eqs. (6) and (18) that the following identity is satisfied

$$\left(\stackrel{}{𝐩}\stackrel{}{𝐒}\right)\stackrel{}{𝚿}=\mathrm{}\times \stackrel{}{𝚿}.$$
(20)   
From Eqs. (10-11) and (18, 20) we obtain

$$\frac{i\mathrm{}}{c}\frac{}{t}\stackrel{}{𝚿}=\mathrm{}\times \stackrel{}{𝚿},$$
(21)   
$$i\mathrm{}\stackrel{}{𝚿}=0.$$
(22)   
The constant $`\mathrm{}`$ can be cancelled out in Eqs. (21-22), and after replacing $`\stackrel{}{𝚿}`$ by Eq. (19), the following equations are obtained

$$\times \left(\stackrel{}{𝐄}i\stackrel{}{𝐁}\right)=i\frac{1}{c}\frac{\left(\stackrel{}{𝐄}i\stackrel{}{𝐁}\right)}{t},$$
(23)   
$$\left(\stackrel{}{𝐄}i\stackrel{}{𝐁}\right)=0.$$
(24)   
If in Eqs. (23-24) the electric and magnetic fields are real, the separation into the real and imaginary parts will lead to the Maxwell equations

$$\times \stackrel{}{𝐄}=\frac{1}{c}\frac{\stackrel{}{𝐁}}{t}$$
(25)   
$$\times \stackrel{}{𝐁}=\frac{1}{c}\frac{\stackrel{}{𝐄}}{t}$$
(26)   
$$\stackrel{}{𝐄}=0$$
(27)   
$$\stackrel{}{𝐁}=0.$$
(28)   
One should note that the Plank constant $`\mathrm{}`$ was cancelled out earlier, in Eqs. (21-22), which explains its absence in the Maxwell equations. Another comment should be made here, starting from equations (16) and (18-19) one can get equations which are equivalent to Maxwell equations (without sources), namely

$$\times \stackrel{}{𝐄}=\frac{1}{c}\frac{\stackrel{}{𝐁}}{t},\times \stackrel{}{𝐁}=\frac{1}{c}\frac{\stackrel{}{𝐄}}{t},\left(\frac{^2}{c^2t^2}^2\right)\stackrel{}{𝐄}=0,\left(\frac{^2}{c^2t^2}^2\right)\stackrel{}{𝐁}=0,$$
(29)   
the Gauss laws (27-28) are satisfied on the basis of Eq. (15).

Dirac  and Wigner , have derived relativistic equations for massless particles of any spin from which the Gauss laws, for the spin one case, can be derived. Moreover, Wigner  has shown that any finite-component massless field has only two possible helicity states. Dirac has derived equations for massless particles with spin $`k`$, which in the ordinary vector notation  are

$$\left\{kp_t+S_xp_x+S_yp_y+S_zp_z\right\}\psi =0,$$
(30)   
$$\left\{kp_x+S_xp_tiS_yp_z+iS_zp_y\right\}\psi =0,$$
(31)   
$$\left\{kp_y+S_yp_tiS_zp_x+iS_xp_z\right\}\psi =0,$$
(32)   
$$\left\{kp_z+S_zp_tiS_xp_y+iS_yp_x\right\}\psi =0,$$
(33)   
where the $`p_n`$ are the momenta, $`p_t=E/c`$, $`E`$ the energy, $`\psi `$ a $`\left(2k+1\right)`$ component wave function and $`S_n`$ are the spin $`\left(2k+1\right)\times \left(2k+1\right)`$ matrices which satisfy

$$[S_x,S_y]=iS_z,[S_z,S_x]=iS_y,[S_y,S_z]=iS_x,S_x^2+S_y^2+S_z^2=k(k+1)I^{\left(k\right)},$$
(34)   
and $`I^{\left(k\right)}`$ is a $`\left(2k+1\right)\times \left(2k+1\right)`$ unit matrix. As we shall see below, for the case $`k=1`$, Eq. (30) will lead to the Faraday and Ampere-Maxwell laws. The Gauss laws can be derived from Eqs. (30-33) in a way which will be described below. Eqs. (30-33) were analyzed extensively by Bacry , who derived them using Wigner’s condition  on the Pauli-Lubanski vector $`W^\mu `$ for massless fields

$$W^\mu =kp^\mu ,\mu =x,y,z,t.$$
(35)   
Let us now demonstrate how Eq. (11) can be derived from Eqs. (30-33). Following Dirac , one replaces Eqs. (30-33), which are linearly dependent, with the Eq. (30) and 3 conditions on the wave function, which are obtained by substituting $`p_t`$ from Eq. (30) into Eqs. (31-33)

$$\left\{kp_t+S_xp_x+S_yp_y+S_zp_z\right\}\psi =0,$$
(36)   
$$\left\{(k^2S_x^2)p_x+(ikS_zS_xS_y)p_y\left(ikS_y+S_xS_z\right)p_z\right\}\psi =0,$$
(37)   
$$\left\{(k^2S_y^2)p_y+(ikS_xS_yS_z)p_z\left(ikS_z+S_yS_x\right)p_x\right\}\psi =0,$$
(38)   
$$\left\{(k^2S_z^2)p_z+(ikS_yS_zS_x)p_x\left(ikS_x+S_zS_y\right)p_y\right\}\psi =0.$$
(39)   
For the case $`k=1`$, $`\psi \stackrel{}{\mathrm{\Psi }}`$, and using the representation (6) for the spin matrices, one obtains for Eq. (37)

$$\left(\begin{array}{ccc}p_x& p_y& p_z\\ 0& 0& 0\\ 0& 0& 0\end{array}\right)\stackrel{}{\mathrm{\Psi }}=0,$$
(40)   
for Eq. (38)

$$\left(\begin{array}{ccc}0& 0& 0\\ p_x& p_y& p_z\\ 0& 0& 0\end{array}\right)\stackrel{}{\mathrm{\Psi }}=0,$$
(41)   
and for Eq. (39)

$$\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& 0\\ p_x& p_y& p_z\end{array}\right)\stackrel{}{\mathrm{\Psi }}=0,$$
(42)   
which all are equivalent to Eq. (11). It is interesting to note that

$`\left(\begin{array}{ccc}p_x^2& p_xp_y& p_xp_z\\ p_yp_x& p_y^2& p_yp_z\\ p_zp_x& p_zp_y& p_z^2\end{array}\right)`$ $`=`$ $`\left(\begin{array}{ccc}p_x& 0& 0\\ p_y& 0& 0\\ p_z& 0& 0\end{array}\right)\left(\begin{array}{ccc}p_x& p_y& p_z\\ 0& 0& 0\\ 0& 0& 0\end{array}\right)`$ (52)   
$`=`$ $`\left(\begin{array}{ccc}0& p_x& 0\\ 0& p_y& 0\\ 0& p_z& 0\end{array}\right)\left(\begin{array}{ccc}0& 0& 0\\ p_x& p_y& p_z\\ 0& 0& 0\end{array}\right)`$ (59)   
$`=`$ $`\left(\begin{array}{ccc}0& 0& p_x\\ 0& 0& p_y\\ 0& 0& p_z\end{array}\right)\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& 0\\ p_x& p_y& p_z\end{array}\right),`$ (66)   
from which we deduce that the equations (10) and (11) and the decomposition (5), can be realized on the basis of Eq. (36) and one of the equations (37),(38) and (39).

Above, we have shown how all Maxwell equations can be derived simultaneously from first principles, similar to those which have been used to derive the Dirac relativistic electron equation. Moreover the wave function $`\stackrel{}{𝚿}`$ has a definite local classical interpretation in terms of the electric and magnetic fields, as given by Eq. (19), which are locally measurable and well understood quantities. Therefore Maxwell equations should be used as a guideline for proper interpretations of quantum theories.