| 查看: 182 | 回复: 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 |
回复此楼» 猜你喜欢
夏威夷大学生物胶体与界面实验室诚招全奖博士生和博士后
已经有1人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有200人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
求助立方晶体In(OH)3的CIF卡片
已经有2人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
[电子教材]天津大学物理化学（第七版）上、下册
已经有1人回复
-
阴离子交换膜电解二氧化碳还原
已经有0人回复
-
分子动力学模拟合作需求
已经有0人回复