Study of the C(3P) + OH(X2Pi) --> CO(X1Sigma(g)+) + H(2S) reaction: a fully global ab initio potential energy surface of the X2A' state

The C((3)P) + OH(X (2)Pi) --> CO(X (1)Sigma(g)(+)) + H((2)S) reaction has been investigated by ab initio electronic structure calculations of the X(2)A' state based on the multireference (MR) internally contracted single and double configuration interaction (SDCI) method plus Davidson correction (+Q) using Dunning aug-cc-pVQZ basis sets. In particular, the multireference space is taken to be a complete active space (CAS). Improvement over previously proposed potential energy surfaces for HCO/COH is obtained in the sense that present surface describes also the potential part where the CO interatomic distance is large. A large number of geometries (around 2000) have been calculated and analytically fitted using the reproducing kernel Hilbert space (RKHS) method of Ho and Rabitz both for the two-body and three-body terms following the many-body decomposition of the total electronic energies. Results show that the global reaction is highly exothermic ( approximately 6.4 eV) and barrierless (relative to the reactant channel), while five potential barriers are located on this surface. The three minima and five saddle points observed are characterized and found to be in good agreement with previous work. The three minima correspond to the formation of HCO and COH complexes and to the CO + H products, with the COH complex being a metastable minimum relative to the product channel. The five saddle points correspond to potential barriers for both the dissociation/formation of HCO and COH into/from CO + H, to barriers for the isomerization of HCO into COH and to barriers for the inversion of HCO and COH through their respective linear configuration.     

Separability of tight and roaming pathways to molecular decomposition

Recent studies have questioned the separability of the tight and roaming mechanisms to molecular decomposition. We explore this issue for a variety of reactions including MgH(2) → Mg + H(2), NCN → CNN, H(2)CO → H(2) + CO, CH(3)CHO → CH(4) + CO, and HNNOH → N(2) + H(2)O. Our analysis focuses on the role of second-order saddle points in defining global dividing surfaces that encompass both tight and roaming first-order saddle points. The second-order saddle points define an energetic criterion for separability of the two mechanisms. Furthermore, plots of the differential contribution to the reactive flux along paths connecting the first- and second-order saddle points provide a dynamic criterion for separability. The minimum in the differential reactive flux in the neighborhood of the second-order saddle point plays the role of a mechanism divider, with the presence of a strong minimum indicating that the roaming and tight mechanisms are dynamically distinct. We show that the mechanism divider is often, but not always, associated with a second-order saddle point. For the formaldehyde and acetaldehyde reactions, we find that the minimum energy geometry on a conical intersection is associated with the mechanism divider for the tight and roaming processes. For HNNOH, we again find that the roaming and tight processes are dynamically separable but we find no intrinsic feature of the potential energy surface associated with the mechanism divider. Overall, our calculations suggest that roaming and tight mechanisms are generally separable over broad ranges of energy covering most kinetically relevant regimes.     

Structure and H(2)-Loss Energies of OsHX(H(2))(CO)L(2) Complexes (L = P(t-Bu)(2)Me, P(i-Pr)(3); X = Cl, I, H): Attempted Correlation of (1)J(H-D), T(1min), and DeltaG()

1J(H-D), T(1min) and k(1) for H(2) dissociation from OsHX(H(2))(CO)L(2) have been measured for X = Cl, I, H (L = P(t-Bu)(2)Me or P(i-Pr)(3)), as well as for OsCl(2)(H(2))(CO)(P(i-Pr)(3))(2). For comparison, new data (including previously unobserved coupling constants) have been reported for W(HD)(CO)(3)(P(i-Pr)(3))(2). A comprehensive consideration of T(1min) data for over 20 dihydrogen complexes containing only 1-2 phosphines cis to H(2), together with a consideration of the shortest "conceivable" H-H distance for H(2) bound to a d(4) or d(6) metal, is used to argue that the "fast spinning" model is not appropriate for determining r(H-H) in such complexes. Regarding OsHX(H(2))(CO)L(2), the stronger electron-donor (lighter) halide, when cis to H(2), facilitates loss of H(2). The complete absence of pi-donor ability when X = H renders H(2) loss most difficult. However, a pi-donor trans to H(2) also makes H(2) loss unobservable. Within the series of isoelectronic, structurally analogous Os complexes, a longer H-H bond shows a larger DeltaG() for H(2) loss. However, this correlation does not continue to W(H(2))(CO)(3)(P(i-Pr)(3))(2), which has r(H-H) comparable to that of OsH(halide)(H(2))(CO)(P(i-Pr)(3))(2), but a significantly higher DeltaG(). This may originate from lack of a pi-donor ligand to compensate as H(2) leaves W.     

