Publications

Abstract New molecular dynamics (MD) simulations and experimental data on a deep eutectic solvent, propeline, composed by choline chloride, ChCl, and propylene glycol, PG, in a molar ratio of 1:2 are reported in this work. The experimental physicochemical properties (density, viscosity and self-diffusion coefficients) were used as support in the development of a new OPLS based force field model (FFM) for propeline. Validation of the new force field was established both through measuring physicochemical properties over a range of temperatures (298.15-373.15 K) and by comparison with experimental and simulated data of ethaline (ChCl:ethylene glycol, at a molar ration of 1:2). Classical MD simulations using the new FFM led to good agreement between experimental and simulated data. Structural properties, namely radial and spatial distribution functions, coordination numbers, and hydrogen bonding were analyzed. Moreover, it was found that the interactions between the anion, Cl-, and the hydrogen bond donor (HBD) form a network that is immutable with increasing temperature. The higher prevalence of anion-HBD hydrogen bonds is likely the major reason for the relatively high viscosity of propeline.

Abstract Despite a growing number of research reports on neat room temperature ionic liquids (RTILs) and their mixtures with molecular solvents in recent years, understanding and rationalising of such systems is still a challenge. In this work, we performed a classical molecular dynamics simulation study of the pure components - 1-ethyl-3-methylimidazolium thiocyanate (C2C1 imSCN), methanol, and ethanol - and their binary mixtures at room temperature. Thermodynamic (density and heats of vaporization), transport (viscosity and self-diffusion coefficients) and structural (in terms of radial, angular and spatial distributions) properties were analysed. It was found, that with the decrease of RTIL content, the ions self-diffusion coefficients notably increase, reaching higher values in the C2C1 imSCN-MeOH system. Density and viscosity follow the opposite trend, reaching their minimum at lower RTIL mole fraction. Negative deviations of excess molar volume from ideality in the studied mixtures with minima at similar to 0.2-03 mole fraction of RTIL suggest the strongest ion-molecular interactions at this mixture composition. A careful analysis at the molecular level revealed that introducing of alcohols to both systems weakens the inter-ionic H-bonding network, particularly, at low RTIL content. The cation-cation arrangement was found to lose its characteristic above/below orientation in neat RTIL and become disordered at low RTIL content. As to the tail length of the selected alcohols, this was found to have an insignificant effect on the structural properties of the addressed systems.

Abstract Mixing of ionic liquids provides new opportunities for their tuning, enabling the applications of ionic liquid mixtures to expand. At the same time, the genesis of the fundamental properties of ionic liquid mixtures is still poorly understood. In this study we carried out a molecular dynamics simulation of binary mixtures of 1-buthyl-3-methylimidazolium hexafluorophosphate, 1-buthyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, and 1-buthyl-3-methylimidazolium tris(perfluoroethyl)trifluorophosphate ([C(4)mim][PF6] + [C(4)mim][NTf2], [C(4)mim][PF6] + [C(4)mim][FAP], [C(4)mim][FAP] + [C(4)mim][NTf2]) in a wide concentration range at 303.15 K and complemented it with quantum mechanical calculations. Three pure ionic liquids underwent the same kind of analysis for comparison purposes. We found that the addition of the [FAP](-)-anion to a mixture enhances the segregation of non-polar domains and weakens the hydrogen-bond network. The H-bonds in the studied mixtures are rather weak, as follows from QTAIM analysis, with the rarest occurrence for the [FAP](-)-anion. The competition of two anions in the mixtures for the most acidic hydrogen of the 1-butyl-3-methylimidazolium cation is reported. In most of the cases, the smaller anion ([PF6](-) or [NTf2](-)) with stronger charge concentration displaces the bigger one ([NTf2](-) or [FAP](-)) from the preferred coordination site. The existing nano-segregation in some mixtures notably slows down ion diffusion. Our results show that the differences in anion size, shape and nature are the main reasons for nano-segregation and the non-ideal behavior of ionic liquid mixtures.

Abstract In this work we developed a new force field model (FFM) for propylene glycol (PG) based on the OPLS all-atom potential. The OPLS potential was refined using quantum chemical calculations, taking into account the densities and self-diffusion coefficients. The validation of this new FFM was carried out based on a wide range of physicochemical properties, such as density, enthalpy of vaporization, self-diffusion coefficients, isothermal compressibility, surface tension, and shear viscosity. The molecular dynamics (MD) simulations were performed over a large range of temperatures (293.15-373.15 K). The comparison with other force field models, such as OPLS, CHARMM27, and GAFF, revealed a large improvement of the results, allowing a better agreement with experimental data. Specific structural properties (radial distribution functions, hydrogen bonding and spatial distribution functions) were then analyzed in order to support the adequacy of the proposed FFM. Pure propylene glycol forms a continuous phase, displaying no microstructures. It is shown that the developed FFM gives rise to suitable results not only for pure propylene glycol but also for mixtures by testing its behavior for a 50 mol % aqueous propylene glycol solution. Furthermore, it is demonstrated that the addition of water to the PG phase produces a homogeneous solution and that the hydration interactions prevail over the propylene glycol self-association interactions.

