Molecular dynamics investigations of ozone on an ab initio potential energy surface with the utilization of pattern-recognition neural network for accurate determination of product formation

The singlet-triplet transformation and molecular dissociation of ozone (O(3)) gas is investigated by performing quasi-classical molecular dynamics (MD) simulations on an ab initio potential energy surface (PES) with visible and near-infrared excitations. MP4(SDQ) level of theory with the 6-311g(2d,2p) basis set is executed for three different electronic spin states (singlet, triplet, and quintet). In order to simplify the potential energy function, an approximation is adopted by ignoring the spin-orbit coupling and allowing the molecule to switch favorably and instantaneously to the spin state that is more energetically stable (lowest in energy among the three spin states). This assumption has previously been utilized to study the SiO(2) system as reported by Agrawal et al. (J. Chem. Phys. 2006, 124 (13), 134306). The use of such assumption in this study probably makes the upper limits of computed rate coefficients the true rate coefficients. The global PES for ozone is constructed by fitting 5906 ab initio data points using a 60-neuron two-layer feed-forward neural network. The mean-absolute error and root-mean-squared error of this fit are 0.0446 eV (1.03 kcal/mol) and 0.0756 eV (1.74 kcal/mol), respectively, which reveal very good fitting accuracy. The parameter coefficients of the global PES are reported in this paper. In order to identify the spin state with high confidence, we propose the use of a pattern-recognition neural network, which is trained to predict the spin state of a given configuration (with a prediction accuracy being 95.6% on a set of testing data points). To enhance the prediction effectiveness, a buffer series of five points are validated to confirm the spin state during the MD process to gain better confidence. Quasi-classical MD simulations from 1.2 to 2.4 eV of total internal energy (including zero-point energy) result in rate coefficients of singlet-triplet transformation in the range of 0.027 ps(-1) to 1.21 ps(-1). Also, we find very low dissociation probability up to 2.4 eV of internal energy during the investigating period (5 ps), which suggests that dissociation does not occur directly from the singlet ground-state, but it involves the excited triplet-state as an intermediate step and requires more reaction time to occur.     

Unimolecular decomposition of chemically activated pentatetraene (H2CCCCCH2) intermediates: A crossed beams study of dicarbon molecule reactions with allene

The reactions dynamics of the dicarbon molecule C2 in the 1Sigma (g)+ singlet ground state and 3Pi(u) first excited triplet state with allene, H2CCCH2(X1A1), was investigated under single collision conditions using the crossed molecular beam approach at four collision energies between 13.6 and 49.4 kJ mol(-1). The experiments were combined with ab initio electronic structure calculations of the relevant stationary points on the singlet and triplet potential energy surfaces. Our investigations imply that the reactions are barrier-less and indirect on both the singlet and the triplet surfaces and proceed through bound C5H4 intermediates via addition of the dicarbon molecule to the carbon-carbon double bond (singlet surface) and to the terminal as well as central carbon atoms of the allene molecule (triplet surface). The initial collision complexes isomerize to form triplet and singlet pentatetraene intermediates (H2CCCCCH2) that decompose via atomic hydrogen loss to yield the 2,4-pentadiynyl-1 radical, HCCCCCH2(X2B1). These channels result in symmetric center-of-mass angular distributions. On the triplet surface, a second channel involves the existence of a nonsymmetric reaction intermediate (HCCCH2CCH) that fragments through atomic hydrogen emission to the 1,4-pentadiynyl-3 radical [C5H3(X2B1)HCCCHCCH]; this pathway was found to account for the backward scattered center-of-mass angular distributions at higher collision energies. The identification of two resonance-stabilized free C5H3 radicals (i.e., 2,4-pentadiynyl-1 and 1,4-pentadiynyl-3) suggests that these molecules can be important transient species in combustion flames and in the chemical evolution of the interstellar medium.     

The lowest-lying electronic singlet and triplet potential energy surfaces for the HNO-NOH system: energetics, unimolecular rate constants, tunneling and kinetic isotope effects for the isomerization and dissociation reactions

