Lithium-coupled electron transfer reactions of nano-confined WOx within Zr-based metal-organic framework

Interfacial charge transfer reactions involving cations and electrons are fundamental to (photo/electro) catalysis, energy storage, and beyond. Lithium-coupled electron transfer (LCET) at the electrode-electrolyte interfaces of lithium-ion batteries (LIBs) is a preeminent example to highlight the importance of charge transfer in modern-day society. The thermodynamics of LCET reactions define the minimal energy for charge/discharge of LIBs, and yet, these parameters are rarely available in the literature. Here, we demonstrate the successful incorporation of tungsten oxides (WOx) within a chemically stable Zr-based metal-organic framework (MOF), MOF-808. Cyclic voltammograms (CVs) of the composite, WOx@MOF-808, in Li+-containing acetonitrile (MeCN)-based electrolytes showed an irreversible, cathodic Faradaic feature that shifted in a Nernstian fashion with respect to the Li+ concentration, i.e., ∼59 mV/log [(Li+)]. The Nernstian dependence established 1:1 stoichiometry of Li+ and e-. Using the standard redox potential of Li+/0, the apparent free energy of lithiation of WOx@MOF-808 (ΔGapp,Li) was calculated to be -36 ± 1 kcal mol-1. ΔGapp,Li is an intrinsic parameter of WOx@MOF-808, and thus by deriving the similar reaction free energies of other metal oxides, their direct comparisons can be achieved. Implications of the reported measurements will be further contrasted to proton-coupled electron transfer (PCET) reactions on metal oxides.

Colorimetric Detection of Sulfur Mustard with 4-(p-Nitrobenzyl)pyridine and Its Derivatives

In this study, we present a simple, highly sensitive, and selective colorimetric method for detecting sulfur mustard (SM) and its simulants. This method relies on a nucleophilic substitution reaction between derivatives of 4-(p-nitrobenzyl)pyridine (NBP) and SM and subsequent treatment with an external base, resulting in a visible response. This reaction exhibits an impressively low detection threshold by the naked eye, as low as 10 ppm at room temperature. In contrast to the conventional use of NBP for detecting other alkylating agents, such as nitrogen mustard, our approach eliminates the need for prolonged heating or intricate extraction processes. Both computational and experimental investigations underscore the significance of water within our detection medium as it stabilizes crucial episulfonium cation intermediates. Furthermore, we demonstrate the practical applicability of this sensor by incorporating it onto cellulose and silica surfaces, which may provide guidance for the design and development of solid-state SM detectors.

The role of charge in microdroplet redox chemistry

In charged water microdroplets, which occur in nature or in the lab upon ultrasonication or in electrospray processes, the thermodynamics for reactive chemistry can be dramatically altered relative to the bulk phase. Here, we provide a theoretical basis for the observation of accelerated chemistry by simulating water droplets of increasing charge imbalance to create redox agents such as hydroxyl and hydrogen radicals and solvated electrons. We compute the hydration enthalpy of OH- and H+ that controls the electron transfer process, and the corresponding changes in vertical ionization energy and vertical electron affinity of the ions, to create OH• and H• reactive species. We find that at ~ 20 - 50% of the Rayleigh limit of droplet charge the hydration enthalpy of both OH- and H+ have decreased by >50 kcal/mol such that electron transfer becomes thermodynamically favorable, in correspondence with the more favorable vertical electron affinity of H+ and the lowered vertical ionization energy of OH-. We provide scaling arguments that show that the nanoscale calculations and conclusions extend to the experimental microdroplet length scale. The relevance of the droplet charge for chemical reactivity is illustrated for the formation of H2O2, and has clear implications for other redox reactions observed to occur with enhanced rates in microdroplets.

Hydration of [Formula: see text]aminobenzoic acid: structures and non-covalent bondings of aminobenzoic acid-water clusters

CONTEXT

