Strong isospin violation and chiral logarithms in the baryon spectrum

TL;DR

This paper reports a precise lattice QCD calculation of the neutron-proton mass difference due to strong isospin breaking, providing evidence for non-analytic chiral behavior and updating the CP-odd pion-nucleon coupling estimate.

Contribution

First conclusive evidence of non-analytic light quark mass dependence in lattice QCD baryon spectrum calculations, with detailed analysis of isospin splittings and chiral logarithms.

Findings

Neutron-proton mass splitting from QCD: 2.32±0.17 MeV

Cascade baryon mass splitting: 5.44±0.31 MeV

Evidence for non-analytic chiral behavior in baryon spectrum

Abstract

We present a precise lattice QCD calculation of the contribution to the neutron-proton mass splitting arising from strong isospin breaking, $m_n-m_p|_{QCD}=2.32\pm0.17$ MeV. We also determine $m_{\Xi^-} - m_{\Xi^0}|_{QCD} = 5.44\pm0.31$ MeV. The calculation is performed at three values of the pion mass, with several values of the quark mass splitting and multiple lattice volumes, but only a single lattice spacing and an estimate of discretization errors. The calculations are performed on the anisotropic clover-Wilson ensembles generated by the Hadron Spectrum Collaboration. The omega-baryon mass is used to set the scale $a_t^{-1}=6111\pm127$ MeV, while the kaon masses are used to determine the value of the light-quark mass spitting. The nucleon mass splitting is then determined as a function of the pion mass. We observe, for the first time, conclusive evidence for non-analytic light quark mass dependence in lattice QCD calculations of the baryon spectrum. When left as a free parameter, the fits prefer a nucleon axial coupling of $g_A=1.24(56)$. To highlight the presence of this chiral logarithm in the nucleon mass splitting, we also compute the isospin splitting in the cascade-baryon system which is less sensitive to chiral dynamics. Finally, we update the best lattice QCD determination of the CP-odd pion-nucleon coupling that would arise from a non-zero QCD theta-term, $\bar{g}_0 / (\sqrt{2}f_\pi) = (14.7\pm1.8\pm1.4) \cdot 10^{-3} \bar{\theta}$. The original lattice QCD correlation functions, analysis results and extrapolated quantities are packaged in HDF5 files made publicly available including a simple Python script to access the numerical results, construct effective mass plots along with our analysis results, and perform the extrapolations of various quantities determined in this work.