计算化学公社

标题: 如何加快dimer计算的收敛速度？跑了一天，仍旧离收敛限比较远 [打印本页]

作者
Author: Ying_Zhang 时间: 2026-3-24 21:33
标题: 如何加快dimer计算的收敛速度？跑了一天，仍旧离收敛限比较远
各位老师，给一下意见，如何加快dimer的计算速度。
1. 先进行了CI-NEB（PBE/DZVP）的计算，提取过渡态结构。
2. 进行过渡态的振动分析（PBE/DZVP），发现仅有一个虚频（-1389.814687 (cm^-1) ），且振动模式对应反应路径，提取虚频下的正则坐标。
3. 进行过渡态的dimer计算（PBE/DZVP），将正则坐标手动输入，下面是dimer的输入文件。

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

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

作者
Author: Ying_Zhang 时间: 2026-3-25 17:50

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

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

感谢sob老师回复

作者
Author: northern_mar 时间: 2026-3-25 18:47
本帖最后由 northern_mar 于 2026-3-25 19:01 编辑

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

作者
Author: Ying_Zhang 时间: 2026-3-25 20:41

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

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

作者
Author: sobereva 时间: 2026-3-26 03:48

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

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

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

作者
Author: Ying_Zhang 时间: 2026-3-26 04:24

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

哦哦好的感谢sob老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)