24小时热门版块排行榜    

查看: 1601  |  回复: 4

梦的翅膀、

新虫 (初入文坛)

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

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

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

ra2ghgzh

木虫 (正式写手)

学痴

【答案】应助回帖

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

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

我永远当不了老板
2楼2017-05-31 09:41:04
已阅   回复此楼   关注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的回帖

梦的翅膀、

新虫 (初入文坛)

送红花一朵
引用回帖:
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的回帖
相关版块跳转 我要订阅楼主 梦的翅膀、 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Z7qfCsdIiD1T 2026-11-10 5/250 2026-11-12 02:45 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 Z7qfCsdIiD1T 2026-11-10 3/150 2026-11-12 02:34 by 9fOTqj074tIg
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 02:23 by 9fOTqj074tIg
[基金申请] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 02:14 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 G2rhSvqvhNzl 2026-11-10 4/200 2026-11-12 01:41 by 9fOTqj074tIg
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 YHiw1gou1Eda 2026-11-10 4/200 2026-11-12 01:17 by 9fOTqj074tIg
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 FChKX8meu1ds 2026-11-09 6/300 2026-11-12 01:06 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 01:06 by 9fOTqj074tIg
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 FChKX8meu1ds 2026-11-09 6/300 2026-11-12 00:55 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 00:55 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 6/300 2026-11-12 00:55 by 9fOTqj074tIg
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +7 cqv1zKnzd7hr 2026-11-09 10/500 2026-11-12 00:35 by 9fOTqj074tIg
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 0TEMgxdT1eIw 2026-11-09 8/400 2026-11-12 00:33 by 9fOTqj074tIg
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 0TEMgxdT1eIw 2026-11-09 8/400 2026-11-12 00:33 by 9fOTqj074tIg
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 zSbvVX9JFl9J 2026-11-09 7/350 2026-11-12 00:30 by 9fOTqj074tIg
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 3/150 2026-11-11 11:41 by gRND4NVAZ8VG
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
信息提示
请填处理意见