Activity coefficient relation using

This example is based on the model description of Sec. 3.3.4, and involves a multicomponent, semi-batch system, with both heating and boiling periods. The compositions and boiling point temperatures will change with time. The water phase will accumulate in the boiler. The system simulated is based on a mixture of n-octane and n-decane, which for simplicity will be assumed to be ideal but which has been simulated using detailed activity coefficient relations by Prenosil (1976). 

Can the species activity coefficients be calculated accurately An activity coefficient relates each dissolved species concentration to its activity. Most commonly, a modeler uses an extended form of the Debye-Hiickel equation to estimate values for the coefficients. Helgeson (1969) correlated the activity coefficients to this equation for dominantly NaCl solutions having concentrations up to 3 molal. The resulting equations are probably reliable for electrolyte solutions of general composition (i.e., those dominated by salts other than NaCl) where ionic strength is less than about 1 molal (Wolery, 1983 see Chapter 8). Calculated activity coefficients are less reliable in more concentrated solutions. As an alternative to the Debye-Hiickel method, the modeler can use virial equations (the Pitzer equations ) designed to predict activity coefficients for electrolyte brines. These equations have their own limitations, however, as discussed in Chapter 8. 

We can estimate the activity coefficients by using the excess Gibbs energy models. Based on the local composition concept, the Wilson, NRTL, and UNIQUAC models for excess Gibbs energy provide relations for activity coefficient... 

This leads to the relation that the distribution coefficient (partition coefficient) is only influenced by the activity coefficients for both phases and depends only on pressure, temperature and concentration. The relationship of the Gibbsche excess enthalpy to the activity coefficient was used for a variety of modem calculations (Wilson, UNI-QUAC equations) ... 

When activity coefficients are used, y is plotted versus (log Ci + log/ ) to obtain T v The area per molecule at the interface provides information on the degree of packing and the orientation of the adsorbed surfactant molecule when compared with the dimensions of the molecule as obtained by use of molecular models. From the surface excess concentration, the area per molecule at the interface a, in square angstroms is calculated from the relation... 

Note Equations like (2.6) and (2.8) can also be put in a form where molarity, molality, or some other concentrative unit is used rather than mole fraction. This means that /ie has another value, but—more important—it also affects the value of the (apparent) activity coefficient. For a very dilute solution, the differences tend to be negligible, but in other cases, the concentrative unit to which the activity coefficient relates should be stated. Naturally, the various kinds of concentration can be recalculated into each other see Appendix A.7. 

It is useful to be aware of these relations since both mole fraction and activity coefficients are used both to define solubility and to explore the relation of solubility to the physical properties of molecules. 

The flash point of a solvent mixture is not identical to that of its most flammable component. When solvents with widely differing hydrogen bond parameters are mixed (e.g., alcohol hydrocarbon), tbe flash point is significantly reduced. On the other hand, the flash point of a mixture of chemically related solvents lies between those of the individual components [14.87]-[14.89]. Methods have been developed for calculating the flash point of solvent mixtures and solutions, activity coefficients are used to account for nonideal behavior [14.90]. 

Similar connections can be made for properties of mixing, excess properties, and activity coefficients. The use of this approach on fluctuation properties is described by Perry and O Connell (1984). The relations for TCFI among pairs of species are extremely complex with the extents of reactions embedded in the TCFI. However, those for DCFI are much more direct and the extents of reaction are contained in the projectors. Defining the desired matrix. 

It what concerns binary systems involving non-electrolytes, we have been measuring mutual diffusion coefficients of some cyclodextrins (a-CD, P-CD, HP-a-CD, and HP-p-CD) [12, 13] and some dmgs (e.g., caffeine and isoniazid [14]) in aqueous solutions. Also, from comparison of these experimental diffusion coefficients with the related calculated values, it is possible to give some stmctural information, such as diffusion coefficients at infinitesimal concentration at different temperatures, estimation of activity coefficients by using equations of Hartley and Gordon, estimation of hydrodynamic radius, and estimation of activation energies, Ea, of the diffusion process at several temperatures. 

