Ionic equilibrium calculations

The molten system KC1 - K2TaF7 was analyzed using the same method and the ionic equilibrium calculated is as follows ... 

Many commercially available and investigational drugs are anionic or cationic salt forms of weak acids or weak bases (undissociated). Their properties (solubility, partition coefficient, bioavailability, etc.) are strongly dependent upon the degree of ionization, the pH of the solution, and other constituents in the solutions of the drugs. In this chapter, ionic equilibrium calculations will be demonstrated in order to facilitate study of their properties. 

This shows that the ionic equilibrium calculation for the salt formed between the weak base and the strong acid is identical for the undissociated acid. Therefore, one may use Equation (2.15), Equation (2.17), and Equation (2.19) for the calculation of H+ in the solution of a salt between a weak base and a strong acid. Equation (2.21a), Equation (2.21b), and Equation (2.21c) may be used for the calculation of OH- in the solution of the salt between a weak acid and a strong base along with Equation (2.20). 

Advances continue in the treatment of detonation mixtures that include explicit polar and ionic contributions. The new formalism places on a solid footing the modeling of polar species, opens the possibility of realistic multiple fluid phase chemical equilibrium calculations (polar—nonpolar phase segregation), extends the validity domain of the EXP6 library,40 and opens the possibility of applications in a wider regime of pressures and temperatures. 

Predictions of high explosive detonation based on the new approach yield excellent results. A similar theory for ionic species model43 compares very well with MD simulations. Nevertheless, high explosive chemical equilibrium calculations that include ionization are beyond the current abilities of the Cheetah code, because of the presence of multiple minima in the free energy surface. Such calculations will require additional algorithmic developments. In addition, the possibility of partial ionization, suggested by first principles simulations of water discussed below, also needs to be added to the Cheetah code framework. 

Equations 8 and 9 can be used for values of I up to 1. M. The second term in these equations accounts for the reversal of slope of activity coefficient versus ionic strength from negative to positive as ionic strength increases. Equations 8 and 9 have been widely used in the equilibrium calculations of the lime or limestone processes. With coals of moderate chloride content and for systems without extensive sludge dewatering, the ionic strength is well below 1.0 M, and equations 8 and 9 reasonable. 

Although one can probably find exceptions, most equilibrium calculations involving flue gas slurries are performed with temperature as a known variable. With temperature known, the numerical values of the appropriate equilibrium constants can be immediately calculated. The remaining unknown variables to be determined are the activities, activity coefficients, molalities, and the gas phase partial pressures. The equations used to determine these variables are formulated from among the equilibrium expressions presented in Table 1, the expressions for the activity coefficients, ionic strength, material balance expressions, and the electroneutrality balance. Although there are occasionally exceptions, the solution sequence generally is an iterative or cyclic sequence. 

Below we present a well-known calculation of membrane potential based on the classical Teorell-Meyer-Sievers (TMS) membrane model [2], [3]. The essence of this model is in treating the ion-selective membrane as a homogeneous layer of electrolyte solution with constant fixed charge density and with local ionic equilibrium at the membrane/solution interfaces. In spite of the obvious idealization involved in the first assumption the TMS model often yields useful results and represents in fact the main tool for practical membrane calculations. We shall return to TMS once again in 4.4 when discussing the electric current effects upon membrane selectivity. In the case of our present interest, the simplest TMS model of membrane potential for a 1,2 valent electrolyte reads... 

For mote about equilibrium calculations, see W. B. Guenther, Unified Equilibrium Calculations (New York Wiley, 1991) J. N. Butler. Ionic Equilibrium Solubility and pH Calculations (New York Wiley, 1998) and M. Meloun, Computation of Solution Equilibria (New York Wiley, 1988). For equilibrium calculation software, see http //www.micromath.com/ and http //www.acadsoft.co.uk/... 

