1.2 选择并审计赝势
这些输入和起始模型为未执行的教学示例。原创配图为示意图,不是计算结果。使用前请验证版本相关语法、授权数据、数值收敛及科学模型。
1.2.1 模型、单位与来源记录
PW 截断能与能量使用 Ry,常见力输出使用 Ry/bohr,压力使用 kbar。几何卡片须明确坐标单位。不同程序具有不同输入语法与时间单位约定。
1.2.2 未执行输入与明确修改片段
按照下方说明确认每个片段的母计算与替换位置;片段不自动构成独立输入。保留空行与文件来源要求。
1.2.2.1 输入块 1
sha256sum pseudo/*.UPF > provenance/pseudopotentials.sha256
# Inspect the readable UPF header; filenames below are local aliases.
head -n 60 pseudo/Si.PBE.UPF
head -n 60 pseudo/H.PBE.UPF
1.2.3 实作研究
1.2.3.1 原理与准备
赝势属于哈密顿量的一部分,并非输入中的装饰附件。价电子组态、交换关联泛函、相对论处理与构建类型都会影响计算。本例为后续练习准备 Si、H,并可加入 Al。请从经过整理的库或原作者获取文件,保留引用及许可信息,得到实际 UPF 文件后再运行模板。
1.2.3.2 原创审计步骤
参见上方输入块 1。
每个文件建立一行记录:元素、z_valence、泛函、模守恒/超软/PAW 类型、标量或全相对论、来源网址与 SHA-256。读取文件确认,不要凭文件名猜测。推荐波函数与电荷密度截断能只能作为测试起点。电荷密度截断能等于波函数的四倍并非普遍规律,尤其不能套用于所有超软或 PAW 赝势。改变赝势后,前后绝对总能不能直接比较。
1.2.3.3 操作与输出
普通案例选用相互一致的 PBE 赝势。可在课程目录使用清晰的别名,但清单必须保留原文件名。读取 1.4 输出电子数,与价电子总数核对。SOC 案例需选用经过验证的全相对论文件,并说明它与标量相对论参考为何可比。特殊响应程序的赝势兼容性要另行检查;基态能运行不代表响应计算一定支持。
1.2.3.4 检查、陷阱与练习
练习:选择两个合规 Si 赝势,分别收敛截断能,再比较平衡晶格常数和各自内部的相对能量;不能比较不同哈密顿量的原始绝对能。更低的推荐截断能不等于更高质量,它只体现部分计算成本。不要把 PBE 与 LDA 文件混在一起,再用 input_dft 强制覆盖来“修复”。