Reaction path construction

There is another useiiil way of depicting the ideas embodied in the variable transition state theory of elimination reactions. This is to construct a three-dimensional potential energy diagram. Suppose that we consider the case of an ethyl halide. The two stepwise reaction paths both require the formation of high-energy intermediates. The El mechanism requires formation of a carbocation whereas the Elcb mechanism proceeds via a caibanion intermediate. 

In this chapter we consider how to construct reactions paths that account for the effects of simple reactants, a name given to reactants that are added to or removed from a system at constant rates. We take on other types of mass transfer in later chapters. Chapter 14 treats the mass transfer implicit in setting a species activity or gas fugacity over a reaction path. In Chapter 16 we develop reaction models in which the rates of mineral precipitation and dissolution are governed by kinetic rate laws. 

As an example, we consider the reaction path traced in the previous section (Fig. 13.1). To extract the overall reaction for each segment of the path, we construct a plot as described above. The result is shown in Figure 13.2. There are three segments in the reaction path the precipitation of kaolinite, the transformation of kaolinite to muscovite, and the continued formation of muscovite once the kaolinite is exhausted. There is a distinct overall reaction for each segment. 

In this chapter we consider how to construct reaction models that are somewhat more sophisticated than those discussed in the previous chapter, including reaction paths over which temperature varies and those in which species activities and gas fugacities are buffered. The latter cases involve the transfer of mass between the equilibrium system and an external buffer. Mass transfer in these cases occurs at rates implicit in solving the governing equations, rather than at rates set explicitly by the modeler. In Chapter 16 we consider the use of kinetic rate laws, a final method for defining mass transfer in reaction models. 

To model the chemical effects of evaporation, we construct a reaction path in which H2O is removed from a solution, thereby progressively concentrating the solutes. We also must account in the model for the exchange of gases such as CO2 and O2 between fluid and atmosphere. In this chapter we construct simulations of this sort, modeling the chemical evolution of water from saline alkaline lakes and the reactions that occur as seawater evaporates to desiccation. 

In this chapter we construct a variety of kinetic reaction paths to explore how this class of model behaves. Our calculations in each case are based on kinetic rate laws determined by laboratory experiment. In considering the calculation results, therefore, it is important to keep in mind the uncertainties entailed in applying laboratory measurements to model reaction processes in nature, as discussed in detail in Section 16.2. 

In Section 4.3.3, it was explained how to construct the reaction (diffusion) path for ternary and higher solid solution systems. In practice, one plots, for example, in a ternary system, the composition variables (measured along the pertinent space coordinate of the reacting solid) into a Gibbs phase triangle, noting that the spatial information is thereby lost. For certain boundary conditions, such a reaction path is independent of reaction time and therefore characterizes the diffusion process. For a one dimensional ternary system with stable interfaces, these boundary conditions are c,-( = oo,f) = c°( oo) q( <0,0) = c (-oo) c,(f>0,0) = c (+oo). 

In the context of the morphological evolution of non-equilibrium systems, let us then ask whether the reaction path, when constructed for a system with stable interfaces, can tell us something about the instability of moving boundaries. For this we... 

A number of empirical tunneling paths have been proposed in order to simplify the two-dimensional problem. Among those are the MEP [Kato et al., 1977], the sudden straight line [Makri and Miller, 1989b], and the so-called expectation value path [Shida et al., 1989]. The results of these papers are hard to compare because somewhat different PES s were used. As to the expectation value path, it was constructed as a parametric line q(Q) on which the vibrational coordinate q takes its expectation value when Q is fixed. Clearly, for the PES at hand this path coincides with MEP, since q is the coordinate of a harmonic oscillator. The results for the tunneling splitting calculated with the use of some of the earlier proposed reaction paths for PES (4.41) are compiled by Bosch et al. 

The methods of constructing different reaction paths are described in numerous papers and reviews (see, for example, Truhlar and Garrett [1984, 1987], Garrett et al. [1988], Ischtwan and Collins [1988], and references therein). In the IRC method proposed by Fukui [1970], the steepest descent path from the saddle point of a multidimensional PES V(X) to the reactant and product valleys is found by numerically solving the equation... 

Three approaches leading to 8 were considered (Figure 3.6.9) Following the biosynthetic pathway directly, polymer-supported thiamine 9 was constructed (path A) and could lead via crossed acyloin couplings to the target structure. Polymer-supported hydrazones 10 were reported to add directly to aldehydes in a non-catalyzed Umpolung reaction (path B) with results reported in due course. Finally, phosphine ylides 11 were investigated as polymer-supported acyl anion equivalents (path C). 

The isolation of these closely related thiolate complexes hints at an important role for 172-vinyl ligands in reactions which lead to net ligand substitution at metal. The SR bridge between Cp and W may resemble a snapshot along a reaction path for alkyne insertion into a M—L bond or for transfer of L from an T 2-vinyl to metal (97). A mechanism for alkyne polymerization based on rj2-vinyl intermediates has also been constructed (186). 

The above calculations provided the electronic ground and the first nine excited energies as well as the corresponding (transition) dipoles, at each point of the above reaction path. Such unperturbed Hamiltonian eigenstates defined the basis set used to construct the perturbed Hamiltonian matrix, Eq. 8-1, which was then diagonalized at each simulation frame, leading to the reaction free energy and related properties. 

---