Effect of Ions on Water Dynamics in Dilute and Concentrated Aqueous Salt Solutions

Structural Study of Aqueous Electrolyte Solution by MeV Liquid Electron Scattering

The impact of ions on water has long been a subject of great interest, as it is closely tied to the hydration structure, dynamics, and properties of electrolyte solutions. Over centuries of investigation, the influence of ions on water's structure remains highly debated. Prevailing techniques, such as neutron and X-ray scattering, primarily focus on the microscopic structure of salt solutions at very high concentrations, mostly above 1 mol/L. In this study, we measured the structure of aqueous potassium iodide (KI) and potassium chloride (KCl) solutions using MeV liquid electron scattering (MeV-LES) across a concentration range of 0.10 to 0.75 mol/L. The obtained results provide detailed insights into the variations in ion-oxygen and oxygen-oxygen correlations as a function of concentration. The observed structural differences between KI and KCl solutions are in line with the structure maker/breaker theory, which suggests that iodide ions exert a more pronounced effect than chloride ions on disrupting the water shell. This work demonstrates the potency of MeV-LES for investigating the atomic structure in liquids, augmenting the modern analytical toolbox.

Effects of Hydrogen Bond Networks on Viscosity in Aqueous Solutions

In aqueous solutions, the impact of ions on hydrogen bond networks plays a crucial role in transport properties. We used molecular dynamics simulations to explain how ions affect viscosity through structural changes. We developed a quantitative model to describe the effect of ions on viscosity. The model comprises two parts: the addition of ions alters hydrogen bond networks, and changes in hydrogen bond networks exponentially lead to changes in viscosity. The influence of ions on hydrogen bond networks involves the following mechanisms: first, ions can disrupt the tetrahedral structures within the first solvation shell into three-coordinated structures through substitution; second, structural changes within the first shells affect the global hydrogen bond network through electrostatic forces and the hindrance of ionic volumes. By analyzing the mechanisms of how hydrogen bond networks determine viscosity through the decomposition of viscosity, we found that the proportion of potential viscosity in aqueous solutions primarily increases due to the enhancement of non-hydrogen bonding interactions, and the proportion of hydrogen bonding viscosity decreases accordingly. Our results demonstrate that hydrogen bond networks are crucial for describing the changes in transport phenomena affected by external factors.

Anionic Effects on Concentrated Aqueous Lithium Ion Dynamics

The dynamics, orientational anisotropy, diffusivity, viscosity, and density were measured for concentrated lithium salt solutions, including lithium chloride (LiCl), lithium bromide (LiBr), lithium nitrite (LiNO2), and lithium nitrate (LiNO3), with methyl thiocyanate as an infrared vibrational probe molecule, using two-dimensional infrared spectroscopy (2D IR), nuclear magnetic resonance (NMR) spectroscopy, and viscometry. The 2D IR, NMR, and viscosity results show that LiNO2 exhibits longer correlation times, lower diffusivity, and nearly 4 times greater viscosity compared to those of the other lithium salt solutions of the same concentration, suggesting that nitrite anions may strongly facilitate structure formation via strengthening water-ion network interactions, directly impacting bulk solution properties at sufficiently high concentrations. Additionally, the LiNO2 and LiNO3 solutions show significantly weakened chemical interactions between the lithium cations and the methyl thiocyanate when compared with those of the lithium halide salts.

Range and sensitivity of 17O nuclear spin-lattice relaxation as a probe of aqueous electrolyte dynamics

The study of electrolytic solutions is of relevance in many research fields, ranging from biophysics, materials, and colloid science to catalysis and electrochemistry. The dependence of solution dynamics on the nature of electrolytes and their concentrations has been the subject of many experimental and computational studies, yet it remains challenging to obtain a full understanding of the factors that govern solution behavior. Here, we provide additional insights into the behavior of aqueous solutions of alkali chlorides by combining 17O relaxation data with diffusion and viscosity data and contrast their behavior with 1H nuclear magnetic resonance relaxation data. The main findings are that 17O relaxation correlates well with viscosity data but not with diffusion data, while 1H relaxation correlates with neither. Certain ionic trends match known ion-specific series behavior, especially at high concentrations. Notably, we also examine the ranges of the interactions and conclude that the majority of the effects are tied to local water reorientation dynamics.

