| 查看: 220 | 回复: 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 |
回复此楼» 猜你喜欢
香港中文大学（深圳）2026年招聘公告
已经有0人回复
-
牛来！米来！面来！
已经有8人回复
-
物理化学论文润色/翻译怎么收费?
已经有299人回复
-
2027申博咨询
已经有11人回复
-
寻求CrC2Al的cif文件，感谢
已经有0人回复
-
求问：LATP煅烧应该使用什么容器？
已经有0人回复
-
ORCA 纯新手向入门指南0.7版
已经有0人回复
-
哈尔滨工业大学（深圳）绿色合成化学团队招收27年春季或秋季入学的博士生
已经有5人回复
-
吉林大学化学学院理论化学研究所刘靖尧教授课题组招收2027级申请考核制博士研究生
已经有0人回复
-
西安交通大学 赵旭课题组 招收研究助理--电催化方向
已经有1人回复
-
西安交通大学 赵旭课题组 2026年招收博士后-电催化方向
已经有0人回复
