液相硝基甲烷振动能量驰豫的二体碰撞模型研究.pdfVIP

赵波等：液相硝基甲烷振动能量驰豫的二体碰掩模世研究 液相硝基甲烷振动能量驰豫的 二体碰撞模型研究1 赵波，崔季平，樊箐 (中国科学院力学研究所高温气体动力学重点实验室，北京100190) C摘要】t本文从能量驰豫的二体碰撞模型、晶体格波振动理论及Dlo仕等人提出的。门槛”模观点研究液相硝 基甲烷振动能量驰豫过程。二体碰掩模型中的单次碰掩能量传输几率用分子动力学打靶法确定，碰撞频率用分子动 力学方法确定并与晶格振动Debye频率关联，以此给出了振动激发速率常数及振动驰豫时间。本文研究将气相及固 相的理论成果向液相进行了推广，并指出高速冲击下化学反应与振动激发产生耦合现象，此外应当考虑集体相互作 用现象。 【关键词】：振动能量驰豫，硝基甲烷，二体碰撞模型，分子动力学 0引言 在含能材料冲击起爆与持续爆轰的传播过程中，冲击波激发声子浴，声子浴能量向分子内部传 递，是含能分子开始起始反应的重要前提。目前研究表明，这一过程的特征时间是标志炸药感度的 一个重要特征量【¨。这一过程的控制步是被称为“门槛振动模”的激发速犁2’31，一旦“门槛振动模” 的激发形成之后，分子内部不同模之间的平衡会在很短的时间内完成14l。这一阶段(ps级)完成的 后果是含能材料的分子内态被高度激发，导致分子失稳发生化学反应峨酬。因此研究“门槛模”振动 激发过程极其重要，它的驰豫过程决定了整个分子振动能量驰豫过程，确定“门槛模”的能量传输 几率既可以确定整个分子振动能量传输几率。 在凝聚态介质中，振动能量传递过程即分子振动能驰豫过程的理论共有两类。一类是二体碰撞 模型，一类是声子浴模型。二体碰撞模型是由气体分子振动驰豫理论引伸出来，探讨向凝聚相的推 White．Millikaa关联，这样的理论和大量的实验，来推定碰撞传能几率。液态二体模型中如何界定“碰 撞”和“碰撞频率”则转用类固模型。例如用液态笼子结构的Debye模型，把Debye频率作为“碰 撞频率”。 声子浴模型是基于晶格动力学中格波的非线性相互作用Il¨，振动传能过程涉及到晶体格波声学 支与光学支之间的能量传输，完全移用于液体也有困难，也需要简化处理【2J。对于简单分子液体， 二体碰撞模型已经成功解释有关了液态含氢双原子分子振动能量驰豫实验现象【l}№】。本文中我们引 入分子动力学方法将二体碰撞模型推广到凝聚态多原子分子的振动能量驰豫研究中。硝基甲烷作为 一种典型的多原子含能材料，有比较多的实验数据可参考，本文也选它为对象，定量研究其“门槛 模”振动能量驰豫过程。 1 理论 1 究方向：凝聚态炸药爆轰微观动力学研究．中国科学院力学研究所高温气体动力学重点实验室，北京100190。 馋十四届中园卒气动力学物理气体动力学学术交流会论文集 詈置—鼍詈暑曼詈鼍詈詈置曼量I ． 一一 一 ． 一． I鼍鼍曼詈詈曼詈詈鼍詈曼詈詈皇曼 1．1二体碰撞模型 分子振动能量驰豫过程是用振动能量驰豫速率常数来定量描述。用二体碰撞理论计算能量驰豫 速率时，它可以表示为两个因子之积：(a)单位时间内碰撞频率，(b)单次碰撞能量传输儿率。对 于气体而言，这两个因子都很容易确定。将二体碰撞理论应用到液体中时，作为一个近似模型一般 方法是单次碰撞能量传输几率采用气相实验结果或者是散射理论计算，用比拟量界定液体中的碰撞 频率。对于液相多原子分子用气相实验结果计算单次碰撞能量传输几率已经不合适【17l，严格理论计 算也很复杂。因此在研究液相硝基甲烷时，我们借用Dlott等人提出的“门槛模”理论把问题简化为 “门槛模”的振动传输几率，转而利用分子动力学的计算机模拟研究。对于碰撞频率，我们用分子 动力学方法直接统计模拟，结合固体Debay理论进行理论计算。 1．2模拟条件 在研究单次碰撞能量驰豫几率时采用的是打靶法，因此系统只需要由两个分子组成即可，分子 内化学键作用能及分子原子间作用势在下节中进行说明。用分子动力学方法直接模拟液相硝基甲烷 分子的碰撞频率时，系统为1000个分子组成立方体系，采用周期性边界条件，分子间作用势为 Lelmard-Jones(L-J)势 删占吟)哆俐] ㈩ 其中D，g为势参数，，．为距离，越为分子问作用势。由等温压缩