A transferable ab initio based force field for aqueous ions

Electrochemical Reduction of N2 to Ammonia Promoted by Hydrated Cation Ions: Mechanistic Insights from a Combined Computational and Experimental Study

Alkali ions, major components at the electrode-electrolyte interface, are crucial to modulating reaction activity and selectivity of catalyst materials. However, the underlying mechanism of how the alkali ions catalyze the N2 reduction reaction (NRR) into ammonia remains elusive, posing challenges for experimentalists to select appropriate electrolyte solutions. In this work, by employing a combined experimental and computational approach, we proposed four essential roles of cation ions at Fe electrodes for N2 fixation: (i) promoting NN bond cleavage; (ii) stabilizing NRR intermediates; (iii) suppressing the competing hydrogen evolution reaction (HER); and (iv) modulating the interfacial charge distribution at the electrode-electrolyte interface. For N2 adsorption on an Fe electrode with cation ions, our constrained ab initio molecular dynamic (c-AIMD) results demonstrate a barrierless process, while an extra 0.52 eV barrier requires to be overcome to adsorb N2 for the pure Fe-water interface. For the formation of *NNH species within the N2 reduction process, the calculated free energy barrier is 0.50 eV at the Li+-Fe-water interface. However, the calculated barrier reaches 0.81 eV in pure Fe-water interface. Furthermore, experiments demonstrate a high Faradaic efficiency for ammonia synthesis on a Li+-Fe-water interface, reaching 27.93% at a working potential of -0.3 V vs RHE and pH = 6.8. These results emphasize how alkali metal cations and local reaction environments on the electrode surface play crucial roles in influencing the kinetics of interfacial reactions.

Impact of hierarchical water dipole orderings on the dynamics of aqueous salt solutions

Ions exhibit highly ion-specific complex behaviours when solvated in water, which remains a mystery despite the fundamental importance of ion solvation in nature, science, and technology. Here we explain these ion-specific properties by the ion-induced hierarchical dipolar, translational, and bond-orientational orderings of ion hydration shell under the competition between ion-water electrostatic interactions and inter-water hydrogen bonding. We first characterise this competition by a new length λHB(q), explaining the ion-specific effects on solution dynamics. Then, by continuously tuning ion size and charge, we find that the bond-orientational order of the ion hydration shell highly develops for specific ion size and charge combinations. This ordering drastically stabilises the hydration shell; its degree changes the water residence time around ions by 11 orders of magnitude for main-group ions. These findings are fundamental to ionic processes in aqueous solutions, providing a physical principle for electrolyte design and application.

Nonpolarizable Force Fields through the Self-Consistent Modeling Scheme with MD and DFT Methods: From Ionic Liquids to Self-Assembled Ionic Liquid Crystals

A key to achieve the accuracy of molecular dynamics (MD) simulation is the set of force fields used to express the atomistic interactions. In particular, the electrostatic interaction remains the main issue for the precise simulation of various ionic soft materials from ionic liquids to their supramolecular compounds. In this study, we test the nonpolarizable force fields of ionic liquids (ILs) and self-assembled ionic liquid crystals (ILCs) for which the intermolecular charge transfer and intramolecular polarization are significant. The self-consistent modeling scheme is adopted to refine the atomic charges of ionic species in a condensed state through the use of density functional theory (DFT) under the periodic boundary condition. The atomic charges of the generalized amber force field (GAFF) are effectively updated to express the electrostatic properties of ionic molecules obtained by the DFT calculation in condensed phase, which improves the prediction accuracy of ionic conductivity with the obtained force field (GAFF-DFT). The derived DFT charges then suggest that the substitution of a hydrophobic liquid-crystalline moiety into IL-based cations enhances the charge localization of ionic groups in the amphiphilic molecules, leading to the amplification of the electrostatic interactions among the hydrophilic/ionic groups in the presence of hydrophobic moieties. In addition, we focus on an ion-conductive pathway hidden in the self-assembled nanostructure. The MD results indicate that the ionic groups of cation and anion interact strongly for keeping the bicontinuous nanosegregation of ionic nanochannel. The partial fractions of hydrophilic/ionic and hydrophobic nanodomains are then quantified with the volume difference from referenced IL systems, while the calculated ionic conductivity decreases in the self-assembled ILCs more than the occupied volume of ionic nanodomains. These analyses suggest that the mobility of ions in the self-assembled ILCs remains quite restricted even with small tetrafluoroborate anions because of strong attractive interaction among ionic moieties.

