Benson group additivity values of phosphines and phosphine oxides: Fast and accurate computational thermochemistry of organophosphorus species

Composite quantum chemical methods W1X-1 and CBS-QB3 are used to calculate the gas phase standard enthalpy of formation, entropy, and heat capacity of 38 phosphines and phosphine oxides for which reliable experimental thermochemical information is limited or simply nonexistent. For alkyl phosphines and phosphine oxides, the W1X-1, and CBS-QB3 results are mutually consistent and in excellent agreement with available G3X values and empirical data. In the case of aryl-substituted species, different computational methods show more variation, with G3X enthalpies being furthest from experimental values. The calculated thermochemical data are subsequently used to determine Benson group additivity contributions for 24 Benson groups and group pairs involving phosphorus, thereby allowing fast and accurate estimations of thermochemical data of many organophosphorus compounds of any complexity. Such data are indispensable, for example, in chemical process design or estimating potential hazards of new chemical compounds. © 2018 Wiley Periodicals, Inc.

Thermochemistry of Guanine Tautomers Re-Examined by Means of High-Level CCSD(T) Composite Ab Initio Methods

We obtained accurate gas-phase tautomerization energies for a set of 14 guanine tautomers by means of high-level thermochemical procedures approximating the CCSD(T) energy at the complete basis set (CBS) limit. For the five low-lying tautomers, we use the computationally demanding W1-F12 composite method for obtaining the tautomerization energies. The relative W1-F12 tautomerization enthalpies at 298K are: 0.00 (1), 2.37 (2), 2.63 (3), 4.03 (3′), and 14.31 (4) kJmol−1. Thus, as many as four tautomers are found within a small energy window of less than 1.0kcalmol−1 (1kcalmol−1=4.184kJmol−1). We use these highly accurate W1-F12 tautomerization energies to evaluate the performance of a wide range of lower-level composite ab initio procedures. The Gn composite procedures (G4, G4(MP2), G4(MP2)-6X, G3, G3B3, G3(MP2), and G3(MP2)B3) predict that the enol tautomer (3) is more stable than the keto tautomer (2) by amounts ranging from 0.36 (G4) to 1.28 (G3(MP2)) kJmol−1. We also find that an approximated CCSD(T)/CBS energy calculated as HF/jul-cc-pV{D,T}Z+CCSD/jul-cc-pVTZ+(T)/jul-cc-pVDZ results in a root-mean-square deviation (RMSD) of merely 0.11kJmol−1 relative to the W1-F12 reference values. We use this approximated CCSD(T)/CBS method to obtain the tautomerization energies of 14 guanine tautomers. The relative tautomerization enthalpies at 298K are: 0.00 (1), 2.20 (2), 2.51 (3), 4.06 (3′), 14.30 (4), 25.65 (5), 43.78 (4′), 53.50 (6′), 61.58 (6), 77.37 (7), 82.52 (8′), 86.02 (9), 100.70 (10), and 121.01 (8) kJmol−1. Using these tautomerization enthalpies, we evaluate the performance of standard and composite methods for the entire set of 14 guanine tautomers. The best-performing procedures emerge as (RMSDs are given in parentheses): G4(MP2)-6X (0.51), CCSD(T)+ΔMP2/CBS (0.52), and G4(MP2) (0.64kJmol−1). The worst performers are CCSD(T)/AVDZ (1.05), CBS-QB3 (1.24), and CBS-APNO (1.38kJmol−1).

A Trip to the Density Functional Theory Zoo: Warnings and Recommendations for the User

