I am studying QFT with Mandl & Shaw, Quantum Field Theory and ran into Problem 2.3 (page 37 in the Second Edition) that says:
Problem 2.3: Show that the Lagrangian density
$$\mathscr{L} = -\frac{1}{2}\big[\partial_\alpha\phi_\beta(x)\big]\big[\partial^\alpha\phi^\beta(x)\big] + \frac{1}{2}\big[\partial_\alpha\phi^\alpha(x)\big]\big[\partial_\beta\phi^\beta(x)\big] + \frac{\mu^2}{2}\phi_\alpha(x)\phi^\alpha(x)$$
for the real vector field $\phi^\alpha(x)$ leads to the field equations
$$\left[g_{\alpha\beta}\left(\Box+\mu^2\right) - \partial_\alpha\partial_\beta\right] \phi^\beta(x) = 0$$
and that the field $\phi^\alpha(x)$ satisfies the Lorenz condition
$$\partial_\alpha \phi^\alpha(x) = 0.$$
I have done the first part, but am not sure how to approach the second part---deriving the Lorenz gauge condition. I now know how to derive both the results.
However, why/how is this possible? The Lagrangian looks very much like a Klein-Gordon Lagrangian for each field component, except for the second-term. However, if the Lorenz condition was indeed derivable from the given Lagrangian, the second term will be identically zero, and the Lagrangian will reduce completely to a Klein-Gordon Lagrangian for each component. Then, what was the point of the second term in the Lagrangian from the beginning? Its contribution to the equation of motion will also be identically zero, if Lorenz condition held.
Is there a meaningful physical interpretation to be learned here?