Self-consistent Reaction Field Methods Cavities

Combined Quantum Mechanical and Molecular Mechanical Potentials Combined Quantum Mechanics and Molecular Mechanics Approaches to Chemical and Biochemical Reactivity Continuum Solvation Density Functional Theory (DFT), Hartree-Fock (HF), and the Self-consistent Field MNDO Monte Carlo Simulations for Liquids Quantum Mechanical/Molecular Mechanical (QM/MM) Coupled Potentials Quantum Mechanics/Molecular Mechanics (QM/MM) Self-consistent Reaction Field Methods Self-consistent Reaction Field Methods Cavities Solvation Modeling TURBOMOLE. 

It may be noted here that the theoretical developments below apply also to other solution chemistry problems of current interest. Electronic structure calculations for species in solution often append to the free energy a contribution to represent cavitation (see Self-consistent Reaction Field Methods Cavities). This is a free energy cost to open space for the solute in the solvent. Theories that solve that problem for water should be carried over to other solvents too. 

SELF-CONSISTENT REACTION FIELD METHODS CAVITIES... 

Self-consistent Reaction Field Methods Cavities... 

In the thermodynamic cycle (1.92), the Gibbs energies of transfer AG° (i) can be calculated by theoretical methods. Two main classes of methods have been developed for modelling solvent effects. Molecular dynamics and Monte Carlo methods use a discrete representation of the solvent molecules whereas, in the second class of so-called SCRF (self-consistent reaction field) methods, the solvent is represented as a dielectric continuum surrounding the solute cavity. These methods are outside the scope of this book. They are described in reviews [105, 107, 108] and books [109, 110]. Examples of their application to Lewis acid/base complexes can be found in the following references hydrogen-bonded complexes [111-113], BF3 and BH3 complexes [114] and diiodine complexes [115]. 

One femily of models for systems in non-aqueous solution are referred to as Self-Consistent Reaction Field (SCRF) methods. These methods all model the solvent as a continuum of uniform dielectric constant e the reaction field. The solute is placed into a cavity within the solvent. SCRF approachs differ in how they define the cavity and the reaction field. Several are illustrated below. 

Conceptually, the self-consistent reaction field (SCRF) model is the simplest method for inclusion of environment implicitly in the semi-empirical Hamiltonian24, and has been the subject of several detailed reviews24,25,66. In SCRF calculations, the QM system of interest (solute) is placed into a cavity within a polarizable medium of dielectric constant e (Fig. 2.2). For ease of computation, the cavity is assumed to be spherical and have a radius ro, although expressions similar to those outlined below have been developed for ellipsoidal cavities67. Using ideas from classical electrostatics, we can show that the interaction potential can be expressed as a function of the charge and multipole moments of the solute. For ease... 

Organometallic systems such as porphyrines have been investigated because of the possibility to fine tune their response by functionalization[105-107]. Systems of increased the dimensionality have been of particular interest [108-111], Concomitant to the large effort to establish useful structure-to-properties relationships, considerable effort has now been put to investigate the environmental effects on TPA[112-114], For example, the solvent effect has been studied for a small linear push-pull chromophore using a self-consistent reaction field (homogeneous solvation) method employing a spherical cavity and an internal force field (IFF) method[l 12] in another study the polarizable continuum model has been employed to calculate the relevant quantities to obtain the TPA cross-section in the limit of a two-state model[113] Woo et al. made a critical study of experimental comparison of TPA cross-sections in different solvents[114]. 

Among the few determinations of of molecular crystals, the CPHF/ INDO smdy of Yamada et al. [25] is unique because, on the one hand, it concerns an open-shell molecule, the p-nitrophenyl-nitronyl-nitroxide radical (p-NPNN) and, on the other hand, it combines in a hybrid way the oriented gas model and the supermolecule approach. Another smdy is due to Luo et al. [26], who calculated the third-order nonlinear susceptibility of amorphous thinmultilayered films of fullerenes by combining the self-consistent reaction field (SCRF) theory with cavity field factors. The amorphous namre of the system justifies the choice of the SCRF method, the removal of the sums in Eq. (3), and the use of the average second hyperpolarizability. They emphasized the differences between the Lorentz Lorenz local field factors and the more general Onsager Bbttcher ones. For Ceo the results differ by 25% but are in similar... 

A short overview of the quantum chemical and statistical physical methods of modelling the solvent effects in condensed disordered media is presented. In particular, the methods for the calculation of the electrostatic, dispersion and cavity formation contributions to the solvation energy of electroneutral solutes are considered. The calculated solvation free energies, proceeding from different geometrical shapes for the solute cavity are compared with the experimental data. The self-consistent reaction field theory has been used for a correct prediction of the tautomeric equilibrium constant of acetylacetone in different dielectric media,. Finally, solvent effects on the molecular geometry and charge distribution in condensed media are discussed. 

Recently, a new category of methods, the cavity model, has been proposed to account for the solvent effect. Molecules or supermolecules are embedded in a cavity surrounded by a dielectric continuum, the solvent being represented by its static dielectric constant. The molecules (supermolecules) polarize the continuum. As a consequence this creates an electrostatic potential in the cavity. This reaction potential interacts with the molecules (supermolecules). This effect can be taken into account through an interaction operator. The usual SCF scheme is modified into a SCRF (self consistent reaction field) scheme, and similar modifications can be implemented beyond the SCF level. Several studies based on this category of methods have been published on protonated hydrates. They account for the solvent effect on the filling of the first solvation shell (53, 69), the charges (69, 76) and the energy barrier to proton transfer (53, 76). 

The quantum Onsager model, which has also been termed the Self-Consistent Reaction Field (SCRF) method, is the simplest of the continuum models used in solvation studies. In this model, which dates from the work of Kirkwood[44] and Onsager[45] in the 1930s, the solvent is represented by a continuous rmifonn dielectric with a static dielectric constant, e, surrounding a solute in a spherical cavity[46] - [48]. 

Reaction field methods model solutions by placing the solute in a cavity of a polarizable medium. The electrostatic potential due to the solute molecule polarizes the surrounding medium which in turn changes the charge distribution of the solute. Hence, the electrostatic interaction has to be evaluated self-consistently (self-consistent reaction field, SCRF). A term for creating the cavity (calculated from the surface of the cavity) has a be added to the solvation energy. Explicit treatment of solvent molecules can be combined with a reaction field method. 

---