Molecular dynamics simulations of aqueous NaCl and KCl solutions: Effects of ion concentration on the single-particle, pair, and collective dynamical properties of ions and water molecules

Spatial-Decomposition Analysis of Electrical Conductivity

Spatial decomposition is conducted for the electrical conductivity. The contribution per ion pair at a certain distance is identified in terms of a two-body velocity time correlation function and is integrated over the whole distance of the ion pair to provide the cross-correlation term of the conductivity. The spatial-decomposition formula is an exact expression at any concentrations of ions and incorporates physically appealing pictures in the space domain into the theory of time correlation functions. Illustrative analyses are presented for 1m NaCl aqueous solution and the [C₄ mim][NTf₂ ] ionic liquid. The contrast between the two systems is discussed for the time and spatial ranges of correlations, and it is shown that the ion-pair contribution to the conductivity for the [C₄ mim][NTf₂ ] system is not localized and extends beyond the first coordination shell of the cation-anion pair.

Development of the ReaxFF Methodology for Electrolyte-Water Systems

A new ReaxFF reactive force field has been developed for water-electrolyte systems including cations Li⁺, Na⁺, K⁺, and Cs⁺ and anions F^-, Cl^-, and I^-. The reactive force field parameters have been trained against quantum mechanical (QM) calculations related to water binding energies, hydration energies and energies of proton transfer. The new force field has been validated by applying it to molecular dynamics (MD) simulations of the ionization of different electrolytes in water and comparison of the results with experimental observations and thermodynamics. Radial distribution functions (RDF) determined for most of the atom pairs (cation or anion with oxygen and hydrogen of water) show a good agreement with the RDF values obtained from DFT calculations. On the basis of the applied force field, the ReaxFF simulations have described the diffusion constants for water and electrolyte ions in alkali metal hydroxide and chloride salt solutions as a function of composition and electrolyte concentration. The obtained results open opportunities to advance ReaxFF methodology to a wide range of applications involving electrolyte ions and solutions.

Slippery and Sticky Graphene in Water

Understanding modulation of water molecule slippage along graphene surfaces is crucial for many promising applications of two-dimensional materials. Here, we examine normal and shear forces on supported single-layer graphene using atomic force microscopy and find that the electrolyte composition affects the molecular slippage of nanometer thick films of aqueous electrolytes along the graphene surface. In light of the shear-assisted thermally activated theory, water molecules along the graphene plane are very mobile when subjected to shear. However, upon addition of an electrolyte, the cations can make water stick to graphene, while ion-specific and concentration effects are present. Recognizing the tribological and tribochemical utility of graphene, we also evaluate the impact of this behavior on its frictional response in the presence of water. It appears that the addition of an electrolyte to pure water causes a reduction of the thermal activation energy and of the shear-activation length at several concentrations, both results conversely affecting the friction force. Further, this work can inspire innovation in research areas where changes of the molecular slippage through the modulation of the doping characteristics of graphene in liquid environment can be of use, including molecular sensing, lubrication, and energy storage.

Effects of Temperature and Ionic Concentration on Nanodroplet Electrocoalescence

Droplet electrocoalescence is of interest for various applications such as petroleum dehydration, electrospray ionization, and surface self-cleaning. Here, the effects of temperature and ionic concentration on nanodroplet electrocoalescence are investigated by molecular dynamics simulation. The results show that low ionic concentration rapidly drives ions towards water clusters and leads to dipole polarization of droplets. With an increase of ionic concentration, the particle-particle interaction is enhanced, but the mobility of free water molecules and salt ions is curbed by hydration and ion pairs, which then slows the electrocoalescence. Low temperature accelerates the rotation of water molecules but does not enhance the mobility of ions. Alternatively, high temperature not only breaks the self-assembly of water molecules along the electric field direction but also helps ions to overcome the electrostatic barrier between particles. The latter effect promotes dipole polarization to compensate for the shortcoming of less orientation polarization. The combined effects of ion concentration and temperature are investigated and unified by droplet conductivity from the microscopic point of view. The conductivity increases with the increase in temperatures and ionic concentrations. We confirm that the accurate control of droplet electrocoalescence can be achieved by a suitable combination of temperature and ionic concentration.

