Biomedical subjects

William L Hase

Publications and source records attributed to William L Hase.

At least 19 recordsLinked to original sources

PubMed2006

Ab initio and analytic intermolecular potentials for Ar-CH(3)OH.

Ab initio calculations at the CCSD(T)/aug-cc-pVTZ level of theory were used to characterize the Ar-CH(3)OH intermolecular potential energy surface (PES). Potential energy curves were calculated for four different Ar + CH(3)OH orientations and used to derive an analytic function for the intermolecular PES. A sum of Ar-C, Ar-O, Ar-H(C), and Ar-H(O) two-body potentials gives an excellent fit to these potential energy curves up to 100 kcal mol(-1), and adding an additional r(-n) term to the Buckingham two-body potential results in only a minor improvement in the fit. Three Ar-CH(3)OH van der Waals minima were found from the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVTZ calculations. The structure of the global minimum is in overall good agreement with experiment (X.-C. Tan, L. Sun and R. L. Kuczkowski, J. Mol. Spectrosc., 1995, 171, 248). It is T-shaped with the hydroxyl H-atom syn with respect to Ar. Extrapolated to the complete basis set (CBS) limit, the global minimum has a well depth of 0.72 kcal mol(-1) with basis set superposition error (BSSE) correction. The aug-cc-pVTZ basis set gives a well depth only 0.10 kcal mol(-1) smaller than this value. The well depths of the other two minima are within 0.16 kcal mol(-1) of the global minimum. The analytic Ar-CH(3)OH intermolecular potential also identifies these three minima as the only van der Waals minima and the structures predicted by the analytic potential are similar to the ab initio structures. The analytic potential identifies the same global minimum and the predicted well depths for the minima are within 0.05 kcal mol(-1) of the ab initio values. Combining this Ar-CH(3)OH intermolecular potential with a potential for a OH-terminated alkylthiolate self-assembled monolayer surface (i.e., HO-SAM) provides a potential to model Ar + HO-SAM collisions.

Argon↗

PubMed2006

Role of projectile and surface temperatures in the energy transfer dynamics of protonated peptide ion collisions with the diamond {111} surface.

The effects of temperature on energy transfer during collisions of protonated diglycine ions, Gly(2)-H(+), with a diamond {111} surface were investigated by chemical dynamics simulations. The simulations were performed for a collision energy of 70 eV and angle of 0 degrees with respect to the surface normal. In one set of simulations the initial surface temperature, T(surf), was varied from 300 to 2000 K, while the Gly(2)-H(+) vibrational and rotational temperatures were maintained at 300 K. For the second set of simulations the Gly(2)-H(+) vibrational temperature, T(vib), was varied from 300 to 2000 K, keeping both the Gly(2)-H(+) rotational and surface temperatures at 300 K. Increasing either the surface temperature or Gly(2)-H(+) vibrational temperature to values as high as 2000 K has, at most, only a negligible effect on the partitioning of the incident collision energy to the surface and to the vibrational and rotational modes of Gly(2)-H(+). To a good approximation, the initial surface and peptide ion energies are nearly adiabatic during the collisional energy transfer. This adiabaticity of the initial peptide ion energy agrees with experiments (J. Phys. Chem. A 2004, 108, 1). A more quantitative analysis of the effects of T(vib) and T(surf) shows there are small, but noticeable, effects on the energy transfer efficiencies. Namely, increasing the vibrational or surface temperature results in a near-linear decrease in the energy transfer to the degrees of freedom associated with this temperature.

Diamond↗

PubMed2006

Post-transition state dynamics for propene ozonolysis: Intramolecular and unimolecular dynamics of molozonide.

