Abstract. The self-consistent decay of mixing (SCDM) semiclassical trajectory method for electronically nonadiabatic dynamics is improved by modifying the switching probability that determines the instantaneous electronic state toward which the system decoheres. The new method is called coherent switching with decay of mixing (CSDM), and it differs from the previously presented SCDM method in that the electronic amplitudes controlling the switching of the decoherent state are treated fully coherently in the electronic equations of motion for each complete passage through a strong interaction region. The new method is tested against accurate quantum mechanical calculations for twelve atom-diatom scattering test cases. Also tested are the SCDM method and the trajectory surface hopping method of Parlant and Gislason that requires coherent passages through each strong interaction region, and which we call the ECP-TSH method. The results are compared with previously presented results for the fewest-switches with time uncertainty and Tully's fewest switches (TFS) surface hopping methods and the semiclassical Ehrenfest method. We find that the CSDM method is the most accurate of the semiclassical trajectory methods tested. Including coherent passages improves the accuracy of the SCDM method (i.e., the CSDM method is more accurate than the SCDM method) but not of the trajectory surface hopping method (i.e., the ECP-TSH method is not more accurate on average than the TFS method).

The CH3 + OH bimolecular reaction and the dissociation of methanol are studied theoretically at conditions relevant to combustion chemistry. Kinetics for the CH3 + OH barrierless association reaction and for the H + CH2OH and H + CH3O product channels are determined in the high-pressure limit using variable reaction coordinate transition state theory and multireference electronic structure calculations to evaluate the fragment interaction energies. The CH3 + OH --> 3CH2 + H2O abstraction reaction and the H2 + HCOH and H2 + H2CO product channels feature localized dynamical bottlenecks and are treated using variational transition state theory and QCISD(T) energies extrapolated to the complete basis set limit. The 1CH2 + H2O product channel has two dynamical regimes, featuring both an inner saddle point and an outer barrierless region, and it is shown that a microcanonical two-state model is necessary to properly describe the association rate for this reaction over a broad temperature range. Experimental channel energies for the methanol system are reevaluated using the Active Thermochemical Tables (ATcT) approach. Pressure dependent, phenomenological rate coefficients for the CH3 + OH bimolecular reaction and for methanol decomposition are determined via master equation simulations. The predicted results agree well with experimental results, including those from a companion high-temperature shock tube determination for the decomposition of methanol.

Stationary points on the C(3)H(7)O potential energy surface relevant to the title reaction are calculated employing RQCISD(T)/cc-pVinfinityZ//B3LYP/6-311++G(d,p) quantum chemical calculations. Rate coefficients at 50-3000 K temperature and from zero to infinite pressure are calculated using an RRKM-based multiwell master equation. Due to the topography of the entrance channel an effective two-transition-state model is used to calculate accurate association rate coefficients. Our calculations are in excellent agreement with the available experimental data. We predict approximately 5% vinyl alcohol branching above 1000 K, the allyl radical formation being the main channel at high temperatures.

The ability to predict the pressure dependence of chemical reaction rates would be a great boon to kinetic modeling of processes such as combustion and atmospheric chemistry. This pressure dependence is intimately related to the rate of collision-induced transitions in energy E and angular momentum J. We present a scheme for predicting this pressure dependence based on coupling trajectory-based determinations of moments of the E,J-resolved collisional transfer rates with the two-dimensional master equation. This completely a priori procedure provides a means for proceeding beyond the empiricism of prior work. The requisite microcanonical dissociation rates are obtained from ab initio transition state theory. Predictions for the CH4 = CH3 + H and C2H3 = C2H2 + H reaction systems are in excellent agreement with experiment.

A semiclassical trajectory method, called the self-consistent decay of mixing (SCDM) method, is presented for the treatment of electronically nonadiabatic dynamics. The SCDM method is a modification of the semiclassical Ehrenfest (SE) method (also called the semiclassical time-dependent self-consistent-field method) that solves the problem of unphysical mixed final states by including decay-of-mixing terms in the equations for the evolution of the electronic state populations. These terms generate a force, called the decoherent force (or dephasing force), that drives the electronic component of each trajectory toward a pure state. Results for several mixed quantum-classical methods, in particular the SCDM, SE, and natural-decay-of-mixing methods and several trajectory surface hopping methods, are compared to the results of accurate quantum mechanical calculations for 12 cases involving five different fully dimensional triatomic model systems. The SCDM method is found to be the most accurate of the methods tested. The method should be useful for the simulation of photochemical reactions.