No Arabic abstract
Using determinant Quantum Monte Carlo, we compare three methods of evaluating the dc Hall coefficient $R_H$ of the Hubbard model: the direct measurement of the off-diagonal current-current correlator $chi_{xy}$ in a system coupled to a finite magnetic field (FF), $chi_{xy}^{text{FF}}$; the three-current linear response to an infinitesimal field as measured in the zero-field (ZF) Hubbard Hamiltonian, $chi_{xy}^{text{ZF}}$; and the leading order of the recurrent expansion $R_H^{(0)}$ in terms of thermodynamic susceptibilities. The two quantities $chi_{xy}^{text{FF}}$ and $chi_{xy}^{text{ZF}}$ can be compared directly in imaginary time. Proxies for $R_H$ constructed from the three-current correlator $chi_{xy}^{text{ZF}}$ can be determined under different simplifying assumptions and compared with $R_H^{(0)}$. We find these different quantities to be consistent with one another, validating previous conclusions about the close correspondence between Fermi surface topology and the sign of $R_H$, even for strongly correlated systems. These various quantities also provide a useful set of numerical tools for testing theoretical predictions about the full behavior of the Hall conductivity for strong correlations.
Since its experimental discovery, many phenomenological theories successfully reproduced the rapid rise from $p$ to $1+p$ found in the Hall number $n_H$ at the critical doping $p^*$ of the pseudogap in superconducting cuprates. Further comparison with experiments is now needed in order to narrow down candidates. In this paper, we consider three previously successful phenomenological theories in a unified formalism---an antiferromagnetic mean field (AF), a spiral incommensurate antiferromagnetic mean field (sAF), and the Yang-Rice-Zhang (YRZ) theory. We find a rapid rise in the specific heat and a rapid drop in the Seebeck coefficient for increasing doping across the transition in each of those models. The predicted rises and drops are locked, not to~$p^*$, but to the doping where anti-nodal electron pockets, characteristic of each model, appear at the Fermi surface shortly before~$p^*$. While such electron pockets are still to be found in experiments, we discuss how they could provide distinctive signatures for each model. We also show that the range of doping where those electron pockets would be found is strongly affected by the position of the van~Hove singularity.
We study the ground state properties of the Hubbard model on a 4-leg cylinder with doped hole concentration per site $deltaleq 12.5%$ using density-matrix renormalization group. By keeping a large number of states for long system sizes, we find that the nature of the ground state is remarkably sensitive to the presence of next-nearest-neighbor hopping $t$. Without $t$ the ground state of the system corresponds with the insulating filled stripe phase with long-range charge-density-wave (CDW) order and short-range incommensurate spin correlations appears. However, for a small negative $t$ a phase characterized by coexisting algebraic d-wave superconducting (SC)- and algebraic CDW correlations. In addition, it shows short range spin- and fermion correlations consistent with a canonical Luther-Emery (LE) liquid, except that the charge- and spin periodicities are consistent with half-filled stripes instead of the $4 k_F$ and $2 k_F$ wavevectors generic for one dimensional chains. For a small positive $t$ yet another phase takes over showing similar SC and CDW correlations. However, the fermions are now characterized by a (near) infinite correlation length while the gapped spin system is characterized by simple staggered antiferromagnetic correlations. We will show that this is consistent with a LE formed from a weakly coupled (BCS like) d-wave superconductor on the ladder where the interactions have only the effect to stabilize a cuprate style magnetic resonance.
The repulsive Hubbard model has been immensely useful in understanding strongly correlated electron systems, and serves as the paradigmatic model of the field. Despite its simplicity, it exhibits a strikingly rich phenomenology which is reminiscent of that observed in quantum materials. Nevertheless, much of its phase diagram remains controversial. Here, we review a subset of what is known about the Hubbard model, based on exact results or controlled approximate solutions in various limits, for which there is a suitable small parameter. Our primary focus is on the ground state properties of the system on various lattices in two spatial dimensions, although both lower and higher dimensions are discussed as well. Finally, we highlight some of the important outstanding open questions.
We present determinant quantum Monte Carlo simulations of the hole-doped single-band Hubbard-Holstein model on a square lattice, to investigate how quasiparticles emerge when doping a Mott insulator (MI) or a Peierls insulator (PI). The MI regime at large Hubbard interaction $U$ and small relative electron-phonon coupling strength $lambda$ is quickly suppressed upon doping, by drawing spectral weight from the upper Hubbard band and shifting the lower Hubbard band towards the Fermi level, leading to a metallic state with emergent quasiparticles at the Fermi level. On the other hand, the PI regime at large $lambda$ and small $U$ persists out to relatively high doping levels. We study the evolution of the $d$-wave superconducting susceptibility with doping, and find that it increases with lowering temperature in a regime of intermediate values of $U$ and $lambda$.
The dualism between superconductivity and charge/spin modulations (the so-called stripes) dominates the phase diagram of many strongly-correlated systems. A prominent example is given by the Hubbard model, where these phases compete and possibly coexist in a wide regime of electron dopings for both weak and strong couplings. Here, we investigate this antagonism within a variational approach that is based upon Jastrow-Slater wave functions, including backflow correlations, which can be treated within a quantum Monte Carlo procedure. We focus on clusters having a ladder geometry with $M$ legs (with $M$ ranging from $2$ to $10$) and a relatively large number of rungs, thus allowing us a detailed analysis in terms of the stripe length. We find that stripe order with periodicity $lambda=8$ in the charge and $2lambda=16$ in the spin can be stabilized at doping $delta=1/8$. Here, there are no sizable superconducting correlations and the ground state has an insulating character. A similar situation, with $lambda=6$, appears at $delta=1/6$. Instead, for smaller values of dopings, stripes can be still stabilized, but they are weakly metallic at $delta=1/12$ and metallic with strong superconducting correlations at $delta=1/10$, as well as for intermediate (incommensurate) dopings. Remarkably, we observe that spin modulation plays a major role in stripe formation, since it is crucial to obtain a stable striped state upon optimization. The relevance of our calculations for previous density-matrix renormalization group results and for the two-dimensional case is also discussed.