Systematic comparison of force fields for microscopic simulations of NaCl in aqueous solutions: Diffusion, free energy of hydration, and structural properties

Scaling Solute-Solvent Distances to Improve Solubility and Ion Paring Predictions in Rigid Ion Models

Force fields based on the rigid ion model (RIM) have been developed to accurately predict the various physical and chemical properties of salts and water. However, the combined use of these models often fails to accurately predict the solubility of salts in water. To address this issue, several approaches, such as charge scaling or reparameterization, have been proposed. Nevertheless, these methods require laborious reparameterization of nonbonded force field parameters. In this article, we propose a scaling solute-solvent distance (SSSD) method to improve force fields in predicting salt solubility without changing the solute-solute and solvent-solvent interactions in the original force fields. This method can also tune the ion pairing of salt in water. One main advantage of the SSSD method is that reparameterization of the crystal and water models is not needed. We use two RIMs for the NaCl-water system (JC-SPC/E and SD-SPC/E) and the CHARMM force field for the KCl-water system to demonstrate the improved accuracy in predicting solubility by the SSSD method. Furthermore, we use the RDG-SPC/Fw force field to show that the SSSD method can also be used to tune the ion pairing of CaCO3 in water. Limitations of this method are also discussed.

Unified polarizable electrode models for open and closed circuits: Revisiting the effects of electrode polarization and different circuit conditions on electrode-electrolyte interfaces

A precise understanding of the interfacial structure and dynamics is essential for the optimal design of various electrochemical devices. Herein, we propose a method for classical molecular dynamics simulations to deal with electrochemical interfaces with polarizable electrodes under the open circuit condition. Less attention has been paid to electrochemical circuit conditions in computation despite being often essential for a proper assessment, especially comparison between different models. The present method is based on the chemical potential equalization principle, as is a method developed previously to deal with systems under the closed circuit condition. These two methods can be interconverted through the Legendre transformation, so that the difference in the circuit conditions can be compared on the same footing. Furthermore, the electrode polarization effect can be correctly studied by comparing the present method with the conventional simulations with the electrodes represented by fixed charges, since both of the methods describe systems under the open circuit condition. The method is applied to a parallel-plate capacitor composed of platinum electrodes and an aqueous electrolyte solution. The electrode polarization effects have an impact on the interfacial structure of the electrolyte solution. We found that the difference in the circuit conditions significantly affects the dynamics of the electrolyte solution. The electric field at the charged electrode surface is poorly screened by the nonequilibrium solution structure in the open circuit condition, which accelerates the motion of the electrolyte solution.

Machine Learning Assisted Screening of Two-Dimensional Materials for Water Desalination

There exists a vast expanse of data in the literature which can be harnessed for accelerated design and discovery of advanced materials for various applications of importance ─ for example, desalination of seawater. Here, we develop a machine learning (ML) model, training it with ∼260 molecular dynamics (MD) computation results, to predict the desalination performance of 2D membranes that exist in the literature. The desalination performance variables of water flux and salt rejection rates are correlated to 49 material features related to the chemistry of the pores and the membranes along with applied pressure, salt concentration, partial charges on the atoms, geometry of the pore, the mechanical properties of the membranes, and the properties of water for the water model used. We used the ML model to screen 3814 structurally optimized 2D materials for maximum water flux and salt rejection rates from the literature. We found some candidates that perform ∼4 times better than the more popularly known 2D materials such as graphene and MoS₂. This result is verified using data obtained from MD simulations performed on several representative 2D membranes for different classes. Such validated statistical frameworks using literature data can be very useful in guiding experiments in the field of functional materials for varied applications.

Molecular dynamics simulation study of adsorption of anionic–nonionic surfactants at oil/water interfaces

