Modeling crystal and molecular deformation in regenerated cellulose fibers

Experimental deformation micromechanics of regenerated cellulose fibers using Raman spectroscopy have been widely reported. Here we report on computer modeling simulations of Raman band shifts in modes close to the experimentally observed 1095 cm(-1) band, which has previously been shown to shift toward a lower wavenumber upon application of external fiber deformation. A molecular mechanics approach is employed using a previously published model structure of cellulose II. Changing the equilibrium c-spacing of this structure and then performing a minimization routine mimics tensile deformation. Normal-mode analysis is then performed on the minimized structure to predict the Raman-intensive vibrations. By using a dot-product analysis on the predicted eigenvectors it is shown that some Raman active modes close to the 1095 cm(-1) band interchange at certain strain levels. Nevertheless, when this is taken into account it is shown that it is possible to find reasonable agreement between theory and experiment. The effect of the experimentally observed broadening of the Raman bands is discussed in terms of crystalline and amorphous regions of cellulose, and this is compared to the lack of X-ray broadening to explain why discrepancies between theory and experiment are present. A hybrid model structure with a series-parallel arrangement of amorphous and misaligned amorphous-crystalline domains is proposed which is shown to agree with what is observed experimentally. Finally, the theoretical crystal modulus for cellulose II is reported as 98 GPa, which is shown to be in agreement with other studies and with an experimental measurement using synchrotron X-ray diffraction.