论文部分内容阅读
含能材料(EMs),包括炸药、推进剂和烟火剂等,已广泛应用于国民经济和国防工业等各个领域。因此,含能材料的基础研究对国民经济的增长和国防工业的发展都有着重要的意义。由于含能材料的实验研究具有一定的危险性,计算机模拟逐渐成为含能材料研究的重要手段之一。本文基于密度泛函理论,并结合分子动力学模拟,分别研究了分子晶体中分子间相互作用,金属有机物骨架聚合物含能材料的晶体结构、电子结构和弹性性质,以及新型钝感炸药FOX-7晶体在常压和高压下的热解行为。目前,大部分含能材料属于分子晶体,分子晶体中分子间的相互作用决定了其晶体结构及其物理化学性质描述的准确性。然而,基于局域密度近似或者广义梯度近似的密度泛函理论由于不能描述长程关联相互作用,不能很好地描述分子晶体中的相互作用。为了能更好地描述分子间相互作用,人们已经发展了多种方案来弥补密度泛函理论框架下这些近似的缺陷。虽然人们对分子二聚体、分子团簇等体系进行了大量的测试,但是这些修正方案在描述分子晶体的表现却知之甚少。本文系统地测试了密度泛函理论及其修正方案在描述分子晶体中相互作用时的表现,为分子晶体的物理化学性质及其动力学性质的研究提供可靠的方法指导。我们首先对甲烷晶体的结构进行了测试,发现通过加入经验对势来修正DFT的PBE-Grimme(PBE-D2)和PBE-D3方案可以很好地重复实验的晶格常数和晶格能,而且和标准DFT计算相比计算成本增加不明显;而vdW-DF方案尽管可以给出合理的晶格常数和体弹性模量,但却高估了甲烷晶体的晶格能,并且增加了计算成本。更进一步,我们采用合作者针对包含C、H、N和O体系发展了DFT+LAP方案,计算了十二种分子晶体的晶格常数和晶格能。通过和实验值及PBE-Grimme的计算结果相对比,评估了此方案在描述分子间相互作用时的表现。结果表明,我们合作者发展的DFT+LAP得到的计算结果与实验值相符合,而且计算量也与DFT方案基本相当;另一方面,PBE-Grimme方案是测试的所有方案中最好的,我们的计算结果与之相近。因此,DFT+LAP方案在研究分子晶体和其它非共价键相互作用体系时是一种高效且精确的计算方法。随后,我们更进一步改进和发展了DFT+LAP方法,研究指出需要更有效的势函数才可以更精确描述分子晶体中分子间的相互作用。和传统炸药相比,人们己经合成出了高含氮量且具有更优异的爆轰性能的金属有机骨架聚合物(MOF),这类材料被证实为潜在的新型含能材料。金属有机骨架聚合物的晶体结构极其复杂,其电子结构和力学性质对这类潜在含能材料的应用有着重要的影响。本文采用第一性原理方法系统地研究了一维、二维和三维MOF-EMs的晶体结构、电子结构和弹性性质。我们首先对MOF-EMs的晶体结构用PBE泛函及色散修正方案进行测试。结果表明,PBE泛函明显地低估了MOF-EMs晶体中的非共价键相互作用,从而高估了晶格常数。所有的修正方法都或多或少地修正了对分子晶体晶格常数的描述,但是一些方案如PBE-D和optB88-vdW却过度修正这种效应。在所有方案中,optPBE-vdW可以获得与实验最一致的晶体结构。电子结构分析表明只有NHN为金属,其余MOF-EMs为半导体或绝缘体。我们进一步讨论了MOF-EMs的电子结构和其感度之间的关系,指出相似结构的带隙越大,MOF的冲击感度越高,则MOF越不敏感,这和实验相符。此外,我们预测的MOF-EMs的体弹性模量介于15.1-35.0 GPa之间,高于传统炸药的体弹性模量。除了含能材料的基本性质,其爆轰机制的研究也是含能材料的研究热点。FOX-7作为一种新型钝感炸药,理解其热解机制不仅对FOX-7的储藏、运输和应用都有着重要的意义,还有助于新型钝感炸药的设计。本文采用第一性原理分子动力学模拟研究了FOX-7晶体在常压和高压下的热解过程。常压热解过程表明FOX-7的起始分解从C—NO2键的断裂开始。我们详细地分析了热解产物及其演化规律。结果表明热解产物包含大量的N2、H2O,还包含少量的CO、CO2等。同时发现中间产物NO、NO2、OH等分子片段都携带电荷,这些带电分子片段可能会促进FOX-7的后续化学反应。进一步模拟压力对热解的影响指出压力的增加并不直接对热解机制产生影响,却加快了热解反应发生的速率。