Four anionic–nonionic surfactants with the same headgroups and different units of oxygen ethyl (EO) and oxygen propyl (PO) were adopted to investigate the influence on oil/water interfacial tensions in this article. Molecular dynamics (MD) simulations were conducted to study the interfacial property of the four surfactants. Four parameters were proposed to reveal the effecting mechanism of molecular structure on interfacial tension, which included the interfacial thickness, order parameter of the hydrophobic chain, radial distribution function, and the solvent accessible surface area. In addition, the electrostatic potential of the four surfactants was calculated. The research results indicated that the interface facial mask formed by the surfactants, which contained three EO or three PO units was more stable, and it was easier for the surfactants of six EO or six PO units to form a microemulsion at higher concentrations. The adsorption mechanism of the anionic–nonionic surfactant systems at the oil/water interfaces was supplemented at a molecular level, which provided fundamental guidance for an in-depth understanding of the optimal selection of the surfactants in enhancing oil recovery.

Four anionic–nonionic surfactants with the same headgroups and different units of oxygen ethyl (EO) and oxygen propyl (PO) were adopted to investigate the influence on oil/water interfacial tensions in this article.

Li+ Selectivity of Carboxylate Graphene Nanopores Inspired by Electric Field and Nanoconfinement

The regulation of the ion selectivity by electric field and ion association on the Li⁺ selectivity of carboxyl functionalized graphene nanopores are investigated by molecular dynamics simulation. Carboxylate graphene nanopores of sub-2 nm exhibit excellent Li⁺ selectivity under the electric field of 1.0 V nm^-1 . The results show that ion association inspired by electric field may be a key factor affecting ion selectivity of sub-2 nm nanopores. The ion association of Mg²⁺ and Cl^- can be promoted obviously near the nanopores under the electric field of 1.0 V nm^-1 . The migrating of Mg²⁺ can be retarded by stable clusters of Mg²⁺ and Cl^- formed near nanopores. The degree of association of Li⁺ with Cl^- is relatively low and the disassociation of the Li⁺ cluster is easier so that Li⁺ can more easily pass through the nanopores. These results gain insight into the effect of ion association inspired by electric field and nanoconfinement of graphene nanopore on Mg²⁺ /Li⁺ separation, and provide helpful information for the application of nanoporous materials in extraction of Li⁺ ion from salt-lake brine.

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.

Benchmark assessment of molecular geometries and energies from small molecule force fields

Background: Force fields are used in a wide variety of contexts for classical molecular simulation, including studies on protein-ligand binding, membrane permeation, and thermophysical property prediction. The quality of these studies relies on the quality of the force fields used to represent the systems. Methods: Focusing on small molecules of fewer than 50 heavy atoms, our aim in this work is to compare nine force fields: GAFF, GAFF2, MMFF94, MMFF94S, OPLS3e, SMIRNOFF99Frosst, and the Open Force Field Parsley, versions 1.0, 1.1, and 1.2. On a dataset comprising 22,675 molecular structures of 3,271 molecules, we analyzed force field-optimized geometries and conformer energies compared to reference quantum mechanical (QM) data. Results: We show that while OPLS3e performs best, the latest Open Force Field Parsley release is approaching a comparable level of accuracy in reproducing QM geometries and energetics for this set of molecules. Meanwhile, the performance of established force fields such as MMFF94S and GAFF2 is generally somewhat worse. We also find that the series of recent Open Force Field versions provide significant increases in accuracy. Conclusions: This study provides an extensive test of the performance of different molecular mechanics force fields on a diverse molecule set, and highlights two (OPLS3e and OpenFF 1.2) that perform better than the others tested on the present comparison. Our molecule set and results are available for other researchers to use in testing.

Optimized utilization of COMB3 reactive potentials in LAMMPS

