α-quartz型GeO2高压相变的分子动力学研究

采用在经典离子晶体作用势中附加Ｍｏｒｓｅ势,并进行必要的量子化修正,对α－ｑｕａｒｔｚ型ＧｅＯ２结构随压力变化特性,进行分子动力学计算模拟,获得了在压力高于６．０ＧＰａ,α－型ＧｅＯ２从晶相向非晶相相变的模拟结果,并利用其摩尔体积变化,键角、径向分布１、配位数等重要信息对模拟结果作了深入的探讨,相变后的非晶相,由占体积６６％的八面体结构和３３％的四面体结构组成的非晶体,其中还有极少量的α－ｑｕａｒｔｚ型ＧｅＯ２存在。     

N2双原子系统经典轨迹的辛算法计算

采用辛算法和Ｒ－Ｋ法计算了Ｎ２双原子系统的经典轨迹，结果表明，辛算法保持系统的辛结构和固有性质不变，与实验和理论一致，Ｒ－Ｋ法则与实验和理论不符。     

PuO2体系的分子反应动力学研究

基于ＰｕＯ２分子（Ｘ＾５Σｇ＾＋）的多体展式分析势能函数,用准经典的Ｍｏｎｔｅ－Ｃａｒｌｏ轨迹法研究了Ｏ（＾３Ｐｇ）＋ＰｕＯ（０,０）和Ｐｕ（＾７Ｆｇ）＋Ｏ２（０,０）的分子反应动力学过程。结果表明：Ｏ（＾３Ｐｇ）＋ＰｕＯ（０,０）碰撞生成少量稳定络合物ＰｕＯ２（Ｘ＾５Σｇ＾＋）,该反应是无阈能反应;而Ｐｕ（＾７Ｆｇ）＋Ｏ２（０,０）碰撞生成稳定络合物的同时,主要还发生了Ｐｕ（＾７Ｆｇ）＋Ｏ２（０     

高温下CO2—N2,CO2—CO2混合气体中CO2分子的谱线压力展宽半宽度

用无限阶突变近似方法计算了ＣＯ２－Ｎ２,ＣＯ２－ＣＯ２混合气体中Ｃ承２９６－１９２９Ｋ温度下转动量子数从１到３０的谱线压力展宽半宽度。计算过程中各系统的相互作用势采用改进的Ｂｕｃｋｉｎｇｈａｍ模型势。得到的结果和现有结果符合较好。     

铍团簇Bem（n—2—55）的基态能量,结构及其拓扑性质

采用一种第一性原理多体势来描述团簇Ｂｅｎ中原子之间的相互作用,用模拟自然界生物进化方法系统地优化计算了Ｂｅｎ团簇（ｎ＝２－５５）的基态能量和几何结构,分析了其拓扑性质,总结了能量和结构随原子数目的变化规律,并将得到Ｂｅｎ的结果和长程Ｍｏｒｓｅ势所获这的金属团簇Ｍｎ,以及Ｌｅｎｎａｒｄ－Ｊｏｎｅｓ势所描述的惰性气体团簇ＬＪｎ进行了比较     

基态HO2根势能面特性分析和HO2束缚态的计算

对ＤＭＢＥⅣ势能面进行了特性分析，图示出了其主要地貌特征。利用该势能面进行ＨＯ２（总角动量Ｊ＝０）束缚态的计算，计算结果与该势能面特性相吻合。将计算结果与实验测量值相对照，说明ＤＭＢＥⅣ势能还有待改进     

C2分子的分析势能函数与垂直电离势

用ａｂｉｎｉｔｉｏ计算方法和实验光谱数据分别导出Ｃ２分子基态Ｘ１Σ＋ｇ和激发态ｄ３Πｇ的Ｍｕｒｅｌ－Ｓｏｒｂｉｅ函数，同时用ＱＣＩＳＤ／６－３１１Ｇ方法优化出Ｃ－２、Ｃ＋２、Ｃ２＋２分子离子基态的平衡核间距，算出Ｃ２分子的垂直电离势，计算数据与实验光谱数据进行了比较。     

UO2体系的分子动力学研究

