Infinite dilution activity coefficient from SMD calculations: accuracy and performance for predicting liquid-liquid equilibria

Single-ion solvation free energy: A new cluster-continuum approach based on the cluster expansion method

Accurate calculation of the solvation free energy of single ions remains an important goal, involving development in the dielectric continuum solvation models, and statistical mechanics with explicit solvent and hybrid discrete-continuum methods. In the last case, many of the research studies involve a quasi-chemical approach using the monomer cycle or the cluster cycle to calculate the solvation free energy of single ions. In this work, a different cluster-continuum approach based on the cluster expansion method was tested for solvation of 16 cations and 32 anions in aqueous solution. The SMD model was used for the dielectric continuum part and three explicit water molecules were introduced in the region of the solute with the highest interaction energy. Harmonic frequency calculations and molecular dynamics sampling of configurations are not required. An empirical γ_N parameter for cations and another for anions is introduced. The method produces a substantial improvement of the SMD model with a mean absolute deviation of 2.3 kcal mol^-1 for cations and 2.9 kcal mol^-1 for anions. The analysis of the correlation between theoretical and experimental data produces a linear regression line with a slope of 1.09 for cations and 1.01 for anions. The good results of this approximated cluster expansion approach suggest that the method could be further improved by including more solvent molecules and sampling the configurations.

Assessment of the SM12, SM8, and SMD Solvation Models for Predicting Limiting Activity Coefficients at 298.15 K

The SMx (x = 12, 8, or D) universal solvent models are implicit solvent models which using electronic structure calculations can compute solvation free energies at 298.15 K. While solvation free energy is an important thermophysical property, within the thermodynamic modeling of phase equilibrium, limiting (or infinite dilution) activity coefficients are preferred since they may be used to parameterize excess Gibbs free energy models to model phase equilibrium. Conveniently, the two quantities are related. Therefore the present study was performed to assess the ability to use the SMx universal solvent models to predict limiting activity coefficients. Two methods of calculating the limiting activity coefficient where compared: (1) the solvation free energy and self-solvation free energy were both predicted and (2) the self-solvation free energy was computed using readily available vapor pressure data. Overall the first method is preferred as it results in a cancellation of errors, specifically for the case in which water is a solute. The SM12 model was compared to both the Universal Quasichemical Functional-group Activity Coefficients (UNIFAC) and modified separation of cohesive energy density (MOSCED) models. MOSCED was the highest performer, yet had the smallest available compound inventory. UNIFAC and SM12 exhibited comparable performance. Therefore further exploration and research should be conducted into the viability of using the SMx models for phase equilibrium calculations.

Infinite Dilution Activity Coefficients as Constraints for Force Field Parametrization and Method Development

Molecular simulations begin with an underlying energy model or force field and from this can predict diverse physical properties. However, force fields were often developed with relatively limited data sets, yet accuracy for diverse properties across a broad chemical space is desirable; therefore, tests of such accuracy are particularly important. Here, to this end, we calculated 237 infinite dilution activity coefficients (IDACs), comparing with experimental values from NIST's ThermoML database. We found that calculated IDAC values correlate strongly with experiment (Pearson R of 0.92 ± 0.01) and allow us to identify specific functional groups that appear to present challenges to the force field employed. One potentially valuable aspect of IDACs, as compared to solvation free energies, which have been frequently employed as force field tests, is that the same molecules serve both as solutes and solvents in different cases, allowing us to ensure that force fields are not overly tuned to one particular environment or solvent.

Potassium fluoride activation for the nucleophilic fluorination reaction using 18-crown-6, [2.2.2]-cryptand, pentaethylene glycol and comparison with the new hydro-crown scaffold: a theoretical analysis

Ion-pair binding in the hydro-crown ether leads to more effective catalytic cycle for nucleophilic fluorination.