Gauge-boson mass matrix (${\mathcal{L}}_{\mathrm{mass}}=\frac{1}{2}{A}_{\mu }^a(M^2{)}_{ab}{A}^{b\mu }$)

With canonically normalized gauge field kinetic terms, a gauge-boson mass matrix is the real symmetric matrix of coefficients in the displayed quadratic mass term. Its eigenvalues are squared physical gauge boson masses. For a complex scalar field vacuum ${\Phi }_0$ in the Higgs mechanism, $(M^2{)}_{ab}=2\mathrm{Re}[(g_a{T}_a{\Phi }_0{)}^{†}(g_b{T}_b{\Phi }_0)]$. Thus for real $x_a$, $x^T{M}^{2}x=2\Vert {\sum }_a{x}_a{g}_a{T}_a{\Phi }_0{\Vert }^2\ge 0$: it is a Gram matrix and a positive semidefinite matrix. An unbroken generator annihilates the vacuum expectation value and lies in its null space. Diagonalizing the matrix identifies the massive and massless field combinations.