Structure and self-diffusivity of mixed-cation electrolytes between neutral and charged graphene sheets.

Graphene-based applications, such as supercapacitors or capacitive deionization, take place in an aqueous environment, and they benefit from molecular-level insights into the behavior of aqueous electrolyte solutions in single-digit graphene nanopores with a size comparable to a few molecular diameters. Under single-digit graphene nanoconfinement (smallest dimension <2 nm), water and ions behave drastically different than in the bulk. Most aqueous electrolytes in the graphene-based applications as well as in nature contain a mix of electrolytes.

Structure and self-diffusivity of alkali-halide electrolytes in neutral and charged graphene nanochannels.

Understanding the microscopic behaviour of aqueous electrolyte solutions in graphene-based ultrathin nanochannels is important in nanofluidic applications such as water purification, fuel cells, and molecular sensing. Under extreme confinement (<2 nm), the properties of water and ions differ drastically from those in the bulk phase. We studied the structural and diffusion behaviour of prototypical aqueous solutions of electrolytes (LiCl, NaCl, and KCl) confined in both neutral and positively-, and negatively-charged graphene nanochannels.

Molecular dynamics of preferential adsorption in mixed alkali-halide electrolytes at graphene electrodes.

Understanding the microscopic behavior of aqueous electrolyte solutions in contact with graphene and related carbon surfaces is important in electrochemical technologies, such as capacitive deionization or supercapacitors. In this work, we focus on preferential adsorption of ions in mixed alkali-halide electrolytes containing different fractions of Li/Na or Li/K and/or Na/K cations with Cl anions dissolved in water. We performed molecular dynamics simulations of the solutions in contact with both neutral and positively and negatively charged graphene surfaces under ambient conditions, using the effectively polarizable force field.

Molecular Force Field Development for Aqueous Electrolytes: 2. Polarizable Models Incorporating Crystalline Chemical Potential and Their Accurate Simulations of Halite, Hydrohalite, Aqueous Solutions of NaCl, and Solubility.

The current state-of-the-art force fields (FFs) for Na and Cl ions are not capable of simultaneously predicting the thermodynamic properties of the aqueous solution and the crystalline phase. This is primarily due to an oversimplification of the interaction models used but partially also due to the insufficient parametrization of the FFs. We have devised a straightforward and simple parametrization procedure for determining the ion-ion interaction parameters in complex molecular models of NaCl electrolytes which involves fitting the density, lattice energy, and chemical potential of crystalline NaCl at ambient conditions.

Chemical potentials of alkaline earth metal halide aqueous electrolytes and solubility of their hydrates by molecular simulation: Application to CaCl, antarcticite, and sinjarite.

We present a molecular-level simulation study of CaCl in water and crystalline hydrates formed by CaCl at ambient (298.15 K, 1 bar) conditions and at a high-temperature high-pressure state (365 K, 275 bars) typical of hydraulic fracturing conditions in natural-gas extraction, at which experimental properties are poorly characterized. We focus on simulations of chemical potentials in both solution and crystalline phases and on the salt solubility, the first time to our knowledge that such properties have been investigated by molecular simulation for divalent aqueous electrolytes.

Modelling aqueous solubility of sodium chloride in clays at thermodynamic conditions of hydraulic fracturing by molecular simulations.

To address the high salinity of flow-back water during hydraulic fracturing, we have studied the equilibrium partitioning of NaCl and water between the bulk phase and clay pores. In shale rocks, such a partitioning can occur between fractures with a bulk-like phase and clay pores. We use an advanced Grand Canonical Monte Carlo (GCMC) technique based on fractional exchanges of dissolved ions and water molecules.

Thermodynamics of supersaturated steam: Molecular simulation results.

