High-accuracy Extrapolated Ab Initio Thermochemistry Research Articles

W4Λ: Leveraging Λ Coupled-Cluster for Accurate Computational Thermochemistry Approaches.

High-accuracy composite wave function methods like Weizmann-4 (W4) theory, high-accuracy extrapolated ab initio thermochemistry (HEAT), and the Feller-Peterson-Dixon (FPD) approach enable sub-kJ/mol accuracy in gas-phase thermochemical properties. Their biggest computational bottleneck is the evaluation of the valence post-CCSD(T) correction term. We demonstrate here, for the W4-17 thermochemistry benchmark and subsets thereof, that the Λ coupled-cluster expansion converges more rapidly and smoothly than the regular coupled-cluster series. By means of CCSDT(Q)Λ and CCSDTQ(5)Λ, we can considerably (up to an order of magnitude) accelerate W4- and W4.3-type calculations without loss in accuracy, leading to the W4Λ and W4.3Λ computational thermochemistry protocols.

Unraveling the electronic structure of LuH, LuN, and LuNH: building blocks of new materials.

Advances in superconductor technology have been pursued for decades, moving towards room temperature models, such as a postulated nitrogen-doped lutetium hydride network. While experimental observations have been contradictory, insight into the building blocks of potential new superconductor materials can be gained theoretically, unravelling the fascinating electronic structure of these compounds at a molecular level. Here, the fundamental building blocks of lutetium materials (LuH, LuN, and LuNH) have been examined. The structures, spectroscopic constants for the ground and excited states, and the potential energy curves have been obtained for these species using complete active self-consistent field (CASSCF) and multireference configuration interaction with Davidson's correction (MRCI+Q) methods. For LuNH, the energetic properties of its isomers are determined. The bond dissociation energies of the three building blocks are calculated with the state-of-the-art f-block ab initio correlation consistent composite approach (f-ccCA) and the high accuracy extrapolated ab initio thermochemistry (HEAT) scheme. As well, an analysis of different formation pathways of LuNH has been provided.

Sub 20 cm-1 computational prediction of the CH bond energy - a case of systematic error in computational thermochemistry.

The bond dissociation energy of methylidyne, D0(CH), is studied using an improved version of the High-Accuracy Extrapolated ab initio Thermochemistry (HEAT) approach as well as the Feller-Peterson-Dixon (FPD) model chemistry. These calculations, which include basis sets up to nonuple (aug-cc-pCV9Z) quality, are expected to be capable of providing results substantially more accurate than the ca. 1 kJ mol-1 level that is characteristic of standard high-accuracy protocols for computational thermochemistry. The calculated 0 K CH bond energy (27 954 ± 15 cm-1 for HEAT and 27 956 ± 15 cm-1 for FPD), along with equivalent treatments of the CH ionization energy and the CH+ dissociation energy (85 829 ± 15 cm-1 and 32 946 ± 15 cm-1, respectively), were compared to the existing benchmarks from Active Thermochemical Tables (ATcT), uncovering an unexpected difference for D0(CH). This has prompted a detailed reexamination of the provenance of the corresponding ATcT benchmark, allowing the discovery and subsequent correction of a systematic error present in several published high-level calculations, ultimately yielding an amended ATcT benchmark for D0(CH). Finally, the current theoretical results were added to the ATcT Thermochemical Network, producing refined ATcT estimates of 27 957.3 ± 6.0 cm-1 for D0(CH), 32 946.7 ± 0.6 cm-1 for D0(CH+), and 85 831.0 ± 6.0 cm-1 for IE(CH).

Elaborated thermochemical treatment of HF, CO, N2, and H2O: Insight into HEAT and its extensions.

