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.     

自由基HO2与C2H2反应势能面的理论研究

应用量子化学从头计算和密度泛函理论(DFT)对HO2+C2H2反应体系的反应机理进行了研究.在B3LYP/6-311G**和CCSD(T)/6-311G**水平上计算了HO2+ C2H2反应的二重态反应势能面.计算结果表明,主要反应方式为自由基HO2的H原子和C2H2分子中的C原子结合,经过一系列异构化,最后分解得到主要产物P1 (CH2O+ HCO).此反应是放热反应,化学反应热为-321.99 kJ·mol-1.次要产物为P2 (CO2 +CH3),也是放热反应.     

Interpolated potential energy surfaces and dynamics for atom exchange between H and H3(+), and D and H3(+)

Two ab initio interpolated potential energy surfaces have been constructed to study the dynamics of atomic hydrogen/deuterium exchange in collisions of H(3)(+) with H (D). One of the surfaces is based on energy calculations using quadratic configuration interaction with single and double excitations. The second includes a perturbative treatment of the triple excitations and an additive correction for basis set deficiency. Results from classical dynamics simulation of the exchange reaction on these surfaces are presented and discussed.     

Quasiclassical trajectory study of formaldehyde unimolecular dissociation: H(2)CO-->H2 + CO, H + HCO

We report quasiclassical trajectory calculations of the dynamics of the two reaction channels of formaldehyde dissociation on a global ab initio potential energy surface: the molecular channel H(2)CO-->H(2) + CO and the radical H(2)CO-->H + HCO. For the molecular channel, it is confirmed that above the threshold of the radical channel a second, intramolecular hydrogen abstraction pathway is opened to produce CO with low rotation and vibrationally hot H(2). The low-j(CO) and high-nu(H(2) ) products from the second pathway increase with the total energy. The competition between the molecular and radical pathways is also studied. It shows that the branching ratio of the molecular products decreases with increasing energy, while the branching ratio of the radical products increases. The results agree well with very recent velocity-map imaging experiments of Suits and co-workers and solves a mystery first posed by Moore and co-workers. For the radical channel, we present the translational energy distributions and HCO rotation distributions at various energies. There is mixed agreement with the experiments of Wittig and co-workers, and this provides an indirect confirmation of their speculation that the triplet surface plays a role in the formation of the radical products.     

The quantum mechanical calculation of rotational spectra. A comparison of methods for C2H2, HCN, HNC, HCO+, N2H+, CO and N2. Predictions for HCNH+, COH+, HBO, HBNH, and HBF+

Rotational frequencies determined with ab initio molecular orbital theory can play an important role in guiding spectroscopic searches for new molecules and in corroborating the assignment of unidentified lines, from the laboratory and from space. In a systematic study of 22 levels of molecular orbital theory, CISD/6-311G** gave rotational frequencies to an accuracy of +/- 0.4 GHz when an empirical correction is applied to the results for C2H2,HCN, HNC, HCO+, N2H+, CO, and N2. Larger errors can be expected when there are large vibrational effects on the rotational constants, as exemplified by COH+. Predicted J = 0--> 1 rotational frequencies using these methods are 73.9 +/- 0.4 GHz for HCNH+, 78.6 +/- 0.4 GHz for HBO, 65.8 +/- 0.4 GHz for HBNH, and 72.1 +/- 0.4 GHz for HBF+.     

Solvated CH5+ in liquid superacid

The transition states for methane activation in liquid superacid have been studied by experimentally determined secondary kinetic deuterium isotope effects (SKIEs) and computational chemistry. For the first time, the SKIEs on hydrogen/deuterium exchange of methane have been measured by using the methane isotopologues in homogeneous liquid superacid (2HF/SbF5). To achieve high accuracy of the SKIEs, the rate constants for pairs of methane isotopologues were simultaneously measured in the same superacid solution by using NMR spectroscopy. Density functional theory (DFT) and high-level ab initio methods have been employed to model possible intermediates and transition states, assuming that the superacids involved in the exchange reactions are H2F+ ions solvated by HF. Only the unsolvated superacid H2F+ is found to be strong enough to protonate methane, yielding the methonium ion solvated by HF as a potential energy minimum. In contrast, the (HF)x-solvated H2F+ superacids (x = 1-4) do not appear to be strong enough to yield stable solvated methonium ions. However, such ions show up as parts of the transition states of the exchange in which the methonium ions are solvated by (HF)x. The calculated DFT activation barrier is in good agreement with that experimentally observed.     

Isomeric differentiation of green tea catechins using gas-phase hydrogen/deuterium exchange reactions