Water dynamics and sum-frequency generation spectra at electrode/aqueous electrolyte interfaces

The dynamics of water at interfaces between an electrode and an electrolyte is essential for the transport of redox species and for the kinetics of charge transfer reactions next to the electrode. However, while the effects of electrode potential and ion concentration on the electric double layer structure have been extensively studied, a comparable understanding of dynamical aspects is missing. Interfacial water dynamics presents challenges since it is expected to result from the complex combination of water-water, water-electrode and water-ion interactions. Here we perform molecular dynamics simulations of aqueous NaCl solutions at the interface with graphene electrodes, and examine the impact of both ion concentration and electrode potential on interfacial water reorientational dynamics. We show that for all salt concentrations water dynamics exhibits strongly asymmetric behavior: it slows down at increasingly positively charged electrodes but it accelerates at increasingly negatively charged electrodes. At negative potentials water dynamics is determined mostly by the electrode potential value, but in contrast at positive potentials it is governed both by ion-water and electrode-water interactions. We show how these strikingly different behaviors are determined by the interfacial hydrogen-bond network structure and by the ions' surface affinity. Finally, we indicate how the structural rearrangements impacting water dynamics can be probed via vibrational sum-frequency generation spectroscopy.

Water-in-Salt: Fast Dynamics, Structure, Thermodynamics, and Bulk Properties

We present concentration-dependent dynamics of highly concentrated LiBr solutions and LiCl temperature-dependent dynamics for two high concentrations and compare the results to those of prior LiCl concentration-dependent data. The dynamical data are obtained using ultrafast optical heterodyne-detected optical Kerr effect (OHD-OKE). The OHD-OKE decays are composed of two pairs of biexponentials, i.e., tetra-exponentials. The fastest decay (t1) is the same as pure water's at all concentrations within error, while the second component (t2) slows slightly with concentration. The slower components (t3 and t4), not present in pure water, slow substantially, and their contributions to the decays increase significantly with increasing concentration, similar to LiCl solutions. Simulations of LiCl solutions from the literature show that the slow components arise from large ion/water clusters, while the fast components are from ion/water structures that are not part of large clusters. Temperature-dependent studies (15-95 °C) of two high LiCl concentrations show that decreasing the temperature is equivalent to increasing the room temperature concentration. The LiBr and LiCl concentration dependences and the two LiCl concentrations' temperature dependences all have bulk viscosities that are linearly dependent on τcslow, the correlation time of the slow dynamics (weighted averages of t3 and t4). Remarkably, all four viscosity vs 1/τCslow plots fall on the same line. Application of transition state theory to the temperature-dependent data yields the activation enthalpies and entropies for the dynamics of the large ion/water clusters, which underpin the bulk viscosity.

Lifting Hofmeister's Curse: Impact of Cations on Diffusion, Hydrogen Bonding, and Clustering of Water

Water plays a role in the stability, reactivity, and dynamics of the solutes that it contains. The presence of ions alters this capacity by changing the dynamics and structure of water. However, our understanding of how and to what extent this occurs is still incomplete. Here, a study of the low-frequency Raman spectra of aqueous solutions of various cations by using optical Kerr-effect spectroscopy is presented. This technique allows for the measurement of the changes that ions cause in both the diffusive dynamics and the vibrations of the hydrogen-bond structure of water. It is found that when salts are added, some of the water molecules become part of the ion solvation layers, while the rest retain the same diffusional properties as those of pure water. The slowing of the dynamics of the water molecules in the solvation shell of each ion was found to depend on its charge density at infinite dilution conditions and on its position in the Hofmeister series at higher concentrations. It is also observed that all cations weaken the hydrogen-bond structure of the solution and that this weakening depends only on the size of the cation. Finally, evidence is found that ions tend to form amorphous aggregates, even at very dilute concentrations. This work provides a novel approach to water dynamics that can be used to better study the mechanisms of solute nucleation and crystallization, the structural stability of biomolecules, and the dynamic properties of complex solutions, such as water-in-salt electrolytes.

Anomalous Concentration Dependence of Viscosity: Hidden Role of Cross-Correlations in Aqueous Electrolyte Solutions

