Barrier heights, distribution

BSS recharging of which may accompany adsorption influences the adsorption-caused transformation of the type of barrier height distribution function. In this case, similarly to the situation which was addressed in section 1.8 one can easily obtain... 

Agmon, N. and Hopfield, J. (1983). CO binding to heme proteins a model for barrier height distributions and slow conformational changes. J. Chem. Phys. 79, 2042-2053... 

Fig. 4. Radial distribution functions between the centre of a test cavity and the (jxygen atom of the surrounding water. The curves correspond to the different barrier heights for the softcore interaction illustrated in Fig. 3...

Application of Eq. (30) corrects the free energies of the endpoints but not those of the intermediate conformations. Therefore, the above approach yields a free energy profile between qp and q-g, that is altered by the restraint(s). In particular, the barrier height is not that of the namral, unrestrained system. It is possible to correct the probability distributions P,. observed all along the pathway (with restraints) to obtain those of the unrestrained system [8,40]. Erom the relation P(q)Z, = P,(q)Z, cxp(UJkT) and Eqs. (6)-(8), one obtains... 

We have seen that the cooperative region, which represents a nominal dynamical unit of liquid, is of rather modest size, resulting in observable fluctuation effects. Xia and Wolynes [45] computed the relaxation barrier distribution. The configurational entropy must fluctuate, with the variance given by the usual expression [77] 5Sc) ) = Cp barrier height for a particular region is directly related to the local density of states, and hence to... 

For the sake of simplicity we can assume that initially the heights of barriers are distributed homogeneously over the interval A = 4o 5) ... 

Fig. 8.3. Histogram of work values for Jarzynski s identity applied to the double-well potential, V(x) = x2(x — a)2 + x, with harmonic guide Vpun(x, t) = k(x — vt)2/2, pulled with velocity v. Using skewed momenta, we can alter the work distribution to include more low-work trajectories. Langevin dynamics on Vtot(x(t),t) = V(x(t)) + Upuii(x(t)yt) with JcbT = 1, k = 100, was run with step size At = 0.001, and friction constant 7 = 0.2 (in arbitrary units). We choose v = 4 and a = 4, so that the barrier height is many times feT and the pulling speed far from reversible. Trajectories were run for a duration t = 1000. Work histograms for 10,000 trajectories, for both equilibrium (Maxwell) initial momenta, with zero average and unit variance, and a skewed distribution with zero average and a variance of 16.0...

Let us consider the case when the diffusion coefficient is small, or, more precisely, when the barrier height A is much larger than kT. As it turns out, one can obtain an analytic expression for the mean escape time in this limiting case, since then the probability current G over the barrier top near xmax is very small, so the probability density W(x,t) almost does not vary in time, representing quasi-stationary distribution. For this quasi-stationary state the small probability current G must be approximately independent of coordinate x and can be presented in the form... 

Isomerizations are important unimolecular reactions that result in the intramolecular rearrangement of atoms, and their rate parameters are of the same order of magnitude as other unimolecular reactions. Consequently, they can have significant impact on product distributions in high-temperature processes. A large number of different types of isomerization reactions seem to be possible, in which stable as well as radical species serve as reactants (Benson, 1976). Unfortunately, with the exception of cis-trans isomerizations, accurate kinetic information is scarce for many of these reactions. This is, in part, caused by experimental difficulties associated with the detection of isomers and with the presence of parallel reactions. However, with computational quantum mechanics theoretical estimations of barrier heights in isomerizations are now possible. 

To retain the analogy with a simple exponential function, it is considered in the cases described by Equations (1.8) and (1.9) that there is a distribution of barrier heights, g(G), each height corresponding to an exponential relaxation (Austin et al. 1975 Nagy et al. 2005). The concentration profile is in this case described by... 

If entropy changes are negligible, one nsnally speaks of enthalpic barrier heights that the reaction has to overcome. The distribntion of relaxation times or rate constants in this case may be eqnivalently considered a distribution of reaction barrier heights. 

One of the widely used methods of analysis of kinetic data is based on extraction of the distribution of relaxation times or, equivalently, enthalpic barrier heights. In this section, we show that this may be done easily by using the distribution function introduced by Raicu (1999 see Equation [1.16] above). To this end, we use the data reported by Walther and coworkers (Walther et al. 2005) from pump-probe as well as the transient phase grating measurements on trehalose-embedded MbCO. Their pump-probe data have been used without modification herein, while the phase grating data (also reproduced in Figure 1.12) have been corrected for thermal diffusion of the grating using the relaxation time reported above, r,, and Equation (1.25). 

Once the best-fit parameters are obtained, the distribution of relaxation times or, equivalently, barrier heights may be easily computed from Equation (1.16). The distribution function corresponding to the data in Figure 1.13 is plotted for the two types of measurements in Figure 1.14. 

---