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. We study several prototypical aqueous mixed alkali-chloride electrolytes containing an equimolar fraction of Li/Na, Li/K, or Na/K cations confined between neutral and positively or negatively charged parallel graphene sheets. The strong hydration shell of small Li+ vs a larger Na+ or large K+ with weaker or weak hydration shells affects the interplay between the ions's propensity to hydrate or dehydrate under the graphene nanoconfinement and the strength of the ion-graphene interactions mediated by confinement-induced layered water. We perform molecular dynamics simulations of the confined mixed-cation electrolytes using the effectively polarizable force field for electrolyte-graphene systems and focused on a relation between the electrochemical adsorption and structural properties of the water molecules and ions and their diffusion behavior. The simulations show that the one-layer nanoslits have the biggest impact on the ions' adsorption and the water and ions' diffusion. The positively charged one-layer nanoslits only allow for Cl- adsorption and strengthen the intermolecular bonding, which along with the ultrathin confinement substantially reduces the water and Cl- diffusion. In contrast, the negatively charged one-layer nanoslits only allow for adsorption of weakly hydrated Na+ or K+ and substantially break up the non-covalent bond network, which leads to the enhancement of the water and Na+ or K+ diffusion up to or even above the bulk diffusion. In wider nanoslits, cations adsorb closer to the graphene surfaces than Cl-'s with preferential adsorption of a weakly hydrated cation over a strongly hydrated cation. The positive graphene charge has an intuitive effect on the adsorption of weakly hydrated Na+'s or K+'s and Cl-'s and a counterintuitive effect on the adsorption of strongly hydrated Li+'s. On the other hand, the negative surface charge has an intuitive effect on the adsorption of both types of cations and only mild intuitive or counterintuitive effects on the Cl- adsorption. The diffusion of water molecules and ions confined in the wider nanoslits is reduced with respect to the bulk diffusion, more for the positive graphene charge, which strengthened the intermolecular bonding, and less for the negative surface charge, which weakened the non-covalent bond network.

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. We performed molecular dynamics simulations of the solutions in the nanochannels with either one, two- or three-layer water structures using the effectively polarisable force field for graphene. We analysed the structure and intermolecular bond network of the confined solutions along with their relation to the self-diffusivity of water and ions. The simulations show that Na and K cations can more easily rearrange their solvation shells under the graphene nanoconfinement and adsorb on the graphene surfaces or dissolve in the confinement-induced layered water than the Li cation. The negative surface charge together with the presence of ions orient water molecules with hydrogens towards the graphene surfaces, which in turn weakens the intermolecular bond network. The one-layer nanochannels have the biggest effect on the water structure and intermolecular bonding as well as on the adsorption of ions with only co-ions entering these nanochannels. The self-diffusivity of confined water is strongly reduced with respect to the bulk water and decreases with diminishing nanochannel heights except for the negatively-charged one-layer nanochannel. The self-diffusivity of ions also decreases with the reducing the nanochannel heights except for the self-diffusivity of cations in the negatively-charged one-layer nanochannel, evidencing cooperative diffusion of confined water and ions. Due to the significant break-up of the intermolecular bond network in the negatively-charged one-layer nanochannel, self-diffusion coefficients of water and cations exceed those for the two- and three-layer nanochannels and become comparable to the bulk values.

Adsorption of amphiphilic grafted polymers as polymer corrosion inhibitors: insights from mesoscopic simulations.