The viscosity of aqueous electrolyte solutions exhibits well-known composition-dependent anomalies that show certain definitive trends and universal features. The viscosity of LiCl and NaCl solutions increases with concentration in a monotonic fashion, while solutions of KCl, RbCl, and CsCl exhibit a more complex behavior. Here, the viscosity first decreases and then increases with increasing concentration, with a rather broad minimum at intermediate concentrations (ca. 1-3 m). To unearth the origin of such puzzling behavior, we carried out detailed molecular-level analyses by interrogating the exact Green-Kubo expression of viscosity in terms of the stress-stress time correlation function (SS-TCF). The total SS-TCF can be decomposed into a collection of three self- and three cross-SS-TCFs arising from the three constituent components (water, cations, and anions). Mode coupling theory (MCT) analysis for the friction on ions and the viscosity of the solution suggests the possible importance of two-particle static and time-dependent cross-correlations between water and the ions. We calculate the viscosity and other dynamical properties for all five electrolyte (LiCl, NaCl, KCl, RbCl, and CsCl) solutions over a range of concentrations, using two models of water (SPC/E and TIP4P/2005). The total viscosity derives non-negligible contributions from all of the terms. The cross-correlations are found to be surprisingly large and seen to play a hidden role in the concentration dependence. However, the importance of cross-correlations is often not discussed. Our study leads to a theoretical understanding of the microscopic origin of the observed anomalies in the composition dependence of viscosity across all five electrolytes.

Ab Initio Molecular Dynamics Investigation of the Solvation States of Hydrated Ions in Confined Water

Ionic transport in nanoscale channels with a critical size comparable to that of ions and solutes exhibits exceptional performance in water desalination, ion separation, electrocatalysts, and supercapacitors. However, the solvation states (SSs), i.e., the hydration structures and probability distribution, of hydrated ions in nanochannels differ from those in the bulk and the perspective of continuum theory. In this work, we conduct ab initio enhanced-sampling atomistic simulations to investigate the ion-specific SSs of monovalent ions (including Li+, Na+, K+, F-, Cl-, and I-) in the graphene channel with a width of 1 nm. Our findings highlight that the SSs of those ions are primarily determined by ion-water hydration, where ion-wall interactions play a minor role. The distribution of ions in layered confined water is a result of ion-specific hydration, which arises from the synergy of entropy and enthalpy. The free energy barriers for transitions between SSs are on the order of 1kBT, allowing for modulation through applying external fields or modifying surface properties. As the ion-wall interaction strengthens, as observed in vermiculite and carbides and nitrides of transition metal channels, the probability of near-wall SSs increases. These results help to improve the performance of nanofluidic devices and provide crucial insights for developing accurate force fields of molecular simulations or advanced theoretical approaches for ion dynamics in confined channels.

Understanding the Anomalous Diffusion of Water in Aqueous Electrolytes Using Machine Learned Potentials

The diffusivity of water in aqueous cesium iodide solutions is larger than that in neat liquid water and vice versa for sodium chloride solutions. Such peculiar ion-specific behavior, called anomalous diffusion, is not reproduced in typical force field based molecular dynamics (MD) simulations due to inadequate treatment of ion-water interactions. Herein, this hurdle is tackled by using machine learned atomic potentials (MLPs) trained on data from density functional theory calculations. MLP based atomistic MD simulations of aqueous salt solutions reproduce experimentally determined thermodynamic, structural, dynamical, and transport properties, including their varied trends in water diffusivities across salt concentration. This enables an examination of their intermolecular structure to unravel the microscopic underpinnings of the differences in their transport properties. While both ions in CsI solutions contribute to the faster diffusion of water molecules, the competition between the heavy retardation by Na ions and the slight acceleration by Cl ions in NaCl solutions reduces their water diffusivity.

Electrical noise in electrolytes: a theoretical perspective