Ionic Conductivity, Diffusion Coefficients, and Degree of Dissociation in Lithium Electrolytes, Ionic Liquids, and Hydrogel Polyelectrolytes

The conductive and diffusional behavior of electrolytes in media with different dielectric and viscoelastic properties is investigated. A revised model to separate the contribution of dissociated and nondissociated species to the diffusion coefficients determined with NMR is proposed. Impedance spectroscopy is used to measure the ionic conductivity of lithium salts in aqueous medium, ionic liquids in aprotic solvents, and hydrogel polyelectrolytes. The diffusion coefficients of the species of interest in those systems are determined with multinuclear pulsed-gradient spin-echo (PGSE) NMR. The results are analyzed using the revised model. It is shown that the degree of ionization could be determined directly from measurements of ionic conductivity and diffusion coefficients in very different types of electrolytes and in a wide range of concentrations. Furthermore, these findings support the original Arrhenius hypothesis about electrolytes and show that the assumption of a complete dissociation is not required to describe their conductive behavior. The reduced conductivity observed in hydrogels, at or near swelling equilibrium, compared to that in solutions could be attributed mainly to the hindered ionic mobility caused by the network structure.

Dielectric constant of ionic solutions: Combined effects of correlations and excluded volume

The dielectric constant of ionic solutions is known to reduce with increasing ionic concentrations. However, the origin of this effect has not been thoroughly explored. In this paper, we study two such possible sources: long-range Coulombic correlations and solvent excluded-volume. Correlations originate from fluctuations of the electrostatic potential beyond the mean-field Poisson-Boltzmann theory, evaluated by employing a field-theoretical loop expansion of the free energy. The solvent excluded-volume, on the other hand, stems from the finite ion size, accounted for via a lattice-gas model. We show that both correlations and excluded volume are required in order to capture the important features of the dielectric behavior. For highly polar solvents, such as water, the dielectric constant is given by the product of the solvent volume fraction and a concentration-dependent susceptibility per volume fraction. The available solvent volume decreases as a function of ionic strength due the increasing volume fraction of ions. A similar decrease occurs for the susceptibility due to the correlations between the ions and solvent, reducing the dielectric response even further. Our predictions for the dielectric constant fit well with experiments for a wide range of concentrations for different salts in different temperatures, using a single fit parameter related to the ion size.

Molecular Dynamics Analysis of Synergistic Effects of Ions and Winter Flounder Antifreeze Protein Adjacent to Ice-Solution Surfaces

The control of freezing saline water at the micrometer level has become very important in cryosurgery and cryopreservation of stem cells and foods. Adding antifreeze protein to saline water is a promising method for controlling the freezing because the protein produces a gap between the melting point and the freezing point. Furthermore, a synergistic effect of the solutes occurs in which the freezing point depression of a mixed solution is more noticeable than the sum of two freezing point depressions of single-solute solutions. However, the mechanism of this effect has not yet been clarified. Thus, we have carried out a molecular dynamics simulation on aqueous solutions of winter flounder antifreeze protein and sodium chloride or calcium chloride with an ice layer. The results show that the cations inhibit the hydrogen bond among water molecules not only in the salt solutions but also in the mixed solutions. This inhibition depends on the local number of ions and the valence of cations. The space for water molecules to form the hydrogen bonds becomes small in the case of the mixed solution of the protein and calcium chloride. These findings are consistent with the synergistic effect. In addition, it is found that the diffusion of ions near positively-charged residues is attenuated. This attenuation causes an increase in the possibility of water molecules staying near or inside the hydration shells of the ions. Furthermore, the first hydration shells of the cations become weak in the vicinity of the arginine, lysine and glutamic-acid residues. These factors can be considered to be possible mechanisms of the synergistic effect.

Effects of nanoparticle heating on the structure of a concentrated aqueous salt solution

