Ab initio studies on the structure of and atomic interactions in cellulose IIII crystals

• Published in: Carbohydrate research 2015-11, Vol.417, p.72-77
• Main Authors: Ishikawa, Tetsuya, Hayakawa, Daichi, Miyamoto, Hitomi, Ozawa, Motoyasu, Ozawa, Tomonaga, Ueda, Kazuyoshi
• Format: Article
• Language: English
• Subjects: CH/O interaction; Crystal structure of cellulose IIII; Density functional theory; FMO4; NBO
• ISSN: 0008-6215; 1873-426X
• DOI: 10.1016/j.carres.2015.09.006

•The crystal structure of cellulose IIII was analyzed using first-principles DFT.•The calculation well reproduced the experimental structure obtained by Wada et al.•NBO analysis suggested the bonding nature of CH/O interaction.•CH/O and hydrogen bonding on the stabilization of crystal structure are comparable. [Display omitted] The crystal structure of cellulose IIII was analyzed using first-principles density functional theory (DFT). The geometry was optimized using variable-cell relaxation, as implemented in Quantum ESPRESSO. The Perdew–Burke–Ernzerhof (PBE) functional with a correction term for long-range van der Waals interactions (PBE-D) reproduced the experimental structure well. By using the optimized crystal structure, the interactions existed among the cellulose chains in the crystal were precisely investigated using the NBO analysis. The results showed that the weak bonding nature of CH/O and the hydrogen bonding occur among glucose molecules in the optimized crystal structure. To investigate the strength of interaction, dimeric and trimeric glucose units were extracted from the crystal, and analyzed using MP2 ab initio counterpoise methods with BSSE correction. The results estimated the strength of the interactions. That is, the packed chains along with a-axis interacts with weak bonding nature of CH/O and dispersion interactions by −7.50 kcal/mol, and two hydrogen bonds of O2HO2…O6 and O6HO6…O2 connect the neighboring packed chains with −11.9 kcal/mol. Moreover, FMO4 calculation was also applied to the optimized crystal structure to estimate the strength of the interactions. These methods can well estimate the interactions existed in the crystal structure of cellulose IIII.