The transition-state approach

The posibility of extending MD to slower diffusion processes has been discussed (98). But applying such algorithms has a tradeoff on the overall quality of the computational approach. To perform calculations at time scales beyond those accessible to MD is possible nowadays only by using the transition state approach (TSA) proposed in (97,115,132). This method will be presented briefly below. 

To follow to actually carry out a TSA simulation a three-dimensional grid, with grid interval of about 0.2 A (5 -106 equispaced points in (132)) is built and the Helmholtz energies at all grid points are computed. Before this can be done in practice, a value for (A2) must be found. Then, local minima and the crest surfaces must be found, using the procedures given in (130,132,165). To study the dynamics of the penetrant molecules on the network of sites a Monte-Carlo procedure is employed, which is presented is some detail in (97). 

In Table 5-2 a comparison between diffusivities obtained with the TSA method and experimental D is presented. From this table one can see that, in all cases computed D agree with experimental data to within an order of magnitude. Moreover most of these D are considerably smaller than the 5 10-7 cm2/s lower threshold assumed to be in reach of nowadays MD simulations Section 5.2.1. This is an encouraging sign that computer simulations of diffusional processes are already able to predict, with a reasonable accuracy and for small and simple penetrants, diffusion coefficients around 10-10 cm2/s. From the point of view the packaging sector it would be interesting to learn if and when further theoretical developments of the TSA method will be able to simulate (predict) such slow diffusional processes for organic penetrants with a much more complex structure, see Chapter 3 and Appendix I. 

Two atomistic approaches have been presented briefly above molecular dynamics and the transition-state approach. They are still not ideal tools for the prediction of diffusion constants because (i) in order to obtain a reliable chain packing with a MD simulation one still needs the experimental density of the polymer and (ii) though TSA does not require classical dynamics it involves a number of simplifying assumptions, i.e. duration of jump mechanism, elastic polymer matrix, size of smearing factor, that impair to a certain degree the ab initio character of the method. However MD and TSA are valuable achievements, they are complementary in several 

Unfortunately, it seems that none of the diffusion models presented in the above sections meets completely these practical goals. 

---

The effect of external pressure on the rates of liquid phase reactions is normally quite small and, unless one goes to pressures of several hundred atmospheres, the effect is difficult to observe. In terms of the transition state approach to reactions in solution, the equilibrium existing between reactants and activated complexes may be analyzed in terms of Le Chatelier s principle or other theorems of moderation. The concentration of activated complex species (and hence the reaction rate) will be increased by an increase in hydrostatic pressure if the volume of the activated complex is less than the sum of the volumes of the reactant molecules. The rate of reaction will be decreased by an increase in external pressure if the volume of the activated complex molecules is greater than the sum of the volumes of the reactant molecules. For a decrease in external pressure, the opposite would be true. In most cases the rates of liquid phase reactions are enhanced by increased pressure, but there are also many cases where the converse situation prevails. 

The effect of pressure on the reaction rate constant can be interpreted by both the collision-, and the transition state or activated complex theories. However, it has generally been found that the role of pressure can be evaluated more clearly by the transition state approach [3]. 

Many reasons may be responsible for the different a )3 ratios observed in the four rings. In particular, the small ar.fi ratios for pyrrole [as well as the low sensitivity of this ring to substituent effects (Section IV, B)] may be due to the fact that in this case the Wheland intermediate is not a good model for the transition state. This hypothesis is in keeping with the Hammond postulate,183 according to which the transition state approaches closer to the unperturbed starting molecule as its reactivity increases. 

The transition-state approach permits us to make a separation of the factors constituting an experimental specific rate constant (for an elementary chemical act) into kinetic and thermodynamic factors. Thus, for the transition state X = A + B + C+ we can write for the rate constant jfc governing the appearance of products (Sec. XII.4) from X... 

Or, let us consider reactivity in the most precise way we know by the transition-state approach (Sec. 2.22). 

The transition state approach leads in a natural way to the Butler-Volmer equation, but is relatively weak in its predictive properties regarding the exchange current, ( 0, which is proportional to the frequency factor kr(., i and to exp(—AG J). The latter is quite closely related to the enthalpy and standard entropy of formation of the adsorbed reduction product or intermediate, and this is one main reason for the very intense modern efforts to develop predictive theoretical tools for the ab initio computation of adsorption energies at... 

The ability to predict accurate potential surfaces means that we are in a position to investigate the nature of the kinetic barriers, including the activated complexes, of key surface processes. In fact, Born-Oppenheimer potential surfaces can be used not only with the transition-state approach to kinetics but also with the much more general and exact collision theory (e.g., scattering S-matrix ilieory). While methods based on collision and scattering theory have pointed out deficiencies in the traditional transition-state theory (TST), they have also served to uphold many of TST s simple claims. In turn, new generalized transition state theories have been born. For complex systems, the transition-dale approach, while admittedly approximate, has been well established. 

By the end of 1972, a second cornerstone of the transition state approach was beginning to crumble significantly, for it was now quite evident that widely different transition states could be assumed for a given reaction, but the Rice-Ramsperger-Kassel-Marcus (RRKM) procedure would give the same result for the shape of the fall-off curve [72.N 72.R 74.F 79.A1]. This, as is now well known, arises through the adjustment of the model after the transition state has been chosen so as to force it to be consistent with the observed high pressure rate constant [72.R 80.P1], Perhaps it should have sounded the knell for the RRKM theory, much as the unsymmetric isotopic replacement experiments did for the Slater theory a decade earlier, but there was no other substitute available. 

Solvent effects on reaction rates - the transition state approach... 

For the calculation of the rate constants of such processes, the transition-state method has been proposed [129, 131, 365, 518]. Its relative simplicity (permitting calculation of the rate constants for many reactions) lies in that it does not attempt taking into account all the dynamical features of the elementary processes. It introduces instead the activated complex concept. However, it does not give unambiguous indication as to how the activated complex properties are connected with those of the reactants, thus leaving aside the dynamical problem. For this reason, the transition-state approach is sometimes opposed to the collision theory, though very often they are correlated. 

In (7.1.15) EcD is the diffusion activation barrier and Be - the pre-exponential factor depending on the tangential energy. The latter has been calculated within the theory of correlation functions (Doll and Voter 1985) and the transition state approach (Voter and Doll 1984) with the account of tunnel effects (Zhdanov 1985). 

Derive an expression for k(T)ER/k(T)LH, i.e., the ratio of the rate constants for the Eley-Rideal and the Langmuir-Hinshelwood rates for reaction of two atoms A + A A2. Assume the atom mass 1 amu, a temperature of 400 K, and vibrational frequencies 100 cm for the surface-bound atom. The gas pressure is 1 torr and the coverage 9a — atom/site. At what binding energy of the atom is the ratio equal to unity Hint Use the transition state approach (see Appendix A). 

A simple expression for the diffusion constant can be obtained using the transition state approach (see Appendix A). This theory gives the following expression for the rate constant for diffusion in the x-direction... 

---