Search

We present an efficient, open-source formulation for coupled-cluster theory through perturbative triples with domain-based local pair natural orbitals [DLPNO-CCSD(T)]. Similar to the implementation of the DLPNO-CCSD(T) method found in the ORCA package, the most expensive integral generation and contraction steps associated with the CCSD(T) method are linear-scaling. In this work, we show that the t1-transformed Hamiltonian allows for a less complex algorithm when evaluating the local CCSD(T) energy without compromising efficiency or accuracy. Our algorithm yields sub-kJ mol−1 deviations for relative energies when compared with canonical CCSD(T), with typical errors being on the order of 0.1 kcal mol−1, using our TightPNO parameters. We extensively tested and optimized our algorithm and parameters for non-covalent interactions, which have been the most difficult interaction to model for orbital (PNO)-based methods historically. To highlight the capabilities of our code, we tested it on large water clusters, as well as insulin (787 atoms). [ABSTRACT FROM AUTHOR]

A consideration of the point group symmetry of molecules is often advantageous from a computational efficiency standpoint and sometimes necessary for the correct treatment of chemical physics problems. Many modern electronic structure software packages include a treatment of symmetry, but these are sometimes incomplete or unusable outside of that program's environment. Therefore, we have developed the MolSym package for handling molecular symmetry and its associated functionalities to provide a platform for including symmetry in the implementation and development of other methods. Features include point group detection, molecule symmetrization, arbitrary generation of symmetry element sets and character tables, and symmetry adapted linear combinations of real spherical harmonic basis functions, Cartesian displacement coordinates, and internal coordinates. We present some of the advantages of using molecular symmetry as achieved by MolSym, particularly with respect to Hartree–Fock theory, and the reduction of finite difference displacements in gradient/Hessian computations. This package is designed to be easily integrated into other software development efforts and may be extended to further symmetry applications. [ABSTRACT FROM AUTHOR]

Over the past half century, computational chemistry has evolved from a niche field to a ubiquitous pillar of modern chemical research. Driven by the increased demand for computational chemistry in research settings, the undergraduate curriculum has evolved alongside to ensure that students are well-equipped for modern research. Toward this end, many excellent computational chemistry exercises have been developed that aim to teach students what kinds of questions computational chemistry can answer and how to properly interpret the results. However, there has been far less attention given to the complexities of determining how reliable computational results are and how constraining computational scaling can be. We present an undergraduate lab exercise that uses ab initio methods to predict the combustion energy of methane. The exercise walks students through the process of benchmarking errors on a small system (methane), estimating the computational cost to perform the same analysis on a larger system (propane), and justifying an affordable yet accurate method for a hypothetical study of the larger system. Furthermore, students are introduced to other cost-saving measures like basis set extrapolation and additive corrections. The entire exercise is intentionally designed to require little technical knowledge of computational chemistry and to be flexibly grafted into a standard undergraduate curriculum. In order to ensure accessibility, the exercise utilizes the free open source software Psi4 (available on any operating system) and provides a detailed installation and use guide for completing this lab. This lab will provide students the understanding of how to properly judge, select, and justify different computational models where cost and accuracy compete, a highly desirable set of skills that generalize to any computational science.

Multifidelity modeling is a technique for fusing the information from two or more datasets into one model. It is particularly advantageous when one dataset contains few accurate results and the other contains many less accurate results. Within the context of modeling potential energy surfaces, the low-fidelity dataset can be made up of a large number of inexpensive energy computations that provide adequate coverage of the N-dimensional space spanned by the molecular internal coordinates. The high-fidelity dataset can provide fewer but more accurate electronic energies for the molecule in question. Here, we compare the performance of several neural network-based approaches to multifidelity modeling. We show that the four methods (dual, Δ-learning, weight transfer, and Meng–Karniadakis neural networks) outperform a traditional implementation of a neural network, given the same amount of training data. We also show that the Δ-learning approach is the most practical and tends to provide the most accurate model. [ABSTRACT FROM AUTHOR]

Lead is heavily utilised in industry, but it bears potential harm due to its toxic nature. Identifying lead is important in diagnosing and preventing lead poisoning. Hydride generation via inductively-coupled plasma mass spectrometer (HG-ICP-MS) can detect small amounts of lead with high accuracy. Here we investigate seven structures of Pb $ _2 $ 2 H $ _2^+ $ 2 + . The geometries of each structure were optimised using coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] method with the cc-pwCVQZ-PP basis set for lead and the cc-pVQZ basis set for hydrogen. Likewise, harmonic vibrational frequencies were obtained using the same method. Anharmonic frequencies were determined using CCSD(T) with cc-pwCVTZ-PP and cc-pVTZ basis sets for lead and hydrogen, respectively. The focal point approach (FPA), including CCSDT and CCSDT(Q) computations, was applied to the seven structures to find the relative energies with respect to the global minimum butterfly structure. Relative to the global minimum butterfly structure, the energies for the planar transition state, monobridged, vinylidene-like, trans, cis, and linear structures are: 1.1, 18.3, 26.5, 27.8, 28.8, and 93.7 kcal mol $ ^{-1} $ − 1 , respectively. Natural bond orders and dipole moments are reported as well. The predicted vibrational frequencies should be helpful in guiding the laboratory observation of Pb $ _2 $ 2 H $ _2^+ $ 2 + . [ABSTRACT FROM AUTHOR]

Criegee intermediates, formed from the ozonolysis of alkenes, are known to have a role in atmospheric chemistry, including the modulation of the oxidizing capacity of the troposphere. Although studies have been conducted since their discovery, the synthesis of these species in the laboratory has ushered in a new wave of investigations of these structures, both theoretically and experimentally. In some of these theoretical studies, high-order corrections for correlation energy are included to account for the mid multi-reference character found in these systems. Many of these studies include a focus on kinetics; therefore, the calculated energies should be accurate (<1 kcal/mol in error). In this research, we compute the enthalpies of formation for a small set of Criegee intermediates, including higher-order coupled cluster corrections for correlation energy up to coupled cluster with perturbative quintuple excitations. The enthalpies of formation for formaldehyde oxide, anti-acetaldehyde oxide, syn-acetaldehyde oxide, and acetone oxide are presented at 0 K as 26.5, 15.6, 12.2, and 0.1 kcal mol−1, respectively. Additionally, we do not recommend the coupled cluster with perturbative quadruple excitations [CCSDT(Q)] energy correction, as it is approximately twice as large as that of the coupled cluster with full quadruple excitations (CCSDTQ). Half of the CCSDT(Q) energy correction may be included as a reliable, cost-effective estimation of CCSDTQ energies for Criegee intermediates. [ABSTRACT FROM AUTHOR]