A direct chemical dynamics simulation, at the B3LYP6-31G(d) level of theory, was used to study the post-transition state intramolecular and unimolecular dynamics for the O3 + propene reaction. Comparisons of B3LYP6-31G(d) with CCSD(T)/cc-pVTZ and other levels of theory show that the former gives accurate structures and energies for the reaction's stationary points. The direct dynamics simulations are initiated at the anti and syn O3 + propene transition states (TSs) and the TS symmetries are preserved in forming the molozonide intermediates. Anti<-->syn molozonide isomerization has a very low barrier of 2-3 kcalmol and its Rice-Ramsperger-Kassel-Marcus (RRKM) lifetime is 0.3 ps. However, the trajectory isomerization is slower and it is unclear whether this anti<-->syn equilibration is complete when the trajectories are terminated at 1.6 ps. The syn (anti) molozonides dissociate to CH3CHO + H2COO and H2CO + syn (anti) CH3CHOO. The kinetics for the latter reactions are in overall good agreement with RRKM theory, but there is a symmetry preserving non-RRKM dynamical constraint for the former. Dissociation of anti molozonide to CH3CHO + H2COO is enhanced and suppressed, respectively, for the trajectory ensembles initiated at the anti and syn O3 + propene TSs. The dissociation of syn molozonide to CH3CHO + H2COO may also be enhanced for trajectories initiated at the syn O3 + propene TS. At the time the trajectories are terminated at 1.6 ps, the ratio of the trajectory and RRKM values of the CH3CHO + H2COO product yield is 1.6 if the symmetries of the initiation and dissociation TSs are the same and 0.6 if their symmetries are different. There are coherences in the intramolecular energy flow, which depend on molozonide's symmetry (i.e., anti or syn). This symmetry related dynamics is not completely understood, but it is clearly related to the non-RRKM dynamics for anti<-->syn isomerization and anti molozonide dissociation to CH3CHO + H2COO. Correlations are found between the stretching motions of molozonide, indicative of nonchaotic and non-RRKM dynamics. The non-RRKM dynamics of molozonide dissociation partitions vibration energy to H2COO that is larger than statistical partitioning. Though the direct dynamics simulations are classical, better agreement is obtained using quantum instead of classical harmonic RRKM theory. This may result from the neglect of anharmonicity in the RRKM calculations, the non-RRKM dynamics of the classical trajectories, or a combination of these two effects. The trajectories suggest that the equilibrium syn/anti molozonide ratio is approximately 1.1-1.2 times larger than that predicted by the harmonic densities of state, indicating an anharmonic correction.

Journal Article↗

PubMed2006

Molecular dynamics simulation of nanoparticle self-assembly at a liquid-liquid interface.

We have used molecular dynamics simulations to investigate the in situ self-assembly of modified hydrocarbon nanoparticles (mean diameter of 1.2 nm) at a water-trichloroethylene (TCE) interface. The nanoparticles were first distributed randomly in the water phase. The MD simulation shows the in situ formation of nanoparticle clusters and the migration of both single particles and clusters from the water phase to the trichloroethylene phase, possibly due to the hydrophobic nature of the nanoparticles. Eventually, the single nanoparticles or clusters equilibrate at the water-TCE interface, and the surrounding liquid molecules pack randomly when in contact with the nanoparticle surfaces. In addition, the simulations show that the water-TCE interfacial thickness analyzed from density profiles is influenced by the presence of nanoparticles either near or in contact with the interface but is independent of the number of nanoparticles present. The nanoparticles, water molecules, and TCE molecules all exhibit diffusion anisotropy.

Journal Article↗

PubMed2006

Dynamics of energy transfer in collisions of O(3P) atoms with a 1-decanethiol self-assembled monolayer surface.