Empirical, highly accurate non-relativistic electronic total atomization energies (eTAEs) are established by combining experimental or computationally converged treatments of the nuclear motion and relativistic contributions with the total atomization energies of HF, CO, N2, and H2O obtained from the Active Thermochemical Tables. These eTAEs, which have estimated (2σ) uncertainties of less than 10 cm-1 (0.12 kJ mol-1), form the basis for an analysis of high-level ab initio quantum chemical calculations that aim at reproducing these eTAEs for the title molecules. The results are then employed to analyze the performance of the high-accuracy extrapolated ab initio thermochemistry, or High-Accuracy Extrapolated Ab Initio Thermochemistry (HEAT), family of theoretical methods. The method known as HEAT-345(Q), in particular, is found to benefit from fortuitous error cancellation between its treatment of the zero-pointenergy, extrapolation errors in the Hartree-Fock and coupled cluster contributions, neglect of post-(T) core-correlation, and the basis-set error involved in higher-level correlation corrections. In addition to shedding light on a longstanding curiosity of the HEAT protocol-where the cheapest HEAT-345(Q) performs comparably to the theoretically more complete HEAT-456QP procedure-this study lays the foundation for extended HEAT variants that offer substantial improvements in accuracy relative to the established approaches.

Semiclassical transition state theory/master equation kinetics of HO + CO: Performance evaluation

AbstractPreviously, master equation (ME) simulations using semiclassical transition state theory (SCTST) and high‐accuracy extrapolated ab initio thermochemistry (HEAT) predicted rate constants in excellent agreement with published experimental data over a wide range of pressure and temperatures ≳250 K, but the agreement was not as good at lower temperatures. Possible reasons for this reduced performance are investigated by (a) critically evaluating the published experimental data and by investigating; (b) three distinct ME treatments of angular momentum, including one that is exact at the zero‐ and infinite‐pressure limits; (c) a hindered‐rotor model for HOCO that implicitly includes the cis‐ and trans‐conformers; (d) possible empirical adjustments of the thermochemistry; (e) possible empirical adjustments to an imaginary frequency controlling tunneling; (f) including or neglecting the prereaction complex PRC1; and (g) its possible bimolecular reactions. Improvements include better approximations to factors in SCTST and using the Hill and van Vleck treatment of angular momentum coupling. Evaluation of literature data does not reveal any specific shortcomings, but the stated uncertainties may be underestimated. All ME treatments give excellent fits to experimental data at T ≥ 250 K, but the discrepancy at T < 250 K persists. Note that each ME model requires individual empirical energy transfer parameters. Thermochemical adjustments were unable to match the experimental H/D kinetic isotope effects. Adjusting an imaginary frequency can achieve good fits, but the adjustments are unacceptably large. Whether PRC1 and its possible bimolecular reactions are included had little effect. We conclude that none of the adjustments is an improvement over the unadjusted theory. Note that only one set of experimental data exists in the regime of the discrepancy with theory, and data for DO + CO are scanty.

High-accuracy extrapolated ab initio thermochemistry. IV. A modified recipe for computational efficiency.

A number of economical modifications to the high-accuracy extrapolated ab initio thermochemistry (HEAT) model chemistry are evaluated. The two resulting schemes, designated as mHEAT and mHEAT+, are designed for efficient and pragmatic evaluation of molecular energies in systems somewhat larger than can be practically studied by the unapproximated HEAT scheme. It is found that mHEAT+ produces heats of formation with nearly subchemical (±1 kJ/mol) accuracy at a substantially reduced cost relative to the full scheme. Total atomization energies calculated using the new thermochemical recipes are compared to the results of the HEAT-345(Q) model chemistry, and enthalpies of formation for the three protocols are also compared to Active Thermochemical Tables. Finally, a small selection of transition states is studied using mHEAT and mHEAT+, which illuminates some interesting features of reaction barriers and serves as an initial benchmark of the performance of these model chemistries for chemical kinetics applications.

A master equation simulation for the •OH + CH3OH reaction.

