计算化学公社

标题: 吸附能计算求助 [打印本页]

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

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

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

作者
Author: 小书童 时间: 2015-10-14 18:39

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

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

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

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

作者
Author: 小书童 时间: 2015-10-14 22:28

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

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

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

作者
Author: 小书童 时间: 2015-10-14 22:50
本帖最后由小书童于 2015-10-14 22:52 编辑

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

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

作者
Author: liyuanhe211 时间: 2015-10-15 00:29

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

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

作者
Author: liyuanhe211 时间: 2015-10-15 00:30
本帖最后由 liyuanhe211 于 2015-10-15 00:32 编辑

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

Gaussian User's Reference：

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

复制代码

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

作者
Author: 小书童 时间: 2015-10-15 08:28

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

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

作者
Author: 小书童 时间: 2015-10-15 08:31

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)