Understanding the physics of the two-dimensional Hubbard model is widely believed to be a key step in achieving a full understanding of high-Tc cuprate superconductors. In recent years, progress has been made by large-scale numerical simulations at finite doping and, on the other hand, by microscopic theories able to capture the physics of individual charge carriers. In this work, we study single pairs of dopants in a cylindrical system using the density-matrix renormalization group algorithm. We identify two coexisting charge configurations that couple to the spin environment in different ways: a tightly bound configuration featuring (next-)nearest-neighbor pairs and a stripelike configuration of dopants on opposite sides of the cylinder, accompanied by a spin domain wall. Thus, we establish that the interplay between stripe order and uniform pairing, central to the models’ phases at finite doping, has its origin at the single-pair level. By interpolating between the Hubbard and the related t−J model, we are able to quantitatively understand discrepancies in the pairing properties of the two models through the three-site hopping term usually omitted from the t−J Hamiltonian. This term is closely related to a next-nearest-neighbor tunneling t′, which we observe to upset the balance between the competing stripe and pair states on the two-dopant level.
Blatz et al. (2025) studied this question.