The homogeneous covering of amphiphillic polymer molecules onto metallic surfaces is of great importance for corrosion inhibitor applications. Lyophillic side chains grafted onto a lyophobic backbone act as anchors that allow the molecule to absorb at the metallic surface preventing the exposure with the solvent. Coarse-grained simulations are used to study the sorption and conformation behaviour of amphiphillic grafted polymers for corrosion inhibition. The backbone insolubility is found to play a key role in the sorption and conformation behaviour in the dilute limit. For finite concentrations, moderate backbone solubility and moderate molecule concentrations achieve optimal surface coverage, while highly a lyophobic backbone leads to bulk-like structures as a consequence of aggregation.

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. The simulations show that large ions are often intuitively attracted to oppositely charged electrodes. In contrast, the adsorption behavior of small ions tends to be counterintuitive. In mixed-cation solutions, one of the cations always supports the adsorption of the other cation, while the other cation weakens the adsorption of the first cation. In mixed-cation solutions containing large and small cations simultaneously, adsorption of the larger cations varies dramatically with the electrode charge in an intuitive way, while adsorption of the smaller cations changes oppositely, i.e., in a counterintuitive way. For (Li/K)Cl mixed-cation solutions, these effects allow the control of Li+ adsorption by varying the electrode charge, whereas, for LiCl single-salt solutions, Li+ adsorption is nearly independent of the electrode charge. We rationalize this cation-cation lever effect as a result of a competition between three driving forces: (i) direct graphene-ion interactions, (ii) the strong tendency of the solutions to saturate the network of non-covalent intermolecular bonds, and (iii) the tendency to suppress local charge accumulation in any region larger than typical interparticle distances. We analyze the driving forces in detail using a general method for intermolecular bonding based on spatial distribution functions and different contributions to the total charge density profiles. The analysis helps to predict whether an ion is more affected by each of the three driving forces, depending on the strength of the ion solvation shells and the compatibility between the contributions of the charge density profiles due to the ion and water molecules. This approach is general and can also be applied to other solutions under different thermodynamic conditions.

Structural properties of cationic surfactant-fatty alcohol bilayers: insights from dissipative particle dynamics.

Bilayers, self-assembled by cationic surfactants and fatty alcohols in water, are the basic units of lamellar gel networks - creamy formulations extensively used in cosmetics and pharmaceutics. Mesoscopic modelling and study of the bilayers formed by single- or double-tail cationic surfactants (CTAC or DHDAC), and fatty alcohols (FAs) in the lamellar fluid and gel phases were employed. Fatty alcohols with alkyl tail equal to or greater than the surfactant alkyl tail, i.e., C16FA or C18FA and C22FA, were considered. A model formulation was explored with the FA concentration greater than that of the surfactant and the structure of the fluid and gel bilayers in tensionless state characterised via the density profiles across the bilayers, orientational order parameters of the surfactant and FA chains, intrinsic analysis of the bilayer interfaces, and bending rigidity. The intrinsic analysis allows identification and quantification of the coexistence of the interdigitated and non-interdigitated phases present within the gel bilayers. The FA chains were found to conform the primary scaffolding of the bilayers while the surfactant chains tessellate bilayer monolayers from their water-hydrophobic interface. Further, the overlap of the FA chains from the apposed monolayers of the fluid bilayers rises with increasing FA length. Finally, the prevalence of the non-interdigitated phase over the interdigitated phase within the gel bilayers becomes enhanced upon the FA length increase with a preference of the surfactant chains to reside in the non-interdigitated phase rather than the interdigitated phase.

Generalised dissipative particle dynamics with energy conservation: density- and temperature-dependent potentials.

We present a generalised, energy-conserving dissipative particle dynamics (DPDE) method appropriate for the non-isothermal simulation of particle interaction force fields that are both density- and temperature-dependent. A detailed derivation is formulated in a bottom-up manner by considering the thermodynamics of small systems with the appropriate consideration of the fluctuations. Connected to the local volume is a local density and corresponding local pressure, which is determined from an equation-of-state based force field that depends also on a particle temperature. Compared to the original DPDE method, the formulation of the generalised DPDE method requires a change in the independent variable from the particle internal energy to the particle entropy. As part of the re-formulation, the terms dressed particle entropy and the corresponding dressed particle temperature are introduced, which depict the many-body contributions in the local volume. The generalised DPDE method has similarities to the energy form of the smoothed dissipative particle dynamics method, yet fundamental differences exist, which are described in the manuscript. The basic dynamic equations are presented along with practical considerations for implementing the generalised DPDE method, including a numerical integration scheme based on the Shardlow-like splitting algorithm. Demonstrations and validation tests are performed using analytical equation-of-states for the van der Waals and Lennard-Jones fluids. Particle probability distributions are analysed, where excellent agreement with theoretical estimates is demonstrated. As further validation of the generalised DPDE method, both equilibrium and non-equilibrium simulation scenarios are considered, including adiabatic flash heating response and vapour-liquid phase separation.