A systematic and critical review of application of molecular dynamics simulation in low salinity water injection

Low Salinity Water Injection (LSWI) has been a well-researched EOR method, with several experimental and theoretical scientific papers reported in the literature over the past few decades. Despite this, there is still an ongoing debate on dominant mechanisms behind this complex EOR process, and some issues remain elusive. Part of the complexity arises from the scale of investigation, which spans from sub-pore scale (atomic and electronic scale) to pore scale, core scale, and reservoir scale. Molecular Dynamics (MD) simulation has been used as a research tool in the past decade to investigate the nano-scale interactions among reservoir rock (e.g., calcite, silica), crude oil, and brine systems in presence of some impurities (e.g., clay minerals) and additives (e.g., nanoparticles). In this paper, fundamental concepts of MD simulation and common analyses driven by MD are briefly reviewed. Then, an overview of molecular models of the most common minerals encountered in petroleum reservoirs: quartz, calcite, and clay, with their most common types of potential function, is provided. Next, a critical review and in depth analysis of application of MD simulations in LSWI process in both sandstone and carbonate reservoirs in terms of sub-pore scale mechanisms, namely electrical double layer (EDL) expansion, multi-ion exchange (MIE), and cation hydration, is presented to scrutinize role of salinity, ionic composition, and rock surface chemistry from an atomic level. Some inconsistencies observed in the literature are also highlighted and the reasons behind them are explained. Finally, a future research guide is provided after critically discussing the challenges and potential of the MD in LSWI to shed more light on governing mechanisms behind LSWI by enhancing the reliability of MD outcomes in future researches. Such insights can be used for design of new MD researches with complementary experimental studies at core scale to capture the main mechanisms behind LSWI.

Rigid Cations Induce Enhancement of Microheterogeneity and Exhibit Anomalous Ion Diffusion in Water-Ethanol Mixtures

Because of the amphiphilic nature of ethanol in the aqueous solution, ions cause an interesting microheterogeneity where the water molecules and the hydroxy groups of ethanol preferentially solvate the ions, while the ethyl groups tend to occupy the intervening space. Using computer simulations, we study the dynamics of rigid monovalent cations (Li⁺, Na⁺, K⁺, and Cs⁺) in aqueous ethanol solutions with chloride as the counterion. We vary both the size of the ions and the composition of the mixture to explore size- and composition-dependent ion diffusion. The relative stability of enhanced microheterogeneous configurations makes ion diffusion slower than what would be surmised by using the bulk properties of the mixture, using the Stokes-Einstein relation. We study the structure through partial radial distribution functions and the stability through coordination number fluctuations. The ion diffusion coefficient exhibits sharp re-entrant behavior when plotted against viscosity varied by composition. Our studies reveal multiple anomalous features of ion motion in this mixture. We formulate a mode-coupling theory (MCT) that takes into account the interaction between different dynamical components; MCT can incorporate the effects of heterogeneous dynamics and nonlinearity in composition dependence that arise from the feedback between mutually dependent ion-solvent dynamics.

NMR Relaxation Rates of Quadrupolar Aqueous Ions from Classical Molecular Dynamics Using Force-Field Specific Sternheimer Factors

The nuclear magnetic resonance (NMR) relaxation of quadrupolar nuclei is governed by the electric field gradient (EFG) fluctuations at their position. In classical molecular dynamics (MD), the electron cloud contribution to the EFG can be included via the Sternheimer approximation, in which the full EFG at the nucleus that can be computed using quantum density functional theory (DFT) is considered to be proportional to that arising from the external, classical charge distribution. In this work, we systematically assess the quality of the Sternheimer approximation as well as the impact of the classical force field (FF) on the NMR relaxation rates of aqueous quadrupolar ions at infinite dilution. In particular, we compare the rates obtained using an ab initio parametrized polarizable FF, a recently developed empirical FF with scaled ionic charges and a simple empirical nonpolarizable FF with formal ionic charges. Surprisingly, all three FFs considered yield good values for the rates of smaller and less polarizable solutes (Li⁺, Na⁺, K⁺, Cl^-), provided that a model-specific Sternheimer parametrization is employed. Yet, the polarizable and scaled charge FFs yield better estimates for divalent and more polarizable species (Mg²⁺, Ca²⁺, Cs⁺). We find that a linear relationship between the quantum and classical EFGs holds well in all of the cases considered; however, such an approximation often leads to quite large errors in the resulting EFG variance, which is directly proportional to the computed rate. We attempted to reduce the errors by including first order nonlinear corrections to the EFG, yet no clear improvement for the resulting variance has been found. The latter result indicates that more refined methods for determining the EFG at the ion position, in particular those that take into account the instantaneous atomic environment around an ion, might be necessary to systematically improve the NMR relaxation rate estimates in classical MD.

