Water structure continuum models

The elucidation of actinide chemistry in solution is important for understanding actinide separation and for predicting actinide transport in the environment, particularly with respect to the safety of nuclear waste disposal.72,73 The uranyl CO + ion, for example, has received considerable interest because of its importance for environmental issues and its role as a computational benchmark system for higher actinides. Direct structural information on the coordination of uranyl in aqueous solution has been obtained mainly by extended X-ray absorption fine structure (EXAFS) measurements,74-76 whereas X-ray scattering studies of uranium and actinide solutions are more rare.77 Various ab initio studies of uranyl and related molecules, with a polarizable continuum model to mimic the solvent environment and/or a number of explicit water molecules, have been performed.78-82 We have performed a structural investigation of the carbonate system of dioxouranyl (VI) and (V), [U02(C03)3]4- and [U02(C03)3]5- in water.83 This study showed that only minor geometrical rearrangements occur upon the one-electron reduction of [U02(C03)3]4- to [U02(C03)3]5-, which supports the reversibility of this reduction. 

For the description of a solution of alanine in water two models were compared and combined with one another (79), namely the continuum model approach and the cluster ansatz approach (148,149). In the cluster approach snapshots along a trajectory are harvested and subsequent quantum chemical analysis is carried out. In order to learn more about the structure and the effects of the solvent shell, the molecular dipole moments were computed. To harvest a trajectory and for comparison AIMD (here CPMD) simulations were carried out (79). The calculations contained one alanine molecule dissolved in 60 water molecules. The average dipole moments for alanine and water were derived by means of maximally localized Wannier functions (MLWF) (67-72). For the water molecules different solvent shells were selected according to the three radial pair distributions between water and the functional groups. An overview about the findings is given in Tables II and III. 

Much like the RISM method, the LD approach is intermediate between a continuum model and an explicit model. In the limit of an infinite dipole density, the uniform continuum model is recovered, but with a density equivalent to, say, the density of water molecules in liquid water, some character of the explicit solvent is present as well, since the magnitude of the dipoles and their polarizability are chosen to mimic the particular solvent (Papazyan and Warshel 1997). Since the QM/MM interaction in this case is purely electrostatic, other non-bonded interaction terms must be included in order to compute, say, solvation free energies. When the same surface-tension approach as that used in many continuum models is adopted (Section 11.3.2), the resulting solvation free energies are as accurate as those from pure continuum models (Florian and Warshel 1997). Unlike atomistic models, however, the use of a fixed grid does not permit any real information about solvent structure to be obtained, and indeed the fixed grid introduces issues of how best to place the solute into the grid, where to draw the solute boundary, etc. These latter limitations have curtailed the application of the LD model. 

Considerable progress has been made in going beyond the simple Debye continuum model. Non-Debye relaxation solvents have been considered. Solvents with nonuniform dielectric properties, and translational diffusion have been analyzed. This is discussed in Section II. Furthermore, models which mimic microscopic solute/solvent structure (such as the linearized mean spherical approximation), but still allow for analytical evaluation have been extensively explored [38, 41-43], Finally, detailed molecular dynamics calculations have been made on the solvation of water [57, 58, 71]. 

DFT was employed to study the mechanism of ammonolysis of phenyl formate in the gas phase, and the effect of various solvents on the title reaction was assessed by the polarizable continuum model (PCM). The calculated results show that the neutral concerted pathway is the most favourable one in the gas phase and in solution.24 The structure and stability of putative zwitterionic complexes in the ammonolysis of phenyl acetate were examined using DFT and ab initio methods by applying the explicit, up to 7H20, and implicit PCM solvation models. The stability of the zwitterionic tetrahedral intermediate required an explicit solvation by at least five water molecules with stabilization energy of approximately 35 kcalmol-1 25... 

Figure 3.17 Comparison between experiment (dashed curve) and calculations combining the polarizable continuum model for solute electronic structure and continuum dielectric theory of solvation dynamics in water. SRF(t) stands for S(t) in our notation. The calculations are for a cavity based on a space-filling model of Cl53, while the experiments are for C343. The two sets of theoretical results correspond to using water e(o>) from simulation (full curve) of SPC/E water and from a fit to experimental data (dash-dotted curve). (Reprinted from F. Ingrosso, A. Tani andJ. Tomasi, J. Mol. Liq., 1117, 85-92. Copyright (2005), with permission from Elsevier).

---