Correspondence of canonical and microcanonical rate constants using variational transition state theory for simple bond fissions

Diffusion in a two-dimensional energy landscape in the presence of dynamical correlations and validity of random walk model

Diffusion in a multidimensional energy surface with minima and barriers is a problem of importance in statistical mechanics and it also has wide applications, such as protein folding. To understand it in such a system, we carry out theory and simulations of a tagged particle moving on a two-dimensional periodic potential energy surface, both in the presence and absence of noise. Langevin dynamics simulations at multiple temperatures are carried out to obtain the diffusion coefficient of a solute particle. Friction is varied from zero to large values. Diffusive motion emerges in the limit of a long time, even in the absence of noise. Noise destroys the correlations and increases diffusion at small friction. Diffusion thus exhibits a nonmonotonic friction dependence at the intermediate value of the damping, ultimately converging to our theoretically predicted value. The latter is obtained using the well-established relationship between diffusion and random walk. An excellent agreement is obtained between theory and simulations in the high-friction limit but not so in the intermediate regime. We explain the deviation in the low- to intermediate-friction regime using the modified random walk theory. The rate of escape from one cell to another is obtained from the multidimensional rate theory of Langer. We find that enhanced dimensionality plays an important role. To quantify the effects of noise on the potential-imposed coherence on the trajectories, we calculate the Lyapunov exponent. At small friction values, the Lyapunov exponent mimics the friction dependence of the rate.

Efficient calculation of low energy statistical rates for gas phase dissociation using umbrella sampling

Monte Carlo (MC) simulations can be used to compute microcanonical statistical rates of gas phase dissociation reactions. Unfortunately, the MC approach may suffer from a slow convergence and large statistical errors for energies just above the dissociation threshold. In this work, umbrella sampling is proposed as a device to reduce the statistical error of MC rate constants. The method is tested by computing the classical dissociation rate for the reaction [H5O2+]* --> H2O + H3O(+) over the range of internal energy 38 < E < or = 100 kcal/mol. Comparing with other literature methods, it is found that umbrella sampling reduces the computational effort by up to two orders of magnitude when used in conjunction with a careful choice of sampling distributions. The comparison between MC rate constants and classical Rice-Ramsperberg-Kassel harmonic theory shows that anharmonicity plays an important role in the dissociation process of the Zundel cation (H5O2+) at all energies.

Monte Carlo-steepest descent (MC-SD) simulations of nanometric cutting

In order to reduce the computational time, Monte Carlo (MC) simulations of nanometric cutting have been modified to include a combination of steepest descent (SD) and Monte Carlo procedures. This MC-SD combination is found to reduce the required computational times by factors of at least two to three over those achieved using MC methods alone. The MC-SD method is applied to the nanometric cutting of single-crystal aluminium along the (100) plane with different rake angle tools at a cutting speed of 5m/s. The results obtained from the MC-SD calculations are found to be almost identical to those resulting from the MC simulations.

Evaluation of canonical and microcanonical nonadiabatic reaction rate constants by using the Zhu–Nakamura formulas

We consider a problem of calculating both thermal and microcanonical rate constants for nonadiabatic chemical reactions. Instead of using the conventional transition state theory, we use a generalized seam surface and introduce a concept of a coordinate dependent effective nonadiabatic transition probability based on the Zhu-Nakamura theory which can treat the nonadiabatic tunneling properly. The present approach can be combined with Monte Carlo method so as to be applicable to chemical reactions in complicated systems. The method is demonstrated to work well in wide energy and temperature range. Numerical tests also show that it is very essential for accurate evaluation of the thermal rate constant to use the generalized seam surface and take into account the nonadiabatic tunneling effect.