The lowest-lying electronic singlet and triplet potential energy surfaces (PES) for the HNO-NOH system have been investigated employing high level ab initio quantum chemical methods. The reaction energies and barriers have been predicted for two isomerization and four dissociation reactions. Total energies are extrapolated to the complete basis set limit applying focal point analyses. Anharmonic zero-point vibrational energies, diagonal Born-Oppenheimer corrections, relativistic effects, and core correlation corrections are also taken into account. On the singlet PES, the (1)HNO → (1)NOH endothermicity including all corrections is predicted to be 42.23 ± 0.2 kcal mol(-1). For the barrierless decomposition of (1)HNO to H + NO, the dissociation energy is estimated to be 47.48 ± 0.2 kcal mol(-1). For (1)NOH → H + NO, the reaction endothermicity and barrier are 5.25 ± 0.2 and 7.88 ± 0.2 kcal mol(-1). On the triplet PES the reaction energy and barrier including all corrections are predicted to be 7.73 ± 0.2 and 39.31 ± 0.2 kcal mol(-1) for the isomerization reaction (3)HNO → (3)NOH. For the triplet dissociation reaction (to H + NO) the corresponding results are 29.03 ± 0.2 and 32.41 ± 0.2 kcal mol(-1). Analogous results are 21.30 ± 0.2 and 33.67 ± 0.2 kcal mol(-1) for the dissociation reaction of (3)NOH (to H + NO). Unimolecular rate constants for the isomerization and dissociation reactions were obtained utilizing kinetic modeling methods. The tunneling and kinetic isotope effects are also investigated for these reactions. The adiabatic singlet-triplet energy splittings are predicted to be 18.45 ± 0.2 and 16.05 ± 0.2 kcal mol(-1) for HNO and NOH, respectively. Kinetic analyses based on solution of simultaneous first-order ordinary-differential rate equations demonstrate that the singlet NOH molecule will be difficult to prepare at room temperature, while the triplet NOH molecule is viable with respect to isomerization and dissociation reactions up to 400 K. Hence, our theoretical findings clearly explain why (1)NOH has not yet been observed experimentally.     

Ab initio energies and product branching ratios for the O+C3H6 reaction

Intermediate and transition-state energies have been calculated for the O+C3H6 (propene) reaction using the compound ab initio CBS-QB3 and G3 methods in combination with density functional theory. The lowest-lying triplet and singlet potential energy surfaces of the O-C3H6 system were investigated. RRKM statistical theory was used to predict product branching fractions over the 300-3000 K temperature and 0.001-760 Torr pressure ranges. The oxygen atom adds to the C3H6 terminal olefinic carbon in the primary step to generate a nascent triplet biradical, CH3CHCH2O. On the triplet surface, unimolecular dissociation of CH3CHCH2O to yield H+CH3CHCHO is favored over the entire temperature range, although the competing H2CO+CH3CH product channel becomes significant at high temperature. Rearrangement of triplet CH3CHCH2O to CH3CH2CHO (propanal) via a 1,2 H-atom shift has a barrier of 122.3 kJ mol(-1), largely blocking this reaction channel and any subsequent dissociation products. Intersystem crossing of triplet CH3CHCH2O to the singlet surface, however, leads to facile rearrangement to singlet CH3CH2CHO, which dissociates via numerous product channels. Pressure was found to have little influence over the branching ratios under most conditions, suggesting that the vibrational self-relaxation rates for p相似文献   

Vibrational cooling of spin-stretched dimer states by He buffer gas: quantum calculations for Li2(a 3Sigma(u)+) at ultralow energies

The interaction between the triplet state of the lithium dimer, (7)Li(2), with (4)He is obtained from accurate ab initio calculations where the vibrational dependence of the potential is newly computed. Vibrational quenching dynamics within a coupled-channel quantum treatment is carried out at ultralow energies, and large differences in efficiency as a function of the initial vibrational state of the targets are found as one compares the triplet results with those of the singlet state of the same target.     

Photodissociation dynamics of nitrobenzene and o-nitrotoluene

Photodissociation of nitrobenzene at 193, 248, and 266 nm and o-nitrotoluene at 193 and 248 nm was investigated separately using multimass ion imaging techniques. Fragments corresponding to NO and NO(2) elimination from both nitrobenzene and o-nitrotoluene were observed. The translational energy distributions for the NO elimination channel show bimodal distributions, indicating two dissociation mechanisms involved in the dissociation process. The branching ratios between NO and NO(2) elimination channels were determined to be NONO(2)=0.32+/-0.12 (193 nm), 0.26+/-0.12 (248 nm), and 0.4+/-0.12(266 nm) for nitrobenzene and 0.42+/-0.12(193 nm) and 0.3+/-0.12 (248 nm) for o-nitrotoluene. Additional dissociation channels, O atom elimination from nitrobenzene, and OH elimination from o-nitrotoluene, were observed. New dissociation mechanisms were proposed, and the results are compared with potential energy surfaces obtained from ab initio calculations. Observed absorption bands of photodissociation are assigned by the assistance of the ab initio calculations for the relative energies of the triplet excited states and the vertical excitation energies of the singlet and triplet excited states of nitrobenzene and o-nitrotoluene. Finally, the dissociation rates and lifetimes of photodissociation of nitrobenzene and o-nitrotoluene were predicted and compared to experimental results.     

