Chemical model determination

The purpose of this chapter is to discuss the task of chemical model determination (CMD) not only with respect to solvent extraction (SX) and liquid membranes (LM) but also to give overview of the fundamentals and also to elucidate the problems of extraction equilibria computation more deeply. There is no intention of providing an exhaustive description of computation in extraction systems due to the limited space and because it is not necessary. Also no exhaustive review of the computer programs will be given here but will be limited to the most important questions. 

FIGURE 3.1 Scheme of classic trial-and-error approach toward the chemical model determination using general regression minimization methods. 

A new scheme for the chemical model determination in equilibria studies might include except statistical tests (1) the application of factor analysis to estimate directly the number of species in solution and (2) the method of direct computation of species stoichiometry via simultaneous calculation of stability constants and stoichiometric indices—the ESI approach. The ESI method might also be used as a new diagnostic tool when searching the best chemical model. The proposed novel scheme of CMD based on PCA and ESI is schematically given in Figure 3.6. 

J. Havel, in Computer Applications in Chemistry VIII. Problem of Chemical Model Determination in Equilibrium Study, Scripta Fac. Sci. Univ. Purk. Brun. 17, 305-320 (1987). 

In the last two decades experimental evidence has been gathered showing that the intrinsic properties of the electrolytes determine both bulk properties of the solution and the reactivity of the solutes at the electrodes. Examples covering various aspects of this field are given in Ref. [16]. Intrinsic properties may be described with the help of local structures caused by ion-ion, ion-solvent, and solvent-solvent interactions. An efficient description of the properties of electrolyte solutions up to salt concentrations significantly larger than 1 mol kg 1 is based on the chemical model of electrolytes. 

The commonly used method for the determination of association constants is by conductivity measurements on symmetrical electrolytes at low salt concentrations. The evaluation may advantageously be based on the low-concentration chemical model (lcCM), which is a Hamiltonian model at the McMillan-Mayer level including short-range nonelectrostatic interactions of cations and anions [89]. It is a feature of the lcCM that the association constants do not depend on the physical... 

Atmospheric aerosols have a direct impact on earth s radiation balance, fog formation and cloud physics, and visibility degradation as well as human health effect[l]. Both natural and anthropogenic sources contribute to the formation of ambient aerosol, which are composed mostly of sulfates, nitrates and ammoniums in either pure or mixed forms[2]. These inorganic salt aerosols are hygroscopic by nature and exhibit the properties of deliquescence and efflorescence in humid air. That is, relative humidity(RH) history and chemical composition determine whether atmospheric aerosols are liquid or solid. Aerosol physical state affects climate and environmental phenomena such as radiative transfer, visibility, and heterogeneous chemistry. Here we present a mathematical model that considers the relative humidity history and chemical composition dependence of deliquescence and efflorescence for describing the dynamic and transport behavior of ambient aerosols[3]. 

Carell and Olin (58) were the first to derive thermodynamic functions relating to beryllium hydrolysis. They determined the enthalpy and entropy of formation of the species Be2(OH)3+ and Be3(OH)3+. Subsequently, Mesmer and Baes determined the enthalpies for these two species from the temperature variation of the respective equilibrium constants. They also determined a value for the species Be5(OH) + (66). Ishiguro and Ohtaki measured the enthalpies of formation of Be2(OH)3+ and Be3(OH)3+ calorimetrically in solution in water and water/dioxan mixtures (99). The agreement between the values is satisfactory considering the fact that they were obtained with different chemical models and ionic media. 

For a first chemical model, we calculate the distribution of species in surface seawater, a problem first undertaken by Garrels and Thompson (1962 see also Thompson, 1992). We base our calculation on the major element composition of seawater (Table 6.2), as determined by chemical analysis. To set pH, we assume equilibrium with CO2 in the atmosphere (Table 6.3). Since the program will determine the HCOJ and water activities, setting the CO2 fugacity (about equal to partial pressure) fixes pH according to the reaction,... 

Thirteen minerals appear supersaturated in the first block of results produced by the chemical model (Table 6.6). These results, therefore, represent an equilibrium achieved internally within the fluid but metastable with respect to mineral precipitation. It is quite common in modeling natural waters, especially when working at low temperature, to find one or more minerals listed as supersaturated. Unfortunately, the error sources in geochemical modeling are large enough that it can be difficult to determine whether or not a water is in fact supersaturated. 

D-Glucose 6-phosphate is converted enzymically into L-wyo-inositol 1-phosphate (20) in a process which requires NAD+. The base-catalysed cyclization of d-xylo-hexos-5-ulose 6-phosphate (21), followed by reduction with borohydride, leads to (20) and epi-inositol 3-phosphate (22) (Scheme 3).59 This has been put forward as a chemical model for the enzymic synthesis. The phosphorylation of inositols with polyphosphoric acid has been described80 and the p-KVs of inositol hexaphosphate have been determined by 31P n.m.r.61... 

For a liquor of known pH and magnesium and chloride concentrations, the degree of gypsum saturation can be determined by measurement of either the total dissolved calcium or the total dissolved SO2 (sulfite plus bisulfite). The chemical model has been used to obtain correlations for gypsum saturation, presented below. The correlations, Equations 7 and 9, are valid for a typical scrubbing temperature of 50 °C, and for the same ranges of pH, magnesium, and chloride as for Equations 1-4. 

Use of Equations 2, 3, 5, and 9 to determine gypsum saturation gives results that agree with the chemical model to within a standard error of estimate of 0.04 fraction gypsum saturation for saturations of 0.5-2.0. 

However, many successful chemical models exist that do not necessarily have obvious connections with quantum mechanics. Typically, these models were developed based on intuitive concepts, i.e., their forms were determined inductively. In principle, any successful model must ultimately find its basis in quantum mechanics, and indeed a posteriori derivations have illustrated this point in select instances, but often the form of a good model is more readily grasped when rationalized on the basis of intuitive chemical concepts rather than on the basis of quantum mechanics (the latter being desperately non-intuitive at first blush). 

We have already determined the following information about the behaviour of the pool chemical model with the exponential approximation. There is a unique stationary-state solution for ass, the concentration of the intermediate A, and 0SS, the temperature rise, for any given combination of the experimental conditions /r and k. If the dimensionless reaction rate constant k is larger than the value e-2, then the stationary state is always stable. If heat transfer is more efficient, so that k Hopf bifurcation points along the stationary-state locus as /r varies (Fig. 4.4). If these bifurcation points are /r and /z (with the stationary state... 

The behaviour we are expecting to emerge from this physico-chemical model is that of a steady wave of reaction moving from left to right in Fig. 11.2 into the region of unreacted A. By steady, we really mean that the wavefront should maintain its shape as it moves with a constant speed. It is this shape and speed which we seek to determine (and express in terms of the rate constant, diffusion coefficient, etc.). 

---