A combined (fixed-J) two-dimensional master-equation/semi-classical transition state theory/variational Rice-Ramsperger-Kassel-Marcus approach has been used to compute reaction rate coefficients of •OH with CH3OH over a wide range of temperatures (10-2500 K) and pressures (10-1-104 Torr) based on a potential energy surface that has been constructed using a modification of the high accuracy extrapolated ab initio thermochemistry (HEAT) protocol. The calculated results show that the title reaction is nearly pressure-independent when T > 250 K but depends strongly on pressure at lower temperatures. In addition, the preferred mechanism and rate constants are found to be very sensitive to temperature. The reaction pathway CH3OH + •OH → CH3O• + H2O proceeds exclusively through tunneling at exceedingly low temperatures (T ≤ 50 K), typical of those established in interstellar environments. In this regime, the rate constant is found to increase with decreasing temperature, which agrees with low-temperature experimental results. The thermodynamically favored reaction pathway CH3OH + •OH → •CH2OH + H2O becomes dominant at higher temperatures (T ≥ 200 K), such as those found in Earth's atmosphere as well as combustion environments. By adjusting the ab initio barrier heights slightly, experimental rate constants from 200 to 1250 K can be satisfactorily reproduced.

Ab initio thermal rate coefficients for H + NH3 ⇌ H2 + NH2

AbstractThe reversible reaction NH3 + H ⇌ H2 + NH2, which plays an important role in NH3 fuel combustion, is studied with a theoretical approach that combines the high‐accuracy extrapolated ab initio thermochemistry (HEAT) protocol with semiclassical transition state theory (SCTST). The calculated forward reaction is endothermic by 11.8 ± 1 kJ/mol, in nearly perfect agreement with the active thermochemical tables (ATcT) value of 11.5 ± 0.2 kJ/mol. Using this improved thermochemistry yields better rate constants, especially at low temperatures. Experimental rate constants available from 400 to 2000 K for the forward and reverse reaction pathways can be reproduced (within 20%) by the calculations from first principles.

Bethe-Salpeter correlation energies of atoms and molecules.

A variety of approaches are presented for the computation of atomic and molecular correlation energies based on the Bethe-Salpeter equation in the framework of the adiabatic-connection fluctuation-dissipation theorem. The performance of the approaches is assessed by computing the total energies of the atoms H-Ne and the atomization energies of the high-accuracy extrapolated ab initio thermochemistry set of small molecules as well as by determining the bond lengths and harmonic vibrational frequencies of the metal monoxides MO with M=Ca-Zn.

High-accuracy extrapolated ab initio thermochemistry of the NCN radical

Enthalpies of formation of the cyanonitrene radical, NCN, were calculated according to the HEAT (High-accuracy Extrapolated Ab initio Thermochemistry) protocol. The obtained, recommended values are 451.7 ± 1.0 kJ mol−1 and 452.2 ± 1.0 kJ mol−1 at 0 K and 298.15 K, respectively (gas phase, 1 bar). The calculated values are in agreement with the most recent Active Thermochemical Tables analysis and some of the previous experimental and computational results. The uncertainties associated with the HEAT values are, however, considerably smaller.

Accurate Determination of Barrier Height and Kinetics for the F + H2O → HF + OH Reaction

The reaction energy and barrier height of the title reaction are investigated using two high-level ab initio protocols, namely Focal Point Analysis (FPA) and modified High Accuracy Extrapolated Ab Initio Thermochemistry (HEAT) methods. It is concluded from these calculations that despite some multireference character, dynamic electron correlation plays a dominant role near the reaction barrier. Thus, the coupled-cluster method with higher excitations than singles and doubles gives a better description than the multireference configuration interaction method for the barrier height. The FPA and HEAT classical barrier heights, including the spin-orbit and other corrections, are 1.919 and 2.007 kcal/mol, respectively. The rate constants and H/D kinetic isotope effect for the title reaction are determined by semiclassical transition-state theory based on the anharmonic potential energy surface near the saddle point, and the agreement with experiment is excellent. The rate constants are also computed using a quasi-classical trajectory method on a global potential energy surface scaled to the FPA barrier height and a similar level of agreement with experimental data is obtained.

