Title: Constraints on the finite volume two-nucleon spectrum at $m_π\approx 806$ MeV

Abstract: The low-energy finite-volume spectrum of the two-nucleon system at a pion mass of $m_\pi \approx 806$ MeV is studied with lattice quantum chromodynamics (LQCD) using variational methods. The interpolating-operator sets used in [Phys.Rev.D 107 (2023) 9, 094508] are extended by including a complete basis of local hexaquark operators, as well as plane-wave dibaryon operators built from products of both positive- and negative-parity nucleon operators. Results are presented for the isosinglet and isotriplet two-nucleon channels. In both channels, these results provide compelling evidence for the presence of an additional state in the low-energy spectrum not present in a non-interacting two-nucleon system. Within the space of operators that are considered in this work, the strongest overlap of this state is onto hexaquark operators, and variational analyses based on operator sets containing only dibaryon operators are insensitive to this state over the temporal extent studied here. The consequences of these studies for the LQCD understanding of the two-nucleon spectrum are investigated.