理论与计算化学第二次上机作业

理论与计算化学第二次上机作业

说明：（1）请将输入文件、输出文件、检查点（.chk）文件压缩成一个文件发到TA邮箱，将压缩文件命名为本人学号_姓名.rar。（2）以下作业中，除了填写相应计算结果外，请同时写上计算所用的方法、基组和关键字，例如HF/6-31G* opt freq等。

1．比较下列不同方法计算得到的N2基态能量。假设N N的键长为1.098Å(单点计算)，采用6-31+g*方法。

方法、基组、关键字：hf/6-31+ g*sp geom=connectivity

SVWN or LsDA/6-31+g*sp geom=connectivity

# MP2/6-31+g*sp geom=connectivity

CISD/6-31+ g*sp geom=connectivity

B3LYP/6-31+ g*sp geom=connectivity

2. 用不同方法计算H2O分子的能量（包含零点振动能校正）、O的Mulliken电荷以及弯曲振动的频率大小，基组为6-311+g*。（优化结构+计算频率）

3. 用B3L YP/6-31++G** 计算下列反应的正反应活化能和逆反应活化能（单位：kcal/mol）

方法、基组、关键字：# b3lyp/6-31+g** opt(maxcycle=100) freq;

# b3lyp/6-31+g** opt(maxcycle=100) freq

# b3lyp/6-31++g** IRC(RCFC,Maxcycle=40,Maxpoint=15) Guess=Read Geom=AllCheck

正反应活化能=HT(TS)-HF(2)=(-114.353393-(-114.484959))* 627.5095=82.55891kcal/mol

逆反应活化能= HT(TS)-HF(1)= (-114.353393-(-114.5116286))* 627.5095 = 99.294342

kcal/mol

1 2

4. 乙炔分子（C2H2）基态的构型如下，请用CASSCF(4,4)方法和6-31+G（d）基组计算该构型下的CASSCF能量，活性空间包含两个成键的 π轨道和两个反键的π*轨道，并写出活性空间中这两组轨道的电子占据数。

CASSCF能量=-76.890920 7

Final one electron symbolic density matrix:

1 2 3 4

1 0.193883D+01

2 0.130956D-27 0.193883D+01

3 0.402247D-13 -0.922285D-15 0.611745D-01

4 -0.495090D-14 0.341570D-13 -0.668175D-28 0.611745D-01

方法、基组、关键字：# CASSCF(4,4)/6-31+G(d) Guess=(Read,Permute)