  Background: The interacting boson model (IBM) has been used extensively to
calculate the matrix elements governing neutrinoless double-beta decay. Studies
within other models indicate that a good description of neutron-proton pairing
is essential for accurate calculations of those matrix elements. The usual
interacting boson model is based only on like-particle pairs, however, and the
extent to which it captures neutron-proton pairing is not clear.
  Purpose: To determine whether neutron-proton pairing should be explicitly
included as neutron-proton bosons in IBM calculations of neutrinoless
double-beta decay matrix elements.
  Method: An isospin-invariant version of the nucleon-pair shell model is
applied to carry out shell-model calculations in a large space and in a
collective subspace, and to define effective operators in the latter. A
democratic mapping is then used to define corresponding boson operators for the
IBM, with and without an isoscalar neutron-proton pair boson.
  Results: IBM calculations with and without the isoscalar boson are carried
out for nuclei near the beginning of the $pf$ shell, with a realistic
shell-model Hamiltonian and neutrinoless double-beta-decay operator as the
starting point. Energy spectra and double-beta matrix elements are compared to
those obtained in the underlying shell model.
  Conclusions: The isoscalar boson does not improve energy spectra but does
improve double-beta matrix elements. To be useful at the level of precision we
need, the mapping procedure must be further developed to better determine the
dependence of the boson Hamiltonian and decay operator on particle number and
isospin. But the benefits provided by the isoscalar boson suggest that through
an appropriate combination of mappings and fitting, it would make IBM matrix
elements more accurate for the heavier nuclei used in experiments.
