Heavy quarkq¯qmatrix elements in the nucleon from perturbative QCD

The scalar heavy quark content of the nucleon, ${\mathcal{M}}_q=⟨N|m_q\overline{q}q|N⟩,$ is relevant for computing the interaction of dark matter candidates with ordinary matter, while ${\mathcal{M}}_s$ is important for predicting the properties of dense matter. We compute ${\mathcal{M}}_q$ in perturbative QCD to $\mathcal{O}({\alpha }^3).$ As one goes from ${\mathcal{M}}_t$ to ${\mathcal{M}}_c$ the leading order contribution decreases as the number of light quarks is dropping, while the radiative corrections grow and are all positive. The leading source of uncertainty in the calculation is due to the poorly known value of ${\mathcal{M}}_s.$ The Higgs-nucleon coupling is estimated to be $g_{hNN}\simeq 1.7\times {10}^{-3}.$