Dissipative particle dynamics with reactions: Application to RDX decomposition.

We present a general, flexible framework for a constant-energy variant of the dissipative particle dynamics method that allows chemical reactions (DPD-RX). In our DPD-RX approach, reaction progress variables are assigned to each particle that monitor the time evolution of an extent-of-reaction associated with the prescribed reaction mechanisms and kinetics assumed to occur within the particle, where chemistry can be modeled using complex or reduced reaction mechanisms. We demonstrate our DPD-RX method by considering thermally initiated unimolecular decomposition of the energetic material, cyclotrimethylene trinitramine (RDX), into a molecular gas mixture. Studies are performed to demonstrate the effect of a spatially averaged particle internal temperature and a local reaction volume term in the chemical kinetics expressions, where both provide implicit mechanisms for capturing condensed phase reactivity. We also present an analysis of the expansion of the product gas mixture during decomposition. Finally, a discussion of other potential applications and extensions of the DPD-RX method is given.

Controlled Transport of Flexible Polymers in Slit and Cylindrical Pores Coated with Polymer Brushes: Insight from Dissipative Particle Dynamics.

We use dissipative particle dynamics to simulate the controlled transport of flexible polymers through coated slit and cylindrical pores. Pores are coated inside with solvent-sensitive polymer brushes. Stretch-to-collapse transition then controls the permeability of the coated pores. We change the solvent quality with respect to the polymer brushes and study the flow of flexible polymers through the pores. We show that stretched brush chains close the pores and compress the polymers in the centre of the pores. The collapsed brush chains relieve compression and rapid change in permeability is observed. In open pore state, polymers partially accommodate on the brush layer and partially migrate towards the center of the pores, where the flux is maximal. Finally, we observe that polymers tend to align in the direction of the flow.

Adsorption of amphiphilic graft copolymers in solvents selective for the grafts on a lyophobic surface: a coarse-grained simulation study.

The sorption of graft copolymers on surfaces attractive only for the backbone and its effect on the conformational behavior of adsorbed/desorbed chains in solvents good for the grafts and poor for the backbone was studied by coarse-grained computer simulations. It was found that the sorption and conformational behavior are very complex and are results of an intricate interplay of solvent quality (polymer-solvent interactions) and solvent strength (polymer-surface vs. solvent-surface interactions). Increasing grafting density and length of grafts protect the backbone against adsorption, but the behavior is non-trivial. A decrease in solvent quality promotes the adsorption, because it lowers the overall solubility, but the backbone collapses and the probability of backbone-surface contacts decreases, which simultaneously hinders the adsorption. The results of simulations are presented in the form of phase diagrams depicting the decisive features of the conformational and sorption behavior.

Concentrated aqueous sodium chloride solution in clays at thermodynamic conditions of hydraulic fracturing: Insight from molecular dynamics simulations.

To address a high salinity of flow-back water during hydraulic fracturing, we use molecular dynamics (MD) simulations and study the thermodynamics, structure, and diffusion of concentrated aqueous salt solution in clay nanopores. The concentrated solution results from the dissolution of a cubic NaCl nanocrystal, immersed in an aqueous NaCl solution of varying salt concentration and confined in clay pores of a width comparable to the crystal size. The size of the nanocrystal equals to about 18 Å which is above a critical nucleus size. We consider a typical shale gas reservoir condition of 365 K and 275 bar, and we represent the clay pores as pyrophyllite and Na-montmorillonite (Na-MMT) slits. We employ the Extended Simple Point Charge (SPC/E) model for water, Joung-Cheatham model for ions, and CLAYFF for the slit walls. We impose the pressure in the normal direction and the resulting slit width varies from about 20 to 25 Å when the salt concentration in the surrounding solution increased from zero to an oversaturated value. By varying the salt concentration, we observe two scenarios. First, the crystal dissolves and its dissolution time increases with increasing salt concentration. We describe the dissolution process in terms of the number of ions in the crystal, and the crystal size and shape. Second, when the salt concentration reaches a system solubility limit, the crystal grows and attains a new equilibrium size; the crystal comes into equilibrium with the surrounding saturated solution. After crystal dissolution, we carry out canonical MD simulations for the concentrated solution. We evaluate the hydration energy, density profiles, orientation distributions, hydrogen-bond network, radial distribution functions, and in-plane diffusion of water and ions to provide insight into the microscopic behaviour of the concentrated aqueous sodium chloride solution in interlayer galleries of the slightly hydrophobic pyrophyllite and hydrophilic Na-MMT pores.

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

