We study the single-orbital Hubbard model on the 1/5-depleted square-lattice geometry, which arises in such diverse systems as the spin-gap magnetic insulator CaV4O9 and ordered-vacancy iron selenides, presenting new issues regarding the origin of both magnetic ordering and superconductivity in these materials. We find a rich phase diagram that includes a plaquette singlet phase, a dimer singlet phase, a Néel and a block-spin antiferromagnetic phase, and stripe phases. Quantum Monte Carlo simulations show that the dominant pairing correlations at half filling change character from d wave in the plaquette phase to extended s wave upon transition to the Néel phase. These findings have intriguing connections to iron-based superconductors, and suggest that some physics of multiorbital systems can be captured by a single-orbital model at different dopings.
We study the single-orbital Hubbard model on the 1/5-depleted square-lattice geometry, which arises in such diverse systems as the spin-gap magnetic insulator CaV4O9 and ordered-vacancy iron selenides, presenting new issues regarding the origin of both magnetic ordering and superconductivity in these materials. We find a rich phase diagram that includes a plaquette singlet phase, a dimer singlet phase, a Néel and a block-spin antiferromagnetic phase, and stripe phases. Quantum Monte Carlo simulations show that the dominant pairing correlations at half filling change character from d wave in the plaquette phase to extended s wave upon transition to the Néel phase. These findings have intriguing connections to iron-based superconductors, and suggest that some physics of multiorbital systems can be captured by a single-orbital model at different dopings.
We consider the one-band Hubbard model on the square lattice by using variational and Green's function Monte Carlo methods, where the variational states contain Jastrow and backflow correlations on top of an uncorrelated wave function that includes BCS pairing and magnetic order. At half-filling, where the ground state is antiferromagnetically ordered for any value of the on-site interaction U, we can identify a hidden critical point UMott, above which a finite BCS pairing is stabilized in the wave function. The existence of this point is reminiscent of the Mott transition in the paramagnetic sector and determines a separation between a Slater insulator (at small values of U), where magnetism induces a potential energy gain, and a Mott insulator (at large values of U), where magnetic correlations drive a kinetic energy gain. Most importantly, the existence of UMott has crucial consequences when doping the system: We observe a tendency for phase separation into hole-rich and hole-poor regions only when doping the Slater insulator, while the system is uniform by doping the Mott insulator. Superconducting correlations are clearly observed above UMott, leading to the characteristic dome structure in doping. Furthermore, we show that the energy gain due to the presence of a finite BCS pairing above UMott shifts from the potential to the kinetic sector by increasing the value of the Coulomb repulsion.
The cluster slave-spin method is employed to investigate systematically the ground state properties of the Hubbard model on a square lattice with doping $\ensuremath{\delta}$ and coupling strength U being its parameters. At half-filling, a relation between the staggered magnetization M and the antiferromagnetic (AFM) gap ${\mathrm{\ensuremath{\Delta}}}_{\text{AFM}}$ is established in the small U limit to compare with that from the Hartree-Fock theory, and a first-order metal-insulator Mott transition in the paramagnetic state is substantiated, which is characterized by discontinuities and hystereses at UMott=10t. The interaction Uc for the crossover in the AFM state, separating the weak- and strong-coupling regimes, is found to remain almost unchanged with large dopings, and smaller than UMott at half-filling because of long range AFM correlations. Finally, an overall phase diagram in the $U\ensuremath{-}\ensuremath{\delta}$ plane is presented, which is composed of four regimes: the AFM insulator at half-filling, the AFM metal with the compressibility $\ensuremath{\kappa}>0$ or $\ensuremath{\kappa}<0$, and the paramagnetic metal, as well as three phase transitions: (i) From the AFM metal to the paramagnetic metal, (ii) between the AFM metal phases with positive and negative $\ensuremath{\kappa}$, and (iii) separating the AFM insulating phase at $\ensuremath{\delta}=0$ from the AFM metal phase for $\ensuremath{\delta}>0$.
We study the half-filled Hubbard model on a one-parameter family of vortex-full square lattices ranging from the isotropic case to weakly coupled Hubbard dimers. The ground-state phase diagram consists of four phases: A semimetal and a band insulator which are connected to the weak-coupling limit, and a magnetically ordered N\'eel phase and a valence bond solid (VBS) which are linked to the strong-coupling Mott limit. The phase diagram is obtained by quantum Monte Carlo (QMC) and continuous unitary transformations (CUTs). The CUT is performed in a two-step process: Nonperturbative graph-based CUTs are used in the Mott insulating phase to integrate out charge fluctuations. The resulting effective spin model is tackled by perturbative CUTs about the isolated dimer limit yielding the breakdown of the VBS by triplon condensation. We find three scenarios when varying the interaction for a fixed anisotropy of hopping amplitudes: (i) one direct phase transition from N\'eel to semimetal, (ii) two phase transitions VBS to N\'eel and N\'eel to semimetal, or (iii) a smooth crossover from VBS to the band insulator. Our results are consistent with the absence of spin-liquid phases in the whole phase diagram.
The recent discovery of superconductivity under high pressure in the ladder compound BaFe2S3 has opened a new field of research in iron-based superconductors with focus on quasi-one-dimensional geometries. In this publication, using the density matrix renormalization group technique, we study a two-orbital Hubbard model defined in one-dimensional chains. Our main result is the presence of hole binding tendencies at intermediate Hubbard U repulsion and robust Hund coupling JH/U=0.25. Binding does not occur either in weak coupling or at very strong coupling. The pair-pair correlations that are dominant near half-filling, or of similar strength as the charge and spin correlation channels, involve hole-pair operators that are spin singlets, use nearest-neighbor sites, and employ different orbitals for each hole. The Hund coupling strength, presence of robust magnetic moments, and antiferromagnetic correlations among them are important for the binding tendencies found here.
Phase diagram of the two-dimensional Hubbard-Holstein model: enhancement of s-wave pairing between charge and magnetic orders
We investigate the role of electron-electron and electron-phonon interactions in strongly correlated systems by performing unbiased quantum Monte Carlo simulations in the square lattice Hubbard-Holstein model at half-filling. We study the competition and interplay between antiferromagnetism (AFM) and charge-density wave (CDW), establishing its very rich phase diagram. In the region between AFM and CDW phases, we have found an enhancement of superconducting pairing correlations, favoring (nonlocal) s-wave pairs. Our study sheds light over past inconsistencies in the literature, in particular the emergence of CDW in the pure Holstein model case.
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.