Quantum chemical study of the mechanism of reaction between NH (X 3sigma-) and H2, H2O, and CO2 under combustion conditions

Reactions of ground-state NH (3sigma-) radicals with H2, H2O, and CO2 have been investigated quantum chemically, whereby the stationary points of the appropriate reaction potential energy surfaces, that is, reactants, products, intermediates, and transition states, have been identified at the G3//B3LYP level of theory. Reaction between NH and H2 takes place via a simple abstraction transition state, and the rate coefficient for this reaction as derived from the quantum chemical calculations, k(NH + H2) = (1.1 x 10(14)) exp(-20.9 kcal mol(-1)/RT) cm3 mol(-1) s(-1) between 1000 and 2000 K, is found to be in good agreement with experiment. For reaction between triplet NH and H2O, no stable intermediates were located on the triplet reaction surface although several stable species were found on the singlet surface. No intersystem crossing seam between triplet NH + H2O and singlet HNO + H2 (the products of lowest energy) was found; hence there is no evidence to support the existence of a low-energy pathway to these products. A rate coefficient of k(NH + H2O) = (6.1 x 10(13)) exp(-32.8 kcal mol(-1)/RT) cm3 mol(-1) s(-1) between 1000 and 2000 K for the reaction NH (3sigma-) + H2O --> NH2 (2B) + OH (2pi) was derived from the quantum chemical results. The reverse rate coefficient, calculated via the equilibrium constant, is in agreement with values used in modeling the thermal de-NO(x) process. For the reaction between triplet NH and CO2, several stable intermediates on both triplet and singlet reaction surfaces were located. Although a pathway from triplet NH + CO2 to singlet HNO + CO involving intersystem crossing in an HN-CO2 adduct was discovered, no pathway of sufficiently low activation energy was discovered to compare with that found in an earlier experiment [Rohrig, M.; Wagner, H. G. Proc. Combust. Inst. 1994, 25, 993.].     

Theoretical study of the C6H3 potential energy surface and rate constants and product branching ratios of the C2H(2Sigma+) + C4H2(1Sigma(g)+) and C4H(2Sigma+) + C2H2(1Sigma(g)+) reactions

Ab initio CCSD(T)cc-pVTZ//B3LYP6-311G(**) and CCSD(T)/complete basis set (CBS) calculations of stationary points on the C(6)H(3) potential energy surface have been performed to investigate the reaction mechanism of C(2)H with diacetylene and C(4)H with acetylene. Totally, 25 different C(6)H(3) isomers and 40 transition states are located and all possible bimolecular decomposition products are also characterized. 1,2,3- and 1,2,4-tridehydrobenzene and H(2)CCCCCCH isomers are found to be the most stable thermodynamically residing 77.2, 75.1, and 75.7 kcal/mol lower in energy than C(2)H + C(4)H(2), respectively, at the CCSD(T)/CBS level of theory. The results show that the most favorable C(2)H + C(4)H(2) entrance channel is C(2)H addition to a terminal carbon of C(4)H(2) producing HCCCHCCCH, 70.2 kcal/mol below the reactants. This adduct loses a hydrogen atom from the nonterminal position to give the HCCCCCCH (triacetylene) product exothermic by 29.7 kcal/mol via an exit barrier of 5.3 kcal/mol. Based on Rice-Ramsperger-Kassel-Marcus calculations under single-collision conditions, triacetylene+H are concluded to be the only reaction products, with more than 98% of them formed directly from HCCCHCCCH. The C(2)H + C(4)H(2) reaction rate constants calculated by employing canonical variational transition state theory are found to be similar to those for the related C(2)H + C(2)H(2) reaction in the order of magnitude of 10(-10) cm(3) molecule(-1) s(-1) for T = 298-63 K, and to show a negative temperature dependence at low T. A general mechanism for the growth of polyyne chains involving C(2)H + H(C[triple bond]C)(n)H --> H(C[triple bond]C)(n+1)H + H reactions has been suggested based on a comparison of the reactions of ethynyl radical with acetylene and diacetylene. The C(4)H + C(2)H(2) reaction is also predicted to readily produce triacetylene + H via barrierless C(4)H addition to acetylene, followed by H elimination.     

Experimental and theoretical studies of the reactions Y (a2D) + H2CO and Y (a2D) + CH3CHO

