幾何優化過程中Displacement 收斂但Force 不收斂應如何解決

六羰基鉬 · 发表于 Post on 2026-7-24 15:57:07

使用 # opt=calcfc freq b3lyp/def2SVP em=gd3bj 計算含金屬-金屬鍵以及碳橋接配體的錯合物出現以下情況

Item             Value    Threshold  Converged?
Maximum Force          0.008604    0.000450    NO
RMS    Force          0.000911    0.000300    NO
Maximum Displacement    0.000000    0.001800    YES
RMS    Displacement    0.000000    0.001200    YES
Predicted change in Energy=       NaN
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

sgwzq · 发表于 Post on 2026-7-24 20:41:16

所以是报错了还是没算完？如果是没算完的话，继续等优化过程结束；如果算完报错提示不收敛，需要重新调初始构型、计算级别等等（参考http://sobereva.com/164）。
另外，金属-金属键用B3LYP是不是不太合适？这个1994年的泛函，拟合参数的时候只用了一些小分子的生成热，可能对过渡金属体系不是太理想。应该考虑一下对过渡金属优化的泛函（参考sobereva老师的文章：http://sobereva.com/272）

阮文升-河内国大 · 发表于 Post on 2026-7-25 08:29:27

同意楼上。另外，这个看起来像是一些原子跑到了接近直线的位置导致了数值问题，可以试试用最后一步的坐标输进去重新开始算（输入文件中写笛卡尔坐标，程序自动转为内坐标，会考虑原子是否共线的情况）

sobereva · 发表于 Post on 2026-7-25 10:57:49

Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（http://bbs.keinsci.com/thread-633-1-1.html）

我不知道你说的是什么金属，如果是两个d族、f族金属之间直接成键，绝对不能用杂化泛函，下文说了
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）
如果是桥联的情况，建议用TPSSh或PBE0

六羰基鉬 · 发表于 Post on 2026-8-5 14:58:54

本帖最后由六羰基鉬于 2026-8-5 15:19 编辑

sobereva 发表于 2026-7-25 10:57
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（http://bbs.kein ...

我想請問一下，如果是4d 金屬的情況下可以先用b3lyp/def2SVP GD3BJ 優化結構然後再用TPSSh/def2TZVP 計算單點能嗎？現在是加了cartesian 關鍵字後能用b3lyp/def2SVP GD3BJ 完成結構優化。但是用了TPSSh/def2SVP 要算很久還不收斂，而且在加了色散校正後要把opt 跟freq分開有點麻煩。

sobereva · 发表于 Post on 2026-8-5 15:32:37

六羰基鉬发表于 2026-8-5 14:58
我想請問一下，如果是4d 金屬的情況下可以先用b3lyp/def2SVP GD3BJ 優化結構然後再用TPSSh/def2TZVP 計算 ...

说不通。见下文
量子化学研究中切换泛函应当注意的问题
http://sobereva.com/415（http://bbs.keinsci.com/thread-9878-1-1.html）
建议都用PBE0（如果牵扯到金属-有机反应、弱相互作用，还应当带DFT-D3校正）

优化def2-SVP、单点def2-TZVP，这点是没问题的

加了色散校正後要把opt 跟freq分開有點麻煩

并不需要分开。opt freq任务带上色散校正的关键词，则opt和freq过程都会带色散校正

六羰基鉬 · 发表于 Post on 2026-8-5 15:49:25

sobereva 发表于 2026-8-5 15:32
说不通。见下文
量子化学研究中切换泛函应当注意的问题
http://sobereva.com/415（http://bbs.keinsci. ...

因爲我在老師的這篇http://sobereva.com/210 看到用TPSSh 要分開opt 跟freq做。所以用PBE0 就只需要寫em=gd3bj 然後直接Opt+freq了對吧？
那我試試PBE0, 先謝過老師撥冗回覆了

sobereva · 发表于 Post on 2026-8-5 15:54:26

六羰基鉬发表于 2026-8-5 15:49
因爲我在老師的這篇http://sobereva.com/210 看到用TPSSh 要分開opt 跟freq做。所以用PBE0 就只需要寫em= ...

TPSSh要分开是因为用了IOp，IOp没法传给之后的子任务
PBE0带D3是那样（记得PBE0的关键词是PBE1PBE）

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 幾何優化過程中Displacement 收斂但Force 不收斂應如何解決

评分 Rate