Properties of the tetravalent actinide series in aqueous phase from a microscopic simulation self-consistent engine

In the context of nuclear fuel recycling and environmental issues, the understanding of the properties of radio-elements with various approaches remains a challenge regarding their dangerousness. Moreover, experimentally, some issues are also of importance; first, it is imperative to work at sufficiently high concentrations to reach the sensitivities of the analytical tools, however this condition often leads to precipitation for some of them; second, stabilizing specific oxidation states of some actinides remains a challenge, thus making it difficult to extract general trends across the actinide series. Complementary to experiments, modeling can be used to unbiasedly probe the actinide's properties in an aquatic environment and offers a predictive tool. We report the first molecular dynamics simulations based on homogeneously built force fields for the whole series of the tetravalent actinides in aqueous phase from ThIV to BkIV and including PuIV. The force fields used to model the interactions among the constituents include polarization and charge donation microscopic effects. They are built from a self-consistent iterative ab initio based engine that can be included in future developments as an element of a potential machine learning procedure devoted to generating accurate force fields. The comparison of our simulated hydrated actinide properties to available experimental data shows the model robustness and the relevance of our parameter assignment engine. Moreover, our simulated structural, dynamical and evolution of the hydration free energy data show that, apart from AmIV and CmIV, the actinide properties change progressively along the series.

Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields

Many applications in chemistry, biology, and energy storage/conversion research rely on molecular simulations to provide fundamental insight into structural and transport properties of materials with high ionic concentrations. Whether the system is comprised entirely of ions, like ionic liquids, or is a mixture of a polar solvent with a salt, e.g., liquid electrolytes for battery applications, the presence of ions in these materials results in strong local electric fields polarizing solvent molecules and large ions. To predict properties of such systems from molecular simulations often requires either explicit or mean-field inclusion of the influence of polarization on electrostatic interactions. In this manuscript, we review the pros and cons of different treatments of polarization ranging from the mean-field approaches to the most popular explicit polarization models in molecular dynamics simulations of ionic materials. For each method, we discuss their advantages and disadvantages and emphasize key assumptions as well as their adjustable parameters. Strategies for the development of polarizable models are presented with a specific focus on extracting atomic polarizabilities. Finally, we compare simulations using polarizable and nonpolarizable models for several classes of ionic systems, discussing the underlying physics that each approach includes or ignores, implications for implementation and computational efficiency, and the accuracy of properties predicted by these methods compared to experiments.

Polarizable Force Fields for Biomolecular Simulations: Recent Advances and Applications

Realistic modeling of biomolecular systems requires an accurate treatment of electrostatics, including electronic polarization. Due to recent advances in physical models, simulation algorithms, and computing hardware, biomolecular simulations with advanced force fields at biologically relevant timescales are becoming increasingly promising. These advancements have not only led to new biophysical insights but also afforded opportunities to advance our understanding of fundamental intermolecular forces. This article describes the recent advances and applications, as well as future directions, of polarizable force fields in biomolecular simulations.

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.

Classical Polarizable Force Field to Study Hydrated Hectorite: Optimization on DFT Calculations and Validation against XRD Data

Following our previous works on dioctahedral clays, we extend the classical Polarizable Ion Model (PIM) to trioctahedral clays, by considering dry Na-, Cs-, Ca- and Sr-hectorites as well as hydrated Na-hectorite. The parameters of the force field are determined by optimizing the atomic forces and dipoles on density functional theory calculations. The simulation results are validated by comparison with experimental X-ray diffraction (XRD) data. The XRD patterns calculated from classical molecular dynamics simulations performed with the PIM force field are in very good agreement with experimental results. In the bihydrated state, the less structured electronic density profile obtained with PIM compared to the one from the state-of-the-art non-polarizable force field clayFF explains the slightly better agreement between the PIM results and experiments.

Underscreening in ionic liquids: a first principles analysis

