计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First-principles) › 针对熔盐体系测试能量随截断能收敛性时，无法收敛怎么办 ...

针对熔盐体系测试能量随截断能收敛性时，无法收敛怎么办？

查看数: 154 | 评论数: 5 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

以玉名诗

发布时间: 2026-1-9 23:29

正文摘要:

本帖最后由以玉名诗于 2026-1-12 22:38 编辑学生使用Packmol进行了周期性体系的构建，由于我不知道掺杂了EuF3之后体系的密度，故将盒子尺寸稍增大了些以构建周期性体系。从VMD来该体系应当是合理的，如图1.现 ...

回复 Reply

以玉名诗 发表于 Post on 2026-1-12 22:34:51

UW_0728. 发表于 2026-1-12 18:23
培训脚本是按照MAX_SCF 1的情况设置的（毕竟CP2K的截断能收敛性测试不是非得要求SCF收敛），这种情况Total ...

谢谢您老师，我修改了脚本现在尝试一下！

UW_0728. 发表于 Post on 2026-1-12 18:23:40

培训脚本是按照MAX_SCF 1的情况设置的（毕竟CP2K的截断能收敛性测试不是非得要求SCF收敛），这种情况Total Energy只会输出一次；你使用了OT且经过了不止一轮Outer SCF使其顺利收敛，导致Total Energy输出了多次，所以CUTOFF.txt的输出是异常的。在提取能量相应指令那行加个“| tail -n 1”即可。

以玉名诗 发表于 Post on 2026-1-12 14:50:31

sobereva 发表于 2026-1-10 01:49
输入文件里有MAX_SCF 1，当然会提示未收敛。培训里明确说了为什么这样设
测试脚本是北京科音CP2K第一性原 ...

Sob老师，我之前的体系测试的很顺利，更换了一个新体系后不知为什么脚本输出的ΔE这么奇怪，手动计算后发现450 Ry时能量变化首次低于了1E-6 Hartree，但按照常理来说含Na的体系截断能应该很高才对，请问老师是不是我的输入文件中存在问题？

cutconv.inp

10.05 KB, 下载次数 Times of downloads: 0

输入文件

test 900 Ry.out

161.91 KB, 下载次数 Times of downloads: 0

900 Ry截断能时输出文件

CUTOFF (1).txt

2.23 KB, 下载次数 Times of downloads: 0

截断能txt文件

以玉名诗 发表于 Post on 2026-1-10 13:33:55

sobereva 发表于 2026-1-10 01:49
输入文件里有MAX_SCF 1，当然会提示未收敛。培训里明确说了为什么这样设
测试脚本是北京科音CP2K第一性原 ...

谢谢老师，我还是得仔细复习才行。好的老师，我先前不了解这个不可以外传，下次不会上传了

sobereva 发表于 Post on 2026-1-10 01:49:49

输入文件里有MAX_SCF 1，当然会提示未收敛。培训里明确说了为什么这样设
测试脚本是北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）我专门给学员提供的，不要随便外传。我已把你的脚本附件删掉

评分 Rate

参与人数 Participants 1	eV +1	收起理由 Reason
Hurono	+ 1	谢谢

查看全部评分 View all ratings

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