Abstract We investigate the hole-doped Hubbard model on a honeycomb lattice using a fermionic projected entangled pair states (fPEPS) method. Our study reveals the presence of quasi-long-range order of Cooper pairs, characterized by power-law decay of correlation functions with exponents K > 1. We further analyze the competing phases of superconductivity, specifically the antiferromagnetic (AFM) order and the charge density wave (CDW) order. Our results show that there are domain wall structures when the hole doping δ is small and the Coulomb parameter U is large. However, these structures disappear as we increase the hole doping δ or decrease U. Furthermore, for small hole doping, the system exhibits AFM order, which diminishes for δ > 0.05. Conversely, as the doping level increases, the CDW order gradually decreases. Notably, a considerable CDW order persists even at higher doping levels. These findings suggest a progressive suppression of the AFM order and a growing prominence of the CDW order with increasing δ.
We have performed large-scale density-matrix renormalization group studies of the lightly doped Hubbard model on the honeycomb lattice on long three- and four-leg cylinders. We find that the ground state of the system upon lightly doping is consistent with that of a superconducting state with coexisting quasi-long-range superconducting and charge density wave orders. Both the superconducting and charge density wave correlations decay as a power law at long distances with corresponding exponents
Ksc<2 and
Kc<2. On the contrary, the spin-spin and single-particle correlations decay exponentially, although with relatively long correlation lengths.
The square-lattice Hubbard and closely related t-J models are considered as basic paradigms for understanding strong correlation effects and unconventional superconductivity (SC). Recent large-scale density matrix renormalization group simulations on the extended t-J model have identified d-wave SC on the electron-doped side (with the next-nearest-neighbor hopping t_{2}>0) but a dominant charge density wave (CDW) order on the hole-doped side (t_{2}<0), which is inconsistent with the SC of hole-doped cuprate compounds. We re-examine the ground-state phase diagram of the extended t-J model by employing the state-of-the-art density matrix renormalization group calculations with much enhanced bond dimensions, allowing more accurate determination of the ground state. On six-leg cylinders, while different CDW phases are identified on the hole-doped side for the doping range δ=1/16-1/8, a SC phase emerges at a lower doping regime, with algebraically decaying pairing correlations and d-wave symmetry. On the wider eight-leg systems, the d-wave SC also emerges on the hole-doped side at the optimal 1/8 doping, demonstrating the winning of SC over CDW by increasing the system width. Our results not only suggest a new path to SC in general t-J model through weakening the competing charge orders, but also provide a unified understanding on the SC of both hole- and electron-doped cuprate superconductors.
Recently, a flurry of works have found strong competition between charge density wave (CDW) and superconductivity (SC) in the doped Hubbard and $t\text{\ensuremath{-}}J$ models on the square lattice. Interestingly, some recent results suggest that the electron-phonon coupling may suppress CDW order and enhance SC. In this work, we consider the square-lattice Hubbard model with the Holstein or Su-Schrieffer-Heeger electron-phonon coupling at the large-
U and antiadiabatic (infinite phonon frequency) limits, which gives an effective $t\text{\ensuremath{-}}J$ model with either a density attractive interaction
V or a
JP term that contributes a larger spin exchange and a density repulsive interaction. To explore how these effective couplings may tune the competition between CDW and SC, we implement the density matrix renormalization group simulations on the $t\text{\ensuremath{-}}{t}^{\ensuremath{'}}\text{\ensuremath{-}}J$ model with
V or
JP coupling. We focus on the six-leg cylinder system with the next-nearest-neighbor hopping ${t}^{\ensuremath{'}}$, which hosts partially filled stripe and
d-wave SC in phase diagram. By tuning ${t}^{\ensuremath{'}}/t>0$ and
V or
JP, we establish two quantum phase diagrams at
121 doping ratio. In the SC phases, the increased
V or
JP coupling can enhance the quasi-long-range SC order, consistent with some previous findings. Nonetheless, no SC emerges when the partially filled stripe phase disappears with increased
V or
JP. Instead, the system has a transition to either a phase-separation-like regime or a filled stripe phase. On the other hand, with increased ${t}^{\ensuremath{'}}/t$, not only the partially filled stripe but the phase separation and filled stripe can also be tuned to SC phase. Our results suggest that although
V and
JP couplings may strengthen hole binding, the hole dynamics controlled by ${t}^{\ensuremath{'}}/t$ appears to play a more crucial role for obtaining a SC in the $t\text{\ensuremath{-}}J$ model.
Inspired by the recent experimental advances in cold atom quantum simulators, we explore the experimentally implemented bosonic t-t^{'}-J model on the square lattice using large-scale density matrix renormalization group simulations. By tuning the doping level δ and hopping ratio t^{'}/t, we uncover six distinct quantum phases, several of which go far beyond the conventional paradigm of phase-coherent superfluidity (SF) expected for bosonic systems. In particular, in the presence of antiferromagnetic (AFM) order, doped holes are tightly bound into pairs, giving rise to a pair density wave phase at low doping and small |t^{'}/t|, which is suppressed on the t^{'}<0 side, resulting in a disordered pair density wave state that lacks coherence of either individual bosons or pairs. Upon further doping, bosons can regain phase coherence and form an SF* state, characterized by condensation at emergent incommensurate momenta concurrent with an incommensurate magnetic order. On the t^{'}>0 side, the sign-induced kinetic frustration inherently disfavors local AFM correlations, leading to a phase separation in which doped holes cluster into ferromagnetic (FM) domains spatially separated by undoped AFM regions. Upon further doping, this inhomogeneous state evolves into a uniform SF+xy-FM phase. Finally, we propose a concrete experimental scheme to realize both signs of t^{'}/t in Rydberg tweezer arrays, with an explicit mapping between model parameters and experimentally accessible regimes. Our results reveal competing and intertwined orders in doped antiferromagnets, which are relevant to central issues in high-T_{c} superconductivity, reflecting the frustrated interplay between doped holes and spin background.
In view of quite recent experimental activities on magnetic and superconducting properties of honeycomb and hexagonal lattice based materials, in this thesis we have used field-theoretic and many-body methods to investigate magnetic and superconducting properties of the large-U Hubbard model on the honeycomb lattice at half-filling and in the hole-doped regime. Within the framework of a functional-integral approach, we obtain the Lagrangian density associated with the charge (Grassmann fields) and spin [SU(2) gauge fields] degrees of freedom. The Hamiltonian related to the charge degrees of freedom is exactly diagonalized. In the strong-coupling regime, we derive a perturbative low-energy theory suitable to describe the (quantum) magnetic and superconducting phases at half-filling and in the hole-doped regime. At half-filling, we deal with the underlying spin degrees of freedom of the quantum antiferromagnetic (AF) Heisenberg model by employing a second-order spin-wave analysis, in which case we have calculated the ground-state energy and the staggered magnetization; the results are in very good agreement with previous studies. Further, in the continuum, we derive a nonlinear σ-model with a topological Hopf term that describes the AF-VBS (valence bond solid) competition. In the challenging hole-doped regime, our approach allows the derivation of a t-J Hamiltonian, and the analysis of the role played by charge and spin quantum fluctuations on the ground-state energy and, particularly, on the breakdown of the AF order at a critical hole doping; the results are benchmarked against recent Grassmann tensor product state simulations. In addition, we have performed an extensive study of the electronic structure of the doped system for each competing phase: AF, ferromagnetic (FM), and (spin-singlet pairing) s-, dx₂₋ʏ₂ – and idxʏ -wave superconducting (SC) state induced by purely electronic effects. In this context, an energetic analysis of the ground state of these phases reveal that the AF order prevails for low hole doping, while a dominantly chiral dx₂₋ʏ₂ + idxʏ superconducting state was found in the vicinity of the Van Hove singularity (high hole doping). We also stress that a thermodynamic analysis of the superconducting phase shows that the critical temperature is directly related to the exchange constant J = 4t²/U, in which t denotes the hopping amplitude and U the on-site Coulomb repulsion of the Hubbard model (purely electronic origin). Remarkably, the competition between the AF and dx₂₋ʏ₂ + idxʏ SC phases takes place by the occurrence of a first-order transition accompanied by a spatial phase separation of the referred phases.
We study the honeycomb lattice $t\text{\ensuremath{-}}J$ model using the fermionic tensor network approach. By examining the ansatz with various unit cells, we discover several different stripe states with different periods that compete strongly with uniform states. At very small doping $\ensuremath{\delta}<0.05$, we find almost degenerate uniform
d-wave superconducting ground states coexisting with antiferromagnetic order. While at larger doping $\ensuremath{\delta}>0.05$, the ground state is an approximately half-filled stripe-ordered state, where the stripe period decreases with increasing hole doping $\ensuremath{\delta}$. Furthermore, the stripe states with the lowest variational energy always display ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave pairing symmetry. The similarity between our results and those on the square lattice contributes to a more comprehensive understanding of doped Mott insulators.