请问老师们，在CP2K计算时，预优化没问题，但进行K点测试不收敛该如何处理

方不惑 · 发表于 Post on 2 hour ago

本帖最后由方不惑于 2026-8-25 18:33 编辑

老师您好，我在计算纯金属体系的时候，元素是Zn，盒子大小我选用的是10.4*9.7*35，使用PBE-D3-DZVP-MOLOPT-SR-GTH基组进行预优化，预优化结束之后进行K点测试，先对XY进行测试，一开始测到10 10 1，K点还是不收敛，于是我这边开始继续加大K点的数据，期间我也提升过截断能，从350提升至700，测试一直到15 15 1才收敛，统一体系下的对比组测到17 17 1也不收敛，感觉这个K点太离谱了，但不知缘由，还望老师赐教，下面是我测试的结果
# k-points |  Energy (Hartree)  | delte E
   1 1 1 -2905.2681981149
   2 2 1 -2905.2988394794 -0.0306413645
   3 3 1 -2905.3103995685 -0.0115600891
   4 4 1 -2905.3146813045 -0.0042817360
   5 5 1 -2905.3151183294 -0.0004370249
   6 6 1 -2905.3139504906    0.0011678388
   7 7 1 -2905.3106934225    0.0032570682
   8 8 1 -2905.3127379852 -0.0020445627
   9 9 1 -2905.3133454338 -0.0006074486
   10 10 1 -2905.3120685642    0.0012768696
   11 11 1 -2905.3127139506 -0.0006453864
   12 12 1 -2905.3127862344 -0.0000722838
   13 13 1 -2905.3124741012    0.0003121332
   14 14 1 -2905.3126615880 -0.0001874868
   15 15 1 -2905.3126453426    0.0000162454
   16 16 1 -2905.3125202476    0.0001250949
   17 17 1 -2905.3125672746 -0.0000470270
   18 18 1 -2905.3121220181    0.0004452565
   19 19 1 -2905.3120740500    0.0000479681

UW_0728. · 发表于 Post on halfhour ago

原子数是多少（看你盒子较大原子数应该不少吧），如果较多（多于10个）最好按1meV/atom的标准来选

方不惑 · 发表于 Post on 10 min ago

UW_0728. 发表于 2026-8-25 19:05
原子数是多少（看你盒子较大原子数应该不少吧），如果较多（多于10个）最好按1meV/atom的标准来选

好的，谢谢老师，老师这个是我的输入文件，大概有40多个原子，其实我还有一个问题想请教老师，咱们算吸附能的时候，做受体表面优化的时候，需要在表面上建立真空层做定胞优化吗，还是直接一开始就变胞优化，变胞优化之后扩展真空层加入吸附分子呢？有些疑惑，我目前这个体系这两种方法都在尝试，不知道哪一种是正确的手法
&FORCE_EVAL
  METHOD Quickstep
  &SUBSYS
&CELL
   A    9.05640000    0.00000000    0.00000000
   B    0.00000000 10.45740000    0.00000000
   C    0.00000000    0.00000000 35.87320000
   PERIODIC XYZ #Direction(s) of applied PBC (geometry aspect)
&END CELL
&COORD
   Zn       6.06967615 10.45740000 0.98827352
   Zn       8.33639444 1.31196221 0.98428187
   Zn       1.54147615 10.45740000 0.98827348
   Zn       3.80819443 1.31196217 0.98428188
   Zn       6.07094652 2.60612049 0.98753335
   Zn       8.33441028 3.92617095 0.98854463
   Zn       1.54274654 2.60612051 0.98753338
   Zn       3.80621025 3.92617097 0.98854468
   Zn       7.52069081 10.45739999 3.40872151
   Zn       5.25198173 1.31395241 3.41061360
   Zn       6.06979716 10.45740000 5.89912355
   Zn       8.33623289 1.31376814 5.90133304
   Zn       2.99249084 10.45739999 3.40872161
   Zn       0.72378172 1.31395239 3.41061361
   Zn       1.54159720 10.45740001 5.89912342
   Zn       3.80803288 1.31376814 5.90133301
   Zn       7.51912332 2.60439421 3.40884201
   Zn       5.25621743 3.92739346 3.40854780
   Zn       6.07093931 2.60506477 5.89954345
   Zn       8.33455255 3.92694567 5.89909910
   Zn       2.99092325 2.60439432 3.40884205
   Zn       0.72801753 3.92739347 3.40854777
   Zn       1.54273936 2.60506485 5.89954345
   Zn       3.80635249 3.92694571 5.89909926
   Zn       6.07332375 5.22869997 0.98212628
   Zn       8.33441026 6.53122905 0.98854463
   Zn       1.54512367 5.22870000 0.98212633
   Zn       3.80621027 6.53122905 0.98854467
   Zn       6.07094654 7.85127948 0.98753336
   Zn       8.33639445 9.14543787 0.98428184
   Zn       1.54274653 7.85127950 0.98753332
   Zn       3.80819444 9.14543786 0.98428186
   Zn       7.51448420 5.22869989 3.41248618
   Zn       5.25621754 6.53000652 3.40854780
   Zn       6.07293991 5.22869993 5.90337970
   Zn       8.33455251 6.53045432 5.89909915
   Zn       2.98628415 5.22870019 3.41248617
   Zn       0.72801747 6.53000656 3.40854778
   Zn       1.54473976 5.22870008 5.90337975
   Zn       3.80635253 6.53045431 5.89909923
   Zn       7.51912330 7.85300571 3.40884203
   Zn       5.25198176 9.14344759 3.41061360
   Zn       6.07093934 7.85233517 5.89954345
   Zn       8.33623288 9.14363190 5.90133304
   Zn       2.99092335 7.85300583 3.40884205
   Zn       0.72378174 9.14344756 3.41061360
   Zn       1.54273932 7.85233521 5.89954345
   Zn       3.80803290 9.14363190 5.90133302
&END COORD
&KIND Zn
   ELEMENT Zn
   BASIS_SET DZVP-MOLOPT-SR-GTH-q12
   POTENTIAL GTH-PBE
&END KIND
  &END SUBSYS

  &DFT
BASIS_SET_FILE_NAME  BASIS_MOLOPT
POTENTIAL_FILE_NAME  POTENTIAL
# WFN_RESTART_FILE_NAME kpconv-RESTART.wfn
CHARGE 0 #Net charge
MULTIPLICITY 1 #Spin multiplicity
&KPOINTS
   SCHEME MONKHORST-PACK kp_test
&END KPOINTS
&QS
   EPS_DEFAULT 1.0E-11 #Set all EPS_xxx to values such that the energy will be correct up to this value
&END QS
&POISSON
   PERIODIC XYZ #Direction(s) of PBC for calculating electrostatics
   PSOLVER PERIODIC #The way to solve Poisson equation
&END POISSON

UW_0728. · 发表于 Post on 59 sec ago

方不惑发表于 2026-8-25 19:37
好的，谢谢老师，老师这个是我的输入文件，大概有40多个原子，其实我还有一个问题想请教老师，咱们算吸附 ...

印象中常用的整个流程是：先对体相做变胞优化，优化完后再切表面、加真空层，固定表面结构最下面2-3层原子，做二维周期性的定胞优化，最后在优化完的表面结构中再加入吸附分子去做定胞优化；算吸附能要用到的能量数据包括上面第二步和第三步对应结构的能量

关于你题目问的问题，40多个原子的话，5 5 1的k点就够了

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 请问老师们，在CP2K计算时，预优化没问题，但进行K点测试不收敛该如何处理