In such cases, only an iterative method can be used to calculate the activity coefficients because, in order to find the molar fractions, we must already know the activity coefficients to use relations such as [4.133],... 

Our laboratory in cooperation with several national and international academic and industrial partners is contributing to these efforts by the establishment of various dedicated characterization techniques (like activity coefficient measurements using GC technology) as well as determination of thermodynamic and physicochemical properties from a continuously growing portfolio of (functionalized) ionic liquids. Based on the received property data we published several papers related to the adjacent prediction of properties (like molar enthalpy of vaporization, parachor, interstice volume, interstice fractions, thermal expansion coefficient, standard entropy etc.). Additionally our laboratory created and launched a new most comprehensive Ionic Liquid property data base—delph-IL.(www.delphil.net). This fast growing collections of IL data will supvport researchers in the field to find and evaluate potential materials for their applications and hence decrease the time for new developments. 

It is strictly for convenience that certain conventions have been adopted in the choice of a standard-state fugacity. These conventions, in turn, result from two important considerations (a) the necessity for an unambiguous thermodynamic treatment of noncondensable components in liquid solutions, and (b) the relation between activity coefficients given by the Gibbs-Duhem equation. The first of these considerations leads to a normalization for activity coefficients for nonoondensable components which is different from that used for condensable components, and the second leads to the definition and use of adjusted or pressure-independent activity coefficients. These considerations and their consequences are discussed in the following paragraphs. 

If we vary the composition of a liquid mixture over all possible composition values at constant temperature, the equilibrium pressure does not remain constant. Therefore, if integrated forms of the Gibbs-Duhem equation [Equation (16)] are used to correlate isothermal activity coefficient data, it is necessary that all activity coefficients be evaluated at the same pressure. Unfortunately, however, experimentally obtained isothermal activity coefficients are not all at the same pressure and therefore they must be corrected from the experimental total pressure P to the same (arbitrary) reference pressure designated P. This may be done by the rigorous thermodynamic relation at constant temperature and composition ... 

Hydrocarbon mixtures can be assumed to be regular solutions it is thus possible to estimate the activity coefficient using a relation published by Hildebrandt (1950) ... 

Whereas the fundamental residual property relation derives its usefulness from its direct relation to experimental PVT data and equations of state, the excess property formulation is useful because and are all experimentally accessible. Activity coefficients are found from vapor—Hquid... 

Evaluation of 9 is usually by Eq. (4-196), based on the two-term virial equation of state, but other equations, such as Eq. (4-200), are also applicable. The activity coefficient Jj is evaluated by Eq. (4-119), which refates In Jj to G /RT as a partial proper. Thus, what is required for the hquid phase is a relation between G /BT and composition. Equations in common use for this purpose have already been described. 

Outlined below are the steps required for of a X T.E calciilation of vapor-phase composition and pressure, given the liquid-phase composition and temperature. A choice must be made of an equation of state. Only the Soave/Redlich/Kwong and Peng/Robinson equations, as represented by Eqs. (4-230) and (4-231), are considered here. These two equations usually give comparable results. A choice must also be made of a two-parameter correlating expression to represent the liquid-phase composition dependence of for each pq binaiy. The Wilson, NRTL (with a fixed), and UNIQUAC equations are of general applicabihty for binary systems, the Margules and van Laar equations may also be used. The equation selected depends on evidence of its suitability to the particular system treated. Reasonable estimates of the parameters in the equation must also be known at the temperature of interest. These parameters are directly related to infinite-dilution values of the activity coefficients for each pq binaiy. 

The proliferation of acidity functions is a consequence of the activity coefficient cancellation assumption. According to Eq. (8-89), a plot of log(cB/cBH+) against Hq should be linear with unit slope. Such plots are usually linear (for bases of closely related structure), but the slopes often differ from unity. - This behavior is an indication that the cancellation assumption (also called the zero-order approximation) is not valid, and several groups have devised alternatives. We will use the symbolism of Cox and Yates. ... 