Thallium chemistry is experiencing unprecedented importance. Therefore, it is valuable to characterize some of the simplest thallium compounds. Stationary points along the singlet and triplet Tl 2H 2 potential energy surface have been characterized. Stationary point geometries were optimized with the CCSD(T)/aug‐cc‐pwCVQZ‐PP method. Harmonic vibrational frequencies were computed at the same level of theory while anharmonic vibrational frequencies were computed at the CCSD(T)/aug‐cc‐pwCVTZ‐PP level of theory. Final energetics were obtained with the CCSDT(Q) method. Basis sets up to augmented quintuple‐zeta cardinality (aug‐cc‐pwCV5Z‐PP) were employed to obtain energetics in order to extrapolate to the complete basis set limits using the focal point approach. Zero‐point vibrational energy corrections were appended to the extrapolated energies in order to determine relative energies at 0 K. It was found that the planar dibridged isomer lies lowest in energy while the linear structure lies highest in energy. The results were compared to other group 13 M 2H 2 (M = B, Al, Ga, In, and Tl) theoretical studies and some interesting variations are found. With respect to experiment, incompatibilities exist. [ABSTRACT FROM AUTHOR]

Five structures of Ge2H2 and Ge2H2+ are investigated in this study. Optimized geometries at the CCSD(T)/cc-pwCVQZ-PP level of theory were obtained. Focal point analyses were performed on these optimized geometries to determine relative energies using the CCSD(T) method with polarized basis sets up to quintuple-zeta. Energy corrections include full T and pertubative (Q) coupled-cluster effects plus anharmonic corrections to the zero-point vibrational energy. Relative ordering in energy from lowest to highest of the five Ge2H2+ structures is butterfly, germylidene, monobridged, trans, then linear. In neutral Ge2H2, the monobridged structure lies lower in energy than the germylidene structure. Fundamental vibrational frequencies and IR intensities were computed for the minima at the CCSD(T)/cc-pwCVTZ-PP level of theory to compare with experimental research. Partial atomic charges and natural bonding orbital analyses indicated that the positive charge of Ge2H2+ is contained in the region of the Ge–Ge bond. [ABSTRACT FROM AUTHOR]

In 2021 Olaru, Mebs, and Beckmann reported the synthesis of remarkable cationic carbene analogues $ \mathrm{PR_2^+} $ P R 2 + and $ \mathrm{AsR_2^+} $ As R 2 + . This work followed the same group's synthesis of $ \mathrm{SbR_2^+} $ Sb R 2 + and $ \mathrm{BiR_2^+} $ Bi R 2 + . To better understand these important systems, $ \mathrm{SbH_2^+} $ Sb H 2 + and $ \mathrm{BiH_2^+} $ Bi H 2 + have been studied with high level ab initio quantum mechanical methods. Geometries were optimised with the CCSDT(Q) method with the cc-pwCVTZ-PP basis set using small core pseudopotentials. Fundamental vibrational frequencies were computed to provide theoretical predictions for future synthetic studies. Relative energies with respect to $ \mathrm{Pn^+} $ P n + + $ \mathrm{H_2} $ H 2 ( $ \mathrm{Pn} $ Pn = $ \mathrm{Sb} $ Sb , $ \mathrm{Bi} $ Bi ) were determined at the CCSDT/CBS level of theory via the Focal Point Analysis method, and anharmonic zero-point vibrational energy and higher order contributions were also computed. For $ \mathrm{SbH_2^+} $ Sb H 2 + , we obtained R = 1.697 Å and $ \theta =90.8^{\circ } $ θ = 90.8 ∘ with $ \mathrm{Sb^+ + H_2 \rightarrow SbH_2^+} $ S b + + H 2 → Sb H 2 + reaction enthalpy of $ \Delta H = -18.71\,{\rm kcal}\ {\rm mol}^{-1} $ ΔH = − 18.71 kcal mol − 1 . For $ \mathrm{BiH_2^+} $ Bi H 2 + , the analogous results were R = 1.774 Å, $ \theta =89.7^{\circ } $ θ = 89.7 ∘ , and $ \Delta H = -7.64\,{\rm kcal}\,{\rm mol}^{-1} $ ΔH = − 7.64 kcal mol − 1 for $ \mathrm{Bi^+ + H_2 \rightarrow BiH_2^+} $ B i + + H 2 → Bi H 2 + . [ABSTRACT FROM AUTHOR]

Preceding even the Hartree-Fock method, molecular integrals are the very foundation upon which quantum chemical molecular modeling depends. Discussions of molecular integrals are normally found only in advanced and technical texts or articles. The objective of the present article is to provide less experienced readers, or students in a physical/computational chemistry course, a thorough understanding of molecular integrals. Through a series of detailed Handouts, the student/reader can participate in the derivation of molecular integrals, and in turn implement them in computer code. Hartree-Fock theory is discussed in enough detail to motivate the molecular integrals and address such topics as the atomic orbital basis. An introduction to the programming language of choice, Python3, is provided, tailored toward developing the essential skills necessary for implementing molecular integrals. The article is intended to be useful not only to instructors of physical/computational chemistry, but also to any reader who has independently sought a primer on this elusive subject.

