Kinetic master equation

The fimdamental kinetic master equations for collisional energy redistribution follow the rules of the kinetic equations for all elementary reactions. Indeed an energy transfer process by inelastic collision, equation (A3.13.5). can be considered as a somewhat special reaction . The kinetic differential equations for these processes have been discussed in the general context of chapter A3.4 on gas kmetics. We discuss here some special aspects related to collisional energy transfer in reactive systems. The general master equation for relaxation and reaction is of the type [H, 12 and 13, 15, 25, 40, 4T ] ... 

In a typical sample with 10 °-coupled proton spins, a full quantum description is impractical and simplifications are needed. To this end, the polarization-transfer process in a many-spin system is modeled by a kinetic master equation for the polarizations p, = (S, ) ... 

The kinetic master equations associated with the reaction mechanism of Eq. (1) can be written as... 

As we have seen in the previous sectimis, the definition of a reasonably accurate kinetic master equation from the microscopic structure can be achieved by combining atomistic simulations, quantum chemical calculatirais, and microelectrostatic... 

Generalized first-order kinetics have been extensively reviewed in relation to teclmical chemical applications [59] and have been discussed in the context of copolymerization [53]. From a theoretical point of view, the general class of coupled kinetic equation (A3.4.138) and equation (A3.4.139) is important, because it allows for a general closed-fomi solution (in matrix fomi) [49]. Important applications include the Pauli master equation for statistical mechanical systems (in particular gas-phase statistical mechanical kinetics) [48] and the investigation of certain simple reaction systems [49, ]. It is the basis of the many-level treatment of... 

An important example for the application of general first-order kinetics in gas-phase reactions is the master equation treatment of the fall-off range of themial unimolecular reactions to describe non-equilibrium effects in the weak collision limit when activation and deactivation cross sections (equation (A3.4.125)) are to be retained in detail [ ]. 

A3.13.3.2 THE MASTER EQUATION FOR COLLISIONAL AND RADIATIVE ENERGY REDISTRIBUTION UNDER CONDITIONS OF GENERALIZED FIRST-ORDER KINETICS... 

Master equation methods are not tire only option for calculating tire kinetics of energy transfer and analytic approaches in general have certain drawbacks in not reflecting, for example, certain statistical aspects of coupled systems. Alternative approaches to tire calculation of energy migration dynamics in molecular ensembles are Monte Carlo calculations [18,19 and 20] and probability matrix iteration [21, 22], amongst otliers. 

Kinetic studies such as these use the master equation to follow the flow of probability between the states of the model. This equation is a basic loss-gain equation that describes the time evolution of the probability pi(t) for finding the system in state i [24]. The basic form of this equation is... 

Solving the master equation for the minimally frustrated random energy model showed that the kinetics depend on the connectivity [23]. Eor the globally connected model it was found that the resulting kinetics vary as a function of the energy gap between the folded and unfolded states and the roughness of the energy landscape. The model... 

Czerminski and Elber [64], who generated an almost complete map of the minima and barriers of an alanine tetrapeptide in vacuum. Using the master equation approach they were able to smdy aspects of this system s kinetics, which involve the crossing of barriers of different heights. 

If it cannot be guaranteed that the adsorbate remains in local equilibrium during its time evolution, then a set of macroscopic variables is not sufficient and an approach based on nonequihbrium statistical mechanics involving time-dependent distribution functions must be invoked. The kinetic lattice gas model is an example of such a theory [56]. It is derived from a Markovian master equation, but is not totally microscopic in that it is based on a phenomenological Hamiltonian. We demonstrate this approach... 

For adsorbates out of local equilibrium, an analytic approach to the kinetic lattice gas model is a powerful theoretical tool by which, in addition to numerical results, explicit formulas can be obtained to elucidate the underlying physics. This allows one to extract simplified pictures of and approximations to complicated processes, as shown above with precursor-mediated adsorption as an example. This task of theory is increasingly overlooked with the trend to using cheaper computer power for numerical simulations. Unfortunately, many of the simulations of adsorbate kinetics are based on unnecessarily oversimplified assumptions (for example, constant sticking coefficients, constant prefactors etc.) which rarely are spelled out because the physics has been introduced in terms of a set of computational instructions rather than formulating the theory rigorously, e.g., based on a master equation. 

Recently " we propo.sed to describe the nonequilibrium alloy kinetics using the fundamental master equation for probability P of finding the occupation number set... 

This reaction was investigated by Klippenstein and Harding [57] using multireference configuration interaction quantum chemistry (CAS + 1 + 2) to define the PES, variable reaction coordinate TST to determine microcanonical rate coefficients, and a one-dimensional (ID) master equation to evaluate the temperature and pressure dependence of the reaction kinetics. There are no experimental investigations of pathway branching in this reaction. 

