Abstract.This work describes the behaviour of water molecules in 1-butyl-3-methylimidazolium tetrafluoroborate ionic liquid under nanoconfinement between graphene sheets. By means of molecular dynamics simulations, an adsorption of water molecules at the graphene surface is studied. A depletion of water molecules in the vicinity of the neutral and negatively charged graphene surfaces and their adsorption at the positively charged surface are observed in line with the preferential hydration of the ionic liquid anions. The findings are appropriately described using a two-level statistical model. The confinement effect on the structure and dynamics of the mixtures is thoroughly analyzed using the density and the potential of mean force profiles, as well as by the vibrational densities of states of water molecules near the graphene surface. The orientation of water molecules and the water-induced structural transitions in the layer closest to the graphene surface are also discussed.

We report the structure and charge transport properties of a novel solid-state proton conductor obtained by acid–base chemistry via proton transfer from 12-tungstophosphoric acid to imidazole. The resulting material (henceforth named Imid3WP) is a solid salt hydrate that, at room temperature, includes four water molecules per structural unit. To our knowledge, this is the first attempt to tune the properties of a heteropolyacid-based solid-state proton conductor by means of a mixture of water and imidazole, interpolating between water-based and ionic liquid-based proton conductors of high thermal and electrochemical stability. The proton conductivity of Imid3WP·4H2O measured at truly anhydrous conditions reads 0.8 × 10–6 S cm–1 at 322 K, which is higher than the conductivity reported for any other related salt hydrate, despite the lower hydration. In the pseudoanhydrous state, that is, for Imid3WP·2H2O, the proton conductivity is still remarkable and, judging from the low activation energy (E a = 0.26 eV), attributed to structural diffusion of protons. From complementary X-ray diffraction data, vibrational spectroscopy, and solid-state NMR experiments, the local structure of this salt hydrate was resolved, with imidazolium cations preferably orienting flat on the surface of the tungstophosphate anions, thus achieving a densely packed solid material, and water molecules of hydration that establish extremely strong hydrogen bonds. Computational results confirm these structural details and also evidence that the path of lowest energy for the proton transfer involves primarily imidazole and water molecules, while the proximate Keggin anion contributes with reducing the energy barrier for this particular pathway.

We perform molecular dynamics simulations of mixtures of a prototypical protic ionic liquid, ethylammonium nitrate, with lithium or magnesium nitrate (LiNO3/Mg(NO3)2) confined between two graphene walls. The structure of the system is analyzed by means of ionic density profiles, angular orientations of ethylammonium cations close to the wall and the lateral structure of the first layer close to the graphene wall. All these results are compared to those of the corresponding aprotic ionic liquid systems, analyzing the influence of the graphene wall charge in the structure of the protic and aprotic mixtures. Moreover, vibrational densities of states are calculated for the salt cations close to the walls. Finally, we investigate the structure of the mixture with Li salt near the interface using ab initio density functional theory, and the results are compared with those obtained by classical molecular dynamics simulations.

The structural and dynamical properties of bulk mixtures of long-chained primary and secondary alcohols (propanol, butanol, and 2-pentanol) with protic ionic liquids (ethylammonium and butylammonium nitrate) were studied by means of molecular dynamics simulations and small angle X-ray scattering (SAXS). Changes in the structure with the alcohol concentration and with the alkyl chain length of the alcohol moieties were found, showing variations in the radial distribution function and in the number of hydrogen bonds in the bulk liquids. Moreover, the structural behaviour of the studied mixtures is further clarified with the spatial distribution functions. The global picture in the local scale is in good agreement with the nanostructured solvation paradigm [T. Méndez-Morales et al. Phys. Chem. B 118, 761 (2014)], according to which alcohols are accommodated into the hydrogen bonds' network of the ionic liquid instead of forming clusters in the bulk. Indeed, our study reveals that the alcohol molecules are placed with their polar heads at the interfaces between polar and nonpolar nanodomains in the ionic liquid, with their alkyl chains inside the nonpolar organic nanodomains. The influence of alcohol chain length in the single-particle dynamics of the mixtures is also reported calculating the velocity autocorrelation function and vibrational densities of states of the different species in the ionic liquid-alcohol mixtures, and a weak caging effect for the ethylammonium cations independent of the chain size of the alcohols was found. However, the SAXS data collected for the studied mixtures show an excess of the scattering intensities which indicates that there are also some structural heterogeneities at the nanoscale.

We perform molecular dynamics simulations of ionic liquids confined between graphene walls under a large variety of conditions (pure ionic liquids, mixtures with water and alcohols, mixtures with lithium salts and defective graphene walls). Our results show that the formation of striped and hexagonal patterns in the Stern layer can be considered as a general feature of ionic liquids at electrochemical interfaces, the transition between patterns being controlled by the net balance of charge in the innermost layer of adsorbed molecules. This explains previously reported experimental and computational results and, for the first time, why these pattern changes are triggered by any perturbation of the charge density at the innermost layer of the electric double layer (voltage and composition changes, and vacancies at the electrode walls, among others), which may help tuning electrode-ionic liquid interfaces. Using Monte Carlo simulations we show that such structures can be reproduced by a simple two-dimensional lattice model with only nearest-neighbour interactions, governed by highly screened ionic interactions and short-range and excluded volume interactions. We also show that the results of our simulations are consistent with those inferred from the Landau-Brazovskii theory of pattern formation in self-assembling systems. The presence of these patterns at the ionic liquid graphene-electrode interfaces may have a strong impact on the process of ionic transfer from the bulk mixtures to the electrodes, on the differential capacitance of the electrode-electrolyte double layer or on the rates of redox reactions at the electrodes, among other physicochemical properties, and is therefore an effect of great technological interest.

We report the first MD simulations of the 3D structure of the electric double layer of ionic liquid–alcohol mixtures.

A new method for switching between structures consisting of equivalent discrete and flexible objects with different particle representation and object configuration, including different resolution levels (number of particles per object), is reported. The method is fully general since it does not require any extra code nor additional database elements for new systems. It is based on a Monte Carlo sampling of the configurational space for each object type of the target system. The sampling is controlled by a Metropolis acceptance criterion of movements (translations, rotations, and relative deformations of the object configuration) that uses the generalized distance between the sets of particles at both representations. For Gaussian distributed distances, such a minimization procedure is equivalent to an optimization of χ in a maximum likelihood method. This provides sound statistical support since the method leads to the most probable configuration of the system at each representation. The configurations obtained in this way are then used to create resolution exchange maps for each object type, which allows the extrapolation of the conversion to every object configuration throughout the whole system. As an example, the method is here tested with several molecular dynamics simulated systems (ionic liquids, cyclodextrins, cell-penetrating peptides, cyclic peptides, lipid bilayers, vesicles, heterogeneous organic molecules, DNA, and solvated proteins) for different resolution force fields (GROMOS, AMBER, OPLS, MARTINI) using GROMACS. In this context, the method can be applied to map structures described by any other pair of force fields, as well as to homogeneous and heterogeneous systems with many different molecules. The method is proved to be highly efficient since the time required for the mapping is practically independent of the number of molecules in the target system.

No abstract