Simone Chiesa

Finite-size error in many-body simulations with long-range interactions.

We discuss the origin of the finite-size error of the energy in many-body simulation of systems of charged particles and we propose a correction based on the random-phase approximation at long wavelengths. The correction is determined mainly by the collective charge oscillations of the interacting system. Finite-size corrections, both on kinetic and potential energy, can be calculated within a single simulation. Results are presented for the electron gas and silicon.

Quantum Monte Carlo study of the two-dimensional fermion Hubbard model

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 noninteracting 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 Green’s function in momentum space, allowing a more quantitative comparison with time-of-flight optical lattice experiments.

Accurate, Efficient, and Simple Forces Computed with Quantum Monte Carlo Methods

Computation of ionic forces using quantum Monte Carlo methods has long been a challenge. We introduce a simple procedure, based on known properties of physical electronic densities, to make the variance of the Hellmann-Feynman estimator finite. We obtain very accurate geometries for the molecules H(2), LiH, CH(4), NH(3), H(2)O, and HF, with a Slater-Jastrow trial wave function. Harmonic frequencies for diatomics are also in good agreement with experiment. An antithetical sampling method is also discussed for additional reduction of the variance.

Magnetism and pairing of two-dimensional trapped fermions.

The emergence of local phases in a trapped two-component Fermi gas in an optical lattice is studied using quantum Monte Carlo simulations. We treat temperatures that are comparable to or lower than those presently achievable in experiments and large enough systems that both magnetic and paired phases can be detected by inspection of the behavior of suitable short-range correlations. We use the latter to suggest the interaction strength and temperature range at which experimental observation of incipient magnetism and d-wave pairing are more likely and evaluate the relation between entropy and temperature in two-dimensional confined fermionic systems.

Disorder-induced stabilization of the pseudogap in strongly correlated systems.

The interplay of strong interaction and strong disorder, as contained in the Anderson-Hubbard model, is addressed using two nonperturbative numerical methods: the Lanczos algorithm in the grand canonical ensemble at zero temperature and quantum Monte Carlo simulations. We find distinctive evidence for a zero-energy anomaly which is robust upon variation of doping, disorder, and interaction strength. Its similarities to, and differences from, pseudogap formation in other contexts, including perturbative treatments of interactions and disorder, classical theories of localized charges, and in the clean Hubbard model, are discussed.

QMCPACK: An open source ab initio quantum Monte Carlo package for the electronic structure of atoms, molecules and solids

QMCPACK is an open source quantum Monte Carlo package for ab initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wavefunctions in conjunction with a sophisticated optimizer capable of optimizing tens of thousands of parameters. The orbital space auxiliary-field quantum Monte Carlo method is also implemented, enabling cross validation between different highly accurate methods. The code is specifically optimized for calculations with large numbers of electrons on the latest high performance computing architectures, including multicore central processing unit and graphical processing unit systems. We detail the programs capabilities, outline its structure, and give examples of its use in current research calculations. The package is available at http://qmcpack.org.

Phases of attractive spin-imbalanced fermions in square lattices

We determine the relative stability of different ground-state phases of spin-imbalanced popula- tions of attractive fermions in square lattices. The phases are systematically characterized by the symmetry of the order parameter and the real- and momentum-space structures using Hartree- Fock-Bogoliubov theory. We find several type of unidirectional Larkin-Ovchinikov-type phases. We discuss the effect of commensuration between the ordering wave vector and the density imbalance, and describe the mechanism of Fermi surface reconstruction and pairing for various orders. A robust supersolid phase is shown to exist when the ordering wave vector is diagonally directed.

Magnetic transition in a correlated band insulator

The effect of on-site electron-electron repulsion U in a band insulator is explored for a bilayer Hubbard Hamiltonian with opposite sign hopping in the two sheets. The ground-state phase diagram is determined at half-filling in the plane of U and the interplanar hybridization V through a computation of the antiferromagnetic (AF) structure factor, local moments, single-particle and spin wave spectra, and spin correlations. Unlike the case of the ionic Hubbard model, no evidence is found for a metallic phase intervening between the Mott and band insulators. Instead, upon increase of U at large V, the behavior of the local moments and of single-particle spectra give quantitative evidence of a crossover to a Mott insulator state preceding the onset of magnetic order. Our conclusions generalize those of single-site dynamical mean-field theory, and show that including interlayer correlations results in an increase of the single-particle gap with U.

Magnetic order in the Hubbard model in three dimensions and the crossover to two dimensions

Systems of fermions described by the three-dimensional (3D) repulsive Hubbard model on a cubic lattice have recently attracted considerable attention due to their possible experimental realization via cold atoms in an optical lattice. Analytical and numerical results are limited away from half-filling. We study the ground state of the doped system from weak to intermediate interaction strengths within the generalized Hartree-Fock approximation. The exact solution to the self-consistent-field equations in the thermodynamic limit is obtained and the ground state is shown to exhibit antiferromagnetic order and incommensurate spin-density waves (SDW). At low interaction strengths, the SDW has unidirectional character with a leading wavevector along the [100]-direction, and the system is metallic. As the interaction increases, the system undergoes a simultaneous structural and metal-to-insulator transition to a unidirectional SDW state along the [111]-direction, with a different wavelength. We systematically determine the real and momentum space properties of these states. The crossover from 3D to two dimensions (2D) is then studied by varying the inter-plane hopping amplitude, which can be experimentally realized by tuning the distance between a stack of square-lattice layers. Detailed comparisons are made between the exact numerical results and predictions from the pairing model, a variational ansatz based on the pairing of spins in the vicinity of the Fermi surface. Most of the numerical results can be understood quantitatively from the ansatz, which provides a simple picture for the nature of the SDW states.

Superconductivity and charge order of confined Fermi systems

The low-temperature properties of the two-dimensional attractive Hubbard model are strongly influenced by the fermion density. Away from half-filling, there is a finite-temperature transition to a phase with s-wave pairing order. However, the critical temperature is suppressed to zero at half-filling, where long-range charge-density-wave order also appears, degenerate with superconductivity. This paper presents Determinant Quantum Monte Carlo simulations of the attractive Hubbard model in the presence of a confining potential V which makes the fermion density \rho{} inhomogeneous across the lattice. Pair correlations are shown to be large at low temperatures in regions of the trapped system with incommensurate filling, and to exhibit a minimum as the local density \rho(i) passes through one fermion per site. In this ring of \rho=1, charge order is enhanced. A comparison is made between treating V within the local-density approximation (LDA) and in an ab initio manner. It is argued that certain sharp features of the LDA result at integer filling do not survive the proximity of doped sites. The critical temperature of confined systems of fixed characteristic density is estimated.