Anisotropy and temperature dependence of structural, thermodynamic, and elastic properties of crystalline cellulose I sub( beta ): a first-principles investigation

• Published in: Modelling and simulation in materials science and engineering 2014-12, Vol.22 (8)
• Main Authors: Dri, Fernando L, Shang, Shunli, Hector, Louis G, Saxe, Paul, Liu, Zi-Kui, Moon, Robert J, Zavattieri, Pablo D
• Format: Article
• Language: English
• Subjects: Anisotropy; Bonding; Cellulose; Computation; Constants; Elastic constants; Temperature dependence; Thermodynamics
• ISSN: 0965-0393; 1361-651X
• DOI: 10.1088/0965-0393/22/8/085012

Anisotropy and temperature dependence of structural, thermodynamic and elastic properties of crystalline cellulose I sub( beta ) were computed with first-principles density functional theory (DFT) and a semi-empirical correction for van der Waals interactions. Specifically, we report the computed temperature variation (up to 500?K) of the monoclinic cellulose I sub( plus or minus ) lattice parameters, constant pressure heat capacity, C sub(p), entropy, S, enthalpy, H, the linear thermal expansion components, xi sub(i), and components of the isentropic and isothermal (single crystal) elastic stiffness matrices, (ProQuest: Formulae and/or non-USASCII text omitted)(T) and (ProQuest: Formulae and/or non-USASCII text omitted)(T), respectively. Thermodynamic quantities from phonon calculations computed with DFT and the supercell method provided necessary inputs to compute the temperature dependence of cellulose I sub( beta ) properties via the quasi-harmonic approach. The notable exceptions were the thermal conductivity components, lambda sub(i) (the prediction of which has proven to be problematic for insulators using DFT) for which the reverse, non-equilibrium molecular dynamics approach with a force field was applied. The extent to which anisotropy of Young's modulus and Poisson's ratio is temperature-dependent was explored in terms of the variations of each with respect to crystallographic directions and preferred planes containing specific bonding characteristics (as revealed quantitatively from phonon force constants for each atomic pair, and qualitatively from charge density difference contours). Comparisons of the predicted quantities with available experimental data revealed reasonable agreement up to 500?K. Computed properties were interpreted in terms of the cellulose I sub( beta ) structure and bonding interactions.