Competing inhomogeneous orders are a central feature of correlated electron materials, including the high-temperature superconductors. The two-dimensional Hubbard model serves as the canonical microscopic physical model for such systems. Multiple orders have been proposed in the underdoped part of the phase diagram, which corresponds to a regime of maximum numerical difficulty. By combining the latest numerical methods in exhaustive simulations, we uncover the ordering in the underdoped ground state. We find a stripe order that has a highly compressible wavelength on an energy scale of a few kelvin, with wavelength fluctuations coupled to pairing order. The favored filled stripe order is different from that seen in real materials. Our results demonstrate the power of modern numerical methods to solve microscopic models, even in challenging settings.
The Hubbard model is an iconic model in quantum many-body physics and has been intensely studied, especially since the discovery of high-temperature cuprate superconductors. Combining the complementary capabilities of two computational methods, we found superconductivity in both the electron- and hole-doped regimes of the two-dimensional Hubbard model with next-nearest-neighbor hopping. In the electron-doped regime, superconductivity was weaker and was accompanied by antiferromagnetic Néel correlations at low doping. The strong superconductivity on the hole-doped side coexisted with stripe order, which persisted into the overdoped region with weaker hole-density modulation. These stripe orders varied in fillings between 0.6 and 0.8. Our results suggest the applicability of the Hubbard model with next-nearest hopping for describing cuprate high–transition temperature ( T c ) superconductivity.
Stripe order is a prominent feature in the phase diagram of the cuprate superconductors and has been established as the lowest-energy state of the two-dimensional Fermi-Hubbard model in relevant regimes. As temperature rises, stripes and superconductivity give way to the interesting strange-metal and pseudogap regimes. Here we investigate a crucial aspect of these regimes using numerical simulations of the square lattice Hubbard model. Using infinite projected entangled pair states with purification in the thermodynamic limit, we show that the thermodynamic charge susceptibility develops a broad maximum above the stripe regime at filling near n = 0.90, that strengthens upon cooling. Minimally entangled typical thermal states on finite cylinders attribute this enhancement to the emergence of large, fluctuating charge clusters, while the charge susceptibility remains finite. Upon further cooling, the cluster sizes stop broadening and lock to the stripe wavelength, indicating that interaction-driven clustering organizes the intermediate-temperature landscape but is ultimately forestalled by stripe order at low temperature.
Understanding the physics of the two-dimensional Hubbard model is widely believed to be a key step in achieving a full understanding of high- T c 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.
Abstract Several state-of-the-art numerical methods have observed static or fluctuating spin and charge stripes in doped two-dimensional Hubbard models, suggesting that these orders play a significant role in shaping the cuprate phase diagram. Many experiments, however, also indicate that the cuprates have strong electron-phonon (e-ph) coupling, and it is unclear how this interaction influences stripe correlations. We study static and fluctuating stripe orders in the doped single-band Hubbard-Holstein model using zero temperature variational Monte Carlo and finite temperature determinant quantum Monte Carlo. We find that the lattice couples more strongly with the charge component of the stripes, leading to an enhancement or suppression of stripe correlations, depending on model parameters like the next-nearest-neighbor hopping
t′ t′ or phonon energy Ω. Our results help elucidate how thee-ph interaction can tip the delicate balance between stripe and superconducting correlations in the Hubbard-Holstein model with implications for our understanding of the high-Tccuprates.
We determine the spin and charge orders in the ground state of the doped two-dimensional (2D) Hubbard model in its simplest form, namely with only nearest-neighbor hopping and on-site repulsion. At half-filling, the ground state is known to be an antiferromagnetic Mott insulator. Doping Mott insulators is believed to be relevant to the superconductivity observed in cuprates. A variety of candidates have been proposed for the ground state of the doped 2D Hubbard model. A recent work employing a combination of several state-of-the-art numerical many-body methods established the stripe order as the ground state near
1/8 doping at strong interactions. In this paper, we apply one of these methods, the cutting-edge constrained-path auxiliary field quantum Monte Carlo (AFQMC) method with self-consistently optimized gauge constraints, to systematically study the model as a function of doping and interaction strength. With careful finite size scaling based on large-scale computations, we map out the ground state phase diagram in terms of its spin and charge order. We find that modulated antiferromagnetic order persists from near half-filling to about
1/5 doping. At lower interaction strengths or larger doping, these ordered states are best described as spin-density waves, with essentially delocalized holes and modest oscillations in charge correlations. When the charge correlations are stronger (large interaction or small doping), they are best described as stripe states, with the holes more localized near the node in the antiferromagnetic spin order. In both cases, we find that the wavelength in the charge correlations is consistent with so-called filled stripes.
The long-studied Hubbard model is one of the simplest models of copper-oxide superconductors. However, the connection between the model and the experimental phase diagram is still under debate, in particular regarding the existence and extent of the
d-wave superconducting phase. Recent rapid progress in improving the accuracy of numerical solvers has opened a way to answer this question reliably. Here, we study the hole-doping concentration ($\ensuremath{\delta}$) dependence of the Hubbard model in the ground states on a square lattice at strong coupling
U/t=10, for the on-site interaction
U and the transfer
t, using a variational Monte Carlo method. The method, which combines tensor network and Lanczos methods on top of Pfaffian wave functions, reveals a rich phase diagram, in which many orders compete severely and degenerate within the energy range of
0.01t. We have identified distinct phases including a uniform
d-wave superconducting phase for $0.17\ensuremath{\lesssim}\ensuremath{\delta}\ensuremath{\lesssim}0.22$ and a stripe charge/spin ordered phase for $\ensuremath{\delta}\ensuremath{\lesssim}0.17$ with the stripe period depending on $\ensuremath{\delta}$, together with presumable spatially coexisting antiferromagnetic and stripe order for $\ensuremath{\delta}\ensuremath{\lesssim}0.07$ and coexisting stripe and
d-wave superconductivity for $0.07\ensuremath{\lesssim}\ensuremath{\delta}\ensuremath{\lesssim}0.17$. The present, improved method revealed a wider region of a charge uniform superconducting phase than the previous studies and shows a qualitative similarity to the phase diagram of the cuprate superconductors. The superconducting order parameter is largest at doping of around $\ensuremath{\delta}=0.17$ in the ground state, which undergoes phase transitions from an inhomogeneous to a uniform state.