Maxwell's Equations with ##\small \rho = q_e\psi^{*}\psi##. Here, ##\small q_e## is the electron charge. So long as there is no transition, magnetic field remains constant and we have ##\small \nabla\cdot E=\rho/\epsilon_0## and ##\small \nabla\times E=0##. The reason all of this works is that E is a linear field, so an electric field from superposition of sources is a superposition of fields from these sources, and that's just a sum of these fields. So probability density of electron is the charge density of electron. If you have entanglement, things can get interesting, but otherwise, it's extremely straight forward. So let's look at the hydrogen atom in ground state.
[tex]\psi_{100}(r) = \frac{1}{\sqrt{\pi a_0^3}}e^{-r/a_0}[/tex]
Here, ##\small a_0## is the Bohr Radius, of course. This gives us a perfectly spherically symmetric distribution, so applying Gauss law, the equation for electric field follows.
[tex]E(r) = \frac{1}{4\pi \epsilon_0}\frac{q_e\hat{r}}{r^2}\int_0^r \psi_{100}^{*}(r)\psi_{100}(r) dr = \frac{1}{4\pi \epsilon_0}\frac{e\hat{r}}{r^2}\left( 1-\frac{e^{-2r/a_0}(a_0^2+2ra_0+2r^2)}{a^2} \right)[/tex]
Of course, that's just the electron charge. If we say that nucleus has charge ##\small q_p = -q_e##, and we assume it to be a point source, which it really close to being, then we have the net electric field in the neighborhood of a hydrogen atom.
[tex]E(r) = \frac{q_p\hat{r}}{4\pi \epsilon_0}\frac{a_0^2+2ra_0+2r^2}{a_0^2r^2}e^{-2r/a_0}[/tex]
Note that because of the exponent, this is basically zero for ##\small r>>a_0##.