Micro-hydration of the aminobenzoic acid is essential to understand its interaction with surrounding water molecules. Understanding the micro-hydration of the aminobenzoic acid is also essential to study its remediation from wastewater. Therefore, we explored the potential energy surfaces (PESs) of the para-aminobenzoic acid-water clusters, ABW[Formula: see text], [Formula: see text], to study the microsolvation of the aminobenzoic acid in water. In addition, we performed a quantum theory of atoms in molecules (QTAIM) analysis to identify the nature of non-covalent bondings in the aminobenzoic acid-water clusters. Furthermore, temperature effects on the stability of the located isomers have been examined. The located structures have been used to calculate the hydration free energy and the hydration enthalpy of the aminobenzoic acid using the cluster continuum solvation model. The hydration free energy and the hydration enthalpy of the aminobenzoic acid at room temperature are evaluated to be -7.0 kcal/mol and -18.1 kcal/mol, respectively. The hydration enthalpy is in perfect agreement with a previous experimental estimate. Besides, temperature effects on the calculated hydration enthalpy and free energy are reported. Finally, we calculated the gas phase binding energies of the most stable structures of the ABW[Formula: see text] clusters using twelve functionals of density functional theory (DFT), including empirical dispersion. The DFT functionals are benchmarked against the DLPNO-CCSD(T)/CBS. We have found that the three most suitable DFT functionals are classified in the following order: PW6B95D3 > MN15 > [Formula: see text]B97XD. Therefore, the PW6B95D3 functional is recommended for further study of the aminobenzoic acid-water clusters and similar systems.

METHODS

The exploration started with classical molecular dynamics simulations followed by complete optimization at the PW6B95D3/def2-TZVP level of theory. Optimizations are performed using Gaussian 16 suite of codes. QTAIM analysis is performed using the AIMAll program.

Experimental Compilation and Computation of Hydration Free Energies for Ionic Solutes

Although charged solutes are common in many chemical systems, traditional solvation models perform poorly in calculating solvation energies of ions. One major obstacle is the scarcity of experimental data for solvated ions. In this study, we release an experiment-based aqueous ionic solvation energy data set, IonSolv-Aq, that contains hydration free energies for 118 anions and 155 cations, more than 2 times larger than the set of hydration free energies for singly charged ions contained in the 2012 Minnesota Solvation Database commonly used in benchmarking studies. We discuss sources of systematic uncertainty in the data set and use the data to examine the accuracy of popular implicit solvation models COSMO-RS and SMD for predicting solvation free energies of singly charged ionic solutes in water. Our results indicate that most SMD and COSMO-RS modeling errors for ionic solutes are systematic and correctable with empirical parameters. We discuss two systematic offsets: one across all ions and one that depends on the functional group of the ionization site. After correcting for these offsets, solvation energies of singly charged ions are predicted using COSMO-RS to 3.1 kcal mol-1 MAE against a challenging test set and 1.7 kcal mol-1 MAE (about 3% relative error) with a filtered test set. The performance of SMD is similar, with MAE against those same test sets of 2.7 and 1.7 kcal mol-1. These results underscore the importance of compiling larger experimental data sets to improve solvation model parametrization and fairly assess performance.

Solv: An Alternative Continuum Model Implementation Based on Fixed Atomic Charges, Scaled Particle Theory, and the Atom-Atom Potential Method

An alternative continuum model implementation is reported. The electrostatic contribution to the solvation Gibbs free energy utilizes the noniterative conductor-like screening model of Vyboishchikov and Voityuk (DOI: 10.1002/jcc.26531) based on the fixed partial atomic charges. The nonelectrostatic solute-solvent dispersion-repulsion energy is calculated through the Caillet-Claverie atom-atom potential method employing the grid-based approach. The nonelectrostatic cavitation energy is computed within the scaled particle theory (SPT) formalism with the solute hard-sphere radius obtained via the Pierotti-Claverie (PC) scheme, from the solute molecular surface (SPT-S) or volume (SPT-V). The solvent hard-sphere radius is derived through the fitting to the experimental total solvation free energies of 2530 neutral species in 92 solvents. Application of the model to reproduce both the absolute and relative (reaction net) solvation free energies indicates that the SPT-V approach based on the CM5 charges is the best performer. The method is suggested for the solvation free energy calculation in the nonaqueous solvents.