We study the combination of orbital-optimized density cumulant theory and a new parameterization of reduced density matrices in which the variables are the particle–hole cumulant elements. We call this combination OλDCT. We find that this new Ansatz solves problems identified in the previous unitary coupled cluster Ansatz for density cumulant theory: the theory is now free of near-zero denominators between occupied and virtual blocks, can correctly describe the dissociation of H2, and is rigorously size-extensive. In addition, the new Ansatz has fewer terms than the previous unitary Ansatz, and the optimal orbitals delivered by the exact theory are the natural orbitals. Numerical studies on systems amenable to full configuration interaction show that the amplitudes from the previous ODC-12 method approximate the exact amplitudes predicted by this Ansatz. Studies on equilibrium properties of diatomic molecules show that even with the new Ansatz, it is necessary to include triples to improve the accuracy of the method compared to orbital-optimized linearized coupled cluster doubles. With a simple iterative triples correction, OλDCT outperforms other orbital-optimized methods truncated at comparable levels in the amplitudes, as well as coupled cluster single and doubles with perturbative triples [CCSD(T)]. By adding four more terms to the cumulant parameterization, OλDCT outperforms CCSDT while having the same O ( V 5 O 3 ) scaling. [ABSTRACT FROM AUTHOR]

The previously proposed Ansatz for density cumulant theory that combines orbital-optimization and a parameterization of the 2-electron reduced density matrix cumulant in terms of unitary coupled cluster amplitudes (OUDCT) is carefully examined. Formally, we elucidate the relationship between OUDCT and orbital-optimized unitary coupled cluster theory and show the existence of near-zero denominators in the stationarity conditions for both the exact and some approximate OUDCT methods. We implement methods of the OUDCT Ansatz restricted to double excitations for numerical study, up to the fifth commutator in the Baker–Campbell–Hausdorff expansion. We find that methods derived from the Ansatz beyond the previously known ODC-12 method tend to be less accurate for equilibrium properties and less reliable when attempting to describe H2 dissociation. New developments are needed to formulate more accurate density cumulant theory variants. [ABSTRACT FROM AUTHOR]

The acetonyl radical (•CH2COCH3) is relevant to atmospheric and combustion chemistry due to its prevalence in many important reaction mechanisms. One such reaction mechanism is the decomposition of Criegee intermediates in the atmosphere that can produce acetonyl radical and OH. In order to understand the fate of the acetonyl radical in these environments and to create more accurate kinetics models, we have examined the reaction system of the acetonyl radical with O2 using highly reliable theoretical methods. Structures were optimized using coupled cluster theory with singles, doubles, and perturbative triples [CCSD(T)] with an atomic natural orbital (ANO0) basis set. Energetics were computed to chemical accuracy using the focal point approach involving perturbative treatment of quadruple excitations [CCSDT(Q)] and basis sets as large as cc-pV5Z. The addition of O2 to the acetonyl radical produces the acetonylperoxy radical, and multireference computations on this reaction suggest it to be barrierless. No submerged pathways were found for the unimolecular isomerization of the acetonylperoxy radical. Besides dissociation to reactants, the lowest energy pathway available for the acetonylperoxy radical is a 1-5 H shift from the methyl group to the peroxy group through a transition state that is 3.3 kcal mol−1 higher in energy than acetonyl radical + O2. The ultimate products from this pathway are the enol tautomer of the acetonyl radical along with O2. Multiple pathways that lead to OH formation are considered; however, all of these pathways are predicted to be energetically inaccessible, except at high temperatures. [ABSTRACT FROM AUTHOR]

Sulfurous acid (H2SO3) is an infamously elusive molecule. Although some theoretical papers have supposed possible roles for it in more complicated systems, it has yet to be experimentally observed. To aid experiment in detecting this molecule, we have examined the H2O + SO2 potential energy surface at the CCSDT(Q)/CBS//CCSD(T)-F12b/cc-pVTZ-F12b level of theory to resolve standing discrepancies in previous reports and predict the gas-phase vibrational spectrum for H2SO3. We find that sulfurous acid has two potentially detectable rotamers, separated by 1.1 kcal mol−1 ΔH0K with a torsional barrier of 1.6 kcal mol−1. The sulfonic acid isomer is only 6.9 kcal mol−1 above the lowest enthalpy sulfurous acid rotamer, but the barrier to form it is 57.2 kcal mol−1. Error in previous reports can be attributed to misidentified stationary points, the use of density functionals that perform poorly for this system, and, most importantly, the basis set sensitivity of sulfur. Using VPT2+K, we determine that the intense S=O stretch fundamental of each species is separated from other intense peaks by at least 25 cm−1, providing a target for identification by infrared spectroscopy. [ABSTRACT FROM AUTHOR]

The reaction between two OH radicals has been studied in the context of combustion chemistry, atmospheric chemistry, and astrochemistry. The global minimum on the triplet reaction surface is the van der Waals complex between water and O( 3 P) atom. We optimised the geometry of the system using ab initio methods up to CCSDT/aug-cc-pVQZ. We use the focal point approach to determine an accurate dissociation energy of 1.35 kcal mol − 1 for the complex dissociating into water and O( 3 P). We also predict the anharmonic vibrational frequencies of the complex in anticipation of future experiments aimed at detecting this system. [ABSTRACT FROM AUTHOR]

Vinoxy radicals are involved in numerous atmospheric and combustion mechanisms. High-level theoretical methods have recently shed new light on the reaction of the unsubstituted vinoxy radical with O2. The reactions of 1-methylvinoxy radical and 2-methylvinoxy radical with molecular oxygen have experimental high pressure limiting rate constants, k∞, 5–7 times higher than that of the vinoxy plus O2 reaction. In this work, high-level ab initio quantum chemical computations are applied to the 2-methylvinoxy radical plus O2 system, namely, the formation and isomerization of the 1-oxo-2-propylperoxy radical, the immediate product of O2 addition to the 2-methylvinoxy radical. Multireference methods were applied to the entrance channel. No barrier to O2 addition could be located, and more sophisticated treatment of dynamic electron correlation shows that the principal difference between O2 addition to the vinoxy and 2-methylvinoxy radicals is a larger steric factor for 2-methylvinoxy + O2. This is attributed to the favorable interaction between the incoming O2 molecule and the methyl group of the 2-methylvinoxy radical. Via the focal point approach, energetics for this reaction were determined, in most cases, to chemical accuracy. The coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] correlation energy and Hartree–Fock energies were independently extrapolated to the complete basis set limit. A correction for the effect of higher excitations was computed at the CCSDT(Q)/6-31G level. Corrections for the frozen-core approximation, the Born–Oppenheimer approximation, the nonrelativistic approximation, and the zero-point vibrational energy were included. From the 1-oxo-2-propylperoxy radical, dissociation to reactants is competitive with the lowest energy isomerization pathway. The lowest energy isomerization pathway ultimately forms acetaldehyde, CO, and ·OH as the final products. [ABSTRACT FROM AUTHOR]

