任务: VASP Slab 计算的 VASP 参数查询
- 模块: computational-modeling
- 深度: L3
- 学科: 化学（表面科学 / 电催化）
- 判定: 通过
- 核心声明: 是
任务描述
查询含 Pt 单原子的 CeO2(111) slab 计算的 VASP ISIF、ENCUT 和 ISMEAR 参数。AI 必须遵循路由协议：先读 INDEX.md，再读 query_router.md，然后执行带 5 层证据标签的融合查询。
输入
User query: "VASP slab 计算中 ISIF 怎么设置？我的体系是 CeO2(111) slab，含 Pt 单原子"
System: 6-layer slab, bottom 3 fixed, top 3 + Pt relaxed
Task Output: computational-modeling/cm-01
Five-Layer Evidence Report
【当前材料】
用户体系：CeO2(111) slab，6 层，底部 3 层固定（selective dynamics: F F F），顶部 3 层 + Pt 单原子弛豫（T T T）。 关键约束：
- slab 计算需要保持底层晶格固定（模拟体相约束）
- 含稀土元素 Ce（f 电子），需要足够高的 ENCUT
- 金属单原子 Pt 可能引入导电性 → 影响 ISMEAR 选择
【论文知识】
未安装 paper-reader-workflow，跳过论文知识层。
【软件知识】
ISIF 参数 (来源：VASP 官方 Wiki, v6.x)
| ISIF 值 | 含义 | 对 slab 的影响 |
|---|---|---|
| 0 | 仅计算能量/力，不更新位置也不改变晶胞 | 单点能计算 |
| 1 | 仅计算能量/力，不更新位置 | 同 0（VASP 6.x 中 0 和 1 行为相同） |
| 2 | 弛豫离子位置，晶胞形状和体积固定 | slab 优化推荐：胞参数保持不变 |
| 3 | 弛豫离子位置 + 晶胞形状和体积 | 不推荐 slab：会改变真空层高度和面内尺寸 |
| 4 | 弛豫离子位置 + 晶胞形状 | 不推荐 |
| 5 | 弛豫离子位置 + 晶胞形状 | 同 4 |
| 6 | 弛豫晶胞形状和体积，离子位置固定 | 不推荐 |
| 7 | 弛豫晶胞体积，离子位置和晶胞形状固定 | 不推荐 |
ISIF 与 selective dynamics 的关系：
- ISIF ≤ 2 时，POSCAR 中的 selective dynamics 标记（F/T）仍然生效
- ISIF = 2 配合 selective dynamics 是 slab 优化的标准做法
ENCUT 参数 (来源：VASP 官方 Wiki, v6.x)
- ENCUT 控制平面波截断能
- 对于含 Ce 的体系（Ce POTCAR 含 f 电子），默认 ENCUT（从 POTCAR 读取）通常需要乘以 1.3-1.5
- Ce 的 POTCAR (Ce_3, Ce_sv 等) 的 ENMAX 通常为 ~275-400 eV
- 推荐 ENCUT = max(ENMAX of all elements) × 1.3 = ~520 eV（保守）
ISMEAR 与 SIGMA (来源：VASP 官方 Wiki)
| 体系类型 | ISMEAR | SIGMA (eV) |
|---|---|---|
| 绝缘体/半导体 | 0 (Gaussian) | 0.05-0.10 |
| 金属 | 1 (Methfessel-Paxton) | 0.10-0.20 |
| 大体系 (>100 atoms) | -1 (Fermi) / 0 | 0.05 |
| 所有体系（准确总能量） | -5 (tetrahedron with Blöchl) | — |
对于 CeO2 slab + Pt（可能引入 gap states → 部分导电），推荐 ISMEAR = 1, SIGMA = 0.10 eV；最终能量计算用 ISMEAR = -5。
【本地经验】
当前未检索到该体系本地验证经验。
【AI 判断】
ISIF 推荐：ISIF = 2
理由：
- slab 计算需要保持面内晶格常数（匹配体相）和真空层高度
- ISIF = 3 会改变晶胞，对 slab 模型有风险（真空层可能被压缩）
- 配合 selective dynamics（底部 3 层 F F F）实现部分弛豫
风险：
- 如果晶格常数来自实验值而非优化后的体相值，可能需要先做 bulk 优化获取准确晶格常数
- 建议先用 ISIF = 3 对 bulk CeO2 做优化，再用优化后的晶格常数构建 slab
ENCUT 推荐：ENCUT = 520 eV
理由：
- 含 Ce 体系（f 电子）需要高于默认值的 ENCUT
- 520 eV = ENMAX(Ce) × 1.3，是含 f 电子元素的标准推荐
- 进一步收敛性测试：ENCUT = 400, 450, 500, 550, 600 eV → 确认 520 eV 已收敛
ISMEAR 推荐：ISMEAR = 1, SIGMA = 0.10
理由：
- CeO2 是绝缘体（带隙 ~3.2 eV），但 Pt 单原子可能引入 mid-gap states
- ISMEAR = 1 对金属/半金属体系更稳定
- 最终能量应使用 ISMEAR = -5 重新计算
创建的记录
Record 1: vasp-param-ISIF-0001
record_id: vasp-param-ISIF-0001
record_type: parameter
scope: single_software
software:
- VASP
topic: ISIF
source_type: official_doc
source_ref: https://www.vasp.at/wiki/index.php/ISIF
version_or_doc_date: VASP 6.x / 2026-06-15
claim: ISIF controls whether the stress tensor is calculated and which degrees of freedom (ionic positions, cell shape, cell volume) are allowed to change during relaxation.
evidence_level: must
confidence: official
related_wiki:
- computer_kb/wiki/concepts/ISIF.md
related_records:
- vasp-param-ENCUT-0001
- vasp-param-ISMEAR-0001
note: For slab calculations, ISIF=2 is recommended to keep cell dimensions fixed. ISIF=3 may compress the vacuum layer. Use selective dynamics (F/T flags in POSCAR) to fix bottom layers.
last_checked: 2026-06-15
Record 2: vasp-param-ENCUT-0001
record_id: vasp-param-ENCUT-0001
record_type: parameter
scope: single_software
software:
- VASP
topic: ENCUT
source_type: official_doc
source_ref: https://www.vasp.at/wiki/index.php/ENCUT
version_or_doc_date: VASP 6.x / 2026-06-15
claim: ENCUT specifies the cutoff energy for the plane-wave basis set in eV. Must be ≥ ENMAX of all POTCAR elements. For f-electron elements (Ce, U, etc.), multiply ENMAX × 1.3.
evidence_level: must
confidence: official
related_wiki:
- computer_kb/wiki/concepts/ENCUT.md
related_records:
- vasp-param-ISIF-0001
note: For Ce-containing systems, recommended ENCUT ≥ 520 eV. Always run convergence test from 1.0×ENMAX to 1.5×ENMAX in steps of 50-100 eV.
last_checked: 2026-06-15