CP2K 输入文件生成器

Skill 概述

本 Skill 用于根据用户的计算要求和结构文件,生成准确可用的 CP2K 输入文件。支持多种计算类型、结构格式和计算方法,并能够根据体系特征自动选择合适的参数。

触发词

当用户提及以下内容时应该使用此 Skill:

• - "生成 CP2K 输入文件"
• "创建 .inp 文件"
• "CP2K calculation input"
• "准备 CP2K 计算任务"
• "写一个 CP2K 输入文件"
• "帮我生成 CP2K 的输入"
• 包含计算类型关键词: "单点计算", "几何优化", "分子动力学", "MD", "频率分析", "NEB", "能带结构", "DOS"
• 提及结构文件格式: cif, pdb, xyz, gaussian
• 从 CP2K 输出文件生成新的输入

执行流程

步骤 1: 理解计算需求

向用户确认或从描述中提取以下信息:

必需参数:

1. 1. 计算类型 (RUN_TYPE):

- ENERGY / ENERGY_FORCE: 单点能量计算
- GEO_OPT: 几何优化
- MD: 分子动力学
- FREQUENCY: 频率分析
- VIBRATIONAL_ANALYSIS: 振动分析
- BAND_STRUCTURE: 能带结构
- NEB: 弹性带路径搜索

1. 2. 结构文件:

- 支持格式: .xyz, .cif, .pdb, .gjf, .com - 或提供直接的坐标数据

重要参数(如果没有提供需要询问):

1. 3. 计算方法:

- Quickstep (DFT): 密度泛函理论
- QMMM: 量子力学/分子力学混合方法
- Classical/MM: 纯分子力学
- xTB: 半经验方法

1. 4. 周期性条件:

- XYZ: 三维周期性(晶体) - XY: 二维周期性(表面/平板) - X/Y/Z: 一维周期性(纳米线) - NONE: 孤立体系(分子/团簇)

1. 5. 电荷和自旋:

- CHARGE: 净电荷 - MULTIPLICITY: 自旋多重度(如果自旋极化)

可选参数(如果没有提供可以使用默认值):

1. 6. 交换关联泛函: PBE, BLYP, PBE0, HSE06 等
2. 基组类型: DZVP-MOLOPT, TZVP-MOLOPT 等
3. k点采样: Monkhorst-Pack 网格
4. CUTOFF: 平面波截止能量
5. SCF 收敛判据
6. MD 参数: 系综,温度,时间步长,步数
7. 几何优化参数: 优化器,收敛判据,最大迭代次数

步骤 2: 读取和解析结构文件

使用 read_file 工具读取结构文件内容,然后:

• - CIF 文件: 提取晶胞参数和原子坐标
• PDB 文件: 提取原子类型和坐标,注意元素符号
• XYZ 文件: 提取原子类型和坐标
• Gaussian 文件: 提取原子类型和坐标
• CP2K 输出文件: 从输出中提取最后的结构坐标

步骤 3: 确定计算参数

根据体系特征和计算需求,智能选择参数:

基组选择:
CODEBLOCK0

CUTOFF 选择:
CODEBLOCK1

k点选择:
CODEBLOCK2

SCF 参数:
CODEBLOCK3

步骤 4: 生成输入文件内容

根据计算类型生成相应的输入文件结构:

通用结构框架:

CODEBLOCK4

不同计算类型的 MOTION 设置:

GEO_OPT (几何优化):
CODEBLOCK5

MD (分子动力学):
CODEBLOCK6

FREQUENCY (频率分析):
CODEBLOCK7

NEB (弹性带):
CODEBLOCK8

步骤 5: 处理特殊情况

自旋极化计算:

• - 添加 UKS 关键字
• 设置正确的 MULTIPLICITY 值
• 对于铁磁性系统,设置所有原子的初始磁矩

金属体系:

• - 使用 Fermi-Dirac smearing
• 设置 INLINECODE3
• 使用较密的 k 点采样

表面/平板模型:

• - 周期性设置为 XY
• 使用适当的真空层(通常 >15 Å)
• k 点设置为 8x8x1 或类似

QMMM 计算:

• - METHOD 设置为 QMMM
• 定义 QM 和 MM 区域
• 设置耦合参数

xTB 方法:

