A potential model for sodium chloride solutions based on the TIP4P/2005 water model

Formation of compounds with diverse polyelectrolyte morphologies and nonlinear ion conductance in a two-dimensional nanofluidic channel

Aqueous electrolytes subjected to angstrom-scale confinement have recently attracted increasing interest because of their distinctive structural and transport properties, as well as their promising applicability in bioinspired nanofluidic iontronics and ion batteries. Here, we performed microsecond-scale molecular dynamics simulations, which provided evidence of nonlinear ionic conductance under an external lateral electric field due to the self-assembly of cations and anions with diverse polyelectrolyte morphologies (e.g., extremely large ion clusters) in aqueous solutions within angstrom-scale slits. Specifically, we found that the cations and anions of Li2SO4 and CaSO4 formed chain-like polyelectrolyte structures, whereas those of Na2SO4 and MgSO4 predominantly formed a monolayer of hydrated salt. Additionally, the cations and anions of K2SO4 assembled into a hexagonal anhydrous ionic crystal. These ion-dependent diverse polyelectrolyte morphologies stemmed from the enhanced Coulomb interactions, weakened hydration and steric constraints within the angstrom-scale slits. More importantly, once the monolayer hydrated salt or ionic crystal structure was formed, the field-induced ion current exhibited an intriguing gating effect at a low field strength. This abnormal ion transport was attributed to the concerted movement of cations and anions within the solid polyelectrolytes, leading to the suppression of ion currents. When the electric field exceeded a critical strength, however, the ion current surged rapidly due to the dissolution of many cations and anions within a few nanoseconds in the aqueous solution.

Accurate Free Energies of Aqueous Electrolyte Solutions from Molecular Simulations with Non-polarizable Force Fields

Non-polarizable force fields fail to accurately predict free energies of aqueous electrolytes without compromising the predictive ability for densities and transport properties. A new approach is presented in which (1) TIP4P/2005 water and scaled charge force fields are used to describe the interactions in the liquid phase and (2) an additional Effective Charge Surface (ECS) is used to compute free energies at zero additional computational expense. The ECS is obtained using a single temperature-independent charge scaling parameter per species. Thereby, the chemical potential of water and the free energies of hydration of various aqueous salts (e.g., NaCl and LiCl) are accurately described (deviations less than 5% from experiments), in sharp contrast to calculations where the ECS is omitted (deviations larger than 20%). This approach enables accurate predictions of free energies of aqueous electrolyte solutions using non-polarizable force fields, without compromising liquid-phase properties.

The temperature of maximum density for aqueous solutions

Experimental and theoretical advances for understanding the temperature of maximum density (TMD) of aqueous solutions are outlined. The main equations that relate the TMD behavior to key thermodynamic properties are stated. The experimental TMD data are classified as a function of the nature of the solute (inorganic electrolytes, non-electrolytes, organic salts and ionic liquids, and amino acids and proteins). In addition, the experimental results that explore the effect of pressure are detailed. These experimental data are rationalized by making use of qualitative and semi-quantitative arguments based on the thermodynamics of aqueous systems. The main theoretical and simulation advances in TMD for aqueous solutions are also shown-including new calculations in the context of the scaled particle theory-and their ability to reproduce the experimental data is evaluated. Finally, new experiments and theoretical and simulation developments, which could give important insights into the problem of TMD for aqueous solutions, are proposed.

Insights into Molecular Dynamics and Oil Extraction Behavior of the Polymeric Surfactant in a Multilayered Heterogeneous Reservoir

Interlayer heterogeneity, an inevitable and complex challenge during water flooding, seriously constrains the spread of the sweep region and oil recovery enhancement in multilayered heterogeneous reservoirs. To overcome this challenge, a novel polymeric surfactant, having an excellent performance in the reduction of interfacial tension (IFT) and the increase of viscosity of displacing fluid, is applied for enlarging the sweep resonance and increasing the oil washing efficiency. Through the molecular dynamics (MD) simulation, the molecular distributing mechanisms of the polymeric surfactant at the oil-water interface are analyzed to provide the theoretical basis for explaining the microscopic mechanism of oil extraction. To directly reflect the microscopic behavior of oil extraction, multiple transparent sand-packed models are designed to investigate the flowing behavior of different fluids and the extracted mechanisms of the remaining oil in both pore and macroscales. The multilayered heterogeneous reservoirs consisting of high-, moderate-, and low-permeability layers are fabricated to represent a heterogeneous characteristic. The recognition from the visual experiment and MD simulation can study the performance control, the extracting performance of the remaining oil, and the expression of the displacing front from different perspectives. The results from MD simulation demonstrate that the polymeric surfactant can promote the disintegration of the remaining oil and enhance its mobility. The experimental results indicate that the sweep efficiency is restricted by viscous fingering and tongue advance. Through the analysis of mathematical models, the rising mobility ratio and the location of the displacing front have a strong positive relationship with viscous fingering and tongue advance, which can reasonably explain the plugging performance of the polymeric surfactant, greatly improving the sweeping effect of the whole reservoir. Moreover, the Marangoni effect generated by the IFT gradient can induce the transformation of interfacial energy to displacement kinetic energy by the emulsification of the oil-water interface so that the remaining oil in the blind-end pore can be effectively extracted. However, by comparing data from image quantification techniques and production dynamic performance, the sweep efficiency (484%) was significantly greater than that of oil recovery (300%), demonstrating that the expanded sweep effect still plays a dominant role in the extraction of remaining oil after polymeric surfactant flooding. This study provides a novel plugging and effective washing agent that is expected to be an excellent and comprehensive method for solving the problem of low oil recovery in multilayered heterogeneous reservoirs.

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.

