Publications

Abstract Heat capacity, a crucial physical property for chemical processes, is often understudied in Deep Eutectic Solvents (DESs), which in turn are promising green alternatives to environmentally hazardous conventional solvents. This work addresses this gap by developing a machine learning model to predict DES heat capacity and identify key structural features influencing it. We employed a dataset of 530 DESs with corresponding experimental heat capacity values. Quantum-chemical COSMO-RS-based descriptors, capturing detailed information about DES structures, were calculated for each data point. Various machine learning algorithms, namely k-Nearest Neighbours (kNN), Random Forests (RF), Neural Network Multilayer Perceptron (MLP), and Support Vector Machines (SVM) were explored alongside a linear model (Multiple Linear Regression, MLR). Hyperparameter optimisation ensured all models were fine-tuned for optimal performance. The most successful model, based on the MLP technique, achieved remarkably low Average Absolute Relative Deviation (AARD) values of 0.500 % and 3.999 % for the training and test sets, respectively. This signifies a significant improvement in prediction accuracy compared to traditional methods. Furthermore, by applying a SHapley Additive exPlanations (SHAP) analysis, we identified the most crucial structural factors within DES components that govern their heat capacity. This comprehensive investigation offers valuable insights that can pave the way for an efficient design of novel DESs in the future. © 2024 The Author(s)

Abstract This study investigates the structural and transport properties of SDS, CTAB, and SB3-12 micelles in three deep eutectic solvents (DESs), Ethaline, Glyceline, and Reline, using molecular dynamics (MD) simulations. The influence of solvent composition on micelle morphology, interactions, and dynamics was explored, revealing key differences driven by the DES environment. Structural analyses, including eccentricity and radius of gyration, demonstrated that micelle shape and compactness vary significantly depending on the solvent. In Ethaline and Reline, larger micelles showed significant deviations from spherical shapes, while micelles in Glyceline became more spherical and compact, particularly those formed by SB3-12. Radial distribution functions highlighted different levels of micelle-solvent interactions, with SDS showing strong interactions with HBD components and SB3-12 exhibiting prominent self-interaction. According to hydrogen bonding analysis, micelles slightly disrupt the DES hydrogen bond network, with SB3-12 establishing the most significant hydrogen bond connections. The transport property analysis revealed that larger micelles have lower diffusion coefficients, whereas smaller micelles enhance DESs' component mobility. These findings advance the understanding of micelle behavior in DESs and also help in the optimization of DES-surfactant systems for applications such as electrodeposition, nanomaterial templating, and drug delivery. Future research will focus on surfactant interactions with surfaces to further improve these applications.

Abstract Magnetic ionic liquids (MILs) stand out as a remarkable subclass of ionic liquids (ILs), combining the desirable features of traditional ILs with the unique ability to respond to external magnetic fields. The incorporation of paramagnetic species into their structures endows them with additional attractive features, including thermochromic behavior and luminescence. These exceptional properties position MILs as highly promising materials for diverse applications, such as gas capture, DNA extractions, and sensing technologies. The present Review synthesizes key experimental findings, offering insights into the structural, thermal, magnetic, and optical properties across various MIL families. Special emphasis is placed on unraveling the influence of different paramagnetic species on MILs' behavior and functionality. Additionally, the Review highlights recent advancements in computational approaches applied to MIL research. By leveraging molecular dynamics (MD) simulations and density functional theory (DFT) calculations, these computational techniques have provided invaluable insights into the underlying mechanisms governing MILs' behavior, facilitating accurate property predictions. In conclusion, this Review provides a comprehensive overview of the current state of research on MILs, showcasing their special properties and potential applications while highlighting the indispensable role of computational methods in unraveling the complexities of these intriguing materials. The Review concludes with a forward-looking perspective on the future directions of research in the field of magnetic ionic liquids.

Abstract Electrolytes play a crucial role in enhancing the performance of energy storage devices, including batteries and supercapacitors. However, traditional electrolytes, such as aqueous solutions, organic solvents, and ionic liquids, exhibit inherent limitations and challenges. Deep eutectic solvents have recently emerged as promising alternatives due to their environmentally friendly nature and favorable properties. Despite their widespread applications in various domains, their potential as electrolytes remains relatively underexplored. This study investigates three distinct types of deep eutectic solvents derived from different isomers of butanediols combined with choline chloride. Ab initio molecular dynamics simulations are employed to analyze the microstructure of these deep eutectic solvents, focusing on non-covalent electrostatic interactions, hydrogen bonding patterns, and vibrational spectra. The results reveal significant differences in the structural configuration of hydrogen bond acceptors and hydrogen bond donors and their interactions within the deep eutectic solvents. Specifically, the positioning of functional groups in hydrogen bond donors significantly impacts the hydrogen bonding network and the interaction with monoatomic ions. Moreover, the vibrational spectra analysis highlights the existence of hydrogen bonds involving stretching modes of the OH group, as evidenced by redshift deviations. Overall, this study provides valuable insights into the unique features of deep eutectic solvents as potential electrolytes for energy storage applications. The comprehensive analysis of their microstructure and vibrational properties enhances our understanding of deep eutectic solvent utilization and opens avenues for further research in sustainable energy storage.