We present a molecular-level simulation study of CaCl2 in water and crystalline hydrates formed by CaCl2 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. We first extend our osmotic ensemble Monte Carlo simulation technique [F. Moucka et al., J. Phys. Chem. B 115, 7849-7861 (2011)] to such solutions. We then describe and apply new methodology for the simulation of the chemical potentials of the experimentally observed crystalline hydrates at ambient conditions (antarcticite, CaCl2·6H2O) and at high-temperature conditions (sinjarite, CaCl2·2H2O). We implement our methodologies using for both phases the CaCl2 transferable force field (FF) based on simple point charge-extended water developed by Mamatkulov et al. [J. Chem. Phys. 138, 024505 (2013)], based on training sets involving single-ion and ion-pair low-concentration solution properties at near-ambient conditions. We find that simulations of the solution chemical potentials at high concentrations are somewhat problematic, exhibiting densities diverging from experimental values and accompanied by dramatically decreasing particle mobility. For the solid phases, the sinjarite crystalline lattice differs from experiment only slightly, whereas the simulations of antarcticite completely fail, due to instability of the crystalline lattice. The FF thus only successfully yields the sinjarite solubility, but its value m = 8.0(7) mol kg-1H2O lies well below the experimentally observed solubility range at 1 bar pressure of (12m, 15m) in the temperature interval (320 K, 400 K). We conclude that the used FF does not provide a good description of the experimental properties considered and suggest that improvement must take into account the crystalline properties.

Adsorption and Diffusion of C1 to C4 Alkanes in Dual-Porosity Zeolites by Molecular Simulations.

We employ grand canonical Monte Carlo and molecular dynamics simulations to systematically study the adsorption and diffusion of C1 to C4 alkanes in hierarchical ZSM-5 zeolite with micropores (∼1 nm) and mesopores (>2 nm). The zeolite is characterized by a large surface area of active sites on the microporous scale with high permeability and access to the active sites, which arises from the enhanced transport at the mesoporous scale. We model this zeolite as a microporous Na+-exchanged alumino-sillicate zeolite ZSM-5/35 (Si/Al = 35) in which cylindrical mesopores with a diameter of 4 nm have been built by deleting atoms accordingly. We use the TraPPE and Vujic-Lyubartsev force fields along with the Lorentz-Berthelot combining rules to describe adsorbate-adsorbate and adsorbate-adsorbent interactions. The performance of the force fields is assessed by comparing against experimental single-component adsorption isotherms of methane and ethane in microporous ZSM-5/35, which we measured as part of this work. We compare the adsorption isotherms and diffusivities of the adsorbed alkanes in the dual-porosity zeolite with those in microporous ZSM-5/35 and discern the specific behavior at each porosity scale on the overall adsorption, self-diffusion, and transport behavior in zeolites with dual micro/mesoporosities.

Flow and aggregation of rod-like proteins in slit and cylindrical pores coated with polymer brushes: an insight from dissipative particle dynamics.