Effects of ion type and concentration on the structure and aggregation of the amyloid peptide A β 16 - 22 $$ {\boldsymbol{\beta}}_{16-22} $$

Among the various factors controlling the amyloid aggregation process, the influences of ions on the aggregation rate and the resulting structures are important aspects to consider, which can be studied by molecular simulations. There is a wide variety of protein force fields and ion models, raising the question of which model to use in such studies. To address this question, we perform molecular dynamics simulations of Aβ16-22 , a fragment of the Alzheimer's amyloid β peptide, using different protein force fields, AMBER99SB-disp (A99-d) and CHARMM36m (C36m), and different ion parameters. The influences of NaCl and CaCl2 at various concentrations are studied and compared with the systems without the addition of ions. Our results indicate a sensitivity of the peptide-ion interactions to the different ion models. In particular, we observe a strong binding of Ca2+ to residue E22 with C36m and also with the Åqvist ion model used together with A99-d, which slightly affects the monomeric Aβ16-22 structures and the aggregation rate, but significantly affects the oligomer structures formed in the aggregation simulations. For example, at high Ca2+ concentrations, there was a switch from an antiparallel to a parallel β-sheet. Such ionic influences are of biological relevance because local ion concentrations can change in vivo and could help explain the polymorphism of amyloid fibrils.

Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property

In this work, we computed electrical conductivities under ambient conditions of aqueous NaCl and KCl solutions by using the Einstein-Helfand equation. Common force fields (charge q = ±1 e) do not reproduce the experimental values of electrical conductivities, viscosities, and diffusion coefficients. Recently, we proposed the idea of using different charges to describe the potential energy surface (PES) and the dipole moment surface (DMS). In this work, we implement this concept. The equilibrium trajectories required to evaluate electrical conductivities (within linear response theory) were obtained by using scaled charges (with the value q = ±0.75 e) to describe the PES. The potential parameters were those of the Madrid-Transport force field, which accurately describe viscosities and diffusion coefficients of these ionic solutions. However, integer charges were used to compute the conductivities (thus describing the DMS). The basic idea is that although the scaled charge describes the ion-water interaction better, the integer charge reflects the value of the charge that is transported due to the electric field. The agreement obtained with experiments is excellent, as for the first time electrical conductivities (and the other transport properties) of NaCl and KCl electrolyte solutions are described with high accuracy for the whole concentration range up to their solubility limit. Finally, we propose an easy way to obtain a rough estimate of the actual electrical conductivity of the potential model under consideration using the approximate Nernst-Einstein equation, which neglects correlations between different ions.

Self-Consistent Implementation of a Solvation Free Energy Framework to Predict the Salt Solubilities of Six Alkali Halides

To assess the salt solubilities of six alkali halides in aqueous systems, we proposed a thermodynamic cycle and an efficient molecular modeling methodology. The Gibbs free energy changes for vaporization, dissociation, and dissolution were calculated using the experimental data of ionic thermodynamic properties obtained from the NBS tables. Additionally, the Marcus' and Tissandier's solvation free energy data for Li+, Na+, K+, Cl-, and Br- ions were compared with the conventional solvation free energies by substituting into our self-consistent thermodynamic cycle. Furthermore, Tissandier's absolute solvation free energy data were used as the training set to refit the Lennard-Jones parameters of OPLS-AA force field for ions. To predict salt solubilities, an assumption of a pseudo-solvent was proposed to characterize the coupling work of a solute with its environment from infinitely diluted to saturated solutions, indicating that the Gibbs energy change of solvation process is a function of ionic strength. Following the self-consistency of the cycle, the newly derived formulas were used to determine the salt solubilities by interpolating the intersection of Gibbs free energy of dissolution and the zero free energy line. The refined ion parameters can also predict the structure and thermodynamic properties of aqueous electrolyte solutions, such as densities, pair correlation functions, hydration numbers, mean activity coefficients, vapor pressures, and the radial dependences of the net charge at 298.15 K and 1 bar. Our method can be used to characterize the solid-liquid equilibria of ions or charged particles in aqueous systems. Furthermore, for highly concentrated strong electrolyte systems, it is essential to introduce accurate water models and polarizable force fields.

Molecular modelling of the thermophysical properties of fluids: expectations, limitations, gaps and opportunities

