CP2K学习笔记

1、在写某一节时，最小节直接用END，而不是用END + 节的名字

2、对于要模拟摩擦时，要控制原子的移动，可以用集合变量DISTANCE来控制，然后在约束部分MOTIOM /CONSTRAINT控制距离的增长，参考代码如下

    &COLVAR
     &DISTANCE
      &POINT
       ATOMS 19..36 55..72  ###注意看是否约束的原子是否多约束或者少约束
       TYPE GEO_CENTER
      &END
      &POINT
       XYZ 3.7  0 13    ###建议对应施加约束的轴设为0，然后在TARGET 位置出写出约束原子的几何中心，比如目前的是y为6.05
       TYPE FIX_POINT
      &END
      POINTS 1 2
      AXIS Y
     &END DISTANCE
    &END COLVAR
####上面是定义DISTANCE集合变量的部分，定义了上面层石墨烯的几何中心为参考位置，以及定义了一个空间点为另一个参考位置

 &CONSTRAINT
  &COLLECTIVE
   &RESTRAINT
   K 0.01
   &END
   INTERMOLECULAR
   COLVAR 1
   TARGET [angstrom] 6.05  ###指定长度单位
   TARGET_GROWTH 0.05
  &END COLLECTIVE
&END CONSTRAINT
###上面定义了实验使用的约束变量是集合变量1  目标距离是6.05(空间点与几何中心的距离) 增长速度是0.05，用弹簧势约束，劲度系数为0.01

3、关于单位的使用

在数值前面用方括号[]指定单位，如果手册中写了[internal_cp2k]对应于单位网页中对应的第一个单位或者衍生单位