Abstract Technologies involving a solvent|surface interface, such as nanotechnology, electrochemistry, and energy storage applications, are actively pursuing ecologically responsible and sustainable development practices. In response to this pressing need, deep eutectic solvents have emerged as a promising solution to bridge the gap between technological requirements and environmental concerns. In this work, we present the results of a molecular dynamics simulation study of the interface between a monocrystalline gold surface and the deep eutectic solvent ethaline, where a molar ratio of 1:2 choline chloride:ethylene glycol was used for ethaline. The simulations covered a range of temperatures from 313 K to 343 K and applied charge values ranging from 0 to +/- 24 mu C cm-2. Several key interfacial properties were thoroughly analyzed, including among others, charge density profiles, radial distribution functions, hydrogen bond close contacts, and molecular orientation. Additionally, we examined how the differential capacitance varied upon the applied potential. Our findings reveal that, at neutral surfaces, all components of the solvent are present in the innermost layer, with ethylene glycol molecules being the most prevalent, followed by choline cations and a residual amount of chloride anions. For lower applied charges, this mixed composition at the boundary layer persists, despite the growing accumulation of ionic species with charges opposite to that of the electrode. As surface polarization increases, unique innermost boundary layers composed exclusively of one of the ionic species and the hydrogen bond donor molecules are observed, forming a multilayer structure, with subsequent layers enriched of paired counterions. Interestingly, even at higher applied charges, choline cations and ethylene glycol molecules tended to orient themselves in a parallel fashion toward the electrodes. Differential capacitance curves exhibited a camel-shaped behavior, suggesting a complex interplay of electrochemical processes at the DES|Au(100) interface. In summary, our study provides valuable insights into the interfacial properties of deep eutectic solvents on gold surfaces and their response to changes in temperature and potential, which are crucial for understanding and optimizing deep eutectic solventbased electrochemical systems.

Abstract Regardless of the many common properties shared by ionic liquids (ILs) and deep eutectic solvents (DES), these are two completely different classes of solvents. Unlike ILs, DES are not composed solely of ions, which adds an extra complexity to the systems. Indeed, DES possess a distinct nanostructure organization which is mainly due to the strong hydrogen network present in these solvents. For this reason ILs and DES cannot be treated as similar mixtures, having their own particularities and different promising applications. Nevertheless, both solvents are highly tunable owing to the infinite number of possible combinations of their components and so it is important to tailor their specific properties according to a specific task, by choosing the right components at the right molar ratios. The growing number of applications involving an IL/surface or DES/surface interface is increasing rapidly, but the upscaling of such technologies needs a profound knowledge of chemical compositions, structure and orientational arrangements at these interfaces, also known as the electrical double layer (EDL). In this work, we provide a summary of the progress made in the interfacial area involving ILs and DES, paying special attention to the EDL arrangements and dynamics. It is clear that classical EDL theories do not apply to neither of these solvents, being the one proposed by Kornyshev the most accepted for ILs. Yet, a multilayer organization of the mixtures components at the solvent|electrode interface is recognized by experimental and computational simulation communities for both electrolytes. Each system has its own characteristic differential capacitance curves, which vary greatly due to the different nature of adsorption/desorption species and their reorientation.

Abstract The applicability of deep eutectic solvents is determined by their physicochemical properties. In turn, the properties of eutectic mixtures are the result of the components' molar ratio and chemical composition. Owing to the relatively low viscosities displayed by alcohol-based deep eutectic solvents (DESs), their application in industry is more appealing. Modeling the composition-property relationships established in polyalcohol-based mixtures is crucial for both understanding and predicting their behavior. In this work, a physicochemical property-structure comparison study is made between four choline chloride polyalcohol-based DESs, namely, ethaline, propeline, propaneline, and glyceline. Physicochemical properties obtained from molecular dynamic simulations are compared to experimental data, whenever possible. The simulations cover the temperature range from 298.15 to 348.15 K. The simulated and literature experimental data are generally in good agreement for all the studied DESs. Structural properties, such as radial and spatial distribution functions, coordination numbers, hydrogen bond donor (HBD)-HBD aggregate formation, and hydrogen bonding are analyzed in detail. The higher prevalence of HBD:HBD and HBD:anion hydrogen bonds is likely to be the major reason for the relatively high density and viscosity of glyceline as well as for lower DES self-diffusions.

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.