Fine tuning the polarizable CL&Pol force field for the deep eutectic solvent ethaline
| Autor: | Rafael Maglia de Souza, Mikko Karttunen, Mauro Carlos Costa Ribeiro |
|---|---|
| Rok vydání: | 2021 |
| Předmět: |
Materials science
General Chemical Engineering Deep Eutectic Solvents Ab initio Ionic Liquids Ionic bonding FOS: Physical sciences LÍQUIDOS IÔNICOS Molecular Dynamics Simulation Library and Information Sciences Condensed Matter - Soft Condensed Matter 010402 general chemistry 01 natural sciences Molecular dynamics chemistry.chemical_compound Chlorides Polarizability 0103 physical sciences Physics::Atomic and Molecular Clusters Soft matter Physics::Chemical Physics 010304 chemical physics Force field (physics) General Chemistry 0104 chemical sciences Computer Science Applications Deep eutectic solvent Dipole chemistry Chemical physics Solvents Soft Condensed Matter (cond-mat.soft) |
| Zdroj: | Repositório Institucional da USP (Biblioteca Digital da Produção Intelectual) Universidade de São Paulo (USP) instacron:USP |
| Popis: | Polarizable force fields are gradually becoming a common choice for ionic soft matter, in particular for molecular dynamics (MD) simulations of ionic liquids (ILs) and deep eutectic solvents (DESs). The CL&Pol force field introduced in 2019 is the first general, transferable and polarizable force field for MD simulations of different types of DESs. The original formulation contains, however, some problems that appear in simulations of ethaline and may also have a broader impact. First, the originally proposed atomic diameter parameters are unbalanced, resulting in too weak interactions between the chlorides and the hydroxyl groups of the ethylene glycol molecules. This, in turn, causes an artificial phase separation in long simulations. Second, there is an overpolarization of chlorides due to strong induced dipoles that give rise to the presence of peaks and antipeaks at very low $q$-vector values (\SI{2.4}{\per\nano\meter}) in the partial components of the structure factors. In physical terms, this is equivalent to overestimated spatial nano-scale heterogeneity. To correct these problems, we adjusted the chloride-hydroxyl radial distribution functions against \textit{ab initio} data and then extended the use of the Tang-Toennis damping function for the chlorides' induced dipoles. These adjustments correct the problems without losing the robustness of the CL\&Pol force field. The results were also compared with the non-polarizable version, the CL&P force field. We expect that the corrections will facilitate reliable use of the CL&Pol force field for other types of DESs. |
| Databáze: | OpenAIRE |
| Externí odkaz: |
|