Molecular dynamics simulations on epsilon-CL-20-based PBXs with added GAP and its derivative polymers

Molecular dynamics simulations have been employed to study the epsilon-CL-20-based PBXs under COMPASS force field. epsilon-CL-20 was chosen as the base explosive due to its higher energy, density and detonation performance than conventional explosives. Four polymers, GAP, GAP-NH2, GAP-NO2 and GAP-NH2-NO2 were added into the epsilon-CL-20(001) crystalline surface to build the PBX models. The cohesive energy densities (CEDs), elastic coefficients, isotropic mechanical properties (Young's moduli, bulk moduli, shear moduli, Poisson's ratio, Cauchy pressure and K/G) and initiation bond length distribution were studied. It turned out that the CEDs order was A1 < A4 < A3 < A2 < A. The mechanical properties of pure epsilon-CL-20(001) were effectively improved by building PBX models. System A3 showed better comprehensive mechanical properties than the other three PBXs. A study on the initiation bond length distribution showed that the L-max and L-ave of N-NO2 increased with increasing temperature and they were related to the sensitivity of the explosives. The order of Lmax was A3 < A4 < A2 < A1 < A, which indicated that the PBXs owned lower sensitivity than system A. These studies are thought to provide guidance for further research on the application of GAP and its derivative polymers. Meanwhile, they are meaningful for the studies on epsilon-CL-20-based PBXs.