An investigation to optimize the application of the third-generation charge optimized many-body (COMB3) interatomic potential and associated input parameters was carried out through the study of solid-liquid interactions in classical molecular dynamics simulations. The rates of these molecular interactions are understood through the wetting rates of water nano-droplets on a bare copper (111) surface. Implementing the Langevin thermostat, the influence of simulation time step, the number of atoms in the system, the frequency at which charge equilibration is performed, and the temperature relaxation rate are all examined. The results indicate that time steps of 0.4 fs are possible when using longer relaxation times for the system temperature, which is almost double the typical time step used for reactive potentials. The use of the charge equilibration allows for a fewer atomic layers to be used in the Cu slab. In addition, charge equilibrium schemes do not need to be performed every time step to ensure accurate charge transfer. Interestingly, the rate of wetting for the nanodroplets is dominantly dependent on the temperature relaxation time, which is predicted to significantly change the viscosity of the water droplets. This work provides a pathway for optimizing simulations using the COMB3 reactive interatomic potential.

Surface Energies and Structure of Salt-Brine Interfaces

Permeability of salt formations is controlled by the equilibrium between the salt-brine and salt-salt interfaces described by the dihedral angle, which can change with the composition of the intergranular brine. Here, classical molecular dynamics (MD) simulations were used to investigate the structure and properties of the salt-brine interface to provide insight into the stability of salt systems. Mixed NaCl-KCl brines were investigated to explore differences in ion size on the surface energy and interface structure. Nonlinearity was noted in the salt-brine surface energy with increasing KCl concentration, and the addition of 10% KCl increased surface energies by 2-3 times (5.0 M systems). Size differences in Na⁺ and K⁺ ions altered the packing of dissolved ions and water molecules at the interface, impacting the surface energy. Additionally, ions at the interface had lower numbers of coordinating water molecules than those in the bulk and increased hydration for ions in systems with 100% NaCl or 100% KCl brines. Ultimately, small changes in brine composition away from pure NaCl altered the structure of the salt-brine interface, impacting the dihedral angle and the predicted equilibrium permeability of salt formations.

Molecular Dynamics and Neutron Scattering Studies of Potassium Chloride in Aqueous Solution

Neutron diffraction with isotopic substitution (NDIS) experiments were done on both natural abundance potassium and isotopically labeled ⁴¹KCl heavy water solutions to characterize the solvent structuring around the potassium ion in water. Preliminary measurements suggested that the literature value for the coherent neutron scattering length (2.69 fm) for ⁴¹K was significantly in error. This value was remeasured using a neutron powder diffractometer and found to be 2.40 fm. This revision increases significantly the contrast between the natural abundance K and ⁴¹K by about 30% (from 1.0 to 1.3 fm). The experimentally determined structure factor of the potassium ion was then compared to that calculated from molecular dynamics (MD) simulations. Previous neutron scattering measurements of potassium gave a solvation number of 5.5 (see below). In this study, the NDIS and MD results are in good agreement and allowed us to derive a coordination number of 6.1 for water molecules and 0.8 for chloride ions around each K⁺ ion in 4 molal aqueous KCl solution.

Glass polymorphism and liquid-liquid phase transition in aqueous solutions: experiments and computer simulations

One of the most intriguing anomalies of water is its ability to exist as distinct amorphous ice forms (glass polymorphism or polyamorphism). This resonates well with the possible first-order liquid-liquid phase transition (LLPT) in the supercooled state, where ice is the stable phase. In this Perspective, we review experiments and computer simulations that search for LLPT and polyamorphism in aqueous solutions containing salts and alcohols. Most studies on ionic solutes are devoted to NaCl and LiCl; studies on alcohols have mainly focused on glycerol. Less attention has been paid to protein solutions and hydrophobic solutes, even though they reveal promising avenues. While all solutions show polyamorphism and an LLPT only in dilute, sub-eutectic mixtures, there are differences regarding the nature of the transition. Isocompositional transitions for varying mole fractions are observed in alcohol but not in ionic solutions. This is because water can surround alcohol molecules either in a low- or high-density configuration whereas for ionic solutes, the water ion hydration shell is forced into high-density structures. Consequently, the polyamorphic transition and the LLPT are prevented near the ions, but take place in patches of water within the solutions. We highlight discrepancies and different interpretations within the experimental community as well as the key challenges that need consideration when comparing experiments and simulations. We point out where reinterpretation of past studies helps to draw a unified, consistent picture. In addition to the literature review, we provide original experimental results. A list of eleven open questions that need further consideration is identified.