This account is written for general users of density functional theory (DFT) methods as well as experimental researchers who are new to the field and would like to conduct such calculations. Its main emphasis lies on how to find a way through the confusing ‘zoo’ of DFT by addressing common misconceptions and highlighting those modern methods that should ideally be used in calculations of energetic properties and geometries. A particular focus is on highly popular methods and the important fact that popularity does not imply accuracy. In this context, we present a new analysis of the openly available data published in Swart and co-workers’ famous annual ‘DFT poll’ (http://www.marcelswart.eu/dft-poll/) to demonstrate the existing communication gap between the DFT user and developer communities. We show that despite considerable methodological advances in the field, the perception of some parts of the user community regarding their favourite approaches has changed little. It is hoped that this account makes a contribution towards changing this status and that users are inspired to adjust their current computational protocols to accommodate strategies that are based on proven robustness, accuracy, and efficiency rather than popularity.

Performance of DFT for C60 Isomerization Energies: A Noticeable Exception to Jacob's Ladder

The ability to accurately calculate relative energies of fullerenes is important in many areas of computational nanotechnology. Because of the large size of fullerenes, their relative energies cannot normally be calculated by means of high-level ab initio procedures, and therefore, density functional theory (DFT) represents a cost-effective alternative. In an extensive benchmark study, we calculate the electronic energies of eight C₆₀ isomers by means of the high-level G4(MP2) composite procedure. G4(MP2) isomerization energies span a wide range between 307.5 and 1074.0 kJ mol^-1. We use these benchmark data to assess the performance of DFT, double-hybrid DFT (DHDFT), and MP2-based ab initio methods. Surprisingly, functionals from the second and third rungs of Jacob's Ladder (i.e., GGA and meta-GGA functionals) significantly and systematically outperform hybrid and hybrid-meta-GGA functionals, which occupy higher rungs of Jacob's Ladder. In addition, DHDFT functionals do not offer a substantial improvement over meta-GGA functionals, with respect to isomerization energies. Overall, the best performing functionals with mean absolute deviations (MADs) below 15.0 kJ mol^-1 are (MADs given in parentheses) the GGA N12 (14.7); meta-GGAs M06-L (10.6), M11-L (10.8), MN15-L (11.9), and TPSS-D3BJ (12.8); and the DHDFT functionals B2T-PLYP (9.3), mPW2-PLYP (9.8), B2K-PLYP (12.1), and B2GP-PLYP (12.3 kJ mol^-1). In light of these results, we recommend the use of meta-GGA functionals for the calculation of fullerene isomerization energies. Finally, we show that inclusion of very small percentages of exact Hartree-Fock exchange (3-5%) slightly improves the performance of the GGA and meta-GGA functionals. However, their performance rapidly deteriorates with the inclusion of larger percentages of exact Hartree-Fock exchange.

Efficient Estimation of Formation Enthalpies for Closed-Shell Organic Compounds with Local Coupled-Cluster Methods

Efficient estimation of the enthalpies of formation for closed-shell organic compounds via atom-equivalent-type computational schemes and with the use of different local coupled-cluster with single, double, and perturbative triple excitation (CCSD(T)) approximations was investigated. Detailed analysis of established sources of uncertainty, inclusive of contributions beyond frozen-core CCSD(T) and errors due to local CCSD(T) approximations and zero-point energy anharmonicity, suggests the lower limit of about 2 kJ·mol-1 for the expanded uncertainty of the proposed estimation framework. Among the tested computational schemes, the best-performing cases demonstrate expanded uncertainty of about 2.5 kJ·mol-1, based on the analysis against 44 critically evaluated experimental values. Computational efficiency, accuracy commensurable with that of a typical experiment, and absence of the need for auxiliary reactions and additional experimental data offer unprecedented advantages for practical use, such as prompt validation of existing measurements and estimation of missing values, as well as resolution of experimental conflicts. The utility of the proposed methodology was demonstrated using a representative sample of the most recent experimental measurements.

Accurate heats of formation of polycyclic saturated hydrocarbons predicted by using the XYG3 type of doubly hybrid functionals

Polycyclic saturated hydrocarbons (PSHs) are attractive candidates as hydrocarbon propellants. To assess their potential values, one of the key factors is to determine their energy contents, such as to calculate their heats of formation (HOF). In this work, we have calculated HOFs for a set of 36 PSHs including exo-Tricyclo[5.2.1.0^(2,6) ] decane, the principal component of the high-energy density hydrocarbon fuel commonly identified as JP-10. The results from B3LYP, B3LYP-D3BJ, M06-2X, B2PLYP, B2PLYP-D3BJ, and the XYG3 type of doubly hybrid (xDH) functionals are presented. It is demonstrated here that the xDH functionals yield accurate HOFs in good agreement with those from experiments or the G4 theory. In particular, XYGJ-OS, a low scaling xDH functional, is shown to hold the promise for accurate prediction of HOFs for PSHs of larger sizes. © 2018 Wiley Periodicals, Inc.

Gas-Phase Heat of Formation Values for Buckminsterfullerene (C60), C70 Fullerene (C70), Corannulene, Coronene, Sumanene, and Other Polycyclic Aromatic Hydrocarbons Calculated Using Density Functional Theory (M06 2X) Coupled with a Versatile Inexpensive Group-Equivalent Approach

A straightforward procedure using density functional theory (M06 2X) coupled with a group-equivalent approach is described that was used to calculate gas-phase heat of formation (Δ_f H°_g,298) values for buckminsterfullerene (C₆₀), C70 fullerene (C₇₀), corannulene, coronene, and sumanene. This procedure was also used to calculate exceptionally accurate Δ_f H°_g,298 values for a variety of single-ring aromatic and 2-7 ring polycyclic aromatic hydrocarbons (PAHs) as well as a large selection of other hydrocarbons and phenols. The approach described herein is internally consistent, and results for C₆₀, C₇₀, corannulene, coronene, and sumanene are in very close agreement with results reported by others who used higher-level computational theory. Statistical analysis of a test set containing benzene and 18 two to seven ring PAHs demonstrated that by using this approach a mean absolute deviation (MAD) and a root-mean-square deviation (RMSD) of 0.8 and 1.3 kJ/mol, respectively, were achieved for reference/experimental Δ_f H°_g,298 values versus calculated/predicted Δ_f H°_g,298 values. For statistical analysis of a larger test set containing 235 aromatic and aliphatic hydrocarbons and phenols, a MAD and a RMSD of 1.2 and 1.9 kJ/mol, respectively, were achieved for reference/experimental Δ_f H°_g,298 values versus calculated/predicted Δ_f H°_g,298 values.

Physicochemical Prediction of Metabolite Fragmentation in Tandem Mass Spectrometry

Current bottleneck of comprehensive non-target metabolite identification is insufficient spectral library. Many research groups have tried to build a theoretical product ion spectral library independent of measurement condition or settings, but mechanisms of metabolite fragmentation are not fully clarified. To achieve the mechanistic prediction of metabolite fragmentation which covers a wide range of metabolites, we will discuss utilization of physicochemical calculation. We introduce bonding patterns, which include two bound atoms and chemical groups adjacent to the bond. Cleavage of each bonding pattern is simulated and its activation energy is precisely calculated with quantum chemistry and assigned on metabolites. By tracing low-energy bond cleavages, fragmentation of a dipeptide molecule is successfully predicted. Prediction on another metabolite requires some additional features to fully reproduce its experimentally observed product ions. Physicochemical calculation shows its promising ability to predict fragmentation pathways only from metabolite structures, while required improvements suggested by comparison between our prediction and standard spectra stored in database are also discussed. Moreover, to construct a prediction strategy which withstands the vast metabolite space, we have to build a comprehensive list of bonding patterns and their activation energy. As theoretically possible bonding patterns are huge in number, proper simplification of the patterns must be implemented. We will discuss how to achieve it in addition to the prediction results.

Can Popular DFT Approximations and Truncated Coupled Cluster Theory Describe the Potential Energy Surface of the Beryllium Dimer?

The potential energy surface (PES) of the ground state of the beryllium dimer poses a significant challenge for high-level ab initio electronic structure methods. Here, we present a systematic study of basis set effects over the entire PES of Be2 calculated at the full configuration interaction (FCI) level. The reference PES is calculated at the valence FCI/cc-pV{5,6}Z level of theory. We find that the FCI/cc-pV{T,Q}Z basis set extrapolation reproduces the shape of the FCI/cc-pV{5,6}Z PES as well as the binding energy and vibrational transition frequencies to within ~10 cm−1. We also use the FCI/cc-pV{5,6}Z PES to evaluate the performance of truncated coupled cluster methods (CCSD, CCSD(T), CCSDT, and CCSDT(Q)) and contemporary density functional theory methods (DFT) methods for the entire PES of Be2. Of the truncated coupled cluster methods, CCSDT(Q)/cc-pV{5,6}Z provides a good representation of the FCI/cc-pV{5,6}Z PES. The GGA functionals, as well as the HGGA and HMGGA functionals with low percentages of exact exchange tend to severely overbind the Be2 dimer, whereas BH&HLYP and M06-HF tend to underbind it. Range-separated DFT functionals tend to underbind the dimer. Double-hybrid DFT functionals show surprisingly good performance, with DSD-PBEP86 being the best performer. Møller–Plesset perturbation theory converges smoothly up to fourth order; however, fifth-order corrections have practically no effect on the PES.

The S66 Non-Covalent Interactions Benchmark Reconsidered Using Explicitly Correlated Methods Near the Basis Set Limit

The S66 benchmark for non-covalent interactions has been re-evaluated using explicitly correlated methods with basis sets near the one-particle basis set limit. It is found that post-MP2 ‘high-level corrections’ are treated adequately well using a combination of CCSD(F12*) with (aug-)cc-pVTZ-F12 basis sets on the one hand, and (T) extrapolated from conventional CCSD(T)/heavy-aug-cc-pV{D,T}Z on the other hand. Implications for earlier benchmarks on the larger S66×8 problem set in particular, and for accurate calculations on non-covalent interactions in general, are discussed. At a slight cost in accuracy, (T) can be considerably accelerated by using sano-V{D,T}Z+ basis sets, whereas half-counterpoise CCSD(F12*)(T)/cc-pVDZ-F12 offers the best compromise between accuracy and computational cost.

Revisiting the thermochemistry of chlorine fluorides

In this work, accurate calculations of standard enthalpies of formation of chlorine fluorides (ClF_n, n = 1-7; Cl₂ F and Cl₃ F₂ ) were performed through the isodesmic reactions scheme. It is argued that, for many chlorine fluorides, the gold standard method of quantum chemistry (CCSD(T)) is not capable to predict enthalpy values nearing chemical accuracy if atomization scheme is used. This is underpinned by a thorough analysis of total atomization energy results and the inspection of multireference features of these compounds. Other thermodynamic quantities were also calculated at different temperatures. To complement the energetic description, elimination curves were studied through density functional theory as a computationally affordable alternative to highly correlated wave function-based methods. © 2017 Wiley Periodicals, Inc.

How to computationally calculate thermochemical properties objectively, accurately, and as economically as possible

We have developed the WnX series of quantum chemistry composite protocols for the computation of highly-accurate thermochemical quantities with advanced efficiency and applicability. The W1X-type methods have a general accuracy of ~3–4 kJ mol−1 and they can currently be applied to systems with ~20–30 atoms. Higher-level methods include W2X, W3X and W3X-L, with the most accurate of these being W3X-L. It can be applied to molecules with ~10–20 atoms and is generally accurate to ~1.5 kJ mol−1. The WnX procedures have opened up new possibilities for computational chemists in pursue of accurate thermochemical values in a highly-productive manner.

Unification of the W1X and G4(MP2)-6X Composite Protocols

We have devised the composite procedures WG and WGh to unify the W1X and the (computationally more economical) G4(MP2)-6X protocols. The WG procedure employs a combination of MP2, MP2-F12, CCSD-F12b, and CCSD(T) to approximate the all-electron scalar-relativistic CCSD(T)/CBS energy. In addition, it incorporates features such as the scaling of the energy components and an empirical "higher-level-correction" term. The WGh protocol represents a somewhat more economical variant of WG with partial removal of diffuse functions. Our benchmark shows that, in general, both WG and WGh have similar performance to that for W1X-2, with WGh (predictably) performing somewhat less well for electron affinities. In terms of computational efficiency, WG is approximately an order of magnitude less costly than W1X-2, while WGh gives not only a further slight savings in computer time but also a notably reduced disk requirement.

Efficient DLPNO–CCSD(T)-Based Estimation of Formation Enthalpies for C-, H-, O-, and N-Containing Closed-Shell Compounds Validated Against Critically Evaluated Experimental Data

An accurate and cost-efficient methodology for the estimation of the enthalpies of formation for closed-shell compounds composed of C, H, O, and N atoms is presented and validated against critically evaluated experimental data. The computational efficiency is achieved through the use of the resolution-of-identity (RI) and domain-based local pair-natural orbital coupled cluster (DLPNO-CCSD(T)) approximations, which results in a drastic reduction in both the computational cost and the number of necessary steps for a composite quantum chemical method. The expanded uncertainty for the proposed methodology evaluated using a data set of 45 thoroughly vetted experimental values for molecules containing up to 12 heavy atoms is about 3 kJ·mol-1, competitive with those of typical calorimetric measurements. For the compounds within the stated scope, the methodology is shown to be superior to a representative, more general, and widely used composite quantum chemical method, G4.

Impact of Hydrogen Bonding on the Susceptibility of Peptides to Oxidation

The tendency of peptides to be oxidized is intimately connected with their function and even their ability to exist in an oxidative environment. Here we report high-level theoretical studies that show that hydrogen bonding can alter the susceptibility of peptides to oxidation, with complexation to a hydrogen-bond acceptor facilitating oxidation, and vice versa, impacting the feasibility of a diverse range of biological processes. It can even provide an energetically viable mechanistic alternative to direct hydrogen-atom abstraction. We find that hydrogen bonding to representative reactive groups leads to a broad (≈400 kJ mol^-1 ) spectrum of ionization energies in the case of model amide, thiol and phenol systems. While some of the oxidative processes at the extreme ends of the spectrum are energetically prohibitive, subtle environmental and solvent effects could potentially mitigate the situation, leading to a balance between hydrogen bonding and oxidative susceptibility.

Thermochemistry of icosahedral closo-dicarboranes: a composite ab initio quantum-chemical perspective

We obtained accurate thermochemical properties for the ortho-, meta-, and para-dicarborane isomers (C2B10H12) by means of explicitly correlated high-level thermochemical procedures. The thermochemical properties include heats of formation, isomerization energies, C–H and B–H bond dissociation energies (BDEs), and ionization potentials. Of these only the ionization potentials are known experimentally. Our best theoretical ionization potentials, obtained by means of the ab initio W1–F12 thermochemical protocol, was 241.50 kcal mol–1 (para-dicarborane), 238.45 kcal mol–1 (meta-dicarborane), and 236.54 kcal mol–1 (ortho-dicarborane). These values agree with the experimental values adopted by the National Institute of Standards and Technology (NIST) thermochemical tables to within overlapping uncertainties. However, they suggest that the experimental values may represent significant underestimations. For all isomers, the C–H BDEs are systematically higher than the B–H BDEs because of the relative stability of the boron-centred radicals. The C–H BDEs for the three isomers cluster within a narrow energetic interval, namely between 110.8 kcal mol–1 (para-dicarborane) and 111.7 kcal mol–1 (meta-dicarborane). The B–H BDEs cluster within a larger interval ranging between 105.8 and 108.1 kcal mol–1 (both obtained for ortho-dicarborane). We used our benchmark W1–F12 data to assess the performance of a number of lower cost composite ab initio methods. We found that the Gaussian-3 procedures (G3(MP2)B3 and G3B3) result in excellent performance with overall root-mean-square deviations (RMSDs) of 0.3–0.4 kcal mol–1 for the isomerization, ionization, and bond dissociation energies. However, the Gaussian-4 procedures (G4, G4(MP2), and G4(MP2)-6X) showed relatively poor performance with overall RMSDs of 1.3–3.7 kcal mol–1.

[Al2O4](-), a Benchmark Gas-Phase Class II Mixed-Valence Radical Anion for the Evaluation of Quantum-Chemical Methods

The radical anion [Al2O4](-) has been identified as a rare example of a small gas-phase mixed-valence system with partially localized, weakly coupled class II character in the Robin/Day classification. It exhibits a low-lying C2v minimum with one terminal oxyl radical ligand and a high-lying D2h minimum at about 70 kJ/mol relative energy with predominantly bridge-localized-hole character. Two identical C2v minima and the D2h minimum are connected by two C2v-symmetrical transition states, which are only ca. 6-10 kJ/mol above the D2h local minimum. The small size of the system and the absence of environmental effects has for the first time enabled the computation of accurate ab initio benchmark energies, at the CCSDT(Q)/CBS level using W3-F12 theory, for a class-II mixed-valence system. These energies have been used to evaluate wave function-based methods [CCSD(T), CCSD, SCS-MP2, MP2, UHF] and density functionals ranging from semilocal (e.g., BLYP, PBE, M06L, M11L, N12) via global hybrids (B3LYP, PBE0, BLYP35, BMK, M06, M062X, M06HF, PW6B95) and range-separated hybrids (CAM-B3LYP, ωB97, ωB97X-D, LC-BLYP, LC-ωPBE, M11, N12SX), the B2PLYP double hybrid, and some local hybrid functionals. Global hybrids with about 35-43% exact-exchange (EXX) admixture (e.g., BLYP35, BMK), several range hybrids (CAM-B3LYP, ωB97X-D, ω-B97), and a local hybrid provide good to excellent agreement with benchmark energetics. In contrast, too low EXX admixture leads to an incorrect delocalized class III picture, while too large EXX overlocalizes and gives too large energy differences. These results provide support for previous method choices for mixed-valence systems in solution and for the treatment of oxyl defect sites in alumosilicates and SiO2. Vibrational gas-phase spectra at various computational levels have been compared directly to experiment and to CCSD(T)/aug-cc-pV(T+d)Z data.