Hydrogen/deuterium exchange reactions in a quadrupole ion trap mass spectrometer are used to differentiate galloylated catechin stereoisomers (catechin gallate and epicatechin gallate; gallocatechin gallate and epigallocatechin gallate) and the nongalloylated analogs (catechin and epicatechin, gallocatechin and epigallocatechin). Significant differences in the hydrogen/deuterium exchange behavior of the four pairs of deprotonated catechin stereoisomers are observed upon reaction with D(2)O. Interestingly, the nongalloylated catechins undergo H/D exchange to a much greater extent than the galloylated species, incorporating deuterium at both aromatic/allylic and active phenolic sites. Nongalloylated catechin isomers are virtually indistinguishable by their H/D exchange kinetics over a wide range of reaction times (0.05 to 10 s). Our experimental results are explained using high-level ab initio calculations to elucidate the subtle structural variations in the catechin stereoisomers that lead to their differing H/D exchange kinetics.     

Classical trajectory study of the reaction between H and HCO

The thermal reaction between H and HCO was studied by classical trajectory calculations on an ab initio potential. The formation of H2+CO, the exchange of hydrogen atoms, and nonreactive encounters, proceeding either via direct or via complex-forming pathways, were separated. The reaction H+HCO-->H2+CO, with direct and complex-forming components, was found to have a low-pressure rate coefficient of 2.0(+/-0.15)x10(-10) cm3 molecule(-1) s(-1), being nearly independent of temperature between 200 and 1000 K. This value is in agreement with the recent experimental value of 1.83(+/-0.4)x10(-10) cm3 molecule(-1) s(-1). Thermal lifetime distributions of H2CO* complexes formed in the reaction are only weakly temperature dependent due to a compensation of energy and angular momentum effects.     

Transition state dynamics of the OH + H2O hydrogen exchange reaction studied by dissociative photodetachment of H3O2-

Dynamics in the transition state region of the bimolecular OH + H2O-->H2O + OH hydrogen exchange reaction have been studied by photoelectron-photofragment coincidence spectroscopy of the H3O2- negative ion and its deuterated analog D3O2-. The data reveal vibrationally resolved product translational energy distributions. The total translational energy distribution shows a vibrational progression indicating excitation of the antisymmetric stretch of the water product. Electronic structure calculations at the QCISD level of theory support this analysis. Examination of the translational energy release between the neutral products reveals a dependence on the product vibrational state. These data should provide a critical test of ab initio potential energy surfaces and dynamics calculations.     

Calculating geochemical reaction pathways--exploration of the inner-sphere water exchange mechanism in Al(H2O)6(3+)(aq) + nH2O with ab Initio calculations and molecular dynamics

We have simulated exchange of inner-sphere and bulk water molecules for different sizes of Al3+(aq) clusters, Al(H2O)63+ + nH2O for n = 0, 1, 6, or 12, with ab initio and molecular dynamics simulations, in order to understand how robust the ab initio method is for identifying hydrolytic reaction pathways of particular importance to geochemistry. In contrast to many interfacial reactions, this particular elementary reaction is particularly simple and well-constrained by experiment. Nevertheless, we find that a rich array of parallel reaction pathways depend sensitively on the details of the solvation sphere and structure and that larger clusters are not necessarily better. Inner-sphere water exchange in Al3+(aq) may occur through two Langford-Gray dissociative pathways, one in which the incoming and outgoing waters are cis, the other in which they are trans to one another. A large majority of exchanges in the molecular dynamics simulations occurred via the trans mechanism, in contrast to the predictions of the ab initio method. In Al(H2O)63+ + H2O, the cis mechanism has a transition state of 84.3 kJ/mol, which is in good agreement with previous experimental and ab initio results, while the trans mechanism has only a saddle point with two negative frequencies, not a transition state, at 89.7 kJ/mol. In addition to the exchange mechanisms, dissociation pathways could be identified that were considerably lower in energy than experiment and varied considerably between 60 and 100 kJ/mol, depending on the particular geometry and cluster size, with no clear relation between the two. Ab initio calculations using large clusters with full second coordination spheres (n = 12) were unable to find dissociation or exchange transition states because the network of hydrogen bonds in the second coordination sphere was too rigid to accommodate the outgoing inner-sphere water. Our results indicate that caution should surround ab initio simulation of complicated dynamic processes such as hydrolysis, ion exchange, and interfacial reactions that involve several steps. Dynamic methods of simulation need to accompany static methods such as ab initio calculation, and it is best to consider simulated pathways as hypotheses to be tested experimentally rather than definitive properties of the reaction.     

Energetics and structural elucidation of mechanisms for gas phase H/D exchange of protonated peptides

