Organic molecule, solvation

Figure 9.4 An organic molecule solvated in a cage inside liquid water.

The oldest of the numerous syntheses of dichlorophosphates dates back to 1911 (S7) it makes use of the reaction of metal oxides (MgO, CaO, MnO) with POCI3 in presence of solvatating organic molecules (e.g. ethyl acetate = E). This reaction can be understood by the access of moisture from the air (26) as well as by the simultaneous formation of metal chloride (27) ... 

Thermochemistry. Chen et al.168 combined the Kohn-Sham formalism with finite difference calculations of the reaction field potential. The effect of mobile ions into on the reaction field potential Poisson-Boltzman equation. The authors used the DFT(B88/P86)/SCRF method to study solvation energies, dipole moments of solvated molecules, and absolute pKa values for a variety of small organic molecules. The list of molecules studied with this approach was subsequently extended182. A simplified version, where the reaction field was calculated only at the end of the SCF cycle, was applied to study redox potentials of several iron-sulphur clusters181. 

To calculate free energies of solvation for several organic molecules, Fortunelli and Tomasi applied the boundary element method for the reaction field in DFT/SCRF framework173. The authors demonstrated that the DFT/SCRF results obtained with the B88 exchange functional and with either the P86 or the LYP correlation functional are significantly closer to the experimental ones than the ones steming from the HF/SCRF calculations. The authors used the same cavity parameters for the HF/SCRF and DFT/SCRF calculations, which makes it possible to attribute the apparent superiority of the DFT/SCRF results to the density functional component of the model. The boundary element method appeared to be very efficient computationally. The DFT/SCRF calculations required only a few percent more CPU time than the corresponding gas-phase SCF calculations. 

As we have recently shown [3], the stability of GIC with Bronstcd acids is governed by ionization potential of the intercalated anion. Solvation of the anion substantially increases ionization potential and hence stabilizes GIC. Such a possibility was predicted by Inagaki [4] who proposed to use organic molecules as additional cointercalants for the purpose of stabilization. Afterwards, we successfully used glacial acetic acid and water to synthesize stable products [5],... 

In recent years the FEP method has fallen into disuse. However, as the studies outlined above show, in many cases the results obtained are in good agreement with experimental measurements. In these cases new information may be obtained, which may be difficult or even impossible to measure. Examples of this are the relative ratios of conformers in the histamine system, a detailed breakdown of the tautomers present in the guanine or cystine systems, or the acidity strengths of organic molecules such as ethane in water. In addition to this thermodynamic data, the simulations then also provide detailed information on the solvation of the species of interest. 

The level of accuracy that can be achieved by these different methods may be viewed as somewhat remarkable, given the approximations that are involved. For relatively small organic molecules, for instance, the calculated AGsoivation is now usually within less than 1 kcal/mole of the experimental value, often considerably less. Appropriate parametrization is of key importance. Applications to biological systems pose greater problems, due to the size and complexity of the molecules,66 156 159 161 and require the use of semiempirical rather than ab initio quantum-mechanical methods. In terms of computational expense, continuum models have the advantage over discrete molecular ones, but the latter are better able to describe solvent structure and handle first-solvation-shell effects. 

QM calculations provide an accurate way to treat strong, long-range electrostatic forces that dominate many solvation phenomena. The errors due to the use of the approximate eontinuum solvation models ean be small enough to allow the quantitative treatment of the solute behavior therefore, this approach is widely used also for evaluating solvent extraetion equilibria of organic molecules [56]. 

It is interesting to note that smaller ions (e.g., Na, Mg, Ca, Cl") form hydration shells larger than bigger ions, which tend to bind water molecules only very weakly. In a simple way, the salting out of nonpolar and weakly polar compounds was explained by Schwarzenbach et al. (2003) by imagining that the dissolved ions compete successfully with the organic compound for solvent molecules. The freedom of some water molecules to solvate an organic molecule depends on the type and concentration of salts. 

Fig. 6. Calculated (y-axis) and experimental (x-axis) molecular properties. Top Left free energy of solvation in kcal/mol for 291 small organic molecules. Top Right boiling points (in Kelvin) of 298 small organic molecules. Bottom Left the log concentration ratio between the blood and brain for 75 compounds. Bottom Right the solubility in water of 1438 small organic compounds (units are log concentration ratios).

What is the exact influence of water and organic molecules on the enzyme structure Could its effects on properties such as selectivity, affinity, binding constants, and catalytic constants be predictable by controlling the hydration/solvation state ... 

Thermochemical data for the solvation of ions as used in the preceding calculations are difficult to measure and even to estimate. Therefore this kind of calculation of AH° for ionic reactions involving organic molecules in solution usually cannot be made. As a result, we have considerably fewer possibilities to assess the thermodynamic feasibility of the individual steps of polar reactions in solution than we do of vapor-phase radical processes. Bond energies are not of much use in predicting or explaining reactivity in ionic reactions unless we have information that can be used to translate gas-phase AH°. values to solution AH° values. Exercise 8-3 will give you a chance to see how this is done. 

Electron transfer reactions have been characterized with much more rigor in inorganic chemistry than with organic molecules. Marcus has provided the principal description relating the kinetics and thermodynamics of electron transfer between metal complexes (1). The Marcus theory, a computationally simple approach with good predictive power, is an empirical treatment which uses thermodynamic parameters and spectroscopic measurements to calculate kinetic data. It assumes that bimolecular electron transfer reactions occur in three stages as shown in Scheme 1 (1) formation of the precursor complex, (2) electron transfer, and (3) solvation of the redox pair. 

---