Supersaturated steam modeled by the Gaussian charge polarizable model [P. Paricaud, M. Předota, and A.

Chemical Potentials, Activity Coefficients, and Solubility in Aqueous NaCl Solutions: Prediction by Polarizable Force Fields.

We describe a computationally efficient molecular simulation methodology for calculating the concentration dependence of the chemical potentials of both solute and solvent in aqueous electrolyte solutions, based on simulations of the salt chemical potential alone. We use our approach to study the predictions for aqueous NaCl solutions at ambient conditions of these properties by the recently developed polarizable force fields (FFs) AH/BK3 of Kiss and Baranyai (J. Chem.

Electrolyte pore/solution partitioning by expanded grand canonical ensemble Monte Carlo simulation.

Using a newly developed grand canonical Monte Carlo approach based on fractional exchanges of dissolved ions and water molecules, we studied equilibrium partitioning of both components between laterally extended apolar confinements and surrounding electrolyte solution. Accurate calculations of the Hamiltonian and tensorial pressure components at anisotropic conditions in the pore required the development of a novel algorithm for a self-consistent correction of nonelectrostatic cut-off effects. At pore widths above the kinetic threshold to capillary evaporation, the molality of the salt inside the confinement grows in parallel with that of the bulk phase, but presents a nonuniform width-dependence, being depleted at some and elevated at other separations.

Vapor-liquid equilibrium and polarization behavior of the GCP water model: Gaussian charge-on-spring versus dipole self-consistent field approaches to induced polarization.

We developed the Gaussian charge-on-spring (GCOS) version of the original self-consistent field implementation of the Gaussian Charge Polarizable water model and test its accuracy to represent the polarization behavior of the original model involving smeared charges and induced dipole moments. For that purpose we adapted the recently proposed multiple-particle-move (MPM) within the Gibbs and isochoric-isothermal ensembles Monte Carlo methods for the efficient simulation of polarizable fluids. We assessed the accuracy of the GCOS representation by a direct comparison of the resulting vapor-liquid phase envelope, microstructure, and relevant microscopic descriptors of water polarization along the orthobaric curve against the corresponding quantities from the actual GCP water model.

Thermodynamics of small alkali metal halide cluster ions: comparison of classical molecular simulations with experiment and quantum chemistry.

We evaluate the ability of selected classical molecular models to describe the thermodynamic and structural aspects of gas-phase hydration of alkali metal halide ions and the formation of small water clusters. To understand the effect of many-body interactions (polarization) and charge penetration effects on the accuracy of a force field, we perform Monte Carlo simulations with three rigid water models using different functional forms to account for these effects: (i) point charge nonpolarizable SPC/E, (ii) Drude point charge polarizable SWM4-DP, and (iii) Drude Gaussian charge polarizable BK3. Model predictions are compared with experimental Gibbs free energies and enthalpies of ion hydration, and with microscopic structural properties obtained from quantum DFT calculations.

Molecular Force Field Development for Aqueous Electrolytes: 1. Incorporating Appropriate Experimental Data and the Inadequacy of Simple Electrolyte Force Fields Based on Lennard-Jones and Point Charge Interactions with Lorentz-Berthelot Rules.

It is known that none of the available simple molecular interaction models of aqueous electrolytes based on SPC/E water and their associated force fields are able to reproduce the concentration dependence of important thermodynamic properties of even the simplest electrolyte, NaCl, at ambient conditions over the entire experimentally accessible concentration range [ Mouc̆ka , F. ; Nezbeda , I. ; Smith , W.

Molecular simulation of aqueous electrolytes: water chemical potential results and Gibbs-Duhem equation consistency tests.

This paper deals with molecular simulation of the chemical potentials in aqueous electrolyte solutions for the water solvent and its relationship to chemical potential simulation results for the electrolyte solute. We use the Gibbs-Duhem equation linking the concentration dependence of these quantities to test the thermodynamic consistency of separate calculations of each quantity. We consider aqueous NaCl solutions at ambient conditions, using the standard SPC/E force field for water and the Joung-Cheatham force field for the electrolyte.

Molecular force fields for aqueous electrolytes: SPC/E-compatible charged LJ sphere models and their limitations.

Thirteen of the most common aqueous NaCl solution force fields based on the SPC/E water solvent are examined with respect to their prediction at ambient conditions of the concentration dependence of the total electrolyte chemical potential and the solution density. We also calculate the salt solubility and the chemical potential and density of the NaCl crystalline solid. We obtain the solution chemical potential in a computationally efficient manner using our recently developed Osmotic Ensemble Monte Carlo method [F.

Molecular simulation of aqueous electrolyte solubility. 3. Alkali-halide salts and their mixtures in water and in hydrochloric acid.

We extend the osmotic ensemble Monte Carlo (OEMC) molecular simulation method (Moučka et al. J. Phys Chem.

Molecular simulation of aqueous electrolyte solubility. 2. Osmotic ensemble Monte Carlo methodology for free energy and solubility calculations and application to NaCl.

We present a new and computationally efficient methodology using osmotic ensemble Monte Carlo (OEMC) simulation to calculate chemical potential-concentration curves and the solubility of aqueous electrolytes. The method avoids calculations for the solid phase, incorporating readily available data from thermochemical tables that are based on well-defined reference states. It performs simulations of the aqueous solution at a fixed number of water molecules, pressure, temperature, and specified overall electrolyte chemical potential.

Novel perturbation approach for the structure factor of the attractive hard-core Yukawa fluid.

A novel perturbation approach for the structure factor S(k) of the Lennard-Jones-type Yukawa fluid with z=1.8 is presented. An approach is based on a new reference system, that is, the short-range Yukawa model with z0>z=1.

Efficient multiparticle sampling in Monte Carlo simulations on fluids: application to polarizable models.

A novel Monte Carlo simulation scheme based on biased simultaneous displacements of all particles of the system has been developed. The method is particularly suited for systems with nonadditive interactions and its efficiency is demonstrated by its implementation for the polarizable Stockmayer fluid. Performance of the method is compared with both the standard one-particle move method and an unbiased multiparticle scheme by computing the mean squared displacements, rotation relaxation, and the speed of equilibration (translational order parameter).

Detection and characterization of structural changes in the hard-disk fluid under freezing and melting conditions.

The fluid of two-dimensional hard disks is investigated over a range of densities by Monte Carlo simulations in order to detect and characterize structural changes which take place when the condition of freezing and melting is approached. A novel method is proposed based on the use of the Voronoi tessellation and a certain shape factor which turns out to be a clear indicator of the presence of different underlying substructures (domains). Close to the freezing condition the probability distribution of the shape factor develops a second distinct maximum corresponding to a predominant presence of near-regular hexagons, whereas the original peak, having its origin primarily in pentagons and distorted hexagons, diminishes and disappears at melting density.