Distribution of species

Equations (2.14) to (2.23) were obtained on the assumption that electroactive molecules are uniformly distributed in the entire volume of the solid. However, this is an unrealistic assumption in cases where bulky guest species are entrapped within the cavities of the porous material. Here, ship-in-a-bottle synthetic procedures most likely yield a nonuniform distribution of guest species in the network of the host porous material. Electrochemical data can then be used for obtaining information on the distribution of electroactive species. 

These concentration profiles, however, must be taken with caution. First of all, it should be noted that the obtained concentration profiles represent the distribution of effectively electroactive molecules within the solid but not necessarily the entire distribution of such molecules. As far as electron hopping and/or ion hopping can 

FIGURE 2.9 Chronocoulometric curve recorded under application of a potential of -0.35 V for PY+ Y deposited on paraffin-impregnated graphite electrode in contact with 0.10 M Et4NClO4/MeCN. 

---

Diflfiisive processes nonnally operate in chemical systems so as to disperse concentration gradients. In a paper in 1952, the mathematician Alan Turing produced a remarkable prediction [37] that if selective diffiision were coupled with chemical feedback, the opposite situation may arise, with a spontaneous development of sustained spatial distributions of species concentrations from initially unifonn systems. Turmg s paper was set in the context of the development of fonn (morphogenesis) in embryos, and has been adopted in some studies of animal coat markings. With the subsequent theoretical work at Brussels [1], it became clear that oscillatory chemical systems should provide a fertile ground for the search for experimental examples of these Turing patterns. 

The distribution of species as a function of anolyte pH is shown in Eigure 8 (21). Chlorine is the primary product at pH < 4. Hypochlorous acid [7790-92-3] HOCl, is the dorninant species from pH 4 to pH 6. Above pH 6 combinations of HOCl and hypochlorite ion, C10 , exist but react to produce CIO- (22) ... 

Let us consider a simple model of a quenched-annealed system which consists of particles belonging to two species species 0 is quenched (matrix) and species 1 is annealed, i.e., the particles are allowed to equlibrate between themselves in the presence of 0 particles. We assume that the subsystem composed of 0 particles has been a usual fluid before quenching. One can characterize it either by the density or by the value of the chemical potential The interparticle interaction Woo(r) does not need to be specified for the moment. It is just assumed that the fluid with interaction woo(r) has reached an equlibrium at certain temperature Tq, and then the fluid has been quenched at this temperature without structural relaxation. Thus, the distribution of species 0 is any one from a set of equihbrium configurations corresponding to canonical or grand canonical ensemble. We denote the interactions between annealed particles by Un r), and the cross fluid-matrix interactions by Wio(r). 

The role of the distribution of species in solution in determining the CdS film composition and structure was studied by Rieke and Bentjen [244], who performed equilibrium analysis of the cadmium-amine-hydroxide system to predict the spe-ciation in solution. The focus was on the formation of Cd(OH)2 and Cd(NH3) species due to their importance in film growfh. If was concluded fhat for deposition of high-quality, adherent, phase-pure CdS films, a surface cafalytically active toward thiourea decomposition is desirable. The Cd(OH)2 film was thought to be responsible for this effect. 

Models describing the adsorption of water-miscible organic compounds on natural materials have not been correlated with field observations under typical injection-zone conditions. Few computer codes contain algorithms for calculating the distribution of species between the adsorbed and aqueous states. 

The clay ion-exchange model assumes that the interactions of the various cations in any one clay type can be generalized and that the amount of exchange will be determined by the empirically determined cation-exchange capacity (CEC) of the clays in the injection zone. The aqueous-phase activity coefficients of the cations can be determined from a distribution-of-species code. The clay-phase activity coefficients are derived by assuming that the clay phase behaves as a regular solution and by applying conventional solution theory to the experimental equilibrium data in the literature.1 2 3... 

Two-step models, which first solve mass momentum and energy balances for each time step and then reequilibrate the chemistry using a distribution-of-species code. 

Fig. 2.5. Measurement of pKas of serotonin by target factor analysis (TFA). (A) 3-D spectrum produced by serotonin in pH gradient experiment (equivalent to A matrix). (B) Molar absorptivity of three serotonin species (equivalent to E matrix). (C) Distribution of species (equivalent to C matrix). In this graph the three sets of data points denote the three...

Garrels and Thompson s calculation, computed by hand, is the basis for a class of geochemical models that predict species distributions, mineral saturation states, and gas fugacities from chemical analyses. This class of models stems from the distinction between a chemical analysis, which reflects a solution s bulk composition, and the actual distribution of species in a solution. Such equilibrium models have become widely applied, thanks in part to the dissemination of reliable computer programs such as SOLMNEQ (Kharaka and Barnes, 1973) and WATEQ (Truesdell and Jones, 1974). 

