Parameter optimization in differential geometry based solvation models

Differential geometry (DG) based solvation models are a new class of variational implicit solvent approaches that are able to avoid unphysical solvent-solute boundary definitions and associated geometric singularities, and dynamically couple polar and nonpolar interactions in a self-consistent framework. Our earlier study indicates that DG based nonpolar solvation model outperforms other methods in nonpolar solvation energy predictions. However, the DG based full solvation model has not shown its superiority in solvation analysis, due to its difficulty in parametrization, which must ensure the stability of the solution of strongly coupled nonlinear Laplace-Beltrami and Poisson-Boltzmann equations. In this work, we introduce new parameter learning algorithms based on perturbation and convex optimization theories to stabilize the numerical solution and thus achieve an optimal parametrization of the DG based solvation models. An interesting feature of the present DG based solvation model is that it provides accurate solvation free energy predictions for both polar and nonploar molecules in a unified formulation. Extensive numerical experiment demonstrates that the present DG based solvation model delivers some of the most accurate predictions of the solvation free energies for a large number of molecules.Comment: 19 pages, 12 figures, convex optimizatio

Mini-grand canonical ensemble: chemical potential in the solvation shell

Quantifying the statistics of occupancy of solvent molecules in the vicinity of solutes is central to our understanding of solvation phenomena. Number fluctuations in small `solvation shells' around solutes cannot be described within the macroscopic grand canonical framework using a single chemical potential that represents the solvent `bath'. In this communication, we hypothesize that molecular-sized observation volumes such as solvation shells are best described by coupling the solvation shell with a mixture of particle baths each with its own chemical potential. We confirm our hypotheses by studying the enhanced fluctuations in the occupancy statistics of hard sphere solvent particles around a distinguished hard sphere solute particle. Connections with established theories of solvation are also discussed

Electrostatic solvation free energies of charged hard spheres using molecular dynamics with density functional theory interactions

Determining the solvation free energies of single ions in water is one of the most fundamental problems in physical chemistry and yet many unresolved questions remain. In particular, the ability to decompose the solvation free energy into simple and intuitive contributions will have important implications for models of electrolyte solution. Here, we provide definitions of the various types of single ion solvation free energies based on different simulation protocols. We calculate solvation free energies of charged hard spheres using density functional theory interaction potentials with molecular dynamics simulation (DFT-MD) and isolate the effects of charge and cavitation, comparing to the Born (linear response) model. We show that using uncorrected Ewald summation leads to unphysical values for the single ion solvation free energy and that charging free energies for cations are approximately linear as a function of charge but that there is a small non-linearity for small anions. The charge hydration asymmetry (CHA) for hard spheres, determined with quantum mechanics, is much larger than for the analogous real ions. This suggests that real ions, particularly anions, are significantly more complex than simple charged hard spheres, a commonly employed representation.Comment: 28 pages, 5 figure

Segue Between Favorable and Unfavorable Solvation

Solvation of small and large clusters are studied by simulation, considering a range of solvent-solute attractive energy strengths. Over a wide range of conditions, both for solvation in the Lennard-Jones liquid and in the SPC model of water, it is shown that the mean solvent density varies linearly with changes in solvent-solute adhesion or attractive energy strength. This behavior is understood from the perspective of Weeks' theory of solvation [Ann. Rev. Phys. Chem. 2002, 53, 533] and supports theories based upon that perspective.Comment: 8 pages, 7 figure

Solvation forces in Ising films with long-range boundary fields: density-matrix renormalization-group study

Using the quasi-exact density-matrix renormalization-group method we calculate the solvation forces in two-dimensional Ising films of thickness L subject to identical algebraically decaying boundary fields with various decay exponents p. At the bulk critical point the solvation force acquires a universal contribution which is long-ranged in L due to the critical fluctuations, a phenomenon known as the critical Casimir effect. For p = 2, 3 and 50, we study the scaling behaviour of the solvation force along the pseudo-phase coexistence and along the critical and sub-critical isotherms.Comment: 9 pages, 6 figures, accepted to Molecular Physic

Theoretical studies of 31P NMR spectral properties of phosphanes and related compounds in solution

Selected theoretical methods, basis sets and solvation models have been tested in their ability to predict 31P NMR chemical shifts of large phosphorous-containing molecular systems in solution. The most efficient strategy was found to involve NMR shift calculations at the GIAO-MPW1K/6-311++G(2d,2p)//MPW1K/6-31G(d) level in combination with a dual solvation model including the explicit consideration of single solvent molecules and a continuum (PCM) solvation model. For larger systems it has also been established that reliable 31P shift predictions require Boltzmann averaging over all accessible conformations in solution

Comparative study on the gas to solution phase solvation free energies of model combustion flue gas compounds (N~2~, O~2~, CO~2~, H~2~O, SO~2~, and CO) in 178 organic solvents using the IEFPCM-UFF, CPCM, and SMD implicit solvent models at the Gaussian-4 (G4) level of theory

Gas to solution phase Gibbs free energies of solvation at 298.15 K for transfer of six representative combustion flue gas compounds (N~2~, O~2~, CO~2~, H~2~O, SO~2~, and CO) were calculated at the Gaussian-4 (G4) level of theory using the IEFPCM-UFF, CPCM, and SMD implicit solvent models for 178 organic solvents. The IEFPCM-UFF and CPCM models yield similar free energies of solvation for all six compounds in each of the solvents considered, having maximum absolute intra-solvent deviations <1.6 kJ mol^-1^. Substantial free energy of solvation differences were observed between the IEFPCM-UFF/CPCM and SMD models, with maximum absolute intra-solvent deviations up to 45.5 kJ mol^-1^. The IEFPCM-UFF and CPCM models displayed strong free energy of solvation correlations with the solvent dielectric constant for each compound, whereas the SMD model exhibits a significantly more variable free energy of solvation relationship with the solvent dielectric constant