Cluster Perturbation Theory (CPT) is a technique for computing the spectral function of fermionic models with local interactions. By combining the solution of the model on a finite cluster with perturbation theory on intra-cluster hoppings, CPT provi
des access to single-particle properties with arbitrary momentum resolution while incurring low computational cost. Here, we introduce Determinantal Quantum Monte Carlo (DQMC) as a solver for CPT. Compared to the standard solver, exact diagonalization (ED), the DQMC solver reduces finite size effects through utilizing larger clusters, allows study of temperature dependence, and enables large-scale simulations of a greater set of models. We discuss the implementation of the DQMC solver for CPT and benchmark the CPT+DQMC method for the attractive and repulsive Hubbard models, showcasing its advantages over standard DQMC and CPT+ED simulations.
In the hole-doped cuprates, the pseudogap refers to a suppression of the density of states at low energies, in the absence of superconducting long-range order. Numerous calculations of the Hubbard model show a pseudogap in the single-particle spectra
, with striking similarities to photoemission and tunneling experiments on cuprates. However, no clear mechanism has been established. Here, we solve the Hubbard model on $2times2$ clusters by exact diagonalization, with integration over twisted boundary conditions. A pseudogap is found in the single-particle density of states with the following characteristics: a decreasing energy scale and onset temperature for increased hole-doping, closure at a critical hole doping near 15%, absence upon electron-doping, particle-hole asymmetry indicated by the location of the gap center, and persistence in the strong-coupling limit of $U/t to infty$. Studying the many-body excitation spectrum reveals that the pseudogap in single-particle spectra is due to orthogonality between bare electrons and the lowest energy excitations for $U/t gtrsim 8$.
We employ Momentum-Resolved Electron Energy Loss Spectroscopy (M-EELS) on Bi2.1Sr1.9CaCu2O8+x to resolve the issue of the kink feature in the electron dispersion widely observed in the cuprates. To this end, we utilize the GW approximation to relate
the density response function measured in in M-EELS to the self-energy, isolating contributions from phonons, electrons, and the momentum dependence of the effective interaction to the decay rates. The phononic contributions, present in the M-EELS spectra due to electron-phonon coupling, lead to kink features in the corresponding single-particle spectra at energies between 40 meV and 80 meV, independent of the doping level. We find that a repulsive interaction constant in momentum space is able to yield the kink attributed to phonons in ARPES. Hence, our analysis of the M-EELS spectra points to local repulsive interactions as a factor that enhances the spectroscopic signatures of electron-phonon coupling in cuprates. We conclude that the strength of the kink feature in cuprates is determined by the combined action of electron-phonon coupling and electron-electron interactions.
Because the cuprate superconductors are doped Mott insulators, it would be advantageous to solve even a toy model that exhibits both Mottness and superconductivity. We consider the Hatsugai-Kohmoto model, an exactly solvable system that is a prototyp
ical Mott insulator above a critical interaction strength at half filling. Upon doping or reducing the interaction strength, our exact calculations show that the system becomes a non-Fermi liquid metal with a superconducting instability. In the presence of a weak pairing interaction, the instability produces a thermal transition to a superconducting phase, which is distinct from the BCS state, as evidenced by a gap-to-transition temperature ratio exceeding the universal BCS limit. The elementary excitations of this superconductor are not Bogoliubov quasiparticles but rather superpositions of doublons and holons, composite excitations signaling that the superconducting ground state of the doped Mott insulator inherits the non-Fermi liquid character of the normal state. An unexpected feature of this model is that it exhibits a superconductivity-induced transfer of spectral weight from high to low energies as seen in the cuprates as well as a suppression of the superfluid density relative to that in BCS theory.
Strange or bad metallic transport, defined by its incompatibility with conventional quasiparticle pictures, is a theme common to strongly correlated materials and ubiquitous in many high temperature superconductors. The Hubbard model represents a min
imal starting point for modeling strongly correlated systems. Here we demonstrate strange metallic transport in the doped two-dimensional Hubbard model using determinantal quantum Monte Carlo calculations. Over a wide range of doping, we observe resistivities exceeding the Mott-Ioffe-Regel limit with linear temperature dependence. The temperatures of our calculations extend to as low as 1/40 the non-interacting bandwidth, placing our findings in the degenerate regime relevant to experimental observations of strange metallicity. Our results provide a foundation for connecting theories of strange metals to models of strongly correlated materials.
A microscopic understanding of the strongly correlated physics of the cuprates must account for the translational and rotational symmetry breaking that is present across all cuprate families, commonly in the form of stripes. Here we investigate emerg
ence of stripes in the Hubbard model, a minimal model believed to be relevant to the cuprate superconductors, using determinant quantum Monte Carlo (DQMC) simulations at finite temperatures and density matrix renormalization group (DMRG) ground state calculations. By varying temperature, doping, and model parameters, we characterize the extent of stripes throughout the phase diagram of the Hubbard model. Our results show that including the often neglected next-nearest-neighbor hopping leads to the absence of spin incommensurability upon electron-doping and nearly half-filled stripes upon hole-doping. The similarities of these findings to experimental results on both electron and hole-doped cuprate families support a unified description across a large portion of the cuprate phase diagram.
Evidence for the presence of high energy magnetic excitations in overdoped La$_{2-x}$Sr$_x$CuO$_4$ (LSCO) has raised questions regarding the role of spin-fluctuations in the pairing mechanism. If they remain present in overdoped LSCO, why does $T_c$
decrease in this doping regime? Here, using results for the dynamic spin susceptibility ${rm Im}chi(q,omega)$ obtained from a determinantal quantum Monte Carlo (DQMC) calculation for the Hubbard model we address this question. We find that while high energy magnetic excitations persist in the overdoped regime, they lack the momentum to scatter pairs between the anti-nodal regions. It is the decrease in the spectral weight at large momentum transfer, not observed by resonant inelastic X-ray scattering (RIXS), which leads to a reduction in the $d$-wave spin-fluctuation pairing strength.
Upon doping, Mott insulators often exhibit symmetry breaking where charge carriers and their spins organize into patterns known as stripes. For high-Tc superconducting cuprates, stripes are widely suspected to exist in a fluctuating form. Here, we us
e numerically exact determinant quantum Monte Carlo calculations to demonstrate dynamical stripe correlations in the three-band Hubbard model, which represents the local electronic structure of the Cu-O plane. Our results, which are robust to varying parameters, cluster size, and boundary condition, strongly support the interpretation of a variety of experimental observations in terms of the physics of fluctuating stripes, including the hourglass magnetic dispersion and the Yamada plot of incommensurability vs. doping. These findings provide a novel perspective on the intertwined orders emerging from the cuprates normal state.
