skip to main content


This content will become publicly available on December 17, 2024

Title: High-level ab initio composites: thermochemical bond dissociation energies for vanadium species
Bond dissociation energies (BDE) are key descriptors for molecules and are among the most sought-after properties in chemistry. Despite their importance, the accurate prediction of BDE’s for transition metal species can be particularly daunting for both experiment and computation. Experimental data has been limited and, when available, often has large error bars, making the critical evaluation and identification of suitable computational methods difficult. However, recent advancements in the experimental determination of BDE’s with techniques such as Velocity Map Imaging and 2 Photon Ionization now provide useful gauges for computational strategies and new methodologies, providing energies with unprecedented accuracies. The vanadium diatomics (VX, X=B, C, N, O, F, Al, Si, P, S, Cl) have been challenging for computational chemistry methods, and, thus, a new experimental gauge enables methods to be reevaluated and developed for these species. Herein, the super-correlation consistent Composite (super-ccCA or s-ccCA), a new thermochemical scheme centered around CCSD(T)/complete basis set (CBS) limit computations with additional contributions that account for scalar-relativistic effects, and coupled cluster contributions beyond CCSD(T) up to quintuple excitations has been considered. The agreement between determinations made by the s-ccCA scheme and by recent experiment is excellent, demonstrating the utility of the new approach in addressing challenging metal systems, even those of multireference nature. In light of recent experimental BDE’s, the longstanding correlation consistent composite approach (ccCA) is also evaluated for the VX species and find that the mean absolute deviation (MAD) is greatly reduced compared to previously used experimental values.  more » « less
Award ID(s):
1900086
NSF-PAR ID:
10486899
Author(s) / Creator(s):
; ;
Publisher / Repository:
Taylor & Francis
Date Published:
Journal Name:
Molecular Physics
Volume:
121
Issue:
24
ISSN:
0026-8976
Subject(s) / Keyword(s):
Electronic structure transition metals thermochemistry ab initio composite s-ccCA
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We have developed a new composite model chemistry method called WMS (Wuhan–Minnesota scaling method) with three characteristics: (1) a composite scheme to approximate the complete configuration interaction valence energy with the affordability condition of requiring no calculation more expensive than CCSD(T)/jul-cc-pV(T+d)Z, (2) low-cost methods for the inner-shell correlation contribution and scalar relativistic correction, and (3) accuracy comparable to methods with post-CCSD(T) components. The new method is shown to be accurate for the W4-17 database of 200 atomization energies with an average mean unsigned error (averaged with equal weight over strongly correlated and weakly correlated subsets of the data) of 0.45 kcal mol −1 , and the performance/cost ratio of these results compares very favorably to previously available methods. We also assess the WMS method against the DBH24-W4 database of diverse barrier heights and the energetics of the reactions of three strongly correlated Criegee intermediates with water. These results demonstrate that higher-order correlation contributions necessary to obtain high accuracy for molecular thermochemistry may be successfully extrapolated from the lower-order components of CCSD(T) calculations, and chemical accuracy can now be obtained for larger and more complex molecules and reactions. 
    more » « less
  2. Abstract

    Thes‐homodesmotic method for computing conventional strain energies (CSE) has been extended for the first time to bicyclic systems and to individual rings within these systems. Unique isodesmic, homodesmotic, and hyperhomodesmotic reactions originate from thes‐homodesmotic method. These are used to investigate 12 bicyclic systems comprising cyclopropane and cyclobutane and how the CSE of each system compares to the sum of the individual rings within each. Equilibrium geometries, harmonic vibrational frequencies, and the corresponding electronic energies and zero point vibrational energy corrections are computed for all relevant molecules using second‐order perturbation theory and density functional theory (B3LYP) with the correlation consistent basis sets cc‐pVDZ and cc‐pVTZ. Single‐point CCSD(T) energies are computed at the MP2/cc‐pVTZ optimized geometries to ascertain the importance of higher order correlation effects. Results indicate that CSEs are additive when the two rings are separated by one or two bonds and somewhat additive in other cases.

     
    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. It is demonstrated that a double hybrid density functional approximation, ωB88PTPSS, that incorporates equipartition of density functional theory and the non-local correlation, however with a meta-generalized gradient approximation correlation functional, as well as with the range-separated exchange of ωB2PLYP, provides accurate excitation energies for conventional systems, as well as correct prescription of negative singlet–triplet gaps for non-conventional systems with inverted gaps, without any necessity for parametric scaling of the same-spin and opposite-spin non-local correlation energies. Examined over “safe” excitations of the QUESTDB set, ωB88PTPSS performs quite well for open-shell systems, correctly and fairly accurately [relative to equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) reference] predicts negative gaps for 50 systems with inverted singlet–triplet gaps, and is one of the leading performers for intramolecular charge-transfer excitations and achieves near-second-order approximate coupled cluster (CC2) and second-order algebraic diagrammatic construction quality for the Q1 and Q2 subsets. Subsequently, we tested ωB88PTPSS on two sets of real-life examples from recent computational chemistry literature–the low energy bands of chlorophyll a (Chl a) and a set of thermally activated delayed fluorescence (TADF) systems. For Chl a, ωB88PTPSS qualitatively and quantitatively achieves DLPNO-STEOM-CCSD-level performance and provides excellent agreement with experiment. For TADF systems, ωB88PTPSS agrees quite well with spin-component-scaled CC2 (SCS-CC2) excitation energies, as well as experimental values, for the gaps between the S1 and T1 excited states. 
    more » « less
  5. To study the contribution of three-body dispersion to crystal lattice energies, we compute the three-body contributions to the lattice energies for crystalline benzene, carbon dioxide, and triazine using various computational methods. We show that these contributions converge quickly as the intermolecular distances between the monomers grow. In particular, the smallest value among the three pairwise intermonomer closest-contact distances, Rmin, shows a strong correlation with the three-body contribution to the lattice energy, and, here, the largest of the closest-contact distances, Rmax, serves as a cutoff criterion to limit the number of trimers to be considered. We considered all trimers up to Rmax=15Å. The trimers with Rmin<4Å contribute 90.4%, 90.6%, and 93.9% of the total three-body contributions for crystalline benzene, carbon dioxide, and triazine, respectively, for the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] method. For trimers with Rmin>4Å, the second-order Møller–Plesset perturbation theory (MP2) supplemented with the Axilrod–Teller–Muto (ATM) three-body dispersion correction reproduces the CCSD(T) values for the cumulative three-body contributions with errors of less than 0.1 kJ mol−1. Moreover, three-body contributions are converged within 0.15 kJ mol−1 by Rmax=10Å. From these results, it appears that in molecular crystals where dispersion dominates the three-body contribution to the lattice energy, the trimers with Rmin>4Å can be computed with the MP2+ATM method to reduce the computational cost, and those with Rmax>10Å appear to be basically negligible. 
    more » « less