Elisabete Sónia do Carmo Ferreira

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.