计算化学公社

标题: ZnNiLDH双层结构优化原子散开变形严重 [打印本页]

作者
Author: PoorChaos 时间: 10 hour ago
标题: ZnNiLDH双层结构优化原子散开变形严重
本帖最后由 PoorChaos 于 2026-4-20 13:24 编辑

请教 ZnNiLDH的结构优化问题
之前是中间插层放入有机物难优化成功
现在优化未吸附有机物的基底结构，原子变化散开
(, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0)
inp文件见附件 (, 下载次数 Times of downloads: 4)

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

作者
Author: PoorChaos 时间: 8 hour ago

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

谢谢！我先试试开smearing

作者
Author: PoorChaos 时间: 8 hour ago

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

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

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

作者
Author: PoorChaos 时间: halfhour 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

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