最近,Capecchi和Werner拟合出了Cl+H<,2>体系的非绝热势能面,被称作CW势能面,并在此势能面上使用量子的非含时方法进行了动力学计算,计算结果与实验有非常大的差别.为了探索理论与实验的差别,有必要对该反应体系进行进一步研究.于是我们参照Alexander等人拟合F+H2耦合势能面的方法以及BW2势能面的拟合方法,使用最小二乘法在高精度的量化计算基础上独立编写程序也拟合了一个势能面,它包含了四个势能项,两个自旋与轨道相互作用项.与不含自旋与轨道相互作用的BW2势能面相比,自旋与轨道相互作用在入口处使基态的能量降低了大约0.84Kcal/mol,而在过渡态附近对基态的能量基本没有影响,因此使势垒高度增加了大约0.84Kcal/mol.同时由于自旋轨道相互作用以及势能面之间的相互作用,使处于激发态的Cl原子可以通过非绝热跃迁到达基态从而参与反应.此外,与CW势能面相比较,我们拟合的势能面中自旋轨道矩阵元更具有合理性.在上述两组势能面基础之上,我们用含时波包法做了进一步的动力学计算.使用三维非绝热含时波包方法对Cl+H<,2>(HD,D<,2>)进行了反应散射计算,给出了反应几率以及积分截面.这三个体系的计算结果均表明,在低能处,自旋轨道激发态由于具有比较大的内能而反应活性比较大;在高能处,自旋轨道激发态比基态高出的那部分内能起着很小的作用,因此自旋轨道基态的反应活性大.对Cl+HD的计算则表明DCl通道的反应活性很高,这与实验结果是一致的.与Alexander等人的理论计算结果相反,Coriolis耦合对我们的计算结果影响不大,并且我们计算出的激发态活性比他们的理论结果要大很多,已经达到了可以和实验结果相比的程度,但是和实验结果仍有一些不同,差别的原因还有待进一步探索.