Consider a ralativistic particle of rest mass $m$ and electric charge $e $ moving in electromagnetic field with four-potential ${\displaystyle A^{\mu}=(\phi ,\mathrm {A} )} $ in vacuum, then the Hamilton–Jacobi equation has the form
$$g^{\mu \nu}\left ( \frac{\partial S}{\partial x^{\mu}} + \frac {e}{c}A_{\mu} \right ) \left ( \frac{\partial S}{\partial x^{\nu}} + \frac {e}{c}A_{\nu} \right ) = m^2 c^2\tag{1}$$
or more compact expressed as Minkowski product
$$ \left( \frac{\partial S}{\partial x^{\mu}} + \frac {e}{c}A_{\mu} \right ) \left ( \frac{\partial S}{\partial x_{\mu}} + \frac {e}{c}A^{\mu} \right ) = m^2 c^2 \tag{2}$$
here we denote $g^{\mu \nu}$ the metric tensor with signature $(+ - - -)$ and $S$ is the action function from Hamilton-Jacobi-theory.
Especially $S$ satisfy the equation
$$p_{\mu}= \nabla_{\mu}S := \frac{\partial S}{\partial x^{\mu}}\tag{3}$$
where $p_{\mu}$ is the four momentum and $\nabla_{\mu}$ the four gradient.
Now I have following two questions:
Does anybody have a reference for a rigorous derivation for $$ \left( \frac{\partial S}{\partial x^{\mu}} + \frac {e}{c}A_{\mu} \right ) \left ( \frac{\partial S}{\partial x_{\mu}} + \frac {e}{c}A^{\mu} \right ) = m^2 c^2 .\tag{4}$$
It is known that applying method of characteristics to the PDE
$$F(S,\frac{\partial S}{\partial x^{\mu}} ,x^{\mu}):= \left( \frac{\partial S}{\partial x^{\mu}} + \frac {e}{c}A_{\mu} \right ) \left ( \frac{\partial S}{\partial x_{\mu}} + \frac {e}{c}A^{\mu} \right ) - m^2 c^2 =0\tag{5}$$
one can derive the relative Lorentz equation
$${\displaystyle {\frac {\mathrm {d} p^{\mu }}{\mathrm {d} \tau }}=eF^{\mu \nu }p_{\nu }}\tag{6}$$
with electromagnetic tensor $$F^{\mu \nu }:= \frac{\partial A_{\mu}}{\partial x^{\nu}}- \frac{\partial A_{\nu}}{\partial x^{\mu}}\tag{7}$$ and four momentum $p_{\mu}$.
Here I'm also looking for an explicit derivation of LE from the HJE using characteristics.
Indeed, the method of characteristics transform a PDE into a system of ODE with respect parametrizing variable $\tau$:
$$\frac{\partial p_{\mu}}{\partial \tau}= -\frac{\partial F}{\partial x^{\mu}} -\frac{\partial F}{\partial S} p_{\mu}\tag{8}$$
$$\frac{\partial x_{\mu}}{\partial \tau}= \frac{\partial F}{\partial p^{\mu}}. \tag{9}$$
Remark: HJ theory says $$p_{\mu}= \frac{\partial S}{\partial x^{\mu}}.\tag{10}$$
The problem is to derive from here the equation for Lorentz force