Chemical dynamics simulations are reported of energy transfer in collisions of O(3P) atoms with a 300 K 1-decanethiol self-assembled monolayer (H-SAM) surface. The simulations are performed with a nonreactive potential energy surface, developed from PMP2/aug-cc-pVTZ calculations of the O(3P) + H-SAM intermolecular potential, and the simulation results represent the energy transfer dynamics in the absence of O(3P) reaction. Collisions energies E(i) of 0.12, 2.30, 11.2, 75.0, and 120.5 kcal/mol and incident angles theta(i) of 15, 30, 45, 60, and 75 degrees were considered in the study (theta(i) = 0 degrees is the surface normal). The translational energy distribution of the scattered O(3P) atoms, P(E(f)), may be deconvoluted into Boltzmann and non-Boltzmann components, with the former fraction identified as f(B). The trajectories are also analyzed in terms of three types; that is, direct scattering from and physisorption on the top of the H-SAM and penetration of the H-SAM. There are three energy regimes in the scattering dynamics. For the low E(i) values of 0.12 and 2.30 kcal/mol, physisorption is important and both f(B) and the average final translational energy of the scattered O(3P) atom, E(f), are nearly independent of the incident angle. The dynamics is much different for hyperthermal energies of 75.0 and 120.5 kcal/mol, where penetration of the surface is important. For hyperthermal collisions, the penetration probability decreases as theta(i) is increased, with a significant transition between theta(i) of 60 and 75 degrees . Hyperthermal penetration occurs upon initial surface impact and is more probable if the impinging O(3P) atom may move down a channel between the chains. For E(i) = 120.5 kcal/mol, 90% of the trajectories penetrate at theta(i) = 15 degrees , while only 3% penetrate at theta(i) = 75 degrees. For the former theta(i), the energy transfer to the surface is efficient with E(f) = 4.04 kcal/mol, but for the latter theta(i), E(f) = 85.3 kcal/mol! Particularly interesting penetrating trajectories are those in which O(3P) is trapped in the H-SAM for times exceeding 60 ps, linger near the Au substrate, and strike the Au substrate and scatter directly. For E(i) = 11.2 kcal/mol, there is a transition between the scattering dynamics for the low and hyperthermal collision energies. Additional detail in the energy transfer dynamics is obtained from the final polar and azimuthal angles, the residence time on/in the H-SAM, the minimum height with respect to the Au substrate, and the number of inner turning points in the O-atom's velocity. Calculated values of E(f) vs the final polar angle, theta(f), are in qualitative agreement with experiment. The O(3P) + H-SAM nonreactive energy transfer dynamics, for E(i) of 11.2 kcal/mol and lower, are very similar to previously reported Ne + H-SAM simulations.

Journal Article↗

PubMed2006

Ab initio and analytic intermolecular potentials for Ar-CF4.

Ab initio calculations at the CCSD(T) level of theory were performed to characterize the Ar + CF4 intermolecular potential. Potential energy curves were calculated with the aug-cc-pVTZ basis set, and with and without a correction for basis set superposition error (BSSE). Additional calculations were performed with other correlation consistent basis sets to extrapolate the Ar-CF4 potential energy minimum to the complete basis set (CBS) limit. Both the size of the basis set and BSSE have substantial effects on the Ar + CF4 potential. Calculations with the aug-cc-pVTZ basis set, and with a BSSE correction, appear to give a good representation of the BSSE corrected potential at the CBS limit. In addition, MP2 theory is found to give potential energies in very good agreement with those determined by the much higher level CCSD(T) theory. Two model analytic potential energy functions were determined for Ar + CF4. One is a fit to the aug-cc-pVTZ calculations with a BSSE correction. The second was derived by fitting an average BSSE corrected potential, which is an average of the CCSD(T)/aug-cc-pVTZ potentials with and without a BSSE correction. These analytic functions are written as a sum of two-body potentials and excellent fits to the ab initio potentials are obtained by representing each two-body interaction as a Buckingham potential.

Journal Article↗

PubMed2006

Use of a single trajectory to study product energy partitioning in unimolecular dissociation: mass effects for halogenated alkanes.

