Solvent behavior of an ionic liquid set around a cellulose Iβ crystallite model through molecular dynamics simulations

A set of imidazolium-based ionic liquids: [C4mim][PF6], [C4mim][BF4], [C4mim][Cl], [C4mim][CF3COO], [C4mim][NTf2], [C4mim][OMs], [C4mim][Br], and [C4mim][OAc], was studied by molecular dynamics simulations to elucidate their solvent behavior around a crystallite model of cellulose Iβ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}, through atomistic interactions and the degree of departure of its thermodynamic properties from their solvent pure phase. These departure changes were correlated with experimental values of the Kamlet-Taft solvent basicity parameter, and it was found that, even at room temperature, density changes, and vaporization enthalpy changes can be correlated with the capacity of ionic liquids for the preconditioning of the cellulose crystallite. Hydrogen bond occupancies indicate that ionic liquids can disrupt external chains of the cellulose crystallite by replacing and reducing the strong O6-H⋯O2/O3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$O6 - H \cdots O2/O3$$\end{document} hydrogen bonds by weak hydrogen bonds such as O6-H⋯O4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$O6 - H \cdots O4$$\end{document} along the interchain network. Also, radial distribution functions indicated that structural changes in the cellulose-ionic liquid mixtures did not depart significantly with respect to the pure IL structure. The results of the free energy of solvation calculations for a cellulose chain, presented the following trend: [C4mim][Cl] > [C4mim][OAc] > [C4mim][CF3COO] > [C4mim][Br] > [C4mim][OMs] > [C4mim][BF4] > [C4mim][PF6] > water > [C4mim][NTf2]. It is important to emphasize, that the focus of this work was not the cellulose dissolution, but instead, the solvent behavior and cellulose preconditioning within each IL at room temperature. Our results can provide insights about the preconditioning stage of cellulose at low temperature, useful in the development of lignocellulosic materials and valuable cellulose derivatives by means of low energy requirements.