Ab initio computations of molecular systems by the auxiliary-field quantum Monte Carlo method

The auxiliary-field quantum Monte Carlo (AFQMC) method provides a computational framework for solving the time-independent Schroedinger equation in atoms, molecules, solids, and a variety of model systems. AFQMC has recently witnessed remarkable growth, especially as a tool for electronic structure computations in real materials. The method has demonstrated excellent accuracy across a variety of correlated electron systems. Taking the form of stochastic evolution in a manifold of non-orthogonal Slater determinants, the method resembles an ensemble of density-functional theory (DFT) calculations in the presence of fluctuating external potentials. Its computational cost scales as a low-power of system size, similar to the corresponding independent-electron calculations. Highly efficient and intrinsically parallel, AFQMC is able to take full advantage of contemporary high-performance computing platforms and numerical libraries. In this review, we provide a self-contained introduction to the exact and constrained variants of AFQMC, with emphasis on its applications to the electronic structure in molecular systems. Representative results are presented, and theoretical foundations and implementation details of the method are discussed.Comment: 22 pages, 11 figure

Calculation of interatomic forces and optimization of molecular geometry with auxiliary-field quantum Monte Carlo

We propose an algorithm for accurate, systematic and scalable computation of interatomic forces within the auxiliary-field Quantum Monte Carlo (AFQMC) method. The algorithm relies on the Hellman-Fenyman theorem, and incorporates Pulay corrections in the presence of atomic orbital basis sets. We benchmark the method for small molecules by comparing the computed forces with the derivatives of the AFQMC potential energy surface, and by direct comparison with other quantum chemistry methods. We then perform geometry optimizations using the steepest descent algorithm in larger molecules. With realistic basis sets, we obtain equilibrium geometries in agreement, within statistical error bars, with experimental values. The increase in computational cost for computing forces in this approach is only a small prefactor over that of calculating the total energy. This paves the way for a general and efficient approach for geometry optimization and molecular dynamics within AFQMC.Comment: 5 pages, 4 figure

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.Comment: 6.5 pages, 5 figure

Ultracold atoms in a square lattice with spin-orbit coupling: Charge order, superfluidity, and topological signatures

We present an $\textit{ab initio}$, numerically exact study of attractive fermions in square lattices with Rashba spin-orbit coupling. The ground state of this system is a supersolid, with co-existing charge and superfluid order. The superfluid is composed of both singlet and triplet pairs induced by spin-orbit coupling. We perform large-scale calculations using auxiliary-field quantum Monte Carlo to provide the first full, quantitative description of the charge, spin, and pairing properties of the system. In addition to characterizing the exotic physics, our results will serve as essential high-accuracy benchmarks for the intense theoretical and especially experimental efforts in ultracold atoms to realize and understand an expanding variety of quantum Hall and topological superconductor systems.Comment: 6 pages, 4 figure

Finite-Temperature Monte Carlo Calculations For Systems With Fermions

We present a quantum Monte Carlo method which allows calculations on many-fermion systems at finite temperatures without any sign decay. This enables simulations of the grand-canonical ensemble at large system sizes and low temperatures. Both diagonal and off-diagonal expectations can be computed straightforwardly. The sign decay is eliminated by a constraint on the fermion determinant. The algorithm is approximate. Tests on the Hubbard model show that accurate results on the energy and correlation functions can be obtained.Comment: 5 pages, RevTex; to appear in Phys. Rev. Let