Cp* non-innocence and the implications of (η4-Cp*H)Rh intermediates in the hydrogenation of CO2, NAD+, amino-borane, and the Cp* framework - a computational study

In hydrogenation mediated by half-sandwich complexes of Rh, Cp*Rh(III)-H intermediates are critical hydride-delivery agents. For bipyridine-supported complexes, a unique transformation named 'Cp* non-innocence' leads to the appearance of (Cp*H)Rh(I) intermediates, which are purported to exhibit enhanced hydride-delivery capabilities. In this work, DFT calculations performed to compare the role of these complexes in hydrogenation reveal that (Cp*H)Rh(I) intermediates do not compete with the conventional pathway (involving Cp*Rh(III)-H); instead they can lead to sequential hydrogenation of the Cp* framework, and potentially, catalyst degradation. Thus, caution is warranted when invoking the truly homogeneous nature of hydrogenation catalysis mediated by this popular class of complexes.

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.

Predicting Solvent Effects on SN2 Reaction Rates: Comparison of QM/MM, Implicit, and MM Explicit Solvent Models

Solvents are one of the key variables in the optimization of a synthesis yield or properties of a synthesis product. In this paper, contemporary solvent models are applied to predict the rates of S_N2 reactions in a range of aqueous and non-aqueous solvents. High-level CCSD(T)/CBS//M06-2X/6-31+G(d) gas phase energies were combined with solvation free energies from SMD, SM12, and ADF-COSMO-RS continuum solvent models, as well as molecular mechanics (MM) explicit solvent models with different atomic charge schemes to predict the rate constants of three S_N2 reactions in eight protic and aprotic solvents. It is revealed that the prediction of rate constants in organic solvents is not necessarily less challenging than in water and popular solvent models struggle to predict their rate constants to within 3 log units of experimental values. Among the continuum solvent models, the ADF-COSMO-RS model performed the best in predicting absolute rate contants while the SM12 model was best at predicting relative rate constants with an average accuracy of about 1.5 and 0.8 log units, respectively. The use of computationally more demanding MM explicit solvent models did not translate to improvements in absolute rate constants but was quite effective at predicting relative rate constants due to systematic error cancellation. Free energy barriers obtained from umbrella sampling with explicit solvent QM/MM simulations led to excellent agreement with experimental values, provided that a validated level of theory is used to treat the QM region.

Universal Trends between Acid Dissociation Constants in Protic and Aprotic Solvents

pKa values in non-aqueous solvents are of critical importance in many areas of chemistry. Our knowledge is, despite their relevance, still limited to the most fundamental properties and few pKa values in the most common solvents. Taking advantage of a recently introduced computationally efficient procedure we computed the pKa values of 182 compounds in 21 solvents. This data set is used to establish for the first time universal trends between all solvents. Our computations indicate, that the total charge of the molecule and the charge of the acidic group combined with the Kamlet-Taft solvatochromic parameters are sufficient to predict pKa values with at least semi- quantitative accuracy. We find, that neutral acids such as alcohols are strongly affected by the solvent properties. This is contrasted by cationic acids like ammonium ions whose pKa is often almost completely independent from the choice of solvent.

Modeling pKa of the Brønsted Bases as an Approach to the Gibbs Energy of the Proton in Acetonitrile

