Computer simulation studies of aqueous sodium chloride solutions at 298 K and 683 K

Cations Do Not Alter the Membrane Structure of POPC—A Lipid With an Intermediate Area

Combining small-angle neutron scattering (SANS), small-angle X-ray scattering (SAXS), and densitometric measurements, we have studied the interactions of the divalent cations Ca²⁺ and Mg²⁺ with the lipid vesicles prepared of a mixed-chain palmitoyl-oleoyl-phosphatidylcholine (POPC) at 25°C. The structural parameters of the POPC bilayer, such as the bilayer thickness, lateral area, and volume per lipid, displayed no changes upon the ion addition at concentrations up to 30 mM and minor changes at > 30 mM Ca²⁺ and Mg²⁺, while some decrease in the vesicle radius was observed over the entire concentration range studied. This examination allows us to validate the concept of lipid–ion interactions governed by the area per lipid suggested previously and to propose the mixed mode of those interactions that emerge in the POPC vesicles. We speculate that the average area per POPC lipid that corresponds to the cutoff length of lipid–ion interactions generates an equal but opposite impact on ion bridges and separate lipid–ion pairs. As a result of the dynamic equilibrium, the overall structural properties of bilayers are not affected. As the molecular mechanism proposed is affected by the structural properties of a particular lipid, it might help us to understand the fundamentals of processes occurring in complex multicomponent membrane systems.

Precipitation Behavior of Salts in Supercritical Water: Experiments and Molecular Dynamics Simulations

Supercritical water desalination (SCWD) shows great potential in the treatment of high-salt wastewater with zero liquid discharge. To investigate the salt precipitation behavior and mechanism in supercritical water, experiments and molecular dynamics simulations (MDs) were used to study the salting-out process of different salts in supercritical water. The equilibrium concentrations of NaCl, KCl, CaCl2, Na2SO4, and Na2CO3 in supercritical water were experimentally measured. When the temperature exceeded 693 K, the salt equilibrium concentration measured in the experiment was less than 130 mg/L. The solubility decreased in the order of KCl > NaCl > CaCl2 > Na2SO4 > Na2CO3. To elucidate the effects of different cations and anions in supercritical water on salt dissolution and precipitation behavior, the potential energy, radial distribution function (RDF) and coordination number in the system were obtained via molecular dynamics simulation. Experimental and MD results showed that salt solubility has significant positive correlation with systemic potential energy and hydration number. MD results indicated that a small ionic radius, large ionic charge, and low hydration coordination number are favorable for inorganic salts to precipitate and crystallize since these factors can strengthen the interaction between free ions and salt clusters. Moreover, due to the formation of multilayer coordination structure, polyatomic ions can achieve a lower equilibrium concentration than that of the corresponding monatomic ions.

How Water-Ion Interactions Control the Formation of Hydrated Electron:Sodium Cation Contact Pairs

Although solvated electrons are a perennial subject of interest, relatively little attention has been paid to the way they behave in aqueous electrolytes. Experimentally, it is known that the hydrated electron's (e_aq^-) absorption spectrum shifts to the blue in the presence of salts, and the magnitude of the shift depends on the ion concentration and the identities of both the cation and anion. Does the blue-shift result from some type of dielectric effect from the bulk electrolyte, or are there specific interactions between the hydrated electron and ions in solution? Previous work has suggested that e_aq^- forms contact pairs with aqueous ions such as Na⁺, leading to the question of what controls the stability of such contact pairs and their possible connection to the observed spectroscopy. In this work, we use mixed quantum/classical simulations to examine the nature of Na⁺:e^- contact pairs in water, using a novel method for quantum umbrella sampling to construct e_aq^--ion potentials of mean force (PMF). We find that the nature of the contact pair PMF depends sensitively on the choice of the classical interactions used to describe the Na⁺-water interactions. When the ion-water interactions are slightly stronger, the corresponding cation:e^- contact pairs form at longer distances and become free energetically less stable. We show that this is because there is a delicate balance between solvation of the cation, solvation of e_aq^- and the direct electronic interaction between the cation and the electron, so that small changes in this balance lead to large changes in the formation and stability of e^--ion contact pairs. In particular, strengthening the ion-water interactions helps to maintain a favorable local solvation environment around Na⁺, which in turn forces water molecules in the first solvation shell of the cation to be unfavorably oriented toward the electron in a contact pair; stronger solvation of the cation also reduces the electronic overlap of e_aq^- with Na⁺. We also find that the calculated spectra of different models of Na⁺:e^- contact pairs do not shift monotonically with cation-electron distance, and that the calculated spectral shifts are about an order of magnitude larger than experiment, suggesting that isolated contact pairs are not the sole explanation for the blue-shift of the hydrated electron's spectrum in the presence of electrolytes.

Effects of Salts on the Solvation of Hydrophobic Objects in Water