This manuscript provides an overview of the current state of the art in terms of the molecular modelling of the thermophysical properties of fluids. It is intended to manage the expectations and serve as guidance to practising physical chemists, chemical physicists and engineers in terms of the scope and accuracy of the more commonly available intermolecular potentials along with the peculiarities of the software and methods employed in molecular simulations while providing insights on the gaps and opportunities available in this field. The discussion is focused around case studies which showcase both the precision and the limitations of frequently used workflows.

Scaled charges for ions: An improvement but not the final word for modeling electrolytes in water

In this work, we discuss the use of scaled charges when developing force fields for NaCl in water. We shall develop force fields for Na⁺ and Cl^- using the following values for the scaled charge (in electron units): ±0.75, ±0.80, ±0.85, and ±0.92 along with the TIP4P/2005 model of water (for which previous force fields were proposed for q = ±0.85 and q = ±1). The properties considered in this work are densities, structural properties, transport properties, surface tension, freezing point depression, and maximum in density. All the developed models were able to describe quite well the experimental values of the densities. Structural properties were well described by models with charges equal to or larger than ±0.85, surface tension by the charge ±0.92, maximum in density by the charge ±0.85, and transport properties by the charge ±0.75. The use of a scaled charge of ±0.75 is able to reproduce with high accuracy the viscosities and diffusion coefficients of NaCl solutions for the first time. We have also considered the case of KCl in water, and the results obtained were fully consistent with those of NaCl. There is no value of the scaled charge able to reproduce all the properties considered in this work. Although certainly scaled charges are not the final word in the development of force fields for electrolytes in water, its use may have some practical advantages. Certain values of the scaled charge could be the best option when the interest is to describe certain experimental properties.

Dynamics of Aqueous Electrolyte Solutions: Challenges for Simulations

This Perspective article focuses on recent simulation work on the dynamics of aqueous electrolytes. It is well-established that full-charge, nonpolarizable models for water and ions generally predict solution dynamics that are too slow in comparison to experiments. Models with reduced (scaled) charges do better for solution diffusivities and viscosities but encounter issues describing other dynamic phenomena such as nucleation rates of crystals from solution. Polarizable models show promise, especially when appropriately parametrized, but may still miss important physical effects such as charge transfer. First-principles calculations are starting to emerge for these properties that are in principle able to capture polarization, charge transfer, and chemical transformations in solution. While direct ab initio simulations are still too slow for simulations of large systems over long time scales, machine-learning models trained on appropriate first-principles data show significant promise for accurate and transferable modeling of electrolyte solution dynamics.

Entropy scaling of viscosity for molecular models of molten salts

Entropy scaling relates dynamic and thermodynamic properties by reducing the viscosity to a function of only the residual entropy. Molecular simulations are used to investigate the entropy scaling of the viscosity of three models of sodium chloride and five monovalent salts. Even though the correlation between the potential energy and the virial is weak, entropy scaling applies at liquid densities for all models and salts investigated. At lower densities, entropy scaling breaks down due to the formation of ion pairs and chains. Entropy scaling can be used to develop more extendable correlations for the dynamic properties of molten salts.

The Contribution of the Ion-Ion and Ion-Solvent Interactions in a Molecular Thermodynamic Treatment of Electrolyte Solutions

Developing molecular equations of state to treat electrolyte solutions is challenging due to the long-range nature of the Coulombic interactions. Seminal approaches commonly used are the mean spherical approximation (MSA) and the Debye-Hückel (DH) theory to account for ion-ion interactions and, often, the Born theory of solvation for ion-solvent interactions. We investigate the accuracy of the MSA and DH approaches using each to calculate the contribution of the ion-ion interactions to the chemical potential of NaCl in water, comparing these with newly computer-generated simulation data; the ion-ion contribution is isolated by selecting an appropriate primitive model with a Lennard-Jones force field to describe the solvent. A study of mixtures with different concentrations and ionic strengths reveals that the calculations from both MSA and DH theories are of similar accuracy, with the MSA approach resulting in marginally better agreement with the simulation data. We also demonstrate that the Born theory provides a good qualitative description of the contribution of the ion-solvent interactions; we employ an explicitly polar water model in these simulations. Quantitative agreement up to moderate salt concentrations and across the relevant range of temperature is achieved by adjusting the Born radius using simulation data of the free energy of solvation. We compute the radial and orientational distribution functions of the systems, thereby providing further insight on the differences observed between the theory and simulation. We thus provide rigorous benchmarks for use of the MSA, DH, and Born theories as perturbation approaches, which will be of value for improving existing models of electrolyte solutions, especially in the context of equations of state.

Obtaining Consistent Free Energies for Ion Binding at Surfaces from Solution: Pathways versus Alchemy for Determining Kink Site Stability

