On the derivation of semiclassical expressions for quantum reaction rate constants in multidimensional systems

Entanglement Effect and Angular Momentum Conservation in a Nonseparable Tunneling Treatment

The important, and often dominant, role of tunneling in low temperature kinetics has resulted in numerous theoretical explorations into the methodology for predicting it. Nevertheless, there are still key aspects of the derivations that are lacking, particularly for nonseparable systems in the low temperature regime, and further explorations of the physical factors affecting the tunneling rate are warranted. In this work we obtain a closed-form rate expression for the tunneling rate constant that is a direct analog of the rigid-rotor-harmonic-oscillator expression. This expression introduces a novel "entanglement factor" that modulates the reaction rate. Furthermore, we are able to extend this expression, which is valid for nonseparable systems at low temperatures, to properly account for the conservation of angular momentum. In contrast, previous calculations have considered only vibrational transverse modes and so effectively employ a decoupled rotational partition function for the orientational modes. We also suggest a simple theoretical model to describe the tunneling effects in the vicinity of the crossover temperature (the temperature where tunneling becomes the dominating mechanism). This model allows one to naturally classify, interpret, and predict experimental data. Among other things, it quantitatively explains in simple terms the so-called "quantum bobsled" effect, also known as the negative centrifugal effect, which is related to curvature of the reaction path. Taken together, the expressions obtained here allow one to predict the thermal and E-resolved rate constants over broad ranges of temperatures and energies.

Development of semiclassical molecular dynamics simulation method

Various quantum mechanical effects such as nonadiabatic transitions, quantum mechanical tunneling and coherence play crucial roles in a variety of chemical and biological systems.

Microcanonical and thermal instanton rate theory for chemical reactions at all temperatures

Semiclassical instanton theory is used to study the quantum effects of tunnelling and delocalization in molecular systems. An analysis of the approximations involved in the method is presented based on a recent first-principles derivation of instanton rate theory [J. Chem. Phys., 2016,144, 114106]. It is known that the standard instanton method is unable to accurately compute thermal rates near the crossover temperature. The causes of this problem are identified and an improved method is proposed, whereby an instanton approximation to the microcanonical rate is defined and integrated numerically to obtain a thermal rate at any temperature. No new computational algorithms are required, but only data analysis of a number of standard instanton calculations.

The direct and precursor mediated dissociation rates of H2 on a Ni(111) surface

The dissociation and recombination rates of physisorbed H2, and the direct and steady state dissociation (i.e., the precursor mediated dissociation) rates of gas phase H2 on Ni(111), as well as the corresponding kinetic isotope effects, are calculated using the quantum instanton method, together with path integral Monte Carlo and adaptive umbrella sampling techniques. All these rates except the recombination one first decrease and then increase with the increasing temperature, and their minimum values appear at about 250, 300 and 250 K, respectively. These non-monotonic behaviors reveal that the quantum effect of H2 should be very remarkable at low temperatures. The steady state rates are smaller than the direct rates at low temperatures, however, they become larger than the direct ones at high temperatures, these two kinds of rates become equal at about 400 and 300 K on the rigid and quantum lattices, respectively. The quantum motion of the lattice can enhance the direct and steady state rates, and it increases the steady state rate much more than the direct one, for instance, the direct and steady state rates on the quantum lattice are 1.30 and 2.08 times larger than that on the rigid one at 300 K. The calculated kinetic isotope effects are much larger than 1, which reveals that H2 always has a larger rate than that of D2, and the direct process predicts much larger kinetic isotope effects than the steady state process at low temperatures. In addition, the kinetic isotope effects are not affected by the lattice motion.