Polarization Effects on the Cellulose Dissolution in Ionic Liquids: Molecular Dynamics Simulations with Polarization Model and Integrated Tempering Enhanced Sampling Method

J Phys Chem B. 2017 May 4;121(17):4319-4332. doi: 10.1021/acs.jpcb.6b12647. Epub 2017 Apr 24.

Abstract

Conformation of cellulose with various degree of polymerization of n = 1-12 in ionic liquid 1,3-dimethylimidazolium chloride ([C1mim]Cl) and the intermolecular interaction between them was studied by means of molecular dynamics (MD) simulations with fixed-charge and charge variable polarizable force fields, respectively. The integrated tempering enhanced sampling method was also employed in the simulations in order to improve the sampling efficiency. Cellulose undergoes significant conformational changes from a gaseous right-hand helical twist along the long axis to a flexible conformation in ionic liquid. The intermolecular interactions between cellulose and ionic liquid were studied by both infrared spectrum measurements and theoretical simulations. Designated by their puckering parameters, the pyranose rings of cellulose oligomers are mainly arranged in a chair conformation. With the increase in the degree of polymerization of cellulose, the boat and skew-boat conformations of cellulose appear in the MD simulations, especially in the simulations with polarization model. The number and population of hydrogen bonds between the cellulose and the chloride anions show that chloride anion is prone to form HBs whenever it approaches the hydroxyl groups of cellulose and, thus, each hydroxyl group is fully hydrogen bonded to the chloride anion. MD simulations with polarization model presented more abundant conformations than that with nonpolarization model. The application of the enhanced sampling method further enlarged the conformational spaces that could be visited by facilitating the system escaping from the local minima. It was found that the electrostatics interactions between the cellulose and ionic liquid contribute more to the total interaction energies than the van der Waals interactions. Although the interaction energy between the cellulose and anion is about 2.9 times that between the cellulose and cation, the role of cation is non-negligible. In contrast, the interaction energy between the cellulose and water is too weak to dissolve cellulose in water.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Cellulose / chemistry*
  • Imidazoles / chemistry*
  • Ionic Liquids / chemistry*
  • Molecular Dynamics Simulation*

Substances

  • 1,3-dimethylimidazolium
  • Imidazoles
  • Ionic Liquids
  • Cellulose