Hydrogen/deuterium exchange reactions involving protonated triglycine and deuterated ammonia (ND(3)) have been examined in the gas phase using a Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. Ab initio and density functional theory (DFT) calculations have been carried out to model the exchanges and to obtain energetics and vibrational frequencies for molecules involved in the proposed exchange mechanisms. Structural optimization and frequency calculations have been performed at the B3LYP level of theory with the 6-311+G(d,p) basis set. Transition states have been calculated at the same level of theory and basis set as above using the QST2 and QST3 methods. Single-point energy calculations have been performed at the MP2/6-311+G(d,p) level. Six labile sites of protonated triglycine were found to undergo H/D exchange. Of these six labile hydrogens, two are amide, three are ammonium, and one is carboxyl. Detailed mechanisms for each of these transfers are proposed. Qualitative onium ion and tautomer mechanisms for the exchanges of ammonium and amide hydrogens, respectively, using semiempirical calculations were suggested in previous studies by Beauchamp et al. As shown by the current ab initio and DFT calculations completed during this study, the mechanisms proposed in that study are notionally correct; however, the tautomer mechanisms are shown here to be the result of the fact that a second stable isomer of protonated triglycine exists in which the amide1 carbonyl oxygen is protonated. The exchange of the carboxyl hydrogen is found to proceed via a transition state resembling an ammonium ion interacting with a carboxylate moiety via two hydrogen bonds. The current work thus provides significant mechanistic and structural detail for a considerably more in-depth understanding of the processes involved in gas phase H/D exchange of peptides.     

Gas phase hydrogen deuterium exchange reactions of a model peptide: FT-ICR and computational analyses of metal induced conformational mutations

We utilized gas phase hydrogen/deuterium (H/D) exchange reactions and ab initio calculations to investigate the complexation between a model peptide (Arg-Gly-AspRGD) with various alkali metal ions. The peptide conformation is drastically altered upon alkali metal ion complexation. The associated conformational changes depend on both the number and type of complexing alkali metal ions. Sodium has a smaller ionic diameter and prefers a multidentate interaction that involves all three amino acids of the peptide. Conversely, potassium and cesium form different types of complexes with the RGD. The [RGD + 2Cs − H]⁺ species exhibit the slowest H/D exchange reactivity (reaction rate constant of 6 × 10⁻¹³ cm³molecule⁻¹s⁻¹ for the fastest exchanging labile hydrogen with ND₃). The reaction rate constant of the protonated RGD is two orders of magnitude faster than that of the [RGD + 2Cs − H]⁺. Addition of the first cesium to the RGD reduces the H/D exchange reaction rate constant (i.e., D₀) by a factor of seven whereas sodium reduces this value by a factor of thirty. Conversely, addition of the second alkali metal ions has the opposite effect; the rate of D₀ disappearance for all [RGD + 2Met − H]⁺ species (MetNa, K, and Cs) decreases with the alkali metal ion size.     

Theoretical study of the dynamics of the H+CH4 and H+C2H6 reactions using a specific-reaction-parameter semiempirical Hamiltonian

We present a theoretical study of the reactions of hydrogen atoms with methane and ethane molecules and isotopomers. High-accuracy electronic-structure calculations have been carried out to characterize representative regions of the potential-energy surface (PES) of various reaction pathways, including H abstraction and H exchange. These ab initio calculations have been subsequently employed to derive an improved set of parameters for the modified symmetrically-orthogonalized intermediate neglect of differential overlap (MSINDO) semiempirical Hamiltonian, which are specific to the H+alkane family of reactions. The specific-reaction-parameter (SRP) Hamiltonian has then been used to perform a quasiclassical-trajectory study of both the H+CH4 and H+C2H6 reactions. The calculated values of dynamics properties of the H+CH4-->H2+CH3 reaction and isotopologues, including alkyl product speed distributions, diatomic product internal-state distributions, and cross sections, are generally in good agreement with experiment and with the results provided by the ZBB3 PES [Z. Xie et al., J. Chem. Phys. 125, 133120 (2006)]. The results of trajectories propagated with the SRP Hamiltonian for the H+C2H6-->H2+C2H5 reaction also agree with experiment. The level of agreement between the results calculated with the SRP Hamiltonian and experiment in both the H+methane and H+ethane reactions indicates that semiempirical Hamiltonians can be improved for not only a specific reaction but also a family of reactions.     

The computation of the potential energy surface for H2 + OH → H2O + H using ab initio and density functional theory methods

The reaction energy profile for H₂ + OH → H + H₂O was computed using HF, MP2, MP4, QCISD, G1, G2, and G2MP2 ab initio methods. In addition, the B3LYP, B3P86, B3PW91, BLYP, BP291, and SVWN density functional theory (DFT) methods were also used. All the ab initio methods, with the exception of the G series, produced much higher activation barriers and heats of reaction than the experimental values. On the other hand, the DFT methods produced negative forward and reverse barriers which were too low, with the exception of the hybrid DFT methods. The G2 ab initio method generated energies which deviated from the experimental values by ∼ 1 kcal/mol and therefore should be considered a very accurate computational method. The hybrid DFT methods produced positive forward reaction barriers with energies that were 2–4 kcal/mol lower than the experimental values. The geometries of the transition state and energies computed by the ab initio and DFT methods were compared. These results suggest that, in the hybrid exchange functional, the portion of the Slater exchange term should be increased. This may be the reason why the computed energies were too low. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 62: 639–644, 1997     

