Cellulose microfibrils are recalcitrant toward dissolution, thus it is difficult to extract and characterize them without modifying their native state. To study the molecular level behavior of microfibrils over 100 sugar residues, we construct a coarse-grained model of solvated cellulose Iβ microfibril using one bead per sugar residue. We derive the coarse-grained force field from atomistic simulation of a 36 chain, 40-residue microfibril by requiring consistency between the chain configuration, intermolecular packing and hydrogen bonding of the two levels of modeling. Coarse-grained force sites are placed at the geometric center of each glucose ring. Intermolecular van der Waals and hydrogen bonding interactions are added sequentially until the microfibril crystal structure in the atomistic simulation is achieved. This requires hydrogen bond potentials for pairs that hydrogen bond in cellulose Iβ, as well as those that can hydrogen bond in other structures, but not in cellulose Iβ. Microfibrils longer than 100 nm form kinks along their longitudinal direction, with an average periodicity of 70 nm. The behavior of kinked regions is similar with a bending angle of approximately 20°. These kinked regions might be linked to observations of periodic disorder from small angle neutron scattering and acid hydrolysis.
All Science Journal Classification (ASJC) codes
- Polymers and Plastics