Seemingly unrelated experiments such as electrolyte transport through nanotubes, nano-scale electrochemistry, NMR relaxometry and surface force balance measurements, all probe electrical fluctuations: of the electric current, the charge and polarization, the field gradient (for quadrupolar nuclei) and the coupled mass/charge densities. The fluctuations of such various observables arise from the same underlying microscopic dynamics of the ions and solvent molecules. In principle, the relevant length and time scales of these dynamics are encoded in the dynamic structure factors. However, modelling the latter for frequencies and wavevectors spanning many orders of magnitude remains a great challenge to interpret the experiments in terms of physical processes such as solvation dynamics, diffusion, electrostatic and hydrodynamic interactions between ions, interactions with solid surfaces, etc. Here, we highlight the central role of the charge-charge dynamic structure factor in the fluctuations of electrical observables in electrolytes and offer a unifying perspective over a variety of complementary experiments. We further analyze this quantity in the special case of an aqueous NaCl electrolyte, using simulations with explicit ions and an explicit or implicit solvent. We discuss the ability of the standard Poisson-Nernst-Planck theory to capture the simulation results, and how the predictions can be improved. We finally discuss the contributions of ions and water to the total charge fluctuations. This work illustrates an ongoing effort towards a comprehensive understanding of electrical fluctuations in bulk and confined electrolytes, in order to enable experimentalists to decipher the microscopic properties encoded in the measured electrical noise.

Dynamics of Concentrated Aqueous Lithium Chloride Solutions Investigated with Optical Kerr Effect Experiments

We report the dynamics of concentrated lithium chloride aqueous solutions over a range of moderate to high concentrations. Concentrations (1-29 to 1-3.3 LiCl-water) were studied in which, at the highest concentrations, there are far too few water molecules to solvate the ions. The measurements were made with optically heterodyne-detected optical Kerr effect experiments, a non-resonant technique able to observe dynamics over a wide range of time scales and signal amplitudes. While the pure water decay is a biexponential, the LiCl-water decays are tetra-exponentials at all concentrations. The faster two decays arise from water dynamics, while the slower two decays reflect the dynamics of the ion-water network. The fastest decay (t₁) is the same as pure water at all concentrations. The second decay (t₂) is also the same as that of pure water at the lower concentrations, and then, it slows with increasing concentration. The slower dynamics (t₃ and t₄), which do not have counterparts in pure water, arise from ion-water complexes and, at the highest concentrations, an extended ion-water network. Comparisons are made between the concentration dependence of the observed dynamics and simulations of structural changes from the literature, which enable the assignment of dynamics to specific ion-water structures. The concentration dependences of the bulk viscosity and the ion-water network dynamics are directly correlated. The correlation provides an atomistic-level understanding of the viscosity.

Shining (Infrared) Light on the Hofmeister Series: Driving Forces for Changes in the Water Vibrational Spectra in Alkali-Halide Salt Solutions

The Hofmeister series is frequently used to rank ions based on their behavior from chaotropes ("structure breakers"), which weaken the surrounding hydrogen-bond network, to kosmotropes ("structure makers"), which enhance it. Here, we use fluctuation theory to investigate the energetic and entropic driving forces underlying the Hofmeister series for aqueous alkali-halide solutions. Specifically, we exploit the OH stretch infrared (IR) spectrum in isotopically dilute HOD/D₂O solutions as a probe of the effect of the salt on the water properties for different concentrations and choice of halide anion. Fluctuation theory is used to calculate the temperature derivative of these IR spectra, including decomposition of the derivative into different energetic contributions. These contributions are used to determine the thermodynamic driving forces in terms of effective internal energy and entropic contributions. This analysis implicates entropic contributions as the key factor in the Hofmeister series behavior of the OH stretch IR spectra, while the effective internal energy is nearly ion-independent.

Probing electrolyte–silica interactions through simulations of the infrared spectroscopy of nanoscale pores

The structural and dynamical properties of nanoconfined solutions can differ dramatically from those of the corresponding bulk systems. Understanding the changes induced by confinement is central to controlling the behavior of synthetic nanostructured materials and predicting the characteristics of biological and geochemical systems. A key outstanding issue is how the molecular-level behavior of nanoconfined electrolyte solutions is reflected in different experimental, particularly spectroscopic, measurements. This is addressed here through molecular dynamics simulations of the OH stretching infrared (IR) spectroscopy of NaCl, NaBr, and NaI solutions in isotopically dilute HOD/D2O confined in hydroxylated amorphous silica slit pores of width 1–6 nm and pH [Formula: see text]. In addition, the water reorientation dynamics and spectral diffusion, accessible by pump–probe anisotropy and two-dimensional IR measurements, are investigated. The aim is to elucidate the effect of salt identity, confinement, and salt concentration on the vibrational spectra. It is found that the IR spectra of the electrolyte solutions are only modestly blue-shifted upon confinement in amorphous silica slit pores, with both the size of the shift and linewidth increasing with the halide size, but these effects are suppressed as the salt concentration is increased. This indicates the limitations of linear IR spectroscopy as a probe of confined water. However, the OH reorientational and spectral diffusion dynamics are significantly slowed by confinement even at the lowest concentrations. The retardation of the dynamics eases with increasing salt concentration and pore width, but it exhibits a more complex behavior as a function of halide.

