NpH-MD-Simulations of the Elastic Moduli of Cellulose II at Room Temperatue

We have used molecular dynamics modeling to investigate the stucture and mechanical properties of regenerated cellulose fibres. This work is motivated by continued interest in replacing the environmentally hazardous viscose process by alternative spinning methods. An important input parameter for any realistic model of the elastic properties is the stiffness tensor of the crystalline constituent, cellulose II. Conventional molecular mechanics techniques can be used to estimate the elastic reaction of a material with respect to small external stresses or strains, i.e. the compliance and stiffness tensors, and the elastic moduli derived therefrom, at zero temperature. In order to access non-zero temperatures, it is necessary to use either the quasi-harmonic approximation for the vibrational free energy or molecular dynamics (MD) simulations. In the present work, Parrinello-Rahman constant-stress MD was performed to generate trajectories in constant particle number (N), constant external stress tensor (p or t) and constant enthalpy H (NpH or HtN) ensemble. This was found to be less time consuming than working with isothermal conditions, as done by other authors. The fluctuations in kinetic energy and MD cell vectors were then used to calculate adiabatic elastic constants, thermal expansion coefficients and heat capacity. The isothermal elastic constants were found by applying a standard thermodynamic relation. The Youngs modulus along the chain direction, E_l, was determined to be 155 GPa, whereas the values in the perpendicular directions vary between 51 and 24 GPa. These results are of the same order of magnitude as those obtained by Tashiro and Kobayashi 1] with the static (T = 0K) method, but our value of E_l is 5% lower and, unexpectedly, the lateral values are up to six times higher. A strong anisotropy is found for shear along the chains in planes containing the chain axis, the shear modulus ranging from 5 to 20 GPa. Convergence was achieved in the simulations, to the extend that the elastic constants become stationary, but significant internal stresses remain, pointing to shortcomings in the software used. Further work is necessary to resolve these problems, although the major conclusions should be unaffected.