24小时热门版块排行榜    

查看: 1605  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

梦的翅膀、

新虫 (初入文坛)

[求助] gaussion09使用 已有1人参与

3原子以上咋消负频,算了好几天了,没有进展论文很着急,望大神指导!!!

发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

避免freq最后出现NO乃至虚频,或者说让优化更精确的做法有二:
(1)几何优化时用更严的收敛限,即opt里面写上tight。对于DFT来讲,建议同时搭配int=ultrafine使用更高精度的积分格点,效果会明显更好(注:从G16开始int=ultrafine已成为默认)。这样来解决出现NO和虚频对计算耗时增加不太多,不过也有可能会反复震荡很难达到这个收敛限。
注意:如果写opt=tight freq,则freq最后也会用tight标准来判断是否收敛。实际上只要能满足默认判断标准即可,不需要非得满足tight的判断标准。

(2)使用精确的Hessian矩阵。在opt里写上calcall,那么几何优化过程中每一步都会精确计算Hessian,此时几何优化和freq所用的Hessian都是一致的,因此对收敛的判断结果也是完全相同的。也就是说,只要opt最后是YES则freq最后也必然是YES,并且多数情况也可以确保无虚频。(很值得一提的是,用calcall的时候,优化任务最后会自动做振动分析,也就是说,此时实际上完全没有必要再单独用freq关键词做振动分析了。)
不过,每一步优化都精确计算Hessian太耗时。如果你之前可以已经在默认情况下优化过,可以取最后的结构重新优化,并在opt里写上calcfc,这样只在优化第一步的时候精确计算Hessian,而之后还是近似计算Hessian。这样收敛后也往往可以避免freq之后出现NO。但是如果初始结构离实际极小点太远则这么做无效,因为等到收敛时Hessian矩阵可能又偏离精确Hessian比较远了。

参考http://sobereva.com/278
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
3楼2017-05-31 10:08:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

ra2ghgzh

木虫 (正式写手)

学痴

【答案】应助回帖

感谢参与,应助指数 +1
沿着负频的振动方向调整一下原子位置不就行了?这个应该不难吧。

» 本帖已获得的红花(最新10朵)

我永远当不了老板
2楼2017-05-31 09:41:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦的翅膀、

新虫 (初入文坛)

送红花一朵
引用回帖:
2楼: Originally posted by ra2ghgzh at 2017-05-31 09:41:04
沿着负频的振动方向调整一下原子位置不就行了?这个应该不难吧。

直接加,直到负频为正吗?加到负频都变大了,但是振动坐标没变,再加振动坐标负频又变小,还继续加吗,用乘倍数在加吗?
初次接触,有些不懂的,望指教!!!

发自小木虫Android客户端
4楼2017-06-05 10:34:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ra2ghgzh

木虫 (正式写手)

学痴

引用回帖:
4楼: Originally posted by 梦的翅膀、 at 2017-06-05 10:34:54
直接加,直到负频为正吗?加到负频都变大了,但是振动坐标没变,再加振动坐标负频又变小,还继续加吗,用乘倍数在加吗?
初次接触,有些不懂的,望指教!!!
...

http://muchong.com/bbs/viewthread.php?tid=10455387
你看一下10楼的方法,我觉得挺靠谱的。
我永远当不了老板
5楼2017-06-05 11:41:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 博士申请 +4 z_l_l_666 2026-11-06 4/200 2026-11-13 06:19 by wlf1982
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 cdqGD6ftHmbl 2026-11-11 5/250 2026-11-12 12:28 by ECLTXnJnRSXf
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 cdqGD6ftHmbl 2026-11-11 4/200 2026-11-12 12:26 by ECLTXnJnRSXf
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:31 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 11:22 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:11 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 11:03 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 czAiolzR0IJC 2026-11-10 9/450 2026-11-12 11:00 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 7/350 2026-11-12 10:49 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 8/400 2026-11-12 10:16 by HajghKOMynmz
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 C6CRn66RYYXS 2026-11-10 5/250 2026-11-12 10:16 by HajghKOMynmz
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +7 zSbvVX9JFl9J 2026-11-09 8/400 2026-11-12 09:19 by HajghKOMynmz
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 C6CRn66RYYXS 2026-11-10 5/250 2026-11-12 06:03 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 6/300 2026-11-12 06:03 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 05:30 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 10/500 2026-11-12 05:19 by 9fOTqj074tIg
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 cqv1zKnzd7hr 2026-11-09 8/400 2026-11-12 05:08 by 9fOTqj074tIg
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-12 04:57 by 9fOTqj074tIg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 17:03 by TgvkbZZl9tqB
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
信息提示
请填处理意见