We report diffusion quantum Monte Carlo calculations of three-dimensional Wigner crystals in the density range r_s=100-150. We have tested different types of orbital for use in the approximate wave functions but none improve upon the simple Gaussian form. The Gaussian exponents are optimized by directly minimizing the diffusion quantum Monte Carlo energy. We have carefully investigated and sought to minimize the potential biases in our Monte Carlo results. We conclude that the uniform electron gas undergoes a transition from a ferromagnetic fluid to a body-centered-cubic Wigner crystal at r_s=106+/-1. The diffusion quantum Monte Carlo results are compared with those from Hartree-Fock and Hartree theory in order to understand the role played by exchange and correlation in Wigner crystals. We also study floating Wigner crystals and give results for their pair-correlation functions.
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.
We probe the superconducting gap in the zero temperature ground state of an attractively interacting spin-imbalanced two-dimensional Fermi gas with Diffusion Monte Carlo. A condensate fraction at nonzero pair momentum evidences a spatially non-uniform superconducting order parameter. Comparison with exact diagonalisation studies confirms that the nonzero condensate fraction across a range of nonzero fermion pair momenta is consistent with non-exclusive pairing between majority and minority fermions, an extension beyond FFLO theory.
According to Landaus Fermi liquid theory, the main properties of the quasiparticle excitations of an electron gas are embodied in the effective mass $m^*$, which determines the energy of a single quasiparticle, and the Landau interaction function, which indicates how the energy of a quasiparticle is modified by the presence of other quasiparticles. This simple paradigm underlies most of our current understanding of the physical and chemical behavior of metallic systems. The quasiparticle effective mass of the three-dimensional homogeneous electron gas has been the subject of theoretical controversy and there is a lack of experimental data. In this work, we deploy diffusion Monte Carlo (DMC) methods to calculate $m^*$ as a function of density for paramagnetic and ferromagnetic three-dimensional homogeneous electron gases. The DMC results indicate that $m^*$ decreases when the density is reduced, especially in the ferromagnetic case. The DMC quasiparticle energy bands exclude the possibility of a reduction in the occupied bandwidth relative to that of the free-electron model at density parameter $r_s=4$, which corresponds to Na metal.
A many-flavor electron gas (MFEG) in a semiconductor with a valley degeneracy ranging between 6 and 24 was analyzed using diffusion Monte Carlo (DMC) calculations. The DMC results compare well with an analytic expression derived by one of us [Phys. Rev. B 78, 035111 (2008)] for the total energy to within 1% over an order of magnitude range of density, which increases with valley degeneracy. For Bi2Te3 (six-fold valley degeneracy) the applicable charge carrier densities are between 7*10^19cm^{-3} and 2*10^20cm^{-3}. DMC calculations distinguished between an exact and a useful approximate expression for the 24-fold degenerate MFEG polarizability for wave numbers 2p_F<q<7p_F. The analytical result for the MFEG is generalized to inhomogeneous systems by means of a gradient correction, the validity range of this approach is obtained. Employed within a density functional theory calculation this approximation compares well with DMC results for a quantum dot.
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.