Ion incorporation or removal from a solid at the interface with solution is a fundamental part of crystal growth. Despite this, there have been few quantitative determinations of the thermodynamics for such processes from atomistic molecular dynamics due to the associated technical challenges. In this study, we compute the free energies for ion removal from kink sites at the interface between NaCl and water as an illustrative example. To examine the influence of the free energy technique used, we compare methods that follow an explicit pathway for dissolution with those that focus on the thermodynamics of the initial and final states using metadynamics and free energy perturbation, respectively. While the initial results of the two approaches are found to be completely different, it is demonstrated that the thermodynamics can be reconciled with appropriate corrections for the standard states, thus illustrating the need for caution in interpreting raw free energy curves for ion binding as widely found in the literature. In addition, a new efficient approach is introduced to correct for the system size dependence of kink site energies both due to the periodic interaction of charges in an inhomogeneous dielectric system and due to the dipolar interactions between pairs of kinks along a row. Ultimately, it is shown that with suitable care, both classes of free energy techniques are capable of producing kink site stabilities that are consistent with the solubility of the underlying bulk solid. However, the precise values for individual kink sites exhibit a small systematic offset, which can be ascribed to the contribution of the interfacial potential to the pathway-based results. For the case of NaCl, the free energies of the kink sites relative to a 1 M aqueous solution for Na⁺ and Cl^- are found to be surprisingly different and of opposite sign, despite the ions having very similar hydration free energies.

Effect of sodium chloride adsorption on the surface premelting of ice

We characterise the structural properties of the quasi-liquid layer (QLL) at two low-index ice surfaces in the presence of sodium chloride (Na⁺/Cl^-) ions by molecular dynamics simulations. We find that the presence of a high surface density of Na⁺/Cl^- pairs changes the surface melting behaviour from step-wise to gradual melting. The ions lead to an overall increase of the thickness and the disorder of the QLL, and to a low-temperature roughening transition of the air-ice interface. The local molecular structure of the QLL is similar to that of liquid water, and the differences between the basal and primary prismatic surface are attenuated by the presence of Na⁺/Cl^- pairs. These changes modify the crystal growth rates of different facets and the solvation environment at the surface of sea-water ice with a potential impact on light scattering and environmental chemical reactions.

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.

The liquidus temperature curve of aqueous methanol mixtures: a numerical simulation study

The liquidus temperature curve that characterizes the boundary between the liquid methanol/water mixture and its coexistence with ice Ih is determined using the direct-coexistence method. Several methanol concentrations and pressures of 0.1 MPa, 50 MPa, and 100 MPa are considered. In this study, we used the TIP4P/Ice model for water and two different models for methanol: OPLS and OPLS/2016, using the geometric rule for the Lennard-Jones cross interactions. We compared our simulation results with available experimental data and found that this combination of models reproduces reasonably well the liquidus curve for methanol mole fractions up to xm=0.3 at p=0.1 MPa. The freezing point depression of these mixtures is calculated and compared to experimental results. We also analyzed the effect of pressure on the liquidus curve, and we found that both models also reproduce qualitatively well the experimental decreasing of the liquidus temperatures as the pressure increases.

Accurate prediction of ice nucleation from room temperature water

Crystal nucleation is one of the most fundamental processes in the physical sciences and almost always occurs heterogeneously with the aid of a nucleating substrate. No example of nucleation is more ubiquitous and impactful than the formation of ice, vital to fields as diverse as geology, biology, aeronautics, and climate science. However, despite considerable effort, we still cannot predict a priori the efficacy of a nucleating agent. Here we utilize deep learning methods to accurately predict nucleation ability from images of room temperature liquid water-generated from molecular dynamics simulations-on a broad range of substrates. The resulting model, named IcePic, can rapidly and accurately infer nucleation ability, eliminating the requirement for either notoriously expensive simulations or direct experimental measurement. In an online poll, IcePic was found to significantly outperform humans in predicting the ice nucleating efficacy of materials. By analyzing the typical errors made by humans, as well as the application of reverse interpretation methods, physical insights into the role the water contact layer plays in ice nucleation have been obtained. Moving forward, we suggest that IcePic can be used as an easy, cheap, and rapid way to discern the nucleation ability of substrates, also with potential for learning other properties related to interfacial water.

Maximum in density of electrolyte solutions: Learning about ion-water interactions and testing the Madrid-2019 force field

In this work, we studied the effect of Li⁺, Na⁺, K⁺, Mg²⁺, and Ca²⁺ chlorides and sulfates on the temperature of maximum density (TMD) of aqueous solutions at room pressure. Experiments at 1 molal salt concentration were carried out to determine the TMD of these solutions. We also performed molecular dynamics simulations to estimate the TMD at 1 and 2 m with the Madrid-2019 force field, which uses the TIP4P/2005 water model and scaled charges for the ions, finding an excellent agreement between experiment and simulation. All the salts studied in this work shift the TMD of the solution to lower temperatures and flatten the density vs temperature curves (when compared to pure water) with increasing salt concentration. The shift in the TMD depends strongly on the nature of the electrolyte. In order to explore this dependence, we have evaluated the contribution of each ion to the shift in the TMD concluding that Na⁺, Ca²⁺, and SO₄ ^2- seem to induce the largest changes among the studied ions. The volume of the system has been analyzed for salts with the same anion and different cations. These curves provide insight into the effect of different ions upon the structure of water. We claim that the TMD of electrolyte solutions entails interesting physics regarding ion-water and water-water interactions and should, therefore, be considered as a test property when developing force fields for electrolytes. This matter has been rather unnoticed for almost a century now and we believe it is time to revisit it.