The effects of a rapidly heated nanoparticle on the structure of a concentrated aqueous salt solution are studied using molecular dynamics simulations. A diamond-like nanoparticle of radius 20 Å is immersed in a sodium-chloride solution at 20% above the experimental saturation concentration and equilibrated at T = 293 K and P = 1 atm. The nanoparticle is then rapidly heated to several thousand degrees Kelvin, and the system is held under isobaric-isoenthalpic conditions. It is observed that after 2-3 ns, the salt ions are depleted far more than water molecules from a proximal zone 15-25 Å from the nanoparticle surface. This leads to a transient reduction in molality in the proximal zone and an increase in ion clustering in the distal zone. At longer times, ions begin to diffuse back into the proximal zone. It is speculated that the formation of proximal and distal zones, and the increase in ion clustering, plays a role in the mechanism of nonphotochemical laser-induced nucleation.

Anomalous cation diffusion in salt-doped confined bilayer ice

The diffusive dynamics of aqueous electrolyte solutions in nanoconfined spaces has attracted considerable attention due to their potential applications in desalination, biosensors and supercapacitors. Here we show by molecular dynamics simulations that lithium and sodium ions diffuse at a rate at least an order of magnitude higher than that of water molecules when the ions are trapped in an ice bilayer confined between two parallel plates. This novel picture is in sharp contrast to the prevailing view that the diffusion rate of ions is comparable to or even lower than that of water in both bulk and confined solutions. The predicted high ion mobility stems from frequent lateral hopping of ions along the coordination sites inside the hydrogen-bonding network connecting the two water layers of the ice bilayer. This anomalous diffusion should provide new insights into the physics of confined aqueous electrolytes.

Molecular Insight into Water Desalination across Multilayer Graphene Oxide Membranes

Transport of ionic solutions through graphene oxide (GO) membranes is a complicated issue because the complex and tortuous structure inside makes it very hard to clarify. Using molecular dynamics (MD) simulations, we investigated the mechanism of water transport and ion movement across multilayer GO. The significant flow rate is considerably influenced by the structural parameters of GO membranes. Because of the size effect on a shrunken real flow area, there is disagreement between the classical continuum model and nanoscaled flow. To eliminate the variance, we obtained modified geometrical parameters from density analysis and used them in the developed hydrodynamic model to give a precise depiction of water flow. Four kinds of solutions (i.e., NaCl, KCl, MgCl₂, and CaCl₂) and different configurational GO sheets were considered to clarify the influence on salt permeation. It is found that the abilities of permeation to ions are not totally up to the hydration radius. Even though the ionic hydration shell is greater than the opening space, the ions can also pass through the split because of the special double-deck hydration structure. In the structure of GO, a smaller layer separation with greater offsetting gaps could substantially enhance the membrane's ability to reject salt. This work establishes molecular insight into the effects of configurational structures and salt species on desalination performance, providing useful guidelines for the design of multilayer GO membranes.

Molecular Simulation of the Ions Ultraconfined in the Nanometer-Channel of Calcium Silicate Hydrate: Hydration Mechanism, Dynamic Properties, and Influence on the Cohesive Strength

Reactive force field molecular dynamics was utilized to investigate the structure, dynamics, and mechanical nature of different cations solvated in the nanometer-channel of highly disordered calcium silicate hydrate. The local structures of different cations bonded with hydroxyl groups are characterized by the long spatial correlation, bond angel distribution preference, and featured coordinated number, resembling those of the tetra-/penta-/octahedron for cation-oxygen structure in the defective region of the silicate glass. Al atoms in the interlayer region play a role in bridging the defective silicate chains and enhance the connectivity of the silicate skeleton. Dynamically, the mobility of ultraconfined water molecules and cations is significantly influenced by the ionic chemistry: the residence time for water molecules in the hydration shell of Al and Mg ions is longer than that in the environment of Na and Ca ions. Furthermore, uniaxial tension simulation provides insight that while both the stiffness and cohesive strength of the C-S-H gels are significantly improved due to the silicate-aluminate branch structure formation, sodium ions with unstable Na-O connection weaken the loading resistance of the C-S-H gels. During the tensile process, the hydrolytic reaction is also affected by the cationic type: water molecules coordinated with Al and Mg cations at high stress state are likely to decompose, but those aggregated with sodium ions are hard to be stretched broken due to the low failure stress.

Temperature dependence of ion diffusion coefficients in NaCl electrolyte confined within graphene nanochannels

