Theoretical Insights into Specific Ion Effects and Strong-Weak Acid-Base Rules for Ions in Solution: Deriving the Law of Matching Solvent Affinities from First Principles

Physicochemical and sweetness behavior of kosmotropic and chaotropic ions in aqueous solutions of polyhydroxy compound

Hydrated ions have essential applications in biochemical processes. The objective of the present work is to study the behaviour of kosmotropic and chaotropic ions in aqueous solutions of polyhydroxy compounds in terms of sweetness response, volumetric and acoustic parameters such as apparent molar volume (∅v), apparent specific volume (ASV), compressibility (∅k0) etc. at different temperatures 20 °C - 45 °C and 101 kPa. The positive ∅v values of arabitol showed more kosmotropic nature for Na+ ion developing more stabilize structure of arabitol with water through strong hydrogen bonding, while the presence of K+ and Cl- form weak hydrogen bonding with water termed as chaotropic ions. The ASV values showed that arabitol remained sweet in the presence of electrolytes. The negative results of ∅k0 for arabitol in aqueous electrolytic solutions give valuable information about interactions occurring among solute, electrolytic ions, and water molecules. The structure-making/breaking behaviour has been observed by Hepler's constant (∂2∅vo/∂T2)p.

Insights into Hildebrand Solubility Parameters - Contributions from Cohesive Energies or Electrophilicity Densities?

We introduce certain concepts and expressions from conceptual density functional theory (DFT) to study the properties of the Hildebrand solubility parameter. The original form of the Hildebrand solubility parameter is used to qualitatively estimate solubilities for various apolar and aprotic substances and solvents and is based on the square root of the cohesive energy density. Our results show that a revised expression allows the replacement of cohesive energy densities by electrophilicity densities, which are numerically accessible by simple DFT calculations. As an extension, the reformulated expression provides a deeper interpretation of the main contributions and, in particular, emphasizes the importance of charge transfer mechanisms. All calculated values of the Hildebrand parameters for a large number of common solvents are compared with experimental values and show good agreement for non- or moderately polar aprotic solvents in agreement with the original formulation of the Hildebrand solubility parameters. The observed deviations for more polar and protic solvents define robust limits from the original formulation which remain valid. Likewise, we show that the use of machine learning methods leads to only slightly better predictability.

Solvation of Nanoions in Aqueous Solutions

In recent years it has been increasingly recognized that different classes of large ions with multiple valency have effects conceptually similar to weakly solvated ions in the Hofmeister series, also labeled by the term chaotropic. The term "superchaotropic effect" has been coined because these effects are much more strongly pronounced for nanometer-sized ions, whose adsorption properties often resemble typical surfactants. Despite this growing interest in these nanometer-sized ions, a simple conceptual extension of the Hofmeister series toward nanoions has not been achieved because an extrapolation of the one-dimensional surface charge density scale does not lead to the superchaotropic regime. In this work, we discuss a generic model that is broadly applicable to ions of nearly spherical shape and thus includes polyoxometalates and boron clusters. We present a qualitative classification scheme in which the ion size appears as a second dimension. Ions of different sizes but the same charge density differ in their bulk solvation free energy. As the ions grow bigger at constant surface charge density, they become more stable in solution, but the adsorption behavior is still governed by the surface charge density. A detailed molecular dynamics simulation study of large ions that is based on a shifted Lennard-Jones potential is presented that supports the presented classification scheme.

Approximating constant potential DFT with canonical DFT and electrostatic corrections

The complexity of electrochemical interfaces has led to the development of several approximate density functional theory (DFT)-based schemes to study reaction thermodynamics and kinetics as a function of electrode potential. While fixed electrode potential conditions can be simulated with grand canonical ensemble DFT (GCE-DFT), various electrostatic corrections on canonical, constant charge DFT are often applied instead. In this work, we present a systematic derivation and analysis of the different electrostatic corrections on canonical DFT to understand their physical validity, implicit assumptions, and scope of applicability. Our work highlights the need to carefully address the suitability of a given model for the problem under study, especially if physical or chemical insight in addition to reaction energetics is sought. In particular, we analytically show that the different corrections cannot differentiate between electrostatic interactions and covalent or charge-transfer interactions. By numerically testing different models for CO₂ adsorption on a single-atom catalyst as a function of the electrode potential, we further show that computed capacitances, dipole moments, and the obtained physical insight depend sensitively on the chosen approximation. These features limit the scope, generality, and physical insight of these corrective schemes despite their proven practicality for specific systems and energetics. Finally, we suggest guidelines for choosing different electrostatic corrections and propose the use of conceptual DFT to develop more general approximations for electrochemical interfaces and reactions using canonical DFT.

