小木虫 --- 600万学术达人喜爱的学术科研平台

首页 >> 量子化学 >>【求助】请各位高手专家看看我这个算是收敛么?

【求助】请各位高手专家看看我这个算是收敛么?

穷的当裤子了,没有金币可以回报啊,只能祝各位身体健康,工作顺利!
其次有几个问题需要讨教!
1)关于四楼的回答,对重原子加弥散,优化好了之后,再去掉弥散是不是还可以优化好
2) 关于六楼,这是一个产物,有2-将冰片烯裂解得到,不是说产物也要优化么?我这个地方还没搞清楚呢,他们之间到底怎么作用?距离该是多远?看到有的文献上,两个分子之间的距离也挺远的,都有3.***的样子。


# opt=calcfc b3lyp/6-31+g(d,p) guess=read geom=(connectivity,allcheck)
c5h6+c2h4 opt
0 1
C
C                  1              B1
C                  2              B2    1              A1
C                  3              B3    2              A2    1              D1
H                  3              B4    2              A3    1              D2
H                  1              B5    2              A4    3              D3
H                  4              B6    3              A5    2              D4
H                  4              B7    3              A6    2              D5
C                  3              B8    2              A7    1              D6
H                  3              B9    2              A8    1              D7
H                  9             B10    3              A9    2              D8
C                  4             B11    3             A10    2              D9
H                 12             B12    4             A11    3             D10
C                 12             B13    4             A12    3             D11
H                 14             B14   12             A13    4             D12
H                  2             B15    1             A14    4             D13
H                  3             B16    2             A15    1             D14
   B1             1.35265980
   B2             2.30246529
   B3             1.50629801
   B4             3.26106574
   B5             1.08361529
   B6             1.10049415
   B7             1.10048936
   B8             1.35273514
   B9             4.44093799
   B10            1.08471020
   B11            4.12794848
   B12            1.08718789
   B13            1.33457914
   B14            1.08709798
   B15            1.08469280
   B16            1.08362044
   A1            75.57961085
   A2            72.04701325
   A3            72.85087912
   A4           126.82490662
   A5           111.91649174
   A6           111.93638862
   A7            37.04708486
   A8           105.91357717
   A9           126.21758000
   A10           82.51552139
   A11          110.93352649
   A12          113.99229221
   A13          121.81749459
   A14          126.22606783
   A15          163.87482060
   D1            -0.01363677
   D2           -84.82860615
   D3          -179.98441922
   D4          -120.54188696
   D5           120.61979512
   D6          -179.99714670
   D7           -88.60872048
   D8          -179.91065264
   D9           -82.83807043
   D10          146.76884301
   D11          -71.65458287
   D12         -137.04183122
   D13          179.96074752
   D14          179.72022832
收敛条件为:我再优化,老是这样!!!
         Item               Value     Threshold  Converged?
Maximum Force            0.000004     0.000450     YES
RMS     Force            0.000001     0.000300     YES
Maximum Displacement     0.133908     0.001800     NO
RMS     Displacement     0.020511     0.001200     NO
Predicted change in Energy=-7.806212D-07
Optimization completed on the basis of negligible forces.
    -- Stationary point found.

今日热帖

用户评论

势能面太平缓了,加gdiis参数

您就是密度泛函小卒啊,呵呵,久仰大名了。谢谢了,不知道能不能QQ加友,方便以后向您请教。
我的QQ:403725301

后两项差的也太多了。
最近我在优化的一个结构,算了几种方法都是除了第三项NO,其它都是YES,后来我对重原子再加一个弥散就全YES了,没有虚频。

Originally posted by zhaooming at 2009-9-10 10:22:
穷的当裤子了,没有金币可以回报啊,只能祝各位身体健康,工作顺利!

# opt=calcfc b3lyp/6-31+g(d,p) guess=read geom=(connectivity,allcheck)
c5h6+c2h4 opt
0 1
C
C                  1        
...
         Item               Value     Threshold  Converged?
