Exploring Ion-Ion Interactions in Aqueous Solutions by a Combination of Molecular Dynamics and Neutron Scattering

On the Elusive Links between Solution Microstructure, Dynamics, and Solvation Thermodynamics: Demystifying the Path through a Bridge over Troubled Conjectures and Misinterpretations

We build a fundamentally based bridge between the solute-induced microstructural perturbation of the species environment and the dynamic as well as thermodynamic responses of the fluid system, regardless of the state conditions, composition, nature of the solvent, and either the magnitude or the type of solute-solvent intermolecular-interaction asymmetries. For that purpose, we advance a fluctuation-based solvation formalism of fluid mixtures to provide meaningful descriptors of solvation phenomena, the microstructural signatures of their solute-solvent intermolecular interaction asymmetry, and the thermodynamic manifestations linked to the solution nonideality. The rigorous foundations afford us to address some crucial issues frequently invoked in the literature including the microstructural perturbation domain, its proper identification and molecular-based meaning toward the interpretation of the solvation process, and the potential impact of the local differential behavior between anions and cations on the actual salt-induced perturbation of the solvent microstructure. Indeed, we link the precisely characterized species solvation behavior to fundamental thermodynamic residual-property relations, and the dynamics associated with either the viscous flow or diffusive behavior of the solvent, to finally illustrate their outcome with experimental data of aqueous electrolyte solutions from the available literature. Ultimately, this effort provides a highly desirable unambiguous identification of the cause-effect connections between the microstructurally perturbed domains and the experimentally measured macroscopic solvation properties, including their effect on the dynamics of the solvent environment. More importantly, it lends a well-established solvation framework to bridge rigorously the microstructural details of the mixture, its dynamics, and its solvation thermodynamics to enhance our understanding of well-defined ranked Hofmeister series, i.e., by avoiding ad hoc conjectures and unsupported microscopic interpretations of solvation phenomena.

Osmotic Force Balance Evaluation of Aqueous Electrolyte Osmotic Pressures and Chemical Potentials

Concentrated aqueous salt solutions are ubiquitous in problems of biological and environmental relevance. The development of accurate force fields that capture the interactions between dissolved species in solution is crucial to simulating these systems to gain molecular insights into the underlying processes under saline conditions. The osmotic pressure is a relatively simple thermodynamic property connecting the experimental and simulation measurements of the associative properties of the ions in solution. Milner [C. Gillespie and S. T. Milner, Soft Matter, 16, 9816 (2020)] proposed a simulation approach to evaluate the osmotic pressures of salts in solution by applying a restraint potential to the ions alone in solution and determining the resulting pressure required to balance that potential, referred to here as the osmotic force balance. Here, we expand Milner's approach, demonstrating that the chemical potentials of the salts in solution as a function of concentration can be fitted to the concentration profiles determined from simulation, additionally providing an analytical expression for the osmotic pressure. This approach is used to determine the osmotic pressures of 15 alkali halide salts in water from simulations. The cross interactions between cations and anions in solution are subsequently optimized to capture their experimental osmotic pressures.

Ion recognition properties of 2,2'-bibenzimidazole regulated by ammonium-modified pillar[5]arenes

With the increasing issues of environmental degradation and health problem, the selective detection of toxic ions has attracted considerable attention from researchers. Chemical fluorescent sensors with the advantages of facile operation, high sensitivity, rapid response, and easy visualization are emerging as powerful detection tools towards ions. However, the selective recognition of ions is always hindered by the presence of other interfering substances. Herein, we show that supramolecular host-guest interaction based on a pillar[5]arene provides a new opportunity to regulate the ionic recognition properties of guest molecules. A pillar[5]arene-based host-guest complex HG was constructed through the host-guest interaction between ammonium functionalized pillar[5]arene (HAP5) and 2,2'-bibenzimidazole (G). The host-gust complex HG can realize the successive, highly selective, and sensitive detection of specific ions. It was found that only in the presence of HAP5, the sensitivity towards cations was evidently enhanced, and selective successive recognition for I^- and HSO₄^- was achieved. Those results indicate that the introduction of HAP5 can effectively improve the ion recognition performance of 2,2'-bibenzimidazole, so it is a feasible strategy using supramolecular host-guest interaction to regulate the ionic recognition properties of guest molecules.

