Higgs boson coupling to Z bosons (${\mathcal{L}}_{HZ}=m_Z^{2}H{Z}_{\mu }{Z}^{\mu }/v+m_Z^2{H}^2{Z}_{\mu }{Z}^{\mu }/(2v^2)$)

The neutral vector mass term is $\frac{1}{2}{m}_Z^2{Z}_{\mu }{Z}^{\mu }(1+H/v{)}^2$. Expanding it gives the cubic $HZZ$ and quartic $HHZZ$ interactions in unitary gauge. Differentiating the action with respect to identical fields gives Feynman rules $2i{m}_Z^2{g}_{\mu \nu }/v$ and $2i{m}_Z^2{g}_{\mu \nu }/v^2$. A single $Z$ does not couple to two copies of the real neutral radial Higgs boson through an elementary vertex.

For $m_H>2m_Z$, the Higgs boson coupling to Z bosons gives $\mathcal{M}=(2m_Z^2/v){ϵ}_1^{\ast }\cdot {ϵ}_2^{\ast }$. The massive vector polarization sum gives $\sum |\mathcal{M}{|}^2=(m_H^4/v^2)(1-4x+12x^2)$, with $x=m_Z^2/m_H^2$. The two-body decay phase space and identical final-state symmetry factor $1/2$ giveThe contraction uses $k_1\cdot k_2=(m_H^2-2m_Z^2)/2$. The on-shell rate vanishes at threshold; below threshold an off-shell calculation is required.