The reactions of ground state Y (a(2)D) with H(2)CO and CH(3)CHO were studied at a range of collision energies in crossed molecular beams. For reaction with H(2)CO, three product channels were observed: formation of YH(2) + CO, YCO + H(2), and YHCO + H. Reaction with CH(3)CHO led to three analogous product channels involving formation of HYCH(3) + CO, YCH(2)CO + H(2), and YCH(3)CO + H. The calculated CCSD(T) energetics and DFT geometries for key intermediates in both reactions, together with RRKM theory, are used to calculate a priori the branching ratios between various product channels. These calculated values are compared to those obtained experimentally.     

Complexation of the triply-bonded dirhenium(II) complex Re(2)Cl(4)(mu-dppm)(2) (dppm = Ph(2)PCH(2)PPh(2)) by up to three acetylene molecules

The triply bonded dirhenium(II) synthons Re(2)X(4)(mu-dppm)(2) (X = Cl, Br; dppm = Ph(2)PCH(2)PPh(2)) react with acetylene at room temperature in CH(2)Cl(2) and acetone to afford the bis(acetylene) complexes Re(2)X(4)(mu-dppm)(2)(mu:eta(2),eta(2)-HCCH)(eta(2)-HCCH) (X = Cl (3), Br(4)). Compound 3 has been derivatized by reaction with RNC ligands in the presence of TlPF(6) to give unsymmetrical complexes of the type [Re(2)Cl(3)(mu-dppm)(2)(mu:eta(2),eta(2)-HCCH)(eta(2)-HCCH)(CNR)]PF(6) (R = Xyl (5), Mes (6), t-Bu (7)), in which the RCN ligand has displaced the chloride ligand cis to the eta(2)-HCCH ligand. The reaction of 3 with an additional 1 equiv of acetylene in the presence of TlPF(6) gives the symmetrical all-cis isomer of [Re(2)Cl(3)(mu-dppm)(2)(mu:eta(2),eta(2)-HCCH)(eta(2)-HCCH)(2)]PF(6) (8). The two terminal eta(2)-HCCH ligands in 8 are very labile and can be displaced by CO and XylNC to give the complexes [Re(2)Cl(3)(mu-dppm)(2)(mu:eta(2),eta(2)-HCCH)(L)(2)]Y (L = CO when Y = PF(6) (9); L = CO when Y = (PF(6))(0.5)/(H(2)PO(4))(0.5) (10); L = XylNC when Y = PF(6) (11)). These substitution reactions proceed with retention of the all-cis stereochemistry. Single-crystal X-ray structure determinations have been carried out on complexes 3, 5, 8, 10, and 11. In no instance have we found that the acetylene ligands undergo reductive coupling reactions.     

Quantum chemical prediction of reaction pathways and rate constants for dissociative adsorption of CO(x) and NO(x) on the graphite (0001) surface

We present predictions of reaction rate constants for dissociative adsorption reactions of CO(x) (x = 1, 2) and NO(x) (x = 1, 2) molecules on the basal graphite (0001) surface based on potential energy surfaces (PES) obtained by the integrated ONIOM(B3LYP:DFTB-D) quantum chemical hybrid approach with dispersion-augmented density functional tight binding (DFTB-D) as low level method. Following an a priori methodology developed in a previous investigation of water dissociative adsorption reactions on graphite, we used a C(94)H(24) dicircumcoronene graphene slab as model system for the graphite surface in finite-size molecular structure investigations, and single adsorbate molecules reacting with the pristine graphene sheet. By employing the ONIOM PES information in RRKM theory we predict reaction rate constants in the temperature range between 1,000 and 5,000 K. We find that among CO(x) and NO(x) adsorbate species, the dissociative adsorption reactions of CO(2) and both radical species NO and NO(2) are likely candidates as a cause for high temperature oxidation and erosion of graphite (0001) surfaces, whereas reaction with CO is not likely to lead to long-lived surface defects. High temperature quantum chemical molecular dynamics simulations (QM/MD) at T = 5,000 K using on-the-fly DFTB-D energies and gradients confirm the results of our PES study.     

碳材料电化学腐蚀的质谱-电化学循环伏安(MSCV)法研究 II:催化剂的影响

