Simulation of melting in crystalline polyethylene

We carry out a molecular dynamics simulation of the first stages of constrained melting in crystalline polyethylene (PE). When heated, the crystal undergoes two structural phase transitions: from the orthorhombic (O) phase to the monoclinic (M) phase, and then to the columnar (C), quasi-hexagonal, phase. The M phase represents the tendency to the parallel packing of planes of PE zigzags, and the C phase proves to be some kind of oriented melt. We follow both the transitions O→M and M→C in real time and establish that, at their beginning, the crystal tries (and fails) to pass into the partially ordered phases similar to the RI and RII phases of linear alkanes, correspondingly. We discuss the molecular mechanisms and driving forces of the observed transitions, as well as the reasons why the M and C phases in PE crystals substitute for the rotator phases in linear alkanes.