Quantitative Solvation Energies from Gas-Phase Calculations: First-Principles Charge Transfer and Perturbation Approaches

We present a first-principles approach for the calculation of solvation energies and enthalpies with respect to different ion pair combinations in various solvents. The method relies on the conceptual density functional theory (DFT) of solvation, from which detailed expressions for the solvation energies can be derived. In addition to fast and straightforward gas phase calculations, we also study the influence of modified chemical reactivity descriptors in terms of electronic perturbations. The corresponding phenomenological changes in molecular energy levels can be interpreted as the influence of continuum solvents. Our approach shows that the introduction of these modified expressions is essential for a quantitative agreement between the calculated and the experimental results.

Specific Ion and Electric Field Controlled Diverse Ion Distribution and Electroosmotic Transport in a Polyelectrolyte Brush Grafted Nanochannel

Controlling ion distribution inside a charged nanochannel is central to using such channels in diverse applications. Here, we show the possibility of using a charged polyelectrolyte (PE) brush-grafted nanochannel for triggering diverse nanoscopic ion distribution and nanofluidic electroosmotic transport by controlling the valence and size of the counterions (that screen the charges of the PE brushes) and the strength of an externally applied axial electric field. We atomistically simulate separate cases of fully charged polyacrylic acid (PAA) brush functionalized nanochannels with Na⁺, Cs⁺, Ca²⁺, Ba²⁺, and Y³⁺ counterions screening the PE charges. Four key findings emerge from our simulations. First, we find that the counterions with a greater valence and a smaller size prefer to remain localized inside the brush layer. Second, for the case where there is an added chloride salt with the same cation (as the screening counterions), there are more coions (Cl^- ions) in the brush-free bulk than counterions (for counterions Na⁺, Ca²⁺, Ba²⁺, Y³⁺): this is a manifestation of the overscreening (OS) of the PE brush layer. Contrastingly, the number of Cs⁺ ions remain higher than the Cl^- ions inside the brush-free bulk, ensuring that there is no OS effect for this case. Third, large applied electric field enables a few Na⁺, Cs⁺, and Ba²⁺ counterions to leave the brush layer and to go to the bulk: this makes the OS of the PE brush layer disappear for the cases of PE brushes being screened by the Na⁺ and Ba²⁺ ions. On the other hand, no such electric-field-mediated disappearance of OS is observed for the cases of Ca²⁺ and Y³⁺ screening counterions; we attribute this to the firm attachment of these counterions to the negatively charged monomers. Free energy associated with a counterion binding to a PE chain corroborates this diversity in the counterion-specific response to the applied electric field. Finally, we demonstrate that such diverse ion distributions, along with specific electric-field-strength-dependent ion properties, lead to (1) electroosmotic (EOS) transport in nanochannels grafted with PAA brushes screened with Cs⁺ ions to be always counterion dominated, (2) EOS transport in nanochannels grafted with PAA brushes screened with Ca²⁺ and Y³⁺ ions to be always coion-dominated, and (3) EOS transport in nanochannels grafted with PAA brushes screened with Na⁺ and Ba²⁺ ions to be coion dominated for smaller electric fields and counterion dominated for larger electric fields.

Application of Fundamental Chemical Principles for Solvation Effects: A Unified Perspective for Interaction Patterns in Solution

