Cluster expansion of the solvation free energy difference: Systematic improvements in the solvation of single ions

Physical and numerical aspects of sodium ion solvation free energies via the cluster-continuum model

Sodium cation solvation Gibbs free energies (ΔG_solv(Na⁺)) have been obtained in water, dimethylformamide, dimethyl sulfoxide, ethanol, acetone, acetonitrile, and methanol through the "monomer cycle" cluster-continuum approach where a solvent reference state is described by infinitely separated molecules. The following steps are vital for obtaining reliable ΔG_solv(Na⁺) values: (a) a meticulous conformational search involving dispersion corrected density functional theory (DFT-D) and the continuum solvation model (CSM); (b) gas-phase DFT-D geometry optimization followed by single-point (SP) domain-based local pair natural orbital coupled clusters including single, double, and partly triple excitation (DLPNO-CCSD(T)) calculations in conjunction with the complete basis set extrapolation; (c) advanced statistical thermodynamic treatment of the low harmonic frequencies (<100 cm^-1) to obtain the robust gas-phase Gibbs free energy correction; (d) gas-phase and dielectric continuum SP with non-electrostatic contributions included in the CSM; (e) an evaluation of the relative thermodynamic stability of the Na⁺(S)_n clusters to identify the number of explicit solvent molecules n to be considered. Our refined computational protocol is promising with a Pearson correlation coefficient between the predicted and experimental data, ρ, of 0.82, and the mean signed and mean unsigned errors of 0.3 and 1.4 kcal mol^-1, respectively.

SN2 versus E2 reactions in a complex microsolvated environment: theoretical analysis of the equilibrium and activation steps of a nucleophilic fluorination

The reactivity of the fluoride ion towards alkyl halides is highly dependent on the solvating environment. In polar aprotic solvents with large counter-ions is highly reactive and produces substantial E2 product, whereas in polar protic solvents leads to slow kinetics and high selectivity for S_N2 reactions. The use of a more complex environment with stoichiometric addition of tert-butanol to acetonitrile solvent is able to module the reactivity and selectivity of tetrabutylammonium fluoride (TBAF). In the present work, we have performed a detailed theoretical analysis of this complex reaction system by density functional theory, continuum solvation model, and including explicit tert-butanol molecules. A kinetic model based on the free energy profile was also used to predict the reactivity and selectivity. The results indicated that the TBAF(tert-butanol) complex plays the key role to increase the S_N2 selectivity, whereas higher aggregates are not relevant. The E2 product is formed exclusively via free TBAF, because the solvating tert-butanol in the TBAF(tert-butanol) complex inhibits the E2 pathway. Our analysis suggests that diols or tetraols could produce an improved selectivity.

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.

Machine Learning-Guided Approach for Studying Solvation Environments

Molecular-level understanding and characterization of solvation environments are often needed across chemistry, biology, and engineering. Toward practical modeling of local solvation effects of any solute in any solvent, we report a static and all-quantum mechanics-based cluster-continuum approach for calculating single-ion solvation free energies. This approach uses a global optimization procedure to identify low-energy molecular clusters with different numbers of explicit solvent molecules and then employs the smooth overlap for atomic positions learning kernel to quantify the similarity between different low-energy solute environments. From these data, we use sketch maps, a nonlinear dimensionality reduction algorithm, to obtain a two-dimensional visual representation of the similarity between solute environments in differently sized microsolvated clusters. After testing this approach on different ions having charges 2+, 1+, 1-, and 2-, we find that the solvation environment around each ion can be seen to usually become more similar in hand with its calculated single-ion solvation free energy. Without needing either dynamics simulations or an a priori knowledge of local solvation structure of the ions, this approach can be used to calculate solvation free energies within 5% of experimental measurements for most cases, and it should be transferable for the study of other systems where dynamics simulations are not easily carried out.