| 查看: 1271 | 回复: 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 这个正常吗？是优化的结构不对还是频率计算有问题呢？请高人指点 |
回复此楼» 猜你喜欢
河南理工大学化学化工学院招收学硕调剂生
已经有0人回复
-
黑龙江科技大学化工 环境工程调剂
已经有22人回复
-
物理化学论文润色/翻译怎么收费?
已经有138人回复
-
图书求助：Fluorescent Probes for Bioactive Species: Design
已经有0人回复
-
莆田学院环境与生物工程学院资源与环境专硕（环境工程方向） 开放调剂招生
已经有13人回复
-
浙江海洋大学石油化工与环境学院招收工科调剂生（学硕、专硕均可，线上面试）
已经有0人回复
-
浙江海洋大学石油化工学院接收工科调剂（学硕、专硕均可，近期线上面试）
已经有8人回复
-
浙江海洋大学石油与天然气工程专业接收工科调剂（13日早上截止，线上面试）
已经有0人回复
-
考研调剂～湖南工程学院材料与化工学院
已经有0人回复
-
同系列相似结构的单层二维COF真空层厚度该如何设置？
已经有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
