| 查看: 1157 | 回复: 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 这个正常吗？是优化的结构不对还是频率计算有问题呢？请高人指点 |
回复此楼» 猜你喜欢
26年秋季博士申请
已经有0人回复
-
推荐给英语教学者的一本单词书《金鱼单词讲义：从26个拉丁字母到106万个英语单词》
已经有59人回复
-
物理化学论文润色/翻译怎么收费?
已经有278人回复
-
推荐给教师的一本单词书《金鱼单词讲义：从26个拉丁字母到106万个英语单词》
已经有33人回复
-
核磁分析软件MestReNova打开文件时报错
已经有0人回复
-
在职博后不能申请博后基金了，那么在职博后意义何在？
已经有2人回复
-
青岛大学化学化工学院分子测量学研究院2026年招收博士研究生
已经有0人回复
-
香港科技大学（广州）诚招电催化方向博士生（2026秋入学）
已经有0人回复
-
求助Cu2+1O的CIF文件(PDF: 05-0667)
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 3D智能数字化与3D打印的一本新书
已经有12人回复
chenxin1992
银虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 1030.1
- 散金: 126
- 红花: 11
- 帖子: 337
- 在线: 103.1小时
- 虫号: 1813050
- 注册: 2012-05-12
- 性别: GG
- 专业: 无机合成和制备化学
2楼2015-08-28 21:53:05
3楼2015-08-31 09:57:35
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