Density functional and Ab initio study of Cr(CO)n (n = 1-6) complexes

Cr(CO)n (n = 1-6) systems were studied for all possible spin states using density functional and high-level ab initio methods to provide a more complete theoretical understanding of the structure of species that may form during ligand dissociation of Cr(CO)6. We carried out geometry optimizations for each system and obtained vibrational frequencies, sequential bond dissociation energies (BDE), and total CO binding energies. We also compared the performance of various DFT functionals. Generally, the ground states of Cr(CO)6, Cr(CO)5, and Cr(CO)4, whose spin multiplicity is a singlet, are in good agreement with both previous theoretical results and currently available experimental data. Calculations on Cr(CO)3, Cr(CO)2, and CrCO provide new findings that the ground state of Cr(CO)3 might be a quintet with C2v symmetry instead of a singlet with C3v symmetry, and the ground state of Cr(CO)2 is not a linear quintet, as suggested by previous DFT calculations, but rather a linear septet. We also found that nonet states of Cr(CO)2 and CrCO display partial C-O bond breakage.     

Ab initio study on the reaction of Y+ + NH3 → Y+ NH + H2

The reaction of Y⁺ + NH₃ → Y⁺ NH + H₂ was theoretically investigated by ab initio MO methods. Two possible pathways (1–1 H₂ loss and 1–2 H₂ loss) on the singlet potential energy surface and reaction mechanism were examined and discussed. The singlet and triplet PESs of this reaction system were compared to confirm the correctness of spin conservation concepts. © 1996 John Wiley & Sons, Inc.     

Theoretical investigation of excited states of C(3)

In this work, we present ab initio calculations for the potential energy surfaces of C(3) in different electronic configurations, including the singlet ground state [X (1)Sigma(g) (+),((1)A(1))], the triplet ground state [a (3)Pi(u),((3)B(1), (3)A(1))], and some higher excited states. The geometries studied include triangular shapes with two identical bond lengths, but different bond angles between them. For the singlet and triplet ground states in the linear geometry, the total energies resulting from the mixed density functional--Hartree-Fock and quadratic configuration interaction methods reproduce the experimental values, i.e., the triplet occurs 2.1 eV above the singlet. In the geometry of an equilateral triangle, we find a low-lying triplet state with an energy of only 0.8 eV above the energy of the singlet in the linear configuration, so that the triangular geometry yields the lowest excited state of C(3). For the higher excited states up to about 8 eV above the ground state, we apply time-dependent density functional theory. Even though the systematic error produced by this approach is of the order of 0.4 eV, the results give different prospective to insight into the potential energy landscape for higher excitation energies.     

Molecular dynamics and metadynamics simulations of [2 + 2] photocycloaddition

Triplet state mechanism of [2 + 2] photocycloaddition forming a cyclobutane ring from two ethylenes is investigated in the context of photocatalysis. High‐level ab initio calculations are combined with ab initio adiabatic molecular dynamics and ab initio metadynamics for rare events modeling. In a photocatalytic scheme, a reactant reaches the triplet state either via intersystem crossing (ISC) or triplet sensitization. The model system adopts a biradical structure, which represents energy intersection with the ground state. The system either completes cyclization or undergoes fragmentation into two olefinic units. The potential and free energy surfaces of the cyclobutane/ethylenes system are mapped with multireference approaches describing possible reaction pathways. To obtain a full picture of a double bond photoreactivity, ab initio adiabatic dynamical calculations were used to estimate reaction yields and to model the effects of excess energy. The potential use of density functional theory based approaches for [2 + 2] photocycloaddition was investigated for future simulations and design of realistic photocatalytic systems. Dynamical aspects of [2 + 2] photocycloaddition via a triplet state manifold are investigated by combining ab initio multireference methods and ab initio molecular dynamics and metadynamics. The reaction pathways are studied for a model system of two ethylenes forming a cyclobutane ring to provide a basis for further studies on design of photocatalytic systems.     