The ion diffusion coefficients of NaCl electrolyte confined within graphene nanochannels at different temperatures are investigated using molecular dynamics simulations.

Dielectric behavior for saline solutions from renormalized diagrammatically proper interaction site model theory

We consider the dielectric response of angularly dependent site-site theories for models of aqueous saline solutions. We find that we can use relatively low order approximations of the angularly dependent correlation functions with correct long ranged behavior to obtain good estimates of the dielectric constant for three site water models and simple 1-1 salts. We find that the solution thermodynamics results for this level of theory, as measured by the Kirkwood G integrals and the excess chemical potentials, are in good quantitative agreement with simulation even when the details of the short ranged structure is not as accurately determined. We find that the dielectric constant predictions of both the pure fluid and the salt-water mixtures are similarly predictive, in comparison to both simulation and experiment.

Consensus on the solubility of NaCl in water from computer simulations using the chemical potential route

The solubility of NaCl in water is evaluated by using three force field models: Joung-Cheatham for NaCl dissolved in two different water models (SPC/E and TIP4P/2005) and Smith Dang NaCl model in SPC/E water. The methodology based on free-energy calculations [E. Sanz and C. Vega, J. Chem. Phys. 126, 014507 (2007)] and [J. L. Aragones et al., J. Chem. Phys. 136, 244508 (2012)] has been used, except, that all calculations for the NaCl in solution were obtained by using molecular dynamics simulations with the GROMACS package instead of homemade MC programs. We have explored new lower molalities and made longer runs to improve the accuracy of the calculations. Exploring the low molality region allowed us to obtain an analytical expression for the chemical potential of the ions in solution as a function of molality valid for a wider range of molalities, including the infinite dilute case. These new results are in better agreement with recent estimations of the solubility obtained with other methodologies. Besides, two empirical simple rules have been obtained to have a rough estimate of the solubility of a certain model, by analyzing the ionic pairs formation as a function of molality and/or by calculating the difference between the NaCl solid chemical potential and the standard chemical potential of the salt in solution.

Formulation Studies During Preclinical Development of Influenza Hemagglutinin and Virus-Like Particle Vaccine Candidates

A critical element of vaccine formulation studies is the identification of chemical and physical degradation pathways that compromise structural integrity, and which may in turn affect the clinical safety and efficacy, of macromolecular antigens. Formulation development helps optimize and maintain the long-term storage stability and viability of vaccine antigens in pharmaceutically relevant dosage forms. The protocols presented in this manuscript highlight the use of accelerated stability studies for the formulation of influenza vaccine candidates including virus-like particles (VLP) and particle forming hemagglutinin (HA) antigens. Three case studies, each targeting a different facet of preclinical vaccine formulation development, are reviewed: (1) excipient screening experiments to mitigate VLP physical degradation, (2) methods for monitoring a specific chemical perturbation of the recombinant HA antigen and elucidating its effect on in vitro potency, and (3) maintaining HA conformational stability in the presence of freeze-thaw and freeze-drying stresses.

Passive fluidic diode for simple fluids using nested nanochannel structures

In this paper, we propose a moving part-free fluidic diode for simple fluids using nested nanochannels, which contain inner and outer channels of different lengths. Molecular dynamics simulations show that the fluidic diode accepts water flows in the forward direction and blocks flows in the backward direction in a wide range of pressure drops. The anisotropic flow rates are generated by the distinct activation pressures in different directions. In the forward direction, the activation pressure is low, which is determined by the infiltration pressure of the inner channel. In the backward direction, the activation pressure is quite high due to the capillary effects when flows are released from the inner to the outer channel. The pressure drop range for the fluidic diode can be varied by changing the channel size or surface wettability. The fluidic diode offers an alternative way for flow control in integrated micro- and nanofluidic devices.

Atomistic Molecular Dynamics Simulations of Charged Latex Particle Surfaces in Aqueous Solution