High-Accuracy Extrapolated Ab Initio Thermochemistry of the Vinyl, Allyl, and Vinoxy Radicals

Enthalpies of formation at both 0 and 298 K were calculated according to the HEAT (High-accuracy Extrapolated Ab initio Thermochemistry) protocol for the title molecules, all of which play important roles in combustion chemistry. At the HEAT345-(Q) level of theory, recommended enthalpies of formation at 0 K are 301.5 ± 1.3, 180.3 ± 1.8, and 23.4 ± 1.5 kJ mol(-1) for vinyl, allyl, and vinoxy, respectively. At 298 K, the corresponding values are 297.3, 168.6, and 16.1 kJ mol(-1), with the same uncertainties. The calculated values for the three radicals are in excellent agreement with the corresponding experimental values, but the uncertainties associated with the HEAT values for vinoxy are considerably smaller than those based on experimental studies.

Towards highly accurate ab initio thermochemistry of larger systems: Benzene

The high accuracy extrapolated ab initio thermochemistry (HEAT) protocol is applied to compute the total atomization energy (TAE) and the heat of formation of benzene. Large-scale coupled-cluster calculations with more than 1500 basis functions and 42 correlated electrons as well as zero-point energies based on full cubic and (semi)diagonal quartic force fields obtained with the coupled-cluster singles and doubles with perturbative treatment of the triples method and atomic natural orbital (ANO) triple- and quadruple-zeta basis sets are presented. The performance of modifications to the HEAT scheme and the scaling properties of its contributions with respect to the system size are investigated. A purely quantum-chemical TAE and associated conservative error bar of 5463.0 ± 3.1 kJ mol(-1) are obtained, while the corresponding 95% confidence interval, based on a statistical analysis of HEAT results for other and related molecules, is ± 1.8 kJ mol(-1). The heat of formation of benzene is determined to be 101.5 ± 2.0 kJ mol(-1) and 83.9 ± 2.1 kJ mol(-1) at 0 K and 298.15 K, respectively.

Tunneling in a Simple Bond Scission: The Surprising Barrier in the H Loss from HCOOH+

The dissociation dynamics of gas phase formic acid ions (HCOOH(+), DCOOD(+), HCOOD(+), DCOOH(+)) are investigated by threshold photoelectron-photoion coincidence (TPEPICO) spectroscopy and high level ab initio calculations. The slow rate constants for this seemingly simple H loss reaction and the large onset energy shifts due to isotopic substitution point to a substantial exit barrier through which the H or D atoms tunnel. Modeling of the HCOOH(+) experimental data using RRKM theory with tunneling through an Eckart potential are best fitted with a barrier of about 17 kJ mol(-1). High level ab initio calculations support the experimental findings with a computed barrier of 15.9 kJ mol(-1), which is associated with the substantial geometry change between the product HOCO(+) cation and the corresponding HCOOH(+) molecular ion. Because of this exit channel barrier, the formic acid ion dissociation does not provide a route for determination of the HOCO(+) heat of formation. Rather, the most accurate value comes from the calculations employing the high accuracy extrapolated ab initio thermochemistry (HEAT) scheme, which yields a Δ(f)H(o)(0K)[HOCO(+)] = 600.3 ± 1.0 kJ mol(-1) (Δ(f)H(o)(298K)[HOCO(+)] = 597.3 ± 1.0 kJ mol(-1)). The calculated proton affinity of CO(2) is thus 534.7 ± 1.0 kJ mol(-1) at 0 K and 539.3 ± 1.0 kJ mol(-1) at 298.15 K.

High-Accuracy Extrapolated ab Initio Thermochemistry of the Propargyl Radical and the Singlet C3H2 Carbenes