For any component i in a liquid phase, the fugacity of f is most conveniently related to the mole fraction xt by use of the activity coefficient, y(, according to... 

At constant temperature, the activity coefficient depends on both pressure and composition. One of the important goals of thermodynamic analysis is to consider separately the effect of each independent variable on the liquid-phase fugacity it is therefore desirable to define and use constant-pressure activity coefficients which at constant temperature are independent of pressure and depend only on composition. The definition of such activity coefficients follows directly from either of the exact thermodynamic relations... 

The difficulties encountered in the Chao-Seader correlation can, at least in part, be overcome by the somewhat different formulation recently developed by Chueh (C2, C3). In Chueh s equations, the partial molar volumes in the liquid phase are functions of composition and temperature, as indicated in Section IV further, the unsymmetric convention is used for the normalization of activity coefficients, thereby avoiding all arbitrary extrapolations to find the properties of hypothetical states finally, a flexible two-parameter model is used for describing the effect of composition and temperature on liquid-phase activity coefficients. The flexibility of the model necessarily requires some binary data over a range of composition and temperature to obtain the desired accuracy, especially in the critical region, more binary data are required for Chueh s method than for that of Chao and Seader (Cl). Fortunately, reliable data for high-pressure equilibria are now available for a variety of binary mixtures of nonpolar fluids, mostly hydrocarbons. Chueh s method, therefore, is primarily applicable to equilibrium problems encountered in the petroleum, natural-gas, and related industries. 

The application of TST to account for the effects on the rate constants of reactions involving ions is one of its notable successes. We shall make use of a previously derived equation [Eq. (9-28)] to relate changes in the rate constant to the activity coefficients ... 

The activity coefficients in the above equation may be determined by obtaining experimental data for D/ m, and relating Ka to the free energy change for the reaction using the equation AGf = -RT In K. 

Note that a number of complicating factors have been left out for clarity For instance, in the EMF equation, activities instead of concentrations should be used. Activities are related to concentrations by a multiplicative activity coefficient that itself is sensitive to the concentrations of all ions in the solution. The reference electrode necessary to close the circuit also generates a (diffusion) potential that is a complex function of activities and ion mobilities. Furthermore, the slope S of the electrode function is an experimentally determined parameter subject to error. The essential point, though, is that the DVM-clipped voltages appear in the exponent and that cheap equipment extracts a heavy price in terms of accuracy and precision (viz. quantization noise such an instrument typically displays the result in a 1 mV, 0.1 mV, 0.01 mV, or 0.001 mV format a two-decimal instrument clips a 345.678. .. mV result to 345.67 mV, that is it does not round up ... 78 to ... 8 ). 

Amongst the various relations given above, the last is the one which is actually used to calculate activity coefficients. For calculating the activity coefficient of the component, A, when the activity coefficient versus composition relationship for the other component, B, is known, the equation is used in the following form ... 

Most measurements include the determination of ions in aqueous solution, but electrodes that employ selective membranes also allow the determination of molecules. The sensitivity is high for certain ions. When specificity causes a problem, more precise complexometric or titri-metric measurements must replace direct potentiometry. According to the Nernst equation, the measured potential difference is a measure of the activity (rather than concentration) of certain ions. Since the concentration is related to the activity through an appropriate activity coefficient, calibration of the electrode with known solution(s) should be carried out under conditions of reasonable agreement of ionic strengths. For quantitation, the standard addition method is used. 

In many drug solutions, it is necessary to use buffer salts in order to maintain the formulation at the optimum pH. These buffer salts can affect the rate of drug degradation in a number of ways. First, a primary salt effect results because of the effect salts have on the activity coefficient of the reactants. At relatively low ionic strengths, the rate constant, k, is related to the ionic strength, p, according to... 

---