Molecular dynamics simulation of polymorphic transformation between ε-CL-20 and β-CL-20 by short range order and orientational order parameter

Phase transition simulation of molecular crystals is a challenge. In order to clarify the polymorphic transformation mechanism of molecular crystal, a molecular dynamics (MD) simulation of the transformation between ε-CL-20 and β-CL-20 was carried out by the short range order g(r) and local orientational order parameter (OPs) P2(cosθ) as well as the spatial orientation correlation function. The effects of the temperature on them were performed. The results indicate that the g(r), orientational OPs and spatial orientation correlation function, composed of the structural parameters from three CL-20 molecules in three different single crystal, could be used to simulate the process of the polymorphic transformation between ε-CL-20 and β-CL-20. Only at a high temperature, could two kinds of the CL-20 supercells realize the transformation of crystal form. Furthermore, there were local and global minimum energy points, indicating that there could be the metastable or stable crystal forms in the process of the transformation from β-CL-20 to ε-CL-20. The values of P2(cosθ) were well related to results of the surface electrostatic potentials (ESPs), suggesting that the ESPs have an influence on polymorphic transformation.