Total atomization energies and enthalpies of formation (0 and 298.15 K) are evaluated using the high-accuracy extrapolated ab initio thermochemistry (HEAT) scheme for the two stable singlet C(3)H(2) carbenes [cyclopropenylidene (c-C(3)H(2)) and propadienylidene (vinylidencarbene, l-C(3)H(2))], as well as for the 2-propynyl (propargyl, C(3)H(3)) radical. In the case of propargyl, the HEAT protocol predicts an enthalpy of formation of 354.9 +/- 1.0 kJ mol(-1) for 0 K; corresponding values of 498.1 +/- 1.0 and 555.6 +/- 1.0 kJ mol(-1) are estimated for c-C(3)H(2) and l-C(3)H(2). Additional consideration of temperature corrections leads to estimates of 352.2 +/- 1.0, 497.1 +/- 1.0, and 556.7 +/- 1.0 kJ mol(-1) for the heats of formation at 298.15 K of the propargyl radical, c-C(3)H(2), and l-C(3)H(2), respectively. Potential strategies for simplifying the HEAT protocol are also investigated and shown to have negligible impact on accuracy.

High-Accuracy Extrapolated Ab Initio Thermochemistry of Vinyl Chloride

High-Accuracy Extrapolated Ab Initio Thermochemistry of Vinyl Chloride

High-accuracy extrapolated ab initio thermochemistry. III. Additional improvements and overview

Effects of increased basis-set size as well as a correlated treatment of the diagonal Born-Oppenheimer approximation are studied within the context of the high-accuracy extrapolated ab initio thermochemistry (HEAT) theoretical model chemistry. It is found that the addition of these ostensible improvements does little to increase the overall accuracy of HEAT for the determination of molecular atomization energies. Fortuitous cancellation of high-level effects is shown to give the overall HEAT strategy an accuracy that is, in fact, higher than most of its individual components. In addition, the issue of core-valence electron correlation separation is explored; it is found that approximate additive treatments of the two effects have limitations that are significant in the realm of <1 kJ mol(-1) theoretical thermochemistry.

High-Accuracy Extrapolated Ab Initio Thermochemistry of Vinyl Chloride

Applying a modified "high accuracy extrapolated ab initio thermochemistry" (HEAT) scheme, the standard heat of formation of vinyl chloride at 0 K is computed to be 29.79 +/- 1 kJ/mol and at 298.15 K to be 20.9 +/- 2 kJ/mol, thus resolving earlier discrepancies among the available experimental values, which span a range from 21 up to 38 kJ/mol. The enthalpies of the reactions C2H4 + Cl2 --> CH2CHCl + HCl and C2H2 + HCl --> CH2CHCl at 298.15 K are determined to be -123.0 and -113.9 +/- 2 kJ/mol, respectively.

High-accuracy extrapolated ab initio thermochemistry. II. Minor improvements to the protocol and a vital simplification

The recently developed high-accuracy extrapolated ab initio thermochemistry method for theoretical thermochemistry, which is intimately related to other high-precision protocols such as the Weizmann-3 and focal-point approaches, is revisited. Some minor improvements in theoretical rigor are introduced which do not lead to any significant additional computational overhead, but are shown to have a negligible overall effect on the accuracy. In addition, the method is extended to completely treat electron correlation effects up to pentuple excitations. The use of an approximate treatment of quadruple and pentuple excitations is suggested; the former as a pragmatic approximation for standard cases and the latter when extremely high accuracy is required. For a test suite of molecules that have rather precisely known enthalpies of formation {as taken from the active thermochemical tables of Ruscic and co-workers [Lecture Notes in Computer Science, edited by M. Parashar (Springer, Berlin, 2002), Vol. 2536, pp. 25-38; J. Phys. Chem. A 108, 9979 (2004)]}, the largest deviations between theory and experiment are 0.52, -0.70, and 0.51 kJ mol(-1) for the latter three methods, respectively. Some perspective is provided on this level of accuracy, and sources of remaining systematic deficiencies in the approaches are discussed.