Alkyl combustion models depend on kinetic parameters derived from reliable experimental or theoretical energetics that are often unavailable for larger species. To this end, we have performed a comprehensive investigation of the tert-butyl radical (R· in this paper) autoignition pathways. CCSD(T)/ANO0 geometries and harmonic vibrational frequencies were obtained for key stationary points for the R· + O2 and QOOH + O2 mechanisms. Relative energies were computed to chemical accuracy (±1 kcal mol-1) via extrapolation of RCCSD(T) energies to the complete basis-set limit, or usage of RCCSD(T)-F12 methods. At 0 K, the minimum energy R· +O2 pathway involves direct elimination of HO. 2 (30.3 kcal mol-1 barrier) from the tert-butyl peroxy radical (ROO·) to give isobutene. This pathway lies well below the competing QOOH-forming intramolecular hydrogen abstraction pathway (36.2 kcal mol-1 barrier) and ROO· dissociation (35.9 kcal mol-1 barrier). The most favorable decomposition channel for QOOH radicals leads to isobutene oxide (12.0 kcal mol-1 barrier) over isobutene (18.6 kcal mol-1 barrier). For the QOOH + O2 pathways, we studied the transition states and initial products along three pathways: (1) α-hydrogen abstraction (42.0 kcal mol-1 barrier), (2) γ-hydrogen abstraction (27.0 kcal mol-1 barrier), and (3) hydrogen transfer to the peroxy moiety (24.4 kcal mol-1 barrier). The barrier is an extensive modification to the previous 18.7 kcal mol-1 value andwarrants further study. However, it is still likely that the lowest energy QOOH+O2 pathway corresponds to pathway (3).We found significant spin contamination and/or multireference character in multiple stationary points, especially for transition states stemming from QOOH. Lastly, we provide evidence for an ... surface crossing at a Cs-symmetric, intramolecular hydrogen abstraction structure. [ABSTRACT FROM AUTHOR]

The n-propylperoxy radical has been described as a molecule of critical importance to studies of low temperature combustion. Ab initio methods were used to study this three-carbon alkylperoxy radical, normal propylperoxy. Reliable CCSD(T) (coupled-cluster theory, incorporating single, double, and perturbative triple)/ANO0 geometries were predicted for the molecule's five rotamers. For each rotamer, energetic predictions were made using basis sets as large as the cc-pV5Z in conjunction with coupled cluster levels of theory up to CCSDT(Q). Along with the extrapolations, corrections for relativistic effects, zero-point vibrational energies, and diagonal Born-Oppenheimer corrections were used to further refine energies. The results indicate that the lowest conformer is the gauche-gauche (GG) rotamer followed by the gauche-trans (0.12 kcal mol-1 above GG), trans-gauche (0.44 kcal mol-1), gauche'-gauche (0.47 kcal mol-1), and trans-trans (0.57 kcal mol-1). Fundamental vibrational frequencies were obtained using second-order vibrational perturbation theory. This is the first time anharmonic frequencies have been computed for this system. The most intense IR features include all but one of the C-H stretches. The O-O fundamental (1063 cm-1 for the GG structure) also has a significant IR intensity, 19.6 km mol-1. The anharmonicity effects on the potential energy surface were also used to compute vibrationally averaged rg,0K bond lengths, accounting for zero-point vibrations present within the molecule. [ABSTRACT FROM AUTHOR]

Four constitutional isomers of diindium dihydride have been studied utilising rigorous quantum mechanical methods. Geometries were optimised with the CCSD(T) method using the aug-cc-pwCVTZ basis set and a small-core pseudopotential for the indium atoms. Relative energetics were determined at the CCSD(T)/CBS level of theory, and the higher order δT(Q) contributions are also computed. The monobridged and vinylidene-like isomers lie 10.8 and 13.6 kcal mol − 1 above the planar dibridged isomer, respectively. At 0 K, the trans isomer is the least favoured energetically, differing from the dibridged isomer by 16.6 kcal mol − 1 . Anharmonic vibrational frequencies were predicted to provide theoretical insight into future experimental studies. Very detailed comparisons are made to matrix isolation experiments from the group of Downs and the group of Andrews. Numerous agreements and disagreements were found. [ABSTRACT FROM AUTHOR]

Pnictinidenes are an increasingly relevant species in main group chemistry and generally exhibit proclivity for the triplet electronic ground state. However, the elusive singlet electronic states are often desired for chemical applications. We predict the singlet‐triplet energy differences (ΔEST=ESinglet−ETriplet) of simple group 15 and 16 substituted pnictinidenes (Pn−R; Pn=P, As, Sb, or Bi) with highly reliable focal‐point analyses targeting the CCSDTQ/CBS level of theory. The only cases we predict to have favorable singlet states are P−PH2 (−3.2 kcal mol−1) and P−NH2 (−0.2 kcal mol−1). ΔEST trends are discussed in light of the geometric predictions as well as qualitative natural bond order analysis to elucidate some of the important electronic structure features. Our work provides a rigorous benchmark for the ΔEST of fundamental Pn−R moieties and provides a firm foundation for the continued study of heavier pnictinidenes. [ABSTRACT FROM AUTHOR]

