Varying the electromagnetic action \begin{equation} S=-m c \int d s\left(\dot{z}^{2}\right)^{\frac{1}{2}}-\frac{e}{c} \int d s A_{\mu} \dot{z}^{\mu}-\frac{1}{16 \pi c} \int d^{4} x F_{\mu \nu} F^{\mu \nu} \end{equation} we get the two equations of motion \begin{equation} m \ddot{z}^{\mu}=\frac{e}{c^{2}} F^{\mu \nu} \dot{z}_{\nu} \end{equation}
\begin{equation} \square A^{\mu}-\partial^{\mu}\left(\partial_{\nu} A^{\nu}\right)=\frac{4 \pi}{c} j^{\mu} \end{equation} If we use the Lorenz gauge, the solutions to the second equation are the Liénard–Wiechert potentials.
Is there a general solution to the second equation without fixing a gauge? By solution I mean a closed expression for $A^{\mu}$ in terms of $j^{\mu}$. Or is it possible to proof that such a solution does not exist.