We demonstrate the utility of basic chemical principles like the "|Δμ| big is good" (DMB) rule for the study of solvation interactions between distinct solutes such as ions and solvents. The corresponding approach allows us to define relevant criteria for maximum solvation energies of ion pairs in different solvents in terms of electronegativities and chemical hardnesses. Our findings reveal that the DMB principle culminates into the strong and weak acids and bases concept as recently derived for specific ion effects in various solvents. The further application of the DMB approach highlights a similar condition for the chemical hardnesses with a reminiscence to the hard/soft acids and bases principle. Comparable conclusions can also be drawn with regard to the change of the solvent. We show that favorable solvent interactions are mainly driven by low chemical hardnesses as well as high electronegativity differences between the ions and the solvent. Our findings highlight that solvation interactions are governed by basic chemical principles, which demonstrates the close similarity between solvation mechanisms and chemical reactions.

Electronic properties of amino acids and nucleobases: similarity classes and pairing principles from chemical reactivity indices

We present a new classification scheme for amino acids and nucleobases based on the electronic properties of the individual molecules. Using chemical reactivity indices such as electronegativity, electrophilicity, and chemical hardness, we can identify similarities and differences between each class of amino acids and nucleobases. Notable differences emerge in particular with regard to high, neutral or low electronegativity as well as different combinations of chemical hardness. Our approach allows us to relate these insights to the properties of the side groups in terms of a unique reference scheme. We further show that hydrophobic differences between amino acids are rather negligible in the context of electronic properties. Our classification scheme also rationalizes the occurrence of distinct stable nucleobase pairs and clearly emphasizes certain differences between individual molecules. The stability and abundant occurrence of Watson-Crick nucleobase pairs is further discussed in the context of the minimum electrophilicity principle.

Probing the distribution of ionic liquid mixtures at charged and neutral interfaces via simulations and lattice-gas theory

Room temperature ionic liquid solutions confined between neutral and charged surfaces are investigated by means of atomistic Molecular Dynamics simulations. We study 1-ethyl-3-methylimidazolium dicyanamide ([EMIm]⁺[DCA]^-) in water or dimethylsulfoxide (DMSO) mixtures in confinement between two interfaces. The analysis is based on the comparison of the molecular species involved and the charged state of the surfaces. Focus is given on the influence of different water/DMSO concentrations on the microstructuring and accumulation of each species. Thermodynamic aspects, such as the entropic contributions in the observed trends are obtained from the simulations using a lattice-gas theory. The results clearly underline the differences in these properties for the water and DMSO mixtures and unravel the underlying mechanisms and inherent details. We were able to pinpoint the importance of the size and the relative permittivity of the molecules in guiding their microstructuring in the vicinity of the surfaces, as well as their interactions with the latter, i.e. the solute-surface interactions. The influence of water and DMSO on the overscreening at charged interfaces is also discussed. The analysis on the molecular accumulation at the interfaces allows us to predict whether the accumulation is entropy or enthalpy driven, which has an impact in the removal of the molecular species from the surfaces. We discuss the impact of this work in providing an essential understanding towards a careful design of electrochemical elements.

Molecular Interactions From the Density Functional Theory for Chemical Reactivity: The Interaction Energy Between Two-Reagents

Reactivity descriptors indicate where a reagent is most reactive and how it is most likely to react. However, a reaction will only occur when the reagent encounters a suitable reaction partner. Determining whether a pair of reagents is well-matched requires developing reactivity rules that depend on both reagents. This can be achieved using the expression for the minimum-interaction-energy obtained from the density functional reactivity theory. Different terms in this expression will be dominant in different circumstances; depending on which terms control the reactivity, different reactivity indicators will be preferred.

Understanding specific ion effects and the Hofmeister series

Specific ion effects (SIE), encompassing the Hofmeister Series, have been known for more than 130 years since Hofmeister and Lewith's foundational work. SIEs are ubiquitous and are observed across the medical, biological, chemical and industrial sciences. Nevertheless, no general predictive theory has yet been able to explain ion specificity across these fields; it remains impossible to predict when, how, and to what magnitude, a SIE will be observed. In part, this is due to the complexity of real systems in which ions, counterions, solvents and cosolutes all play varying roles, which give rise to anomalies and reversals in anticipated SIEs. Herein we review the historical explanations for SIE in water and the key ion properties that have been attributed to them. Systems where the Hofmeister series is perturbed or reversed are explored, as is the behaviour of ions at the liquid-vapour interface. We discuss SIEs in mixed electrolytes, nonaqueous solvents, and in highly concentrated electrolyte solutions - exciting frontiers in this field with particular relevance to biological and electrochemical applications. We conclude the perspective by summarising the challenges and opportunities facing this SIE research that highlight potential pathways towards a general predictive theory of SIE.