Parametrization of analytic interatomic potential functions using neural networks

A generalized method that permits the parameters of an arbitrary empirical potential to be efficiently and accurately fitted to a database is presented. The method permits the values of a subset of the potential parameters to be considered as general functions of the internal coordinates that define the instantaneous configuration of the system. The parameters in this subset are computed by a generalized neural network (NN) with one or more hidden layers and an input vector with at least 3n-6 elements, where n is the number of atoms in the system. The Levenberg-Marquardt algorithm is employed to efficiently affect the optimization of the weights and biases of the NN as well as all other potential parameters being treated as constants rather than as functions of the input coordinates. In order to effect this minimization, the usual Jacobian employed in NN operations is modified to include the Jacobian of the computed errors with respect to the parameters of the potential function. The total Jacobian employed in each epoch of minimization is the concatenation of two Jacobians, one containing derivatives of the errors with respect to the weights and biases of the network, and the other with respect to the constant parameters of the potential function. The method provides three principal advantages. First, it obviates the problem of selecting the form of the functional dependence of the parameters upon the system's coordinates by employing a NN. If this network contains a sufficient number of neurons, it will automatically find something close to the best functional form. This is the case since Hornik et al., [Neural Networks 2, 359 (1989)] have shown that two-layer NNs with sigmoid transfer functions in the first hidden layer and linear functions in the output layer are universal approximators for analytic functions. Second, the entire fitting procedure is automated so that excellent fits are obtained rapidly with little human effort. Third, the method provides a procedure to avoid local minima in the multidimensional parameter hyperspace. As an illustrative example, the general method has been applied to the specific case of fitting the ab initio energies of Si(5) clusters that are observed in a molecular dynamics (MD) simulation of the machining of a silicon workpiece. The energies of the Si(5) configurations obtained in the MD calculations are computed using the B3LYP procedure with a 6-31G(**) basis set. The final ab initio database, which comprises the density functional theory energies of 10 202 Si(5) clusters, is fitted to an empirical Tersoff potential containing nine adjustable parameters, two of which are allowed to be the functions of the Si(5) configuration. The fitting error averaged over all 10 202 points is 0.0148 eV (1.43 kJ mol(-1)). This result is comparable to the accuracy achieved by more general fitting methods that do not rely on an assumed functional form for the potential surface.     

Ultrafast Decay of the Excited Singlet States of Thioxanthone by Internal Conversion and Intersystem Crossing

The experimental ultrafast photophysics of thioxanthone in several aprotic organic solvents at room temperature is presented, measured using femtosecond transient absorption together with high‐level ab initio CASPT2 calculations of the singlet‐ and triplet‐state manifolds in the gas phase, including computed state minima and conical intersections, transition energies, oscillator strengths, and spin–orbit coupling terms. The initially populated singlet ππ* state is shown to decay through internal conversion and intersystem crossing processes via intermediate nπ* singlet and triplet states, respectively. Two easily accessible conical intersections explain the favorable internal conversion rates and low fluorescence quantum yields in nonpolar media. The presence of a singlet–triplet crossing near the singlet ππ* minimum and the large spin–orbit coupling terms also rationalize the high intersystem crossing rates. A phenomenological kinetic scheme is proposed that accounts for the decrease in internal conversion and intersystem crossing (i.e. the very large experimental crescendo of the fluorescence quantum yield) with the increase of solvent polarity.     

Prediction of product branching ratios in the C(3P)+C2H2-->l-C3H+H/c-C3H+H/C3+H2 reaction using ab initio coupled clusters calculations extrapolated to the complete basis set combined with Rice-Ramsperger-Kassel-Marcus and radiationless transition theories

Ab initio CCSD(T) calculations of intermediates and transition states on the singlet and triplet C3H2 potential energy surfaces extrapolated to the complete basis set limit are combined with statistical computations of energy-dependent rate constants of the C(3P)+C2H2 reaction under crossed molecular beam conditions. Rice-Ramsperger-Kassel-Marcus theory is applied for isomerization and dissociation steps within the same multiplicity and radiationless transition and nonadiabatic transition state theories are used for singlet-triplet intersystem crossing rates. The calculated rate constants are utilized to predict product branching ratios. The results demonstrate that, in qualitative agreement with available experimental data, c-C3H+H and C3+H2 are the most probable products at low collision energies, whereas l-C3H+H becomes dominant at higher Ec above approximately 25 kJ/mol.     

