Correlation Energy Extrapolation by Many-Body Expansion.

Accounting for electron correlation is required for high accuracy calculations of molecular energies. The full configuration interaction (CI) approach can fully capture the electron correlation within a given basis, but it does so at a computational expense that is impractical for all but the smallest chemical systems. In this work, a new methodology is presented to approximate configuration interaction calculations at a reduced computational expense and memory requirement, namely, the correlation energy extrapolation by many-body expansion (CEEMBE). This method combines a MBE approximation of the CI energy with an extrapolated correction obtained from CI calculations using subsets of the virtual orbitals. The extrapolation approach is inspired by, and analogous to, the method of correlation energy extrapolation by intrinsic scaling. Benchmark calculations of the new method are performed on diatomic fluorine and ozone. The method consistently achieves agreement with CI calculations to within a few mhartree and often achieves agreement to within ∼1 millihartree or less, while requiring significantly less computational resources.

The transition from the open minimum to the ring minimum on the ground state and on the lowest excited state of like symmetry in ozone: A configuration interaction study.

The metastable ring structure of the ozone 1(1)A1 ground state, which theoretical calculations have shown to exist, has so far eluded experimental detection. An accurate prediction for the energy difference between this isomer and the lower open structure is therefore of interest, as is a prediction for the isomerization barrier between them, which results from interactions between the lowest two (1)A1 states. In the present work, valence correlated energies of the 1(1)A1 state and the 2(1)A1 state were calculated at the 1(1)A1 open minimum, the 1(1)A1 ring minimum, the transition state between these two minima, the minimum of the 2(1)A1 state, and the conical intersection between the two states. The geometries were determined at the full-valence multi-configuration self-consistent-field level. Configuration interaction (CI) expansions up to quadruple excitations were calculated with triple-zeta atomic basis sets. The CI expansions based on eight different reference configuration spaces were explored. To obtain some of the quadruple excitation energies, the method of Correlation Energy Extrapolation by Intrinsic Scaling was generalized to the simultaneous extrapolation for two states. This extrapolation method was shown to be very accurate. On the other hand, none of the CI expansions were found to have converged to millihartree (mh) accuracy at the quadruple excitation level. The data suggest that convergence to mh accuracy is probably attained at the sextuple excitation level. On the 1(1)A1 state, the present calculations yield the estimates of (ring minimum-open minimum) ∼45-50 mh and (transition state-open minimum) ∼85-90 mh. For the (2(1)A1-(1)A1) excitation energy, the estimate of ∼130-170 mh is found at the open minimum and 270-310 mh at the ring minimum. At the transition state, the difference (2(1)A1-(1)A1) is found to be between 1 and 10 mh. The geometry of the transition state on the 1(1)A1 surface and that of the minimum on the 2(1)A1 surface nearly coincide. More accurate predictions of the energy differences also require CI expansions to at least sextuple excitations with respect to the valence space. For every wave function considered, the omission of the correlations of the 2s oxygen orbitals, which is a widely used approximation, was found to cause errors of about ±10 mh with respect to the energy differences.

GVVPT2 multireference perturbation theory description of diatomic scandium, chromium, and manganese.

With relatively simple model spaces derived from valence bond models, a straightforward zero-order Hamiltonian, and the use of moderate-sized Dunning-type correlation consistent basis sets (cc-pVTZ, aug-cc-pVTZ, and cc-pVQZ), the second order generalized Van Vleck perturbation theory (GVVPT2) method is shown to produce potential energy curves (PECs) and spectroscopic constants close to experimental results for both ground and low-lying excited electronic states of Sc(2), Cr(2) and Mn(2). In spite of multiple quasidegeneracies (particularly for the cases of Sc(2) and Mn(2)), the GVVPT2 PECs are smooth with no discontinuities. Since these molecules have been identified as ones that widely used perturbative methods are inadequate for describing well, due to intruder state problems, unless shift parameters are introduced that can obfuscate the physics, this study suggests that the conclusion about the inadequacy of multireference perturbation theory be re-evaluated. The ground state of Sc(2) is predicted to be X(5)∑(u)(-), and its spectroscopic constants are close to the ones at the MRCISD level. Near equilibrium geometries, the 1(3)∑(u)(-) electronic state of Sc(2) is found to be less stable than the quintet ground state by 0.23 eV. The Cr(2) PEC has several features of the Rydberg-Klein-Rees (RKR) experimental curve (e.g., the pronounced shelf at elongated bond lengths), although the predicted bond length is slightly long (R(e) = 1.80 Å with cc-pVQZ compared to the experimental value of 1.68 Å). The X(1)∑(g)(+) ground state of Mn(2) is predicted to be a van der Waals molecule with a long bond length, R(e), of 3.83 Å using a cc-pVQZ basis set (experimental value = 3.40 Å) and a binding energy, D(e), of only 0.05 eV (experimental value = 0.1 eV). We obtained R(e) = 3.40 Å and D(e) = 0.09 eV at the complete basis set (CBS) limit for ground state Mn(2). Low lying excited state curves have also been characterized for all three cases (Cr(2), Mn(2), and Sc(2)) and show similar mathematical robustness as the ground states. These results suggest that the GVVPT2 multireference perturbation theory method is more broadly applicable than previously documented.

