Accelerating the Convergence of Auxiliary-Field Quantum Monte Carlo inSolids with Optimized Gaussian Basis Sets

Abstract

We investigate the use of optimized correlation consistent gaussian basis sets for the study of insulating solids with auxiliary-field quantum Monte Carlo (AFQMC). The exponents of the basis set are optimized through the minimization of the second order Møller--Plesset perturbation theory (MP2) energy in a small unit cell of the solid. We compare against other alternative basis sets proposed in the literature, namely calculations in the Kohn--Sham basis and in the natural orbitals of an MP2 calculation. We find that our optimized basis sets accelerate the convergence of the AFQMC correlation energy compared to a Kohn--Sham basis, and offer similar convergence to MP2 natural orbitals at a fraction of the cost needed to generate them. We also suggest the use of an improved, method independent, MP2-based basis set correction that significantly reduces the required basis set sizes needed to converge the correlation energy. With these developments, we study the relative performance of these basis sets in LiH, Si and MgO, and determine that our optimized basis sets yield the most consistent results as a function of volume. Using these optimized basis sets, we systematically converge the AFQMC calculations to the complete basis set and thermodynamic limit and find excellent agreement with experiment for systems studied. Although we focus on AFQMC, our basis set generation procedure is independent of the subsequent correlated wavefunction method used.

arXiv commentAdd up to date supplementary material

References in corpus (28)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Optimized norm-conserving Vanderbilt pseudopotentials
• Optimization Algorithm for the Generation of ONCV Pseudopotentials
• Recent developments in the PySCF program package
• Random-phase approximation and its applications in computational chemistry and materials science
• A Constrained Path Monte Carlo Method for Fermion Ground States
• Quantum Monte Carlo method using phase-free random walks with Slater determinants
• Daubechies wavelets as a basis set for density functional pseudopotential calculations
• QMCPACK : An open source ab initio Quantum Monte Carlo package for the electronic structure of atoms, molecules, and solids
• Assessment of correlation energies based on the random-phase approximation
• QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion Quantum Monte Carlo
• Convergence of many-body wavefunction expansions using a plane wave basis: from the homogeneous electron gas to the solid state
• Applying the Coupled-Cluster Ansatz to Solids and Surfaces in the Thermodynamic Limit
• Phaseless auxiliary-field quantum Monte Carlo calculations with planewaves and pseudopotentials--applications to atoms and molecules
• Explicitly correlated plane waves: Accelerating convergence in periodic wavefunction expansions
• Pressure-induced diamond to beta-tin transition in bulk silicon: a near-exact quantum Monte Carlo study
• From plane waves to local Gaussians for the simulation of correlated periodic systems
• Quantum Monte Carlo calculations in solids with downfolded Hamiltonians
• Overcoming the Memory Bottleneck in Auxiliary Field Quantum Monte Carlo Simulations with Interpolative Separable Density Fitting
• Excited state calculations in solids by auxiliary-field quantum Monte Carlo
• Accelerating Auxiliary-Field Quantum Monte Carlo Simulations of Solids with Graphical Processing Unit
• Frozen-orbital and downfolding calculations with auxiliary-field quantum Monte Carlo
• Auxiliary-field quantum Monte Carlo calculations of the structural properties of nickel oxide
• On the duality of ring and ladder diagrams and its importance for many-electron perturbation theories
• Auxiliary-field quantum Monte Carlo calculations with multiple-projector pseudopotentials
• RPA natural orbitals and their application to post-Hartree-Fock electronic structure methods
• An optimized twist angle to find the twist-averaged correlation energy applied to the uniform electron gas
• Particle-particle ladder based basis-set corrections applied to atoms and molecules using coupled-cluster theory

Cited by in corpus (6)

• Correlation-consistent Gaussian basis sets for solids made simple
• Fault-tolerant quantum simulation of materials using Bloch orbitals
• Low-Scaling Algorithm for the Random Phase Approximation using Tensor Hypercontraction with k-point Sampling
• Focal-point approach with pair-specific cusp correction for coupled-cluster theory
• Toward an accurate equation of state and B1-B2 phase boundary for magnesium oxide to TPa pressures and eV temperatures
• Approaching the Basis Set Limit in Gaussian-Orbital-Based Periodic Calculations with Transferability: Performance of Pure Density Functionals for Simple Semiconductors