Ab initio quantum chemical investigation of the first steps of the photocycle of phototropin: a model study

Phototropin is a blue light-activated photoreceptor that plays a dominant role in the phototropism of plants. The protein contains two subunits that bind flavin mononucleotide (FMN), which are responsible for the initial steps of the light-induced reaction. It has been proposed that the photoexcited flavin molecule adds a cysteine residue of the protein backbone, thus activating autophosphorylation of the enzyme. In this study, the electronic properties of several FMN-related compounds in different charge and spin states are characterized by means of ab initio quantum mechanical calculations. The model compounds serve as idealized model chromophores for phototropism. Reaction energies are estimated for simple model reactions, roughly representing the addition of a cysteine residue to the flavin molecule. Excitation energies were calculated with the help of time-dependent density functional theory. On the basis of these calculations we propose the following mechanism for the addition reaction: (1) after photoexcitation of FMN out of the singlet ground state S0, excited singlet state(s) are populated; these relax to the lowest excited singlet state S1, and subsequently by intersystem crossing FMN in the lowest triplet state, T1 is formed; (2) the triplet easily removes the neutral hydrogen atom from the H-S group of the cysteine residue; and (3) the resulting thio radical is added.     

A diabatic representation including both valence nonadiabatic interactions and spin-orbit effects for reaction dynamics

A diabatic representation is convenient in the study of electronically nonadiabatic chemical reactions because the diabatic energies and couplings are smooth functions of the nuclear coordinates and the couplings are scalar quantities. A method called the fourfold way was devised in our group to generate diabatic representations for spin-free electronic states. One drawback of diabatic states computed from the spin-free Hamiltonian, called a valence diabatic representation, for systems in which spin-orbit coupling cannot be ignored is that the couplings between the states are not zero in asymptotic regions, leading to difficulties in the calculation of reaction probabilities and other properties by semiclassical dynamics methods. Here we report an extension of the fourfold way to construct diabatic representations suitable for spin-coupled systems. In this article we formulate the method for the case of even-electron systems that yield pairs of fragments with doublet spin multiplicity. For this type of system, we introduce the further simplification of calculating the triplet diabatic energies in terms of the singlet diabatic energies via Slater's rules and assuming constant ratios of Coulomb to exchange integrals. Furthermore, the valence diabatic couplings in the triplet manifold are taken equal to the singlet ones. An important feature of the method is the introduction of scaling functions, as they allow one to deal with multibond reactions without having to include high-energy diabatic states. The global transformation matrix to the new diabatic representation, called the spin-valence diabatic representation, is constructed as the product of channel-specific transformation matrices, each one taken as the product of an asymptotic transformation matrix and a scaling function that depends on ratios of the spin-orbit splitting and the valence splittings. Thus the underlying basis functions are recoupled into suitable diabatic basis functions in a manner that provides a multibond generalization of the switch between Hund's cases in diatomic spectroscopy. The spin-orbit matrix elements in this representation are taken equal to their atomic values times a scaling function that depends on the internuclear distances. The spin-valence diabatic potential energy matrix is suitable for semiclassical dynamics simulations. Diagonalization of this matrix produces the spin-coupled adiabatic energies. For the sake of illustration, diabatic potential energy matrices are constructed along bond-fission coordinates for the HBr and the BrCH(2)Cl molecules. Comparison of the spin-coupled adiabatic energies obtained from the spin-valence diabatics with those obtained by ab initio calculations with geometry-dependent spin-orbit matrix elements shows that the new method is sufficiently accurate for practical purposes. The method formulated here should be most useful for systems with a large number of atoms, especially heavy atoms, and/or a large number of spin-coupled electronic states.     

High‐level ab initio study of the N+(3P)+SH2 reactions in the gas phase: Role of spin‐forbidden pathways