用MSCV法研究了Pt催化剂作用下, 碳材料在PH0.3～11.8范围内的阳极腐蚀.Pt催化剂明显加速了活性炭及乙炔黑的腐蚀速度, 使这两种电极上CO2开始析出的电位ΦO2/H2O。因此两种碳电极无论作为阳极还是氧阴极工作时均会被腐蚀。Pt对乙炔黑的催化活性明显强于活性炭, 对于活性炭, Pt并不改变CO2的ΦS-PH曲线的形状(ηS-CO2如析出的超电势)(即无论有无Pt它们均为折点在PH～7折线)。然而乙炔黑的情况更为复杂, Pt使CO2的ΦS-PH及ηS-PH形状从折线变为直线; 在CO2的IM-Φ曲线上出现两个波。它们似与乙炔黑上两种表面氧化物相对应。     

Dihydrogen complexes of electrophilic metal centers: observation of Cr(CO)5(H2), W(CO)5(H2) and [Re(CO)5(H2)]+

Photolysis of dichloromethane solutions of M(CO)6 (M = Cr, W) at low temperature in the presence of hydrogen gas affords W(CO)5(H2) (1) and Cr(CO)5(H2) (2). Complexes 1 and 2 are characterized as dihydrogen complexes based on short T1 values for the hydride resonances and a large HD coupling of 35.3 Hz (W) and 35.8 Hz (Cr) in the HD derivatives. A cationic analogue, [Re(CO)5(H2)]+ (3), was prepared by reaction of Re(CO)5Cl with [Et3Si][B(C6F5)4] in fluorobenzene under hydrogen. Complex 3-d1 exhibits JHD = 33.9 Hz. Complex 3 is strongly acidic, with complete deprotonation by diethyl ether; complexes 1 and 2 are moderately acidic. Deprotonation of 1 is complete in the presence of one equivalent of triethylamine.     

Adsorption of CO and O2 on W2C(0001)

CO, O(2), and H(2) adsorption on a clean W(2)C(0001)√13×√13 R ± 13.9° reconstructed surface at room temperature (RT) were investigated using high-resolution electron energy loss spectroscopy (HREELS). The W(2)C(0001) adsorbs CO molecularly and adsorbs O(2) dissociatively, but does not adsorb H(2) at RT. In the CO adsorption system, two C-O stretching (antisymmetric CCO stretching) modes were found at 242.3 meV (1954 cm(-1)) and at 253.0 meV (2041 cm(-1)). The low-frequency site is occupied at first with subsequent conversion to the high-frequency site with increasing coverage. Additionally, a small peak was apparent at 104.5 meV (843 cm(-1)), and a middle peak at 50-51 meV (400-410 cm(-1)), which are assignable to a symmetric stretching mode and a hindered translational mode, respectively, of a CCO (ketenylidene) species. These observations are consistent with the CO adsorption model on top of the surface carbon. For oxygen adsorption, two adsorption states were found at 65.2-68.1 meV (526-549 cm(-1)) and 73.6 meV (594 cm(-1)): typical frequencies to oxygen adsorption on metal surfaces. Results suggest that atomic oxygen adsorption occurred on a threefold hollow site of the second W layer.     

Coupling of carbon monoxide molecules over oxygen-defected UO2(111) single crystal and thin film surfaces

While coupling reactions of carbon-containing compounds are numerous in organometallic chemistry, they are very rare on well-defined solid surfaces. In this work we show that the reductive coupling of two molecules of carbon monoxide to C2 compounds (acetylene and ethylene) could be achieved on oxygen-defected UO2(111) single crystal and thin film surfaces. This result allows in situ electron spectroscopic investigation of a typical organometallic reaction such as carbon coupling and extends it to heterogeneous catalysis and solids. By using high-resolution photoelectron spectroscopy (HRXPS) it was possible to track the changes in surface states of the U and O atoms as well as identify the intermediate of the reaction. Upon CO adsorption U cations in low oxidation states are oxidized to U4+ ions; this was accompanied by an increase of the O-to-U surface ratios. The HRXPS C 1s lines show the presence of adsorbed species assigned to diolate species (-OCH=CHO-) that are most likely the reaction intermediate in the coupling of two CO molecules to acetylene and ethylene.     

Quasiclassical trajectory calculations of the reaction C+C2H2-->l-C3H, c-C3H+H, C3+H2 using full-dimensional triplet and singlet potential energy surfaces

