Lattice QCD determination of the radiative decay rates $h_{c}\to \eta_{c}\, \gamma$ and $h_{b}\to \eta_{b}\, \gamma$

We present the results of our lattice QCD computation of the hadronic matrix elements relevant to the $h_{c}\to \eta_{c}\gamma$ and $h_{b}\to \eta_{b}\gamma$ decays by using the gauge configurations produced by the Extended Twisted Mass Collaboration with $N_{f}=2+1+1$ dynamical Wilson-Clover twisted mass fermions at five different lattice spacings with physical dynamical $u$ , $d$, $s$ and $c$ quark masses (except for the the coarsest lattice for which the lightest sea quark corresponds to a pion with $m_{\pi}\simeq 175~\mathrm{MeV}$). While the hadronic matrix element for $h_{c}\to \eta_{c}\gamma$ is obtained directly, the one relevant to $h_{b}\to\eta_{b}\gamma$ is reached by working with heavy quark masses $m^{(n)}_{H} = \lambda^{n-1} m_{c}$, with $\lambda \sim 1.24$ and $n=1,2, \ldots ,6$, and then extrapolated to $m_{b}$ by several judicious ans\"atze. In the continuum limit we obtain $\Gamma( h_{c}\to \eta_{c} \gamma ) = 0.604(24)~\mathrm{MeV}$, which is by a factor of $2.3$ more accurate than the previous lattice estimates, and in good agreement with the experimental measurement. In the $b$-quark case we obtain $\Gamma( h_{b}\to \eta_{b} \gamma) =46.0(4.8)~\mathrm{keV}$.