The cyano radical (CN) is an abundant, open-shell molecule found in a variety of environments, including the atmosphere, the interstellar medium and combustion processes. In these environments, it often reacts with small, closed-shell molecules via hydrogen abstraction. Both carbon and nitrogen atoms of the cyano radical are reactive sites, however the carbon is more reactive with reaction barrier heights generally between 2–15 kcal mol−1 lower than those of the analogous nitrogen. The CN + HX → HCN/HNC + X, with X = H, CH3, NH2, OH, F, SiH3, PH2, SH, Cl, C2H, CN reactions have been studied at a high-level of theory, including CCSD(T)-F12a. Furthermore, kinetics were obtained over the 100–1000 K temperature range, showing excellent agreement with those rate constants that have been determined experimentally. [ABSTRACT FROM AUTHOR]

Oxiranes are a class of cyclic ethers formed in abundance during low‐temperature combustion of hydrocarbons and biofuels, either via chain‐propagating steps that occur from unimolecular decomposition of β‐hydroperoxyalkyl radicals (β‐̇QOOH) or from reactions of HȮO with alkenes. The cis‐ and trans‐isomers of 2,3‐dimethyloxirane are intermediates of n‐butane oxidation, and while rate coefficients for β‐̇QOOH → 2,3‐dimethyloxirane + ̇OH are reported extensively, subsequent reaction mechanisms of the cyclic ethers are not. As a result, chemical kinetics mechanisms commonly adopt simplified chemistry to describe the consumption of 2,3‐dimethyloxirane by convoluting several elementary reactions into a single step, which may introduce mechanism truncation error—uncertainty derived from missing or incomplete chemistry. The present research examines the isomer dependence of 2,3‐dimethyloxirane reaction mechanisms in support of ongoing efforts to minimize mechanism truncation error. Reaction mechanisms are inferred via the detection of products from Cl‐initiated oxidation of both cis‐2,3‐dimethyloxirane and trans‐2,3‐dimethyloxirane using multiplexed photoionization mass spectrometry (MPIMS). The experiments were conducted at 10 Torr and temperatures of 650 K and 800 K. To complement the experiments, the enthalpies of stationary points on the ̇R + O2 surfaces were computed at the ccCA‐PS3 level of theory. In total, 28 barrier heights were computed on the 2,3‐dimethyloxiranylperoxy surfaces. Two notable aspects are low‐lying pathways that form resonance‐stabilized ketohydroperoxide‐type radicals caused by ̇QOOH ring‐opening when the unpaired electron is localized adjacent to the ether group, and cis‐trans isomerization of ̇R and ̇QOOH radicals, via inversion, which enable reaction pathways otherwise restricted by stereochemistry. Several species were identified in the MPIMS experiments from ring opening of 2,3‐dimethyloxiranyl radicals. Neither of the two conjugate alkene isomers prototypical of ̇R + O2 reactions were detected. Products were also identified from decomposition of ketohydroperoxide‐type radicals. The present work provides the first analysis of 2,3‐dimethyloxirane oxidation chemistry and reveals that consumption pathways are complex and require the expansion of submechanisms in chemical kinetics mechanisms. [ABSTRACT FROM AUTHOR]

Oxiranes are a class of cyclic ethers formed in abundance during low‐temperature combustion of hydrocarbons and biofuels, either via chain‐propagating steps that occur from unimolecular decomposition of β‐hydroperoxyalkyl radicals (β‐̇QOOH) or from reactions of HOȮ with alkenes. Ethyloxirane is one of four alkyl‐substituted cyclic ether isomers produced as an intermediate from n‐butane oxidation. While rate coefficients for β‐̇QOOH → ethyloxirane + ȮH are reported extensively, subsequent reaction mechanisms of the cyclic ether are not. As a result, chemical kinetics mechanisms commonly adopt simplified chemistry to describe ethyloxirane consumption by convoluting several elementary reactions into a single step, which may introduce mechanism truncation error—uncertainty derived from missing or incomplete chemistry. The present work provides fundamental insight on reaction mechanisms of ethyloxirane in support of ongoing efforts to minimize mechanism truncation error. Reaction mechanisms are inferred from the detection of products during chlorine atom‐initiated oxidation experiments using multiplexed photoionization mass spectrometry conducted at 10 Torr and temperatures of 650 K and 800 K. To complement the experiments, calculations of stationary point energies were conducted using the ccCA‐PS3 composite method on ̇R + O2 potential energy surfaces for the four ethyloxiranyl radical isomers, which produced barrier heights for 24 reaction pathways. In addition to products from ̇QOOH → cyclic ether + ȮH and ̇R + O2 → conjugate alkene + HOȮ, both of which were significant pathways and are prototypical to alkane oxidation, other species were identified from ring‐opening of both ethyloxiranyl and ̇QOOH radicals. The latter occurs when the unpaired electron is localized on the ether group, causing the initial ̇QOOH structure to ring‐open and form a resonance‐stabilized ketohydroperoxide‐type radical. The present work provides the first analysis of ethyloxirane oxidation chemistry, which reveals that consumption pathways are complex and may require an expansion of submechanisms to increase the fidelity of chemical kinetics mechanisms. [ABSTRACT FROM AUTHOR]

Formic acid (FA) is a ubiquitous molecule found in the atmosphere, and is relevant to many important processes. The FA molecule generally exists as the trans isomer, which can decompose into H2O and CO (dehydration). It can also exist in the less favorable cis isomer which can decompose into H2 and CO2 (decarboxylation). Our work examines the complexes formed between each isomer of FA with water. We present geometries and vibrational frequencies obtained at the reliable CCSD(T)/aug-cc-pVTZ level of theory for seven FA⋯water complexes. We utilize the focal point method to determine CCSDT(Q)/CBS plus corrections binding energies of 7.37, 3.36, and 2.02 kcal mol−1 plus 6.07, 3.79, 2.60, and 2.55 kcal mol−1 for the trans-FA⋯water and cis-FA⋯water complexes, respectively. Natural bond orbital analysis is used to further decompose the interactions in each complex and gain insight into their relative strengths. Furthermore, we examine the effect that a single water molecule has on the barrier heights to each decomposition pathway by optimizing the transition states and verifying their connectivity with intrinsic reaction coordinate computations as well as utilizing a kinetic model. Water lowers the barrier to dehydration by at most 15.78 kcal mol−1 and the barrier to decarboxylation by up to 15.90 kcal mol−1. Our research also examines for the first time the effect of one water molecule on the interconversion barrier and we find that the barrier from trans to cis is not catalyzed by water due to the strong FA and water interactions. Our results highlight some instances where different binary complexes result in different decomposition pathways and even a case where one binary complex can form the same decomposition products via two distinct mechanisms. Our results provide a reliable benchmark of the FA⋯H2O system as well as provide insight into future studies of similar atmospheric systems. [ABSTRACT FROM AUTHOR]

