I am trying to get the formula for energy of EM waves:
$$W = \frac{E^2 + B^2}{2}$$ calculating the work done on a test charge by the force: $$\mathbf F = q(\mathbf E + v \times \mathbf B)$$ $\mathbf E$ and $\mathbf B$ are vectors of the type $\mathbf F(u)$, $u = (\mathbf {k.x} - \omega t)$ and $\omega = \frac{c}{|k|}$ solutions of Maxwell wave equation. It seems go well until I get $$\frac{\partial E_v}{\partial t} = \mathbf {j.E}$$ where the left side is power per unit of volume and $\mathbf j$ is density of current.
But if I try to get rid of $\mathbf j$, using the Maxwell equation:
$$\mathbf j = \nabla \times \mathbf B - \frac{\partial \mathbf E}{\partial t}$$ the right side vanishes. And it is not a surprise, because the wave equation, from which $\mathbf E$ and $\mathbf B$ are solutions requires no charges or currents.
Searching the web, the energy formula comes from circuits, inductors and capacitors storing energy. Energy of EM waves simply use that results.
The other approach is from Lagrangian, but in this case, as I understand, it is the opposite way: the expression for the energy is postulated, and Maxwell equations are derived from it.
Is it possible to derive the quadratic energy expression from the wave equation and Lorentz force?