| 查看: 1270 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
关于优化收敛和频率计算不收敛问题 已有2人参与
|
||
|
我刚刚接触高斯，最近做了一个25个原子的化合物，结构优化收敛是4个yes,如下： Maximum Force 0.000428 0.000450 YES RMS Force 0.000103 0.000300 YES Maximum Displacement 0.000127 0.001800 YES RMS Displacement 0.000032 0.001200 YES Predicted change in Energy=-2.977102D-07 Optimization completed. 但优化好的结构计算频率只有一个yes,如下： Maximum Force 0.000914 0.000450 NO RMS Force 0.000294 0.000300 YES Maximum Displacement 0.041760 0.001800 NO RMS Displacement 0.008580 0.001200 NO Predicted change in Energy=-4.117907D-05 这个正常吗？是优化的结构不对还是频率计算有问题呢？请高人指点 |
回复此楼» 猜你喜欢
奥尔堡Xinxin Xiao课题组招收一名岗位制博士（酶基生物电化学）
已经有0人回复
-
招调剂
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有235人回复
-
26年电池方向博士申请
已经有4人回复
-
运动康复求调剂
已经有0人回复
-
中原工学院现有材料科学与工程和纳米科学与工程学术型、材料工程研究生调剂名额
已经有0人回复
-
北方工业大学储能技术02电化学储能（电池）方向第二次调剂，4月12日22:00调剂系统关闭
已经有2人回复
-
北方民族大学085602化学工程有调剂名额34个，欢迎调剂！4月13日20:00开系统
已经有0人回复
-
中南林科大化工院第三批调剂
已经有1人回复
-
华南理工大学电子显微中心博士招生（仪器、电子、机械、材料、软件/AI方向）
已经有0人回复
-
重庆文理学院2026年硕士研究生招生调剂信息 化学 、材料工程
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 3D智能数字化与3D打印的一本新书
已经有12人回复
3楼2015-08-31 09:57:35
chenxin1992
银虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 1030.1
- 散金: 126
- 红花: 11
- 帖子: 337
- 在线: 103.1小时
- 虫号: 1813050
- 注册: 2012-05-12
- 性别: GG
- 专业: 无机合成和制备化学
2楼2015-08-28 21:53:05
zhou2009
版主 (著名写手)
- QC强帖: 15
- 应助: 337 (大学生)
- 贵宾: 4.053
- 金币: 42654.6
- 红花: 153
- 帖子: 1729
- 在线: 1274.7小时
- 虫号: 676800
- 注册: 2008-12-18
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
【答案】应助回帖
感谢参与，应助指数 +1
|
Sobereva在帖子《量子化学计算中帮助几何优化收敛的常用方法》中说： RFO、GDIIS、GEDIIS这三种优化方法在走步时都需要Hessian矩阵（力常数矩阵），然而计算Hessian矩阵是很耗时的。因此，在默认情况下，Gaussian会通过价力场的方法近似估计出初始的Hessian矩阵，在每一步优化中只精确计算梯度，利用梯度对原先Hessian矩阵进行修正不断得到新的Hessian矩阵。所以，在默认情况下，优化从头到尾使用的Hessian矩阵都只是近似的（这也是为什么在优化出来的结构上做Freq往往会显示还没收敛，因为Freq用的是精确的Hessian矩阵）。当Hessian矩阵离精确值偏差较大，就会造成收敛缓慢，或者始终不收敛。为解决这个问题，可以用opt=calcfc，这会在优化的第一步使用精确计算的Hessian矩阵，但是仍可能后续优化过程中Hessian矩阵逐渐变得越来越不精确而依然收敛失败。opt=calcall则不仅在第一步，在后续的每一步中也都精确计算Hessian矩阵，这使得很多优化失败的情况都能得到解决，优化所需步数通常也会减少很多，而且能够保证最终优化结果准确（因为最终判断是否收敛时是基于精确的Hessian矩阵所得结果），但代价是每一步计算量会很大。 |
4楼2015-08-31 10:51:16
20