The solvation of large, hydrophobic objects in water is facilitated by the formation of a low-density region surrounding the solute that is separated from the bulk liquid by an interface, which has a structure that resembles that between a liquid and its vapor. We study the effect of dissolved sodium chloride on the thermodynamics of solvation and on the solvent structure surrounding hydrophobic solutes in the size regime where this interface is not yet fully formed. Using biased Molecular Dynamics computer simulations, we calculate solvation free energies and orientational distributions of water molecules at different salt concentrations and solute sizes. We find that while the effects of sodium chloride on thermodynamic properties are small, the ions' response to the presence of a hydrophobic solute differs significantly from that of the water. Our findings provide mechanistic insight into how our understanding of hydrophobic solvation in water can be extended to electrolyte solutions.

Failure of molecular dynamics to provide appropriate structures for quantum mechanical description of the aqueous chloride ion charge-transfer-to-solvent ultraviolet spectrum

The lowest band in the charge-transfer-to-solvent ultraviolet absorption spectrum of aqueous chloride ion is studied by experiment and computation. Interestingly, the experiments indicate that at concentrations up to at least 0.25 M, where calculations indicate ion pairing to be significant, there is no notable effect of ionic strength on the spectrum. The experimental spectra are fitted to aid comparison with computations. Classical molecular dynamic simulations are carried out on dilute aqueous Cl-, Na+, and NaCl, producing radial distribution functions in reasonable agreement with experiment and, for NaCl, clear evidence of ion pairing. Clusters are extracted from the simulations for quantum mechanical excited state calculations. Accurate ab initio coupled-cluster benchmark calculations on a small number of representative clusters are carried out and used to identify and validate an efficient protocol based on time-dependent density functional theory. The latter is used to carry out quantum mechanical calculations on thousands of clusters. The resulting computed spectrum is in excellent agreement with experiment for the peak position, with little influence from ion pairing, but is in qualitative disagreement on the width, being only about half as wide. It is concluded that simulation by classical molecular dynamics fails to provide an adequate variety of structures to explain the experimental CTTS spectrum of aqueous Cl-.

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.

An ab initio molecular dynamics study of benzene in water at supercritical conditions: Structure, dynamics, and polarity of hydration shell water and the solute

Anisotropic structure and dynamics of the hydration shell of a benzene solute in supercritical water are investigated by means of ab initio molecular dynamics simulations. The polarity and structural distortion of the benzene solute in supercritical water are also investigated in this study. Calculations are done at 673 K for three different densities of the solvent. The simulations are carried out using the Becke-Lee-Yang-Parr (BLYP) and also the Becke-Lee-Yang-Parr functional including dispersion corrections of Grimme (BYLP-D). The structural anisotropy is found to exist even at supercritical conditions as elucidated by the radial distribution functions of different conical regions and also by angular and spatial distribution functions. The benzene-water πH-bond and also the water-water hydrogen bonds are found to exist even at the supercritical temperature of 673 K. However, the numbers of these hydrogen bonds are reduced substantially with a decrease in water density. The water molecules in the axial region of benzene are found to be preferably oriented with one OH vector pointing toward the benzene ring, whereas the water molecules located in the equatorial region are found to orient their dipoles mostly parallel to the ring plane. The orientational distributions, however, are found to be rather broad at the supercritical temperature due to thermal fluctuations. Although the water molecules have faster dynamics at these supercritical conditions, a slight difference is observed in the dynamics of the solvation shell and bulk molecules. The conformational flexibility of the ring is found to be enhanced which causes an increase in polarity of the benzene solute in water under supercritical conditions.

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.

Morphology and ion diffusion in PEDOT:Tos. A coarse grained molecular dynamics simulation

A Martini coarse-grained Molecular Dynamics (MD) model for the doped conducting polymer poly(3,4-ethylenedioxythiophene) (PEDOT) is developed. It is shown that the diffusion coefficients decrease exponentially as the hydration level is reduced.

Understanding the temperature effect on transport dynamics and structures in polyamide reverse osmosis system via molecular dynamics simulations

Molecular simulations could disclose the transport dynamics, membrane structures and temperature effect on reverse osmosis process.

An Activity-Based Dissolution Model for Solute-Containing Microdroplets

When a solute is present in an aqueous droplet, the water activity in the droplet and the rate of droplet dissolution are both decreased (as compared to a pure water droplet). One of the main parameters that controls this effect is the dynamically changing solute concentration, and therefore water activity and chemical potential, at the droplet interface. This work addresses the importance of understanding how water activity changes during solution droplet dissolution. A model for dissolution rate is presented that accounts for the kinetic effects of changing water activity at the droplet interface during the dissolution of an aqueous salt solution microdroplet into a second immiscible liquid phase. The important underlying question in this model is whether the dissolving component can be considered in local equilibrium on both sides of the droplet interface and whether this assumption is sufficient to account for the kinetics of dissolution. The dissolution model is based on the Epstein-Plesset equation, which has previously been applied to pure gas (bubble) and liquid (droplet) dissolution into liquid phases, but not to salt solutions. The model is tested by using the micropipet technique to form and observe the dehydration of single NaCl solution microdroplets in octanol or butyl acetate. The model successfully predicts the droplet diameter as a function of time in both organic solvents. The NaCl concentration in water is measured well into the supersaturated area >5.4 M, and the supersaturation limit at which NaCl nucleation happens is reported to be 10.24 ± 0.31 M.

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.

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.

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.