Charged particles in aqueous suspension form an electrical double layer at their surfaces, which plays a key role in suspension properties. For example, binder particles in latex paint remain suspended in the can because of repulsive forces between overlapping double layers. Existing models of the double layer assume sharp interfaces bearing fixed uniform charge, and so cannot describe aqueous binder particle surfaces, which are soft and diffuse, and bear mobile charge from ionic surfactants as well as grafted multivalent oligomers. To treat this industrially important system, we use atomistic molecular dynamics simulations to investigate a structurally realistic model of commercial binder particle surfaces, informed by extensive characterization of particle synthesis and surface properties. We determine the interfacial profiles of polymer, water, bound and free ions, from which the charge density and electrostatic potential can be calculated. We extend the traditional definitions of the inner and outer Helmholtz planes to our diffuse interfaces. Beyond the Stern layer, the simulated electrostatic potential is well described by the Poisson-Boltzmann equation. The potential at the outer Helmholtz plane compares well to the experimental zeta potential. We compare particle surfaces bearing two types of charge groups, ionic surfactant and multivalent oligomers, with and without added salt. Although the bare charge density of a surface bearing multivalent oligomers is much higher than that of a surfactant-bearing surface at realistic coverage, greater counterion condensation leads to similar zeta potentials for the two systems.

Confinement and surface effects of aqueous solutions within charged carbon nanotubes

Size-charge effects of brines in charged nanotubes from a molecular dynamics investigation of ion hydration, water coordination, and hydrogen bonding.

Molecular structure and dynamics of an aqueous sodium chloride solution in nano-pores between portlandite surfaces: a molecular dynamics study

Portlandite plays an important role in the hydration phase of cement-based materials and influences the strength and durability of such materials. This study describes a molecular dynamics study of the structure and dynamics of water and ions confined at ambient temperature in calcium hydroxyl nanopores with widths of 35 Å. Strong layering of water in the vicinity of the (001) surface of portlandite demonstrates special structural features such as large density, good orientation preference, ordered interfacial organization and low diffusion rate. Due to the fixed vibration and rotation of the hydroxyl groups at the interface, water molecules within the first adsorbed layer adopt both H-downward and H-upward orientations by donating H-bonds and accepting H-bonds from the OH groups in the solid surface. Regarding the interaction of the ions and portlandite, Na(+) ions, deeply rooted in spaces in the surface hydroxyl groups, are significantly slowed and remain near the surface for long periods of time. On the other hand, due to the weak H-bonds formed by chloride ions and hydroxyl groups, adsorbed chloride ions near the surface cannot remain for longer times. In addition, when water and ions are confined in the nano-pores, the residence time for the ion-water and ion-ion clusters is lengthened so that the ion adsorption capability of the porlandite surface is enhanced due to the stable Na-Cl connections in the electrolyte solution.

Spontaneous Formation of KCl Aggregates in Biomolecular Simulations: A Force Field Issue?

Realistic all-atom simulation of biological systems requires accurate modeling of both the biomolecules and their ionic environment. Recently, ion nucleation phenomena leading to the rapid growth of KCl or NaCl clusters in the vicinity of biomolecular systems have been reported. To better understand this phenomenon, molecular dynamics simulations of KCl aqueous solutions at three (1.0, 0.25, and 0.10 M) concentrations were performed. Two popular water models (TIP3P and SPC/E) and two Lennard-Jones parameter sets (AMBER and Dang) were combined to produce a total of 80 ns of molecular dynamics trajectories. Results suggest that the use of the Dang cation Lennard-Jones parameters instead of those adopted by the AMBER force-field produces a more accurate description of the ionic solution. In the later case, formation of salt aggregates is probably indicative of an artifact resulting from misbalanced force-field parameters. Because similar results were obtained with two different water parameter sets, the simulations exclude a water model dependency in the formation of anomalous ionic clusters. Overall, the results strongly suggest that for accurate modeling of ions in biomolecular systems, great care should be taken in choosing balanced ionic parameters even when using the most popular force-fields. These results invite a reexamination of older data obtained using available force-fields and a thorough check of the quality of current parameters sets by performing simulations at finite (>0.25 M) instead of minimal salt conditions.

Molecular Insights into Aqueous NaCl Electrolytes Confined within Vertically-oriented Graphenes