Maximum Force            0.000004     0.000450     YES
RMS     Force            0.000001     0.000300     YES
Maximum Displacement     0.133908     0.001800     NO
RMS     Displacement     0.020511     0.001200     NO
Predicted change in Energy=-7.806212D-07
Optimization completed on the basis of negligible forces.
    -- Stationary point found. ...
不是要求tight标准的优化的话,可以算是收敛了。
on the basis of negligible forces:这句话是关键的提示,

我复制了你的gjf,放到我的GV里一看,我晕!你想算啥?
如果你想算它们之间的相互作用,距离这么远,他们之间会有相互作用吗?
如果你觉得他们之间没有相互作用,为什么要放到一起优化?

其次有几个问题需要讨教!
(1)关于四楼的回答,对重原子加弥散,优化好了之后,再去掉弥散是不是还可以优化好
(2) 关于六楼,这是一个产物,有2-将冰片烯裂解得到,不是说产物也要优化么?我这个地方还没搞清楚呢,他们之间到底怎么作用?距离该是多远?看到有的文献上,两个分子之间的距离也挺远的,都有3.***的样子。
在这上面发帖交流经验不足,不知道这样他们能不能看到,我在这祈祷了,呵呵呵
:):):):):):)

产物优化,你把两个分子分别优化就行了吧!

Originally posted by hairan at 2009-9-11 09:57:
产物优化,你把两个分子分别优化就行了吧!
楼上正解!没有你这样蛮干的。

两个分子之间的距离也挺远的,都有3.***的样子。
单位是埃?

>3.0 埃的优化俺也做过,:)

Originally posted by yjcmwgk at 2009-9-10 19:00:
我复制了你的gjf,放到我的GV里一看,我晕!你想算啥?
如果你想算它们之间的相互作用,距离这么远,他们之间会有相互作用吗?
如果你觉得他们之间没有相互作用,为什么要放到一起优化?
找原因就要像这样从源头开始,赞.

(1)两个产物分子分别优化是不是比放在一起优化时的能量大啊?
(2)有的文献上两个分子放在一起优化,有的是分开,到底该什么时候放在一起,什么时候分开啊?
(3)貌似放在一起优化也没什么好啊,算活化能什么的算放在一起的还是分开后的,貌似效果差不多啊

请问楼上你说的涉及基组叠加误差么?
师兄说基组比较大的时候基组叠加误差是可以不考虑的

我贴一个关于基组叠加误差的东东。不知和你的问题有没有关系
基组重叠误差(Basis Set Superposition Error)的校正主要是用于结合能计算。如两个水分子的二聚物的结合能,一般是先计算二聚物的能量(产物),再减去两个单个水分子(反应物)的能量,但这样计算一般会对实际的结合能进行过高的估计,这个现象就是BSSE。因为两个分子接近时,除了两个分子的相互作用导致能量降低外,每个分子的基函数对另一个分子的电子结构提供了更好的描述而导致能量的降低(根据变分原理,描述一个体系,波函数越适宜其能量越低)。很显然,使用较小的和不合适的基组(对离核较远的电子分布不能提供好的描述,特别是在非键相互作用最强的区域),BSSE会更明显。估算BSSE的一种方法是Boys and Bernardi提出的抵消校正法,从头计算中由Counterpoise关键词实现。校正时,单个组分A的计算要具有和B一样的幽灵轨道,单组分B计算也一样,这样就可以得到扣除BSSE的结合能。用公式可表述为:
E(int,BSSE)=E(AB,rc)AB-E(A,rc)AB-E(B,rc)AB+E(def)
E(def)=+
  E(int,BSSE)为修正后的结合能;E(def)为产物A、B结合成AB结构变形后的变形能;上标AB代表产物AB的平衡结构;rc、re分别表示产物AB结构和孤立A、B的结构。

猜你喜欢

学术必备
与600万学术达人在线互动!


扫描下载送金币