Criegee intermediates are important targets for study in atmospheric chemistry because of their capacity to oxidize airborne species. Among these species, ammonia has received critical attention for its presence in polluted agricultural or industrial areas and its role in forming particulate matter and condensation nuclei. Although methylamine has been given less attention than ammonia, both theoretical and experimental studies have demonstrated that the additional methyl substitution on the ammonia derivatives increases the rate constants for some systems. This suggests that the methylamine addition to Criegee intermediates could be more significant to atmospheric processes. In this work, geometries are optimized at the DF-CCSD(T)/ANO1 level for the methylamine addition reactions to the simplest Criegee intermediate and the anti- and syn-methylated Criegee intermediates. Energies for each stationary point were computed at the CCSD(T)/CBS level with corrections from the CCSDT(Q) method. Rate constants are obtained for each reaction using canonical transition state theory. Although methylamine addition proved to be a more favorable reaction relative to ammonia addition, the significantly lower concentration of atmospheric methylamine limits the prevalence of these reactions, even in the most optimal conditions. It is unlikely that the methylamine addition to Criegee intermediates will contribute significantly to the consumption of Criegee intermediates in the atmosphere. [ABSTRACT FROM AUTHOR]

The ethynyl radical (C 2 H) is ubiquitous, appearing in a variety of environments from interstellar space to oxy-acetylene flames. Under these diverse conditions, ethynyl has a high affinity to abstract hydrogen atoms from nearby molecules. In this study, C 2 H + HX → C 2 H 2 + X , X = CH 3 , NH 2 , OH, F, PH 2 , SH, Cl, C 2 H, CN, and NC reactions have been examined at a rigorously high level of theory. Various additive corrections have been appended to the coupled-cluster framework to converge the relative enthalpy of stationary points on the C 2 H + HX surfaces to subchemical accuracy (≤0.5 kcal mol − 1 ). The barriers for the abstraction of a hydrogen from NH 3 , PH 3 , and H 2 S are submerged below the relative enthalpy of their reactant asymptotes by 2.75, 2.48, and 1.39 kcal mol − 1 , respectively. Abstractions from H 2 , CH 4 , H 2 O, and HCl are predicted to have moderate barriers to the abstraction of 2.64, 0.77, 3.82, and 2.19 kcal mol − 1 , respectively, suggesting these reactions will be significant to various systems at moderate to high temperatures. Accurate kinetics were obtained for the C 2 H + H 2 , CH 4 , NH 3 , and H 2 O reactions over a wide range of temperatures with excellent agreement to present experimental rate constants. [ABSTRACT FROM AUTHOR]

The ethyl radical has been isolated and spectroscopically characterized in 4He nanodroplets. The band origins of the five CH stretch fundamentals are shifted by < 2 cm-1 from those reported for the gas phase species [S. Davis, D. Uy, and D. J. Nesbitt, J. Chem. Phys. 112, 1823 (2000); T. Häber, A. C. Blair, D. J. Nesbitt, and M. D. Schuder, J. Chem. Phys. 124, 054316 (2006)]. The symmetric CH2 stretching band (v1) is rotationally resolved, revealing nuclear spin statistical weights predicted by G12 permutation-inversion group theory. A permanent electric dipole moment of 0.28 (2) D is obtained via the Stark spectrum of the v1 band. The four other CH stretch fundamental bands are significantly broadened in He droplets and lack rotational fine structure. This broadening is attributed to symmetry dependent vibration-to-vibration relaxation facilitated by the He droplet environment. In addition to the five fundamentals, three a1′ overtone/combination bands are observed, and each of these have resolved rotational substructure. These are assigned to the 2v12, v4 + v6, and 2v6 bands through comparisons to anharmonic frequency computations at the CCSD(T)/cc-pVTZ level of theory. [ABSTRACT FROM AUTHOR]

This research presents an ab initio characterization of the potential energy surface for the methylamine plus 1D oxygen atom reaction, which may be relevant to interstellar chemistry. Geometries and harmonic vibrational frequencies were determined for all stationary points at the CCSD(T)/aug-cc-pVTZ level of theory. The focal point method along with several additive corrections was used to obtain reliable CCSDT(Q)/CBS potential energy surface features. Extensive conformational analysis and intrinsic reaction coordinate computations were performed to ensure accurate chemical connectivity of the stationary points. Five minima were determined to be possible products of this reaction and three novel transition states were found that were previously unreported or mislabeled in the literature. The pathways we present can be used to guide further searches for NH2 containing species in the interstellar medium. [ABSTRACT FROM AUTHOR]

Symmetry-adapted perturbation theory (SAPT) provides a means of probing the fundamental nature of intermolecular interactions. Low-orders of SAPT (here, SAPT0) are especially attractive since they provide qualitative (sometimes quantitative) results while remaining tractable for large systems. The application of density fitting and Laplace transformation techniques to SAPT0 can significantly reduce the expense associated with these computations and make even larger systems accessible. We present new factorizations of the SAPT0 equations with density-fitted two-electron integrals and the first application of Laplace transformations of energy denominators to SAPT. The improved scalability of the DF-SAPT0 implementation allows it to be applied to systems with more than 200 atoms and 2800 basis functions. The Laplace-transformed energy denominators are compared to analogous partial Cholesky decompositions of the energy denominator tensor. Application of our new DF-SAPT0 program to the intercalation of DNA by proflavine has allowed us to determine the nature of the proflavine-DNA interaction. Overall, the proflavine-DNA interaction contains important contributions from both electrostatics and dispersion. The energetics of the intercalator interaction are are dominated by the stacking interactions (two-thirds of the total), but contain important contributions from the intercalator-backbone interactions. It is hypothesized that the geometry of the complex will be determined by the interactions of the intercalator with the backbone, because by shifting toward one side of the backbone, the intercalator can form two long hydrogen-bonding type interactions. The long-range interactions between the intercalator and the next-nearest base pairs appear to be negligible, justifying the use of truncated DNA models in computational studies of intercalation interaction energies. [ABSTRACT FROM AUTHOR]