基于多体展式方法所导出的ＵＯ２（Ｚ＾３Σu＾＋）分析势能函数,用准经典的Ｍｏｎｔ ｅ－Ｃａｒｌｏ轨迹法研究了Ｏ（Ｘ＾３Ｐｇ）＋ＵＯ（０,０）的原子与分子反应动力学过程。研究结果表明：Ｏ（Ｘ＾３Ｐｇ）与ＵＯ（０,０）碰境生成稳定络合物ＵＯ２（Ｘ＾３Σｕ＾＋）,该反应是无阈能反应;而Ｕ（Ｘ＾５Lu）与Ｏ２（０,０）碰撞不能生成稳定络合物,主要发生交换反应：Ｕ（Ｘ＾５Ｌu）＋Ｏ２（０,０）→ＵＯ（Ｖ＾１,Ｊ     

NeH2^＋体系的势能函数与动力学研究（二）Ne＋H2^＋（v）与NeH^…

基于多体项展式方法所导出的ＮｅＨ２＾＋（＾２Σ＾＋）的分析势能函数，用准经典的Ｍｏｎｔｅ－Ｃａｒｌｏ轨迹法研究了Ｎｅ＋Ｈ２＾＋（ｖ）→ＮｅＨ＾＋＋Ｈ、Ｎｅ＋Ｈ２＾＋（ｖ）→Ｎｅ＋Ｈ＋Ｈ＾＋和ＮｅＨ＾＋＋Ｈ→Ｈ２＾＋＋Ｎｅ的反应动力学。研究表明，前两个反应在Ｈ２＾＋的初始振动基态时的反应阈能分别为１６ｋｃａｌ·ｍｏｌ＾－１和６２ｋｃａｌ·ｍｏｌ＾－１，反应阈能随振动量子数增大而减小。最后一个反应是无     

氦原子与H2、D2、T2分子碰撞角分布的理论计算

采用Ｔａｎｇ－Ｔｏｅｎｎｉｅｓ势模型,用密耦法求解散射方程,计算了入射能量Ｅ＝０．０５ｅＶ时,００－００,００－０２转动激发截面,进而计算了Ｈｅ－Ｈ２,Ｈｅ－Ｄ２和Ｈｅ－Ｔ２碰撞的角分布。     

NeH~＋_2体系的势能函数与动力学研究（二）Ne＋H~＋_2（v）与NeH~＋（v＝o）＋H的反应截面计算

基于多体项展式方法所导出的ＮｅＨ＋２（２Σ＋）的分析势能函数〔６〕，用准经典的Ｍｏｎｔｅ－Ｃａｒｌｏ轨迹法研究了Ｎｅ＋Ｈ＋２（ｖ）→ＮｅＨ＋＋Ｈ、Ｎｅ＋Ｈ＋２（ｖ）→Ｎｅ＋Ｈ＋Ｈ＋和ＮｅＨ＋＋Ｈ→Ｈ＋２＋Ｎｅ的反应动力学。研究表明，前两个反应在Ｈ＋２的初始振动基态时的反应阈能分别为１６ｋｃａｌｍｏｌ－１和６２ｋｃａｌｍｏｌ－１，反应阈能随振动量子数增大而减小。最后一个反应是无阈能的     

双原子分子振动能级半量子化及量子理论的数值计算

利用半经典量子化理论和量子理论，分别采用Ｌｅｎｎａｒｄ－Ｊｏｎｅｓ势和Ｍｏｒｓｅ势，对双原子分子的振动能级给予了较为详细的数值计算。并将所得结果与部分实验数据进行比较，峄双原子分子的振动能级作了较为详细的 讨论。     

低能He—H2（D2,T2）碰撞分波截面的计算

采用Ｔａｎｇ－Ｔｏｅｎｎｅｓ势模型，当入射氨原子能量是Ｅ＝０．０５ｅＶ，计算了Ｈｅ－Ｈｅ－Ｈ２（Ｄ２，Ｔ２）弹性分波截面和非弹性激发分波截面随量子数的变化。     

NO2体系的分子动力学研究（二）

