CO在高温高压下的反应反应动力学的研究.pdfVIP

第八届全国冲击动力学学术讨论会 Co在高温高压下的反应反应动力学研究 姬广富‘李晓风1·3张艳丽1，2赵峰1孟川民1宋振飞1刘永刚1，3 (1中国工程物理研究院流体物理研究所冲击波物理与爆轰物理国防重点实验室，621900 2西南交通大学高温高压物理研究所．61003l 3四川大学原子与分子物理研究所．四川成都，610065) 摘要：用理论模拟方法研究炸药降解产物的动力学过程，可以为超高速燃烧和高速空气动力学等研究提供参考，也 可以对炸药降解产物在高温高压下动态化学变化过程提供更直接的理解。本文用从头算分子动力学方法结合平面波赝势 应过程。结果发现，在此压力和温度条件下c0发生了歧化分解形成了更复杂的混合物：C0、c02、c：O和以及c】02， 没有形成在更高压力下才能产生的cO聚合物。这些结果为理解炸药爆轰分解过程和极端条件下原子分子结构变化以及 粒子间相互作用规律有着重要的意义。 关键词： cO炸药分解产物从头算分子动力学 对炸药降解产物的状态方程进行研究，可以为超高速燃烧以及高速空气动力学等研究提供数据和理 论方法。通常，炸药在爆轰过程中，cJ状态时产物气体的密度可以高达2幽m3以上，温度可以到3000K 以上。而在聚心爆轰或核爆加载期间，其密度和温度将会更高。为了全面了解和认识它们在此种状态时 的行为特性，科学家们借助二级轻气炮、激光等技术手段，分别对爆轰产物各成份气体进行冲击波加载， 立相应的物态方程。在分子水平上理解小分子物质在高温高压下的物理与化学性质可以帮助我们了解炸 药爆炸过程以及炸药在爆轰条件下产生的小分子爆轰产物的性质具有非常重要的意义。 一氧化碳(cO)是碳原子的二价化合物，在常温常压下可以稳定存在又具有一定的还原性，它是凝 聚态炸药的主要爆轰产物。因此．对co高温高压状态方程及辐射、光谱特性研究无论是在科学或是技术 80GPa下仍是 方面都有很重要的意义。这主要归因于CO不像N2属于等电子配对(N2的分子晶体至少在1 稳定的[R cO在5GPa及室温下就已经变得不稳定了【A．I_Katz等，1984，以及R．LMills等，1986】。通常认为此时会 za_n等，1990)、(F Cansell等，1993)、(M．NicoI等，1986)以及(H．G 击波实验和静高压实验的验证。这些过程的微观机理也有人做过研究，例如有人用从头算方法研究了乙 压下c0的化学反应，他们分析了产物的可见光谱、拉曼光谱和傅立叶变换红外光谱等。通过分析吸收 光谱和拉曼光谱确定了该物质包含c02、类似石墨结构的碳以及一种可确认的重复单元的聚合物。但对 CO气体而言是否也会形成聚合物．至今少见报道。通过对其高压物态方程及辐射、光谱特性的研究， 我们可以深入地了解该物质在极端条件下的分子间及分子内的物理、化学过程．特别是在高温高压下cO 会分解形成复杂的混合反应物．并且可能有碳的析出。能否定性甚至定量地描述其化学组分，对校验现 有理论和计算结果的正确性、准确性以及对爆轰产物化学反应过程的理解都是非常重要的。况且，通过 对物质在极端条件下的宏现状态测量和计算到对其辐射、光谱性质等这些能反应物质原子、分子性质的 物理量的实时测量，可以深刻理解组成混合物的原子分子结构变化以及粒子间相互作用规律。这是高压 凝聚态物理学延伸到极端条件和微观层次的必然趋势。 本文用从头算分子动力学方法结合平面波赝势基组以及广义梯度近似(PP．GGA)模型模拟研究了 cO在高温高压下的反应反应动力学研究 爆轰分解过程和极端条件下原子分子结构变化以及粒子间相互作用规律有着重要的意义。 计算方法 计算中，我们采用广义梯度近似方法[(A．D．Becke，1988)；(J．PPcrdew1986)】。C和O均采用 p4’c和。的1s 态包含在赝核内，价电子波函数以及缀加电子密度在平面波基组上展开以截断能表示时分别取470eV和 拟中，波函数内插方案采用变系数一阶外推方法。模拟温度1000K，压力050Pa。模拟时间为2ps