We use a meso-scale dissipative particle dynamics method to simulate the flow and aggregation of rod-like protein solutions through pores grafted with a solvent-sensitive polymer brush. The coated pores can control protein permeability and aggregation by a stretch-to-collapse conformational transition of the brush polymers in response to changes in the solvent quality. The protein solutions mimic aqueous glycoprotein solutions and proteins are represented as rod-like objects formed by coarse-grain beads. The model further employs two types of beads to represent the existence of cystein-like terminal groups in real glycoproteins and mimic the aggregation of real glycoproteins in aqueous solutions. We vary the solvent quality with respect to the brush chains and study the flow and aggregation of rod-like proteins in the slit and cylindrical pores as the brush polymers undergo the stretch-to-collapse transition. The results show that stretched brush chains close the pore, hamper proteins' flow and promote proteins' aggregation. The collapsed brush chains open the pores for proteins' flow and suppress their aggregation. Therefore, we observe more than a ten-fold reduction in the permeation rate of proteins in both pore geometries. Finally, due to pore confinement, larger proteins' aggregates are formed in the slit pore than in the cylindrical pore, while more pronounced orientation of proteins in the flow direction is seen in the cylindrical pore than in the slit pore.

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. We consider a typical shale gas reservoir condition of a temperature of 365 K and pressure of 275 bar, and we represent clay pores by pyrophyllite and Na-montmorillonite slits of a width ranging from about 7 to 28 Å, covering clay pores from dry clay to clay pores with a bulk-like layer in the middle of the pore. We employ the Joung-Cheatham model for ions, SPC/E model for water and CLAYFF for the clay pores. We first determine the chemical potentials for NaCl and water in the bulk phase using Osmotic Ensemble Monte Carlo simulations. The chemical potentials are then used in GCMC to simulate the adsorption of ions and water molecules in the clay pores, and in turn to predict the salt solubility in confined solutions. Besides the thermodynamic properties, we evaluate the structure and in-plane diffusion of the adsorbed fluids, and ion conductivities.

Highly grafted polystyrene/polyvinylpyridine polymer gold nanoparticles in a good solvent: effects of chain length and composition.

In this work, the structural features of spherical gold nanoparticles (NPs) decorated with highly grafted poly(styrene) (PS), poly(vinylpyridine) (PVP) and PS-PVP diblock copolymer brushes immersed in a good solvent are investigated by means of Dissipative Particle Dynamics (DPD) simulations as a function of grafted chain length and of homopolymer and copolymer chain composition. For NPs grafted either by PS or PVP homopolymer brushes (selected as a proof of concept), good agreement between the Daoud-Cotton theory, experimental evidence, and our DPD simulations is observed in the scaling behavior of single chain properties, especially for longer grafted chains, and in brush thickness prediction. On the other hand, for grafted chain lengths comparable to NP dimensions parabolic-like profiles of the end-monomer distributions are obtained. Furthermore, a region of high concentration of polymer segments is observed in the monomer density distribution for long homopolymers. In the case of copolymer-decorated NPs, the repulsion between PS and PVP blocks is found to substantially influence the radius of gyration and the shape of the end-monomer distribution of the relevant polymer shell. Moreover, for diblock chains, the un-swollen region is observed to be thinner (and, correspondingly, the swollen layer thicker) than that of a NP modified with a homopolymer of the same length. Finally, the lateral segregation of PS and PVP blocks is evidenced by our calculations and a detailed analysis of the corona behavior is reported, thus revealing the key parameters in controlling the surface properties and the response of diblock copolymer modified nanoparticles.

The self-assembly of copolymers with one hydrophobic and one polyelectrolyte block in aqueous media: a dissipative particle dynamics study.

The reversible self-assembly of symmetrical block copolymers consisting of one hydrophobic block and one ionizable polyelectrolyte block of the same length has been studied in aqueous solutions by dissipative particle dynamics simulations. In addition to three standard dissipative particle dynamics forces (conservative soft repulsion, dissipative and stochastic forces), explicit interaction between smeared charges on ions and on ionized polymer beads described by the electrostatic potential with appropriately localized charges was taken into account. The self-assembly and properties of formed core-shell micelles were investigated as functions of the degree of ionization for systems differing in the hydrophobicity of the non-ionized polyelectrolyte block and in the compatibility of the polymer blocks. This study shows that micelles undergo massive dissociation with increasing degree of ionization. The simulation data compare well with the predictions of scaling theories for systems with soluble polyelectrolytes on a semi-quantitative level and broaden the knowledge of systems in poor solvents.

