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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【讨论】为何结果不同
172/2返回列表上一页12
查看: 623  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

wolfzhong

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
不知道你这个和收敛标准有没有关系。

之前我见过类似的例子，就是同一构型得到的轨道是不一样的。准确来说，你的收敛标准就那么高，你还指望它能到10**-10啊。
赞一下(1人) 回复此楼
11楼2009-12-31 17:31:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gnli

铜虫 (小有名气)
引用回帖:
Originally posted by wolfzhong at 2009-12-31 17:31:
不知道你这个和收敛标准有没有关系。

之前我见过类似的例子，就是同一构型得到的轨道是不一样的。准确来说，你的收敛标准就那么高，你还指望它能到10**-10啊。

这和收敛标准有啥关系呢？
不太明白
赞一下 回复此楼
争取早日成为合格的菜鸟
12楼2010-01-01 18:49:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gnli

铜虫 (小有名气)
引用回帖:
Originally posted by yflchx at 2009-12-31 08:58:
你仔细比较一下两次计算任务的分子结构有没有区别，往小数点后多看几位。

我对比了结构 没有不同之处 具体数据对比：

opt+freq
最后算频率时的标准输入结构
1          7             0        0.735517   -1.035369    0.532954
    2          1             0        0.226481   -1.821177    0.088264
    3          1             0        0.935641   -1.328758    1.506344
    4          6             0        2.023680   -0.718352   -0.212125
    5          6             0       -0.135194    0.205193    0.538263
    6          1             0        1.898421   -1.089523   -1.236151
    7          1             0        2.847929   -1.247469    0.276713
    8          6             0        2.103957    0.804363   -0.164013
    9          6             0        0.638829    1.238722   -0.315319
   10          1             0       -0.228193    0.540905    1.583235
   11          6             0       -1.503529   -0.147082   -0.025120
   12          1             0        2.739639    1.195758   -0.967249
   13          1             0        2.511925    1.139478    0.800014
   14          1             0        0.326012    1.173213   -1.367074
   15          1             0        0.451019    2.257303    0.042673
   16          8             0       -1.714720   -1.193470   -0.607573
   17          8             0       -2.461681    0.767821    0.098945
   18          1             0       -2.212747    1.535985    0.681465

我从gview另存得到的结构
N                  0.73551700   -1.03536900    0.53295400
H                  0.22648100   -1.82117700    0.08826400
H                  0.93564100   -1.32875800    1.50634400
C                  2.02368000   -0.71835200   -0.21212500
C                 -0.13519400    0.20519300    0.53826300
H                  1.89842100   -1.08952300   -1.23615100
H                  2.84792900   -1.24746900    0.27671300
C                  2.10395700    0.80436300   -0.16401300
C                  0.63882900    1.23872200   -0.31531900
H                 -0.22819300    0.54090500    1.58323500
C                 -1.50352900   -0.14708200   -0.02512000
H                  2.73963900    1.19575800   -0.96724900
H                  2.51192500    1.13947800    0.80001400
H                  0.32601200    1.17321300   -1.36707400
H                  0.45101900    2.25730300    0.04267300
O                 -1.71472000   -1.19347000   -0.60757300
O                 -2.46168100    0.76782100    0.09894500
H                 -2.21274700    1.53598500    0.68146500
赞一下 回复此楼
争取早日成为合格的菜鸟
13楼2010-01-01 19:03:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yflchx

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
再检查一下自恰场收敛？

结构没问题，分开算用GUESS=READ了么？
赞一下(1人) 回复此楼
14楼2010-01-01 19:52:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gnli

铜虫 (小有名气)
引用回帖:
Originally posted by yflchx at 2010-1-1 19:52:
再检查一下自恰场收敛？

结构没问题，分开算用GUESS=READ了么？

没有用 guess 没有保存 opt 的chk文件
用gview倒的结构直接算的
赞一下 回复此楼
争取早日成为合格的菜鸟
15楼2010-01-01 20:40:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gnli

铜虫 (小有名气)
引用回帖:
Originally posted by yflchx at 2010-1-1 19:52:
再检查一下自恰场收敛？

结构没问题，分开算用GUESS=READ了么？

我检查了一下
发现SCF能量后三位有差别
难道是因为初始猜测的问题造成的？
我用的是同样的泛函和机组 不应该有问题吧？

opt + freq 的最后一步 SCF 能量
SCF Done:  E(RB+HF-LYP) =  -401.692030655     A.U. after    1 cycles
freq=raman的最后一步SCF能量
SCF Done:  E(RB+HF-LYP) =  -401.692030525     A.U. after    7 cycles
赞一下 回复此楼
争取早日成为合格的菜鸟
16楼2010-01-01 20:44:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yflchx

金虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+1,VIP+0): 1-8 14:11
我想问题也许就出在这里。

因为其余的原因已经被排除了。

对于Gaussian或其它量子化学计算软件来讲，键长一般保留到小数点后三位，键角两位，能量(kJ/mol,kcal/mol)一位，频率取整。
赞一下(11人) 回复此楼
17楼2010-01-02 19:55:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gnli 的主题更新
172/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