Nuclear matrix elements of neutrinoless double-beta decay in covariant density functional theory with different mechanisms

Nuclear matrix elements (NMEs) for neutrinoless double-beta (0νββ) decay in candidate nuclei play a crucial role in interpreting results from current experiments and in designing future ones. Accurate NME values serve as important nuclear inputs for constraining parameters in new physics, such as neutrino mass and the Wilson coefficients of lepton-number-violating (LNV) operators. In this study, we present a comprehensive calculation of NMEs for 0νββ decay in ^76Ge, ^82Se, ^100Mo, ^130Te, and ^136Xe, using nuclear wave functions obtained from multi-reference covariant density functional theory (MR-CDFT). We employ three types of transition potentials at the leading order in chiral effective field theory. Our results, along with recent data, are utilized to constrain the coefficients of LNV operators. We find that NMEs for the standard and short-range mechanisms are significantly larger than those for the non-standard long-range mechanism. The use of NMEs from various nuclear models does not notably change the parameter space intervals for the coefficients, although MR-CDFT yields the most stringent constraint. Furthermore, our NMEs can also be used to perform a more comprehensive analysis with multiple isotopes.