Atomistic simulation framework for molten salt vapor–liquid equilibrium prediction and its application to NaCl

Knowledge of the vapor–liquid equilibrium (VLE) properties of molten salts is important in the design of thermal energy storage systems for solar power and nuclear energy production applications. The high temperatures involved make their experimental determination problematic, and the development of both macroscopic thermodynamic correlations and predictive molecular-based methodologies are complicated by the requirement to appropriately incorporate the chemically reacting vapor-phase species. We derive a general thermodynamic-based atomistic simulation framework for molten salt VLE prediction and show its application to NaCl. Its input quantities are temperature-dependent ideal-gas free energy data for the vapor phase reactions and density and residual chemical potential data for the liquid. If these are not available experimentally, the former may be predicted using standard electronic structure software, and the latter may be predicted by means of classical atomistic simulation methodology. The framework predicts the temperature dependence of vapor pressure, coexisting phase densities, vapor phase composition, and vaporization enthalpy. It also predicts the concentrations of vapor phase species present in minor amounts (such as the free ions), quantities that are extremely difficult to measure experimentally. We furthermore use the results to obtain an approximation to the complete VLE binodal dome and the critical properties. We verify the framework for molten NaCl, for which experimentally based density and chemical potential data are available in the literature. We then apply it to the analysis of NaCl simulation data for two commonly used atomistic force fields. The framework can be readily extended to molten salt mixtures and to ionic liquids.

Freezing point depression of salt aqueous solutions using the Madrid-2019 model

Salt aqueous solutions are relevant in many fields, ranging from biological systems to seawater. Thus, the availability of a force-field that is able to reproduce the thermodynamic and dynamic behavior of salt aqueous solutions would be of great interest. Unfortunately, this has been proven challenging, and most of the existing force-fields fail to reproduce much of their behavior. In particular, the diffusion of water or the salt solubility are often not well reproduced by most of the existing force-fields. Recently, the Madrid-2019 model was proposed, and it was shown that this force-field, which uses the TIP4P/2005 model for water and non-integer charges for the ions, provides a good description of a large number of properties, including the solution densities, viscosities, and the diffusion of water. In this work, we assess the performance of this force-field on the evaluation of the freezing point depression. Although the freezing point depression is a colligative property that at low salt concentrations depends solely on properties of pure water, a good model for the electrolytes is needed to accurately predict the freezing point depression at moderate and high salt concentrations. The coexistence line between ice and several salt aqueous solutions (NaCl, KCl, LiCl, MgCl₂, and Li₂SO₄) up to the eutectic point is estimated from direct coexistence molecular dynamics simulations. Our results show that this force-field reproduces fairly well the experimentally measured freezing point depression with respect to pure water freezing for all the salts and at all the compositions considered.

Anomalous Phase Behaviors of Monolayer NaCl Aqueous Solutions Induced by Effective Coulombic Interactions within Angstrom-Scale Slits

Interests in subnanofluidic devices have called for molecular dynamics (MD) simulation studies of the thermodynamic behavior of monolayer salt solution within angstrom-scale slits. However, it still remains a grand challenge to accurately describe the Coulombic interactions by incorporating the effects of charge transfer and electronic dielectric screening. Herein, by using the electronic continuum model, where the effective ion charges are fine-tuned with a scaling factor of λ, we present simulation evidence that the effective Coulombic interactions among Na⁺/Cl^- ions can strongly affect the behavior of monolayer ionic aqueous solution. Our microsecond-scale MD simulations show that only the counterions with moderate effective charges (0.3 ≤ λ ≤ 0.8) can dissolve in monolayer water, whereas the high effective charges (λ ≥ 0.85) induce ions to assemble into monolayer nanocrystals, and ions with the low effective charges (λ ≤ 0.2) exhibit gas-like nanobubble. These findings could provide deeper insights into the physical chemistry behind subnanofluidic iontronic devices.

Effects of Externally Applied Electric Fields on the Manipulation of Solvated-Chignolin Folding: Static- versus Alternating-Field Dichotomy at Play

The interaction between a protein and external electric field (EF) can alter its structure and dynamical behavior, which has a potential impact on the biological function of proteins and cause uncertain health consequences. Conversely, the application of EFs of judiciously selected intensity and frequency can help to treat disease, and optimization of this requires a greater understanding of EF-induced effects underpinning basic protein biophysics. In the present study, chignolin─an artificial protein sufficiently small to undergo fast-folding events and transitions─was selected as an ideal prototype to investigate how, and to what extent, externally applied electric fields may manipulate or influence protein-folding phenomena. Nonequilibrium molecular dynamics (NEMD) simulations have been performed of solvated chignolin to determine the distribution of folding states and their underlying transition dynamics, in the absence and presence of externally applied electric fields (both static and alternating); a key focus has been to ascertain how folding pathways are altered in an athermal sense by external fields. Compared to zero-field conditions, a dramatically different─indeed, bifurcated─behavior of chignolin-folding processes emerges between static- and alternating-field scenarios, especially vis-à-vis incipient stages of hydrophobic-core formation: in alternating fields, fold-state populations diversified, with an attendant acceleration of state-hopping folding kinetics, featuring the concomitant emergence of a new, quasi-stable structure compared to the native structure, in field-shifted energy landscapes.