A simple but efficient computational approach to calculate pK_a in acetonitrile for a set of phosphorus, nitrogen, and carbon bases was established. A linear function that describes relations between the calculated ΔG’_a.sol(BH⁺) and pK_a values was determined for each group of bases. The best model was obtained through the variations in the basis set, in the level of theory (density functionals or MP2), and in the continuum solvation model (IPCM, CPCM, or SMD). The combination of the IPCM/B3LYP/6-311+G(d,p) solvation approach with MP2/6-311+G(2df,p)//B3LYP/6-31G(d) gas-phase energies provided very good results for all three groups of bases with R² values close to or above 0.99. Interestingly, the slopes and the intercepts of the obtained linear functions showed significant deviations from the theoretical values. We made a linear plot utilizing all the conducted calculations and all the structural variations and employed methods to prove the systematic nature of the intercept/slope dependence. The interpolation of the intercept to the ideal slope value enabled us to determine the Gibbs energy of the proton in acetonitrile, which amounted to −258.8 kcal mol⁻¹. The obtained value was in excellent agreement with previously published results.

Conformational energies of microsolvated Na+ clusters with protic and aprotic solvents from GFNn-xTB methods

Performance of contemporary tight-binding semiempirical GFNn-xTB methods for the conformational energies of singly charged sodium clusters Na⁺ (S)_n (n = 4-8) with 3 protic and 8 aprotic solvents is examined against the reference RI-MP2/CBS method. The median Pearson correlation coefficients of ρ = 0.84 (GFN2-xTB) and ρ = 0.82 (GFN1-xTB) do not give the clear preference to any tested approach. GFN1-xTB method demonstrates more stable performance than its GFN2-xTB successor with the average mean absolute errors (MAEs)/mean signed errors (MSEs) of 1.2/0.2 and 2.3/1.6 kcal mol^-1 , respectively. Conformational energies produced by the computationally efficient DFT functional PBE and double-ζ basis set complemented with -D3(BJ) dispersion correction are suitable for the preliminary sampling (median ρ = 0.93), but should be used with a caution for the calculations of the average ensemble properties (MAE/MSE = 1.7/1.1 kcal mol^-1 ). Higher-ranking PBE0-D3(BJ) and ωB97M-V with triple-ζ basis sets yield significantly lower MAEs/MSEs of 0.55/0.20 and 0.51/0.23 kcal mol^-1 , respectively.

How to Predict the pK a of Any Compound in Any Solvent

Acid-base properties of molecules in nonaqueous solvents are of critical importance for almost all areas of chemistry. Despite this very high relevance, our knowledge is still mostly limited to the pK _a of rather few compounds in the most common solvents, and a simple yet truly general computational procedure to predict pK _a's of any compound in any solvent is still missing. In this contribution, we describe such a procedure. Our method requires only the experimental pK _a of a reference compound in water and a few standard quantum-chemical calculations. This method is tested through computing the proton solvation energy in 39 solvents and by comparing the pK _a of 142 simple compounds in 12 solvents. Our computations indicate that the method to compute the proton solvation energy is robust with respect to the detailed computational setup and the construction of the solvation model. The unscaled pK _a's computed using an implicit solvation model on the other hand differ significantly from the experimental data. These differences are partly associated with the poor quality of the experimental data and the well-known shortcomings of implicit solvation models. General linear scaling relationships to correct this error are suggested for protic and aprotic media. Using these relationships, the deviations between experiment and computations drop to a level comparable to that observed in water, which highlights the efficiency of our method.

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.

Can Anions Be Inserted into MXene?

Despite the continuous progress in the research and development of Ti₃C₂T_x (MXene) electrodes for high-power batteries and supercapacitor applications, the role of the anions in the electrochemical energy storage and their ability to intercalate between the MXene sheets upon application of positive voltage have not been clarified. A decade after the discovery of MXenes, the information about the possibility of anion insertion into the restacked MXene electrode is still being questioned. Since the positive potential stability range in diluted aqueous electrolytes is severely limited by anodic oxidation of the Ti, the possibility of anion insertion was evaluated in concentrated aqueous electrolyte solutions and aprotic electrolytes as well. To address this issue, we have conducted in situ gravimetric electrochemical quartz crystal microbalance with dissipation monitoring (EQCM-D) measurements in highly concentrated LiCl and LiBr electrolytes, which enable a significant extension of the operation range of the MXene electrodes toward positive potentials. Also, halogens are among the smallest anions and should be easier to intercalate between MXene layers, in comparison to multiatomic anions. On the basis of mass change variations in the positive voltage range and complementary density functional theory calculations, it was demonstrated that insertion of anionic species into MXene, within the range of potentials of interest for capacitive energy storage, is not likely to occur. This can be explained by the strong negative charge on Ti₃C₂T_x sheets terminated by functional groups.