Dielectric Properties of Water in Charged Nanopores

In this study, we examine the spectral dielectric properties of liquid water in charged nanopores over a wide range of frequencies (0.3 GHz to 30 THz) and pore widths (0.3 to 5 nm). This has been achieved using classical molecular dynamics simulations of hydrated Na-smectite, the prototypical swelling clay mineral. We observe a drastic (20-fold) and anisotropic decrease in the static relative permittivity of the system as the pore width decreases. This large decrement in static permittivity reflects a strong attenuation of the main Debye relaxation mode of liquid water. Remarkably, this strong attenuation entails very little change in the time scale of the collective relaxation. Our results indicate that water confined in charged nanopores is a distinct solvent with a much weaker collective nature than bulk liquid water, in agreement with recent observations of water in uncharged nanopores. Finally, we observe remarkable agreement between the dielectric properties of the simulated clay system against a compiled set of soil samples at various volumetric water contents. This implies that saturation may not be the sole property dictating the dielectric properties of soil samples, rather that the pore-size distribution of fully saturated nanopores may also play a critically important role.

Bridging atomistic simulations and thermodynamic hydration models of aqueous electrolyte solutions

Chemical thermodynamic models of solvent and solute activities predict the equilibrium behavior of aqueous solutions. However, these models are semi-empirical. They represent micro-scale ion and solvent behaviors controlling the macroscopic properties using small numbers of parameters whose values are obtained by fitting to activities and other partial derivatives of the Gibbs energy measured for the bulk solutions. We have conducted atomistic simulations of aqueous electrolyte solutions (MgCl₂ and CaCl₂) to determine the parameters of thermodynamic hydration models. We have implemented a cooperative hydration model to categorize the water molecules in electrolyte solutions into different subpopulations. The value of the electrolyte-specific parameter, k, was determined from the ion-affected subpopulation with the lowest absolute value of the free energy of removing the water molecule. The other equilibrium constant parameter, K₁, associated with the first degree of hydration, was computed from the free energy of hydration of hydrated clusters. The hydration number, h, was determined from a reorientation dynamic analysis of the water subpopulations compared to bulk-like behavior. The reparameterized models [R. H. Stokes and R. H. Robinson, J. Solution Chem. 2, 173 (1973) and Balomenos et al., Fluid Phase Equilib. 243, 29 (2006)] using the computed values of the parameters lead to the osmotic coefficients of MgCl₂ solutions that are consistent with measurements. Such an approach removes the dependence on the availability of experimental data and could lead to aqueous thermodynamic models capable of estimating the values of solute and solvent activities as well as thermal and volumetric properties for a wide range of compositions and concentrations.

Bridging Structure, Dynamics, and Thermodynamics: An Example Study on Aqueous Potassium Halides

Aqueous salt systems are ubiquitous in all areas of life. The ions in these solutions impose important structural and dynamic perturbations to water. In this study, we employ a combined neutron scattering, nuclear magnetic resonance, and computational modeling approach to deconstruct ion-specific perturbations to water structure and dynamics and shed light on the molecular origins of bulk thermodynamic properties of the solutions. Our approach uses the atomistic scale resolution offered to us by neutron scattering and computational modeling to investigate how the properties of particular short-ranged microenvironments within aqueous systems can be related to bulk properties of the system. We find that by considering only the water molecules in the first hydration shell of the ions that the enthalpy of hydration can be determined. We also quantify the range over which ions perturb water structure by calculating the average enthalpic interaction between a central halide anion and the surrounding water molecules as a function of distance and find that the favorable anion-water enthalpic interactions only extend to ∼4 Å. We further validate this by showing that ions induce structure in their solvating water molecules by examining the distribution of dipole angles in the first hydration shell of the ions but that this perturbation does not extend into the bulk water. We then use these structural findings to justify mathematical models that allow us to examine perturbations to rotational and diffusive dynamics in the first hydration shell around the potassium halide ions from NMR measurements. This shows that as one moves down the halide series from fluorine to iodine, and ionic charge density is therefore reduced, that the enthalpy of hydration becomes less negative. The first hydration shell also becomes less well structured, and rotational and diffusive motions of the hydrating water molecules are increased. This reduction in structure and increase in dynamics are likely the origin of the previously observed increased entropy of hydration as one moves down the halide series. These results also suggest that simple monovalent potassium halide ions induce mostly local perturbations to water structure and dynamics.

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.

