ZnNiLDH双层结构优化原子散开变形严重

PoorChaos · 发表于 Post on 3 day ago

本帖最后由 PoorChaos 于 2026-4-22 15:01 编辑

请教 ZnNiLDH的结构优化问题
之前是中间插层放入有机物难优化成功
现在优化未吸附有机物的基底结构，原子变化散开

opt.7z (263.49 KB, 下载次数 Times of downloads: 2) inp文件见附件
部分inp，restart，优化完的out和未优化完的都在压缩包

ZnNiLDH.inp (17.95 KB, 下载次数 Times of downloads: 5)

Uus/pMeC6H4-/キ · 发表于 Post on 3 day ago

先尝试算单点能并按http://sobereva.com/665解决SCF不收敛问题（此处可以考虑开smearing，以及把截断能从400 Ry提高），确定把IGNORE_CONVERGENCE_FAILURE给去掉也能正常收敛了再说几何优化的问题。

PoorChaos · 发表于 Post on 3 day ago

Uus/pMeC6H4-/キ发表于 2026-4-20 14:35
先尝试算单点能并按http://sobereva.com/665解决SCF不收敛问题（此处可以考虑开smearing，以及把截断能从40 ...

谢谢！我先试试开smearing

PoorChaos · 发表于 Post on 3 day ago

Uus/pMeC6H4-/キ发表于 2026-4-20 14:35
先尝试算单点能并按http://sobereva.com/665解决SCF不收敛问题（此处可以考虑开smearing，以及把截断能从40 ...

谢谢
不开smearing SCF收敛没问题
我试试开smearing ，看看结构是否正常

northern_mar · 发表于 Post on 3 day ago

个人感觉还有几个小地方也许可以一起讨论一下。
1. x方向和y方向没设置k点，但是晶胞又不够大；
2. 不知道这个例子的收敛情况如何，自旋多重度不一定设置合理；
3. 如同上一位老师说的，含Ni的体系CUTOFF设置太低了。如果耗时太长可以考虑OT+全电子基组；
4. 除非特别大的体系，不然一开始一般不用LBFGS，而是优先使用BFGS和CG；

PoorChaos · 发表于 Post on 3 day ago

本帖最后由 PoorChaos 于 2026-4-20 23:13 编辑

northern_mar 发表于 2026-4-20 21:34
个人感觉还有几个小地方也许可以一起讨论一下。
1. x方向和y方向没设置k点，但是晶胞又不够大；
2. 不知 ...

谢谢答复
不加smearing，SCF容易收敛，结构也可以收敛，只是会散的比例厉害
我会尝试加K 点

如果加221 K点，全pop dzvp rev2  ，收敛不太好
是tep    Update method    Time Convergence       Total energy Change
  ------------------------------------------------------------------------------
   1 NoMix/Diag. 0.40E+00 12.3    4.15358412 -94736.7063089387 -9.47E+04
   2 Broy./Diag. 0.40E+00 13.9    1.32952109 -94619.6904807796  1.17E+02
   3 Broy./Diag. 0.40E+00 13.7 17.51187137 -94727.3492622163 -1.08E+02
   4 Broy./Diag. 0.40E+00 13.5 120.26863828 -94467.8404698225  2.60E+02
   5 Broy./Diag. 0.40E+00 13.5 162.11694664 -93633.4100162796  8.34E+02
   6 Broy./Diag. 0.40E+00 13.7 231.03799664 -92704.3994921492  9.29E+02
   7 Broy./Diag. 0.40E+00 13.6 231.06615429 -92242.9112017088  4.61E+02
   8 Broy./Diag. 0.40E+00 13.7 163.87664050 -90074.2296181857  2.17E+03
   9 Broy./Diag. 0.40E+00 13.7 122.28343960 -89881.3651945907  1.93E+02
10 Broy./Diag. 0.40E+00 13.4 175.46421865 -85390.7740296135  4.49E+03
11 Broy./Diag. 0.40E+00 13.6 227.75519169 -70859.0561701974  1.45E+04
12 Broy./Diag. 0.40E+00 13.7 227.75536503 -65671.5757085649  5.19E+03
13 Broy./Diag. 0.40E+00 13.8 217.70651379 -65908.6503469109 -2.37E+02

northern_mar · 发表于 Post on yesterday 09:53

PoorChaos 发表于 2026-4-20 22:59
谢谢答复
不加smearing，SCF容易收敛，结构也可以收敛，只是会散的比例厉害
我会尝试加K 点

如果不介意的话可以上传一下输入输出文件，大家一起看一下。
然后就是不知道你在结构优化前有没有做过cellrelax呢？

PoorChaos · 发表于 Post on yesterday 14:42

本帖最后由 PoorChaos 于 2026-4-22 15:01 编辑

northern_mar 发表于 2026-4-22 09:53
如果不介意的话可以上传一下输入输出文件，大家一起看一下。
然后就是不知道你在结构优化前有没有做过ce ...

没做过cell opt
已上传压缩包

阿夸夸 · 发表于 Post on yesterday 15:29

PoorChaos 发表于 2026-4-22 14:42
没做过cell opt
已上传压缩包

我跑了一下这个双层结构的变胞优化，优化了18步目前来看结构还正常，今天我看看能不能跑完，晚些我会把文件上传

PoorChaos · 发表于 Post on 4 hour ago

阿夸夸发表于 2026-4-22 15:29
我跑了一下这个双层结构的变胞优化，优化了18步目前来看结构还正常，今天我看看能不能跑完，晚些我会把文 ...

谢谢！
我跑完后是变形的
前面二十多步也很小变形

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

[CP2K] ZnNiLDH双层结构优化原子散开变形严重

浏览过的版块