Uncovering Differences in Hydration Free Energies and Structures for Model Compound Mimics of Charged Side Chains of Amino Acids

Free energies of hydration are of fundamental interest for modeling and understanding conformational and phase equilibria of macromolecular solutes in aqueous phases. Of particular relevance to systems such as intrinsically disordered proteins are the free energies of hydration and hydration structures of model compounds that mimic charged side chains of Arg, Lys, Asp, and Glu. Here, we deploy a Thermodynamic Cycle-based Proton Dissociation (TCPD) approach in conjunction with data from direct measurements to obtain estimates for the free energies of hydration for model compounds that mimic the side chains of Arg+, Lys+, Asp-, and Glu-. Irrespective of the choice made for the hydration free energy of the proton, the TCPD approach reveals clear trends regarding the free energies of hydration for Arg+, Lys+, Asp-, and Glu-. These trends include asymmetries between the hydration free energies of acidic (Asp- and Glu-) and basic (Arg+ and Lys+) residues. Further, the TCPD analysis, which relies on a combination of experimental data, shows that the free energy of hydration of Arg+ is less favorable than that of Lys+. We sought a physical explanation for the TCPD-derived trends in free energies of hydration. To this end, we performed temperature-dependent calculations of free energies of hydration and analyzed hydration structures from simulations that use the polarizable Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field and water model. At 298 K, the AMOEBA model generates estimates of free energies of hydration that are consistent with TCPD values with a free energy of hydration for the proton of ca. -259 kcal/mol. Analysis of temperature-dependent simulations leads to a structural explanation for the observed differences in free energies of hydration of ionizable residues and reveals that the heat capacity of hydration is positive for Arg+ and Lys+ and negative for Asp- and Glu-.

The role of intermolecular forces in ionic reactions: the solvent effect, ion-pairing, aggregates and structured environment

The environment enclosing an ionic species has a critical effect on its reactivity. In a more general sense, medium effects are not limited to the solvent, but involve the counter ion effect (ion pairing), formation of larger aggregates and structured environment as provided by the host in the case of host-guest complexes. In this review, a general view of the medium effect on anion-molecule reactions is presented. Nucleophilic substitution reactions of aliphatic (SN2) and aromatic (SNAr) systems, as well as elimination reactions (E2), are the focus of the discussion. In particular, nucleophilic fluorination with KF, CsF and tetraalkylammonium fluoride was used as the main model, because of the importance of this kind of reaction and the recent advances in the study of these systems. The solvent effect, ion pairing, formation of aggregates and formation of complexes with crown ethers, cryptands and calixarenes are discussed. For a deeper insight into the medium effect, many results of reliable theoretical calculations in close agreement with experiments were chosen as examples.

Predictive QM/MM Modeling of Modulations in Protein-Protein Binding by Lysine Methylation