基于多体展式方法所导出的ＮＯ２（Ｘ２Ａ１）的分析势能函数［１］，用准经典的Ｍｏｎｔｅ－Ｃａｒｌｏ轨迹法研究了Ｏ（３Ｐｇ）＋ＮＯ（Ｘ２ＩＩ）和Ｎ（３Ｄｕ）＋Ｏ２（Ｘ３∑－ｇ）的原子与分子反应动力学。研究结果指出，反应Ｏ（３Ｐｇ）＋ＮＯ（Ｘ３ＩＩ）→ＮＯ２（Ｘ２Ａ１）的阈能值约为２ｋｃａｌ／ｍｏｌ，而反应Ｎ（２Ｄｕ）＋Ｏ２（Ｘ３∑－ｇ）→ＮＯ（Ｘ３ＩＩ）＋Ｏ（３Ｐｇ）是无阈能的。对后者，随初始平动能的增加，反应产物的向前散射减少，在１６０ｋｃａｌ／ｍｏｌ时达到极小，之后又缓慢增加。     

BS分子α带系A~2Π→X~2Σ~+及γ带系C~2Π→X~2Σ~+法兰克-康登因子的计算

在双原子分子核运动的薛定谔方程中，计入分子的振转相互作用项，在Ｍｏｒｓｅ势近似下得出的波函数除与振动量子数有关外，还与转动量子数有关。本文用该波函数编程计算了ＢＳ分子α带系Ａ２Π→Ｘ２Σ＋及γ带系Ｃ２Π→Ｘ２Σ＋法兰克－康登（Ｆｒａｎｃｋ－Ｃｏｎｄｏｎ）因子（简称Ｆ－Ｃ因子）。计算中转动量子数的取值由Ｊ＝０至Ｊ＝２００，其结果适用于低温、高温和强激波条件。     

BS分子a带系A^2П→X^2∑＋及γ带系C^2П→X^2∑^＋法兰克—康登因子？…

在双原子分子核运动的薛定谔方程中，计入分子的振转相互作用项，在Ｍ势近似下得出的波函数除与振动量子数有关外，不写转动量子数有关。本和该波函数编程了ＢＳα带系Ａ＾２П→Ｘ＾２∑＾＋法兰克－康登因子。计算中转动量子数的取值由Ｊ＝０至Ｊ＝２００，其结果适用于低温、高温和强激波条件。     

HeH^＋2体系的势能函数与动力学研究（二）He＋H^＋2（V）与HeH^…

基于多体项展式方法所导出的ＨｅＨ（２∑＋）的分析势能函数［８］，用准经典的Ｍｏｎｔｅｃａｒｌｏ轨迹法研究了Ｈｅ＋Ｈ（Ｖ）→ＨｅＨ＋＋Ｈ、Ｈｅ＋Ｈ（Ｖ）→Ｈｅ＋Ｈ＋Ｈ＋和ＨｅＨ＋＋Ｈ→Ｈ＋Ｈｅ的反应动力学。研究结果指出，前两个反应在Ｈ的初始振动基态时的反应阈能分别为２０ｋｃａｌ·ｍｏｌ－１和６５ｋｃａｌ·ｍｏｌ－１，反应阈能随振动量子数增大而减小。最后一个反应是无阈能的。     

S_2分子的从头算研究

采用了微机版的Ｇａｕｓｓｉａｎ９４程序以及ＣＩＳ方法和ＣＩＳ－ＭＰ２方法对Ｓ２分子的部分单重激发态和三重激发态进行了研究。计算表明：ＣＩＳ方法和ＣＩＳ－ＭＰ２方法对Ｓ２分子的２低单重激发态基本适用,可以给出较为满意的结果。但对Ｓ２分子的三重激发态和高激发态,ＣＩＳ方法还需进上步改进。     

利用半经典微扰理论研究H—H2O的非弹性碰撞—H2O的转动,振动…

利用半经典微扰理论计算了Ｈ－Ｈ２Ｏ的非弹性碰撞中的Ｈ２Ｏ分子的转动，振动激发，并将计算结果与实验结果进行了比较。计算中采用了Ｓｃｈａｔｚ和Ｅｌｇｅｒｓｍａ的半经验势能面，水分子的势函数，包括了谐振和非谐振力函数，在考虑了振动和转动耦合的情况焉，通过半经典微扰理论来确定水分子末态的振动量子数。     

Ar—O2混合气体输运系数的量子力学计算

基于物理力学理论，利用量子力学的ＩＯＳＡ方法和Ａｒ－Ｏ２等效势模型计算了Ａｒ在Ｏ２气体中的扩散系数和粘滞系数，并将计算值与实验值进行了比较。