GVVPT2 energy gradient using a Lagrangian formulation.

A Lagrangian based approach was used to obtain analytic formulas for GVVPT2 energy nuclear gradients. The formalism can use either complete or incomplete model (or reference) spaces, and is limited, in this regard, only by the capabilities of the MCSCF program. An efficient means of evaluating the gradient equations is described. Demonstrative calculations were performed and compared with finite difference calculations on several molecules and show that the GVVPT2 gradients are accurate. Of particular interest, the suggested formalism can straightforwardly use state-averaged MCSCF descriptions of the reference space in which the states have arbitrary weights. This capability is demonstrated by some calculations on the ground and first excited singlet states of LiH, including calculations near an avoided crossing. The accuracy and usefulness of the GVVPT2 method and its gradient are highlighted by comparing the geometry of the near-C(2v) minimum on the conical intersection seam between the 1 (1)A(1) and 2 (1)A(1) surfaces of O(3) with values that were calculated at the multireference configuration interaction, including single and double excitations (MRCISD), level of theory.

Reactions of hydrogen atom with hydrogen peroxide.

Rate coefficients are calculated using canonical variational transition state theory with multidimensional tunneling (CVT/SCT) for the reactions H + H2O2 --> H2O + OH (1a) and H + H2O2 --> HO2 + H2 (1b). Reaction barrier heights are determined using two theoretical approaches: (i) comparison of parametrized rate coefficient calculations employing CVT/SCT to experiment and (ii) high-level ab initio methods. The evaluated experimental data reveal considerable variations of the barrier height for the first reaction: although the zero-point-exclusive barrier for (1a) derived from the data by Klemm et al. (First Int. Chem. Kinet. Symposium 1975, 61) is 4.6 kcal/mol, other available measurements result in a higher barrier of 6.2 kcal/mol. The empirically derived zero-point-exclusive barrier for (1b) is 10.4 kcal/mol. The electronic structure of the system at transition state geometries in both reactions was found to have "multireference" character; therefore special care was taken when analyzing electronic structure calculations. Transition state geometries are optimized by multireference perturbation theory (MRMP2) with a variety of one-electron basis sets, and by a multireference coupled cluster (MR-AQCCSD) method. A variety of single-reference benchmark-level calculations have also been carried out; included among them are BMC-CCSD, G3SX(MP3), G3SX, G3, G2, MCG3, CBS-APNO, CBS-Q, CBS-QB3, and CCSD(T). Our data obtained at the MRMP2 level are the most complete; the barrier height for (1a) using MRMP2 at the infinite basis set limit is 4.8 kcal/mol. Results are also obtained with midlevel single-reference multicoefficient correlation methods, such as MC3BB, MC3MPW, MC-QCISD/3, and MC-QCISD-MPWB, and with a variety of hybrid density functional methods, which are compared with high-level theory. On the basis of the evaluated experimental values and the benchmark calculations, two possible recommended values are given for the rate coefficients.

A portable automated system for trace gas sampling in the field and stable isotope analysis in the laboratory.

A computer-controllable mobile system is presented which enables the automatic collection of 33 air samples in the field and the subsequent analysis for delta13C and delta18O stable isotope ratios of a carbon-containing trace gas in the laboratory, e.g. CO2, CO or CH4. The system includes a manifold gas source input for profile sampling and an infrared gas analyzer for in situ CO2 concentration measurements. Measurements of delta13C and delta18O of all 33 samples can run unattended and take less than six hours for CO2. Laboratory tests with three gases (compressed air with different pCO2 and stable isotope compositions) showed a measurement precision of 0.03 per thousand for delta13C and 0.02 per thousand for delta18O of CO2 (standard error (SE), n = 11). A field test of our system, in which 66 air samples were collected within a 24-hour period above grassland, showed a correlation of 0.99 (r2) between the inverse of pCO2 and delta13C of CO2. Storage of samples until analysis is possible for about 1 week; this can be an important factor for sampling in remote areas. A wider range of applications in the field is open with our system, since sampling and analysis of CO and CH4 for stable isotope composition is also possible. Samples of compressed air had a measurement precision (SE, n = 33) of 0.03 per thousand for delta13C and of 0.04 per thousand for delta18O on CO and of 0.07 per thousand for delta13C on CH4. Our system should therefore further facilitate research of trace gases in the context of the carbon cycle in the field, and opens many other possible applications with carbon- and possibly non-carbon-containing trace gases.