Since, my similar question was closed, I will answer here.
The gauge transformation
$$ \mathbf {A} \rightarrow \mathbf {A} +\nabla \lambda, $$
$$ \varphi \rightarrow \varphi - \frac {\partial \lambda}{\partial t}, $$
(where $ \lambda=\lambda(\vec {r},t) $ is an arbitrary scalar function of coordinates $ \mathbf{r} $ and time $ t $) do not change the form of Maxwell's equations, and hence are admissible from a physical point of view.
In practice, no one chooses a special function $\lambda(\vec{r}, t)$ per se, although one is always implicitly assumed. But the described ambiguity of the potentials from mathematical point of view tells us that one can always be chosen to satisfy one arbitrary additional condition. One, since we can arbitrarily choose only one function $\lambda(\vec{r}, t)$.
For example, one can always choose the field potentials so that the scalar potential $\varphi = 0$ (is equal to zero). To make the vector potential equal to zero, is impossible, since the condition $\mathbf{A} = 0$ is three additional conditions (for the three components of $\mathbf{A}$).
Another possible way is to choose one arbitrary additional is
- Coulomb gauge: $\mathrm{div}\mathbf{A}' = 0 $. Where an arbitrary function can be chosen so that it satisfies the condition $\nabla^2 \lambda= -\nabla \cdot \mathbf{A}$. By solving the equation, one can get the following function.
In practice, they proceed as follows: imagine that we know that a given $\mathbf{A}$ solves the Maxwell's Equation. We can always find a gauge transformation that converts it to a new solution $\mathbf{A}'$ that satisfies the Lorentz Gauge condition $$\mathrm{div}\mathbf{A}' = 0. $$ Since we can always do this, rather than imposing this condition after solving Maxwell's Equations, we can require that we are looking for solutions of this type $$\mathrm{div}\mathbf{A} = 0. $$ before we solve that equation.
- For Lorenz gauge $\mathrm{div}\,\mathbf{A}' + {1 \over c^2}{\partial \mathbf{\varphi}' \over \partial t} = 0$, an arbitrary function can be chosen so that it satisfies the condition $\Box^2 \lambda= - \nabla \mathbf{A} - \frac{1}{c^2}\frac{\partial \varphi}{\partial t}$.
Just as a clever choice of coordinates can make solving a problem easier, so
we will find that a clever choice of gauge can make finding solutions easier .