24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhaooming

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 36.2
红花: 1
帖子: 44
在线: 3.7小时
虫号: 614127
注册: 2008-09-27

[交流] 【求助】请各位高手专家看看我这个算是收敛么？

穷的当裤子了，没有金币可以回报啊，只能祝各位身体健康，工作顺利！

其次有几个问题需要讨教！
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.

[ Last edited by yjcmwgk on 2009-9-11 at 10:31 ]

回复此楼

» 猜你喜欢

祈福论文小修已经有21人回复
香港理工大学计算材料与物理课题组招聘博后已经有15人回复
物理化学论文润色/翻译怎么收费? 已经有204人回复
香港理工大学计算材料与物理课题组招聘博后已经有15人回复
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340 已经有2人回复
西班牙Institut Català de Nanociència i Nanotecnologia, CSIC-BIST招收公派博士生已经有0人回复
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料已经有0人回复
英国利物浦大学屠昕和雍雪课题组招收CSC博士生已经有26人回复
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）已经有0人回复
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建已经有0人回复
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人已经有45人回复

» 本主题相关商家推荐: (我也要在这里推广)

1楼2009-09-10 10:22:53

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

qianshan

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 940.4
帖子: 55
在线: 38.9小时
虫号: 478358
注册: 2007-12-14
专业: 理论和计算化学

★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+3,VIP+0):谢谢 9-22 08:43

我贴一个关于基组叠加误差的东东。不知和你的问题有没有关系
基组重叠误差(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(A,rc)-E(A,re)]+[E(B,rc)-E(B,re)]

E(int,BSSE)为修正后的结合能；E(def)为产物A、B结合成AB结构变形后的变形能；上标AB代表产物AB的平衡结构；rc、re分别表示产物AB结构和孤立A、B的结构。