A single trajectory (ST) direct dynamics approach is compared with quasiclassical trajectory (QCT) direct dynamics calculations for determining product energy partitioning in unimolecular dissociation. Three comparisons are made by simulating C(2)H(5)F-->HF + C(2)H(4) product energy partitioning for the MP26-31G(*) and MP26-311 + + G(**) potential energy surfaces (PESs) and using the MP26-31G(*) PES for C(2)H(5)F dissociation as a model to simulate CHCl(2)CCl(3)-->HCl + C(2)Cl(4) dissociation and its product energy partitioning. The trajectories are initiated at the transition state with fixed energy in reaction-coordinate translation E(t) (double dagger). The QCT simulations have zero-point energy (ZPE) in the vibrational modes orthogonal to the reaction coordinate, while there is no ZPE for the STs. A semiquantitative agreement is obtained between the ST and QCT average percent product energy partitionings. The ST approach is used to study mass effects for product energy partitioning in HX(X = F or Cl) elimination from halogenated alkanes by using the MP26-31G(*) PES for C(2)H(5)F dissociation and varying the masses of the C, H, and F atoms. There is, at most, only a small mass effect for partitioning of energy to HX vibration and rotation. In contrast, there are substantial mass effects for partitioning to relative translation and the polyatomic product's vibration and rotation. If the center of mass of the polyatomic product is located away from the C atom from which HX recoils, the polyatomic has substantial rotation energy. Polyatomic products, with heavy atoms such as Cl atoms replacing the H atoms, receive substantial vibration energy that is primarily transferred to the wag-bend motions. For E(t) (double dagger) of 1.0 kcalmol, the ST calculations give average percent partitionings to relative translation, polyatomic vibration, polyatomic rotation, HX vibration, and HX rotation of 74.9%, 6.8%, 1.5%, 14.4%, and 2.4% for C(2)H(5)F dissociation and 39.7%, 38.1%, 0.2%, 16.1%, and 5.9% for a model of CHCl(2)CCl(3) dissociation.

Journal Article↗

PubMed2006

Comparison of levels of electronic structure theory in direct dynamics simulations of C2H5F --> HF + C2H4 product energy partitioning.

Direct dynamics simulations at the MP2/6-311++G** level of theory were performed to study C(2)H(5)F --> HF + C(2)H(4) product energy partitioning. The simulation results are compared with experiment and a previous MP2/6-31G* simulation. The current simulation with the larger basis set releases more energy to HF vibration and less to HF + C(2)H(4) relative translation as compared to the previous simulation with the 6-31G* basis set. The HF rotation and vibration energy distributions determined from the current simulation are in overall very good agreement with previous experimental studies of C(2)H(5)F dissociation by chemical activation and IRMPA. A comparison of the simulations with experiments suggests there may be important mass effects for energy partitioning in HX elimination from haloalkanes. The transition state (TS) structures and energies calculated with MP2 and the 6-31G* and 6-311++G** basis sets are compared with those calculated using CCD, CCSD, CCSD(T), and the 6-311++G** basis set.

Journal Article↗

PubMed2006

Dynamics and kinetics of heat transfer at the interface of model diamond [111] nanosurfaces.

A molecular dynamics simulation was performed to study the effect of an applied force on heat transfer at the interface of model diamond [111] nanosurfaces. The force was applied to a small, hot nanosurface at 800, 1000, or 1200 K brought into contact with a larger, colder nanosurface at 300 K. The relaxation of the initial nonequilibrium interfacial force occurs on a subpicosecond time scale, much shorter than that required for heat transfer. Heat transfer occurs with exponential kinetics and a rate constant that increases linearly with the interfacial force according to 7 x 10(-4) ps(-1)/nN. This rate constant only increases by at most 10% as the temperature of the hot surface is increased from 800 to 1200 K. Replacing the interfacial H-atoms on both surfaces by D atoms also has a very small effect on the heat transfer. However, if one nanosurface has H atoms on its interface and the other nanosurface's interface has D atoms, then there is a marked 25% decrease in the rate constant for heat transfer. Increasing the size of the hot surface, and, thus, the interfacial contact area, increases the rate of heat transfer but not the rate constant. For the same interfacial force, different anharmonic models for the nanosurfaces' potential energy function give the same heat transfer rate constant. The possibility of quantum effects for heat transfer across the diamond interface is considered.

Journal Article↗

PubMed2005

Direct dynamics trajectory study of the reaction of formaldehyde cation with D2: vibrational and zero-point energy effects on quasiclassical trajectories.