Lysine methylation is a key regulator of protein-protein binding. The amine group of lysine can accept up to three methyl groups, and experiments show that protein-protein binding free energies are sensitive to the extent of methylation. These sensitivities have been rationalized in terms of chemical and structural features present in the binding pockets of methyllysine binding domains. However, understanding their specific roles requires an energetic analysis. Here we propose a theoretical framework to combine quantum and molecular mechanics methods, and compute the effect of methylation on protein-protein binding free energies. The advantages of this approach are that it derives contributions from all local non-trivial effects of methylation on induction, polarizability and dispersion directly from self-consistent electron densities, and at the same time determines contributions from well-characterized hydration effects using a computationally efficient classical mean field method. Limitations of the approach are discussed, and we note that predicted free energies of fourteen out of the sixteen cases agree with experiment. Critical assessment of these cases leads to the following overarching principles that drive methylation-state recognition by protein domains. Methylation typically reduces the pairwise interaction between proteins. This biases binding toward lower methylated states. Simultaneously, however, methylation also makes it easier to partially dehydrate proteins and place them in protein-protein complexes. This latter effect biases binding in favor of higher methylated states. The overall effect of methylation on protein-protein binding depends ultimately on the balance between these two effects, which is observed to be tuned via several combinations of local features.

Reversible PtII-CH3 deuteration without methane loss: metal-ligand cooperation vs. ligand-assisted PtII-protonation

Di(2-pyridyl)ketone dimethylplatinum(ii), (dpk)Pt^II(CH₃)₂, reacts with CD₃OD at 25 °C to undergo complete deuteration of Pt-CH₃ fragments in ∼5 h without loss of methane to form (dpk)Pt^II(CD₃)₂ in virtually quantitative yield. The deuteration can be reversed by dissolution in CH₃OH or CD₃OH. Kinetic analysis and isotope effects, together with support from density functional theory calculations indicate a metal-ligand cooperative mechanism wherein DPK enables Pt-CH₃ deuteration by allowing non-rate-limiting protonation of Pt^II by CD₃OD. In contrast, other model di(2-pyridyl) ligands enable rate-limiting protonation of Pt^II, resulting in non-rate-limiting C-H(D) reductive coupling. Owing to its electron-poor nature, following complete deuteration, DPK can be dissociated from the Pt^II-centre, furnishing [(CD₃)₂Pt^II(μ-SMe₂)]₂ as the perdeutero analogue of [(CH₃)₂Pt^II(μ-SMe₂)]₂, a commonly used Pt^II-precursor.

Ambiguities in solvation free energies from cluster-continuum quasichemical theory: lithium cation in protic and aprotic solvents

Gibbs free energies for Li+ solvation in water, methanol, acetonitrile, DMSO, dimethylacetamide, dimethoxyethane, dimethylformamide, gamma-butyrolactone, pyridine, and sulfolane have been calculated using the cluster-continuum quasichemical theory. With n independent solvent molecules S initial state forming the "monomer" thermodynamic cycle, Li+ solvation free energies are found to be on average 14 kcal mol-1 more positive compared to those from the "cluster" thermodynamic cycle where the initial state is the cluster Sn. We ascribe the inconsistency between the "monomer" and "cluster" cycles mainly to the incorrectly predicted solvation free energies of solvent clusters Sn from the SMD and CPCM continuum solvation models, which is in line with the earlier study of Bryantsev et al., J. Phys. Chem. B, 2008, 112, 9709-9719. When experimental-based solvation free energies of individual solvent molecules and solvent clusters are employed, the "monomer" and "cluster" cycles result in identical numbers. The best overall agreement with experimental-based "bulk" scale lithium cation solvation free energies was obtained for the "monomer" scale, and we recommend this set of values. We expect that further progress in the field is possible if (i) consensus on the accuracy of experimental reference values is achieved; (ii) the most recent continuum solvation models are properly parameterized for all solute-solvent combinations and become widely accessible for testing.

Absolute ion hydration free energy scale and the surface potential of water via quantum simulation

With a goal of determining an absolute free energy scale for ion hydration, quasi-chemical theory and ab initio quantum mechanical simulations are employed to obtain an accurate value for the bulk hydration free energy of the Na⁺ ion. The free energy is partitioned into three parts: 1) the inner-shell or chemical contribution that includes direct interactions of the ion with nearby waters, 2) the packing free energy that is the work to produce a cavity of size λ in water, and 3) the long-range contribution that involves all interactions outside the inner shell. The interfacial potential contribution to the free energy resides in the long-range term. By averaging cation and anion data for that contribution, cumulant terms of all odd orders in the electrostatic potential are removed. The computed total is then the bulk hydration free energy. Comparison with the experimentally derived real hydration free energy produces an effective surface potential of water in the range -0.4 to -0.5 V. The result is consistent with a variety of experiments concerning acid-base chemistry, ion distributions near hydrophobic interfaces, and electric fields near the surface of water droplets.