• - 在 &QS 中设置 INLINECODE5
• 添加 &xTB 部分
• 包含必要的参数文件

步骤 6: 使用辅助脚本(可选)

对于复杂的输入文件生成,可以使用 scripts/generate_cp2k_input.py 脚本:

CODEBLOCK9

支持的命令:

• - energy: 单点能量计算
• INLINECODE9: 几何优化
• INLINECODE10: 分子动力学

步骤 7: 输入文件验证

生成输入文件后,检查:

1. 1. 语法正确性:

- 所有 &SECTION 都有对应的 &END SECTION - 关键字拼写正确 - 参数在合理范围内

1. 2. 参数一致性:

- CHARGE 和 MULTIPLICITY 匹配 - CUTOFF 和基组匹配 - k 点和周期性匹配

1. 3. 文件路径:

- 结构文件路径正确 - 基组文件和赝势文件路径正确

1. 4. 资源预估:

- 估算计算所需的内存和 CPU 核心数 - 提醒用户计算可能需要的时间

步骤 8: 输出文件

将生成的输入文件写入 function.inp:

CODEBLOCK10

参考资源

内部参考文档

在 references/cp2k_input_reference.md 中包含:

• - 完整的 CP2K 输入文件结构说明
• 各个关键字的详细解释
• 不同计算类型的示例
• 参数选择指南
• 常见错误和解决方案

外部参考

当遇到复杂情况或特殊需求时,参考:

1. 1. CP2K 官方文档: https://www.cp2k.org/

- Input Reference: 详细的关键字说明 - Howtos: 特定计算类型的教程

1. 2. 用户工作区的示例文件:

- /file/SP/: 单点计算示例 - /file/opt/: 几何优化示例 - /file/MD/: 分子动力学示例 - /file/NEB/: NEB 计算示例 - /file/DOS/: 态密度示例 - /file/freq/: 频率分析示例 - /file/QMMM/: QMMM 计算示例

加载参考文档

在生成输入文件时,如果需要详细参考:

CODEBLOCK11

计算类型快速参考

计算类型	RUNTYPE	主要用途	关键参数
单点能量	ENERGY	计算体系总能量	CUTOFF, k-points
能量+力

ENERGYFORCE | MD前准备 | 同上 | | 几何优化 | GEO_OPT | 优化结构 | OPTIMIZER, 收敛判据 | | 分子动力学 | MD | 模拟时间演化 | ENSEMBLE, TIMESTEP, STEPS | | 频率分析 | FREQUENCY | 计算振动频率 | - | | 能带结构 | BAND_STRUCTURE | 固体能带 | k点路径 | | 态密度 | ENERGY + 特殊k点 | 电子态密度 | 高密度k点 | | NEB | NEB | 反应路径 | NIMAGE, K_SPRING | | QMMM | ENERGY/QMMM | 大体系精确计算 | QM/MM分区 |

常见问题处理

问题 1: SCF 不收敛
解决方案:

• - 降低 EPSSCF 到 1.0E-05
• 增加混合参数 ALPHA
• 使用不同的混合方法 (BROYDENMIXING → PULAY_MIXING)
• 增加网格精度 (提高 CUTOFF)

问题 2: 计算太慢
解决方案:

• - 降低基组精度 (TZVP → DZVP)
• 减少k点数
• 使用 xTB 方法
• 使用 OT 方法代替 DIAGONALIZATION

问题 3: 内存不足
解决方案:

• - 降低 CUTOFF
• 使用 DAVIDSON 对角化
• 减少k点数
• 使用 OT 方法

问题 4: 结果不准确
解决方案:

• - 提高 CUTOFF
• 使用更高精度的基组
• 增加k点采样
• 使用更精确的泛函 (PBE → PBE0 → HSE06)

注意事项

1. 1. 始终使用绝对路径 当引用文件时
2. 保留原始文件 不要直接修改用户的结构文件
3. 添加有意义的注释 在输入文件中解释关键参数的选择
4. 提供多种选择 当有多种可能方案时
5. 询问重要参数 当不确定时,不要随意设置
6. 参考已有示例 查看工作区中相似计算的输入文件

输出文件命名

默认输出文件名为 function.inp,但如果用户指定了项目名称,可以使用:

• - INLINECODE23
• INLINECODE24
• INLINECODE25

始终确保输出文件在用户指定的位置。