Formal Connection Between Symmetry and Gauss's Law In the standard undergraduate treatment of E&M, Gauss's Law is loosely stated as "the electric flux through a closed surface is proportional to the enclosed charge". Equivalently, in differential form, and in terms of the potential (in the static case):
$$\nabla^2 \phi = -\frac{\rho}{\epsilon_0}$$
Now, when using the integral form, one typically uses the symmetries of a known charge distribution to deduce related symmetries in the electric field, allowing the magnitude of the field to be factored out of the integral. To do so, one usually relies on intuitive, heuristic arguments about how the field in question "ought to" behave$^1$.
I'm wondering how one goes about formalizing this notion in precise mathematical terms. In particular, it seems that there ought to be an equivalent statement for Gauss's law in differential form, along the lines of "symmetries in $\rho$ induce related symmetries in $\phi$". Is there a way to formally state this claim? In particular:


*

*How would one formulate a proof of the conditions (necessary and sufficient) under which a symmetry of $\rho$ induces a symmetry in $\phi$?

*When it exists, how does one explicitly state the induced symmetry in terms of the known symmetry?

*Can such a result be generalized for arbitrary linear PDEs whose source terms exhibit some symmetry?


It seems to me like there must exist a concise, elegant, and general way to state and prove the above, but I can't quite seem to connect all the dots right now.

$^1$ See, for example, in Griffiths, Example 2.3, p. 72: "Suppose, say, that it points due east, at the 'equator.' But the orientation of the equator is perfectly arbitrary—nothing is spinning here, so there is no natural "north-south" axis—any argument purporting to show that $\mathbf{E}$ points east could just as well be used to show it points west, or north, or any other direction. The only unique direction on a sphere is radial."
 A: Let $\mathcal{D}$ be the operator corresponding to your equation ($-\epsilon_0 \nabla^2$ in this case). Let $U$ be some operator corresponding to the symmetry. It might be a rotation or parity transformation, etc.
If $U f = f$, we say the function $f$ is symmetric.
If $U\mathcal D=\mathcal D U$ as operators, then we say your equation is symmetric.
Let, $\mathcal {D} f= g$. If $\mathcal{D}$ is symmetric and $g$ is symmetric, then we can easily show $\mathcal{D}f = \mathcal{D}\mathcal U f$. If we can take an inverse of $\mathcal{D}$ we've proven $f$ is symmetric.
Taking an inverse of $\mathcal{D}$ is the the same thing as being able to solve the equation uniquely. In your particular case we can solve the equation uniquely if we restrict our function space to have some boundary conditions, say vanishing at infinity.
So this is the example you had in mind and this a formalization of the argument that $\phi$ must be symmetric.
Now if we can't solve the equation uniquely then there may be a loophole in the argument. A particular case I have in mind is a magnetic monopole which is rotationally symmetric, but the vector potential solution has a Dirac string and is not. But any two solutions $f$ and $Uf$ in this case are connected by a gauge transformation.