Prediction of Solvation Free Energies of Ionic Solutes in Neutral Solvents

The prediction of solvation free energies is essential for a variety of applications. Solvation free energies of neutral systems can be predicted quite accurately. The accuracy of predictions for solvation free energies of ionic solutes dissolved in neutral solvents, however, has been reported to be worse by at least 1 order of magnitude. In this study, the performance of three approaches for solvation free energy prediction of several hundred ions dissolved in neutral solvents is evaluated. The applied methods are COSMO-RS, cluster continuum model (CCM) together with COSMO-RS, and COSMO-RS-ES. It is emphasized that the reference data for model evaluation are subject to large uncertainties stemming from the impossibility to measure the so-called elusive absolute free energies of solvation of a single ion. Consequently, such uncertainty must be considered during the evaluation of prediction methods. Therefore, a straightforward approach to account for the underlying uncertainty is applied here. Hereby, it is revealed that the true performance of the method is better than what is often reported. The average absolute deviation (AAD) of COSMO-RS is calculated to be 2.3 kcal mol^-1, while applying the CCM and COSMO-RS-ES each results in AADs of 2.0 kcal mol^-1. This accuracy allows for qualitative assessment of solvation free energy-dependent quantities, such as reaction rate constants.

Solvation energies of ions with ensemble cluster-continuum approach

An alternative cluster-continuum approach for the calculation of solvation free energies of ions.

Ion hydration free energies and water surface potential in water nano drops: The cluster pair approximation and the proton hydration Gibbs free energy in solution

We estimate both single ion hydration Gibbs free energies in water droplets, comprising from 50 to 1000 molecules, and water/vacuum surface potentials in pure water droplets comprising up to 10 000 molecules. We consider four ions, namely, Li⁺, NH₄ ⁺, F^-, and Cl^-, and we model their hydration process and water/water interactions using polarizable force fields based on an induced point dipole approach. We show both ion hydration Gibbs free energies and water surface potentials to obey linear functions of the droplet radius as soon as droplets comprising a few hundred water molecules. Moreover, we also show that the differences in anion/cation hydration Gibbs free energies in droplets obey a different regime in large droplets than in small clusters comprising no more than six water molecules, in line with the earlier results computed from standard additive point charge force fields. Hence, both point charge and more sophisticated induced point dipole molecular modeling approaches suggest that methods considering only the thermodynamical properties of small ion/water clusters to estimate the absolute proton hydration Gibbs free energy in solution are questionable. In particular, taking into account the data of large ion/water droplets may yield a proton hydration Gibbs free energy in solution value to be shifted by several k_BT units compared to small clusters-based approaches.

Effect of Cation Size on Solvation and Association with Superoxide Anion in Aprotic Solvents

Influence of cation size on solvation strength, diffusion, and kinetics of the association reaction with anions O₂ ^- in aprotic solvents, such as acetonitrile and dimethyl sulfoxide, has been investigated by means of molecular dynamics simulations. The work is motivated by the need to understand the molecular nature of the solvent-induced changes in capacity of Li-air batteries. We have shown that the dependence of the solvation shell stability on the cation size has a maximum at a particular ion radius that corresponds to a solvent coordination number of 4. The shell stability maximum coincides with the diffusion coefficient minimum. The variation of the cation shell stability has a crucial impact on the kinetics of the cation-O₂ ^- association. We have demonstrated that profound inhibition of the association reaction for Li⁺ in dimethyl sulfoxide is a result of the lock-and-key effect that cannot be described in the framework of Hard Soft Acid Base theory.