An attempt is made to understand the underscreening effect, observed in concentrated electrolyte solutions or melts, on the basis of simple, admittedly crude models involving charged (for the ions) and neutral (for the solvent molecules) hard spheres. The thermodynamic and structural properties of these 'primitive' and 'semi-primitive' models are calculated within mean spherical approximation, which provides the basic input required to determine the partial density response functions. The screening length [Formula: see text], which is unambiguously defined in terms of the wave-number-dependent response functions, exhibits a cross-over from a low density, Debye-like regime, to a regime where [Formula: see text] increases with density beyond a critical density at which the Debye length [Formula: see text] becomes comparable to the ion diameter. In this high density regime the ratio [Formula: see text] increases according to a power law, in qualitative agreement with experimental measurements, albeit at a much slower rate.

Ca2+ -Cl- Association in Water Revisited: the Role of Cation Hydration

We investigate the dissociation of a Ca²⁺ -Cl^- pair in water using classical molecular dynamics simulations with a polarizable interaction potential, parameterized from ab initio calculations. By computing the potential of mean force as a function not only of the interionic distance but also of the coordination numbers by water molecules, we show that it is necessary to use a collective variable describing the cation hydration in order to capture the dissociation mechanism. In the contact ion pair, the Ca²⁺ cation has a first coordination sphere containing 5 or 6 water molecules. The minimum free-energy path for dissociation involves a two-step process: First one or two additional water molecules enter the cation coordination shell, increasing the coordination number up to 7 with an almost fixed interionic distance. Then the dissociation of the ionic pair occurs at this fixed coordination number.

Collective water dynamics in the first solvation shell drive the NMR relaxation of aqueous quadrupolar cations

Using molecular simulations, we analyze the microscopic processes driving the Nuclear Magnetic Resonance (NMR) relaxation of quadrupolar cations in water. The fluctuations of the Electric Field Gradient (EFG) experienced by alkaline and magnesium cations, which determine the NMR relaxation time, are mainly due to the dynamics of water molecules in their solvation shell. The dynamics of the ion plays a less important role, with the exception of the short-time dynamics in the lighter Li⁺ case, for which rattling in the solvent cage results in oscillations of the EFG autocorrelation function (ACF). Several microscopic mechanisms that may a priori contribute to the decay of the EFG-ACF occur in fact over too long time scales: entrance/exit of individual water molecules into/from the solvation shell, rotation of a molecule around the ion, or reorientation of the molecule. In contrast, the fluctuations of the ion-water distance are clearly correlated to that of the EFG. Nevertheless, it is not sufficient to consider a single molecule due to the cancellations arising from the symmetry of the solvation shell. The decay of the EFG-ACF, hence NMR relaxation, is in fact governed by the collective symmetry-breaking fluctuations of water in the first solvation shell.

Probing ice VII crystallization from amorphous NaCl–D2O solutions at gigapascal pressures

The high density amorphous solution NaCl·10.2D₂O crystallises at 260 K as almost pure ice VII during annealing at gigapascal pressures.

Multidimensional materials and device architectures for future hybrid energy storage

Electrical energy storage plays a vital role in daily life due to our dependence on numerous portable electronic devices. Moreover, with the continued miniaturization of electronics, integration of wireless devices into our homes and clothes and the widely anticipated ‘Internet of Things', there are intensive efforts to develop miniature yet powerful electrical energy storage devices. This review addresses the cutting edge of electrical energy storage technology, outlining approaches to overcome current limitations and providing future research directions towards the next generation of electrical energy storage devices whose characteristics represent a true hybridization of batteries and electrochemical capacitors.

With the continued miniaturization of electronics, there are increasing efforts to engineer small, powerful energy storage devices. Here the authors review the cutting edge of this rapidly developing field, highlighting the most promising materials and architectures for our future energy storage requirements.

On the microscopic fluctuations driving the NMR relaxation of quadrupolar ions in water