Mineralization of phosphorylated cellulose: crucial role of surface structure and monovalent ions for optimizing calcium content

Cellulose can be phosphorylated to produce organic matrices with highly adsorptive properties for, e.g., biocompatible materials for biomedical applications. We focus on the effects of phosphorylation of surfaces of crystalline nanocellulose and, in particular, on the competitive adsorption of mono- and divalent cations (Na+ and Ca2+) typically contained in mineralizing salt mixtures using all-atom molecular dynamics (MD) simulations. Phosphorylation was applied at 12% and 25% both in water and CaCl2 solutions. Our main result shows that Na+ and Ca2+ cations are concentrated in different interfacial layers with Na+ ions penetrating much closer to the surface. This behavior cannot be described by the Poisson-Boltzmann theory or implicit solvent simulations. Our analysis shows that the physical origin of this observation is due to a balance between the electrostatic interactions and hydration free energy associated with the ions. Adsorption levels of the different ions also respond differently to changes in the degree of phosphorylation. We show that the number of adsorbed Na+ ions per phosphate group increases whereas the number of adsorbed Ca2+ ions decreases with an increasing degree of phosphorylation (or when the number of binding sites increases). The decrease in the number of adsorbed Ca2+ ions can be explained by an increasing "charge-charge" repulsion between the Ca2+ ions attracted by the charged surface. Importantly, our results demonstrate the existence of an optimum degree of phosphorylation in terms of adsorbed Ca2+ ions and can be used as a guideline in materials design, for example, when choosing the cellulose matrix or with other similarly structured biomolecular and polymer surfaces.

Charged nanoporous graphene membranes for water desalination

Water desalination using positively and negatively charged single-layer nanoporous graphene membranes.

Single-Ion Thermodynamics from First Principles: Calculation of the Absolute Hydration Free Energy and Single-Electrode Potential of Aqueous Li+ Using ab Initio Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations

A recently proposed thermodynamic integration (TI) approach formulated in the framework of quantum mechanical/molecular mechanical molecular dynamics (QM/MM MD) simulations is applied to study the structure, dynamics, and absolute intrinsic hydration free energy Δ_s G_M⁺,wat^◦ of the Li⁺ ion at a correlated ab initio level of theory. Based on the results, standard values (298.15 K, ideal gas at 1 bar, ideal solute at 1 molal) for the absolute intrinsic hydration free energy [Formula: see text] of the proton, the surface electric potential jump χ_wat^◦ upon entering bulk water, and the absolute single-electrode potential [Formula: see text] of the reference hydrogen electrode are calculated to be -1099.9 ± 4.2 kJ·mol^-1, 0.13 ± 0.08 V, and 4.28 ± 0.04 V, respectively, in excellent agreement with the standard values recommended by Hünenberger and Reif on the basis of an extensive evaluation of the available experimental data (-1100 ± 5 kJ·mol^-1, 0.13 ± 0.10 V, and 4.28 ± 0.13 V). The simulation results for Li⁺ are also compared to those for Na⁺ and K⁺ from a previous study in terms of relative hydration free energies ΔΔ_s G_M⁺,wat^◦ and relative electrode potentials [Formula: see text]. The calculated values are found to agree extremely well with the experimental differences in standard conventional hydration free energies ΔΔ_s G_M⁺,wat^• and redox potentials [Formula: see text]. The level of agreement between simulation and experiment, which is quantitative within error bars, underlines the substantial accuracy improvement achieved by applying a highly demanding QM/MM approach at the resolution-of-identity second-order Møller-Plesset perturbation (RIMP2) level over calculations relying on purely molecular mechanical or density functional theory (DFT) descriptions. A detailed analysis of the structural and dynamical properties of the Li⁺ hydrate indicates that a correct description of the solvation structure and dynamics is achieved as well at this level of theory. Consideration of the QM/MM potential-energy components also shows that the partitioning into QM and MM zones does not induce any significant energetic artifact for the system considered.

