Thermal Response in Cellulose Iβ Based on Molecular Dynamics

Abstract The structural details of cellulose I β were discussed according to molecular dynamics simulations with the GLYCAM-06 force field. The simulation outcomes were in agreement with previous experimental data, including structural parameters and hydrogen bond pattern at 298 K. We found a new conformation of cellulose Iβ existed at the intermediate temperature that is between the low and high temperatures. Partial chain rotations along the backbone direction were found and conformations of hydroxymethyl groups that alternated from tg to either gt or gg were observed when the temperature increased from 298 K to 400 K. In addition, the gg conformation is preferred than gt. For the structure adopted at high temperature of 500 K, major chains were twisted and two chains detached from each plain. In contrast to the observation under intermediate temperature, the population of hydroxymethyl groups in gt exceeded that in gg conformation at high temperature. In addition, three patterns of hydrogen bonding were identified at low, intermediate and high temperatures in the simulations. The provided structural information indicated the transitions occurred around 350 K and 450 K, considered as the transitional temperatures of cellulose Iβ in this work.