I'll notate things in terms of clifford algebra and geometric calculus instead; it's a little cleaner.
Let $k = \omega e_0 + \kappa v$, where $v$ is a unit spatial vector.
In an EM plane wave, we take the magnitudes of the electric and magnetic fields to be the same, and we say that they're orthogonal to each other. In terms of the EM bivector $F$, this imposes a certain structure: Let $p$ be the direction of the electric field. Then we can write the EM field as
$$F = \alpha (e_0 p + vp) \exp(i k \cdot x)$$
where $e_0, p, v$ are all orthogonal, so these geometric products could be replaced with wedge products instead, if you prefer. The $e_0 p$ part corresponds to the electric field; the $pv$ part corresponds to the magnetic field.
Check: the magnitudes of the electric and magnetic fields should be equal. We can check this by taking $F^2$ and looking at the scalar part. Write $F = e_0 E + \epsilon_3 B$ and square it, looking at the scalar part:
$$\langle F^2 \rangle_0 = (e_0 E + \epsilon_3 B)^2 = e_0 E e_0 E + \epsilon_3 B \epsilon_3 B = E^2 - B^2$$
This is a well-known invariant of the EM tensor.
Looking at our posited form of the EM bivector, we get
$$\langle F^2 \rangle_0 = \alpha^2 \cos^2(k \cdot x) (e_0 + v) p (e_0 + v)p = -\alpha^2 p^2 (e_0 + v)^2 \cos^2(k \cdot x)$$
But $e_0$ and $v$ are unit vectors, so $e_0 + v$ is null, and $\langle F^2 \rangle_0 = 0$. This implies $E^2 = B^2$, which is exactly what we should get.
Now let's do some calculus. Maxwell's equations in vacuum reduce to $\nabla F = 0$, but $\nabla F = -ikF$ using our established form of the Faraday bivector. Let's write out those terms:
$$-i kF = -i \alpha [-\omega p + \kappa v e_0 p + \omega e_0 vp + \kappa p] \exp(i k \cdot x)$$
More clearly, let's separate out the vector terms:
$$\langle -ikF \rangle_1 = -i \alpha [-\omega p + \kappa p] \exp(i k \cdot x)$$
For this to be zero, $\omega = \kappa$ always. This means $k = \omega(e_0 + v)$, and we already established that vectors of the form $e_0 + v$ are null, so $k$ is null. The equation for the trivector terms gives the same result.
Incidentally, you might notice that $F = \ldots \alpha (e_0 + v)p \ldots$ also. That is, $F \propto kp \exp(i k \cdot x)$. This actually gives us a stronger result: That $F^2 = 0$, for $F^2 = kpkp = -pkkp = 0$. That also tells us about the pseudoscalar invariant of the Faraday tensor, which is proportional to $E \cdot B$, which we also constructed to be zero.