Examining the Hofmeister Series through Activation Energies: Water Diffusion in Aqueous Alkali-Halide Solutions

The effect of ions on the properties of aqueous solutions is often categorized in terms of the Hofmeister series that ranks them from chaotropes ("structure-breakers"), which weaken the surrounding hydrogen-bond network to kosmotropes ("structure-makers"), which enhance it. Here, we investigate the Hofmeister series in ∼1 M sodium-halide solutions using molecular dynamics simulations to calculate the effect of the identity and proximity of the halide anion on both the water diffusion coefficient and its activation energy. A recently developed method for calculating the activation energy from a single-temperature simulation is used, which also permits a rigorous decomposition into contributions from different interactions and motions. The mechanisms of the salt effects on the water dynamics are explored by separately considering water molecules based on their location relative to the ions. The results show that water diffusion is accelerated moving down the halide group from F^- to I^-. The behavior of the diffusion activation energy, E_a, is more complex, indicating a significant role for entropic effects. However, water molecules in the first or second solvation shell of an ion exhibit a decrease in E_a moving down the halide series and Na⁺ exhibits a larger effect than any of the anions. The E_a for water molecules within the second solvation shell of an ion are modest, indicating a short-ranged nature of the ion influence.

Correlated atomic dynamics in liquid seen in real space and time

In liquids, the timescales for structure, diffusion, and phonon are all similar, of the order of a pico-second. This not only makes characterization of liquid dynamics difficult but also renders it highly questionable to describe liquids in these terms. In particular, the current definition of the structure of liquids by the instantaneous structure may need to be expanded because the liquid structure is inherently dynamic. Here, we advocate describing the liquid structure through the distinct-part of the Van Hove function, which can be determined by inelastic neutron and x-ray scattering measurements as well as by simulation. It depicts the dynamic correlation between atoms in space and time, starting with the instantaneous correlation function at t = 0. The observed Van Hove functions show that the atomic dynamics is strongly correlated in some liquids, such as water. The effect of atomic correlation on various transport properties of fluid, including viscosity and diffusivity, is discussed.

Molecular modeling of aqueous electrolytes at interfaces: Effects of long-range dispersion forces and of ionic charge rescaling

Molecular dynamics simulations of aqueous electrolytes generally rely on empirical force fields, combining dispersion interactions-described by a truncated Lennard-Jones (LJ) potential-and electrostatic interactions-described by a Coulomb potential computed with a long-range solver. Recently, force fields using rescaled ionic charges [electronic continuum correction (ECC)], possibly complemented with rescaling of LJ parameters [ECC rescaled (ECCR)], have shown promising results in bulk, but their performance at interfaces has been less explored. Here, we started by exploring the impact of the LJ potential truncation on the surface tension of a sodium chloride aqueous solution. We show a discrepancy between the numerical predictions for truncated LJ interactions with a large cutoff and for untruncated LJ interactions computed with a long-range solver, which can bias comparison of force field predictions with experiments. Using a long-range solver for LJ interactions, we then show that an ionic charge rescaling factor chosen to correct long-range electrostatic interactions in bulk accurately describes image charge repulsion at the liquid-vapor interface, and the rescaling of LJ parameters in ECCR models-aimed at capturing local ion-ion and ion-water interactions in bulk- describes well the formation of an ionic double layer at the liquid-vapor interface. Overall, these results suggest that the molecular modeling of aqueous electrolytes at interfaces would benefit from using long-range solvers for dispersion forces and from using ECCR models, where the charge rescaling factor should be chosen to correct long-range electrostatic interactions.

Influence of Calcium Binding on Conformations and Motions of Anionic Polyamino Acids. Effect of Side Chain Length