Nuclear Magnetic Resonance (NMR) relaxation is sensitive to the local structure and dynamics around the probed nuclei. The Electric Field Gradient (EFG) is the key microscopic quantity to understand the NMR relaxation of quadrupolar ions, such as (7)Li(+), (23)Na(+), (25)Mg(2+), (35)Cl(-), (39)K(+), or (133)Cs(+). Using molecular dynamics simulations, we investigate the statistical and dynamical properties of the EFG experienced by alkaline, alkaline Earth, and chloride ions at infinite dilution in water. Specifically, we analyze the effect of the ionic charge and size on the distribution of the EFG tensor and on the multi-step decay of its auto-correlation function. The main contribution to the NMR relaxation time arises from the slowest mode, with a characteristic time on the picosecond time scale. The first solvation shell of the ion plays a dominant role in the fluctuations of the EFG, all the more that the ion radius is small and its charge is large. We propose an analysis based on a simplified charge distribution around the ion, which demonstrates that the auto-correlation of the EFG, hence the NMR relaxation time, reflects primarily the collective translational motion of water molecules in the first solvation shell of the cations. Our findings provide a microscopic route to the quantitative interpretation of NMR relaxation measurements and open the way to the design of improved analytical theories for NMR relaxation for small ionic solutes, which should focus on water density fluctuations around the ion.

Force Field for Mg(2+), Mn(2+), Zn(2+), and Cd(2+) Ions That Have Balanced Interactions with Nucleic Acids

Divalent metal ions are of fundamental importance to the function and folding of nucleic acids. Divalent metal ion-nucleic acid interactions are complex in nature and include both territorial and site specific binding. Commonly employed nonbonded divalent ion models, however, are often parametrized against bulk ion properties and are subsequently utilized in biomolecular simulations without considering any data related to interactions at specific nucleic acid sites. Previously, we assessed the ability of 17 different nonbonded Mg(2+) ion models to reproduce different properties of Mg(2+) in aqueous solution including radial distribution functions, solvation free energies, water exchange rates, and translational diffusion coefficients. In the present work, we depart from the recently developed 12-6-4 potential models for divalent metal ions developed by Li and Merz and tune the pairwise parameters for Mg(2+), Mn(2+), Zn(2+), and Cd(2+) binding dimethyl phosphate, adenosine, and guanosine in order to reproduce experimental site specific binding free energies derived from potentiometric pH titration data. We further apply these parameters to investigate a metal ion migration previously proposed to occur during the catalytic reaction of the hammerhead ribozyme. The new parameters are shown to be accurate and balanced for nucleic acid binding in comparison with available experimental data and provide an important tool for molecular dynamics and free energy simulations of nucleic acids where these ions may exhibit different binding modes.

Capacitive Energy Storage: Current and Future Challenges

Capacitive energy storage devices are receiving increasing experimental and theoretical attention due to their enormous potential for energy applications. Current research in this field is focused on the improvement of both the energy and the power density of supercapacitors by optimizing the nanostructure of porous electrodes and the chemical structure/composition of the electrolytes. However, the understanding of the underlying correlations and the mechanisms of electric double layer formation near charged surfaces and inside nanoporous electrodes is complicated by the complex interplay of several molecular scale phenomena. This Perspective presents several aspects regarding the experimental and theoretical research in the field, discusses the current atomistic and molecular scale understanding of the mechanisms of energy and charge storage, and provides a brief outlook to the future developments and applications of these devices.

Tuning the Electronic Structure of Anatase Through Fluorination

A highly fluorinated anatase lattice has been recently reported, providing a new class of materials whose general chemical formula is Ti_1−x□_xX_4xO_2−4x (X⁻ = F⁻ or OH⁻). To characterise the complex structural features of the material and the different F environments, we here apply a computational screening procedure. After deriving a polarisable force—field from DFT simulations, we screen in a step-wise fashion a large number of possible configurations differing in the positioning of the titanium vacancies (□) and of the fluorine atoms. At each step only 10% of the configurations are retained. At the end of the screening procedure, a configuration is selected and simulated using DFT-based molecular dynamics. This allows us to analyse the atomic structure of the material, which is strongly disordered, leading to a strong decrease (by 0.8 eV) of the band gap compared to conventional anatase.

Direct computation of parameters for accurate polarizable force fields

We present an improved electronic linear response model to incorporate polarization and charge-transfer effects in polarizable force fields. This model is a generalization of the Atom-Condensed Kohn-Sham Density Functional Theory (DFT), approximated to second order (ACKS2): it can now be defined with any underlying variational theory (next to KS-DFT) and it can include atomic multipoles and off-center basis functions. Parameters in this model are computed efficiently as expectation values of an electronic wavefunction, obviating the need for their calibration, regularization, and manual tuning. In the limit of a complete density and potential basis set in the ACKS2 model, the linear response properties of the underlying theory for a given molecular geometry are reproduced exactly. A numerical validation with a test set of 110 molecules shows that very accurate models can already be obtained with fluctuating charges and dipoles. These features greatly facilitate the development of polarizable force fields.

