Computing rates reversible work

Stable states via small fluctuations. One can often infer the mechanism of a reaction by comparing stable and transition states. Transition rates can subsequently be calculated by computing the reversible work to reach the transition state from a stable state, and then initiating many fleeting trajectories from the transition state [3]. 

In an ergodic system, every possible trajectory of a particular duration occurs with a unique probability. This fact may be used to define a distribution functional for dynamical paths, upon which the statistical mechanics of trajectories is based. For example, with this functional one can construct partition functions for ensembles of trajectories satisfying specific constraints, and compute the reversible work to convert between these ensembles. In later sections, we will show that such manipulations may be used to compute transition rate constants. In this section, we derive the appropriate path distribution functionals for several types of microscopic dynamics, focusing on the constraint that paths are reactive, that is, that they begin in a particular stable state. A, and end in a different stable state, B. 

Umbrella sampling for ensembles of trajectories can be carried out in close analogy to the procedure described above. For the purpose of computing rate constants (in this case the rate of conformational transitions), we focus on computing the reversible work to confine trajectories endpoints to state B, given that these trajectories begin in state A. P (A, t) is defined to be the distribution of A at the endpoints of trajectories initiated in A,... 

Consideration of Eqs. (29) and (32) shows that the mechanical energy equation involves only the recoverable or reversible work. In order to calculate this term on the average, however, it is necessary to compute the total work done W and subtract from it the part lost due to friction or the irreversible work F. If Eq. (60) is applied to steady flow in a pipe and divided by the mass flow rate, the following per unit mass form is obtained,... 

The factorization of C(t) introduced in Section IV.D is reminiscent of the factorization used in conventional rate constant calculations [3]. In both cases, a reversible work calculation is needed to compute the correlation of... 

Brusatori and Van Tassel [20] presented a kinetic model of protein adsorption/surface-induced transition kinetics evaluated by the scale particle theory (SPT). Assuming that proteins (or, more generally, particles ) on the surface are at all times in an equilibrium distribution, they could express the probability functions that an incoming protein finds a space available for adsorption to the surface and an adsorbed protein has sufficient space to spread in terms of the reversible work required to create cavities in a binary system of reversibly and irreversibly adsorbed states. They foimd that the scale particle theory compared well with the computer simulation in the limit of a lower spreading rate (i.e., smaller surface-induced unfolding rate constant) and a relatively faster rate of surface filling. 

The earliest estimate of kB was by Johnston.224 On the basis of work he had performed on the N205 system,313 he computed a lower limit for k5 of 1010 Af-1 sec-1 at room temperature. Hisatsune, Crawford, and Ogg,202 using a rapid-scanning infrared spectrometer, studied the decomposition of N205 in the presence of NO. Relevant to the NO + N03 reaction is their determination of k6 and k6k5/k-6, where kg and k 6 are forward- and reverse-rate constants for the reaction... 

More recently Andrieux et. al. (5a,5b) have described a procedure for computer simulation of a second-order ec catalytic mechanism. In their work cyclic voltammetric data were calculated while changing the rate and reversibility of the follow-up reaction. Using the implicit finite-difference method... 

One cannot divorce the computational studies from all that has been done in analytic theory or in experiment (much of which predates the significant increase in the number of computational studies that occurred in the mid-1980s). We will therefore discuss some aspects of the analytic theories that shed light on the interaction between theory and simulation. A number of reviews have concentrated on analytic theories of chemical reactions and reaction rates in solution. In particular, we commend to the reader those of Hynes, Berne et al., and Hanggi et al. These reviews usually contain some discussion of computer simulations. However, here we reverse the priority and concentrate primarily on simulation. In addition, we will describe much of the work that has been done on how reactions climb barriers and what happens as they come off a barrier and return to equilibrium (or in the case of nonthermally activated reactions, how the energy placed into the reaction coordinate by outside means is dissipated into the solvent). Some of these areas have recently been discussed in a review by Ohmine and Sasai of the computer simulation of the dynamics of liquid water and this solvent s effect on chemical reactions. 

A typical set of the required five variables would be the temperatures and pressures of the inlet and outlet water streams, T, Tg, P, and Pg, plus the inlet water flow rate N. With values for these five variables, we can solve the steady-state material balance for the outlet water flow rate (the inlet and outlet mass flow rates are equal here) and we can solve the steady-state energy balance for Q. In this example the value computed for the heat duty is the actual value for the real process, regardless of reversibility, because the process is workfree. However, in the general case, when heat and work both cross a system boundary, the energy balance gives only their sum. Variations on this problem are also possible for example, if we knew values for the five variables T, Tg, P, Pg and Q, then we could solve the energy balance for the required water flow rate. Or, if we knew T, P, P , Q, and N, then we could solve for the outlet water temperature Tg. 

The kinetic propagation reaction equations for radical copolymerization of two monomers, Mj and M2, are written in Fig. 3.45. The complications due to additional, different monomers, as well as transfer, reverse, and termination reactions increase the numbers of equations and rate constants, so that the resulting reaction equations are unhandy, and it needs much experimental work to establish the rate constants. Also, the computational effort to solve the many equations becomes excessive. 

---