The electrostatic co-assembly in non-stoichiometric aqueous mixtures of copolymers composed of one neutral water-soluble and one polyelectrolyte (either positively or negatively charged) block: a dissipative particle dynamics study.

The electrostatic co-assembly in non-stoichiometric aqueous mixtures of diblock copolymers composed of a neutral water-soluble block and an either positively or negatively charged polyelectrolyte (PE) block has been studied by dissipative particle dynamics (DPD) simulations. The employed DPD variant includes explicit electrostatics and enables the investigation of the role of small ions in the co-assembly. The properties of core-shell associates containing insoluble interpolyelectrolyte complex cores and protective neutral shells were investigated as functions of the ratio of positive-to-negative charges in the system. This ratio was varied by increasing the number of positively charged PE chains of the same length as those of negatively charged chains, and by changing the PE length and charge density. The simulation results show that the associates formed in non-stoichiometric mixtures differ from those formed in stoichiometric mixtures: their association numbers are lower, their cores are charged and a fraction of excess chains remain free in the non-associated state. The study demonstrates the important role of the compatibility of the counterions with the polymer blocks. It simultaneously emphasizes the necessity of including the electrostatic interaction of all the charged species in the DPD computational scheme.

Free-energy calculations using classical molecular simulation: application to the determination of the melting point and chemical potential of a flexible RDX model.

We present an extension of various free-energy methodologies to determine the chemical potential of the solid and liquid phases of a fully-flexible molecule using classical simulation. The methods are applied to the Smith-Bharadwaj atomistic potential representation of cyclotrimethylene trinitramine (RDX), a well-studied energetic material, to accurately determine the solid and liquid phase Gibbs free energies, and the melting point (Tm). We outline an efficient technique to find the absolute chemical potential and melting point of a fully-flexible molecule using one set of simulations to compute the solid absolute chemical potential and one set of simulations to compute the solid-liquid free energy difference. With this combination, only a handful of simulations are needed, whereby the absolute quantities of the chemical potentials are obtained, for use in other property calculations, such as the characterization of crystal polymorphs or the determination of the entropy. Using the LAMMPS molecular simulator, the Frenkel and Ladd and pseudo-supercritical path techniques are adapted to generate 3rd order fits of the solid and liquid chemical potentials. Results yield the thermodynamic melting point Tm = 488.75 K at 1.0 atm. We also validate these calculations and compare this melting point to one obtained from a typical superheated simulation technique.

Surface segregation in a binary mixture of ionic liquids: Comparison between high-resolution RBS measurements and moleculardynamics simulations.

Surface structure of equimolar mixture of 1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide ([C2C1Im][Tf2N]) and 1-ethyl-3-methylimidazolium tetrafluoroborate ([C2C1Im][BF4]) is studied using high-resolution Rutherford backscattering spectroscopy (HRBS) and molecular dynamics (MD) simulations. Both HRBS and MD simulations show enrichment of [Tf2N] in the first molecular layer although the degree of enrichment observed by HRBS is more pronounced than that predicted by the MD simulation. In the subsurface region, MD simulation shows a small depletion of [Tf2N] while HRBS shows a small enrichment here. This discrepancy is partially attributed to the artifact of the MD simulations. Since the number of each ion is fixed in a finite-size simulation box, surface enrichment of particular ion results in its artificial depletion in the subsurface region.

Surface structure of imidazolium-based ionic liquids: Quantitative comparison between simulations and high-resolution RBS measurements.

Elemental depth profiles of 1-alkyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide ([CnMIM][TFSI], n = 4, 6, 8) are measured using high-resolution Rutherford backscattering spectroscopy (HRBS). The profiles are compared with the results of molecular dynamics (MD) simulations. Both MD simulations and HRBS measurements show that the depth profiles deviate from the uniform stoichiometric composition in the surface region, showing preferential orientations of ions at the surface. The MD simulations qualitatively reproduce the observed HRBS profiles but the agreement is not satisfactory. The observed discrepancy is ascribed to the capillary waves. By taking account of the surface roughness induced by the capillary waves, the agreement becomes almost perfect.