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

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » 量化图形可视化软件 » 高斯优化结构震荡严重
【调剂】北京石油化工学院2024年16个专业接受调剂
101/1返回列表
查看: 2525  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

节能灯希望

禁虫 (小有名气)
本帖内容被屏蔽

» 猜你喜欢
1楼2020-05-07 14:13:34
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
节能灯希望: 金币+12, ★★★★★最佳答案 2020-05-11 22:28:01
从你这个电子能量来看，分子也不是很大。
那过渡态这个样子大概率是初始结构不够好导致的。要不可能就是该加溶剂模型的没加或者不该加溶剂模型的加了。之类的问题。
赞一下 回复此楼
2楼2020-05-07 14:51:47
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

节能灯希望

禁虫 (小有名气)
本帖内容被屏蔽
3楼2020-05-07 17:59:56
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

节能灯希望

禁虫 (小有名气)
本帖内容被屏蔽
4楼2020-05-07 23:39:06
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

★ ★ ★
zhou2009: 金币+3 2020-05-08 15:21:49
溶剂模型肯定会有影响。不过90%以上的过渡态不收敛都是因为初始结构不够好导致的。 可以尝试调整初始结构。或者你先把这个过渡态两侧的中间体找到。根据中间体结构再去猜测过渡态结构。
赞一下 回复此楼
5楼2020-05-08 14:09:06
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

节能灯希望

禁虫 (小有名气)
本帖内容被屏蔽
6楼2020-05-08 15:34:27
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)
引用回帖:
6楼: Originally posted by 节能灯希望 at 2020-05-08 02:34:27
老师，不是过渡态振荡，而是前一步优化就出现了振荡，气相中优化正常，加了溶剂模型，优化不出来。
...

这个状况也不少见。可以试试别的溶剂模型看看。
赞一下 回复此楼
7楼2020-05-08 16:04:34
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

节能灯希望

禁虫 (小有名气)
本帖内容被屏蔽
8楼2020-05-08 16:39:42
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

对付几何优化不收敛。而你不想换溶剂模型那就只能用常规手段来一套。

1        OPT=GDIIS (或者OPT=RFO)

高斯的默认优化算法是OPT=GEDIIS ，在优化比较“soft” 比较巨大的体系时候可能不好用。没准换成别的优化算法可以帮助收敛。
2        OPT=CALCALL (或者 OPT=RECALC=3)
OPT=CACALL这个关键词帮助高斯在优化的步骤中进行严格的HESSIAN矩阵计算。这样会帮助收敛。但是慎用，对于原子很多的体系每一步都计算HESSIAN其实就是每一步都算一次频率一样。会让计算总时长翻十倍都不止可能。只有在非常难收敛的时候再考虑使用。OPT=RECALC=3 意思是每三步计算一次HESSIAN，会比CALCALL快很多。G16中支持。
3        OPT=MAXCYC=200
这个就是增加优化步数的上限。这个仅仅对于原子数小于10个的体系使用才有效。因为原子数目比较少的时候，高斯的默认步数太小了 （最大步数=2倍的冗余内坐标数目）
4        OPT(MAXSTEP=10, NOTRUST)   (默认值是30 * 0.01Bohr)
这个关键词是缩短优化过程每一步的最大步长。对于优化中反复出现震荡的体系可能会有效。（查看优化进程是否出现振荡，查看本文GaussView使用部分3.1）NOTRUST的意思是让优化中每一步的步长上限都是10.不写这个的话从第二个点开始步长又回到30.
5        OPT=LOOSE
对于有的很大的体系，就是很难收敛这时候可以释放放宽收敛限。但是如果要进行频率计算。则这个做法无效。因为频率需要很准的结构，一般至少是OPT=TIGHT 甚至要OPT=VERYTIGHT才行
6    以上提到的方法可以各种组合，帮助收敛。比如
        OPT(MAXSTEP=5, NOTRUST, RECALC=3,GDIIS)
  而且如果分子又高对称性也最好打破对称性去算。要不然也容易振荡。
赞一下(2人) 回复此楼
9楼2020-05-08 17:08:29
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wzhsun

新虫 (小有名气)
opt=gdiis

发自小木虫IOS客户端
回复此楼
10楼2020-07-08 22:44:01
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 节能灯希望 的主题更新
101/1返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 科研好难啊，想退学去打工。 +7 byywnyl 2024-04-27 8/400 2024-04-28 06:49 by lizhiweiouc
[教师之家] 博士论文被抄袭 +27 和尚敲小木鱼 2024-04-22 44/2200 2024-04-27 21:57 by ay_guobo
[硕博家园] 博士白读了 +49 Da_Meng_Zi 2024-04-21 54/2700 2024-04-27 20:12 by 慕容家族
[基金申请] 优青已经开始了么？ +3 majunge000 2024-04-27 3/150 2024-04-27 18:51 by spring_1987
[考博] 真的好想读博！ +16 wangzhe_bs 2024-04-22 23/1150 2024-04-27 17:10 by @tanzelin
[留学DIY] 10043以后咋出去啊 +4 RuiLove 2024-04-22 4/200 2024-04-27 16:03 by 鱼翔浅底1
[论文投稿] Nature一直在编辑手里，考虑好几天了，是悬了吗 +13 彩虹初见 2024-04-24 13/650 2024-04-27 15:24 by cuan118
[教师之家] 大学直属学院卸任的副院长退休后还享受副处级的养老待遇吗？ +5 苏东坡二世 2024-04-27 5/250 2024-04-27 14:54 by zpv2011
[考博] 25年博士申请 +7 Changzixuan 2024-04-25 13/650 2024-04-27 14:11 by 我属驴核动力驴
[有机交流] 环肽的合成 +3 徐来不惊 2024-04-25 6/300 2024-04-27 13:54 by xxzhangci
[基金申请] 两类问题算是白选了~ +8 jurkat.1640 2024-04-23 13/650 2024-04-27 12:03 by 淀粉搬运工
[有机交流] 如何分离原料与产品 5+4 qwerasdf587 2024-04-24 13/650 2024-04-26 23:18 by 119966984
[论文投稿] LWT投 +3 AChen92 2024-04-26 3/150 2024-04-26 22:16 by hizifu
[考研] 学硕专硕 +5 小蜗牛* 2024-04-26 5/250 2024-04-26 16:43 by 鱼翔浅底1
[考研] 0854－0855调剂 +8 shangannum1 2024-04-21 12/600 2024-04-26 16:42 by yz仔
[硕博家园] 考研，求职还是考编? +15 xizj 2024-04-21 24/1200 2024-04-26 11:49 by Kan客
[考研] 381求调剂 +4 小刺猬987654321 2024-04-25 6/300 2024-04-26 10:57 by czl12138
[教师之家] 某种做法不行。说过几遍了。同学还那样做。再那样做就给低分 +4 河西夜郎 2024-04-24 4/200 2024-04-26 08:51 by Quakerbird
[教师之家] 大家访学都是怎么找的啊？ +3 luokereng 2024-04-22 3/150 2024-04-24 11:40 by xuechenli
[论文投稿] 期刊推荐 20+4 木颜尘ip 2024-04-22 7/350 2024-04-24 10:06 by bobvan
信息提示
关闭
请填处理意见
关闭