CP2K模拟计算

CP2K在大体系DFT中的定位

CP2K的核心优势在于其高斯平面波混合基组（Gaussian and Plane Wave, GPW）方法。传统平面波代码（如VASP、Quantum ESPRESSO）的计算量随体系增大呈三次方增长，而CP2K的GPW方法利用高斯函数的局域性，使得计算量更接近线性标度。在300原子以上的大体系中，CP2K的效率优势开始显现。

本项目中CP2K主要用于两类任务：大体系AIMD模拟和金属/氧化物界面的第一性原理计算。关于[DFT计算](https://www.keyanxueshu.com/category/dft/)的软件选型对比，站内有详细分析。

输入文件结构

CP2K的输入文件（通常命名为.inp或.input）采用嵌套的&关键字结构：

&GLOBAL
PROJECT_NAME water_cluster_aimd
RUN_TYPE MD
&END GLOBAL

&FORCE_EVAL
METHOD Quickstep
&DFT
…DFT设置…
&END DFT
&SUBSYS
…体系定义…
&END SUBSYS
&END FORCE_EVAL

&MOTION
…MD设置…
&END MOTION

三个顶层模块各司其职：GLOBAL定义项目信息，FORCE_EVAL定义如何计算力（DFT/MM/混合），MOTION定义如何运动（MD/优化/NEB）。

DFT模块关键参数

基组和赝势

CP2K使用双高斯基组（double-zeta）作为默认，精度选项通过BASIS_SET指定：

赝势选择：

&KIND H
BASIS_SET DZVP-MOLOPT-SR-GTH
POTENTIAL GTH-PBE-q1
&END KIND

&KIND O
BASIS_SET DZVP-MOLOPT-SR-GTH
POTENTIAL GTH-PBE-q6
&END KIND

关键提醒：基组和赝势必须匹配——PBE泛函用GTH-PBE赝势，BLYP泛函用GTH-BLYP赝势。混用会导致严重错误。CP2K的错误信息中常见”basis set not found”，通常是因为基组名称拼写错误或该元素在选定基组中不可用。

SCF收敛：OT优化器

CP2K的核心特色之一是Orbital Transformation (OT)方法，它避免了直接对角化哈密顿量，而是通过迭代优化轨道变换来达到SCF收敛。对于绝缘体/半导体体系，OT比传统DIIS方法快5-10倍。

&SCF
SCF_GUESS ATOMIC
EPS_SCF 1.0E-6
MAX_SCF 50

&OT
MINIMIZER DIIS
PRECONDITIONER FULL_SINGLE_INVERSE
&END OT

&OUTER_SCF
EPS_SCF 1.0E-6
MAX_SCF 10
&END OUTER_SCF
&END SCF

金属体系处理：OT默认不支持部分占据的能带，因此金属体系需要额外设置：

&SCF
&OT
MINIMIZER CG
PRECONDITIONER FULL_ALL
STEPSIZE 0.08
&END OT

&SMEAR
METHOD FERMI_DIRAC
ELECTRONIC_TEMPERATURE 300
&END SMEAR

&DIAGONALIZATION
ALGORITHM STANDARD
&END DIAGONALIZATION
&END SCF

金属体系建议关闭OT，改用传统DIIS+对角化方法，虽然慢但更稳定。

密度拟合（Auxiliary Density Matrix Method）

ADMM是CP2K的加速技术，用一个辅助基组拟合密度矩阵，减少交换相关泛函的计算量：

&AUXILIARY_DENSITY_MATRIX_METHOD
METHOD BASIS_PROJECTION
ADMM_PURIFICATION_METHOD NONE
&END AUXILIARY_DENSITY_MATRIX_METHOD

配合辅助基组：

&KIND O
BASIS_SET DZVP-MOLOPT-SR-GTH
POTENTIAL GTH-PBE-q6
BASIS_SET ADMM    # 使用默认ADMM辅助基组
&END KIND

ADMM对含杂化泛函（HSE06/PBE0）的计算加速效果最明显（3-5倍），对纯GGA也有1.5-2倍加速。

AIMD设置

CP2K做AIMD的核心在&MOTION模块：

&MD
ENSEMBLE NVT
STEPS 10000
TIMESTEP 0.5
TEMPERATURE 300.0

&THERMOSTAT
TYPE NOSE
NOSE_LENGTH 4
NOSE_YOSHIDA 3
NOSE_TIME 100.0
&END THERMOSTAT
&END MD

时间步长的选择：含氢原子的体系用0.5 fs是标准推荐。如果所有氢原子都替换成氘（同位素替代），可以用1.0 fs。本项目曾用1.0 fs跑含氢体系，5 ps后能量漂移超过0.01 Hartree——这是不可接受的。降到0.5 fs后，能量漂移<0.001 Hartree/ps。

QM/MM耦合

CP2K内置了QM/MM功能，可以对体系的一部分做DFT计算，其余部分用力场：

&FORCE_EVAL
METHOD QMMM

&QM_KIND O
…DFT参数…
&END QM_KIND

&MM_KIND C
…力场参数…
&END MM_KIND

&QMMM
ECOUPL NONE    # 或用ELECTROSTATIC
QM_CELL 30 30 30
&END QMMM
&END FORCE_EVAL

QM/MM适合处理溶液中的化学反应——溶质和第一溶剂壳层用QM，远端溶剂用MM。本项目在水溶液中有机反应的QM/MM模拟中，QM区域约50-100原子（反应中心+6-8个水分子），MM区域约500-2000原子（其余溶剂），总计算量比全QM减少5-10倍。

性能对比：CP2K vs VASP

以300原子的TiO₂-水界面AIMD为例：

CP2K在这个体系中快了4倍，精度偏差在可接受范围内。但需要注意，CP2K的高斯基组在描述真空区域时不如平面波准确——对于表面/分子体系，真空层中的密度尾部可能出现非物理行为。关于[科研学术网](https://www.keyanxueshu.com/)中CP2K的更多应用案例，站内有系列教程。

常见报错排查

CP2K在中小体系的精度略逊于VASP，但在大体系和AIMD场景中的效率优势是压倒性的。随着机器学习势函数（MLP）的兴起，”CP2K生成训练数据→训练MLP→LAMMPS大规模模拟”的流程正在成为新的标准范式，CP2K作为数据生成工具的角色越来越重要。