Molecular Dynamics Simulation on CL-20/TNT Cocrystal Explosive

Liu Qiang,Xiao Jijun,Zhang Jiang,Zhao Feng,He Zhenghua,Xiao Heming
DOI: https://doi.org/10.7503/cjcu20150605
2016-01-01
Abstract:Molecular dynamics( MD) simulation was conducted for 2,4,6,8,10,12-hexanitrohexaazaiso-wurtzitane (epsilon-CL-20) and 2,4,6-trinitrotoluene( TNT) crystal, CL-20/TNT cocrystal with equal molar ratio and CL- 20/ TNT composite with equal molar ratio using COMPASS force field with isothermal. isobaric( NPT) ensembles at different temperatures. Their structures, interactions between different components and mechani. cal properties were investigated. The cohesive energy density( CED) of the cocrystal decreases with rising temperature, and is greater than that of the composite. The binding energy( E-bind) between two components of the cocrystal also decreases gradually with rising temperature, and is twice times as strong as that of the CL-20/TNT composite. In addition, pair correlation function analysis and enlarged local structure show that hydrogen bonds exist between the two components in the cocrystal and these bonds are mainly formed between H atoms in TNT and O atoms in CL-20, and H atoms in CL-20 and O atoms in TNT. The values of the ten. sile, bulk and shear modulus( E, K, G) for CL-20/TNT cocrystal fall in between those for.. CL-20 and TNT crystals and decrease with increasing temperature, which agrees with general expectation. However, the Cauchy pressure( C12-C44), K/G and Poisson's ratio(v) for CL-20/TNT cocrystal are much larger than those for the two crystals, which predicts the cocrystal has much higher ductility and elastic elongation.
What problem does this paper attempt to address?