Homogeneous nucleation of NaCl in supersaturated solutions

The seeding method is an approximate approach to investigate nucleation that combines molecular dynamics simulations with classical nucleation theory. Recently, this technique has been successfully implemented in a broad range of nucleation studies. However, its accuracy is subject to the arbitrary choice of the order parameter threshold used to distinguish liquid-like from solid-like molecules. We revisit here the crystallization of NaCl from a supersaturated brine solution and show that consistency between seeding and rigorous methods, like Forward Flux Sampling (from previous work) or spontaneous crystallization (from this work), is achieved by following a mislabelling criterion to select such threshold (i.e. equaling the fraction of the mislabelled particles in the bulk parent and nucleating phases). This work supports the use of seeding to obtain fast and reasonably accurate nucleation rate estimates and the mislabelling criterion as one giving the relevant cluster size for classical nucleation theory in crystallization studies.

Structure and dynamics of aqueous NaCl solutions at high temperatures and pressures

The structure of a concentrated solution of NaCl in D₂O was investigated by in situ high-pressure neutron diffraction with chlorine isotope substitution to give site-specific information on the coordination environment of the chloride ion. A broad range of densities was explored by first increasing the temperature from 323 to 423 K at 0.1 kbar and then increasing the pressure from 0.1 to 33.8 kbar at 423 K, thus mapping a cyclic variation in the static dielectric constant of the pure solvent. The experimental work was complemented by molecular dynamics simulations using the TIP4P/2005 model for water, which were validated against the measured equation of state and diffraction results. Pressure-induced anion ordering is observed, which is accompanied by a dramatic increase in the Cl-O and O-O coordination numbers. With the aid of bond-distance resolved bond-angle maps, it is found that the increased coordination numbers do not originate from a sizable alteration to the number of either Cl⋯D-O or O⋯D-O hydrogen bonds but from the appearance of non-hydrogen-bonded configurations. Increased pressure leads to a marked decrease in the self-diffusion coefficients but has only a moderate effect on the ion-water residence times. Contact ion pairs are observed under all conditions, mostly in the form of charge-neutral NaCl⁰ units, and coexist with solvent-separated Na⁺-Na⁺ and Cl^--Cl^- ion pairs. The exchange of water molecules with Na⁺ adopts a concerted mechanism under ambient conditions but becomes non-concerted as the state conditions are changed. Our findings are important for understanding the role of extreme conditions in geochemical processes.

Spontaneous NaCl-doped ices Ih, Ic, III, V and VI. Understanding the mechanism of ion inclusion and its dependence on the crystalline structure of ice

Direct coexistence simulations on a microsecond time scale have been performed for different types of ice (I_h, I_c, III, V, and VI) in contact with a NaCl aqueous solution at different pressures. In line with the previous results obtained for ice I_h [Conde et al., Phys. Chem. Chem. Phys., 2017, 19, 9566-9574], our results reveal the spontaneous growth of a new ice doped phase and the formation of a brine rejection phase in all ices studied. However, both the preferential incorporation of ions into the ice lattice and the inclusion mechanisms depend on the crystalline structure of each ice. This work shows the inclusion of Cl^- and Na⁺ ions in ice from salt using molecular dynamics simulation, in agreement with the experimental evidence found in the literature. The model used for water is TIP4P/2005. For NaCl we employ a set of potential parameters that uses unit charges for the ions.

Structural and dynamic properties of some aqueous salt solutions

Aqueous salt solutions are utilized and encountered in wide-ranging technological applications and natural settings. Towards improved understanding of the effect of salts on the dynamic properties of such systems, dilute aqueous salt solutions (up to 1 molar concentration) are investigated here, via experiments and molecular simulations. Four salts are considered: sodium chloride, for which published results are readily available for comparison, ammonium acetate, barium acetate and barium nitrate, for which published data are scarce. In the present work, molecular dynamics (MD) simulations are conducted to quantify viscosity and water self-diffusion coefficients, together with rheometry and Pulsed Field Gradient Spin Echo (PFGSE)-NMR experiments for validation. Simulation predictions are consistent with experimental observations in terms of trend and magnitude of salt-specific effects. Combining insights from the approaches considered, an interpretation of the results is proposed whereby the capacity of salts to influence bulk dynamics arises from their molecular interfacial area and strength of interaction with first hydration-shell water molecules. For the concentration range investigated, the interpretation could be useful in formulating aqueous systems for applications including the manufacturing of advanced catalysts.

Role of water model on ion dissociation at ambient conditions

