Cellulose Iβ Behaviors in Non-solvent Liquid Media: Molecular Dynamic Simulations

The structural changes of cellulose in non-solvent liquid media can provide insights into the high-value utilization of cellulose. This study includes molecular dynamics simulations of 36-chain cellulose Iβ microfibril model (Iβ-MF) behavior in 16 non-solvent liquids with different polarities at room temperature using two carbohydrate force fields (CHARMM36, GLYCAM06). Iβ-MF in CHARMM36 retains more than 70% of the tg conformation in 16 liquids, and the retention of the tg conformation increased with decreasing liquid polarity. Liquid polarity can affect the hydroxymethyl conformation of cellulose, which is only an appearance, and the real driving force behind is the electrostatic interaction between liquid molecules and cellulose. Furthermore, changing the 1,4 electrostatic scaling factor of GLYCAM06 can effectively affect the structural convergence of Iβ-MF. The Iβ-MF forms an alternating layer structure in the gg/gt conformation in a medium to high polarity non-solvent liquid, while the model undergoes untwisting. Model untwisting is inextricably linked to the degree of alternate layer structure formation. This paper provides a theoretical basis for the molecular study of nanocellulose structures from an energy-structure-property perspective.