| 查看: 228 | 回复: 1 | |||
| 【奖励】 本帖被评价1次，作者waterdp510增加金币 1 个 | |||
| 当前主题已经存档。 | |||
[资源]
Gaussian-1 (G1)方法的处理步骤
|
|||
|
第一步：采用HF/6-31G(d)产生初始的几何构型和频率分析得到零点能ZPE，矫正因子0.8929 第二步：从上一步的优化结果开始，采用MP2(Full)/6-31G(d)进行几何优化。所得几何构型用于后面的计算 第三步：计算基态能量Ebase，在上一步得到的几何结构上采用MP4/6-311G(d，p)计算。得到的数值在后面进一步矫正 第四步：增加弥散函数，采用MP4/6-311G+(d，p)计算基态能量，与上一个数值比较得到dE+ 第五步：增加高级极化函数，采用MP4/6-311G(2df，p)计算基态能量，与第三步的数值比较得到dE2df。如果该数值为正，则设该项为零 第六步：采用QCISD(T)/6-311G(d，p)计算基态能量，差值为dEQCI 第七步：矫正第六步的结果，dEHLC=-0.00019na + -0.00595nb，其中na，nb是处于alpha和beta自旋状态的电子的数目 这样，就得到的G1能量 EG1 = Ebase + dE+ + dE2df + dEQCI + dEHLC + ZPE 这样得到的EG1和QCISD(T)/6-311+G(2df，p)得到的结果近似，但速度要快的多。 Gaussian-2 (G2)方法的处理步骤 在G1的基础上，增加处理步骤 第八步：运行MP2/6-311+G(3df，2p)能量计算， dEG2 = dE+2df - dE+ - dE2df + dE3d2p 将G1方法中的2df项进行修正，由于所需要的MP2计算可以在前面找到，最终的dEG2的计算可以表示为：dEG2 = E(8) - E(5，MP2) - E(4，MP2) + E(3，MP2) 其中数字代表进行的步骤，后面的方法为该步骤中该理论的能量值。 第九步：将G1中的dEHLC修正，增加0.00114nb，记为dHLC G2能量为 EG2 = EG1 + dEG2 + dHLC |
回复此楼» 猜你喜欢
求贤纳士 | 西交利物浦大学丁理峰课题组招收计算化学博士生
已经有55人回复
-
欢迎一志愿报考华北电力大学数理学院 物理/光电方向研究生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有85人回复
-
英国利物浦大学屠昕和雍雪课题组招收CSC博士生
已经有58人回复
-
英国利物浦大学龙超老师招收电气工程专业智能电网方向CSC博士生
已经有32人回复
-
催化剂的作用原理或者催化剂失活机理相关的课件
已经有2人回复
-
中科院大连化物所潘秀莲团队研发的OXZEO双功能催化剂
已经有0人回复
-
Materials Studio 版权求助
已经有0人回复
-
中英对照版ORCA6.1.1版本manual
已经有0人回复
-
综述写作方法交流｜从选题到投稿前自查的学习路径
已经有3人回复
50