It is known that the Coulomb potential can be obtained by Fourier transform of the propagator from E&M. Is this because one of Maxwell's equations have the form $\nabla \cdot \mathbf{E}=\rho$?
1 Answers
To obtain the classical (in the sense of non quantum mechanical) potential due to a pair of charged particles, says, two electrons, you may calculate the scattering amplitude and compare it to
$$\langle p'|iT|p \rangle=-i2\pi\tilde{V}( \textbf{q} )\delta^{3}(E_{p'}-E_{p}), \hspace{2cm} \bf{q}=\bf{p'}-\bf{p'} $$
this is equation $(4.123)$ in Peskin's & Schroeder - An Introduction to Quantum Field Theory, and you in the first chapter of this lecture notes (link) that expression is justified. In the book this is first done to the Yukawa potential (scalar-fermion theory) and then generalised to Quantum Electrodynamics (vector-fermion theory).
The form of the Hamiltonian interaction for QED, which is the relevant piece of information to do the scattering amplitude computation (in particular, the lowest order contribution) has nothing to do with the classical form of Maxwell's equations but with the imposition of the $U(1)$ symmetry of the theory. Pretty much that underlying symmetry dictates the form of the quantum amplitudes and classical field theory equations.
In particular, Gauss law is pretty much the same thing (contains the same information) as the Coulomb equation, the latter can be derived from the former and vice versa.
You may find the following question interesting Alternative methods to derive the static potential in the NR limit of QED