PubMed Health⌕ Search

Biomedical subjects

N J M Horing

Publications and source records attributed to N J M Horing.

2 recordsLinked to original sources

Transition rate prefactors for systems of many degrees of freedom.

When a minimum on the potential energy surface is surrounded by multiple saddle points with similar energy barriers, the transition pathways with greater prefactors are more important than those that have similar energy barriers but smaller prefactors. In this paper, we present a theoretical formulation for the prefactors, computing the probabilities for transition paths from a minimum to its surrounding saddle points. We apply this formulation to a system of 2 degrees of freedom and a system of 14 degrees of freedom. The first is Brownian motion in a two-dimensional potential whose global anharmonicities play a dominant role in determining the transition rates. The second is a Lennard-Jones (LJ) cluster of seven particles in two dimensions. Low lying transition states of the LJ cluster, which can be reached directly from a minimum without passing through another minimum, are identified without any presumption of their characteristics nor of the product states they lead to. The probabilities are computed for paths going from an equilibrium ensemble of states near a given minimum to the surrounding transition states. These probabilities are directly related to the prefactors in the rate formula. This determination of the rate prefactors includes all anharmonicities, near or far from transition states, which are pertinent in the very sophisticated energy landscape of LJ clusters and in many other complex systems.

Algorithms↗

Efficient transition path sampling for systems with multiple reaction pathways.

A new algorithm is developed for sampling transition paths and computing reaction rates. To illustrate the use of this method, we study a two-dimensional system that has two reaction pathways: one pathway is straight with a relatively high barrier and the other is roundabout with a lower barrier. The transition rate and the ratio between the numbers of the straight and roundabout transition paths are computed for a wide range of temperatures. Our study shows that the harmonic approximation for fluctuations about the steepest-descent paths is not valid even at relatively low temperatures and, furthermore, that factors related to entropy have to be determined by the global geometry of the potential-energy surface (rather than just the local curvatures alone) for complex reaction systems. It is reasonable to expect that this algorithm is also applicable to higher dimensional systems.

Algorithms↗