2,6-二氨基-3,5-二硝基吡啶-1-氧化物晶体形貌的MD模拟

为探讨溶剂对炸药晶体形貌的影响机制和溶剂的选择依据, 采用附着能(AE)模型预测2, 6-二氨基-3, 5-二硝基吡啶-1-氧化物(ANPyO)在真空的晶体形貌, 确定其主要生长晶面; 并运用分子动力学(MD)模拟研究ANPyO晶面与溶剂N, N-二甲基甲酰胺(DMF)的相互作用及其本质, 进而通过修正的AE模型预测其在DMF中的晶体形貌。结果表明：ANPyO在真空中的晶体形状接近椭圆, 主要生长晶面为(1 1 0)、(1 0 0)、(1 1 -2)及(1 0 -1)面。溶剂与晶面间存在较强相互作用, 径向分布函数分析表明相互作用能主要包括范德华作用, 库仑作用和氢键。在溶剂DMF中, 晶面的修正附着能绝对值顺序为(1 1 0)＜(1 1 -2)＜(1 0 -1)＜(1 0 0), ANPyO晶体形貌接近片状, 与已有实验结果一致。此外, DMF分子在ANPyO晶面的扩散系数研究结果表明, 扩散系数与修正附着能绝对值成线性关系, 晶体生长形貌亦受溶剂扩散能力的影响