Influence of Magnetic Field on Aqueous NaCl Solutions: A Foundational Research on the Desalination Method Based on the Rotating Electromagnetic Effect

In order to investigate whether magnetic field can improve the seawater desalination efficiency, the influence of magnetic field on the aqueous solutions of NaCl is studied by molecular dynamics (MD) simulations. The results show that, under the influence of magnetic field, the hydration number of Cl- ions increases and the mobility of hydrated Cl- ions is weakened, while the mean size of water clusters decreases and the mobility of water molecules is enhanced. These may lead to higher salt rejection and water flux by the usage of magnetic field in our novel desalination method based on rotating electromagnetic effect.

SINGLE-PARTICLE AND PAIR DYNAMICAL PROPERTIES OF ACETONE–METHANOL MIXTURES CONTAINING CHARGED AND NEUTRAL SOLUTES: A MOLECULAR DYNAMICS STUDY

The dynamical properties of acetone–methanol mixtures containing either an ionic or a neutral hydrophobic solute are investigated by means of a series of molecular dynamics simulations. The primary goal has been to study how the solute and solvent dynamical properties change with variation of composition of the mixture ranging from pure acetone to pure methanol. The variations of structure and energetics of the mixture with composition are also calculated. The diffusion coefficients of both ionic and neutral solutes are found to show nonlinear variation with composition of the mixture, although the extent of nonlinearity in the diffusion of the neutral solute is much weaker. Calculations of appropriate solute-solvent distribution functions reveal the extent and nature of selective solvation of these solute species which play a role in determining the nonideal dynamical characteristics of these solutes. The free energies of solvation of the ionic solutes are also calculated and the results are discussed in the context of their dynamical behavior. The hydrogen bond statistics and dynamics of these mixtures are also calculated over their entire composition range. The energies and lifetimes of hydrogen bonds between an acetone and a methanol molecule or between two methanol molecules are found to increase with increase of acetone mole fraction of the mixture. Residence times of methanol molecules in solvation shells of acetone and methanol are also found to follow the same trend as relaxation times. However, these pair dynamical properties show essentially linear dependence on composition, thus behave almost ideally with respect to changes in composition of the mixture.

Aqueous NaCl and CsCl solutions confined in crystalline slit-shaped silica nanopores of varying degree of protonation

All-atom molecular dynamics simulations were conducted to study the dynamics of aqueous electrolyte solutions confined in slit-shaped silica nanopores of various degrees of protonation. Five degrees of protonation were prepared by randomly removing surface hydrogen atoms from fully protonated crystalline silica surfaces. Aqueous electrolyte solutions containing NaCl or CsCl salt were simulated at ambient conditions. In all cases, the ionic concentration was 1 M. The results were quantified in terms of atomic density distributions within the pores, and the self-diffusion coefficient along the direction parallel to the pore surface. We found evidence for ion-specific properties that depend on ion-surface, water-ion, and only in some cases ion-ion correlations. The degree of protonation strongly affects the structure, distribution, and the dynamic behavior of confined water and electrolytes. Cl(-) ions adsorb on the surface at large degrees of protonation, and their behavior does not depend significantly on the cation type (either Na(+) or Cs(+) ions are present in the systems considered). The cations show significant ion-specific behavior. Na(+) ions occupy different positions within the pore as the degree of protonation changes, while Cs(+) ions mainly remain near the pore center at all conditions considered. For a given degree of protonation, the planar self-diffusion coefficient of Cs(+) is always greater than that of Na(+) ions. The results are useful for better understanding transport under confinement, including brine behavior in the subsurface, with important applications such as environmental remediation.

Understanding the role of ion interactions in soluble salt flotation with alkylammonium and alkylsulfate collectors

There is anecdotal evidence for the significant effects of salt ions on the flotation separation of minerals using process water of high salt content. Examples include flotation of soluble salt minerals such as potash, trona and borax in brine solutions using alkylammonium and alkylsulfate collectors such as dodecylamine hydrochloride and sodium dodecylsulfate. Although some of the effects are expected, some do not seem to be encompassed by classical theories of colloid science. Several experimental and modeling techniques for determining solution viscosity, surface tension, bubble-particle attachment time, contact angle, and molecular dynamics simulation have been used to provide further information on air-solution and solid-solution interfacial phenomena, especially with respect to the interfacial water structure due to the presence of dissolved ions. In addition atomic force microscopy, and sum frequency generation vibrational spectroscopy have been used to provide further information on surface states. These studies indicate that the ion specificity effect is the most significant factor influencing flotation in brine solutions.