The singlet and triplet potential energy surfaces involved in N⁺+SH₂ reactions have been explored using high‐level ab initio techniques. The geometries of the stationary points were optimized at the QCISD/6‐311G(df,p) level. The final energies were obtained in CCSD(T)/6‐311+G(3df,2p) single‐point calculations. The results obtained show that, although the N⁺(¹D)+SH₂ entrance channel is higher in energy than the N⁺(³P)+SH₂ one, most of the [H₂, S, N]⁺ singlet state cations are lower in energy than the corresponding triplets, due to their different bonding characteristics. Both singlet and triplet potential energy surfaces are quite close each other, and crossover between them can occur. The minimum energy crossing points were located by means of CASSCF(6,5) calculations. The spin‐orbit couplings show that the transition probability from the triplet to the singlet potential energy surface is significantly large. One of the most important consequences is that some of the products of the reaction, such as SH⁺, can be formed in typical spin‐forbidden processes. Since all the relevant structures along these pathways are much lower in energy than the reactants, this mechanism should be accessible even at low impact energies and therefore could be important in processes taking place in interstellar media. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Ab initio potential-energy surfaces for complex, multichannel systems using modified novelty sampling and feedforward neural networks

A neural network/trajectory approach is presented for the development of accurate potential-energy hypersurfaces that can be utilized to conduct ab initio molecular dynamics (AIMD) and Monte Carlo studies of gas-phase chemical reactions, nanometric cutting, and nanotribology, and of a variety of mechanical properties of importance in potential microelectromechanical systems applications. The method is sufficiently robust that it can be applied to a wide range of polyatomic systems. The overall method integrates ab initio electronic structure calculations with importance sampling techniques that permit the critical regions of configuration space to be determined. The computed ab initio energies and gradients are then accurately interpolated using neural networks (NN) rather than arbitrary parametrized analytical functional forms, moving interpolation or least-squares methods. The sampling method involves a tight integration of molecular dynamics calculations with neural networks that employ early stopping and regularization procedures to improve network performance and test for convergence. The procedure can be initiated using an empirical potential surface or direct dynamics. The accuracy and interpolation power of the method has been tested for two cases, the global potential surface for vinyl bromide undergoing unimolecular decomposition via four different reaction channels and nanometric cutting of silicon. The results show that the sampling methods permit the important regions of configuration space to be easily and rapidly identified, that convergence of the NN fit to the ab initio electronic structure database can be easily monitored, and that the interpolation accuracy of the NN fits is excellent, even for systems involving five atoms or more. The method permits a substantial computational speed and accuracy advantage over existing methods, is robust, and relatively easy to implement.     

Cis-->trans, trans-->cis isomerizations and N-O bond dissociation of nitrous acid (HONO) on an ab initio potential surface obtained by novelty sampling and feed-forward neural network fitting

The isomerization and dissociation dynamics of HONO are investigated on an ab initio potential surface obtained by fitting the results of electronic structure calculations at 21 584 configurations by using previously described novelty sampling and feed-forward neural network (NN) methods. The electronic structure calculations are executed by using GAUSSIAN 98 with a 6-311G(d) basis set at the MP4(SDQ) level of accuracy. The average absolute error of the NN fits varies from 0.012 eV (1.22 kJ mol(-1)) to 0.017 eV (1.64 kJ mol(-1)). The average computation time for a HONO trajectory using a single NN surface is approximately 4.8 s. These computation times compare very favorably with those required by other methods primarily because the NN fitting needs to be executed only one time rather than at every integration point. If the average result obtained from a committee of NNs is employed at each point rather than a single NN, increased fitting accuracy can be achieved at the expense of increased computational requirements. In the present investigation, we find that a committee comprising five NN potentials reduces the average absolute interpolation error to 0.0111 eV (1.07 kJ mol(-1)). Cis-trans isomerization rates with total energy of 1.70 eV (including zero point energy) have been computed for a variety of different initial distributions of the internal energy. In contrast to results previously reported by using an empirical potential, where cis-->trans to trans-->cis rate coefficient ratios at 1.70 eV total energy were found to lie in the range of 2.0-12.9 depending on the vibration mode excited, these ratios on the ab initio NN potential lie in the range of 0.63-1.94. It is suggested that this result is a reflection of much larger intramode coupling terms present in the ab initio potential surface. A direct consequence of this increased coupling is a significant decrease in the mode specific rate enhancement when compared to results obtained by using empirical surfaces. All isomerizations are found to be first order in accordance with the results reported by using empirical potentials. The dissociation rate to NO+OH has been investigated at internal HONO energies of 3.10 and 3.30 eV for different distributions of this energy among the six vibrational modes of HONO. These dissociations are also found to be first order. The computed dissociation rate coefficients exhibit only modest mode specific rate enhancement that is significantly smaller than that obtained on an empirical surface because of the much larger mode couplings present on the ab initio surface.