Structure, interactions and action of Mycobacterium tuberculosis 3-hydroxyisobutyric acid dehydrogenase

Biochemical and crystallographic studies on Mycobacterium tuberculosis 3-hydroxyisobutyric acid dehydrogenase (MtHIBADH), a member of the 3-hydroxyacid dehydrogenase superfamily, have been carried out. Gel filtration and blue native PAGE of MtHIBADH show that the enzyme is a dimer. The enzyme preferentially uses NAD+ as the cofactor and is specific to S-hydroxyisobutyric acid (HIBA). It can also use R-HIBA, l-serine and 3-hydroxypropanoic acid (3-HP) as substrates, but with much less efficiency. The pH optimum for activity is ∼11. Structures of the native enzyme, the holoenzyme, binary complexes with NAD+, S-HIBA, R-HIBA, l-serine and 3-HP and ternary complexes involving the substrates and NAD+ have been determined. None of the already known structures of HIBADH contain a substrate molecule at the binding site. The structures reported here provide for the first time, among other things, a clear indication of the location and interactions of the substrates at the active site. They also define the entrance of the substrates to the active site region. The structures provide information on the role of specific residues at the active site and the entrance. The results obtained from crystal structures are consistent with solution studies including mutational analysis. They lead to the proposal of a plausible mechanism of the action of the enzyme.

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.

Quantum Chemical Modeling of Hydrogen Bonding in Ionic Liquids

Hydrogen bonding (H-bonding) is an important and very general phenomenon. H-bonding is part of the basis of life in DNA, key in controlling the properties of water and ice, and critical to modern applications such as crystal engineering, catalysis applications, pharmaceutical and agrochemical development. H-bonding also plays a significant role for many ionic liquids (IL), determining the secondary structuring and affecting key physical parameters. ILs exhibit a particularly diverse and wide range of traditional as well as non-standard forms of H-bonding, in particular the doubly ionic H-bond is important. Understanding the fundamental nature of the H-bonds that form within ILs is critical, and one way of accessing this information, that cannot be recovered by any other computational method, is through quantum chemical electronic structure calculations. However, an appropriate method and basis set must be employed, and a robust procedure for determining key structures is essential. Modern generalised solvation models have recently been extended to ILs, bringing both advantages and disadvantages. QC can provide a range of information on geometry, IR and Raman spectra, NMR spectra and at a more fundamental level through analysis of the electronic structure.

Correlation of natural muscovite exfoliation with interlayer and solvation forces

Natural muscovite exfoliation was correlated with interlayer and solvation forces respectively in this work.

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.

Molecular dynamics simulation of water permeation through the alpha-hemolysin channel

The alpha-hemolysin (AHL) nanochannel is a non-selective channel that allows for uncontrolled transport of small molecules across membranes leading to cell death. Although it is a bacterial toxin, it has promising applications, ranging from drug delivery systems to nano-sensing devices. This study focuses on the transport of water molecules through an AHL nanochannel using molecular dynamics (MD) simulations. Our results show that AHL can quickly transport water across membranes. The first-passage time approach was used to estimate the diffusion coefficient and the mean exit time. To study the energetics of transport, the potential of mean force (PMF) of a water molecule along the AHL nanochannel was calculated. The results show that the energy barriers of water permeation across a nanopore are always positive along the channel and the values are close to thermal energy (kBT). These findings suggest that the observed quick permeation of water is due to small energy barriers and a hydrophobic inner channel surface resulting in smaller friction. We speculate that these physical mechanisms are important in how AHL causes cell death.

