1
|
Quiton SJ, Wu H, Xing X, Lin L, Head-Gordon M. The Staggered Mesh Method: Accurate Exact Exchange Toward the Thermodynamic Limit for Solids. J Chem Theory Comput 2024. [PMID: 39213528 DOI: 10.1021/acs.jctc.4c00771] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
In periodic systems, the Hartree-Fock (HF) exchange energy exhibits the slowest convergence of all HF energy components as the system size approaches the thermodynamic limit. We demonstrate that the recently proposed staggered mesh method for Fock exchange energy [Xing, Li, and Lin, Math. Comp., 2024], which is specifically designed to sidestep certain singularities in exchange energy evaluation, can expedite the finite-size convergence rate for the exact exchange energy across a range of insulators and semiconductors when compared to the regular and truncated Coulomb methods. This remains true even for two computationally cheaper versions of this new method, which we call non-SCF and split-SCF staggered mesh. Additionally, a sequence of numerical tests on simple solids showcases the staggered mesh method's ability to improve convergence toward the thermodynamic limit for band gaps, bulk moduli, equilibrium lattice dimensions, energies, and phonon force constants.
Collapse
Affiliation(s)
- Stephen Jon Quiton
- College of Chemistry, University of California, Berkeley, California 94720, United States
| | - Hamlin Wu
- College of Chemistry, University of California, Berkeley, California 94720, United States
| | - Xin Xing
- Department of Mathematics, University of California, Berkeley, California 94720, United States
| | - Lin Lin
- Department of Mathematics, University of California, Berkeley, California 94720, United States
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California, 94720, United States
| | - Martin Head-Gordon
- College of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
2
|
Harsha G, Abraham V, Zgid D. Challenges with relativistic GW calculations in solids and molecules. Faraday Discuss 2024. [PMID: 39101408 DOI: 10.1039/d4fd00043a] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/06/2024]
Abstract
For molecules and solids containing heavy elements, accurate electronic-structure calculations require accounting not only for electronic correlations but also for relativistic effects. In molecules, relativity can lead to severe changes in the ground-state description. In solids, the interplay between both correlation and relativity can change the stability of phases or it can lead to an emergence of completely new phases. Traditionally, the simplest illustration of relativistic effects can be done either by including pseudopotentials in non-relativistic calculations or alternatively by employing large all-electron basis sets in relativistic methods. By analyzing different electronic properties (band structure, equilibrium lattice constant and bulk modulus) in semiconductors and insulators, we show that capturing the interplay of relativity and electron correlation can be rather challenging in Green's function methods. For molecular problems with heavy elements, we also observe that similar problems persist. We trace these challenges to three major problems: deficiencies in pseudopotential treatment as applied to Green's function methods, the scarcity of accurate and compact all-electron basis sets that can be converged with respect to the basis-set size, and linear dependencies arising in all-electron basis sets, particularly when employing Gaussian orbitals. Our analysis provides detailed insight into these problems and opens a discussion about potential approaches to mitigate them.
Collapse
Affiliation(s)
- Gaurav Harsha
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA.
| | - Vibin Abraham
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA.
| | - Dominika Zgid
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA.
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
3
|
Liang YH, Ye HZ, Berkelbach TC. Can Spin-Component Scaled MP2 Achieve kJ/mol Accuracy for Cohesive Energies of Molecular Crystals? J Phys Chem Lett 2023; 14:10435-10441. [PMID: 37956873 DOI: 10.1021/acs.jpclett.3c02411] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2023]
Abstract
Attaining kJ/mol accuracy in cohesive energy for molecular crystals is a persistent challenge in computational materials science. In this study, we evaluate second-order Møller-Plesset perturbation theory (MP2) and its spin-component scaled models for calculating cohesive energies for 23 molecular crystals (X23 data set). Using periodic boundary conditions and Brillouin zone sampling, we converge results to the thermodynamic and complete basis set limits, achieving an accuracy of about 2 kJ/mol (0.5 kcal/mol), which is rarely achieved in previous MP2 calculations for molecular crystals. When compared to experimental data, our results have a mean absolute error of 12.9 kJ/mol, comparable to Density Functional Theory with the PBE functional and TS dispersion correction. By separately scaling the opposite-spin and same-spin correlation energy components, using predetermined parameters, we reduce the mean absolute error to 9.5 kJ/mol. Further fine-tuning of these scaling parameters specifically for the X23 data set brings the mean absolute error down to 7.5 kJ/mol.
Collapse
Affiliation(s)
- Yu Hsuan Liang
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| |
Collapse
|
4
|
Rettig A, Lee J, Head-Gordon M. Even Faster Exact Exchange for Solids via Tensor Hypercontraction. J Chem Theory Comput 2023; 19:5773-5784. [PMID: 37586065 DOI: 10.1021/acs.jctc.3c00407] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/18/2023]
Abstract
Hybrid density functional theory (DFT) remains intractable for large periodic systems due to the demanding computational cost of exact exchange. We apply the tensor hypercontraction (THC) (or interpolative separable density fitting) approximation to periodic hybrid DFT calculations with Gaussian-type orbitals using the Gaussian plane wave approach. This is done to lower the computational scaling with respect to the number of basis functions (N) and k-points (Nk) at a fixed system size. Additionally, we propose an algorithm to fit only occupied orbital products via THC (i.e., a set of points, NISDF) to further reduce computation time and memory usage. This algorithm has linear scaling cost with k-points, no explicit dependence of NISDF on basis set size, and overall cubic scaling with unit cell size. Significant speedups and reduced memory usage may be obtained for moderately sized k-point meshes, with additional gains for large k-point meshes. Adequate accuracy can be obtained using THC-oo-K for self-consistent calculations. We perform illustrative hybrid density function theory calculations on the benzene crystal in the basis set and thermodynamic limits to highlight the utility of this algorithm.
Collapse
Affiliation(s)
- Adam Rettig
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
5
|
Gygi F. All-Electron Plane-Wave Electronic Structure Calculations. J Chem Theory Comput 2023; 19:1300-1309. [PMID: 36757291 PMCID: PMC9979607 DOI: 10.1021/acs.jctc.2c01191] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
Abstract
We demonstrate the use of the plane wave basis for all-electron electronic structure calculations. The approach relies on the definition of an analytic, norm-conserving, regularized Coulomb potential, and a scalable implementation of the plane wave method capable of handling large energy cutoffs (up to 80 kRy in the examples shown). The method is applied to the computation of electronic properties of isolated atoms as well as the diamond and silicon crystals, MgO, solid argon, and a configuration of 64 water molecules extracted from a first-principles molecular dynamics simulation. The computed energies, band gaps, ionic forces, and stress tensors provide reference results for the validation of pseudopotentials and/or localized basis sets. A calculation of the all-electron band structure of diamond and silicon using the SCAN meta-GGA density functional allows for a validation of calculations based on pseudopotentials derived using the PBE exchange-correlation functional. In the case of (H2O)64, the computed ionic forces provide a reference from which the errors incurred in pseudopotential calculations and in localized Gaussian basis sets calculations can be estimated.
Collapse
|
6
|
Zhang C, Verma P, Wang J, Liu Y, He X, Wang Y, Truhlar DG, Liu Z. Performance of Screened-Exchange Functionals for Band Gaps and Lattice Constants of Crystals. J Chem Theory Comput 2023; 19:311-323. [PMID: 36520598 DOI: 10.1021/acs.jctc.2c00822] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Kohn-Sham density functional theory is the most widely used method for electronic structure calculations of solid-state systems. The screened-exchange functionals developed following the influential work of Scuseria and co-workers in 2003-2006 have significantly improved the accuracy of the predictions of solid-state properties. This work assesses six screened-exchange density functionals for the prediction of 60 band gaps (database BG60) and 68 lattice constants (database LC68). The band gaps are calculated with both consistently calculated lattice constants and experimental lattice constants. Results for the nonlocal screened-exchange functionals are compared with those for six widely used or recently developed local functionals. The results show that all the screened-exchange functionals have smaller mean absolute errors (MAEs) than any of the local functionals. All the functionals except HLE17 overestimate (on average) the lattice constants, and M06-SX gives the best performance among the compared functionals, with a MAE of 0.051 Å. All the functionals underestimate (on average) the band gaps, and M06-SX outperforms all other functionals, with a MAE of 0.47 eV. M06-SX also has the lowest root-mean-squared error for both LC68 and BG60. For the subdatabases of BG60, M06-SX shows better performance for ionic crystals and systems with large band gaps, while HSE12s gives better results for semiconductors and systems with small band gaps. Overall, M06-SX shows the best performance for solid-state systems, followed by N12-SX and HSE12s. The best-performing local functionals are M06-L, revM06-L, and HLE17 for band gaps and M06-L and revM06-L for lattice constants. We found that M06-SX, revM06-L, and N12-SX not only are well optimized for a broad array of chemical properties but also have very good performance for the databases in this paper, making them well-suited for applications involving heterogeneous chemistry.
Collapse
Affiliation(s)
- Cheng Zhang
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha, Hunan 410006, China
| | - Pragya Verma
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Jiaxu Wang
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha, Hunan 410006, China
| | - Yiwei Liu
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.,New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200062, China
| | - Ying Wang
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha, Hunan 410006, China
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Zhonghua Liu
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha, Hunan 410006, China
| |
Collapse
|
7
|
Lee J, Rettig A, Feng X, Epifanovsky E, Head-Gordon M. Faster Exact Exchange for Solids via occ-RI-K: Application to Combinatorially Optimized Range-Separated Hybrid Functionals for Simple Solids with Pseudopotentials Near the Basis Set Limit. J Chem Theory Comput 2022; 18:7336-7349. [PMID: 36459992 PMCID: PMC10441520 DOI: 10.1021/acs.jctc.2c00742] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/04/2022]
Abstract
In this work, we developed and showcased the occ-RI-K algorithm to compute the exact exchange contribution in density functional calculations of solids near the basis set limit. Within the Gaussian planewave (GPW) density fitting, our algorithm achieves a 1-2 orders of magnitude speedup compared to conventional GPW algorithms. Since our algorithm is well suited for simulations with large basis sets, we applied it to 12 hybrid density functionals with pseudopotentials and a large uncontracted basis set to assess their performance on band gaps of 25 simple solids near the basis set limit. The largest calculation performed in this work involves 16 electrons and 350 basis functions in the unit cell utilizing a 6 × 6 × 6 k-mesh. With 20-27% exact exchange, global hybrid functionals (B3LYP, PBE0, revPBE0, B97-3, SCAN0) perform similarly with a root-mean-square deviation (RMSD) of 0.61-0.77 eV, while other global hybrid functionals such as M06-2X (2.02 eV) and MN15 (1.05 eV) show higher RMSD due to their increased fraction of exact exchange. A short-range hybrid functional, HSE achieves a similar RMSD (0.76 eV) but shows a notable underestimation of band gaps due to the complete lack of long-range exchange. We found that two combinatorially optimized range-separated hybrid functionals, ωB97X-rV (3.94 eV) and ωB97M-rV (3.40 eV), and the two other range-separated hybrid functionals, CAM-B3LYP (2.41 eV) and CAM-QTP01 (4.16 eV), significantly overestimate the band gap because of their high fraction of long-range exact exchange. Given the failure of ωB97X-rV and ωB97M-rV, we have yet to find a density functional that offers consistent performance for both molecules and solids. Our algorithm development and density functional assessment will serve as a stepping stone toward developing more accurate hybrid functionals and applying them to practical applications.
Collapse
Affiliation(s)
- Joonho Lee
- Department of Chemistry, Columbia University, New York, NY, USA
| | - Adam Rettig
- Department of Chemistry, University of California, Berkeley, CA, USA
| | | | | | | |
Collapse
|
8
|
Goldzak T, Wang X, Ye HZ, Berkelbach TC. Accurate thermochemistry of covalent and ionic solids from spin-component-scaled MP2. J Chem Phys 2022; 157:174112. [PMID: 36347707 PMCID: PMC9637026 DOI: 10.1063/5.0119633] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2022] [Accepted: 10/11/2022] [Indexed: 11/14/2022] Open
Abstract
We study the performance of spin-component-scaled second-order Møller-Plesset perturbation theory (SCS-MP2) for the prediction of the lattice constant, bulk modulus, and cohesive energy of 12 simple, three-dimensional covalent and ionic semiconductors and insulators. We find that SCS-MP2 and the simpler scaled opposite-spin MP2 (SOS-MP2) yield predictions that are significantly improved over the already good performance of MP2. Specifically, when compared to experimental values with zero-point vibrational corrections, SCS-MP2 (SOS-MP2) yields mean absolute errors of 0.015 (0.017) Å for the lattice constant, 3.8 (3.7) GPa for the bulk modulus, and 0.06 (0.08) eV for the cohesive energy, which are smaller than those of leading density functionals by about a factor of two or more. We consider a reparameterization of the spin-scaling parameters and find that the optimal parameters for these solids are very similar to those already in common use in molecular quantum chemistry, suggesting good transferability and reliable future applications to surface chemistry on insulators.
Collapse
Affiliation(s)
- Tamar Goldzak
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Xiao Wang
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, USA
| | - Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | |
Collapse
|
9
|
Bintrim SJ, Berkelbach TC, Ye HZ. Integral-Direct Hartree-Fock and Møller-Plesset Perturbation Theory for Periodic Systems with Density Fitting: Application to the Benzene Crystal. J Chem Theory Comput 2022; 18:5374-5381. [PMID: 35969856 DOI: 10.1021/acs.jctc.2c00640] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present an algorithm and implementation of integral-direct, density-fitted Hartree-Fock (HF) and second-order Møller-Plesset perturbation theory (MP2) for periodic systems. The new code eliminates the formerly prohibitive storage requirements and allows us to study systems 1 order of magnitude larger than before at the periodic MP2 level. We demonstrate the significance of the development by studying the benzene crystal in both the thermodynamic limit and the complete basis set limit, for which we predict an MP2 cohesive energy of -72.8 kJ/mol, which is about 10-15 kJ/mol larger in magnitude than all previously reported MP2 calculations. Compared to the best theoretical estimate from literature, several modified MP2 models approach chemical accuracy in the predicted cohesive energy of the benzene crystal and hence may be promising cost-effective choices for future applications on molecular crystals.
Collapse
Affiliation(s)
- Sylvia J Bintrim
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, New York 10027, United States.,Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| | - Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| |
Collapse
|
10
|
Zen A, Grüneis A, Alfè D, Rossi M. Beyond GGA total energies for solids and surfaces. J Chem Phys 2022; 157:050401. [PMID: 35933211 DOI: 10.1063/5.0107716] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andrea Zen
- Dipartimento di Fisica Ettore Pancini, Università di Napoli Federico II, Monte S. Angelo, I-80126 Napoli, Italy
| | - Andreas Grüneis
- Institute for Theoretical Physics, Vienna University of Technology (TU Wien), Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria
| | - Dario Alfè
- Dipartimento di Fisica Ettore Pancini, Università di Napoli Federico II, Monte S. Angelo, I-80126 Napoli, Italy
| | - Mariana Rossi
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|
11
|
Villot C, Ballesteros F, Wang D, Lao KU. Coupled Cluster Benchmarking of Large Noncovalent Complexes in L7 and S12L as Well as the C 60 Dimer, DNA-Ellipticine, and HIV-Indinavir. J Phys Chem A 2022; 126:4326-4341. [PMID: 35766331 DOI: 10.1021/acs.jpca.2c01421] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this work, we report the benchmark binding energies of the seven complexes within the L7 data set, six host-guest complexes from the S12L data set, a C60 dimer, the DNA-ellipticine intercalation complex, and the largest system of the study, the HIV-indinavir system, which contained 343 atoms or 139 heavy atoms. The high-quality values reported were obtained via a focal point method that relies on the canonical form of second-order Møller-Plesset theory and the domain-based local pair natural orbital scheme for the coupled cluster with single double and perturbative triple excitations [DLPNO-CCSD(T)] extrapolated to the complete basis set (CBS) limit. The results in this work not only corroborate but also improve upon some previous benchmark values for large noncovalent complexes albeit at a relatively steep cost. Although local CCSD(T) and the largely successful fixed-node diffusion Monte Carlo (FN-DMC) have been shown to generally agree for small- to medium-size systems, a discrepancy in their reported binding energy values arises for large complexes, where the magnitude of the disagreement is a definite cause for concern. For example, the largest deviation in the L7 data set was 2.8 kcal/mol (∼10%) on the low end in C3GC. Such a deviation only grows worse in the S12L set, which showed a difference of up to 10.4 kcal/mol (∼25%) by a conservative estimation in buckycatcher-C60. The DNA-ellipticine complex also generated a disagreement of 4.4 kcal/mol (∼10%) between both state-of-the-art methods. The disagreement between local CCSD(T) and FN-DMC in large noncovalent complexes shows that it is urgently needed to have the canonical CCSD(T), the Monte Carlo CCSD(T), or the full configuration interaction quantum Monte Carlo approaches available to large systems on the hundred-atom scale to solve this dilemma. In addition, the performances of cheaper popular computational methods were assessed for the studied complexes with respect to DLPNO-CCSD(T)/CBS. r2SCAN-3c, B97M-V, and PBE0+D4 work well in large noncovalent complexes in this work, and GFN2-xTB performs well in π-π stacking complexes. B97M-V is the most reliable computationally efficient approach to predicting noncovalent interactions for large complexes, being the only one to have binding errors within the so-called 1 kcal/mol "chemical accuracy". The benchmark interaction energies of these host-guest complexes, molecular materials, and biological systems with electronic and medicinal implications provide crucial reference data for the improvement of current and future lower-cost methods.
Collapse
Affiliation(s)
- Corentin Villot
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Francisco Ballesteros
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Danyang Wang
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Ka Un Lao
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| |
Collapse
|
12
|
Laun J, Bredow T. BSSE-corrected consistent Gaussian basis sets of triple-zeta valence with polarization quality of the fifth period for solid-state calculations. J Comput Chem 2022; 43:839-846. [PMID: 35302265 DOI: 10.1002/jcc.26839] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2022] [Revised: 03/03/2022] [Accepted: 03/04/2022] [Indexed: 01/28/2023]
Abstract
Revised versions of our published pob-TZVP basis sets [Laun, J.; Vilela Oliveira, D. and Bredow, T., J. Comput. Chem., 2018, 39 (19), 1285-1290] have been derived for periodic quantum-chemical solid-state calculations. They complete our pob-TZVP-rev2 series [Vilela Oliveira, D.; Laun, J.; Peintinger, M. F. and Bredow, T., J. Comput. Chem., 2019, 40 (27), 2364-2376 and Laun, J. and Bredow, J. Comput. Chem. 2021; 42 (15), 1064-1072] for the elements of the fifth period and are based on the fully relativistic effective core potentials (ECPs) of the Stuttgart/Cologne group and the def2-TZVP valence basis of the Ahlrichs group. The pob-TZVP-rev2 basis sets are developed to minimize the basis set superposition error (BSSE) in crystalline systems. For the applied PW1PW hybrid functional, the overall performance, transferability, and SCF stability of the resulting pob-TZVP-rev2 basis sets are significantly improved compared to the original pob-TZVP basis sets. After augmentation with single diffuse s- and p-functions, reference plane-wave band structures of metals can be accurately reproduced.
Collapse
Affiliation(s)
- Joachim Laun
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, University of Bonn, Bonn, Germany
| | - Thomas Bredow
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, University of Bonn, Bonn, Germany
| |
Collapse
|
13
|
Ye HZ, Berkelbach TC. Correlation-Consistent Gaussian Basis Sets for Solids Made Simple. J Chem Theory Comput 2022; 18:1595-1606. [PMID: 35192359 DOI: 10.1021/acs.jctc.1c01245] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The rapidly growing interest in simulating condensed-phase materials using quantum chemistry methods calls for a library of high-quality Gaussian basis sets suitable for periodic calculations. Unfortunately, most standard Gaussian basis sets commonly used in molecular simulation show significant linear dependencies when used in close-packed solids, leading to severe numerical issues that hamper the convergence to the complete basis set (CBS) limit, especially in correlated calculations. In this work, we revisit Dunning's strategy for construction of correlation-consistent basis sets and examine the relationship between accuracy and numerical stability in periodic settings. We find that limiting the number of primitive functions avoids the appearance of problematic small exponents while still providing smooth convergence to the CBS limit. As an example, we generate double-, triple-, and quadruple-ζ correlation-consistent Gaussian basis sets for periodic calculations with Goedecker-Teter-Hutter (GTH) pseudopotentials. Our basis sets cover the main-group elements from the first three rows of the periodic table. Especially for atoms on the left side of the periodic table, our basis sets are less diffuse than those used in molecular calculations. We verify the fast and reliable convergence to the CBS limit in both Hartree-Fock and post-Hartree-Fock (MP2) calculations, using a diverse test set of 19 semiconductors and insulators.
Collapse
Affiliation(s)
- Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, New York 10027, United States.,Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| |
Collapse
|