# Gauge invariance and Gauss's law in $U(1)$ E&M theory

I am reading an article: Stable Gapless Bose Liquid Phases without Any Symmetry (and also see Pretko's paper on the same subject). There, the authors wrote that Gauss's laws are associated with different gauge transformations as summarized in table 1 of Ref. 1. I am not quite sure how the authors obtained the Gauss's laws only from the gauge transformation of the vector potential. Take the simplest case rank-1 (the classical E&M) theory for example. The gauge transformation

$$$$A_i \to A_i + \partial_i \lambda,$$$$

where $$\lambda$$ is a scalar function of spacetime is accompanied by a corresponding change in the electrostatic potential $$\phi \to \phi - \partial_t \lambda$$ for the electric field to be invariant. How does this imply Gauss's law $$\partial_i E^i =\nabla \cdot \vec{E} =0$$? Are there some hidden assumptions? (e.g. for the classical pure E&M theory $$L = \int d^4 x E^2+B^2$$, it is known that in the Coulomb gauge, $$\nabla \cdot \vec{E} = 0$$) Could someone please help?

Take the variational derivative of the Maxwell action (the one you have written down) w.r.t. $$A_0$$ to get the Gauss law.
The relation to gauge transformations lies within the theory of constrained systems. Pass to the Hamiltonian formulation by applying the field-theoretic Legendre transform. In the temporal gauge $$A_0 = 0$$, your canonical coordinates will be the spatial components of the vector potential $$\vec{A}$$, and canonical momenta will be the components of the electric field $$\vec{E}$$. Gauge symmetry induces a first-class constraint on the canonical momenta, which are fixed by the Lagrange multiplier $$A_0$$ (which is why the variational derivative w.r.t. $$A_0$$ gives you the Gauss law). Also, the integrated Gauss law constraint generates gauge transformations on the phase space through Poisson brackets.
Another interesting and related observation is that in the naive quantum theory where the canonical momenta is replaced with (variational) derivatives $$E^i = - i \hbar \frac{\delta}{\delta A_i},$$
The quantum constraint equation $$\left( \nabla_i E^i \right) \Psi[\vec{A}] = 0$$
implies that $$\Psi[\vec{A}]$$ is a gauge invariant functional.