Quasiclassical, direct dynamics trajectories have been used to study the reaction of formaldehyde cation with molecular hydrogen, simulating the conditions in an experimental study of H2CO+ vibrational effects on this reaction. Effects of five different H2CO+ modes were probed, and we also examined different approaches to treating zero-point energy in quasiclassical trajectories. The calculated absolute cross-sections are in excellent agreement with experiments, and the results provide insight into the reaction mechanism, product scattering behavior, and energy disposal, and how they vary with impact parameter and reactant state. The reaction is sharply orientation-dependent, even at high collision energies, and both trajectories and experiment find that H2CO+ vibration inhibits reaction. On the other hand, the trajectories do not reproduce the anomalously strong effect of nu2(+) (the CO stretch). The origin of the discrepancy and approaches for minimizing such problems in quasiclassical trajectories are discussed.

Journal Article↗

PubMed2005

Reaction products with internal energy beyond the kinematic limit result from trajectories far from the minimum energy path: an example from H + HBr --> H2 + Br.

The importance of reactive trajectories straying far from the minimum energy path is demonstrated for the bimolecular reaction H + HBr --> H2(v', j') + Br at 53 kcal/mol collision energy. Product quantum state distributions are measured and calculated using the quasi-classical trajectory technique, and the calculations indicate that highly internally excited H2 products result from indirect reactive trajectories with bent transition states. A general argument is made suggesting that reaction products with internal energy exceeding a kinematic constraint can, in general, be attributed to reactive collisions straying far from the minimum energy path.

Journal Article↗

PubMed2005

Quasiclassical dynamics simulation of the collision-induced dissociation of Cr(CO)6 + with Xe.

Quasiclassical trajectory calculations are employed to investigate the dynamics of collision-induced dissociation (CID) of Cr(CO)6 + with Xe atoms at collision energies ranging from 1.3 to 5.0 eV. The trajectory simulations show that direct elimination of CO ligands, during the collision, becomes increasingly important as the collision energy increases. In a significant number of cases, this shattering mechanism is accompanied with a concomitant formation of a transient Xe-Cr(CO)x +(x<6) complex. The calculated results are in very good agreement with the experimental results presented previously [F. Muntean and P. B. Armentrout, J. Chem. Phys. 115, 1213 (2001)]. In particular, the computed cross sections and scattering maps for the product ions Cr(CO)x +(x=3-5) compare very favorably with the reported experimental data. However, in contrast with the conclusions of the previous study, the present calculations suggest that CID dynamics for this system exhibits a significant impulsive character rather than proceeding via a complex surviving more than a rotational period.

Journal Article↗

PubMed2005

Intermolecular potential to represent collisions of protonated peptide ions with fluorinated alkane surfaces.

The MP2/6-311++G(2df,2pd) level of theory was used to calculate intermolecular potential curves between CF(4), as a model for the C and F atoms of a fluorinated alkane surface, and CH(4), NH(3), NH(4)(+), H(2)CO, and H(2)O as models for different types of atoms and functional groups comprising protonated peptide ions. This level of theory was tested by comparisons with the MP2/aug-cc-pVTZ and CCSD(T)/aug-cc-pVTZ theories. Explicit-atom (EA) analytic potential energy functions were then derived by fitting these potential energy curves with two-body potentials between the atoms of the two interacting molecules. An intermolecular potential for the interaction of a protonated peptide ion with a fluorinated alkane surface may be constructed from these two-body potentials. Intermolecular potentials, for which CF(4) is treated as a united atom (UA), were developed by isotropically averaging the CF(4) orientation for each of the EA potential energy curves. The intermolecular potential energy curves calculated for CF(4) are compared with curves calculated previously for CH(4) interacting with the same molecules, to consider the relative efficiency of energy transfer for protonated peptide ion collisions with hydrogenated and fluorinated alkane surfaces.

Alkanes↗

PubMed2005

Nonequilibrium energy dissipation at the interface of sliding model hydroxylated alpha-alumina surfaces.