Full-dimensional, density functional theory (B3LYP/6-311g(d,p))-based potential energy surfaces (PESs) are reported and used in quasi-classical calculations of the reaction of C with C(2)H(2). For the triplet case, the PES spans the region of the reactants, the complex region (with numerous minima and saddle points) and the products, linear(l)-C(3)H+H, cyclic(c)-C(3)H+H and c-(3)C(3)+H(2). For the singlet case, the PES describes the complex region and products l-C(3)H+H, c-C(3)H+H and l-(1)C(3)+H(2). The PESs are invariant under permutation of like nuclei and are fit to tens of thousands of electronic energies. Energies and harmonic frequencies of the PESs agree well the DFT ones for all stationary points and for the reactant and the products. Dynamics calculations on the triplet PES find both l-C(3)H and c-C(3)H products, with l-C(3)H being dominant at the energies considered. Limited unimolecular reaction dynamics on the singlet PES find both products in comparable amounts as well as the C(3)+H(2) product.     

Computational study on the reactions of H2O2 on TiO2 anatase (101) and rutile (110) surfaces

This study investigates the adsorption and reactions of H(2)O(2) on TiO(2) anatase (101) and rutile (110) surfaces by first-principles calculations based on the density functional theory in conjunction with the projected augmented wave approach, using PW91, PBE, and revPBE functionals. Adsorption mechanisms of H(2)O(2) and its fragments on both surfaces are analyzed. It is found that H(2)O(2) , H(2)O, and HO preferentially adsorb at the Ti(5c) site, meanwhile HOO, O, and H preferentially adsorb at the (O(2c))(Ti(5c)), (Ti(5c))(2), and O(2c) sites, respectively. Potential energy profiles of the adsorption processes on both surfaces have been constructed using the nudged elastic band method. The two restructured surfaces, the 1/3 ML oxygen covered TiO(2) and the hydroxylated TiO(2), are produced with the H(2)O(2) dehydration and deoxidation, respectively. The formation of main products, H(2)O(g) and the 1/3 ML oxygen covered TiO(2) surface, is exothermic by 2.8 and 5.0 kcal/mol, requiring energy barriers of 0.8 and 1.1 kcal/mol on the rutile (110) and anatase (101) surface, respectively. The rate constants for the H(2)O(2) dehydration processes have been predicted to be 6.65 × 10(-27) T(4.38) exp(-0.14 kcal mol(-1)/RT) and 3.18 × 10(-23) T(5.60) exp(-2.92 kcal mol(-1)/RT) respectively, in units of cm(3) molecule(-1) s(-1).     

Unsaturated Ru(0) species with a constrained bis-phosphine ligand: [Ru(CO)2(tBu2PCH2CH2PtBu2)]2. Comparison to [Ru(CO)2(PtBu2Me)2

The synthesis of Ru(C2H4)(CO)2(dtbpe) (dtbpe = tBu2PC2H4PtBu2), then green [Ru(CO)2(dtbpe)]n is described. In solution, n = 1, while in the solid state, n = 2; the dimer has two carbonyl bridges. DFTPW91, MP2, and CCSD(T) calculations show that the potential energy surface for bending one carbonyl out of the RuP2C(O) plane is essentially flat. Ru(CO)2(dtbpe) reacts rapidly in benzene solution to oxidatively add the H-E bond of H2, HCl, HCCR (R = H, Ph), [HOEt2]BF4, and HSiEt3. The H-C bond of C6HF5 oxidatively adds at 80 degrees C. CO adds, as does the C=C bond of H2C=CHX (X = H, F, Me). The following do not add: N2, THF, acetone, H3COH, and H2O.     

Quantum chemical and theoretical kinetics study of the O(3P) + C2H2 reaction: a multistate process

The potential energy surfaces of the two lowest-lying triplet electronic surfaces 3A' and 3A' for the O(3P) + C2H2 reaction were theoretically reinvestigated, using various quantum chemical methods including CCSD(T), QCISD, CBS-QCI/APNO, CBS-QB3, G2M(CC,MP2), DFT-B3LYP and CASSCF. An efficient reaction pathway on the electronically excited 3A' surface resulting in H(2S) + HCCO(A2A') was newly identified and is predicted to play an important role at higher temperatures. The primary product distribution for the multistate multiwell reaction was then determined by RRKM statistical rate theory and weak-collision master equation analysis using the exact stochastic simulation method. Allowing for nonstatistical behavior of the internal rotation mode of the initial 3A' adducts, our computed primary-product distributions agree well with the available experimental results, i.e., ca. 80% H(2S) + HCCO(X2A' + A2A') and 20% CH2(X3B1) + CO(X1sigma+) independent of temperature and pressure over the wide 300-2000 K and 0-10 atm ranges. The thermal rate coefficient k(O + C2H2) at 200-2000 K was computed using multistate transition state theory: k(T) = 6.14 x 10(-15)T (1.28) exp(-1244 K/T) cm3 molecule(-1) s(-1); this expression, obtained after reducing the CBS-QCI/APNO ab initio entrance barriers by 0.5 kcal/mol, quasi-perfectly matches the experimental k(T) data over the entire 200-2000 K range, spanning 3 orders of magnitude.     