The "good," the "bad," and the "hidden" in neutron scattering and molecular dynamics of ionic aqueous solutions

We characterize a concentrated 7.3 m CaCl₂ solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca²⁺-Cl^- interactions were primarily "hidden" in the Ca-NDIS signal due to overlapping with Ca²⁺-O_w and Ca²⁺-H_w contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models "on average," the question "which model is the best" has not been answered but replaced by the question "which model is better for a given research." An overall "good" model can still be inappropriate or, in some instances, "bad" and, unfortunately, produce erroneous results. The accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models.

Moderate molecular recognitions on ZnO m-plane and their selective capture/release of bio-related phosphoric acids

Herein, we explore the hidden molecular recognition abilities of ZnO nanowires uniformly grown on the inner surface of an open tubular fused silica capillary via liquid chromatography. Chromatographic evaluation revealed that ZnO nanowires showed a stronger intermolecular interaction with phenylphosphoric acid than any other monosubstituted benzene. Furthermore, ZnO nanowires specifically recognized the phosphate groups present in nucleotides even in the aqueous mobile phase, and the intermolecular interaction increased with the number of phosphate groups. This discrimination of phosphate groups in nucleotides was unique to the rich (101̄0) m-plane of ZnO nanowires with a moderate hydrophilicity and negative charge. The discrimination could be evidenced by the changes in the infrared bands of the phosphate groups on nucleotides on ZnO nanowires. Finally, as an application of the molecular recognition, nucleotides were separated by the number of phosphate groups, utilizing optimized gradient elution on ZnO nanowire column. Thus, the present results elucidate the unique and versatile molecular selectivity of well-known ZnO nanostructures for the capture and separation of biomolecules.

Surfactant-induced wettability reversal on oil-wet calcite surfaces: Experimentation and molecular dynamics simulations with scaled-charges

HYPOTHESIS

Surfactant flooding is the leading approach for reversing the wettability of oil-wet carbonate reservoirs, which is critical for the recovery of the remaining oil. Combination of molecular dynamics (MD) simulations with experiments on simplified model systems can uncover the molecular mechanisms of wettability reversal and identify key molecular properties for systematic design of new, effective chemical formulations for the enhanced oil recovery.

EXPERIMENTS/SIMULATIONS

Wettability reversal by a series of surfactant solutions was studied experimentally using contact angle measurements on aged calcite chips, and a novel MD simulation methodology with scaled-charges that provides superior description of the ionic interactions in aqueous solutions.

FINDINGS

The MD simulation results were in excellent agreement with the experiments. Cationic surfactants were the most effective in reversing the calcite wettability, resulting in complete detachment of the oil from the surface. Some nonionic surfactants also altered the wettability, but to a lesser degree, while the amphoteric and anionic surfactants had no effect. From the tested cationic surfactants, the double-tailed one was the least effective, but the experiments were inconclusive due to its poor solubility. Contributions of specific interactions to the wettability reversal process and implications for the design and optimization of surfactants for the enhanced oil recovery are discussed.

A surface-promoted redox reaction occurs spontaneously on solvating inorganic aerosol surfaces

[Figure: see text].

Examining the Role of Different Molecular Interactions on Activation Energies and Activation Volumes in Liquid Water