We study ion pair dissociation in water at ambient conditions using a combination of classical and ab initio approaches. The goal of this study is to disentangle the sources of discrepancy observed in computed potentials of mean force. In particular, we aim to understand why some models favor the stability of solvent-separated ion pairs vs contact ion pairs. We found that some observed differences can be explained by non-converged simulation parameters. However, we also unveil that for some models, small changes in the solution density can have significant effects on modifying the equilibrium balance between the two configurations. We conclude that the thermodynamic stability of contact and solvent-separated ion pairs is very sensitive to the dielectric properties of the underlying simulation model. In general, classical models are very robust in providing a similar estimation of the contact ion pair stability, while this is much more variable in density functional theory-based models. The barrier to transition from the solvent-separated to contact ion pair is fundamentally dependent on the balance between electrostatic potential energy and entropy. This reflects the importance of water intra- and inter-molecular polarizability in obtaining an accurate description of the screened ion-ion interactions.

Monolayer Structures of Supramolecular Antagonistic Salt Aggregates

The speculated presence of monomolecular lamellae of antagonistic salts in oil-water mixtures has left several open questions besides their hypothetical existence, including their microscopic structure and stabilization mechanism. Here, we simulate the spontaneous formation of supramolecular aggregates of the antagonistic salt sodium tetraphenylborate (NaBPh₄) in water and 3-methylpyridine (3-MP) at the atomistic level. We show that, indeed, the lamellae are formed by a monomolecular layer of the anion, enveloped by 3-MP and hydrated sodium counterions. To understand which thermodynamic forces drive the aggregation, we compare the full-atomistic model with a simplified one for the salt and show that the strong hydrophobic effect granted by the large excluded volume of the anion, together with electrostatic repulsion, suffice to explain the stability of the monomolecular lamellae.

"In Silico" Seawater

Many important processes affecting the earth's climate are determined by the physical properties of seawater. In addition, desalination of seawater is a significant source of drinking water for the human population living in coastal areas. Since the physical properties of seawater governing these processes depend on the molecular interactions among its components, a deeper knowledge of seawater at the molecular level would contribute to a better understanding of these phenomena. However, in strong contrast with the situation in other areas such as biomolecules or materials science, molecular simulation studies reporting the physical properties of seawater are currently lacking. This is probably due to the usual perception of the seawater composition being too complex to approach. This point of view ignores the fact that physical properties of seawater are dependent on a single parameter representing the composition, namely the salinity. This is because the relative proportions of any two major constituents of seasalt are always the same. Another obstacle to performing molecular simulations of seawater could have been the unavailability of a satisfactory force field representing the interactions between water molecules and dissolved substances. However, this drawback has recently been overcome with the proposal of the Madrid-2019 force field. In this work we show for the first time that molecular simulation of seawater is feasible. We have performed molecular dynamics simulations of a system, the composition of which is close to the average composition of standard seawater and with the molecular interactions given by the Madrid-2019 force field. In this way we are able to provide quantitative or semiquantitative predictions for a number of relevant physical properties of seawater for temperatures and salinities from the oceanographic range to those relevant to desalination processes. The computed magnitudes include static (density), dynamical (viscosity and diffusion coefficients), structural (ionic hydration, ion-ion distribution functions), and interfacial (surface tension) properties.

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.

Dynamic and Static Properties of Aqueous NaCl Solutions at 25°C as a Function of NaCl Concentration: A Molecular Dynamics Simulation Study

We present the result of molecular dynamics (MD) simulations to calculate the molar conductivity ${\mathrm{\Lambda }}_m\left(={\lambda }_{N{a}^{+}}+{\lambda }_{C{l}^{-}}\right)$ of NaCl in SPC/E water at 25°C as a function of NaCl concentration (c) using Ewald sums employing a velocity Verlet algorithm. It is found that the MD result for Λm with Ewald sum parameter κ = 0.10 Å−1 gives the closest one to the experimental data and that the obtained radial distribution functions $g_{ii}$ (r) with κ = 0.10 Å−1 show a dramatic change with a very deep minimum of $g_{\text{NaCl}}$ (r) and, as a result, sharp maxima of $g_{\text{NaNa}}$ (r) and $g_{\text{ClCl}}$ (r) at the distance 9.95 Å, which indicates a characteristic of ionic atmosphere, the basis of the Debye–Hückel theory of ionic solutions. The static and dynamic properties of NaCl (aq) solutions are analyzed in terms of radial distribution functions, hydration numbers, coordination numbers around Na+ and Cl−, residence times of water around Na+ and Cl−, water diffusion, and ion-ion electrostatic energies to explain the behavior of the molar conductivity Λm of NaCl obtained from our MD simulations.

Effect of temperature on the coupling transport of water and ions through a carbon nanotube in an electric field