Nonequilibrium molecular dynamics simulations were performed to study the dynamics of energy transfer at the interface of a small nanoscale hydroxylated alpha-alumina surface sliding across a much larger surface of the same material. Sliding velocities of 0.05, 0.5, 5, and 50 ms and loads of 0, 0.0625, 5, 15, 25, and 100 nN were considered. Nonequilibrium energy distributions were found at the interface for each of these conditions. The velocity distribution P(v) for the atoms in a sublayer of the smaller surface oscillates during the sliding, reflecting the periodicity of the interfacial intermolecular potential. When averaged over the sliding, this P(v) for each of the sublayers is bimodal with Boltzmann and non-Boltzmann components. The non-Boltzmann component, with temperatures in excess of 1000 K and as high as 2500 K, is most important for the interfacial H-atom sublayer and becomes less important in moving to a sublayer further from the interface. Similarly, the temperature of the Boltzmann component decreases for sublayers further from the interface and approaches the 300 K temperature of the boundary. The temperature of the Boltzmann component decreases, but the importance of the non-Boltzmann component increases, as the sliding velocity is decreased. The temperature of the non-Boltzmann component is relatively insensitive to the sliding velocity. Friction forces are determined by calculating the energy dissipation during the sliding, and different regimes are found for variation in the friction force versus sliding velocity v(s) and applied load. For v(s) of 0.05, 0.5, and 5 ms, the friction force is inversely proportional to v(s) reflecting the increased time for energy dissipation as v(s) is decreased.

Journal Article↗

PubMed2005

Effect of the Ar-Ni(s) potential on the cross section for Ar+CH4/Ni{111} collision-induced desorption and the need for a more accurate CH4/Ni{111} potential.

In a previous paper [L. Sun, P. de Sainte Claire, O. Meroueh, and W. L Hase, J. Chem. Phys. 114, 535 (2001)], a classical trajectory simulation was reported of CH(4) desorption from Ni{111} by Ar-atom collisions. At an incident angle theta(i) of 60 degrees (with respect to the surface normal), the calculated collision-induced desorption (CID) cross sections are in excellent agreement with experiment. However, for smaller incident angles the calculated cross sections are larger than the experimental values and for normal collisions, theta(i)=0 degrees , the calculated cross sections are approximately a factor of 2 larger. This trajectory study used an analytic function for the Ar+Ni(s) intermolecular potential which gives an Ar-Ni{111} potential energy minimum which is an order of magnitude too deep. In the work reported here, the previous trajectory study is repeated with an Ar+Ni(s) analytic intermolecular potential which gives an accurate Ar-Ni{111} potential energy minimum and also has a different surface corrugation than the previous potential. Though there are significant differences between the two Ar+Ni(s) analytic potentials, they have no important effects on the CID dynamics and the cross sections reported here are nearly identical to the previous values. Zero-point energy motions of the surface and the CH(4)-Ni(s) intermolecular modes are considered in the simulation and they are found to have a negligible effect on the CID cross sections. Calculations of the intermolecular potential between CH(4) and a Ni atom, at various levels of theory, suggest that there are substantial approximations in the ab initio calculation used to develop the CH(4)+Ni{111} potential. The implication is that the differences between the trajectory and experimental CID cross sections may arise from an inaccurate CH(4)+Ni{111} potential used in the trajectory simulation.

Journal Article↗

PubMed2004

Ab initio direct dynamics trajectory simulation of C2H5F-->C2H4 + HF product energy partitioning.

