No Arabic abstract
Fixed-node Greens function Monte Carlo calculations have been performed for very large 16x6 2D Hubbard lattices, large interaction strengths U=10,20, and 40, and many (15-20) densities between empty and half filling. The nodes were fixed by a simple Slater-Gutzwiller trial wavefunction. For each value of U we obtained a sequence of ground-state energies which is consistent with the possibility of a phase separation close to half-filling, with a hole density in the hole-rich phase which is a decreasing function of U. The energies suffer, however, from a fixed-node bias: more accurate nodes are needed to confirm this picture. Our extensive numerical results and their test against size, shell, shape and boundary condition effects also suggest that phase separation is quite a delicate issue, on which simulations based on smaller lattices than considered here are unlikely to give reliable predictions.
We explore the Matsubara quasiparticle fraction and the pseudogap of the two-dimensional Hubbard model with the dynamical cluster quantum Monte Carlo method. The character of the quasiparticle fraction changes from non-Fermi liquid, to marginal Fermi liquid to Fermi liquid as a function of doping, indicating the presence of a quantum critical point separating non-Fermi liquid from Fermi liquid character. Marginal Fermi liquid character is found at low temperatures at a very narrow range of doping where the single-particle density of states is also symmetric. At higher doping the character of the quasiparticle fraction is seen to cross over from Fermi Liquid to Marginal Fermi liquid as the temperature increases.
We study a two-dimensional model of an isolated narrow topological band at partial filling with local attractive interactions. Numerically exact quantum Monte Carlo calculations show that the ground state is a superconductor with a critical temperature that scales nearly linearly with the interaction strength. We also find a broad pseudogap regime at temperatures above the superconducting phase that exhibits strong pairing fluctuations and a tendency towards electronic phase separation; introducing a small nearest neighbor attraction suppresses superconductivity entirely and results in phase separation. We discuss the possible relevance of superconductivity in this unusual regime to the physics of flat band moir{e} materials, and as a route to designing higher temperature superconductors.
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.
We report large scale determinant Quantum Monte Carlo calculations of the effective bandwidth, momentum distribution, and magnetic correlations of the square lattice fermion Hubbard Hamiltonian at half-filling. The sharp Fermi surface of the non-interacting limit is significantly broadened by the electronic correlations, but retains signatures of the approach to the edges of the first Brillouin zone as the density increases. Finite size scaling of simulations on large lattices allows us to extract the interaction dependence of the antiferromagnetic order parameter, exhibiting its evolution from weak coupling to the strong coupling Heisenberg limit. Our lattices provide improved resolution of the Greens function in momentum space, allowing a more quantitative comparison with time-of-flight optical lattice experiments.
Spontaneous phase separation instabilities with the formation of various types of charge and spin pairing (pseudo)gaps in $U>0$ Hubbard model including the {it next nearest neighbor coupling} are calculated with the emphasis on the two-dimensional (square) lattices generated by 8- and 10-site Betts unit cells. The exact theory yields insights into the nature of quantum critical points, continuous transitions, dramatic phase separation instabilities and electron condensation in spatially inhomogeneous systems. The picture of coupled anti-parallel (singlet) spins and paired charged holes suggests full Bose condensation and coherent pairing in real space at zero temperature of electrons complied with the Bose-Einstein statistics. Separate pairing of charge and spin degrees at distinct condensation temperatures offers a new route to superconductivity different from the BCS scenario. The conditions for spin liquid behavior coexisting with unsaturated and saturated Nagaoka ferromagnetism due to spin-charge separation are established. The phase separation critical points and classical criticality found at zero and finite temperatures resemble a number of inhomogeneous, coherent and incoherent nanoscale phases seen near optimally doped high-$T_c$ cuprates, pnictides and CMR nanomaterials.