No Arabic abstract
The Kondo and Periodic Anderson models describe many of the qualitative features of local moments coupled to a conduction band, and thereby the physics of materials such as the heavy fermions. In particular, when the exchange coupling $J$ or hybridization $V$ between the moments and the electrons of the metallic band is large, singlets form, quenching the magnetism. In the opposite, small $J$ or $V$, limit, the moments survive, and the conduction electrons mediate an effective interaction which can trigger long range, often antiferromagnetic, order. In the case of the Kondo model, where the moments are described by local spins, Nozi`eres considered the possibility that the available conduction electrons within the Kondo temperature of the Fermi surface would be insufficient in number to accomplish the screening. Much effort in the literature has been devoted to the study of the temperature scales in the resulting `exhaustion problem, and how the `coherence temperature where a heavy Fermi liquid forms is related to the Kondo temperature. In this paper, we study a version of the Periodic Anderson model in which some of the conduction electrons are removed in a way which avoids the fermion sign problem and hence allows low temperature Quantum Monte Carlo simulations which can access both singlet formation and magnetic ordering temperature scales. We are then able to focus on a somewhat different aspect of exhaustion physics than previously considered: the effect of dilution on the critical $V$ for the singlet-antiferromagnetic transition.
We have performed numerical studies of the Hubbard-Holstein model in two dimensions using determinant quantum Monte Carlo (DQMC). Here we present details of the method, emphasizing the treatment of the lattice degrees of freedom, and then study the filling and behavior of the fermion sign as a function of model parameters. We find a region of parameter space with large Holstein coupling where the fermion sign recovers despite large values of the Hubbard interaction. This indicates that studies of correlated polarons at finite carrier concentrations are likely accessible to DQMC simulations. We then restrict ourselves to the half-filled model and examine the evolution of the antiferromagnetic structure factor, other metrics for antiferromagnetic and charge-density-wave order, and energetics of the electronic and lattice degrees of freedom as a function of electron-phonon coupling. From this we find further evidence for a competition between charge-density-wave and antiferromagnetic order at half-filling.
The superconducting (SC) and charge-density-wave (CDW) susceptibilities of the two dimensional Holstein model are computed using determinant quantum Monte Carlo (DQMC), and compared with results computed using the Migdal-Eliashberg (ME) approach. We access temperatures as low as 25 times less than the Fermi energy, $E_F$, which are still above the SC transition. We find that the SC susceptibility at low $T$ agrees quantitatively with the ME theory up to a dimensionless electron-phonon coupling $lambda_0 approx 0.4$ but deviates dramatically for larger $lambda_0$. We find that for large $lambda_0$ and small phonon frequency $omega_0 ll E_F$ CDW ordering is favored and the preferred CDW ordering vector is uncorrelated with any obvious feature of the Fermi surface.
We characterize the three-orbital Hubbard model using state-of-the-art determinant quantum Monte Carlo (DQMC) simulations with parameters relevant to the cuprate high-temperature superconductors. The simulations find that doped holes preferentially reside on oxygen orbitals and that the ({pi},{pi}) antiferromagnetic ordering vector dominates in the vicinity of the undoped system, as known from experiments. The orbitally-resolved spectral functions agree well with photoemission spectroscopy studies and enable identification of orbital content in the bands. A comparison of DQMC results with exact diagonalization and cluster perturbation theory studies elucidates how these different numerical techniques complement one another to produce a more complete understanding of the model and the cuprates. Interestingly, our DQMC simulations predict a charge-transfer gap that is significantly smaller than the direct (optical) gap measured in experiment. Most likely, it corresponds to the indirect gap that has recently been suggested to be on the order of 0.8 eV, and demonstrates the subtlety in identifying charge gaps.
Determinant Quantum Monte Carlo (DQMC) is used to determine the pairing and magnetic response for a Hubbard model built up from four-site clusters -a two-dimensional square lattice consisting of elemental 2x2 plaquettes with hopping $t$ and on-site repulsion $U$ coupled by an inter-plaquette hopping $t leq t$. Superconductivity in this geometry has previously been studied by a variety of analytic and numeric methods, with differing conclusions concerning whether the pairing correlations and transition temperature are raised near half-filling by the inhomogeneous hopping or not. For $U/t=4$, DQMC indicates an optimal $t/t approx 0.4$ at which the pairing vertex is most attractive. The optimal $t/t$ increases with $U/t$. We then contrast our results for this plaquette model with a Hamiltonian which instead involves a regular pattern of site energies whose large site energy limit is the three band CuO$_2$ model; we show that there the inhomogeneity rapidly, and monotonically, suppresses pairing.
Layered antiferromagnetic spin density wave (LAF) state is one of the plausible ground states of charge neutral Bernal stacked bilayer graphene. In this paper, we use determinant quantum Monte Carlo method to study the effect of the electric field on the magnetic order in bilayer Hubbard model on a honeycomb lattice. Our results qualitatively support the LAF ground state found in the mean field theory. The obtained magnetic moments, however, are much smaller than what are estimated in the mean field theory. As electric field increases, the magnetic order parameter rapidly decreases.