- Our Capabilities
- Specialty Builders
- Analysis and Properties
- Mechanical/Thermal Properties
- Phonon - Thermodynamic Properties
- Transition State Search
- P3C Polymer Property Prediction using Correlations
- LAMMPS-Thermal conductivity
- Embedded Atom Potentials
- MedeA Surface Tension
- Compute Engines
- Climbing Length Scales
- Application Notes
- Atomistic Simulations of Multi-Phase Systems
- MedeA ICME seminar
- Users Group Meeting 2016
Anisotropy and temperature dependence of structural, thermodynamic, and elastic properties of crystalline cellulose Iβ: a first-principles investigation.
Modelling and Simulation in Materials Science and Engineering, 2014 vol. 22 (8) p. 085012.
Anisotropy and temperature dependence of structural, thermodynamic and elastic properties of crystalline cellulose Iβ 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 500K) of the monoclinic cellulose Iβ lattice parameters, constant pressure heat capacity, Cp, entropy, S, enthalpy, H, the linear thermal expansion components, ξᵢ, and components of the isentropic and isothermal (single crystal) elastic stiffness matrices, respectively.
Thermodynamic quantities from phonon calculations computed with DFT and the supercell method provided necessary inputs to compute the temperature dependence of cellulose properties via the quasi-harmonic approach. The notable exceptions were the thermal conductivity components, λᵢ (the prediction of which has proven to be problematic for insulators using DFT) for which the reverse, non-equilibrium molecular dynamics approach with a forcefield 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β structure and bonding interactions.