There are a large number of force fields available to model water in molecular dynamics simulations, which each have their own strengths and weaknesses in describing the behavior of the liquid. One particular weakness in many of these models is their description of dynamics away from ambient conditions, where their ability to reproduce measurements is mixed. To investigate this issue, we use the recently developed fluctuation theory for dynamics to directly evaluate measures of the local temperature and pressure dependence: the activation energy and the activation volume. We examine these activation parameters for hydrogen-bond jump exchange times, OH reorientation times, and diffusion coefficients calculated from the SPC/E, SPC/Fw, TIP3P-PME, TIP3P-PME/Fw, OPC3, TIP4P/2005, TIP4P/Ew, E3B2, and E3B3 water models. Activation energy decompositions available through the fluctuation theory approach provide mechanistic insight into the origins of different temperature dependences between the various models, as well as the influence of three-body effects and flexibility.

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.

A practical guide to biologically relevant molecular simulations with charge scaling for electronic polarization

Molecular simulations can elucidate atomistic-level mechanisms of key biological processes, which are often hardly accessible to experiment. However, the results of the simulations can only be as trustworthy as the underlying simulation model. In many of these processes, interactions between charged moieties play a critical role. Current empirical force fields tend to overestimate such interactions, often in a dramatic way, when polyvalent ions are involved. The source of this shortcoming is the missing electronic polarization in these models. Given the importance of such biomolecular systems, there is great interest in fixing this deficiency in a computationally inexpensive way without employing explicitly polarizable force fields. Here, we review the electronic continuum correction approach, which accounts for electronic polarization in a mean-field way, focusing on its charge scaling variant. We show that by pragmatically scaling only the charged molecular groups, we qualitatively improve the charge-charge interactions without extra computational costs and benefit from decades of force field development on biomolecular force fields.

Charge Scaling Manifesto: A Way of Reconciling the Inherently Macroscopic and Microscopic Natures of Molecular Simulations

Electronic polarization effects play an important role in the interactions of charged species in biologically relevant aqueous solutions, such as those involving salt ions, proteins, nucleic acids, or phospholipid membranes. Explicit inclusion of electronic polarization in molecular modeling is tedious both from the point of view of force field parametrization and actual performance of the simulations. Therefore, the vast majority of biomolecular simulations is performed using nonpolarizable force fields, which can lead to artifacts such as dramatically overestimated ion pairing, particularly when polyvalent ions are involved. Here, we show that many of these issues can be remedied without extra computational costs by including electronic polarization in a mean field way via charge rescaling. We also lay the solid physical foundations of this approach and reconcile from this perspective the microscopic versus macroscopic natures of nonpolarizable force fields.

Accurate Biomolecular Simulations Account for Electronic Polarization

In this perspective, we discuss where and how accounting for electronic many-body polarization affects the accuracy of classical molecular dynamics simulations of biomolecules. While the effects of electronic polarization are highly pronounced for molecules with an opposite total charge, they are also non-negligible for interactions with overall neutral molecules. For instance, neglecting these effects in important biomolecules like amino acids and phospholipids affects the structure of proteins and membranes having a large impact on interpreting experimental data as well as building coarse grained models. With the combined advances in theory, algorithms and computational power it is currently realistic to perform simulations with explicit polarizable dipoles on systems with relevant sizes and complexity. Alternatively, the effects of electronic polarization can also be included at zero additional computational cost compared to standard fixed-charge force fields using the electronic continuum correction, as was recently demonstrated for several classes of biomolecules.

Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields

Many applications in chemistry, biology, and energy storage/conversion research rely on molecular simulations to provide fundamental insight into structural and transport properties of materials with high ionic concentrations. Whether the system is comprised entirely of ions, like ionic liquids, or is a mixture of a polar solvent with a salt, e.g., liquid electrolytes for battery applications, the presence of ions in these materials results in strong local electric fields polarizing solvent molecules and large ions. To predict properties of such systems from molecular simulations often requires either explicit or mean-field inclusion of the influence of polarization on electrostatic interactions. In this manuscript, we review the pros and cons of different treatments of polarization ranging from the mean-field approaches to the most popular explicit polarization models in molecular dynamics simulations of ionic materials. For each method, we discuss their advantages and disadvantages and emphasize key assumptions as well as their adjustable parameters. Strategies for the development of polarizable models are presented with a specific focus on extracting atomic polarizabilities. Finally, we compare simulations using polarizable and nonpolarizable models for several classes of ionic systems, discussing the underlying physics that each approach includes or ignores, implications for implementation and computational efficiency, and the accuracy of properties predicted by these methods compared to experiments.

