Partial crystalline transformation of solvated cellulose IIII crystals, reproduced by theoretical calculations

In hot-water molecular dynamics simulation at 370 K, four cellulose III_I crystal models, with different lattice planes and dimensions, exhibited partial crystalline transformations of (1 ?1 0) chain sheets, in which hydroxymethyl groups were irreversibly rotated from gt into tg conformations, accompanied by hydrogen-bond exchange from the original O3–O6 to cellulose-I-like O2–O6 bonds. The final hydrogen-bond exchange ratio was about 95 % for some of the crystal models after 50 ns simulation. The corrugated (1 ?1 0) chain sheet was converted to a cellulose-I-like flat chain sheet with a slightly right-handed twist. The 3D structures of the three types of isolated chain sheet models were optimized using density functional theory calculations to compare their stabilities without crystal packing forces. The cellulose Iβ (1 0 0) models were more stable than the cellulose III_I (1 ?1 0) models. The optimized structure of cellulose III_I (1 0 0) models deviated largely from the initial sheet form. It was proposed to the crystalline transformation from cellulose III_I to Iβ that conversion of the chain sheet structure first take place, followed by sliding of the chain sheet along the fiber axis.