应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】 优化结构后计算红外的问题
81/1返回列表
查看: 925  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

longyun80

木虫 (小有名气)

[交流] 【求助】 优化结构后计算红外的问题

我把优化后结构用GV导出后计算频率，在最后的结果中怎么会出现这句话？高手帮忙解释下，我刚开始学习，不太懂，谢谢！
  Item               Value     Threshold  Converged?
Maximum Force            0.000541     0.000450     NO
RMS     Force            0.000088     0.000300     YES
Maximum Displacement     0.197210     0.001800     NO
RMS     Displacement     0.043308     0.001200     NO
Predicted change in Energy=-7.253577D-05
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

[ Last edited by lei0736 on 2009-3-29 at 09:15 ]
回复此楼

» 猜你喜欢
1楼 2009-03-28 15:25:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamddba

铁杆木虫 (正式写手)
★ ★ ★
longyun80(金币+1):谢谢参与
lei0736(金币+2,VIP+0):谢谢 3-29 09:15
优化没收敛吧，优化需要的收敛最少前两项都是Yes，应该重新计算
赞一下(11人) 回复此楼
2楼2009-03-28 16:00:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

longyun80

木虫 (小有名气)
引用回帖:
Originally posted by iamddba at 2009-3-28 16:00:
优化没收敛吧，优化需要的收敛最少前两项都是Yes，应该重新计算

优化最后已经收敛了。四项都是YES，可是用GV导出来后再进行频率计算最后的结果就这样了！
赞一下 回复此楼
3楼2009-03-28 16:03:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nyzhaoyin

新虫 (正式写手)
★ ★
longyun80(金币+1):谢谢参与
lei0736(金币+1,VIP+0):谢谢 3-29 09:15
引用回帖:
Originally posted by longyun80 at 2009-3-28 16:03:

优化最后已经收敛了。四项都是YES，可是用GV导出来后再进行频率计算最后的结果就这样了！

你是不是在打开优化后的文件后，动了图像啊？  不要动图像  直接导出
赞一下(11人) 回复此楼
4楼2009-03-28 17:30:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

longyun80

木虫 (小有名气)
引用回帖:
Originally posted by nyzhaoyin at 2009-3-28 17:30:

你是不是在打开优化后的文件后，动了图像啊？  不要动图像  直接导出

我没有动过图像是直接导出的，而且导出的跟最后优化的不一样啊，
导出的文件：
B1             5.84968305
   B2             1.68199593
   B3             1.84555824
   B4             1.71934003
   B5             1.72004643
   B6             1.93328494
   B7             1.70705285
   B8             1.68199587
   B9             1.48156639
   B10            1.48091548
   B11            1.48206373
   B12            1.48091770
   B13            1.48156085
   B14            1.48206557
   B15            0.95727458
   B16            0.95594492
   A1             5.41312099
   A2           134.30399312
   A3           110.73677025
   A4           115.13606586
   A5           128.08156939
   A6           139.68553755
   A7             5.41083569
   A8           107.39383297
   A9           109.68774667
   A10          108.94025217
   A11          104.32809223
   A12          109.35295796
   A13          112.24966032
   A14          121.39099509
   A15          126.24013892
   D1          -105.33358002
   D2            62.49084447
   D3           -71.49841179
   D4           169.19022883
   D5             7.39914168
   D6             0.00905043
   D7             6.72768525
   D8          -113.57763628
   D9           126.23221987
   D10         -171.94574084
   D11           69.66113059
   D12          -53.13419615
   D13         -133.54395050
   D14          -48.37114041
而优化的结果是：
Variables:
B1=5.84968288
B2=1.68199543
B3=1.84555929
B4=1.71934004
B5=1.72004623
B6=1.93328534
B7=1.70705296
B8=1.68199541
B9=1.48156589
B10=1.48091543
B11=1.48206339
B12=1.48091803
B13=1.48156102
B14=1.48206581
B15=0.95727505
B16=0.95594455
A1=5.41311053
A2=134.30399324
A3=110.73671073
A4=115.13605892
A5=128.08160563
A6=139.68550775
A7=5.41082467
A8=107.393842
A9=109.68772791
A10=108.94024537
A11=104.32810347
A12=109.35294536
A13=112.24966951
A14=121.39097284
A15=126.24009325
D1=-105.3335538
D2=62.49087201
D3=-71.49835638
D4=169.19027961
D5=7.3991061
D6=0.00907922
D7=175.9179047
D8=55.61260817
D9=-64.57751007
D10=-171.94569022
D11=69.66117654
D12=-53.13412297
D13=-133.54399572
D14=-48.37117686
为什么会这样呢？？？
赞一下 回复此楼
5楼2009-03-28 18:00:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

北京一叶567

木虫 (小有名气)
★ ★ ★
longyun80(金币+1):谢谢参与
yuhuobuku(金币+2):欢迎参加讨论
优化的结构没有收敛
赞一下(2人) 回复此楼
6楼2009-03-31 00:01:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alwens

铁杆木虫 (正式写手)

老木虫
★ ★ ★
longyun80(金币+1):谢谢参与
yuhuobuku(金币+2):欢迎参加讨论
gv直接导出的肯定有问题，
需要从chk文件直接读，或者计算最后输出的自己手工导出来

更好的办法是读的文件再优化一遍
赞一下(2人) 回复此楼
万里夕阳锦背高 翻身犹恨东洋小 太公怎钓？
7楼2009-03-31 00:08:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)
★ ★ ★ ★
longyun80(金币+1):谢谢参与
yuhuobuku(金币+3):欢迎参加讨论
一方面，用gv导出时，由于对有效数字的截断，结构有细微的变化，这个比较优化和频率计算时算出的能量就知道，后几位是不一样的，解决方法是用geom=checkpoint从chk文件读入结构。
另一方面，优化时对力常数采用的是估算的方法，在一些情况下与用频率计算得到的力常数有较明显的差异。尽管如此，大多数情况下仍可以完成优化。但这时就会发现，原本4个yes，频率计算完再评价收敛结果，yes个数变少了。通常如果前两项标准还是yes的话，我们认为优化结构仍是可以接受的。但像你这样差别很大的，应重新优化，可以将所得结构用geom=checkpoint继续opt，同时加readfc参数从频率计算的chk中读取力常数。

[ Last edited by Jasminer on 2009-3-31 at 13:38 ]
赞一下(2人) 回复此楼
8楼2009-03-31 13:37:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 longyun80 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