Equilibrium, chemical Estimation methods

Two standard estimation methods for heat of reaction and CART are Chetah 7.2 and NASA CET 89. Chetah Version 7.2 is a computer program capable of predicting both thermochemical properties and certain reactive chemical hazards of pure chemicals, mixtures or reactions. Available from ASTM, Chetah 7.2 uses Benson s method of group additivity to estimate ideal gas heat of formation and heat of decomposition. NASA CET 89 is a computer program that calculates the adiabatic decomposition temperature (maximum attainable temperature in a chemical system) and the equilibrium decomposition products formed at that temperature. It is capable of calculating CART values for any combination of materials, including reactants, products, solvents, etc. Melhem and Shanley (1997) describe the use of CART values in thermal hazard analysis. 

There are numerous sources of phase equilibrium data available that serve as a database to those developing or improving equations of state. References to these databases are widely available. In addition, new data are added mainly through the Journal of Chemical and Engineering Data and the journal Fluid Phase Equilibria. Next, we give data for two systems so that the reader may practice the estimation methods discussed in this chapter. 

F. Van Zeggeten and S. H. Storey, The Computation of Chemical Equilibria, Cambridge University Press, Cambridge, UK, 1970 W. R. Smith and R. W. Missen, Chemical Reaction Equilibrium Analysis Theory and Algorithms, Wiley-Interscience, New York, 1982 W. J. Lyman, W. F. Reehl, and D. H. Rosenblatt, eds.. Handbook of Chemical Property Estimation Methods, McGraw-HiU, New York, 1982 C. M. Wal and S. G. Hutchison, J. Chem. Educ. 66, 546 (1989) F. G. Heherich, Chemical Engineering Education, 1989. 

Proper answers are rather complex, because different properties and conditions of a chemical system affect both equilibrium and reaction rate. Although the questions are related, no unified quantitative treatment yet exists, and to a large extent they are handled separately by the sciences of thermodynamics and reaction kinetics. Fortunately, with the help of thermodynamic and kinetics, the questions can be answered for many reactions with the aid of data and generalizations obtained by thermal, spectroscopic, and chromatographic measurements, and/or experimental computer chemistry, and the estimation methods of Benson [15]. 

Furthermore, this book is about equilibrium, mostly for mixtures. The steam tables and other such tables are for pure substances. Few if any such detailed tables exist for mixtures. We will see that much of the rest of this book is devoted to ways of measuring and/or estimating the thermodynamic properties of mixtures. Once we know those properties, we can carry out the material and energy balances on which to base our chemical plant designs. These estimating methods generally begin with the assumption that we have or can estimate the properties of the individual pure substances that make up the mixtures. The next few sections describe where pure substance tables like the steam tables come from, and how we would make up such a table (or the part we need of it) for some new substance. 

The absence of chemical equilibrium between media phases is conveniently expressed by unequal fugacities. This condition dramatically influences the concentration gradients in the vicinity of interfaces and the chemical transport rate or flux. This forms the thermodynamic pillar, commonly termed chemical equilibrium partitioning between the environmental media. A huge literature exists on partition coefficients, and estimation methods are well developed at least for conventional contaminants. 

These models are semiempirical and are based on the concept that intermolecular forces will cause nonrandom arrangement of molecules in the mixture. The models account for the arrangement of molecules of different sizes and the preferred orientation of molecules. In each case, the models are fitted to experimental binary vapor-liquid equilibrium data. This gives binary interaction parameters that can be used to predict multicomponent vapor-liquid equilibrium. In the case of the UNIQUAC equation, if experimentally determined vapor-liquid equilibrium data are not available, the Universal Quasi-chemical Functional Group Activity Coefficients (UNIFAC) method can be used to estimate UNIQUAC parameters from the molecular structures of the components in the mixture3. 

As we will see further in the book, almost all methods for calculating free energies in chemical and biological problems by means of computer simulations of equilibrium systems rely on one of the three approaches that we have just outlined, or on their possible combination. These methods can be applied not only in the context of the canonical ensemble, but also in other ensembles. As will be discussed in Chap. 5, AA can be also estimated from nonequilibrium simulations, to such extent that FEP and TI methods can be considered as limiting cases of this approach. 

Very few generalized computer-based techniques for calculating chemical equilibria in electrolyte systems have been reported. Crerar (47) describes a method for calculating multicomponent equilibria based on equilibrium constants and activity coefficients estimated from the Debye Huckel equation. It is not clear, however, if this technique has beep applied in general to the solubility of minerals and solids. A second generalized approach has been developed by OIL Systems, Inc. (48). It also operates on specified equilibrium constants and incorporates activity coefficient corrections for ions, non-electrolytes and water. This technique has been applied to a variety of electrolyte equilibrium problems including vapor-liquid equilibria and solubility of solids. 

In the second approach, the chemical equilibrium between the reactant(s) and the transition state is expressed in terms of conventional thermodynamic functions, i.e., enthalpy and entropy changes. This method is easier to implement and provides useful insights for estimating both the preexponential factors and the activation energies. Consequently, we shall utilize the thermodynamic formulation of the TST in this paper. 

In the application of Anet s equations (see Section II,B,5 below) to the estimation of AG and AG° by low temperature 13C-NMR spectroscopy, the magnitude of Av (the chemical shift difference between the exchanging sites) is required. Because this cannot always be observed, resort has to be made to some indirect method of estimation of Av. This has been done, for example, in the case of the l,2,4-trimethylhexahydro-l,2,4-triazine equilibrium (Section III,F,3) by estimating the chemical shifts in the various conformers from chemical shift effects based on model systems (Table VIII). Utilization of... 

With the discussion of the free-energy function G in this chapter, all of the thermodynamic functions needed for chemical equilibrium and kinetic calculations have been introduced. Chapter 8 discussed methods for estimating the internal energy E, entropy S, heat capacity Cv, and enthalpy H. These techniques are very useful when the needed information is not available from experiment. 

Determination of T y. In the formulation of the phase equilibrium problem presented earlier, component chemical potentials were separated into three terms (1) 0, which expresses the primary temperature dependence, (2) solution mole fractions, which represent the primary composition dependence (ideal entropic contribution), and (3) 1, which accounts for relative mixture nonidealities. Because little data about the experimental properties of solutions exist, Tg is usually evaluated by imposing a model to describe the behavior of the liquid and solid mixtures and estimating model parameters by semiempirical methods or fitting limited segments of the phase diagram. Various solution models used to describe the liquid and solid mixtures are discussed in the following sections, and the behavior of T % is presented. 

Usually, this method is applied to enzymatic reactions, and the equilibrium IEs are obtained along with kinetic IEs that are of greater interest. An example is the deuterium IE on the reaction of acetone-c/6 with NADH, to form 2-propanol-fi 6 + NAD+. A mixture of acetone-c/6 and 2-propanol is prepared along with coreactants NADH and NAD+ at concentrations such that the reaction is at chemical equilibrium. Isotopic equilibration is initiated by adding enzyme. In this case the spectral signature lies in the NADH, but the measured maximum or minimum of absorbance provides the right-hand side of Equation (25) or (26) and thus a for each mixture. An estimate of AThh is needed to solve for each R in Equation (23) in order to fit the data to Equation (27), but after successive iterations the values of R and XEIE converge. 

---