Vertically-oriented graphenes (VGs) are promising active materials for electric double layer capacitors (EDLCs) due to their unique morphological and structural features. This study, for the first time, reports the molecular dynamics (MD) simulations on aqueous NaCl electrolytes confined within VG channels with different surface charge densities and channel widths. Simulation results show that the accessibility of ions and the structure of EDLCs are determined by the ion type/size, surface charging, and VG channel width. For relatively narrow VG channels with the same width, the threshold charge density (to compensate the energy penalty for shedding hydration shell) and the dehydration rate of Cl(-) ions are larger than those of Na(+) ions. To achieve the highest ion concentration coefficient, the effective VG channel width should be between the crystal and hydration diameters of the ions. The results are further quantified and elucidated by calculating the electrolyte density profiles. The molecular insights obtained in the current work are useful in guiding the design and fabrication of VGs for advancing their EDLC applications.

Dissecting ion-specific dielectric spectra of sodium-halide solutions into solvation water and ionic contributions

Using extensive equilibrium molecular dynamics simulations we determine the dielectric spectra of aqueous solutions of NaF, NaCl, NaBr, and NaI. The ion-specific and concentration-dependent shifts of the static dielectric constants and the dielectric relaxation times match experimental results very well, which serves as a validation of the classical and non-polarizable ionic force fields used. The purely ionic contribution to the dielectric response is negligible, but determines the conductivity of the salt solutions. The ion-water cross correlation contribution is negative and reduces the total dielectric response by about 5%-10% for 1 M solutions. The dominating water dielectric response is decomposed into different water solvation shells and ion-pair configurations, by this the spectral blue shift and the dielectric decrement of salt solutions with increasing salt concentration is demonstrated to be primarily caused by first-solvation shell water. With rising salt concentration the simulated spectra show more pronounced deviations from a single-Debye form and can be well described by a Cole-Cole fit, in quantitative agreement with experiments. Our spectral decomposition into ionic and different water solvation shell contributions does not render the individual contributions more Debye-like, this suggests the non-Debye-like character of the dielectric spectra of salt solutions not to be due to the superposition of different elementary relaxation processes with different relaxation times. Rather, the non-Debye-like character is likely to be an inherent spectral signature of solvation water around ions.

The effect of concentration- and temperature-dependent dielectric constant on the activity coefficient of NaCl electrolyte solutions

Our implicit-solvent model for the estimation of the excess chemical potential (or, equivalently, the activity coefficient) of electrolytes is based on using a dielectric constant that depends on the thermodynamic state, namely, the temperature and concentration of the electrolyte, ε(c, T). As a consequence, the excess chemical potential is split into two terms corresponding to ion-ion (II) and ion-water (IW) interactions. The II term is obtained from computer simulation using the Primitive Model of electrolytes, while the IW term is estimated from the Born treatment. In our previous work [J. Vincze, M. Valiskó, and D. Boda, "The nonmonotonic concentration dependence of the mean activity coefficient of electrolytes is a result of a balance between solvation and ion-ion correlations," J. Chem. Phys. 133, 154507 (2010)], we showed that the nonmonotonic concentration dependence of the activity coefficient can be reproduced qualitatively with this II+IW model without using any adjustable parameter. The Pauling radii were used in the calculation of the II term, while experimental solvation free energies were used in the calculation of the IW term. In this work, we analyze the effect of the parameters (dielectric constant, ionic radii, solvation free energy) on the concentration and temperature dependence of the mean activity coefficient of NaCl. We conclude that the II+IW model can explain the experimental behavior using a concentration-dependent dielectric constant and that we do not need the artificial concept of "solvated ionic radius" assumed by earlier studies.

Fluid infiltration pressure for hydrophobic nanochannels

We investigate water infiltration pressure for hydrophobic nanochannels through molecular dynamics simulations. It is found that the entrance energy barrier significantly raises the infiltration pressure, which makes the classic Young-Laplace equation invalid for nanochannels. As the channel surface is tuned from superhydrophobic to hydrophobic, the infiltration pressure is greatly reduced mainly due to the decrease of the capillary pressure (Young-Laplace equation) caused by the contact angle change, while the contribution of the entrance energy barrier to the infiltration pressure, which is termed entrance barrier pressure, increases from 25% to 60%.