Structural characterization of eutectic aqueous NaCl solutions under variable temperature and pressure conditions

The structure of amorphous NaCl solutions produced by fast quenching is studied as a function of pressure, up to 4 GPa, by combined neutron diffraction experiments and classical molecular dynamics simulations. Similarly to LiCl solutions the system amorphizes at ambient pressure in a dense phase structurally similar to the e-HDA phase in pure water. The measurement of the static structure factor as a function of pressure allowed us to validate a new polarizable force field developed by Tazi et al., 2012, never tested under non-ambient conditions. We infer from simulations that the hydration shells of Na(+) cations form well defined octahedra composed of both H2O molecules and Cl(-) anions at low pressure. These octahedra are gradually broken by the seventh neighbour moving into the shell of first neighbours yielding an irregular geometry. In contrast to LiCl solutions and pure water, the system does not show a polyamorphic transition under pressure. This confirms that the existence of polyamorphism relies on the tetrahedral structure of water molecules, which is broken here.

Dissolution of NaCl nanocrystals: an ab initio molecular dynamics study

NaCl nanocrystal dissolution was investigated in atomistic detail revealing a difference in the solvation of two different ionic species.

Rational Design of Particle Mesh Ewald Compatible Lennard-Jones Parameters for +2 Metal Cations in Explicit Solvent

Metal ions play significant roles in biological systems. Accurate molecular dynamics (MD) simulations on these systems require a validated set of parameters. Although there are more detailed ways to model metal ions, the nonbonded model, which employs a 12-6 Lennard-Jones (LJ) term plus an electrostatic potential is still widely used in MD simulations today due to its simple form. However, LJ parameters have limited transferability due to different combining rules, various water models and diverse simulation methods. Recently, simulations employing a Particle Mesh Ewald (PME) treatment for long-range electrostatics have become more and more popular owing to their speed and accuracy. In the present work we have systematically designed LJ parameters for 24 +2 metal (M(II)) cations to reproduce different experimental properties appropriate for the Lorentz-Berthelot combining rules and PME simulations. We began by testing the transferability of currently available M(II) ion LJ parameters. The results showed that there are differences between simulations employing Ewald summation with other simulation methods and that it was necessary to design new parameters specific for PME based simulations. Employing the thermodynamic integration (TI) method and performing periodic boundary MD simulations employing PME, allowed for the systematic investigation of the LJ parameter space. Hydration free energies (HFEs), the ion-oxygen distance in the first solvation shell (IOD) and coordination numbers (CNs) were obtained for various combinations of the parameters of the LJ potential for four widely used water models (TIP3P, SPC/E, TIP4P and TIP4PEW). Results showed that the three simulated properties were highly correlated. Meanwhile, M(II) ions with the same parameters in different water models produce remarkably different HFEs but similar structural properties. It is difficult to reproduce various experimental values simultaneously because the nonbonded model underestimates the interaction between the metal ions and water molecules at short range. Moreover, the extent of underestimation increases successively for the TIP3P, SPC/E, TIP4PEW and TIP4P water models. Nonetheless, we fitted a curve to describe the relationship between ε (the well depth) and radius (Rmin/2) from experimental data on noble gases to facilitate the generation of the best possible compromise models. Hence, by targeting different experimental values, we developed three sets of parameters for M(II) cations for three different water models (TIP3P, SPC/E and TIP4PEW). These parameters we feel represent the best possible compromise that can be achieved using the nonbonded model for the ions in combination with simple water models. From a computational uncertainty analysis we estimate that the uncertainty in our computed HFEs is on the order of ±1kcal/mol. Further improvements will require more advanced non-bonded models likely with inclusion of polarization.

Interacting ions in biophysics: real is not ideal