Investigation of the effect of CaCl2 salt on conformations of two anionic poly(amino acids) with different side chain lengths, poly-(α-l glutamic acid) (PGA) and poly-(α-l aspartic acid) (PASA), was performed by atomistic molecular dynamics (MD) simulations. The simulations were performed using both unbiased MD and the Hamiltonian replica exchange (HRE) method. The results show that at low CaCl2 concentration adsorption of Ca2+ ions lead to a significant chain size reduction for both PGA and PASA. With the increase in concentration, the chains sizes partially recover due to electrostatic repulsion between the adsorbed Ca2+ ions. Here, the side chain length becomes important. Due to the longer side chain and its ability to distance the charged groups with adsorbed ions from both each other and the backbone, PGA remains longer in the collapsed state as the CaCl2 concentration is increased. The analysis of the distribution of the mineral ions suggests that both poly(amino acids) should induce the formation of mineral with the same structure of the crystal cell.

Water dynamics at electrified graphene interfaces: a jump model perspective

The reorientation dynamics of water at electrified graphene interfaces was recently shown [J. Phys. Chem. Lett., 2020, 11, 624-631] to exhibit a surprising and strongly asymmetric behavior: positive electrode potentials slow down interfacial water reorientation, while for increasingly negative potentials water dynamics first accelerates before reaching an extremum and then being retarded for larger potentials. Here we use classical molecular dynamics simulations to determine the molecular mechanisms governing water dynamics at electrified interfaces. We show that changes in water reorientation dynamics with electrode potential arise from the electrified interfaces' impacts on water hydrogen-bond jump exchanges, and can be quantitatively described by the extended jump model. Finally, our simulations indicate that no significant dynamical heterogeneity occurs within the water interfacial layer next to the weakly interacting graphene electrode.

Trimethylamine N-oxide (TMAO) resists the compression of water structure by magnesium perchlorate: terrestrial kosmotrope vs. Martian chaotrope

Neutron diffraction and computational modelling provide insight into water structure.

A force field of Li+, Na+, K+, Mg2+, Ca2+, Cl-, and SO4 2- in aqueous solution based on the TIP4P/2005 water model and scaled charges for the ions

In this work, a force field for several ions in water is proposed. In particular, we consider the cations Li⁺, Na⁺, K⁺, Mg²⁺, and Ca²⁺ and the anions Cl^- and SO₄ ^2-. These ions were selected as they appear in the composition of seawater, and they are also found in biological systems. The force field proposed (denoted as Madrid-2019) is nonpolarizable, and both water molecules and sulfate anions are rigid. For water, we use the TIP4P/2005 model. The main idea behind this work is to further explore the possibility of using scaled charges for describing ionic solutions. Monovalent and divalent ions are modeled using charges of 0.85 and 1.7, respectively (in electron units). The model allows a very accurate description of the densities of the solutions up to high concentrations. It also gives good predictions of viscosities up to 3 m concentrations. Calculated structural properties are also in reasonable agreement with the experiment. We have checked that no crystallization occurred in the simulations at concentrations similar to the solubility limit. A test for ternary mixtures shows that the force field provides excellent performance at an affordable computer cost. In summary, the use of scaled charges, which could be regarded as an effective and simple way of accounting for polarization (at least to a certain extend), improves the overall description of ionic systems in water. However, for purely ionic systems, scaled charges will not adequately describe neither the solid nor the melt.

The role of water in amyloid aggregation kinetics

The role of water in protein function and aggregation is highly important and may hold some answers to understanding initiation of misfolding diseases such as Parkinson's, Alzheimer's and Huntington's where soluble intrinsically disordered proteins (IDPs) aggregate into fibrillar structures. IDPs are highly dynamic and have larger solvent exposed areas compared to globular proteins, meaning they make and break hydrogen bonds with the surrounding water more frequently. The mobility of water can be altered by presence of ions, sugars, osmolytes, proteins and membranes which differ in different cell types, cell compartments and also as cells age. A reduction in water mobility and thus protein mobility enhances the probability that IDPs can associate to form intermolecular bonds and propagate into aggregates. This poses an interesting question as to whether localised water mobility inside cells can influence the propensity of an IDP to aggregate and furthermore whether it can influence fibril polymorphism and disease outcome.