Temperature governs the motion of molecules at the nanoscale and thus should play an essential role in determining the transport of water and ions through a nanochannel, which is still poorly understood. This work devotes to revealing the temperature effect on the coupling transport of water and ions through a carbon nanotube by molecular dynamics simulations. A fascinating finding is that the ion flux order changes from cation > anion to anion > cation with the increase in field strength, leading to the same direction change of water flux. The competition between ion hydration strength and mobility should be a partial reason for this ion flux order transition. High temperatures significantly promote the transport of water and ions, stabilize the water flux direction, and enhance the critical field strength. The ion translocation time exhibits an excellent Arrhenius relation with the temperature and a power law relation with the field strength, yielding to the Langevin dynamics. However, because of self-diffusion, the water translocation time displays different behaviors without following the ions. The high temperature also leads to an abnormal maximum behavior of the ion flux, deciphered by the massive increase in water flow that inversely hinders the ion flux, suggesting the coexistence of water-ion coupling transport and competition. Our results shed deep light on the temperature dependence of coupling transport of water and ions, answering a fundamental question on the water flux direction during the ionic transport, and thus should have great implications in the design of high flux nanofluidic devices.

Pressure in Molecular Simulations with Scaled Charges. 1. Ionic Systems

Charge scaling, rationalized as MDEC (molecular dynamics in electronic continuum) or ECC (electronic continuum correction), has become a widely used simple approach to how to avoid self-consistent induced dipoles yet approximately take into account the effects of electronic polarizability. It has been assumed that the continuum permittivity does not depend on density; in turn, pressure is calculated by standard formulas. In this work, we elaborate a complementary approximation of density-independent molecular polarizability and derive formulas for pressure corrections within the MDEC framework; real behavior lies between these two extremes. The pressure corrections for test ionic systems are huge and negative, leading to sizable densities in constant-pressure MDEC simulations. A comparison of MDEC results with equivalent polarizable systems gives a good pressure match for a crystal but very low MDEC pressures for ionic liquids. These results witness about the importance of a correct density dependence not only of continuum permittivity in MDEC simulations but also of polarizability in polarizable simulations.

Molecular simulation of liquid-vapor coexistence for NaCl: Full-charge vs scaled-charge interaction models

Scaled-charge models have been recently introduced for molecular simulations of electrolyte solutions and molten salts to attempt to implicitly represent polarizability. Although these models have been found to accurately predict electrolyte solution dynamic properties, they have not been tested for coexistence properties, such as the vapor pressure of the melt. In this work, we evaluate the vapor pressure of a scaled-charge sodium chloride (NaCl) force field and compare the results against experiments and a non-polarizable full-charge force field. The scaled-charge force field predicts a higher vapor pressure than found in experiments, due to its overprediction of the liquid-phase chemical potential. Reanalyzing the trajectories generated from the scaled-charge model with full charges improves the estimation of the liquid-phase chemical potential but not the vapor pressure.

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.

Effect of dissolved salt on the anomalies of water at negative pressure

Adding salt to water at ambient pressure affects its thermodynamic properties. At low salt concentration, anomalies such as the density maximum are shifted to lower temperature, while at large enough salt concentration, they cannot be observed any more. Here, we investigate the effect of salt on an anomaly recently observed in pure water at negative pressure: the existence of a sound velocity minimum along isochores. We compare experiments and simulations for an aqueous solution of sodium chloride with molality around 1.2 mol kg^-1, reaching pressures beyond -100 MPa. We also discuss the origin of the minima in the sound velocity and emphasize the importance of the relative position of the temperatures of sound velocity and density anomalies.

Hierarchical phenomena in multicomponent liquids: simulation methods, analysis, chemistry

Complex, multicomponent, solutions have often been studied solely through the lens of specific applications of interest. Yet advances to both simulation methodologies (enhanced sampling, etc.) and analysis techniques (network analysis algorithms and others), are creating a trove of data that reveal transcending characteristics across vast compositional phase space. This perspective discusses technical considerations of the reliable and accurate simulations of complex solutions, followed by the advances to analysis algorithms that elucidate coupling of different length and timescale behavior (hierarchical phenomena). The different manifestations of hierarchical phenomena are presented across an array of solution environments, emphasizing fundamental and ongoing science questions. With a more advanced molecular understanding in hand, a quintessential application (solvent extraction) is discussed, where significant opportunities exist to re-imagine the technical scope of an established technology.

Dielectric constant and density of aqueous alkali halide solutions by molecular dynamics: A force field assessment

The concentration dependence of the dielectric constant and the density of 11 aqueous alkali halide solutions (LiCl, NaCl, KCl, RbCl, CsCl, LiI, NaI, KI, CsI, KF, and CsF) is investigated by molecular simulation. Predictions using eight non-polarizable ion force fields combined with the TIP4P/ε water model are compared to experimental data. The influence of the water model and the temperature on the results for the NaCl brine are also addressed. The TIP4P/ε water model improves the accuracy of dielectric constant predictions compared to the SPC/E water model. The solution density is predicted well by most ion models. Almost all ion force fields qualitatively capture the decline of the dielectric constant with the increase of concentration for all solutions and with the increase of temperature for NaCl brine. However, the sampled dielectric constant is mostly in poor quantitative agreement with experimental data. These results are related to the microscopic solution structure, ion pairing, and ultimately the force field parameters. Ion force fields with excessive contact ion pairing and precipitation below the experimental solubility limit generally yield higher dielectric constant values. An adequate reproduction of the experimental solubility limit should therefore be a prerequisite for further investigations of the dielectric constant of aqueous electrolyte solutions by molecular simulation.