Ab initio potential energy and dipole moment surfaces of (H2O)2

A full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) are reported for the water dimer, (H2O)2. The CCSD(T)-PES is a very precise fit to 19,805 ab initio energies obtained with the coupled-cluster (CCSD(T)) method, using an aug-cc-pVTZ basis. The standard counterpoise correction was applied to approximately eliminate basis set superposition errors. The fit is based on an approach that incorporates the permutational symmetry of identical atoms [Huang, X.; Braams, B.; Bowman, J. M. J. Chem.Phys. 2005, 122, 044308]. The DMS is a fit to the dipole moment obtained with M?ller-Plesset (MP2) theory, using an aug-cc-pVTZ basis. The PES has an RMS fitting error of 31 cm(-1) for energies below 20,000 cm(-1), relative to the global minimum. This surface can describe various internal floppy motions, including various monomer inversions, and isomerization pathways. Ten characteristic stationary points have been located on the surface, four of which are transition structures and the rest are higher order saddle points. Their geometrical and vibrational properties are presented and compared with available previous theoretical work. The CCSD(T)-PES and MP2-DMS dissociate correctly (and symmetrically) to two H2O monomers, with D(e) = 1665.7 cm(-1) (19.93 kJ/mol). Accurate quantum calculations of the zero-point energy of the dimer (using diffusion Monte Carlo) and the monomers (using a vibrational configuration interaction approach) are reported, and these together with D(e) give a value of D0 of 1042 cm(-1) (12.44 kJ/mol). A best estimated value is 1130 cm(-1) (13.5 kJ/mol).     

Structure and binding energies of unsaturated hydrocarbons on Si(001) and Ge(001)

The adsorption of acetylene, ethylene, and benzene on the Si(001) and Ge(001) surfaces is investigated by first-principles density-functional calculations within the generalized-gradient approximation. We find that the adsorption energies of the three hydrocarbons containing a triple bond, a double bond, and a pi-conjugated aromatic ring decrease as the sequence of C2H2>C2H4>C6H6. We also find that the bondings of acetylene, ethylene, and benzene to Ge(001) are much weaker than those to Si(001). As a result, benzene is weakly bound to Ge(001) while it is chemisorbed on Si(001), consistent with temperature-programmed desorption data.     

Explicit correlation and intermolecular interactions: investigating carbon dioxide complexes with the CCSD(T)-F12 method

We have optimized the lowest energy structures and calculated interaction energies for the CO(2)-Ar, CO(2)-N(2), CO(2)-CO, CO(2)-H(2)O, and CO(2)-NH(3) dimers with the recently developed explicitly correlated coupled cluster singles doubles and perturbative triples [CCSD(T)]-F12 methods and the associated VXZ-F12 (where X = D,T,Q) basis sets. For a given cardinal number, we find that results obtained with the CCSD(T)-F12 methods are much closer to the CCSD(T) complete basis set limit than the conventional CCSD(T) results. The relatively modest increase in the computational cost between explicit and conventional CCSD(T) is more than compensated for by the impressive accuracy of the CCSD(T)-F12 method. We recommend use of the CCSD(T)-F12 methods in combination with the VXZ-F12 basis sets for the accurate determination of equilibrium geometries and interaction energies of weakly bound electron donor acceptor complexes.     

Direct dynamics simulations of O(3P) + HCl at hyperthermal collision energies

The dynamics of the O(3P) + HCl reaction at hyperthermal collision energies were investigated using the quasiclassical trajectory method. Stationary points on the OClH 3A" and 3A' potential energy surfaces (PESs) were also examined. The lowest transition state leading to OCl + H on the 3A" surface is 2.26 eV above the reagents at the CCSD(T)/cc-pVTZ level of theory. This saddle point is bent and product-like. Direct dynamics calculations at the MP2/cc-pVTZ level of theory were used to investigate the excitation functions for OH + Cl, OCl + H, and O + H + Cl formation. OCl is formed mainly from small-impact-parameter collisions, and the OCl + H excitation function peaks around 5 eV, where it is similar in magnitude to the OH + Cl excitation function. The shape of the OCl + H excitation function is discussed, and features are identified that should be general to hyperthermal collision dynamics.