Ions in water are important throughout biology, from molecules to organs. Classically, ions in water were treated as ideal noninteracting particles in a perfect gas. Excess free energy of each ion was zero. Mathematics was not available to deal consistently with flows, or interactions with other ions or boundaries. Nonclassical approaches are needed because ions in biological conditions flow and interact. The concentration gradient of one ion can drive the flow of another, even in a bulk solution. A variational multiscale approach is needed to deal with interactions and flow. The recently developed energetic variational approach to dissipative systems allows mathematically consistent treatment of the bio-ions Na(+), K(+), Ca(2+), and Cl(-) as they interact and flow. Interactions produce large excess free energy that dominate the properties of the high concentration of ions in and near protein active sites, ion channels, and nucleic acids: the number density of ions is often >10 M. Ions in such crowded quarters interact strongly with each other as well as with the surrounding protein. Nonideal behavior found in many experiments has classically been ascribed to allosteric interactions mediated by the protein and its conformation changes. The ion-ion interactions present in crowded solutions-independent of conformation changes of the protein-are likely to change the interpretation of many allosteric phenomena. Computation of all atoms is a popular alternative to the multiscale approach. Such computations involve formidable challenges. Biological systems exist on very different scales from atomic motion. Biological systems exist in ionic mixtures (like extracellular and intracellular solutions), and usually involve flow and trace concentrations of messenger ions (e.g., 10(-7) M Ca(2+)). Energetic variational methods can deal with these characteristic properties of biological systems as we await the maturation and calibration of all-atom simulations of ionic mixtures and divalents.

Driving force for the hydration of the swelling clays: case of montmorillonites saturated with alkaline-earth cations

Important structural modifications occur in swelling clays upon water adsorption. The multi-scale evolution of the swelling clay structure is usually evidenced by various experimental techniques. However, the driving force behind such phenomena is still not thoroughly understood. It appears strongly dependent on the nature of the interlayer cation. In the case of montmorillonites saturated with alkaline cations, it was inferred that the compensating cation or the layer surface could control the hydration process and thus the opening of the interlayer space, depending on the nature of the interlayer cation. In the present study, emphasis is put on the impact of divalent alkaline-earth cations compensating the layer charge in montmorillonites. Since no experimental technique offers the possibility of directly determining the hydration contributions related to interlayer cations and layer surfaces, an approach based on the combination of electrostatic calculations and immersion data is developed here, as already validated in the case of montmorillonites saturated by alkaline cations. This methodology allows to estimate the hydration energy for divalent interlayer cations and therefore to shed a new light on the driving force for hydration process occurring in montmorillonites saturated with alkaline-earth cations. Firstly, the surface energy values obtained from the electrostatic calculations based on the Electronegativity Equalization Method vary from 450 mJ m(-2) for Mg-montmorillonite to 1100 mJ m(-2) for Ba-montmorillonite. Secondly, considering both the hydration energy for cations and layer surfaces, the driving force for the hydration of alkaline-earth saturated montmorillonites can be attributed to the interlayer cation in the case of Mg-, Ca-, Sr-montmorillonites and to the interlayer surface in the case of Ba-montmorillonites. These results explain the differences in behaviour upon water adsorption as a function of the nature of the interlayer cation, thereby allowing the macroscopic swelling trends to be better understood. The knowledge of hydration processes occurring in homoionic montmorillonites saturated with both the alkaline and the alkaline-earth cations may be of great importance to explain the behaviour of natural clay samples where mixtures of the two types of interlayer cation are present and also provides valuable information on the cation exchange occurring in the swelling clays.

Structure-dependent interatomic dispersion coefficients in oxides with maximally localized Wannier functions

The interatomic C(6) dispersion coefficients in crystalline and amorphous SiO(2) and ZrO(2) structures were obtained with the approach proposed by Silvestrelli (2008 Phys. Rev. Lett. 100 053002) and based on the use of maximally localized Wannier functions (MLWFs) for partitioning the electron density. Localization of Wannier functions close to the nuclei in oxide systems makes it possible to assign the MLWFs to the atoms in an unambiguous way and then to compute the C(6) coefficients in an atom pairwise manner. A modification of the method is suggested in which the MLWFs are condensed to effective orbitals centred on the atoms and parameters of these effective orbitals are used for computing the interatomic dispersion coefficients. The obtained values of the dispersion coefficients were found to vary not only from one oxide to another, but also between different modifications of the same compound. The oxygen-oxygen coefficient C6(OO) reveals the largest variation and its value in ZrO(2) structures is twice as large as that in SiO(2) ones. Atomic characteristics obtained in the frame of the effective orbital method, such as the self-atom dispersion coefficient, and the oxide ion polarizability were found to correlate with the metal-oxygen bond length and the oxygen coordination number in the systems. This behaviour is attributed to the confinement of electrons by the electrostatic potential. The values of the coefficient and of the polarizability were related to charges of the oxygen atoms. In all studied systems the oxygen atoms having larger absolute values of charge were found to be less polarizable because of a stronger confinement effect. The obtained results can be used in the development of polarizable force fields for the atomistic modelling of oxide materials.