skip to main content


Title: Explaining and Fixing DFT Failures for Torsional Barrier
Most torsional barriers are predicted with high accuracies (about 1 kJ/mol) by standard semilocal functionals, but a small subset was found to have much larger errors. We created a database of almost 300 carbon–carbon torsional barriers, including 12 poorly behaved barriers, that stem from the Y═C—X group, where Y is O or S and X is a halide. Functionals with enhanced exchange mixing (about 50%) worked well for all barriers. We found that poor actors have delocalization errors caused by hyperconjugation. These problematic calculations are density-sensitive (i.e., DFT predictions change noticeably with the density), and using HF densities (HF-DFT) fixes these issues. For example, conventional B3LYP performs as accurately as exchange-enhanced functionals if the HF density is used. For long-chain conjugated molecules, HF-DFT can be much better than exchange-enhanced functionals. We suggest that HF-PBE0 has the best overall performance.  more » « less
Award ID(s):
1856165
PAR ID:
10220443
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
The journal of physical chemistry letters
Issue:
12
ISSN:
1948-7185
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Empirical fitting of parameters in approximate density functionals is common. Such fits conflate errors in the self-consistent density with errors in the energy functional, but density-corrected DFT (DC-DFT) separates these two. We illustrate with catastrophic failures of a toy functional applied to H2+ at varying bond lengths, where the standard fitting procedure misses the exact functional; Grimme’s D3 fit to noncovalent interactions, which can be contaminated by large density errors such as in the WATER27 and B30 data sets; and double-hybrids trained on self-consistent densities, which can perform poorly on systems with density-driven errors. In these cases, more accurate results are found at no additional cost by using Hartree–Fock (HF) densities instead of self-consistent densities. For binding energies of small water clusters, errors are greatly reduced. Range-separated hybrids with 100% HF at large distances suffer much less from this effect. 
    more » « less
  2. Kohn–Sham (KS) inversion, that is, the finding of the exact KS potential for a given density, is difficult in localized basis sets. We study the precision and reliability of several inversion schemes, finding estimates of density-driven errors at a useful level of accuracy. In typical cases of substantial density-driven errors, Hartree–Fock density functional theory (HF-DFT) is almost as accurate as DFT evaluated on CCSD(T) densities. A simple approximation in practical HF-DFT also makes errors much smaller than the density-driven errors being calculated. Two paradigm examples, stretched NaCl and the HO·Cl– radical, illustrate just how accurate HF-DFT is. 
    more » « less
  3. Abstract

    The hafnium‐rich portion of the of the hafnium‐nitrogen phase diagram is dominated by a substoichiometric rocksalt HfN1‐x, the ζ‐Hf4N3−x, the η‐Hf3N2−x, and the elemental Hf phase. The zeta and eta nitride phases have a close packed metal atom stacking sequence but their nitrogen atom ordering has yet to be concretely identified. With respect to the composition of these phases, recent computational studies of their phase stability using density functional theory (DFT) are not in agreement with reported experimental observations. In this work, we re‐examine the phase stability of the zeta and eta phases using DFT combined with enumerated searches using the known metal atom stacking sequences of these phases but with variable carbon concentration and ordering. We have found new structures for the zeta and eta phases that are now in better agreement with experimental findings. Furthermore, we report a new eta phase,‐Hf12N7, which lies on the convex hull and has a nitrogen atom ordering that is substantially different from the zeta phase. This work also demonstrates the importance of configurational entropy in dictating the finite temperature phase diagrams in this system.

     
    more » « less
  4.  
    more » « less
  5. We report the implementation of a symmetry-adapted perturbation theory algorithm based on a density functional theory [SAPT(DFT)] description of monomers. The implementation adopts a density-fitting treatment of hybrid exchange–correlation kernels to enable the description of monomers with hybrid functionals, as in the algorithm by Bukowski, Podeszwa, and Szalewicz [Chem. Phys. Lett. 414, 111 (2005)]. We have improved the algorithm by increasing numerical stability with QR factorization and optimized the computation of the exchange–correlation kernel with its 2-index density-fitted representation. The algorithm scales as O( N 5 ) formally and is usable for systems with up to ∼3000 basis functions, as demonstrated for the C 60 –buckycatcher complex with the aug-cc-pVDZ basis set. The hybrid-kernel-based SAPT(DFT) algorithm is shown to be as accurate as SAPT(DFT) implementations based on local effective exact exchange potentials obtained from the local Hartree–Fock (LHF) method while avoiding the lower-scaling [ O( N 4 )] but iterative and sometimes hard-to-converge LHF process. The hybrid-kernel algorithm outperforms Hartree–Fock-based SAPT (SAPT0) for the S66 test set, and its accuracy is comparable to the many-body perturbation theory based SAPT2+ approach, which scales as O( N 7 ), although SAPT2+ exhibits a more narrow distribution of errors. 
    more » « less