skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Correlation consistent basis sets designed for density functional theory: Third-row atoms (Ga–Br)
The correlation consistent basis sets (cc-pVnZ with n = D, T, Q, 5) for the Ga–Br elements have been redesigned, tuning the sets for use for density functional approximations. Steps to redesign these basis sets for an improved correlation energy recovery and efficiency include truncation of higher angular momentum functions, recontraction of basis set coefficients, and reoptimization of basis set exponents. These redesigned basis sets are compared with conventional cc-pVnZ basis sets and other basis sets, which are, in principle, designed to achieve systematic improvement with respect to increasing basis set size. The convergence of atomic energies, bond lengths, bond dissociation energies, and enthalpies of formation to the Kohn–Sham limit is improved relative to other basis sets where convergence to the Kohn–Sham limit is typically not observed.  more » « less
Award ID(s):
2154526
PAR ID:
10554468
Author(s) / Creator(s):
;
Publisher / Repository:
American Physical Society
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
160
Issue:
8
ISSN:
0021-9606
Page Range / eLocation ID:
084105
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Conditional-probability density functional theory (CP-DFT) is a formally exact method for finding correlation energies from Kohn-Sham DFT without evaluating an explicit energy functional. We present details on how to generate accurate exchange-correlation energies for the ground-state uniform gas. We also use the exchange hole in a CP antiparallel spin calculation to extract the high-density limit. We give a highly accurate analytic solution to the Thomas-Fermi model for this problem, showing its performance relative to Kohn-Sham and it may be useful at high temperatures. We explore several approximations to the CP potential. Results are compared to accurate parametrizations for both exchange-correlation energies and holes. 
    more » « less
  2. A new reference state for density functional theory (DFT), termed the independent atom ansatz, is introduced in this work. This ansatz allows for the formally exact representation of electron density in terms of atom-localized orbitals. Self-consistent equations for such states are derived in general and asymptotic forms. The resultant total energy functional is found to closely resemble tight-binding theory. The independent atom ansatz facilitates partial cancellation of inter-atomic electron–electron and nucleus–electron interactions, which allows for the derivation of analytical tight-binding Hamiltonian matrix elements in a weak interaction limit. The formalism provides energy decomposition and charge analyses at no additional cost and links tight-binding, localized orbital, and electronegativity concepts. Numerical accuracy of the total energy functional has been previously reported for hydrogenic systems [Mironenko, J. Phys. Chem. A 127, 7836 (2023)] and is demonstrated here for He2, Li2, Be2, B2, N2, O2, F2, and Ne2. The method accurately reproduces the shapes of potential energy curves, capturing large-basis CCSD(T)-level bond lengths and bond dissociation energies for N2, O2, and F2 using only a minimal basis set. It outperforms both CCSD(T) and some mainstream approximate restricted Kohn–Sham DFT functionals in describing bond dissociation behavior away from equilibrium geometries. 
    more » « less
  3. Abstract The global minima of urea and thiourea were characterized along with other low‐lying stationary points. Each structure was optimized with the CCSD(T) method and triple‐ζcorrelation consistent basis sets followed by harmonic vibrational frequency computations. Relative energies evaluated near the complete basis set limit with both canonical and explicitly correlated CCSD(T) techniques reveal several subtle but important details about both systems. These computations resolve a discrepancy by demonstrating that the electronic energy of the C2vsecond‐order saddle point of urea lies at least 1.5 kcal mol−1above the C2global minimum regardless of whether the structures were optimized with MP2, CCSD, or CCSD(T). Additionally, urea effectively has one minimum instead of two because the electronic barrier for inversion at one amino group in the Cslocal minimum vanishes at the CCSD(T) CBS limit. Characterization of both systems with the same ab initio methods and large basis sets conclusively establishes that the electronic barriers to inversion at one or both NH2groups in thiourea are appreciably smaller than in urea. CCSDT(Q)/cc‐pVTZ computations show higher‐order electron correlation effects have little impact on the relative energies and are consistently offset by core correlation effects of opposite sign and comparable magnitude. 
    more » « less
  4. Hydrazoic acid (HN3) is used as a case study for investigating the accuracy and precision by which a molecular structure—specifically, a semi-experimental equilibrium structure (reSE)—may be determined using current state-of-the-art methodology. The influence of the theoretical corrections for effects of vibration–rotation coupling and electron-mass distribution that are employed in the analysis is explored in detail. The small size of HN3 allowed us to deploy considerable computational resources to probe the basis-set dependence of these corrections using a series of coupled-cluster single, double, perturbative triple [CCSD(T)] calculations with cc-pCVXZ (X = D, T, Q, 5) basis sets. We extrapolated the resulting corrections to the complete basis set (CBS) limit to obtain CCSD(T)/CBS corrections, which were used in a subsequent reSE structure determination. The reSE parameters obtained using the CCSD(T)/cc-pCV5Z corrections are nearly identical to those obtained using the CCSD(T)/CBS corrections, with uncertainties in the bond distances and angles of less than 0.0006 Å and 0.08°, respectively. The previously obtained reSE structure using CCSD(T)/ANO2 agrees with that using CCSD(T)/cc-pCV5Z to within 0.000 08 Å and 0.016° for bond distances and angles, respectively, and with only 25% larger uncertainties, validating the idea that reSE structure determinations can be carried out with significantly smaller basis sets than those needed for similarly accurate, strictly ab initio determinations. Although the purely computational re structural parameters [CCSD(T)/cc-pCV6Z] fall outside of the statistical uncertainties (2σ) of the corresponding reSE structural parameters, the discrepancy is rectified by applying corrections to address the theoretical limitations of the CCSD(T)/cc-pCV6Z geometry with respect to basis set, electron correlation, relativity, and the Born–Oppenheimer approximation, thereby supporting the contention that the semi-experimental approach is both an accurate and vastly more efficient method for structure determinations than is brute-force computation. 
    more » « less
  5. Abstract Local correlation methods rely on the assumption that electron correlation is nearsighted. In this work, we develop a method to alleviate this assumption. This new method is demonstrated by calculating the random phase approximation (RPA) correlation energies in several one‐dimensional model systems. In this new method, the first step is to approximately decompose the RPA correlation energy to the nearsighted and farsighted components based on the wavelength decomposition of electron correlation developed by Langreth and Perdew. The short‐wavelength (SW) component of the RPA correlation energy is then considered to be nearsighted, and the long‐wavelength (LW) component of the RPA correlation energy is considered to be farsighted. The SW RPA correlation energy is calculated using a recently developed local correlation method: the embedded cluster density approximation (ECDA). The LW RPA correlation energy is calculated globally based on the system's Kohn‐Sham orbitals. This new method is termedλ‐ECDA, whereλindicates the wavelength decomposition. The performance ofλ‐ECDA is examined on a one‐dimensional model system: aH24chain, in which the RPA correlation energy is highly nonlocal. In this model system, a softened Coulomb interaction is used to describe the electron‐electron and electron‐ion interactions, and slightly stronger nuclear charges (1.2e) are assigned to the pseudo‐H atoms. Bond stretching energies, RPA correlation potentials, and Kohn‐Sham eigenvalues predicted byλ‐ECDA are in good agreement with the benchmarks when the clusters are made reasonably large. We find that the LW RPA correlation energy is critical for obtaining accurate prediction of the RPA correlation potential, even though the LW RPA correlation energy contributes to only a few percent of the total RPA correlation energy. 
    more » « less