Once we have calculated the distribution of species in the fluid, we can determine the degree to which it is undersaturated or supersaturated with respect to the many minerals in the thermodynamic database. Only a few of the minerals can exist in equilibrium with the fluid, which is therefore undersaturated or supersaturated with respect to each of the rest. For any mineral A , we can write a reaction,... 

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,... 

Here, we set oxidation state in the model using the dissolved oxygen content and calculate the distribution of species assuming redox equilibrium. 

Geochemical modelers currently employ two types of methods to estimate activity coefficients (Plummer, 1992 Wolery, 1992b). The first type consists of applying variants of the Debye-Hiickel equation, a simple relationship that treats a species activity coefficient as a function of the species size and the solution s ionic strength. Methods of this type take into account the distribution of species in solution and are easy to use, but can be applied with accuracy to modeling only relatively dilute fluids. 

Virial methods, the second type, employ coefficients that account for interactions among the individual components (rather than species) in solution. The virial methods are less general, rather complicated to apply, require considerable amounts of data, and allow little insight into the distribution of species in solution. They can, however, reliably predict mineral solubilities even in concentrated brines. 

The virial methods differ conceptually from other techniques in that they take little or no explicit account of the distribution of species in solution. In their simplest form, the equations recognize only free ions, as though each salt has fully dissociated in solution. The molality m/ of the Na+ ion, then, is taken to be the analytical concentration of sodium. All of the calcium in solution is represented by Ca++, the chlorine by Cl-, the sulfate by SO4-, and so on. In many chemical systems, however, it is desirable to include some complex species in the virial formulation. Species that protonate and deprotonate with pH, such as those in the series COg -HCOJ-C02(aq) and A1+++-A10H++-A1(0H), typically need to be included, and incorporating strong ion pairs such as CaSO aq) may improve the model s accuracy at high temperatures. Weare (1987, pp. 148-153) discusses the criteria for selecting complex species to include in a virial formulation. 

In the virial methods, therefore, the activity coefficients account implicitly for the reduction in the free ion s activity due to the formation of whatever ion pairs and complex species are not included in the formulation. As such, they describe not only the factors traditionally accounted for by activity coefficient models, such as the effects of electrostatic interaction and ion hydration, but also the distribution of species in solution. There is no provision in the method for separating the traditional part of the coefficients from the portion attributable to speciation. For this reason, the coefficients differ (even in the absence of error) in meaning and value from activity coefficients given by other methods. It might be more accurate and less confusing to refer to the virial methods as activity models rather than as activity coefficient models. 

Unlike the Debye-Hiickel equations, the virial methods provide little or no information about the distribution of species in solution. Geochemists like to identify the dominant species in solution in order to write the reactions that control a system s behavior. In the virial methods, this information is hidden within the complexities of the virial equations and coefficients. Many geochemists, therefore, find the virial methods to be less satisfying than methods that predict the species distribution. The information given by Debye-Hiickel methods about species distributions in concentrated solutions, however, is not necessarily reliable and should be used with caution. 

In this case (Fig. 14.10), we observe a more complicated distribution of species. At low pH, the H+ activity and positive surface potential drive SOJ to sorb according to the reaction,... 

Alternative ways to describe the distribution of species in the mixture are by mass fraction,... 

Pyridinecarboxaldehyde, 3. Possible hydration of the aldehyde group makes the aqueous solution chemistry of 3 potentially more complex and interesting than the other compounds. Hydration is less extensive with 3 than 4-pyridinecarboxaldehyde but upon protonation, about 80% will exist as the hydrate (gem-diol). The calculated distribution of species as a function of pH is given in Figure 4 based on the equilibrium constants determined by Laviron (9). 

Fig. 17. Distribution of species as a function of pH, calculated from constants in Table IX. Tungsten(VI) concentration 0.001 M (A) and 0.20 M (B).

The purpose of this chapter is to outline the simplest methods of arriving at a description of the distribution of species in mixtures of liquids, gases and solids. Homogeneous equilibrium deals with single phase systems, such as electrolyte solutions (e.g., seawater) or gas mixtures (e.g., a volcanic gas). Heterogeneous equilibrium involves coexisting gaseous, liquid and solid phases. 

Figure 8.11 Evidence of cross-diffusional effects. The homogeneous distribution of species 2 (dashed line, top) is perturbed by a coexisting gradient of species 1 (bottom).

More sophisticated calculations (14,20), using either stochastic Monte Carlo or deterministic methods, are able to consider not only different Irradiating particles but also reactant diffusion and variations In the concentration of dissolved solutes, giving the evolution of both transient and stable products as a function of time. The distribution of species within the tracks necessitates the use of nonhomogeneous kinetics (21,22) or of time dependent kinetics (23). The results agree quite well with experimental data. 

Although it is difficult to determine the spatial distribution of species experimentally, it provides an illustrative view of the electrode reaction. Simulations usually provide values of c = /(x, f) for each species as the primary result. The space dependence of c is termed a concentration... 

---