Coupling approaches solvation free energy

Most of these extensions have involved electron correlation methods based on variational approaches (DFT, MCSCF, CI,VB). These methods can be easily formulated by optimizing the free energy functional (1.117), expressed as a function of the appropriate variational parameters, as in the case of the HF approximation. In contrast, for nonva-riational methods such as the Moller-Plesset theory or Coupled-Cluster, the parallel extension to solvation model is less straightforward. 

Density Functional Theory does not require specific modifications, in relation to the solvation terms [9], with respect to the Hartree-Fock formalism presented in the previous section. DFT also absorbs all the properties of the HF approach concerning the analytical derivatives of the free energy functional (see also the contribution by Cossi and Rega), and as a matter of fact continuum solvation methods coupled to DFT are becoming the routine approach for studies of solvated systems. 

Although the correlative methods based on the coupled-cluster (CC) ansatz are among the most accurate approaches for molecules in vacuum, their extension to introduce the interactions between a molecule and a surrounding solvent have not yet reached a satisfactory stage. The main complexity in coupling CC to solvation methods comes from the evaluation of the electronic density, or of the related observables, needed for the calculation of the reaction field. Within the CC scheme the electronic density can only be evaluated by a relaxed approach, which implies the evaluation of the first derivative of the free energy functional. As discussed previously for the cases of the Cl and MPn approaches, this leads to a more involved formalism. 

Solvent effects can significantly influence the function and reactivity of organic molecules.1 Because of the complexity and size of the molecular system, it presents a great challenge in theoretical chemistry to accurately calculate the rates for complex reactions in solution. Although continuum solvation models that treat the solvent as a structureless medium with a characteristic dielectric constant have been successfully used for studying solvent effects,2,3 these methods do not provide detailed information on specific intermolecular interactions. An alternative approach is to use statistical mechanical Monte Carlo and molecular dynamics simulation to model solute-solvent interactions explicitly.4 8 In this article, we review a combined quantum mechanical and molecular mechanical (QM/MM) method that couples molecular orbital and valence bond theories, called the MOVB method, to determine the free energy reaction profiles, or potentials of mean force (PMF), for chemical reactions in solution. We apply the combined QM-MOVB/MM method to... 

Several problems are encountered when potentials in different solvents are sought compared to the potential scale in water. A variety of approaches [186-194] have been followed to attack this problem usually the approach has been to introduce some kind of nonthermodynamic assumption, such as the supposition that certain large, monovalent ions (Rb, CS ) [191] or redox systems [186-188] of the charge type n/n + 1 (preferably 0/ + 1 [186,187]) have a nearly equal free energy of solvation in the two solvents so that the free energy of a transfer of the reference ion is small. The redox couples [194] ferro-cenium/ferrocene and bis(biphenyl)chromium(I)/bis(biphenyl)chromium(0) (BCr" "/BCr) have been recommended as reference redox systems for measurements in nonaqueous solvents however, an investigation concluded [195] that the electrochemistry of ferrocene in MeCN at microelectrodes was far from ideal, as some film formation may occur. 

For distances larger than q, the usual D-H treatment applies to ions that are considered free. The cut-off distance is reasonable, because at that separation distance, the thermal energy is half the work necessary to separate the ion-pair. The ions in the couple cannot get closer than a distance of closest approach, a, which is at least the sum of the crystal radii, but it is not larger than the sum of the solvated ion radii. Since the neat charge of the duplex is zero, the ion-pair is not acted upon by coulombic fields, it does not migrate, and it does therefore not contribute to the electrical conductivity of the solution. 

---