Enhancement of thermal stability of proteinase K by biocompatible cholinium-based ionic liquids

Proteinase K (PK) is a proteolytic enzyme that has been widely used in nucleic acid purification, leather production, environmental protection, and other industrial applications. However, this biocatalyst cannot tolerate high temperatures which has severely restricted its wider application. As reported in previous studies, cholinium-based ionic liquids (ILs) have gained tremendous attention serving as a promising media to stabilize and preserve proteins, DNA, and other biomolecules due to their environmentally benign nature and biocompatibility. In this work, we chose 13 different kinds of cholinium-based ILs to examine their effects on the thermal stability and enzymatic activity of PK. We found that biocompatible cholinium-based ions with appropriately chosen anions can greatly improve the thermal stability of PK, whose melting temperature (T_m) is increased from ∼74.4 °C to 87.7 °C. However, the enzymatic activity is slightly reduced in the presence of ILs. Further comparison of our results with other literature findings suggests that kosmotropic anions of cholinium-based ILs are crucial to maintain the thermal stability of proteins. However, to achieve the best performance, the choice of IL anions is protein specific.

Specific Ion Effects in Different Media: Current Status and Future Challenges

We discuss the current state of research as well as the future challenges for a deeper understanding of specific ion effects in protic and aprotic solvents as well as various additional media. Despite recent interest in solute or interfacial effects, we focus exclusively on the specific properties of ions in bulk electrolyte solutions. Corresponding results show that many mechanisms remain unknown for these simple media, although theoretical, computational, and experimental studies have provided some insights into explaining individual observations. In particular, the importance of local interactions and electronic properties is emphasized, which enabled a more consistent interpretation of specific ion effects over the past years. Despite current insufficient knowledge, we also discuss future challenges in relation to dynamic properties as well as the influence of different concentrations, different solvents, and solute contributions to gain a deeper understanding of specific ion effects for technological applications.

Specific Ion Solvation and Pairing Effects in Glycerol Carbonate

Identifying the driving forces behind the solvation of inorganic salts by nonaqueous solvents is an important step in the development of green solvents. Here we focus on one promising solvent: glycerol carbonate (GC). Using ab initio molecular dynamics simulations, we build upon our previous work by detailing glycerol carbonate's interactions with a series of anions, a lithium ion, and the LiF ion pair. Through these investigations, we highlight the changes in solvation behavior as the anion size increases, the competition of binding shown by lithium for the oxygens of GC, and the behavior of the LiF ion pair in a GC solution. These results indicate the importance of the cation's identity in ion-pairing structure and dynamics and lend insight into the key factors behind the specific ion effects seen in GC.

Electronic Properties of Protein Destabilizers and Stabilizers: Implications for Preferential Binding and Exclusion Mechanisms

We study the electronic properties of low-weight organic co-solutes by means of conceptual density functional theory calculations. Our results highlight the important role of certain chemical reactivity descriptors such as chemical hardness, electronegativity, nucleofugality, and the electrofugality as important criteria to classify protein stabilizers and destabilizers. Our results imply Lewis basic properties with lower chemical hardness for stabilizers, while destabilizers show higher Lewis acidity with higher chemical hardness. Further consideration of analytical calculations in terms of transfer energies reveals the crucial role of co-solute-protein interactions which significantly change the interaction pattern of the stabilizing or destabilizing species. The corresponding outcomes connect statistical thermodynamics with the electronic properties of co-solutes and also allow us to define general principles for strong stabilizers and destabilizers.

Stabilization of DPPC lipid bilayers in the presence of co-solutes: molecular mechanisms and interaction patterns

