| 查看: 221 | 回复: 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 |
回复此楼» 猜你喜欢
厦门大学侯旭团队招聘青年教师、博士后及博硕士研究生
已经有2人回复
-
北京邮电大学物理学院老师招收博士生
已经有3人回复
-
物理化学论文润色/翻译怎么收费?
已经有175人回复
-
哈工大深圳 表界面物理化学方向（胶体与界面，光谱探测研究） 博士招生
已经有0人回复
-
元素硫的化合价和氧化态
已经有1人回复
-
江苏科技大学2026年招聘启事
已经有1人回复
-
Comsol 6.0 桶内肥料产热产湿过程模拟求助
已经有1人回复
-
International Journal of Nanoscience属于中科院4区吗？
已经有0人回复
-
澳大利亚南昆士兰大学量子点课题组招收博士生
已经有0人回复
-
香港中文大学（深圳）2026年招聘公告
已经有0人回复
-
大湾区大学刘天辉课题组招收2027级博士生（与中山大学物理学院联合培养）
已经有65人回复
