| 查看: 1304 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
gaussion09使用已有1人参与
|
||
|
3原子以上咋消负频，算了好几天了，没有进展论文很着急，望大神指导！！！ 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
研究发现一个新的临界系数，独立于临界压缩因子
已经有7人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有213人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
gaussian激发态计算出错，求帮助
已经有4人回复- 请教高斯09算溶剂效应及能量校正问题，多送金币+小红花
已经有19人回复
- gaussian 09 做好输入文件后总是关联不上，总说找不到文件，请高手指点！
已经有16人回复
ra2ghgzh
木虫 (正式写手)
学痴
- 应助: 113 (高中生)
- 金币: 2136.9
- 散金: 5
- 红花: 33
- 帖子: 773
- 在线: 444.6小时
- 虫号: 3739328
- 注册: 2015-03-15
- 性别: GG
- 专业: 原子和分子物理
5楼2017-06-05 11:41:57
ra2ghgzh
木虫 (正式写手)
学痴
- 应助: 113 (高中生)
- 金币: 2136.9
- 散金: 5
- 红花: 33
- 帖子: 773
- 在线: 444.6小时
- 虫号: 3739328
- 注册: 2015-03-15
- 性别: GG
- 专业: 原子和分子物理
梦的翅膀、2楼2017-05-31 09:41:04
小范范1989
木虫 (著名写手)
- 应助: 123 (高中生)
- 金币: 3821.3
- 散金: 968
- 红花: 243
- 帖子: 1979
- 在线: 452.8小时
- 虫号: 2013403
- 注册: 2012-09-19
- 性别: GG
- 专业: 原子和分子物理
|
避免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 |
3楼2017-05-31 10:08:21
4楼2017-06-05 10:34:54