Using a Lagrangian-based approach, we present a more elegant derivation of the equations necessary for the variational optimization of the molecular orbitals (MOs) for the coupled-cluster doubles (CCD) method and second-order Mo\ller-Plesset perturbation theory (MP2). These orbital-optimized theories are referred to as OO-CCD and OO-MP2 (or simply 'OD' and 'OMP2' for short), respectively. We also present an improved algorithm for orbital optimization in these methods. Explicit equations for response density matrices, the MO gradient, and the MO Hessian are reported both in spin-orbital and closed-shell spin-adapted forms. The Newton-Raphson algorithm is used for the optimization procedure using the MO gradient and Hessian. Further, orbital stability analyses are also carried out at correlated levels. The OD and OMP2 approaches are compared with the standard MP2, CCD, CCSD, and CCSD(T) methods. All these methods are applied to H2O, three diatomics, and the O4+ molecule. Results demonstrate that the CCSD and OD methods give nearly identical results for H2O and diatomics; however, in symmetry-breaking problems as exemplified by O4+, the OD method provides better results for vibrational frequencies. The OD method has further advantages over CCSD: its analytic gradients are easier to compute since there is no need to solve the coupled-perturbed equations for the orbital response, the computation of one-electron properties are easier because there is no response contribution to the particle density matrices, the variational optimized orbitals can be readily extended to allow inactive orbitals, it avoids spurious second-order poles in its response function, and its transition dipole moments are gauge invariant. The OMP2 has these same advantages over canonical MP2, making it promising for excited state properties via linear response theory. The quadratically convergent orbital-optimization procedure converges quickly for OMP2, and provides molecular properties that are somewhat different than those of MP2 for most of the test cases considered (although they are similar for H2O). Bond lengths are somewhat longer, and vibrational frequencies somewhat smaller, for OMP2 compared to MP2. In the difficult case of O4+, results for several vibrational frequencies are significantly improved in going from MP2 to OMP2. [ABSTRACT FROM AUTHOR]

Although never spectroscopically identified in the laboratory, hydrogenated nitrogen (HN2) is thought to be an important species in combustion chemistry. The classical barrier height (10.6±0.2 kcal mol-1) and exothermicity (3.6±0.2 kcal mol-1) for the HN2→N2+H reaction are predicted by high level ab initio quantum mechanical methods [up to CCSDT(Q)]. Total energies are extrapolated to the complete basis set limit applying the focal point analysis. Zero-point vibrational energies are computed using fundamental (anharmonic) frequencies obtained from a quartic force field. Relativistic and diagonal Born–Oppenheimer corrections are also taken into account. The quantum mechanical barrier with these corrections is predicted to be 6.4±0.2 kcal mol-1 and the reaction exothermicity to be 8.8±0.2 kcal mol-1. The importance of these parameters for the thermal NOx decomposition (De-NOx) process is discussed. The unimolecular rate constant for dissociation of the HN2 molecule and its lifetime are estimated by canonical transition-state theory and Rice–Ramsperger–Kassel–Marcus theory. The lifetime of the HN2 molecule is here estimated to be 2.8×10-10 s at room temperature. Our result is in marginal agreement with the latest experimental kinetic modeling studies (τ=1.5×10-8 s), albeit consistent with the very rough experimental upper limit (τ<0.5 μs). For the dissociation reaction, kinetic isotope effects are investigated. Our analysis demonstrates that the DN2 molecule has a longer lifetime than the HN2 molecule. Thus, DN2 might be more readily identified experimentally. The ionization potential of the HN2 molecule is determined by analogous high level ab initio methods and focal point analysis. The adiabatic IP of HN2 is predicted to be 8.19±0.05 eV, in only fair agreement with the experimental upper limit of 7.92 eV deduced from sychrothon-radiation-based photoionization mass spectrometry. [ABSTRACT FROM AUTHOR]

Although never spectroscopically identified in the laboratory, hydrogenated nitrogen (HN2) is thought to be an important species in combustion chemistry. The classical barrier height (10.6±0.2 kcal mol-1) and exothermicity (3.6±0.2 kcal mol-1) for the HN2→N2+H reaction are predicted by high level ab initio quantum mechanical methods [up to CCSDT(Q)]. Total energies are extrapolated to the complete basis set limit applying the focal point analysis. Zero-point vibrational energies are computed using fundamental (anharmonic) frequencies obtained from a quartic force field. Relativistic and diagonal Born–Oppenheimer corrections are also taken into account. The quantum mechanical barrier with these corrections is predicted to be 6.4±0.2 kcal mol-1 and the reaction exothermicity to be 8.8±0.2 kcal mol-1. The importance of these parameters for the thermal NOx decomposition (De-NOx) process is discussed. The unimolecular rate constant for dissociation of the HN2 molecule and its lifetime are estimated by canonical transition-state theory and Rice–Ramsperger–Kassel–Marcus theory. The lifetime of the HN2 molecule is here estimated to be 2.8×10-10 s at room temperature. Our result is in marginal agreement with the latest experimental kinetic modeling studies (τ=1.5×10-8 s), albeit consistent with the very rough experimental upper limit (τ<0.5 μs). For the dissociation reaction, kinetic isotope effects are investigated. Our analysis demonstrates that the DN2 molecule has a longer lifetime than the HN2 molecule. Thus, DN2 might be more readily identified experimentally. The ionization potential of the HN2 molecule is determined by analogous high level ab initio methods and focal point analysis. The adiabatic IP of HN2 is predicted to be 8.19±0.05 eV, in only fair agreement with the experimental upper limit of 7.92 eV deduced from sychrothon-radiation-based photoionization mass spectrometry. [ABSTRACT FROM AUTHOR]

