The effect of electron-electron interactions on Dirac fermions, and the possibility of an intervening spin-liquid phase between the semimetal and antiferromagnetic (AF) regimes, has been a focus of intense quantum simulation effort over the last five years. We use determinant quantum Monte Carlo simulations to study the Holstein model on a honeycomb lattice and explore the role of electron- phonon interactions on Dirac fermions. We show that they give rise to charge-density-wave (CDW) order and present evidence that this occurs only above a finite critical interaction strength. We evaluate the temperature for the transition into the CDW which, unlike the AF transition, can occur at finite values owing to the discrete nature of the broken symmetry.

The Holstein model describes the interaction between fermions and a collection of local (dispersionless) phonon modes. In the dilute limit, the phonon degrees of freedom dress the fermions, giving rise to polaron and bipolaron formation. At higher densities, the phonons mediate collective superconducting (SC) and charge-density wave (CDW) phases. Quantum Monte Carlo (QMC) simulations have considered both these limits but have not yet focused on the physics of more general phonon spectra. Here we report QMC studies of the role of phonon dispersion on SC and CDW order in such models. We quantify the effect of finite phonon bandwidth and curvature on the critical temperature $T_{\mathrm{cdw}}$ for CDW order and also uncover several novel features of diagonal long-range order in the phase diagram, including a competition between charge patterns at momenta $\mathbf{q}=(\pi ,\pi )$ and $\mathbf{q}=(0,\pi )$ which lends insight into the relationship between Fermi surface nesting and the wave vector at which charge order occurs. We also demonstrate SC order at half filling in situations where a nonzero bandwidth sufficiently suppresses $T_{\mathrm{cdw}}$.

The electron–electron and electron–phonon interactions play an important role in correlated materials, being key features for spin, charge and pair correlations. Thus, here we investigate their effects 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, favouring (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.

We use determinant quantum Monte Carlo (DQMC), in combination with the principal component analysis (PCA) approach to unsupervised learning, to extract information about phase transitions in several of the most fundamental Hamiltonians describing strongly correlated materials. We first explore the zero-temperature antiferromagnet to singlet transition in the periodic Anderson model, the Mott insulating transition in the Hubbard model on a honeycomb lattice, and the magnetic transition in the 1/6-filled Lieb lattice. We then discuss the prospects for learning finite temperature superconducting transitions in the attractive Hubbard model, for which there is no sign problem. Finally, we investigate finite temperature charge density wave (CDW) transitions in the Holstein model, where the electrons are coupled to phonon degrees of freedom, and carry out a finite size scaling analysis to determine $T_c$. We examine the different behaviors associated with Hubbard-Stratonovich auxiliary field configurations on both the entire space-time lattice and on a single imaginary time slice, or other quantities, such as equal-time Green's and pair-pair correlation functions.

Recent experiments with ultracold fermionic atoms in optical lattices have provided a tuneable and clean realization of the attractive Hubbard model (AHM). In view of this, several physical properties may be thoroughly studied across the crossover between weak [Bardeen-Cooper-Schrieffer, (BCS)] and strong [Bose-Einstein condensation (BEC)] couplings. Here we report on extensive determinant Quantum Monte Carlo (DQMC) studies of the AHM on a square lattice from which several different quantities have been calculated and should be useful as a road map to experiments. We have obtained a detailed phase diagram for the critical superconducting temperature $T_c$ in terms of the band filling $\langle n\rangle$ and interaction strength $U$ from which we pinpoint a somewhat wide region $\left|U\right|/t\approx 5\pm 1$ ($t$ is the hopping amplitude) and $\langle n\rangle \approx 0.79\pm 0.09$ leading to a maximum $T_c\approx 0.16t$. Two additional temperature scales, namely, pairing $T_p$ and degeneracy $T_d$ have been highlighted: The former sets the scale for pair formation (believed to be closely related to the scale for the gap of spin excitations in cuprates), whereas the latter sets the scale for dominant quantum effects. Our DQMC data for the distribution of doubly occupied sites for the momentum distribution function, and for the quasiparticle weight show distinctive features on both sides of the BCS-BEC crossover, being also suggestive of an underlying crossover between Fermi- and non-Fermi liquid behaviors.

The noninteracting electronic structures of tight-binding models on bipartite lattices with unequal numbers of sites in the two sublattices have a number of unique features, including the presence of spatially localized eigenstates and flat bands. When a uniform on-site Hubbard interaction $U$ is turned on, Lieb proved rigorously that at half-filling ($\rho =1$) the ground state has a nonzero spin. In this paper we consider a “${\mathrm{CuO}}_2$ lattice” (also known as “Lieb lattice,” or as a decorated square lattice), in which “$d$ orbitals” occupy the vertices of the squares, while “$p$ orbitals” lie halfway between two $d$ orbitals; both $d$ and $p$ orbitals can accommodate only up to two electrons. We use exact determinant quantum Monte Carlo (DQMC) simulations to quantify the nature of magnetic order through the behavior of correlation functions and sublattice magnetizations in the different orbitals as a function of $U$ and temperature; we have also calculated the projected density of states, and the compressibility. We study both the homogeneous (H) case, $U_d=U_p$, originally considered by Lieb, and the inhomogeneous (IH) case, $U_d\ne U_p$. For the H case at half-filling, we found that the global magnetization rises sharply at weak coupling, and then stabilizes towards the strong-coupling (Heisenberg) value, as a result of the interplay between the ferromagnetism of like sites and the antiferromagnetism between unlike sites; we verified that the system is an insulator for all $U$. For the IH system at half-filling, we argue that the case $U_p\ne U_d$ falls under Lieb's theorem, provided they are positive definite, so we used DQMC to probe the cases $U_p=0,U_d=U$ and $U_p=U,U_d=0$. We found that the different environments of $d$ and $p$ sites lead to a ferromagnetic insulator when $U_d=0$; by contrast, $U_p=0$ leads to to a metal without any magnetic ordering. In addition, we have also established that at density $\rho =1/3$, strong antiferromagnetic correlations set in, caused by the presence of one fermion on each $d$ site.