磁电子学器件应用原理近代物理实验 (75).doc

实验6-4二氟二氯甲烷分子轨道的从头算法

操作流程

用Chem3D建立分子构型

打开Chem3D程序，新建一个文件，在空白处键入分子式，回车后能看到出现了分子的初始构型。在View菜单下的setting中可以设置分子结构的显示方式，一般用BallStick。在Tools菜单下的showmodeltable下可以选择显示分子中各原子的初始坐标。将此坐标记录下来。

对分子构型做能量的初步优化

在Mechanics菜单下选择MinimizeEnergy，参数设置见下方Summary。参数设置好以后点选Run，则进行分子构型的能量的初步优化。该过程至少重复2次。构型优化结束后，记录分子中各原子坐标。

3.对分子构型做进一步优化

选择Gaussian菜单下的CreateInputFile,各参数设置如下方Summary所示。设置好以后点击Create生成.gjf文件（可取名为input.gjf），这是Gaussian09程序的输入文件。

打开Gaussian09程序，从file菜单下打开刚保存的.gjf文件，检查一下是否有问题。没有问题了就可以run了。计算结果为input.out文件。

打开input.out文件可以看到优化结果，找到最后一组坐标，这是最终优化结果，将此坐标记录。然后将此坐标填写到input.gjf中（替换原来的坐标），修改方法和基组为HF/6-311++G（d，p）进行第二次优化，优化结束后在输出文件中找到最后一组坐标，得到最终优化结果。在实验报告中对比所有记录的坐标，比较分析计算结果。

4.计算分子轨道波函数和电离能

按input.gjf格式建一个文档（可以直接用Gaussian09程序打开input.gjf进行修改）。如利用HF方法6-31G基组做它的位置空间波函数，则各命令行如下:

%Chk=HF631G.chk（生成.chk文件，用来计算电子密度分布）

#HF/6-31G6D10FSCF=TIGHTPOP=FULLGFINPUT

（计算分子轨道波函数和轨道电离能的命令行）

CF2Cl2（分子名称）

01（电荷和自旋）

C0.0000000.0000000.000000

F0.0000001.0807940.788312

F0.000000-1.0807940.788312

Cl1.4758330.000000-0.993804（原子坐标）

Cl-1.4758330.000000-0.993804

其他方法基组类似。输入文件改好以后就可以直接利用Gaussian09运行了，计算结束后给出一个相应的输出文件(可命名为hf631g.out),在最后一组结构下面给出了HF/6-31G方法的基组，如下：

S61.00

0.3047524880D+040.1834737130D-02

0.4573695180D+030.1403732280D-01

0.1039486850D+030.6884262220D-01

0.2921015530D+020.2321844430D+00

0.9286662960D+010.4679413480D+00

0.3163926960D+010.3623119850D+00

SP31.00

0.7868272350D+01-0.1193324200D+000.6899906660D-01

0.1881288540D+01-0.1608541520D+000.3164239610D+00

0.5442492580D+000.1143456440D+010.7443082910D+00

SP11.00

0.1687144782D+000.1000000000D+010.1000000000D+01

****

......

紧接其后给出了CF2Cl2分子占据和非占据轨道及其能量特征值（电离能），如下：

OrbitalSymmetries：

Occupied(B1)(A1)(B2)(A1)(A1)(B1)(A1)(A1)(B1)(A2)

(B1)(B2)(A1)(A1)(B2)(A1)(B1)(A1)(B2)(B1)

(A1)