We use atomistic molecular dynamics simulations to study the crystal structure of cellulose nanofibrils, whose sizes are comparable with the crystalline parts in commercial nanocellulose. The simulations show twisting, whose rate of relaxation is strongly temperature dependent. Meanwhile, no significant bending or stretching of nanocellulose is discovered. Considerations of atomic-scale interaction patterns bring about that the twisting arises from hydrogen bonding within and between the chains in a fibril.