如何加快dimer计算的收敛速度？跑了一天，仍旧离收敛限比较远

Ying_Zhang

各位老师，给一下意见，如何加快dimer的计算速度。
1. 先进行了CI-NEB（PBE/DZVP）的计算，提取过渡态结构。
2. 进行过渡态的振动分析（PBE/DZVP），发现仅有一个虚频（-1389.814687 (cm^-1) ），且振动模式对应反应路径，提取虚频下的正则坐标。
3. 进行过渡态的dimer计算（PBE/DZVP），将正则坐标手动输入，下面是dimer的输入文件。

请问大家有什么建议，可以提升收敛速度，跑了一天，仍旧离收敛限比较远。

sobereva

DZVP-MOLOPT-SR-GTH绝对不要简单写成DZVP，有极大歧义
你的盒子和体系完全不对应，Multiwfn载入你的inp文件后进入主功能0一看便知，这是最显而易见的硬伤

Ying_Zhang

sobereva 发表于 2026-3-25 02:59
DZVP-MOLOPT-SR-GTH绝对不要简单写成DZVP，有极大歧义
你的盒子和体系完全不对应，Multiwfn载入你的inp文 ...

感谢sob老师回复，这个inp中的盒子，是我基于CI-NEB计算，提取的过渡态结构，拿Multiwfn生成dimer的输入文件时，它自动生成的。下面是过渡态的结构文件(TS.xyz)。那请问一下老师：
1）这种体系和盒子不一致对计算有什么影响？
2）如何修改这个盒子的尺寸，与计算体系适配？

感谢sob老师回复

northern_mar

Ying_Zhang 发表于 2026-3-25 17:50
感谢sob老师回复，这个inp中的盒子，是我基于CI-NEB计算，提取的过渡态结构，拿Multiwfn生成dimer的输入 ...

2）手动打开xyz文件，输入晶胞信息。下面的是我之前做的，在第二行加入晶胞信息即可。
   166
Tv_1:  16.246300   0.000000   0.000000 Tv_2:   0.000000  20.821800   0.000000 Tv_3:   0.000000   0.000000  23.751900
Zn      1.62463000      0.00000000      0.00000000
Zn      0.00000000      0.00000000      2.81388759
O       1.62463000      1.99160517      0.00000000
O       0.00000000      1.99160517      2.81388759

Ying_Zhang

northern_mar 发表于 2026-3-25 18:47
2）手动打开xyz文件，输入晶胞信息。下面的是我之前做的，在第二行加入晶胞信息即可。
   166
Tv_1:   ...

哦哦 感谢 我直接在inp文件修改了

sobereva

Ying_Zhang 发表于 2026-3-25 17:50
感谢sob老师回复，这个inp中的盒子，是我基于CI-NEB计算，提取的过渡态结构，拿Multiwfn生成dimer的输入 ...

由于CI-NEB计算过程中晶胞不会变，所以最简单的做法就是把CI-NEB的inp文件里的晶胞信息直接复制到Multiwfn新产生的dimer任务的输入文件里

盒子信息不对，则体系的周期性、所处状态完全都是错的，结果毫无意义，并且此时体系处于奇怪的电子结构，往往SCF也极难收敛

Ying_Zhang

sobereva 发表于 2026-3-26 03:48
由于CI-NEB计算过程中晶胞不会变，所以最简单的做法就是把CI-NEB的inp文件里的晶胞信息直接复制到Multiwf ...

哦哦 好的 感谢sob老师