So I understand that we select the divergence of A (magnetic vector potential) to be:
$$\frac{1}{c^2}\dot{\phi} + \nabla\cdot\vec{A} = 0.$$
The Lorenz gauge (1).
because of the symmetries in Maxwell's equations. Also, I know that the A and the scalar potential ($\varphi$) are not unique. We can write:
$$\phi \rightarrow \phi' = \phi-\dot{\Psi},$$
$$ \vec{A} \rightarrow \vec{A}' = \vec{A} +\nabla\Psi.$$
Which give the same E and B fields. Until here everything is general. So we have some freedom to select A and $\varphi$.
The question is: how does this imply the Lorenz gauge (Equation 1)? If there is further information needed to prove Equation 1 please let me know!