Gas-phase hydrogen/deuterium exchange and conformations of deprotonated flavonoids and gas-phase acidities of flavonoids

Gas-phase hydrogen/deuterium (H/D) exchange was used to probe the conformations, gas-phase acidities, and sites of deprotonation of isomeric flavonoid aglycons and glycosides. The flavonoids in each isomer series were differentiated on the basis of their relative rate constants and total numbers of exchanges. For example, flavonoids that possess neohesperidose-type disaccharides may undergo faster and far more extensive exchange than isomeric rutinoside flavonoids. The structural factors that promote or prevent H/D exchange were identified and correlated with collisionally activated dissociation (CAD) patterns and/or molecular modeling data (both high-level ab initio acidity calculations and conformational analysis with molecular dynamics (MD) simulations), thus providing a framework for the use of H/D exchange reactions in the structural elucidation of new flavonoids.     

Quantum and quasiclassical state-to-state dynamics of the NH + H reaction: competition between abstraction and exchange channels

Quantum and quasiclassical state-to-state dynamics for the NH + H' reaction at high collision energies up to 1.6 eV was studied on an accurate ab initio potential energy surface. Both of the endothermic abstraction (NH + H' → N + HH') and thermoneutral exchange (NH + H' → H + NH') channels were investigated from the same set of wave packets using an efficient coordinate transformation method. It is found that the abstraction represents a minor reaction channel in the energy range studied, primarily due to endothermicity. The cross section for the abstraction reaction increases monotonically with the collision energy, while that for the exchange reaction is relatively energy insensitive. As a result, the thermal rate constant for the abstraction reaction follows the Arrhenius law, where that for the exchange reaction is nearly temperature independent. Finally, it is shown that the quantum mechanical results can be reasonably reproduced by the Gaussian-binning quasiclassical trajectory method and to a lesser extent by a quantum statistical model.     

Quantum scattering calculation for reaction Br + H2 on two potential energy surfaces

Three‐dimensional time‐dependent quantum wave packet calculations have been carried out for Br + H₂ on a new global ab initio and a semi‐empirical extended London–Eyring–Polanyi–Sato potential energy surface. It is shown that on the ab initio surface, the threshold energy is much lower, and the reaction probabilities, cross sections, and rate constants are much larger. The effects of the initial rovibrational excitation have also been studied. Comparison of rate constants with experimental measurement implies that the ab initio surface is more suitable for quantum dynamic calculation. The possible reasons and mechanism for the dynamical difference on the two PES are analyzed and discussed. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Modeling the global potential energy surface of the N + N2 reaction from ab initio data

A new global potential energy surface for the N + N2 exchange reaction has been built from ab initio data. To overcome the difficulty of carrying out ab initio calculations for a large set of geometries the alternative strategy of fitting the minimum energy paths of the surface and their angular dependence using a modified LAGROBO functional form has been adopted. In this way we have been able to reproduce all the main features of the potential using a fairly small set of ab initio values.     

Dissociation dynamics of isotopologs of CH₅ studied by charge exchange of CH₅+ with Cs and quasiclassical trajectory calculations

The product branching fractions and dissociation dynamics for the dissociation of CH?, CH?D?, CHD?, and CD? have been experimentally measured following charge exchange of the parent cations with Cs. These results are compared with quasiclassical trajectory calculations run on an ab initio potential energy surface for CH?. In all cases it is found that dissociation channels involving the ejection of a single atom (H or D) are dominant. The statistically weighted branching fractions show that in the mixed isotopologs, ejection of a hydrogen atom is favored over a deuterium atom. The results are consistent with the fluxional nature of CH?(+) isotopologs--no evidence was found for quantum localization upon deuteration in the comparison of the experimental and theoretical branching fractions, possibly as a result of the vibrational temperature of the cations (estimated to be as high as 1660 K for high frequency modes). This comparison of experiment and theory provides a test of the CH? potential energy surface at high levels of excitation.     

Interpolated potential energy surface for abstraction and exchange reactions of NH 3 + H and deuterated analogues

An ab initio interpolated potential energy surface for the hydrogen abstraction and exchange reactions between ammonia and a hydrogen atom is reported. The interpolation is constructed over a set of data points calculated at the unrestricted coupled cluster approximation, using single and double excitations, and including the triple excitations non-iteratively. New data point selection methods were used to improve the convergence and accuracy of the interpolated surface.