Abstract In this work, we combined various parameters found in the literature for the choline cation, chloride anion, and ethylene glycol to set up force field models (FFMs) for a eutectic mixture, namely, ethaline (1:2 choline chloride/ethylene glycol (ChCl:2EG)). The validation of these models was carried out on the basis of physical and chemical properties, such as the density, expansion coefficient, enthalpy of vaporization, self-diffusion coefficients, isothermal compressibility, surface tension, and shear viscosity. After the initial evaluation of the FFMs, a refinement was found necessary and accomplished by taking into account polarization effects in a mean-field manner. This was achieved by rescaling the electrostatic charges of the ions based on partial charges derived from ab initio molecular dynamics (MD) simulations of the bulk system. Classical all-atom MD simulations performed over a large range of temperatures (298.15-373.15 K) using the refined FFMs clearly showed improved results, allowing a better prediction of experimental properties. Specific structural properties (radial distribution functions and hydrogen bonding) were then analyzed in order to support the adequacy of the proposed refinement. The final selected FFM leads to excellent agreement between simulated and experimental data on dynamic and structural properties. Moreover, compared to the previously reported force field model (Perkins, S. L.; Painter, P.; Colina, C. M. Experimental and Computational Studies of Choline Chloride-Based Deep Eutectic Solvents. J. Chem. Eng. Data 2014, 59, 3652-3662), a 10% improvement in simulated transport properties, i.e., self-diffusion coefficients, was achieved. The isothermal compressibility, surface tension, and shear viscosity for ethaline are accessed in MD simulations for the first time.

Abstract The results of a. systematic molecular dynamics study of the interfacial structure between the gold (100) surface and two room-temperature ionic liquids, namely, 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIm][PF6]) and 1-butyl-3-methylimadazolium bis(trifluoromethylsulfonyl)imide ([BMIm][NTf2]), are herein reported. It is found that near an uncharged surface the IL structure differs from its bulk, having an enhanced density extended until the two first layers. Interfacial layering is clearly observed at the gold surface, with a higher effect for the [BMIm][NTf2] IL but a higher packing for [BMIm][PF6]. In both ILs the alkyl side chains are oriented parallel to the interface while the imidazolium rings tend to be parallel to the interface in about 60% of the cases. The presence of the interface has a higher impact on the orientation of the cations than on the chemical properties of the counterion. The surface potential drop across the interface is more pronounced toward a negative value for ([BMIm][PF6]) than for ([BMIm][NTf2]), due to relatively larger local density of the anions for ([BMIm][PF6]) near the gold surface.

Abstract In this work we describe the electrochemical behavior of an ionic liquid formed between choline chloride (ChCl) and chromium chloride hexahidrated (CrCl3 center dot 6H(2)O) in the ratio of 2.5:1 containing 20 wt.% of added ultrapure water, using cyclic voltammetry and chronoamperometry. The reduction of Cr(III) occurs in two steps, Cr(III) to Cr(II) and Cr(II) to Cr(0), respectively. The first step is controlled by diffusion. For the first time a thorough study of the nucleation of chromium from Cr(III) ionic liquid solutions is described. Furthermore for less cathodic potentials there is a diffusion control of nuclei growth. For more cathodic potentials the lattice incorporation of adatoms to the growing nuclei is the limiting process. The additive free bright chromium deposit obtained is formed by semi-spherical nuclei. The effect of hydrogen evolution is only observed for E < - 1.90 V.

Abstract In this work we present a fundamental study of the electrodeposition of tin from Deep Eutectic Solvents (DES) formed by a mixture of choline chloride and different hydrogen bond donors (HBD). Results shows that choline chloride based solvents can be successfully used for the electrodeposition of tin. Furthermore we demonstrate that the choice of hydrogen bond donor does not affect, significantly, the chemistry of tin in solution and we characterize the first stages of tin deposits at glassy-carbon (GC) electrode. The electrochemical characterization of tin deposits is carried out using cyclic voltammetry and chronoamperometry. The comparison of the theoretically and experimentally obtained current transients via dimensionless plots based on Bewick-Fleischman-Thirsk (BFT) theory, Scharifker and Hills (SH) and Scharifker and Mostany (SM) models and a non-linear fitting method showed that tin nucleation on GC surface occurs though a 3D instantaneous process with growth controlled by diffusion.

Abstract Mathematical models under conditions of cyclic staircase voltammetry and electrochemical impedance spectroscopy (EIS), which,consider the kinetic effects due to the complexation reaction by the facilitated transfer of metal ions at polarized interfaces, are presented. Criteria for qualitative recognition of these kinetic effects from the features of simulated cyclic voltammograms are given. In case of the existence of these effects, only the EIS can bring access to the thermodynamics and kinetics of the complexation chemical reaction. Analytical equations for estimating the thermodynamic parameters by such systems under EIS conditions are evaluated. The theoretical results are compared with the experimental results of the facilitated Cu2+ transfer at the polarized water(-)1,2-dichlorethane interface, assisted by two phenanthroline-containing macrocycles. In the experimental case where kinetic effects due to the complexation step exist, we show how elegantly EIS can be used as a too] for estimation of the complexation constant of Cu2+ and 5-oxo-2,8-dithia [9],(2,9)-1,10-phenanthrolinophane (PhenOS(2)).