吸附能计算求助

小书童

sob老师，我对一个吸附的构型对不同自旋多重度1  3  5进行能量计算，想找出最低能量的稳定体，图中左边是利用OPT的结果输出的，右边的是能量计算的结果输出的，计算的方法和基组分别为B3PW91 6-31G(d)（OPT）  B2PLYP def2tzvp（能量）    从OPT结果中看3最低，从能量结果看1最低。3多重度的构型是想要的 ，迷茫了不知道那个对了！

sobereva

你弄得太混乱。
你先用B3PW91/6-31G*分别对1、3、5重态优化。每个优化完了，保存成新的单点任务的.gjf，然后用B2PLYP/def2TZVP算单点。比较单点能。

如果要考虑ZPE，B3PW91/6-31G*优化的时候写成opt freq，到时候把freq部分输出的ZPE加到B2PLYP/def2TZVP单点能上。

小书童

sobereva 发表于 2015-10-14 15:36
你弄得太混乱。
你先用B3PW91/6-31G*分别对1、3、5重态优化。每个优化完了，保存成新的单点任务的.gjf，然 ...

老师我表述的不太好，我就是您第二种做法做的。OPT的时候也做了Freq，最后计算了能量。也直接用OPT计算的能量进行了计算，结果显示这两中方法计算出来的结果不同，指向了不同的自旋多重度，而且从结构上看1不是我要的结构，3和5是我要的结构。所以我茫然了

sobereva

在三种自旋多重度已优化好的结构下，用B3PW91/def2TZVP算单点能看看。

从你的图上，B3PW91/6-31G*和B2PLYP/def2TZVP的能量趋势相差太大，感觉可能有问题。
注意B2PLYP的单点能别取错。

小书童

sobereva 发表于 2015-10-14 22:22
在三种自旋多重度已优化好的结构下，用B3PW91/def2TZVP算单点能看看。

从你的图上，B3PW91/6-31G*和B2PL ...

对于能量计算，B2PLYP/def2TZVP 我选取的是MP2= 后面的数值  ；对于OPT计算   B3PW91/6-31G*  我选取的是HF= 后的数值

小书童

小书童 发表于 2015-10-14 22:28
对于能量计算，B2PLYP/def2TZVP 我选取的是MP2= 后面的数值  ；对于OPT计算   B3PW91/6-31G*  我选取的是 ...

老师，对于算吸附这个问题，我需要对吸附上的构型对不同的多重度计算吗？还是高斯默认的就可以了？
您说能量差距大，是由于Hg的存在吗？我对汞采用的是SDD

liyuanhe211

小书童 发表于 2015-10-14 22:50
老师，对于算吸附这个问题，我需要对吸附上的构型对不同的多重度计算吗？还是高斯默认的就可以了？
您说 ...

Gaussian 默认肯定不行，Gaussian默认（实际上是GaussView默认）就是偶数个电子为1，奇数为2。
可以用一些经验的理论预测一下有几个单电子，要么严格的就算一下各个多重态的能量。

liyuanhe211

小书童 发表于 2015-10-14 22:28
对于能量计算，B2PLYP/def2TZVP 我选取的是MP2= 后面的数值  ；对于OPT计算   B3PW91/6-31G*  我选取的是 ...

Gaussian User's Reference：

1. Example B2PLYP Output. The B2PLYP energy appears in the output labeled as E(B2PLYP):
   
2. E2(B2PLYP) =    -0.3262340664D-01 E(B2PLYP) =    -0.39113226645200D+02

复制代码

以及“能量计算”是什么意思，如果是单点能计算，那和OPT是一回事，都是SCF后的E_total，如果是频率计算，不应该去找Free Energy的输出吗？

小书童

liyuanhe211 发表于 2015-10-15 00:29
Gaussian 默认肯定不行，Gaussian默认（实际上是GaussView默认）就是偶数个电子为1，奇数为2。
可以用一 ...

我计算的过程，用1  3  5三个自选多重度对结构进行优化和频率计算（OPT+Freq），结构优化后分别对三个优化好的结构进行了能量计算。

小书童

liyuanhe211 发表于 2015-10-15 00:30
Gaussian User's Reference：

我说的能量计算就是单点能的计算，我计算频率只是为了找到ZPE零点能修正。B2PLYP计算的单点能不是MP2后面的数值吗？