In general, the equations for the density operator should be solved to describe the kinetics of the process. However, if the nondiagonal matrix elements of the density operator (with respect to electron states) do not play an essential role (or if they may be expressed through the diagonal matrix elements), the problem is reduced to the solution of the master equations for the diagonal matrix elements. Equations of two types may be considered. One of them is the equation for the reduced density matrix which is obtained after the calculation of the trace over the states of the nuclear subsystem. We will consider the other type of equation, which describes the change with time of the densities of the probability to find the system in a given electron state as a function of the coordinates of heavy particles Pt(R, q, Q, s,...) and Pf(R, q, ( , s,... ).74,77 80... 

The methodology of stochastic treatment of e-ion recombination kinetics is basically the same as for neutrals, except that the appropriate electrostatic field term must be included (see Sect. 7.3.1). This means the coulombic field in the dielectric for an isolated pair and, in the multiple ion-pair case, the field due to all unrecombined charges on each electron and ion. All the three methods of stochastic analysis—random flight Monte Carlo (MC), independent reaction time (IRT), and the master equation (ME)—have been used (Pimblott and Green, 1995). 

Immediately when the dynamic interpretation of Monte Carlo sampling in terms of the master equation, Eq. (31), was realized an application to study the critical divergence of the relaxation time in the two-dimensional Ising nearest-neighbor ferromagnet was attempted . For kinetic Ising and Potts models without any conservation laws, the consideration of dynamic universality classespredicts where z is the dynamic exponent , but the... 

The magnitude of the off-diagonal Hamiltonian (i.e. the energy transfer rate) thus depends on the strengths of the electron-phonon and Coulombic couplings and also the overlap of the two exciton wavefunctions[53]. Energy transfer rates from states to state jx, are calculated via the golden rule [54] and used as inputs to a master equation calculation of the excitation transfer kinetics in PSI, in which the dynamical information is included in the matrix K. 

An interesting, but probably incorrect, application of the probabilistic master equation is the description of chemical kinetics in a dilute gas.5 Instead of using the classical deterministic theory, several investigators have introduced single time functions of the form P(n1,n2,t) where P(nu n2, t) is the probability that there are nl particles of type 1 and n2 particles of type 2 in the system at time t. They use the transition rate A(nt, n2 n2, n2, t) from the state with particles of type 1 and n2 particles of type 2 to the state with nt and n2 particles of types 1 and 2, respectively, at time t. The rates that are used are obtained by assuming that only uncorrelated binary collisions occur in the system. These rates, however, are only correct in the thermodynamic limit for a low density system. In this limit, the Boltzmann equation is valid from which the deterministic theory follows. Thus, there is no reason to attach any physical significance to the differences between the results of the stochastic theory and the deterministic theory.6... 

The first one is that this particular form of H can also be used to prove the approach to equilibrium in the case of Boltzmann s kinetic equation for dilute gases. The Boltzmann equation is nonlinear and a different technique is needed to prove that all solutions tend to equilibrium. This technique is based on (5.6) other convex functions cannot be used. Incidentally, the Boltzmann equation is not a master equation for a probability density, but an evolution equation for the particle density in the six-dimensional one-particle phase space ( /i-space ). The linearized Boltzmann equation, however, has the same structure as a master equation (compare XIV.5). 

Equation (5.1) described the vibrational response of a single particle to an applied forceF(t). In a (crystalline) system of many mobile particles (ensemble), the problem is analogous but the question now is how the whole system responds to an external force or perturbation Let us define the system s state (a) as a particular configuration of its particles and the probability of this state as pa. In a thermodynamic system, transitions from an a to a p configuration occur as thermally activated events. If the transition frequency a- /5 is copa and depends only on a and / (Markovian), the time evolution of the system is given by a master equation which links atomic and macroscopic parameters (dynamics and kinetics)... 

One can derive Eqn. (12,12) in a more fundamental way by starting the statistical approach with the (Markovian) master equation, assuming that the jump probabilities obey Boltzmann statistics on the activation saddle points. Salje [E. Salje (1988)] has discussed the following general form of a kinetic equation for solid state processes... 

Since the formal chemical kinetics operates with large numbers of particles participating in reaction, they could be considered as continuous variables. However, taking into account the atomistic nature of defects, consider hereafter these numbers N as random integer variables. The chemical reaction can be treated now as the birth-death process with individual reaction events accompanied by creation and disappearance of several particles, in a line with the actual reaction scheme [16, 21, 27, 64, 65], Describing the state of a system by a vector N = TV),..., Ns, we can use the Chapmen-Kolmogorov master equation [27] for the distribution function P(N, t)... 

In terms of the master equation for the Markov process the formal kinetics is nothing but the mean-field theory where the fluctuation terms like that on the r.h.s. of equation (2.2.43) are neglected. Strictly speaking, the macroscopic description, equation (2.1.2), were correct if the fluctuation terms vanished as V —> oo. In a general case the function P(N, t) does not satisfy the Poisson distribution [16, 27] in particular, °N (N> ... 

Note here that the relation between mesoscopic and microscopic approaches is not trivial. In fact, the former is closer to the macroscopic treatment (Section 2.1.1) which neglects the structural characteristics of a system. Passing from the micro- to meso- and, finally, to macroscopic level we loose also the initial statement of a stochastic model of the Markov process. Indeed, the disadvantages of deterministic equations used for rather simplified treatment of bimolecular kinetics (Section 2.1) lead to the macro- and mesoscopic models (Section 2.2) where the stochasticity is kept either by adding the stochastic external forces (Section 2.2.1) or by postulating the master equation itself for the relevant Markov process (Section 2.2.2). In the former case the fluctuation source is assumed to be external, whereas in the latter kinetics of bimolecular reaction and fluctuations are coupled and mutually related. Section 2.3.1.2 is aimed to consider the relation between these three levels as well as to discuss problem of how determinicity and stochasticity can coexist. 

---