We examine the phase diagram of the extended Hubbard model on a square lattice, for both attractive and repulsive nearest-neighbor interactions, using CDMFT+HFD, a combination of Cluster Dynamical Mean Field theory (CDMFT) and a Hartree-Fock mean-field decoupling of the inter-cluster extended interaction. For attractive non-local interactions, this model exhibits a region of phase separation near half-filling, in the vicinity of which we find islands of d-wave superconductivity, decaying rapidly as a function of doping, with disconnected regions of extended s-wave order at smaller (higher) electron densities. On the other hand, when the extended interaction is repulsive, a Mott insulating state at half-filling is destabilized by hole doping, in the strong-coupling limit, in favor of d-wave superconductivity. At the particle-hole invariant chemical potential, we find a first-order phase transition from antiferromagnetism (AF) to d-wave superconductivity as a function of the attractive nearest-neighbor interaction, along with a deviation of the density from the half-filled limit. A repulsive extended interaction instead favors charge-density wave (CDW) order at half-filling.
We examine the phase diagram of the extended Hubbard model on a square lattice, for both attractive and repulsive nearest-neighbor interactions, using CDMFT+HFD, a combination of Cluster Dynamical Mean Field theory (CDMFT) and a Hartree-Fock mean-field decoupling of the inter-cluster extended interaction. For attractive non-local interactions, this model exhibits a region of phase separation near half-filling, in the vicinity of which we find islands of d-wave superconductivity, decaying rapidly as a function of doping, with disconnected regions of extended s-wave order at smaller (higher) electron densities. On the other hand, when the extended interaction is repulsive, a Mott insulating state at half-filling is destabilized by hole doping, in the strong-coupling limit, in favor of d-wave superconductivity. At the particle-hole invariant chemical potential, we find a first-order phase transition from antiferromagnetism (AF) to d-wave superconductivity as a function of the attractive nearest-neighbor interaction, along with a deviation of the density from the half-filled limit. A repulsive extended interaction instead favors charge-density wave (CDW) order at half-filling.
The nearest-neighbor superexchange-mediated mechanism for ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$ superconductivity in the one-band Hubbard model faces the challenge that nearest-neighbor Coulomb repulsion can be larger than superexchange. To answer this question, we use cellular dynamical mean-field theory (CDMFT) with a continuous-time quantum Monte Carlo solver to determine the superconducting phase diagram as a function of temperature and doping for on-site repulsion U=9t and nearest-neighbor repulsion V=0,2t,4t. In the underdoped regime, V increases the CDMFT superconducting transition temperature Tcd even though it decreases the superconducting order parameter at low temperature for all dopings. However, in the overdoped regime V decreases Tcd. We gain insight into these paradoxical results through a detailed study of the frequency dependence of the anomalous spectral function, extracted at finite temperature via the MaxEntAux method for analytic continuation. A systematic study of dynamical positive and negative contributions to pairing reveals that even though V has a high-frequency depairing contribution, it also has a low frequency pairing contribution since it can reinforce superexchange through $J=4{t}^{2}/(U\ensuremath{-}V)$. Retardation is thus crucial to understanding pairing in doped Mott insulators, as suggested by previous zero-temperature studies. We also comment on the tendency to charge order for large V and on the persistence of d-wave superconductivity over extended-s or s+d wave.
We consider an efficient scheme to simulate fermionic Hubbard models with nonlocal density-density interactions in two dimensions, based on bond-centered auxiliary-field quantum Monte Carlo. The simulations are shown to be sign-problem-free within a finite, restricted parameter range. Using this approach, we first study the Hubbard model on the half-filled square lattice bilayer, including an interlayer repulsion term in addition to the local repulsion, and present the ground-state phase diagram within the accessible parameter region. Starting from the antiferromagnetically ordered state in the absence of interlayer repulsion, the interlayer interactions are found to destabilize the antiferromagnetic order, leading to a band insulator state. Moreover, for sufficiently strong interlayer tunneling, we also observe the emergence of a direct dimer product state of mixed D-Mott and S-Mott character along the equal coupling line. We discuss the stability range of this state within strong-coupling perturbation theory. Furthermore, we consider the Hubbard model on the honeycomb lattice with next-nearest-neighbor interactions. Such an interaction is found to enhance both charge density and spin-current correlations within the semimetallic region. However, inside the accessible parameter region, they do not stabilize long-ranged charge density wave order nor a quantum spin Hall state, and the only insulating state that we observe exhibits long-range antiferromagnetism.
High temperature superconductivity in copper oxides (cuprates) represents one of the most striking phenomena in condensed matter physics. These materials host superconducting states with the highest known critical temperatures, emerging upon doping antiferromagnetic insulator compounds. The antiferromagnetic (AF) insulating and superconducting (SC) phases defy standard descriptions in terms of band picture and conventional phonon-mediated superconductivity. In fact, these phenomena are genuine signatures of strong electronic correlations which hinder effective single-particle descriptions of the many-electron system. After 40 years since their discovery, it is widely accepted that the Hubbard model on the square lattice provides a minimal description of the low-energy electronic properties of cuprates. The model describes the competition between the electronic kinetic energy, modelled by nearest-neighbours hopping t, and the Coulomb repulsion, modelled by a purely local interaction potential U. At half-filling, this competition results in the so-called "Mott localisation" phenomenon, which distinguishes the metallic phase for U/t « 1 from an insulating phase for U/t » 1. The latter is characterised by the almost complete localisation of electrons due to the high energetic cost of hopping towards occupied sites. Interestingly, this basic competition is intimately related to AF ordering and, upon doping, to the formation of the SC state with characteristic d-wave symmetry. The fact that a purely repulsive interaction U is responsible for the formation of a SC state is remarkable, being superconductivity typically associated with an effective attractive interaction. In fact, a d-wave SC state can also be obtained, in general, by considering an attractive potential between electrons of neighbouring sites. Motivated by this observation, in this thesis, we study a modified version of the Hubbard model, dubbed as the "extended Hubbard model", in which we include a nearest-neighbours attractive potential -V. The general goal of this work is to explore the interplay between the purely local repulsion and the non-local attraction. This is not only motivated by the fact that both interactions are, in principle, responsible for d-wave superconductivity, but also by recent experimental observations which strongly suggest that the term -V needs to be included in order to accurately describe the magnetic spectra of some specific compounds. In this thesis, we study the extended Hubbard model using the Hartree-Fock method. Namely, we approximate the ground state by implementing a variational search in the subspace of the many-body Hilbert space of gaussian states. Despite its simplicity, the Hartree-Fock method provides a powerful framework to study competition among different phases, originating from different symmetry breaking mechanisms. In particular, we study the competition between three different phases: the normal phase, characterised by no symmetry breaking, the AF phase, characterised by a lowering of SU(2) symmetry, and the SC phase, characterised by the breaking of U(1) charge conservation. This thesis work led to a number of original results, which can be summarised by the general observation that the non-local attraction -V acts cooperatively with the local repulsion U for the stabilisation of the various phases hosted by the Hubbard model. For the normal phase, we find that the non-local attraction acts as a source of Mott localisation, leading to an effective reduction of the electronic bandwidth. In particular, for values of V larger than a critical value and close to half-filling, we find a strong Mott localisation characterised by the complete flattening of the electronic bands. Interestingly, this result shows an example of Mott localisation that can be described at the level of the Hartree-Fock approximation, whereas Mott localisation in the conventional Hubbard model usually requires more sophisticated approximation schemes. For the AF phase, we restrict our analysis to the half-filled case. We find that the pairing-induced Mott localisation leads to a significant enhancement of the magnetisation, primarily originating from the local repulsion U. Finally, we focus on the superconducting phase and derive self-consistency equations in both singlet and triplet Cooper pairing channels, showing how Cooper pairing can be mediated simultaneously by both the local repulsion and the non-local attraction. Therefore, we numerically investigate the singlet channel and find evidence of "symmetry mixing", namely, the formation of a superconducting gap that is neither symmetric nor anti-symmetric under rotations of π/2 in the plane. We show that symmetry mixing is connected to the emergence of a finite nematic order parameter, that signals the breaking of planar rotational symmetry. This nematic symmetry breaking results in the fact that the renormalisation of bands becomes anisotropic, producing a difference in electronic kinetics in the x and y directions. To summarise, the original results obtained in this thesis show that the non-local pairing provides a different route towards the formation of key strongly correlated phenomena, including Mott localisation and nematic superconductivity, and call for future investigation using more sophisticated methods such as, for example, the Dynamical Mean Field Theory.
We investigate the phase diagram of the extended Hubbard model on a half-filled honeycomb lattice through determinant quantum Monte Carlo simulations. At zero temperature, the repulsive Hubbard model undergoes a transition to an antiferromagnetic state only for sufficiently strong on-site coupling U. However, the inclusion of a repulsive extended interaction V between fermions in nearest-neighbor orbitals induces a charge-ordered phase. For attractive interactions, we also identify superconducting and phase-separated phases. We have determined a complete phase diagram with the aid of structure factors associated with various orderings, combined with measurements of double occupancy and the average sign of fermionic determinants. We have established that the semimetallic phase forms a zone from which ordered phases are excluded, preventing the stabilization of a superconducting phase with d-wave symmetry, i.e., only s-wave pairing is allowed.
A slave-spin representation of fermion operators has recently been proposed for the half-filled single and multiband Hubbard model. We show that with the addition of a gauge variable, the formalism can be extended to finite doping. We solve the resulting spin problem using the cluster mean-field approximation. This approximation takes short-range correlations into account by exact diagonalization on the cluster, whereas long-range correlations beyond the size of clusters are treated at the mean-field level. In the limit where the cluster has only one site and the interaction strength U is infinite, this approach reduces to the Gutzwiller approximation. There are some qualitative differences when the size of the cluster is finite. We first compute the critical U for the Mott transition as a function of a frustrating nearest-neighbor interaction on lattices relevant for various correlated systems, namely, the cobaltates, the layered organic superconductors and the high-temperature superconductors. For the triangular lattice, we also study the extended Hubbard model with nearest-neighbor repulsion. In addition to a uniform metallic state, we find a $\sqrt{(3)}\ifmmode\times\else\texttimes\fi{}\sqrt{(3)}$ charge density wave in a broad doping regime, including commensurate ones. We find that in the large U limit, intersite Coulomb repulsion V strongly suppresses the single-particle weight of the metallic state.
Using a cluster extension of the dynamical mean-field theory (CDMFT) we map out the magnetic phase diagram of the anisotropic square lattice Hubbard model with nearest-neighbor intrachain t and interchain ${t}_{\ensuremath{\perp}}$ hopping amplitudes at half filling. A fixed value of the next-nearest-neighbor hopping ${t}^{\ensuremath{'}}=\ensuremath{-}{t}_{\ensuremath{\perp}}/2$ removes the nesting property of the Fermi surface and stabilizes a paramagnetic metal phase in the weak-coupling regime. In the isotropic and moderately anisotropic regions, a growing spin entropy in the metal phase is quenched out at a critical interaction strength by the onset of long-range antiferromagnetic (AF) order of preformed local moments. It gives rise to a first-order metal-insulator transition consistent with the Mott-Heisenberg picture. In contrast, a strongly anisotropic regime ${t}_{\ensuremath{\perp}}/t\ensuremath{\lesssim}0.3$ displays a quantum critical behavior related to the continuous transition between an AF metal phase and the AF insulator. Hence, within the present framework of CDMFT, the opening of the charge gap is magnetically driven as advocated in the Slater picture. We also discuss how the lattice-anisotropy-induced evolution of the electronic structure on a metallic side of the phase diagram is tied to the emergence of quantum criticality.