The influence of temperature and component proportion on stability, sensitivity, and mechanical properties of LLM-105/HMX co-crystals via molecular dynamics simulation

Based on molecular dynamics (MD) simulation, the binding energy, cohesive energy density (CED), bond length, and mechanical parameters were calculated for 2,6-diamino-3,5-dinitropyrazine-l-oxide (LLM-105) crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) crystal, and their co-crystals under different temperatures. Three LLM-105/HMX patterns were constructed to investigate the influence of component proportion on structures and properties of co-crystals, in which the mole ratios of LLM-105 and HMX are 1:1, 1:2, and 2:1. The effect of temperature and components on the stability and sensitivity were investigated as well. The results show that the binding energies, CED and mechanical parameters of all the co-crystals, decrease when the temperature increases from 248 to 398 K, while their maximum N-NO2 bond length (L-max) increases with rising temperature, indicating that the sensitivities increase and stabilities decrease when temperature rises. At all temperatures, co-crystals exhibit larger CED and shorter bond length than that of single explosive, demonstrating that they are more stable and less sensitive than single crystal, where the stability of co-crystals was ordered as 2:1>1:1>1:2. Moreover, the bulk modulus (K) and shear modulus (G) of co-crystals are lower than that of HMX, conversely, the Cauchy pressure and K/G are higher than that of HMX, implying co-crystals have better ductility. Finally, the 2:1 ratio of LLM-105/HMX co-crystal was identified as the excellent one, owning to the highest binding energy, highest CED, shortest L-max, and greatest ductility. Models of LLM-105/HMX and one of the properties