Lattice Computation of the Nucleon Scalar Quark Contents at the Physical Point

We present a calculation of the $u$, $d$ and $s$ scalar quark contents of nucleons in QCD. Our calculation is based on $47$ lattice QCD ensembles with $N_f=2+1$ dynamical sea quarks, $5$ lattice spacings down to $0.054\,\text{fm}$, lattice sizes up to $6\,\text{fm}$ and pion masses down to $120\,\text{MeV}$. Using the Feynman-Hellmann theorem, we obtain $f_{ud}^N = 0.0405(40)(35)$ and $f_{s}^N=0.113(45)(40)$, which translates into $\sigma_{\pi N}=38(3)(3)\,\text{MeV}$, $\sigma_{sN}=105(41)(37)\,\text{MeV}$ and $y_N=0.20(8)(8)$ for the sigma terms and the related ratio, where the first errors are statistical and the second an estimate of the systematics. Using isospin relations and our recent result for the QCD contribution to the $n$-$p$ mass difference, we also compute the individual up and down quark contents of the proton and neutron. We obtain $f_u^p=0.0139(13)(12)$, $f_d^p=0.0253(28)(24)$, $f_u^n=0.0116(13)(11)$ and $f_d^n=0.0302(28)(25)$.