Zhang Yuqing, Wang Jingtao
School of Chemical Engineering and Technology, Tianjin University, Tianjin, 300350, People's Republic of China.
J Mol Model. 2025 May 30;31(6):178. doi: 10.1007/s00894-025-06396-7.
Cocrystals are recognized as an effective strategy to mitigate the high sensitivity of energetic materials. In this paper, we use molecular dynamics (MD) to simulate the binding energies and XRD of CL-20/DNAN cocrystals with different ratios to prove that there is a new type of cocrystal structure between CL-20 and DNAN. To further investigate the effect of moisture absorption on the storage safety of cocrystals, this paper simulates the adsorption process of water molecules on the (0 0 1), (0 1 0), (1 0 0) crystal surfaces, respectively. By analyzing the radial distribution function plots, it is found that there are hydrogen bonding interactions between water molecules and molecules on the surface of the cocrystal, and the (1 0 0) face is most likely to adsorb water molecules. By analyzing the length of the triggering bond and the cohesive energy density of the cocrystal, it was found that when water molecules are absorbed on the surface, the N-NO2 bond of the cocrystal explosive is more likely to be broken and the sensitivity of the explosive is higher on the surface, leading to increased sensitivity.
The MD simulation is conducted utilizing the Materials Studio software, operating under the NPT set with a temperature of 298 K, a pressure of 0.0001 GPa, temperature control of Andersen, and pressure control of Berendsen. The simulation spans a duration of 150 ps, with samples being recorded at an interval of 1 fs throughout the computational process. Similar simulations are conducted for all systems.
共晶体被认为是缓解含能材料高敏感性的一种有效策略。在本文中,我们使用分子动力学(MD)模拟不同比例的CL-20/DNAN共晶体的结合能和XRD,以证明CL-20和DNAN之间存在一种新型共晶体结构。为进一步研究吸湿对共晶体储存安全性的影响,本文分别模拟了水分子在(0 0 1)、(0 1 0)、(1 0 0)晶体表面的吸附过程。通过分析径向分布函数图,发现水分子与共晶体表面分子之间存在氢键相互作用,且(1 0 0)面最容易吸附水分子。通过分析引发键的长度和共晶体的内聚能密度,发现当水分子吸附在表面时,共晶体炸药的N-NO2键更容易断裂,炸药在表面的敏感性更高,导致敏感度增加。
使用Materials Studio软件进行MD模拟,在NPT系综下运行,温度为298 K,压力为0.0001 GPa,采用Andersen温度控制和Berendsen压力控制。模拟时长为150 ps,在整个计算过程中以1 fs的间隔记录样本。对所有系统进行类似模拟。