Fractionation in reacting systems

Integrating isotope fractionation into the reaction path calculation is a matter of applying the mass balance equations while tracing over the course of the reaction path the system s total isotopic composition. Much of the effort in programming an isotope model consists of devising a careful accounting of the mass of each isotope. 

The calculation begins with an initial fluid of specified isotopic composition. The model, by mass balance, assigns the initial compositions of the solvent and each aqueous species. The model then sets the composition of each unsegregated mineral to be in equilibrium with the fluid. The modeler specifies the composition of each segregated mineral as well as that of each reactant to be added to the system. 

The model traces the reaction path by taking a series of steps along reaction progress, moving forward each step from i to 2- Over a step, the system s isotopic composition can change in two ways reactants can be added or removed, and segregated minerals can dissolve. 

The parallel equations for hydrogen, carbon, and sulfur isotopes are, 

Once we have computed the total isotopic compositions, we calculate the compositions of the reference species using the mass balance equations (Eqns. 19.13, 19.20, 19.21, 19.22). We can then use the isotopic compositions of the reference species to calculate the compositions of the other species (Eqns. 19.10, 19.14, 19.15, 19.16) and the unsegregated minerals (Eqns. 19.11,19.17,19.18, 19.19). 

Using 18O as an example, we can calculate the composition 818OT( ,2) at the end of the step according to the equation 

---

In this chapter, we develop a mass balance model of the fractionation in reacting systems of the stable isotopes of hydrogen, carbon, oxygen, and sulfur. We then demonstrate application of the model by simulating the isotopic effects of the dolomitization reaction of calcite. 

Lee, M.-K. and C. M. Bethke, 1996, A model of isotope fractionation in reacting geochemical systems. American Journal of Science 296,965-988. 

Any property of a reacting system that changes regularly as the reaction proceeds can be formulated as a rate equation which should be convertible to the fundamental form in terms of concentration, Eq. (7-4). Examples are the rates of change of electrical conductivity, of pH, or of optical rotation. The most common other variables are partial pressure p and mole fraction Ni. The relations between these units... 

To evaluate the rate constant in a more rigorous way (as is done in detail in Ref. 1), let us consider first the average behavior of many reacting systems with a one-dimensional surface of the type described in Fig. 2.1. We will try to determine what fraction of the systems that pass the reactant region toward the barrier would react. The number of systems in a path length Ax, which have a momentum between p and p + Ap, is given by (see Ref. lb)... 

The fraction conversion is based on the limiting reagent, in this case assumed to be species A. The parameter SA is the fraction change in the volume of a closed reacting system between zero conversion and complete conversion. As such it may take on both positive and negative values. Hence... 

For a catalyst-reactant system in which the effectiveness factors for the first and second reactions both approach unity, the analysis presented in Sections 9.2 and 9.3 is appropriate. If yields are based on the initial concentration of reactant A and if no V or W is present in the original feed, equations 9.3.9 and 9.3.10 can be rearranged to express the yield of V as a function of the fraction A reacted ... 

All the described properties of such a s-fraction of poly(NVCl-co-NVIAz) synthesized at the temperature above the PST of the reacting system allowed us to draw the conclusion that the chains of this type had the comonomer sequence, which at the temperatures above the conformation transition facilitated the formation of polymer particles, where H-blocks are in the interior shielded by the P-blocks against additional intermolecular association. Such a behaviour of this copolymer in aqueous media is close to that of oligomeric proteins similar to casein [46] possessing a rather hydrophobic core surrounded by the polar segments. 

The process design of a batch reactor may involve determining the time (t) required to achieve a specified fractional conversion (/A, for limiting reactant A, say) in a single batch, or the volume (V) of reacting system required to achieve a specified rate of production on a continual basis. The phrase continual basis means an ongoing operation,... 

Likewise, the matrix M can be found using the method presented in Section 5.3 applied in the limiting case of a non-reacting system (i.e S = 0). In the simplest case (binary mixing), only one mixture-fraction component is required, and Mj is easily found from the species concentrations in the inlet streams. 

The simplest are statistical theories, where the input information is reduced to the distribution of units in different reaction states. The reaction state of a unit is defined by the number and type of bonds issuing from the unit. In a reacting system, the distribution fraction of units in different reaction states is a function of the reaction time (conversion) (cf. e.g. [7, 8, 29, 30] and can be obtained either experimentally (e.g. by NMR) or calculated by solution of a few simple kinetic differential equations. An example of reaction state distribution of an AB2 unit is... 

If kinetic data are to be used, it is necessary to transform the variables to conform with those of the partial equilibrium model. The units used in the model equations for and nj are moles formed/kg of solution. Thus the mass of solution in the reacting system from which the kinetic data comes must be known. Frequently, one will know the volume and have to approximate the density. A relation between and t is also needed. For this, the mass of solid originally present must be known. The amount of solid reacing, -ANg, for a time interval At can be obtained from rate curves or calculated from an integrated rate equation. The fraction of the original mass reacting in the time interval gives an approximate value of 5, e.g.,... 

In combustion problems, one is interested in the rate of energy conversion or utilization. Thus it is more convenient to deal with the fractional change of a particular substance rather than the absolute concentration. If (M) is used to denote the concentrations in a chemical reacting system of arbitrary order n, the rate expression is... 

A clear avenue of future research is to explore the S-Fe redox couple in biologic systems. Bacterial sulfate reduction and DIR may be spatially decoupled, dependent upon the distribution of poorly crystalline ferric hydroxides and sulfate (e.g., Canfield et al. 1993 Thamdrup and Canfield 1996), or may be closely associated in low-suUate environments. Production of FIjS from bacterial sulfate reduction may quickly react with Fefll) to form iron sulfides (e.g., Sorensen and Jeorgensen 1987 Thamdrup et al. 1994). In addition to these reactions, Fe(III) reduchon may be coupled to oxidation of reduced S (e.g., Thamdrup and Canfield 1996), where the net result is that S and Fe may be cycled extensively before they find themselves in the inventory of sedimentary rocks (e.g., Canfield et al. 1993). Investigation of both S and Fe isotope fractionations produced during biochemical cycling of these elements will be an important future avenue of research that will bear on our understanding of the isotopic variations of these elements in both modem and ancient environments. 

In a system containing equivalent numbers of A and B groups, the number of monomer molecules present initially is No and the corresponding total number of functional groups is Af0/avg. If IV is the number of molecules after reaction has occurred, then 2(A/ — N) is the number of functional groups that have reacted. The extent of reaction p is the fraction of functional groups lost ... 

This parameter can be termed as fractional mole change. Note that in the calculation of 8, inerts that are involved in the gas phase are not taken into account. They are taken into account only in the determination of the total moles of the reacting system. 

Using appropriate software, simulate the chemically reacting flow in this system. Plot the calculated Si mole fraction versus height above the susceptor. What is the maximum density, number of atoms/cm3 At what height does this occur What accounts for the sharp drop in Si number density at the surface ... 

A unimolecular reaction is one in which the absolute rate of change is proportional to the first power of the concentration of the reacting substance. The fraction of the total number of the molecules in the system which change in unit time is therefore independent of the concentration, and thus, in gaseous systems, cannot be proportional to the nuinber of collisions undergone in unit time by the molecules. It must therefore be concluded that, whether or not previously received collisions have done anything to put the molecule into an abnormal condition, the actual chemical transformation is an event happening to the isolated molecule. 

---