Evaluation of the Pre-Exponential Factor

Taking the malonaldehyde molecule as an example, we explain how to efficiently evaluate the pre-exponential factor B given by [see Equation (6.100)] 

FIGURE 6.1 (a) One-dimensional interpolation of the potential V[qo(z)] as a function of the parameter z along the instanton path qo(z) for three ab initio methods in the case of malon-aldehyde. The points represent the ab initio data. The lines are obtained by piecewise cubic interpolation, (b) One-dimensional interpolation of the elements of Hessian 9 V [q(z)]/9 along the instanton path for MP2/cc-pVDZ ab initio method. Two examples are shown. (Taken from Reference [122] with permission.) 

This BF frame is defined by the two conditions. The first one is to put the origin in the center of mass of the molecule. 

FIGURE 6.2 Malonaldehyde molecule. (Taken from Reference [123] with permission.) 

Construction of the Hamiltonian for zero total angular momentum (7 = 0) is now straightforward. We consider the kinetic metric 

---

In order to construct the expression for the equilibrium number of nuclei in a unit volume (the dimension of b(x)dx is cm-3, the dimension of b(x), when x is defined as the radius, is cm-2), we must multiply the exponent exp(- /fcT), where is determined by (17), by a quantity of dimension cm-2. Exact evaluation of a pre-exponential factor is presently an unsolved problem of statistical mechanics. Erom dimensional considerations we may propose d 2 or x 2, where d is the linear size of a molecule of liquid and x is the radius of a bubble. In the present problem of evaluating the critical (i.e., minimum) value of the equilibrium concentration, we are dealing with a region where the factor in the exponent is large and exact evaluation of the pre-exponential factor is not actually necessary. 

Briefly, in the Eyring theory for any reaction, homogeneous or heterogeneous, there is an evaluation of the pre-exponential factor A in the Arrhenius equation... 

An increase in the number of ways to store energy increases the entropy of a system. Thus, an estimate of the pre-exponential factor A in TST requires an estimate of the ratio g /gr. A common approximation in evaluating a partition function is to separate it into contributions from the various modes of energy storage, translational (tr), rotational (rot), and vibrational (vib) ... 

An exception to this approach is to be found in the work of Cohen [44] who has used the transition state methods developed by Benson and his colleagues [38] as a framework for evaluating experimental data for bimolecular metathetical reactions such as those of H, O and OH with alkanes. Using a model of the transition state, a theoretical value of the pre-exponential factor in the rate expression is derived which may be combined with experimental data at one temperature to give the exponential term. The rate expression so derived may be used to calculate values at other temperatures. By treating families of reactions adjustments can be made to the transition states to make them compatible with the experimental data on the whole range of reactions considered. 

To evaluate the Al-dimensional analog of the pre-exponential factor B [see Equation (2.135)1, we consider the path integral over harmonic fluctuations along the instanton. 

In most uncertainty studies published so far (see e.g. Brown et al. (1999), Turanyi et al. (2002), Zsely et al. (2005), Zador et al. (2005a, b, 2006a) and Zsely et al. (2008)), where the uncertainties of the rate coefficients were utilised, the uncertainty of k was considered to be equal to the xmcertainty of the pre-exponential factor A. This implies that the uncertainty of parameters E and n is zero, which is an unrealistic assumption. In a global sensitivity analysis study of a turbulent reacting atmospheric plume, Ziehn et al. (2009a) demonstrated the importance of uncertainties in EIR for the reaction N0 + 03 = N02 + 02. In this case for the prediction of mean plume centre line O3 concentratiOTis, the sensitivity to the assumed value for EIR was almost a factor of 20 higher than that of the A-factor, based on input parameter uncertainty factors provided by the evaluation of Androulakis (2004, 2004). However, in this case the parameters of the Anhenius expression for the chemical reactions considered were allowed to vary independently. In fact, the characterisation of the joint uncertainty of the Arrhenius parameters is important for the reahstic calculatiOTi of the uncertainty of chemical kinetic simulation results as will be discussed in the next section. 

The value of reaction rate Eq. (43) can be negative when N02 present in the mixture is transformed to NO via backward reaction, typically at higher temperatures. A comparison of measured and simulated outlet N02 concentrations in dependence on temperature can be seen for two different space velocities in Fig. 13. The pre-exponential factor k j and activation energy Ej of the kinetic constant no/no2 in the global rate law were evaluated by the weighted least squares method, Eq. (35). 

For the fresh and the specifically aged catalyst materials, the dependence of the normalized NOx storage capacity on temperature could be kept the same (Giithenke et al, 2007b). This minimized the number of parameters to be re-adapted for two catalysts with different ageing level. Thus, only the maximum NOx storage capacity and the pre-exponential factors for the reactions R1-R22 had to be re-evaluated, cf. Table III and Eq. (36). 

As discussed in Chapter 1 (Sections III and TV), the kinetics of drug degradation has been the topic of numerous books and articles. The Arrhenius relationship is probably the most commonly used expression for evaluating the relationship between rates of reaction and temperature for a given order of reaction (For a more thorough treatment of the Arrhenius equation and prediction of chemical stability, see Ref. 13). If the decomposition of a drug obeys the Arrhenius relationship [i.e., k = A exp(—Ea/RT), where k is the degree of rate constant, A is the pre-exponential factor ... 

The two-site reaction kinetics model proposed by Bonn [1] was used to evaluate the kinetic parameters. Activation energies and pre-exponential factors were determined from experiments between 570-630 K at 10 MPa. In order to decrease the strong inter-correlation between pre-exponential factors and activation energies, the reparametrisadon method of Kittrell [4] was used. Values for the pre-exponential factors at a reference temperature and activation energies are presented in Table 2. Experimental and theoretical details on HDM reaction kinetics will be published elsewhere [5]. 

Let us consider the last statement in more detail. Keeping in mind serious doubts stated in Section III.B, we nevertheless assume that pre-exponential factors for rate constants of collision-type reactions between gas species and surface sites can be evaluated as a gas-surface collision frequency (v). The dimension (or cross-section) of the surface sites (surface area (3), as well as a steric factor (c) should also be taken into account. The latter value shows the probability of reciprocal orientation of two reactants optimal for the reaction to proceed during the collision. As a result, the pre-exponential factor can be calculated as follows ... 

The evaluation of the desorption energy for multilayer desorption, i.e. the heat of evaporation, is straightforward. The plot of In 7 des vs. 1 IT for the leading edge of the spectrum yields a straight line, where the slope is equal to -EdeJk. If the desorption rate is known quantitatively then the intercept of this straight line with the 7-axis yields the pre-exponential factor. 

The deviation in the pre-exponential factor, represented here by the Damkoehler number Da , seems to be quite severe. But this is a superficial judgement only. For numbers of this magnitude, deviations in percentages are always very large although the realistic assessment comes up with moderate differences. It becomes easier to evaluate the quality of the data determined, if they are used in a new simulation and if this result is compared with the original simulation. This is presented in Figures 4-76 to 4-78. 

---