Search

Partial charges are a central concept in general chemistry and chemical biology, yet dozens of different computational definitions exist. In prior work [Cho et al., ChemPhysChem 21, 688-696 (2020)], we showed that these can be reduced to at most three 'principal components of ionicity'. The present study addressed the dependence of computed partial charges q on 1-particle basis set and (for WFT methods) n -particle correlation treatment or (for DFT methods) exchange-correlation functional, for several representative partial charge definitions such as QTAIM, Hirshfeld, Hirshfeld-I, HLY (electrostatic), NPA, and GAPT. Our findings show that semi-empirical double hybrids can closely approach the CCSD(T) 'gold standard' for this property. In fact, owing to an error compensation in MP2, CCSD partial charges are further away from CCSD(T) than is MP2. The nonlocal correlation is important, especially when there is a substantial amount of nonlocal exchange. Employing range separation proves to be "mostly" not advantageous, while global hybrids perform optimally for 20%-30% Hartree-Fock exchange across all charge types. Basis set convergence analysis shows that an augmented triple-zeta heavy-aug-cc-pV(T+d)Z basis set or a partially augmented jun-cc-pV(T+d)Z basis set is sufficient for Hirshfeld, Hirshfeld-I, HLY, and GAPT charges. In contrast, QTAIM and NPA display slower basis set convergence. It is noteworthy that for both NPA and QTAIM, HF exhibits markedly slower basis set convergence than the correlation components of MP2 and CCSD. Triples corrections in CCSD(T), denoted as CCSD(T)-CCSD, exhibit even faster basis set convergence., (© 2024 The Authors. Journal of Computational Chemistry published by Wiley Periodicals LLC.)

Hydrocarbon cages are key reference materials for the validation and parameterization of computationally cost-effective procedures such as density functional theory (DFT), semiempirical molecular orbital theory, and molecular mechanics. We obtain accurate total atomization energies (TAEs) and heats of formation (Δf H°298) for platonic and prismatic hydrocarbon cages by means of the Wn-F12 explicitly correlated thermochemical protocols. We consider the following kinetically stable (CH)n polycyclic hydrocarbon cages: (i) platonic hydrocarbons (tetrahedrane, cubane, and dodecahedrane), (ii) prismatic hydrocarbons (triprismane, cubane, and pentaprismane), and (iii) one truncated tetrahedrane (octahedrane). Our best theoretical heat of formation for cubane (144.8 kcal mol(-1)) suggests that the experimental value adopted by the NIST thermochemical database (142.7 ± 1.2 kcal mol(-1)) should be revised upwards by ∼2 kcal mol(-1). Our best heat of formation for dodecahedrane (20.2 kcal mol(-1)) suggests that the semiexperimental value (22.4 ± 1 kcal mol(-1)) should be revised downward by ∼2 kcal mol(-1). We use our benchmark Wn-F12 TAEs to evaluate the performance of a variety of computationally less demanding composite thermochemical procedures. These include the Gaussian-n (Gn) and the complete basis set (CBS) methods. The CBS-QB3 and CBS-APNO procedures show relatively poor performance with root-mean-squared deviations (RMSDs) of 4.2 and 2.5 kcal mol(-1), respectively. The best performers of the Gn procedures are G4 and G3(MP2)B3 (RMSD = 0.5 and 0.6 kcal mol(-1), respectively), while the worst performers are G3 and G4(MP2)-6X (RMSD = 2.1 and 2.9 kcal mol(-1), respectively). Isodesmic and even homodesmotic reactions involving these species are surprisingly challenging targets for DFT computations., (© 2015 Wiley Periodicals, Inc.)

Ab initio and density functional theory (DFT) calculations were performed for obtaining fundamental vibrational frequencies of methylamine, CH3NH2, and its deuterated variants CH3ND2, CD3NH2, and CD3ND2. The calculations were carried out using the CCSD(T) coupled cluster approximation with cc-pVTZ and cc-pVQZ basis sets, and by the DFT method with the semiempirical hybrid functional B97-1 with polarization consistent pc-2 and pc-3 basis sets. Reasonable performance of the DFT harmonic and ab initio harmonic calculations was found, which improved considerably upon combination of the harmonic fundamental frequencies with anharmonic corrections from the smaller, pc-2, basis. The computed anharmonic fundamental frequencies of methylamine isotopologues agree very well with the experimental values and represent a useful tool for assignment and analysis of the dominant resonances.