Direct dynamics classical trajectory simulations were performed to study product energy partitioning in C(2)H(5)F-->C(2)H(4)+HF dissociation. The intrinsic reaction coordinate potential energy curve, reaction energetics, and transition state (TS) properties were calculated for this reaction at different levels of electronic structure theory, and MP2/6-31G( *) was chosen as a meaningful and practical method for performing the direct dynamics. The trajectories show that the HF bond, uncoupled from the other degrees of freedom, is formed within the first 10 fs as the system moves from the TS towards products. The populations of the HF vibration states, determined from the simulations, decrease monotonically as found from experiments. However, the simulation's populations for the low and high energy vibration states are larger and smaller, respectively, than the experimental results. The HF rotational temperature found from the simulations is in agreement with experiment. Increasing the TS's excess energy gives higher rotational temperatures for both C(2)H(4) and HF. Energy is partitioned to the products from both the excess energy in the TS and the potential energy release in the exit channel. Partitioning from these two energy sources is distinguished by varying the TS's excess energy. An analysis of the simulation's energy disposal shows that the fractions of the excess energy partitioned to relative translation, C(2)H(4) vibration, C(2)H(4) rotation, HF vibration, and HF rotation, are 0.17, 0.64, 0.076, 0.067, and 0.046, respectively, and are in good agreement with previous simulations on empirical potentials and experiments. The partitioning found for the potential energy release is 81%, <0.05%, 5%, 11%, and 3% to relative translation, C(2)H(4) vibration, C(2)H(4) rotation, HF vibration, and HF rotation. This result is substantially different than the deduction from experiments, which summarizes the partitioning as 20%, 45%, 24%, and <12% to relative translation, C(2)H(4) vibration+rotation, HF vibration, and HF rotation. Possible origins of the difference between the simulations and experiments in the release of the potential energy is discussed.

Journal Article↗

PubMed2004

Direct dynamics trajectory study of vibrational effects: can polanyi rules be generalized to a polyatomic system?

We present an ab initio direct dynamics trajectory study of the hydrogen abstraction reaction: H2CO+ + CD4 --> H2COD+ + CD3, with methane excited in two different distortion modes (nu4 and nu2). The trajectory simulations were able to reproduce experimental results and for the first time show how vibrational enhancement originates in reaction of small polyatomic species. Roughly equal contributions from two vibrational enhancement mechanisms were found. The "distortion" mechanism correlates the vibrational effects with vibration-induced reactant distortions, and the "velocity" mechanism correlates vibrational effects with vibrational velocities of the constituent atoms. This reaction has a reactant-like transition state and, thus, would correspond to an "early" barrier system in the context of the well-known Polanyi rules for predicting effects of vibration and collision energy. Straightforward application of these rules would predict that vibration should be ineffective in driving reaction, in disagreement with both experiment and trajectory results. Using the trajectories for guidance, we were able to construct a two-dimensional cut through the reaction potential energy surface that does suggest a predictive, Polanyi-type rule.

Journal Article↗

PubMed2004

Energetics, transition states, and intrinsic reaction coordinates for reactions associated with O(3P) processing of hydrocarbon materials.

Electronic structure calculations based on multiconfiguration wave functions are used to investigate a set of archetypal reactions relevant to O(3P) processing of hydrocarbon molecules and surfaces. These include O(3P) reactions with methane and ethane to give OH plus methyl or ethyl radicals, O(3P) + ethane to give CH3O + CH3, and secondary reactions of the OH product radical with ethane and the ethyl radical. Geometry optimization is carried out with CASSCF/cc-pVTZ for all reactions, and with CASPT2/cc-pVTZ for O(3P) + methane/ethane. Single-point energy corrections are applied with CASPT2, CASPT3, and MRCI + Q with the cc-pVTZ and cc-pVQZ basis sets, and the energies extrapolated to the complete basis set limit (CBL). Where comparison of computed barriers and energies of reaction with experiment is possible, the agreement is good to excellent. The best agreement (within experimental error) is found for MRCI + Q/CBL applied to O(3P) + methane. For the other reactions, CASPT2/CBL and MRCI + Q/CBL predictions differ from experiment by 1-5 kcal/mol for 0 K enthalpies of reaction, and are within 1 kcal/mol of the best-estimate experimental range of 0 K barriers for O(3P) + ethane and OH + ethane. The accuracy of MRCI + Q/CBL is limited mainly by the quality of the active space. CASPT2/CBL barriers are consistently lower than MRCI + Q/CBL barriers with identical reference spaces.

Journal Article↗