The Hubbard model has attracted considerable interest due to its prototypical role in describing strongly interacting electronic systems, such as high-critical-temperature superconductors as well as many novel quantum materials. By introducing next-nearest-neighbor (NNN) hoppings to the Hubbard model, the phase diagram becomes richer, and fascinating phenomena arise in both, one-dimensional chains and square lattices, such as: antiferromagnetism, ferromagnetism, superconductivity (SC), as well as charge orders, among others. Moreover, NNN hoppings play a fundamental role in understanding effects of doping on magnetism and pairing orders in strongly interacting regimes. In this article, we review the recent progress in understanding the different competing phases of this model in one and two dimensions from a computational perspective. We comment on the pressing technical challenges, illustrate the controversial results concerning the emergence of the SC phase, and conclude with our perspectives on future explorations.
The Hubbard model has attracted considerable interest due to its prototypical role in describing strongly interacting electronic systems, such as high-critical-temperature superconductors as well as many novel quantum materials. By introducing next-nearest-neighbor (NNN) hoppings to the Hubbard model, the phase diagram becomes richer, and fascinating phenomena arise in both, one-dimensional chains and square lattices, such as: antiferromagnetism, ferromagnetism, superconductivity (SC), as well as charge orders, among others. Moreover, NNN hoppings play a fundamental role in understanding effects of doping on magnetism and pairing orders in strongly interacting regimes. In this article, we review the recent progress in understanding the different competing phases of this model in one and two dimensions from a computational perspective. We comment on the pressing technical challenges, illustrate the controversial results concerning the emergence of the SC phase, and conclude with our perspectives on future explorations.
The emergence of superconductivity (SC) in lattice models, such as the attractive Hubbard one, has renewed interest since the realization of cold-atom experiments. However, reducing the temperature in these experiments is a bottleneck; therefore, investigating how to increase the energy scale for SC is crucial to cold atoms. In view of this, we examine the effects of next-nearest-neighbor hoppings (t′) on the pairing properties of the attractive Hubbard model in a square lattice. To this end, we analyze the model through unbiased Quantum Monte Carlo simulations for fixed density n=0.87, and perform finite-size scaling analysis to the thermodynamic limit. As our main result, we notice that the existence of further hopping channels leads to an enhancement of the pairing correlations, which, in turn, increases the ground-state order parameter Δ. Finally, at finite temperatures, for t′/t≠0, this enhancement of pairing correlations leads to an increase in the critical temperature Tc. That is, the fine-tuning of second-neighbor hoppings increases the energy scales for SC, and may be a route by which cold-atom experiments can achieve such a phase and to help us further understand the nature of this phenomenon.
Using the finite-temperature determinant quantum Monte Carlo (DQMC) algorithm, we study the pairing symmetries of the Hubbard Hamiltonian with next-nearest-neighbor (NNN) hopping t' on square lattices. By varying the value of t', we find that the d-wave pairing is suppressed by the onset of t', while the p + ip-wave pairing tends to emerge for low electron density and t' around -0.7. Together with the calculation of the anti-ferromagnetic and ferromagnetic spin correlation function, we explore the relationship between anti-ferromagnetic order and the d-wave pairing symmetry, and the relationship between ferromagnetic order and the p + ip-wave pairing symmetry. Our results may be useful for the exploration of the mechanism of the electron pairing symmetries, and for the realization of the exotic p + ip-wave superconductivity.
Correlated electron systems with competing interactions provide a valuable platform for examining exotic magnetic phases. Theoretical models often focus on nearest-neighbor interactions, although long-range interactions can have a considerable impact on the behavior of the system, creating distinct and complicated magnetic phases. We investigate the consequences of competing interactions in a half-filled one-band Hubbard model on a simple cubic lattice, incorporating hopping processes up to third-nearest-neighbors, to explore the underlying magnetotransport properties. Our magnetic phase diagrams at low temperatures, obtained using semiclassical Monte Carlo analysis, reveal that the long-range interactions can disrupt one form of magnetic phase while creating a new type of magnetic order. For the nonperturbative regime (on-site Hubbard repulsive strength $U\ensuremath{\sim}$ bandwidth), the C-type antiferromagnetic ground state is preferred over the G-type antiferromagnetic phase when the interaction between second-nearest-neighbor sites becomes significant relative to the nearest-neighbor interactions. However, interactions beyond the second-nearest-neighbors are required to stabilize the A-type antiferromagnetic ground state. Remarkably, at low temperatures, a highly correlated paramagnetic insulating phase develops at the intersection between the antiferromagnetic phases, which might promote a three-dimensional spin-liquid state.
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.
The attractive Hubbard model plays a paradigmatic role in the study of superconductivity (superfluidity) and has become directly realizable in ultracold atom experiments on optical lattices. However, the critical temperatures, Tc's, remain lower than the lowest temperatures currently achievable in experiments. Here, we explore a possible route to enhance Tc by introducing an additional next-nearest-neighbor (NNN) hopping, t′, in a two-dimensional square lattice. We perform sign-problem-free determinant quantum Monte Carlo simulations to compute response functions such as pairing correlation functions, superfluid density, and uniform spin susceptibility. Our results show that a judicious choice of t′ can increase Tc by up to 50% compared to the case with only nearest-neighbor hopping. In contrast, the preformed pairs temperature scale, named pairing temperature, Tp, decreases with increasing ∣t′/t∣, which should represent a reduction of the pseudogap region, favoring a more BCS-like behavior at intermediate coupling. We further analyze the interacting density of states to characterize the transition from a pseudogap regime to a fully gapped superconducting state. These findings suggest that NNN hopping could be a viable route to increase Tc to values closer to experimentally accessible temperature scales.
Abstract Using the strong coupling diagram technique, we investigate the extended Hubbard model on a two-dimensional square lattice. This approach allows for charge and spin fluctuations and a short-range antiferromagnetic order at nonzero temperatures. The model features the first-order phase transition to states with alternating site occupations (SAO) at the intersite repulsion v = v c . In this work, we show that doping decreases v c . For a nonzero next-nearest neighbor hopping t ′ , less mobile carriers produce a stronger fall in v c . For half-filling and t ′ = 0 , we consider phase diagrams for a fixed temperature T, on-site U, and intersite repulsions. The diagrams contain regions of SAO, Mott insulator, and several metallic states distinguished by their densities of states. Two of them are characterized by a dip and a peak at the Fermi level. The dip originates from the Slater mechanism for itinerant electrons, while the peak to bound states of electrons with localized magnetic moments. The existence of these two metallic regions in the phase diagram is a manifestation of the Pomeranchuk effect. The boundary between these regions reveals itself as a kink in the curve v c (T) and as a maximum in the T dependence of the double occupancy D. Our calculated D for different values of v, U, and T are in semiquantitative agreement with Monte Carlo results.