Tensor-Network Analysis of Root Patterns in the XXX Model with Open Boundaries

The string hypothesis of Bethe roots is a cornerstone in the thermodynamic analysis of quantum integrable systems, since it connects root configurations with physical quantities such as the ground-state energy, surface energy and excitation spectra. For integrable models with $U(1)$ symmetry, this connection is well established. When the $U(1)$ symmetry is broken by generic non-diagonal boundary fields, however, the off-diagonal Bethe Ansatz leads to an inhomogeneous $T\text{--}Q$ relation whose Bethe roots have highly nontrivial distributions. This raises two fundamental questions: whether the zero roots and the ODBA Bethe roots still possess regular and classifiable structures in the large-size limit, and whether such structures can be used to extract physical quantities. In this work, we address these two questions for the isotropic Heisenberg spin chain with non-diagonal open boundaries. By combining tensor-network algorithms with Bethe-Ansatz techniques, we determine the zero-root and Bethe-root configurations associated with the $\Lambda\text{--}\theta$ relation and the inhomogeneous Bethe Ansatz equations for large system sizes, up to $N\simeq 60$ and $100$. We find that, despite the absence of $U(1)$ symmetry, the roots exhibit well-organized patterns. The zero roots form bulk strings, boundary strings and additional roots, while the ODBA Bethe roots split into four geometric classes: regular roots, line roots, arc roots and paired-line roots.