Quantifying the Strength of a Salt Bridge by Neutron Scattering and Molecular Dynamics

The molecular structure and strength of a model salt bridge between a guanidinium cation (side chain group of arginine) and the acetate carboxylic group in an aqueous solution is characterized by a combination of neutron diffraction with isotopic substitution and molecular dynamics simulations. The present neutron scattering experiments provide direct information about ion pairing in the solution. At the same time, these measurements are used to assess the quality of the force field employed in the simulation. We show that a standard nonpolarizable force field overestimates the strength of salt bridges. In contrast, accounting for electronic polarization effects via charge scaling allows to quantitatively reproduce the experiment. Such simulations are used to quantify the weak character of a fully hydrated salt bridge. Finally, on top of the canonical hydrogen-bonding binding mode, we uncover another interaction motif involving an out-of-plane hydrophobic contact of the acetate methyl group with the guanidinium cation.

The role of pH, metal ions and their hydroxides in charge reversal of protein-coated nanoparticles

In this study, we investigated charge inversion of protein-coated Au nanoparticles caused by the addition of metal ions. Adsorbed metal hydroxides were identified to cause the charge inversion of the NPs by using a combination of cryo-TEM, EFTEM andζ-potential measurements.

Ion Hydration and Ion Pairing as Probed by THz Spectroscopy

Ion hydration is of pivotal importance for many fundamental processes. Various spectroscopic methods are used to study the retardation of the hydration bond dynamics in the vicinity of anions and cations. Here we introduce THz-FTIR spectroscopy as a powerful method to answer the open questions. We show through dissection of THz spectra that we can pinpoint characteristic absorption features that can be attributed to the rattling modes of strongly hydrating ions within their hydration cages as well as vibrationally induced charge fluctuations in the case of weakly hydrating ions. Further analysis yields information on anion-cation cooperativity, the size of the dynamic hydration shell, as well as the lifetimes of these collective ion-hydration water modes and their connecting thermal bath states. Our study provides evidence for a non-additive behavior, thus questioning the simplified Hofmeister model. THz spectroscopy enables ion pairing to be observed and quantified at a high salt concentration.

Decoding Oxyanion Aqueous Solvation Structure: A Potassium Nitrate Example at Saturation

The ability to probe the structure of a salt solution at the atomic scale is fundamentally important for our understanding of many chemical reactions and their mechanisms. The capability of neutron diffraction to "see" hydrogen (or deuterium) and other light isotopes is exceptional for resolving the structural complexity around the dissolved solutes in aqueous electrolytes. We have made measurements using oxygen isotopes on aqueous nitrate to reveal a small hydrogen-bonded water coordination number (3.9 ± 1.2) around a nitrate oxyanion. This is compared to estimates made using the existing method of nitrogen isotope substitution and those of computational simulations (>5-6 water molecules). The low water coordination number, combined with a comparison to classical molecular dynamics simulations, suggests that ion-pair formation is significant. This insight demonstrates the utility of experimental diffraction data for benchmarking atomistic computer simulations, enabling the development of more accurate intermolecular potentials.

On the mechanism of spontaneous thiol–disulfide exchange in proteins

Free energy surfaces from QM/MM metadynamics on model peptides will be a basis to build reactive force field for reaction in proteins.

Ca2+ -Cl- Association in Water Revisited: the Role of Cation Hydration

We investigate the dissociation of a Ca²⁺ -Cl^- pair in water using classical molecular dynamics simulations with a polarizable interaction potential, parameterized from ab initio calculations. By computing the potential of mean force as a function not only of the interionic distance but also of the coordination numbers by water molecules, we show that it is necessary to use a collective variable describing the cation hydration in order to capture the dissociation mechanism. In the contact ion pair, the Ca²⁺ cation has a first coordination sphere containing 5 or 6 water molecules. The minimum free-energy path for dissociation involves a two-step process: First one or two additional water molecules enter the cation coordination shell, increasing the coordination number up to 7 with an almost fixed interionic distance. Then the dissociation of the ionic pair occurs at this fixed coordination number.

