Stability and Equilibrium Calculations

Substances that are not complexes, for example, methylmercury, HgCHa, will not undergo ligand replacement reactions. 

For the same reaction, if Fe + is the metal and Cl is the ligand several hours are required, but if Fe has an OH attached to it, that is, FeOH, only minutes are required. The reaction between Fe and S04 or SON , on the other hand, requires only minutes. 

Notable exceptions to the statement that coordination reactions occur rapidly are the reactions involving replacement of H2O from CrfHaOle by ligands such as SO4 , which take days to years. Changes in the structure of polymeric metal hydroxide species may take weeks. 

Equilibrium constants for-complexes are usually stated for reactions written in the direction of complex formation, that is. 

Ligand + central metal ion complex For example, for the ammine copper complex. 

---

The results of ab initio calculations are compatible with the conclusion that -SH stabilizes an adjacent carbanion more than —OH363 Furthermore, the above discussion represents a viable alternative to the d-orbital model for explaining the enhanced kinetic and equilibrium acidities of molecules containing sulfur groups. 

In the preceding discussion we considered equilibrium void stability however, actual processing conditions involve changing temperature and pressure with time. Whereas equilibrium calculations provide bounds on void growth, it is the time-dependent growth process that is most important from a product quality viewpoint. 

Where I is the initial amount (g) of excipient in the solution, Wc is the weight of the container (g), Fs is the solution volume (L), and E, is the equilibrium partitioning constant, the ratio of the concentration of solute in the film to that in water, at equilibrium (66). This can be calculated from the more familiar, and referenced, solvent solvent partition coefficients. Plastics and rubber stoppers can also leach stabilizers and plasticizers into the contained injection volume. The extent of this can be calculated by considering the same factors described above. 

The charge density, Volta potential, etc., are calculated for the diffuse double layer formed by adsorption of a strong 1 1 electrolyte from aqueous solution onto solid particles. The experimental isotherm can be resolved into individual isotherms without the common monolayer assumption. That for the electrolyte permits relating Guggenheim-Adam surface excess, double layer properties, and equilibrium concentrations. The ratio u0/T2N declines from two at zero potential toward unity with rising potential. Unity is closely reached near kT/e = 10 for spheres of 1000 A. radius but is still about 1.3 for plates. In dispersions of Sterling FTG in aqueous sodium ff-naphthalene sulfonate a maximum potential of kT/e = 7 (170 mv.) is reached at 4 X 10 3M electrolyte. The results are useful in interpretation of the stability of the dispersions. 

More refined continuum models—for example, the well-known Fumi-Tosi potential with a soft core and a term for attractive van der Waals interactions [172]—have received little attention in phase equilibrium calculations [51]. Refined potentials are, however, vital when specific ion-ion or ion-solvent interactions in electrolyte solutions affect the phase stability. One can retain the continuum picture in these cases by using modified solvent-averaged potentials—for example, the so-called Friedman-Gumey potentials [81, 168, 173]. Specific interactions are then represented by additional terms in (pap(r) that modify the ion distribution in the desired way. Finally, there are models that account for the discrete molecular nature of the solvent—for example, by modeling the solvent as dipolar hard spheres [174, 175]. 

Ginsburg and Soloviev (1998, pp. 150-151) state that the BSR is the most widely used indirect indication of gas hydrates. The most important evidence of the hydrate caused nature of the BSR is the coincidence of temperature and pressure calculated at it s depth with the equilibrium temperatures and pressure of gas hydrate stability. The association with the base of the hydrate stability zone is beyond question. ... 

The keto-enol equilibrium of the 1,3-diketones has been the subject of intensive studies using various physical techniques and theoretical calculations [78-80], Recently, X-ray crystal analysis of acetylacetone (83) was carried out at 110 K, and it was found that it exists as an equilibrium mixture of the two enol forms 83b and 83c [81]. Room-temperature studies show an acetylacetone molecule with the enolic H-atom centrally positioned, which can be attributed to the dynamically averaged structure 83d. Application of a crystal engineering technique showed that a 1 1 inclusion complex of83 can be formed with l,l/-binaphthyl-2,2/-dicarboxylic acid in which the enol form is stabilized by a notably short intramolecular hydrogen bond [82],... 

Bursi et al. (2001) reported two methods to calculate the stability of testosterone-like steroids. These were the use of a decision tree and molecular descriptors or quantum mechanical methods. For satisfactory accuracy, Bursi and colleagues (2001) had to use a 3-21G basis set with spin correction and equilibrium geometries. This required 12 hours of computation for reactants and products. Optimization of transition state geometry was also required. The simpler decision tree analysis approach indicated that descriptors such as the volumes and, to a lesser degree, the shape were important. Correlations of calculated and experimental rates of metabolism were reported. 

Christiansen et al. (54) applied the Naphtali-Sandholm method to natural gas mixtures. They replaced the equilibrium relationships and component vapor rates with the bubble-point equation and total liquid rate to get practically half the number of functions and variables [to iV(C + 2)]. By exclusively using the Soave-Redlich-Kwong equation of state, they were able to use analytical derivatives of revalues and enthalpies with respect to composition and temperature. To improve stability in the calculation, they limited the changes in the independent variables between trials to where each change did not exceed a preset maximum. There is a Naphtali-Sandholm method in the FraChem program of OLI Systems, Florham Park, New Jersey CHEMCAD of Coade Inc, of Houston, Texas PRO/II of Simulation Sciences of Fullerton, California and Distil-R of TECS Software, Houston, Texas. Variations of the Naphtali-Sandholm method are used in other methods such as the homotopy methods (Sec. 4,2.12) and the nonequilibrium methods (Sec. 4.2.13). 

---