最新量子计算化学教学内容.doc

《最新量子计算化学教学内容.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、量子计算化学报告姓名：李丽学号：2012101217此报告，我将从赖老师布置的两道计算题来简单阐述我在这个课程中所学到的知识：1.过渡态的优化；2.在气相和液相中分子构象的变化和能量的变化。一、TSsearchingN-3-丁烯基硝酮分子内环加成反应产物选择性（反应条件：甲苯中回流，甲苯的沸点是：T=383.75K）。研究目的：化学反应的过渡态优化计算方法：3-21G(d)1.在GaussView5.0软件中绘制出合适的过渡态结构，并找到合适的键长，保存成.gif格式之后；用Gaussian09W软件打开.gif

2、文件之后，输入以下命令对过渡态进行优化：%chk=F:GaussC-N-O1.chk#opt(modredundant)hf/6-31g(d)TitleCardRequired01B216FB13F我找到的优化的键长是C—C:2.28;C—O:2.322.把优化之后的过渡态结构的.gif文件用Gaussian09W软件打开之后，输入以下命令：%chk=F:GaussC-N-O1-beautiful.chk#opt(modredundant,ts,noeigentest,calcfc)hf/3-21gfre

3、qTitleCardRequired01B216AB13A打开相应的.out文件，找到过渡态的焓值（H），熵值（S）和Cv填入下面的表格。3.在GaussView5.0软件中打开优化后的过渡态的.gif后，将两个键长调到比较大的状态，这样就到了反应物了，保存好。再用Gaussian09W软件打开并输入以下命令：%chk=F:GaussC-N-O1-reactor.chk#opthf/3-21gfreqTitleCardRequired01打开相应的.out文件，找到反应物的焓值（H），熵值（S）和Cv填入下面

4、的表格。4.重复以上的步骤计算第二种过渡态即可。（1a.u.=627.51kcal/mol）计算结果：如下热力学量(T=273.15K)RaTSaRbTSbH(a.u.)-0.51282-0.51264-0.51089-0.51079S(cal/mol/K)84.28978.60273.31176.483CV(cal/mol/K)25.57121.77319.92721.3225.由于采用软件计算得出的结论是在T=273K的温度下的数据，因此，我们要采用下面的公式来计算T=383.75K温度下的数据(R=8.31

5、4J/mol·K=1.986cal/mol·K;Cp=Cv+R)：热力学量(T=383.75K)RaTSaRbTSbH(a.u.)4.344203.674903.351303.59730S(cal/mol/K)93.6584086.680080.761084.4080CV(cal/mol/K)25.571021.773019.927021.3220速率比：计算得：分析讨论：由此可以看出，生成a产物的比率稍大于生成b产物。从表格中的数据可以看到对于焓值，我所得到的结果跟老师的有比较大的出入，不知道是什么原因？二、T

6、heenergydifferencebetween1,2-Dichloroethaneconformersinsolution研究目的：1,2-二氯乙烷两种不同的构象在液相（环己烷）和气相中能量和键长的变化。计算方法：HF/6-31+G(d)；MP2/6-31+G(d)1.在GaussView5.0软件中绘制出合适的对位交叉式和邻位交叉式结构，保存成.gif之后，用Gaussian09W软件打开.gif文件之后，输入以下命令（气相中）：HF方法%chk=F:GaussC2H4Cl2(cis).chk#opth

7、f/6-31+G(d)TitleCardRequired01MP2方法%chk=F:GaussC2H4Cl2(cis).chk#optMP2/6-31+G(d)TitleCardRequired01输入以下命令（液相中）HF方法%chk=F:GaussC2H4Cl2(cis).chk#opthf/6-31GSCRF=(dipole,Solvent=cyclohexane)TitleCardRequired013.65MP2方法%chk=F:GaussC2H4Cl2(cis)2MP2.chk#optMP

8、2/6-31+G(d)SCRF=(Solvent=cyclohexane)TitleCardRequired013.65计算结果：以上计算出的结果，打开.out文件之后即可找到自己所需要的数据，填入上面的表格即可（键长的变化和能量的变化）。6-31+G(d)键长能量C-CC-HC-Cl气相中计算方法对位交叉式HF1.515801.077581.79272-997.0328