A form of Jastrow factor is introduced for use in quantum Monte Carlo simulations of finite and periodic systems. Test data are presented for atoms, molecules, and solids, including both all-electron and pseudopotential atoms. We demonstrate that our Jastrow factor is able to retrieve a large fraction of the correlation energy.
We propose a Jastrow factor for electron-electron correlations that interpolates between the radial symmetry of the Coulomb interaction at short inter-particle distance and the space-group symmetry of the simulation cell at large separation. The proposed Jastrow factor captures comparable levels of the correlation energy to current formalisms, is 40% quicker to evaluate, and offers benefits in ease of use, as we demonstrate in quantum Monte Carlo simulations.
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 wave functions 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 (CPU) and graphical processing unit (GPU) 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://www.qmcpack.org .
One of the standard generalized-gradient approximations (GGAs) in use in modern electronic-structure theory, PBE, and a recently proposed modification designed specifically for solids, PBEsol, are identified as particular members of a family of functionals taking their parameters from different properties of homogeneous or inhomogeneous electron liquids. Three further members of this family are constructed and tested, together with the original PBE and PBEsol, for atoms, molecules and solids. We find that PBE, in spite of its popularity in solid-state physics and quantum chemistry, is not always the best performing member of the family, and that PBEsol, in spite of having been constructed specifically for solids, is not the best for solids. The performance of GGAs for finite systems is found to sensitively depend on the choice of constraints steaming from infinite systems. Guidelines both for users and for developers of density functionals emerge from this work.
We demonstrate a trap that confines polarizable particles around the antinode of a standing-wave microwave field. The trap relies only on the polarizability of the particles far from any resonances, so can trap a wide variety of atoms and molecules in a wide range of internal states, including the ground state. The trap has a volume of about 10 cm$^3$, and a depth approaching 1 K for many polar molecules. We measure the trap properties using $^{7}$Li atoms, showing that when the input microwave power is 610 W, the atoms remain trapped with a $1/e$ lifetime of 1.76(12) s, oscillating with an axial frequency of 28.55(5) Hz and a radial frequency of 8.81(8) Hz. The trap could be loaded with slow molecules from a range of available sources, and is particularly well suited to sympathetic cooling and evaporative cooling of molecules.
The calculation of the MP2 correlation energy for extended systems can be viewed as a multi-dimensional integral in the thermodynamic limit, and the standard method for evaluating the MP2 energy can be viewed as a trapezoidal quadrature scheme. We demonstrate that existing analysis neglects certain contributions due to the non-smoothness of the integrand, and may significantly underestimate finite-size errors. We propose a new staggered mesh method, which uses two staggered Monkhorst-Pack meshes for occupied and virtual orbitals, respectively, to compute the MP2 energy. The staggered mesh method circumvents a significant error source in the standard method, in which certain quadrature nodes are always placed on points where the integrand is discontinuous. One significant advantage of the proposed method is that there are no tunable parameters, and the additional numerical effort needed can be negligible compared to the standard MP2 calculation. Numerical results indicate that the staggered mesh method can be particularly advantageous for quasi-1D systems, as well as quasi-2D and 3D systems with certain symmetries.