We study the interactions between dipalmitoylphosphatidylcholine (DPPC) lipid bilayers in the gel and the fluid phase with ectoine, amino ectoine and water molecules by means of atomistic molecular dynamics (MD) simulations and conceptual density functional theory (DFT) calculations. Our results reveal a pronounced preferential exclusion of both co-solutes from the DPPC lipid bilayer which is stronger for the fluid phase. The corresponding outcomes can be brought into relation with the Kirkwood-Buff theory of solutions in order to provide a thermodynamic rationale for the experimentally observed stabilization of the gel phase. Closely related to preferential exclusion of both co-solutes, our simulations also highlight a preferential hydration behavior as manifested by an increased number of hydrogen bonds between water and DPPC molecules. All results are rationalized by conceptual DFT calculations with regard to differences in the electronic properties between ectoine and amino ectoine.

Long-Range Electrostatic Colloidal Interactions and Specific Ion Effects in Deep Eutectic Solvents

While the traditional consensus dictates that high ion concentrations lead to negligible long-range electrostatic interactions, we demonstrate that electrostatic correlations prevail in deep eutectic solvents where intrinsic ion concentrations often surpass 2.5 M. Here we present an investigation of intermicellar interactions in 1:2 choline chloride:glycerol and 1:2 choline bromide:glycerol using small-angle neutron scattering. Our results show that long-range electrostatic repulsions between charged colloidal particles occur in these solvents. Interestingly, micelle morphology and electrostatic interactions are modulated by specific counterion condensation at the micelle interface despite the exceedingly high concentration of the native halide from the solvent. This modulation follows the trends described by the Hofmeister series for specific ion effects. The results are rationalized in terms of predominant ion-ion correlations, which explain the reduction in the effective ionic strength of the continuum and the observed specific ion effects.

Extending conceptual DFT to include additional variables: oriented external electric field

The extension of the E = E[N, v] functional for exploring chemical reactivity in a conceptual DFT context to include external electric fields is discussed. Concentrating on the case of a homogeneous field the corresponding response functions are identified and integrated, together with the conventional response functions such as permanent dipole moment and polarizability, in an extended response function tree associated with the E = E[N, v, ε] functional. In a case study on the dihalogens F2, Cl2, Br2, I2 the sensitivity of condensed atomic charges (∂q/∂ε) is linked to the polarizability of the halogen atoms. The non-integrated (∂ρ(r)/∂ε) response function, directly related to the field induced density change, is at the basis of these features. It reveals symmetry breaking for a perpendicular field, not detectable in its atom condensed counterpart, and accounts for the induced dipole moment directly related to the molecular polarizability. The much higher sensitivity of the electronic chemical potential/electronegativity as compared to the chemical hardness is highlighted. The response of the condensed Fukui functions to a parallel electric field increases when going down in the periodic table and is interpreted in terms of the extension of the outer contours in the non-condensed Fukui function. In the case of a perpendicular field the (∂f(r)/∂ε) response function hints at stereoselectivity with a preferential side of attack which is not retrieved in its condensed form. In an application the nucleophilic attack on the carbonyl group in H2CO is discussed. Similar to the dihalogens, stereoselectivity is displayed in the Fukui function for nucleophilic attack (f+) in the case of a perpendicular electric field, and opposite to the one that would arise based on the induced density. Disentangling the expression for the evolution of the Fukui function in the presence of an electric field reveals that this difference can be traced back to local differences in the polarization or induced density between the anionic and the neutral system. This difference may be exploited, e.g. for an appropriately substituted H2CO, to generate enantioselectivity.

Artificial neural networks for the prediction of solvation energies based on experimental and computational data

The knowledge of thermodynamic properties for novel electrolyte formulations is of fundamental interest for industrial applications as well as academic research. Herewith, we present an artificial neural networks (ANN) approach for the prediction of solvation energies and entropies for distinct ion pairs in various protic and aprotic solvents. The considered feed-forward ANN is trained either by experimental data or computational results from conceptual density functional theory calculations. The proposed concept of mapping computed values to experimental data lowers the amount of time-consuming and costly experiments and helps to overcome certain limitations. Our findings reveal high correlation coefficients between predicted and experimental values which demonstrate the validity of our approach.