Approaching the Basis Set Limit in Gaussian-Orbital-Based Periodic Calculations with Transferability: Performance of Pure Density Functionals for Simple Semiconductors

Abstract

Simulating solids with quantum chemistry methods and Gaussian-type orbitals (GTOs) has been gaining popularity. Nonetheless, there are few systematic studies that assess the basis set incompleteness error (BSIE) in these GTO-based simulations over a variety of solids. In this work, we report a GTO-based implementation for solids, and apply it to address the basis set convergence issue. We employ a simple strategy to generate large uncontracted (unc) GTO basis sets, that we call the unc-def2-GTH sets. These basis sets exhibit systematic improvement towards the basis set limit as well as good transferability based on application to a total of 43 simple semiconductors. Most notably, we found the BSIE of unc-def2-QZVP-GTH to be smaller than 0.7 m$E_h$ per atom in total energies and 20 meV in band gaps for all systems considered here. Using unc-def2-QZVP-GTH, we report band gap benchmarks of a combinatorially designed meta generalized gradient functional (mGGA), B97M-rV, and show that B97M-rV performs similarly (a root-mean-square-deviation (RMSD) of 1.18 eV) to other modern mGGA functionals, M06-L (1.26 eV), MN15-L (1.29 eV), and SCAN (1.20 eV). This represents a clear improvement over older pure functionals such as LDA (1.71 eV) and PBE (1.49 eV) though all these mGGAs are still far from being quantitatively accurate. We also provide several cautionary notes on the use of our uncontracted bases and on future research on GTO basis set development for solids.

References in corpus (22)

• Restoring the density-gradient expansion for exchange in solids and surfaces
• Generalized gradient approximation for solids and their surfaces
• The SIESTA method for ab initio order-N materials simulation
• Efficient index handling of multidimensional periodic boundary conditions
• Separable Dual Space Gaussian Pseudo-potentials
• Relativistic separable dual-space Gaussian Pseudopotentials from H to Rn
• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• Quantum ESPRESSO toward the exascale
• The GW method
• Efficient implementation of a van der Waals density functional: Application to double-wall carbon nanotubes
• Nonlocal van der Waals density functional: The simpler the better
• Recent developments in the PySCF program package
• Understanding Band Gaps of Solids in Generalized Kohn-Sham Theory
• Predicting Band Gaps with Hybrid Density Functionals
• An exact Coulomb cutoff technique for supercell calculations
• More-Realistic Band Gaps from Meta-Generalized Gradient Approximations: Only in a Generalized Kohn-Sham Scheme
• Comparative study of hybrid functionals applied to structural and electronic properties of semiconductors and insulators
• A DFT-based Molecular Transport Implementation in ADF/BAND
• Accurate Hartree-Fock energy of extended systems using large Gaussian basis sets
• RPA natural orbitals and their application to post-Hartree-Fock electronic structure methods
• Accelerating the Convergence of Auxiliary-Field Quantum Monte Carlo inSolids with Optimized Gaussian Basis Sets
• LiH as a Li${}^{+}$ and H${}^{-}$ ion provider

Cited by in corpus (3)

• Roadmap on Electronic Structure Codes in the Exascale Era
• Correlation-consistent Gaussian basis sets for solids made simple
• Integral-direct Hartree-Fock and Møller-Plesset Perturbation Theory for Periodic Systems with Density Fitting: Application to the Benzene Crystal