计算化学公社

标题: 针对熔盐体系测试能量随截断能收敛性时，无法收敛怎么办？ [打印本页]

作者
Author: 以玉名诗 时间: 2026-1-9 23:29
标题: 针对熔盐体系测试能量随截断能收敛性时，无法收敛怎么办？
本帖最后由以玉名诗于 2026-1-12 22:38 编辑

学生使用Packmol进行了周期性体系的构建，由于我不知道掺杂了EuF3之后体系的密度，故将盒子尺寸稍增大了些以构建周期性体系。从VMD来该体系应当是合理的，如图1.现在使用Sob老师写的进行截断能测试的脚本，设置参考截断能为60 Ry进行从100 Ry到900 Ry的计算。并在最前面加入了几行新命令：
odule purge
module load compiler/gnu/13.2
module load mpi/openmpi/4.1.1

export OMPI_MCA_btl="^openib"
export OMPI_MCA_pml="ucx"
export UCX_TLS="all"

之后提交至超算平台进行计算。输出文件的结果中显示所有截断能测试点均未收敛，且能量非常之大，学生调整了多次参数无果，请教老师提出宝贵的建议。

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

作者
Author: 以玉名诗 时间: 2026-1-10 13:33

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

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

作者
Author: 以玉名诗 时间: 2026-1-12 14:50

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

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

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

作者
Author: 以玉名诗 时间: 2026-1-12 22:34

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

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

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