Calcium Directly Regulates Phosphatidylinositol 4,5-Bisphosphate Headgroup Conformation and Recognition

The orchestrated recognition of phosphoinositides and concomitant intracellular release of Ca²⁺ is pivotal to almost every aspect of cellular processes, including membrane homeostasis, cell division and growth, vesicle trafficking, as well as secretion. Although Ca²⁺ is known to directly impact phosphoinositide clustering, little is known about the molecular basis for this or its significance in cellular signaling. Here, we study the direct interaction of Ca²⁺ with phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂), the main lipid marker of the plasma membrane. Electrokinetic potential measurements of PI(4,5)P₂ containing liposomes reveal that Ca²⁺ as well as Mg²⁺ reduce the zeta potential of liposomes to nearly background levels of pure phosphatidylcholine membranes. Strikingly, lipid recognition by the default PI(4,5)P₂ lipid sensor, phospholipase C delta 1 pleckstrin homology domain (PLC δ1-PH), is completely inhibited in the presence of Ca²⁺, while Mg²⁺ has no effect with 100 nm liposomes and modest effect with giant unilamellar vesicles. Consistent with biochemical data, vibrational sum frequency spectroscopy and atomistic molecular dynamics simulations reveal how Ca²⁺ binding to the PI(4,5)P₂ headgroup and carbonyl regions leads to confined lipid headgroup tilting and conformational rearrangements. We rationalize these findings by the ability of calcium to block a highly specific interaction between PLC δ1-PH and PI(4,5)P₂, encoded within the conformational properties of the lipid itself. Our studies demonstrate the possibility that switchable phosphoinositide conformational states can serve as lipid recognition and controlled cell signaling mechanisms.

Beyond the Hofmeister Series: Ion-Specific Effects on Proteins and Their Biological Functions

Ions differ in their ability to salt out proteins from solution as expressed in the lyotropic or Hofmeister series of cations and anions. Since its first formulation in 1888, this series has been invoked in a plethora of effects, going beyond the original salting out/salting in idea to include enzyme activities and the crystallization of proteins, as well as to processes not involving proteins like ion exchange, the surface tension of electrolytes, or bubble coalescence. Although it has been clear that the Hofmeister series is intimately connected to ion hydration in homogeneous and heterogeneous environments and to ion pairing, its molecular origin has not been fully understood. This situation could have been summarized as follows: Many chemists used the Hofmeister series as a mantra to put a label on ion-specific behavior in various environments, rather than to reach a molecular level understanding and, consequently, an ability to predict a particular effect of a given salt ion on proteins in solutions. In this Feature Article we show that the cationic and anionic Hofmeister series can now be rationalized primarily in terms of specific interactions of salt ions with the backbone and charged side chain groups at the protein surface in solution. At the same time, we demonstrate the limitations of separating Hofmeister effects into independent cationic and anionic contributions due to the electroneutrality condition, as well as specific ion pairing, leading to interactions of ions of opposite polarity. Finally, we outline the route beyond Hofmeister chemistry in the direction of understanding specific roles of ions in various biological functionalities, where generic Hofmeister-type interactions can be complemented or even overruled by particular steric arrangements in various ion binding sites.

Effect of interionic interactions on the structure and dynamics of ionic solvation shells in aqueous electrolyte solutions

Molecular dynamics (MD) simulations to explore the structure and dynamics of the ionic solvation shell of alkali ions and halide ions in aqueous solution.

The i-TTM model for ab initio-based ion–water interaction potentials. II. Alkali metal ion–water potential energy functions

A new set of i-TTM potential energy functions describing the interactions between alkali metal ions and water molecules is reported.