Since the discovery of ozone depletion, the doublet electronic states of the ozone radical cation have received much attention in experimental and theoretical investigations, while the low-lying quartet states have not. In the present research, viable pathways to the quartet states from the lowest three triplet states of ozone, 3A2, 3B2, and 3B1, and excitations from the 2A1 and 2B2 states of the ozone radical cation have been studied in detail. The potential energy surfaces, structural optimizations, and vibrational frequencies for several states of ozone and its radical cation have been thoroughly investigated using the complete active space self-consistent field, unrestricted coupled cluster theory from a restricted open-shell Hartree-Fock reference including all single and double excitations (UCCSD), UCCSD method with the effects of connected triple excitations included perturbatively, and unrestricted coupled cluster including all single, double, and triple excitations with the effects of connected quadruple excitations included perturbatively. These methods used Dunning’s correlation-consistent polarized core-valence basis sets, cc-pCVXZ (X=D, T, Q, and 5). The most feasible pathways (symmetry and spin allowed transitions) to the quartet states are 4A1←3A2, 4A2←3A2, 4A1←3B2, 4A2←3B1, 4B2←3B1, 4A2←1A1, 4B2←1A1, and 4A1←1A1 with vertical ionization potentials of 12.46, 12.85, 12.82, 12.46, 12.65, 13.43, 13.93, and 14.90 eV, respectively. [ABSTRACT FROM AUTHOR]

The singlet ground (X 1Σ+) and excited (1Σ-,1Δ) states of HCP and HPC have been systematically investigated using ab initio molecular electronic structure theory. For the ground state, geometries of the two linear stationary points have been optimized and physical properties have been predicted utilizing restricted self-consistent field theory, coupled cluster theory with single and double excitations (CCSD), CCSD with perturbative triple corrections [CCSD(T)], and CCSD with partial iterative triple excitations (CCSDT-3 and CC3). Physical properties computed for the global minimum (X 1Σ+HCP) include harmonic vibrational frequencies with the cc-pV5Z CCSD(T) method of ω1=3344 cm-1, ω2=689 cm-1, and ω3=1298 cm-1. Linear HPC, a stationary point of Hessian index 2, is predicted to lie 75.2 kcal mol-1 above the global minimum HCP. The dissociation energy D0[HCP(X 1Σ+)→H(2S)+CP(X 2Σ+)] of HCP is predicted to be 119.0 kcal mol-1, which is very close to the experimental lower limit of 119.1 kcal mol-1. Eight singlet excited states were examined and their physical properties were determined employing three equation-of-motion coupled cluster methods (EOM-CCSD, EOM-CCSDT-3, and EOM-CC3). Four stationary points were located on the lowest-lying excited state potential energy surface, 1Σ-→1A″, with excitation energies Te of 101.4 kcal mol-1(1A″ HCP), 104.6 kcal mol-1(1Σ- HCP), 122.3 kcal mol-1(1A″ HPC), and 171.6 kcal mol-1(1Σ- HPC) at the cc-pVQZ EOM-CCSDT-3 level of theory. The physical properties of the 1A″ state with a predicted bond angle of 129.5° compare well with the experimentally reported first singlet state (Ã 1A″). The excitation energy predicted for this excitation is T0=99.4 kcal mol-1(34 800 cm-1,4.31 eV), in essentially perfect agreement with the experimental value of T0=99.3 kcal mol-1(34 746 cm-1,4.308 eV). For the second lowest-lying excited singlet surface, 1Δ→1A′, four stationary points were found with Te values of 111.2 kcal mol-1 (21A′ HCP), 112.4 kcal mol-1 (1Δ HPC), 125.6 kcal mol-1(2 1A′ HCP), and 177.8 kcal mol-1(1Δ HPC). The predicted CP bond length and frequencies of the 2 1A′ state with a bond angle of 89.8° (1.707 Å, 666 and 979 cm-1) compare reasonably well with those for the experimentally reported C 1A′ state (1.69 Å, 615 and 969 cm-1). However, the excitation energy and bond angle do not agree well: theoretical values of 108.7 kcal mol-1 and 89.8° versus experimental values of 115.1 kcal mol-1 and 113°. [ABSTRACT FROM AUTHOR]

The existence or nonexistence of GaH5 has been widely discussed [N. M. Mitzel, Angew. Chem. Int. Ed. 42, 3856 (2003)]. Seven possible structures for gallium pentahydride have been systematically investigated using ab initio electronic structure theory. Structures and vibrational frequencies have been determined employing self-consistent field, coupled cluster including all single and double excitations (CCSD), and CCSD with perturbative triples levels of theory, with at least three correlation-consistent polarized-valence-(cc-pVXZ and aug-cc-pVXZ) type basis sets. The X 1A′ state for GaH5 is predicted to be weakly bound complex 1 between gallane and molecular hydrogen, with Cs symmetry. The dissociation energy corresponding to GaH5→GaH3+H2 is predicted to be De=2.05 kcal mol-1. The H–H stretching fundamental is predicted to be v=4060 cm-1, compared to the tentatively assigned experimental feature of Wang and Andrews [J. Phys. Chem. A 107, 11371 (2003)] at 4087 cm-1. A second Cs structure 2 with nearly equal energy is predicted to be a transition state, corresponding to a 90° rotation of the H2 bond. Thus the rotation of the hydrogen molecule is essentially free. However, hydrogen scrambling through the C2v structure 3 seems unlikely, as the activation barrier for scrambling is at least 30 kcal mol-1 higher in energy than that for the dissociation of GaH5 to GaH3 and H2. Two additional structures consisting of GaH3 with a dihydrogen bond perpendicular to gallane (C3v structure 4) and an in-plane dihydrogen bond [Cs(III) structure 5] were also examined. A C3v symmetry second-order saddle point has nearly the same energy as the GaH3+H2 dissociation limit, while the Cs(III) structure 5 is a transition structure to the C3v structure. The C4v structure 6 and the D3h structure 7 are much higher in energy than GaH3+H2 by 88 and 103 kcal mol-1, respectively. [ABSTRACT FROM AUTHOR]