From the various possible closures, the mean spherical approximation (MSA) [189] has found particularly wide attention in phase equilibrium calculations of ionic fluids. The Percus-Yevick (PY) closure is unsatisfactory for long-range potentials [173, 187, 190]. The hypemetted chain approximation (HNC), widely used in electrolyte thermodynamics [168, 173], leads to an increasing instability of the numerical algorithm as the phase boundary is approached [191]. There seems to be no decisive relation between the location of this numerical instability and phase transition lines [192-194]. Attempts were made to extrapolate phase transition lines from results far away, where the HNC is soluble [81, 194]. 

All reactions of this type must be studied at constant ionic strength because of the considerable non-ideality of the solutions, and the equilibrium calculations should be given in terms of activities rather than concentrations. 

The goal of this research was to improve activity coefficient prediction, and hence, equilibrium calculations in flue gas desulfurization (FGD) processes of both low and high ionic strength. A data base and methods were developed to use the local composition model by Chen et al. (MIT/Aspen Technology). The model was used to predict solubilities in various multicomponent systems for gypsum, magnesium sulfite, calcium sulfite, calcium carbonate, and magnesium carbonate SCU vapor pressure over sulfite/ bisulfite solutions and, C02 vapor pressure over car-bonate/bicarbonate solutions. 

Because K, depends on concentrations and the product KyKx is concentration independent, Kx must also depend on concentration. This shows that the simple equilibrium calculations usually carried out in first courses in chemistry are approximations. Actually such calculations are often rather poor approximations when applied to solutions of ionic species, where deviations from ideality are quite large. We shall see that calculations using Eq. (47) can present some computational difficulties. Concentrations are needed in order to obtain activity coefficients, but activity coefficients are needed before an equilibrium constant for calculating concentrations can be obtained. Such problems are usually handled by the method of successive approximations, whereby concentrations are initially calculated assuming ideal behavior and these concentrations are used for a first estimate of activity coefficients, which are then used for a better estimate of concentrations, and so forth. A G is calculated with the standard state used to define the activity. If molality-based activity coefficients are used, the relevant equation is... 

From Eqn. (14) it follows that with an exothermic reaction - and this is the case for most reactions in reactive absorption processes - decreases with increasing temperature. The electrolyte solution chemistry involves a variety of chemical reactions in the liquid phase, for example, complete dissociation of strong electrolytes, partial dissociation of weak electrolytes, reactions among ionic species, and complex ion formation. These reactions occur very rapidly, and hence, chemical equilibrium conditions are often assumed. Therefore, for electrolyte systems, chemical equilibrium calculations are of special importance. Concentration or activity-based reaction equilibrium constants as functions of temperature can be found in the literature [50]. 

J. N. Butler, Ionic Equilibrium Solubility and pH Calculations, Wiley Interscience, New York, 1998. 

The following example shows how this can be modeled in PHREEQC. First of all, a master- ami a solution species tritium T or T1 have to be defined. Since the input of data for log k und -gamma within the key word SOLUTION SPECIES is required, but unknown, any value can be entered here as a free parameter ( dummy , e g. 0.0). This value is not used for kinetic calculations and thus, does not cause any problems. However, all results based on equilibrium calculations (e.g. the calculation of the saturation index) are nonsense for this species . The tritium values have to be entered in tritium units. However, in order not to have to define or convert them in an extra step, they are entered fictitiously with the unit umol/kgw instead of TU in PHREEQC. As no interactions of tritium with any other species are defined, the unit is eventually irrelevant. After modeling, remember that the result is displayed in mol/kgw as always in PHREEQC and has to be recalculated to the fictitious tritium unit umol/kgw. Entering mol/kgw in the input file, the solution algorithm quits due to problems with too high total ionic strengths. 

The calculation of the change in HC03 concentration with pH, done according to the ionic equilibrium in the solution indicate that lgC pH) is linear for pH = 1 - 6. Probably... 

This chapter is concerned primarily with the computation of potentials of a cell using the hydrogen electrode as a probe for studying ionic equilibrium processes in mixed-organic-aqueous solvent systems. Computation of a number of other thermodynamic functions of the ionic process under investigation or of the solvent used is rather straightforward once the standard potential of the measuring cell has been calculated. 

---