Solvation energy surface density

Smooth COSMO solvation model. We have recently extended our smooth COSMO solvation model with analytical gradients [71] to work with semiempirical QM and QM/MM methods within the CHARMM and MNDO programs [72, 73], The method is a considerably more stable implementation of the conventional COSMO method for geometry optimizations, transition state searches and potential energy surfaces [72], The method was applied to study dissociative phosphoryl transfer reactions [40], and native and thio-substituted transphosphorylation reactions [73] and compared with density-functional and hybrid QM/MM calculation results. The smooth COSMO method can be formulated as a linear-scaling Green s function approach [72] and was applied to ascertain the contribution of phosphate-phosphate repulsions in linear and bent-form DNA models based on the crystallographic structure of a full turn of DNA in a nucleosome core particle [74],... 

From the chemical point of view, the solvent in which the CL experiment is carried out can have a dramatic influence on the efficiency of the CL reaction as solvation can alter the shapes, the depths, and the densities of the vibrational states of the potential surfaces representing the ground states of products and reactants and the lowest excited singlet state of the potential fluorophore. The alteration of the intersections of these potential energy surfaces can affect the enthalpies of reaction and the enthalpies of activation for dark and lumigenic reactions. In some cases, these changes will favor CL (if AH decreases relative to AHa) and in some cases, they will make it thermodynamically unfavorable for CL to occur. 

The quantities defined by Eqs. (2)—(7) plus Vs max, Vs min, and the positive and negative areas, A and, enable detailed characterization of the electrostatic potential on a molecular surface. Over the past ten years, we have shown that subsets of these quantities can be used to represent analytically a variety of liquid-, solid-, and solution-phase properties that depend on noncovalent interactions [14-17, 84] these include boiling points and critical constants, heats of vaporization, sublimation and fusion, solubilities and solvation energies, partition coefficients, diffusion constants, viscosities, surface tensions, and liquid and crystal densities. 

It has been our experience that 7s(r) and Vs(r) play different but complementary roles with respect to molecular reactivity [71,83-85], Vs(r) is effective for treating noncova-lent interactions, which are primarily electrostatic in nature [74,86-89], For instance, a variety of condensed-phase physical properties - boiling points, critical constants, heats of phase transitions, solubilities and solvation energies, partition coefficients, surface tensions, viscosities, diffusion constants and densities - can be expressed quantitatively in terms of one or more key features of Vs(r), such as its maximum and minimum, average deviation, positive and negative variances, etc. [80,90-92], Hydrogen bond donating... 

Equation (11.20) is the primary PCM equation. It must be discretized for actual computation (see Section 11.2.2), but then given the solute s electrostatic potential evaluated at the surface discretization points, this equation can be solved for the induced surface charge at those points (i.e.,the discretized a). In an MM/PCM calculation, the electrostatic solvation energy is then immediately available via a discretized version of Eq. (11.3), although in QM applications the surface charge must be included in the next self-consistent field (SCF) iteration, and the SCF procedure is iterated until both the electron density and the surface charge have reached mutual self-consistency. 

Since an SAS is computationally more expensive to generate than a van der Waals surface, and since the difference is often small, a van der Waals surface is often used in practice. Furthermore, a very small displacement of an atom may alter the SAS in a discontinuous fashion, as a cavity suddenly becomes too small to allow a solvent molecule to enter. Alternatively, the cavity may be calculated directly from the wave function, for example by taking a surface corresponding to an electron density of 0.001. It is generally found that the shape of the hole is importan, and that molecular shaped cavities are necessary to be able to obtain good agreement with experimental data (such as solvation energies). It should be emphasized, however, that reaction field... 

The main outcome of this theory is the stable bundle size k that is obtained by minimization oifik). The bundle size can be evaluated as a function of geometric parameters of single rods, vr and Lr, rod density, p, surface charge density, a, surface tension, /, solvation energy, Esoiv and the dissociation parameter, rj. 

---