Reduced systems initial value problem

Equation (4.2) is to be solved as an initial value problem under the conditions Cl(—00) = 0 and C2(—oo) = 1, so that the nonadiabatic transition probability (or, equivalently, the probability of the system remains as the diabatic state 2) is given by P = C2(- -oo)p. Eliminating C2, the problem is reduced to a second order differential equation... 

Reduced systems modelling based on the initial value problem in Eq. (7.88) requires the application of three functions. Function d=g(a) defines the time derivative of a, function Y = h(a) calculates the concentrations from the parameters of the manifold (mapping whilst function a = h(Y) (mapping... 

The number of equations, M5C + 1), for a large number of trays and components, can be excessive. The global Newton method will suffer from the same problem of requiring initial values near the answer. This problem is aggravated with nonequilibrium models because of difficulties due to nonideal if-values and enthalpies then compounded by the addition of mass transfer coefficients to the thermodynamic properties and by the large number of equations. Taylor et al. (80) found that the number of sections of packing does not have to be great to properly model the column, and so the number of equations can be reduced. Also, since a system is seldom mass-transfer-limited in the vapor phase, the rate equations for the vapor can be eliminated. To force a solution, a combination of this technique with a homotopy method may be required. 

A closer look at the system, however, does pique curiosity. The initial pH within the chamber is not 7 but 2-3, and the reactions are non-equilibrium, often irreversible, and involve other intermediates that can become important end products. The acidic pH represents a problem in that thiolates, not thiols, are the operative reductants, thus cannot reduce at pH values below their typical i.e. 8-9. This is resolved by proteins, including mfp-6, by sequence specific effects such as flanking cationic groups that reduce the Cys pK, e.g. redox active Cys-59 in DsB-A has a p Tg of 3.5. Several Cys residues in mfp-6 are acidic, but specific p Tg values have yet to be measured. The non-equilibrium, irreversible nature of the oxidation reactions is a particular problem with Dopa and other catechols. Indeed, the chemical fate of catechols in mussel byssus is highly dependent on their location. In the cuticle, the fate of Dopa appears to be tris catecholato-Fe complexes in the thread and plaque core, Dopa forms covalent cross-links after oxidation to quinones, whereas at the plaque-substratum interface, it is some combination of metal chelates and reduced H-bonded Dopa on metal oxide surfaces. The reducing capacity of mfp-6 plays a role in maximizing the latter and is astonishingly sustained, i.e. >21 days. ... 

As a result of AEP, the initial system of the set of Eqs. (81) is reduced to the equation describing the diffusional motion of a Brownian particle which undergoes the action of an additive and a multiplicative noise (with intensities D and Q, respectively) in the presence of a renormalized bounding potential, Eq. (90). The Markovian l t corresponds to X 00. If we take such a limit at a ed value of y, = d, and the case studied by Htoggi is recovered. Of course, having neglected the condition X c y we have reduced the problem to a trivial diffusional (lowest-order) approximation. 

In general, in the above considerations the coordinate x is presumed to describe nuclear motion normal to the intersection line L of the diabatic.potential energy surfaces of reactants and products. In particular cases, however, the coordinate x can coincide with a dynamically separable reaction coordinate. Then, the whole manydimensional problem of calculating the transition probability for any energy value is simply reduced to a one-dimensional one. Such is, for instance, the situation in a system of oscillators making harmonic vibrations with the same frequency in both the initial and final state /67/ which we considered in Sec.3.1.1. The diabatic surfaces (50.1) then represent two similar (N+1>dimensional rotational paraboloids which intersect in a N-dimensional plane S, and the intersection... 

---