Binding of Disordered Peptides to Kelch: Insights from Enhanced Sampling Simulations

Keap1 protein plays an essential role in regulating cellular oxidative stress response and is a crucial binding hub for multiple proteins, several of which are intrinsically disordered proteins (IDP). Among Kelch's IDP binding partners, NRF2 and PTMA are the two most interesting cases. They share a highly similar binding motif; however, NRF2 binds to Kelch with a binding affinity of approximately 100-fold higher than that of PTMA. In this study, we perform an exhaustive sampling composed of 6 μs well-tempered metadynamics and 2 μs unbiased molecular dynamics (MD) simulations aiming at characterizing the binding mechanisms and structural properties of these two peptides. Our results agree with previous experimental observations that PTMA is remarkably more disordered than NRF2 in both the free and bound states. This explains PTMA's lower binding affinity. Our extensive sampling also provides valuable insights into the vast conformational ensembles of both NRF2 and PTMA, supports the hypothesis of coupled folding-binding, and confirms the essential role of linear motifs in IDP binding.

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.

Halide, Ammonium, and Alkali Metal Ion Parameters for Modeling Aqueous Solutions

A complete set of Lennard Jones parameters for the halide ions, F(-), Cl(-), Br(-), and I(-), ammonium ion, and the alkali metal ions is reported. The parameters have been optimized using Monte Carlo simulations and free energy perturbation theory with the TIP4P water model to reproduce experimental free energies of hydration and locations of the first maxima of the ion-oxygen radial distribution functions, to provide water coordination numbers consistent with experimental ranges, and to exhibit gas-phase monohydrate energies in reasonable agreement with ab initio values. Average errors for absolute and relative free energies of hydration for the ions are ca. 1 kcal/mol. For the halides, this is the first self-consistent set of parameters that has been optimized for aqueous-phase performance. The good results for relative free energies of hydration are particularly auspicious for use of the new parameters in a wide variety of liquid-phase simulations where halide and alkali cations are systematically varied.

Effect of Force Field Parameters on Sodium and Potassium Ion Binding to Dipalmitoyl Phosphatidylcholine Bilayers

The behavior of electrolytes in molecular dynamics simulations of zwitterionic phospholipid bilayers is very sensitive to the force field parameters used. Here, several 200 ns molecular dynamics of simulations of dipalmitoyl phosphotidylcholine (PC) bilayers in 0.2 M sodium or potassium chloride using various common force field parameters for the cations are presented. All employed parameter sets give a larger number of Na(+) ions than K(+) ions that bind to the lipid heads, but depending on the parameter choice quite different results are seen. A wide range of coordination numbers for the Na(+) and K(+) ions is also observed. These findings have been analyzed and compared to published experimental data. Some simulations produce aggregates of potassium chloride, indicating (in accordance with published simulations) that these force fields do not reproduce the delicate balance between salt and solvated ions. The differences between the force fields can be characterized by one single parameter, the electrostatic radius of the ion, which is correlated to σMO (M represents Na(+)/K(+)), the Lennard-Jones radius. When this parameter exceeds a certain threshold, binding to the lipid heads is no longer observed. One would, however, need more accurate experimental data to judge or rank the different force fields precisely. Still, reasons for the poor performance of some of the parameter sets are clearly demonstrated, and a quality control procedure is provided.

Accounting for Electronic Polarization Effects in Aqueous Sodium Chloride via Molecular Dynamics Aided by Neutron Scattering

Modeled ions, described by nonpolarizable force fields, can suffer from unphysical ion pairing and clustering in aqueous solutions well below their solubility limit. The electronic continuum correction takes electronic polarization effects of the solvent into account in an effective way by scaling the charges on the ions, resulting in a much better description of the ionic behavior. Here, we present parameters for the sodium ion consistent with this effective polarizability approach and in agreement with experimental data from neutron scattering, which could be used for simulations of complex aqueous systems where polarization effects are important.