计算化学公社

标题: 求助：CP2K晶胞优化报错 [打印本页]
作者
Author: wangyueda    时间: 2022-6-22 18:30
标题: 求助：CP2K晶胞优化报错
本帖最后由 wangyueda 于 2022-6-22 18:41 编辑

请问下各位老师，我用CP2K优化Cu单胞时报错，报错信息及输入文件如下：

cp2k.out:
  1. *******************************************************************************
  2. *   ___                                                                       *
  3. *  /   \                                                                      *
  4. * [ABORT]                                                                     *
  5. *  \___/    found unexpected extra argument BETA at , File: 'cp2k.inp', Line: *
  6. *    |                         44, Column: 2, Chunk: <0.1>                    *
  7. *  O/|                                                                        *
  8. * /| |                                                                        *
  9. * / \                                               input/input_parsing.F:629 *
  10. *******************************************************************************


  13. ===== Routine Calling Stack =====

  15.             7 val_create_parsing
  16.             6 section_vals_parse
  17.             5 section_vals_parse
  18.             4 section_vals_parse
  19.             3 section_vals_parse
  20.             2 section_vals_parse
  21.             1 read_input
复制代码



cp2k.inp:

  2. DATAPATH /home/gengzi/install-packages/cp2k-9.1/data

  4. #The cp2k inp is consist of SECTIONs and KEYWORDs
  5. #The SECTION is green (pink) with '&' ahead
  6. #The KEYWORD is yellow without '&'
  7. #The red and black words the parameters matching with the KEYWORD
  8. #The writing order : SECTION => KEYSORDS => SUBSECTION

  10. #######################  1. GLOBAL SECTION   ########################
  11. #This section defines the calculation type
  12. &GLOBAL
  13.    PROJECT cp2k
  14.    #  RUN_TYPE GEO_OPT #the lattice constant unchanged
  15.    #  RUN_TYPE ENERGY  #single point calculation
  16.    #  RUN_TYPE MD      # MD
  17.    RUN_TYPE CELL_OPT # the lattice constant changed
  18.    PRINT_LEVEL LOW
  19. &END GLOBAL
  20. #######################  END GLOBAL SECTION  ########################

  22. #######################  2. FORCE_EVAL SECTION  #####################
  23. #This section defines the contents for electron step
  24. &FORCE_EVAL
  25.    STRESS_TENSOR ANALYTICAL
  26.    METHOD Quickstep

  28.    ###################   2.1 DFT SECTION  ###################
  29.    &DFT
  30.       BASIS_SET_FILE_NAME ${DATAPATH}/BASIS_MOLOPT
  31.       POTENTIAL_FILE_NAME ${DATAPATH}/GTH_POTENTIALS
  32.       WFN_RESTART_FILE_NAME ./cp2k-RESTART.wfn
  33.       UKS T

  35.       ###############    2.1.1 QS SECTION  #################
  36.       &QS
  37.          EPS_DEFAULT 1.0E-12
  38.          WF_INTERPOLATION ASPC
  39.          EXTRAPOLATION_ORDER 3
  40.       &END QS

  42.       ###############    2.1.2 MGRID SECTION  #############
  43.       &MGRID
  44.          CUTOFF 400 #unit Rydberg
  45.          NGRIDS 4
  46.          REL_CUTOFF 60
  47.       &END MGRID

  49.       ###############    2.1.3 XC SECTION  ###############
  50.       &XC
  51.          &XC_FUNCTIONAL PBE
  52.          &END XC_FUNCTIONAL
  53.       &END XC

  55.       #################  2.1.4 SCF SECTION  ##############
  56.       &SCF
  57.          SCF_GUESS RESTART
  58.          EPS_SCF 1.0E-6
  59.          MAX_SCF 500
  60.          ADDED_MOS 500
  61.          CHOLESKY INVERSE

  63.          &DIAGONALIZATION
  64.             ALGORITHM STANDARD
  65.             EPS_ADAPT 0.01
  66.          &END DIAGONALIZATION
  67.                                     
  68.          &SMEAR ON
  69.             METHOD FERMI_DIRAC
  70.             ELECTRONIC_TEMPERATURE [K] 300
  71.          &END SMEAR

  73.          &MIXING
  74.             METHOD BROYDEN_MIXING
  75.             ALPHA 0.1 BETA 1.5
  76.             NBROYDEN 8
  77.          &END MIXING

  79.          &PRINT
  80.             &RESTART OFF
  81.             &END RESTART
  82.          &END PRINT
  83.       &END SCF

  85.       &KPOINTS
  86.          SCHEME MONKHORST-PACK 3 3 3
  87.       &END KPOINTS
  88.    &END DFT   
  89.    ###################   END DFT SECTION  ###################

  91.    ###################   2.2 SUBSYS SECTION  ################
  92.    #This section mainly contains the CELL ,TOPOLOGY (or COORD) and KIND sections

  94.    #the CELL section contains the basic cell info: 1. lattice constant 2. the angle between lattice
  95.    #the TOPOLOGY (COORD) section specified the structral models. COORD for calculation that need to fixed atoms
  96.    #the KIND section point out the specified the BASIS and POTENTIAL for each elements

  98.    &SUBSYS
  99.       #################  2.2.1 CELL SECTION  ##############
  100.       &CELL
  101.          ABC 3.6147 3.6147 3.6147
  102.          ALPHA_BETA_GAMMA 90 90 90
  103.       &END CELL

  105.       #################  2.2.2 TOPOLOGY(COORD) SECTION  ##############
  106.       #usage: no atoms need to be fixed
  107.        &TOPOLOGY
  108.          COORD_FILE_NAME Cu.cif
  109.          COORD_FILE_FORMAT CIF
  110.        &END TOPOLOGY
  111.       
  112.       ##usage: some atoms need to be fixed
  113.       ##the fixed atoms list is writed in MOTION=>CONSTRAINT

  115.       #################  2.2.3 KIND SECTION  ##################
  116.       &KIND C
  117.          ELEMENT C     
  118.          BASIS_SET DZVP-MOLOPT-SR-GTH
  119.          POTENTIAL GTH-PBE-q4
  120.       &END KIND

  122.       &KIND H
  123.          ELEMENT H
  124.          BASIS_SET DZVP-MOLOPT-SR-GTH
  125.          POTENTIAL GTH-PBE-q1
  126.          MASS 2.014
  127.       &END KIND

  129.       &KIND O
  130.          ELEMENT O
  131.          BASIS_SET DZVP-MOLOPT-SR-GTH
  132.          POTENTIAL GTH-PBE-q6
  133.       &END KIND

  135.       &KIND F
  136.          ELEMENT F
  137.          BASIS_SET DZVP-MOLOPT-SR-GTH
  138.          POTENTIAL GTH-PBE-q7
  139.       &END KIND

  141.       &KIND Cu
  142.          ELEMENT Cu
  143.          BASIS_SET DZVP-MOLOPT-SR-GTH
  144.          POTENTIAL GTH-PBE-q11
  145.       &END KIND

  147.       &KIND Li
  148.          ELEMENT Li
  149.          BASIS_SET DZVP-MOLOPT-SR-GTH
  150.          POTENTIAL GTH-PBE-q3
  151.       &END KIND

  153.       &KIND Na
  154.          ELEMENT Na
  155.          BASIS_SET DZVP-MOLOPT-SR-GTH
  156.          POTENTIAL GTH-PBE-q9
  157.       &END KIND

  159.    &END SUBSYS
  160.    ###################   END SUBSYS SECTION  ################
  161. &END FORCE_EVAL
  162. ######################   END FORCE_EVAL SECTION #####################

  164. #######################  3. MOTION SECTION  #########################
  165. #This section defines the contents for ion step
  166. &MOTION

  168.    &CELL_OPT
  169.       EXTERNAL_PRESSURE 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0
  170.       KEEP_ANGLES T
  171.       KEEP_SYMMETRY T
  172.       OPTIMIZER CG
  173.       &CG
  174.          &LINE_SEARCH
  175.             TYPE 2PNT
  176.          &END LINE_SEARCH
  177.       &END CG
  178.   &END CELL_OPT

  180.    &GEO_OPT
  181.       MAX_ITER 400
  182.       OPTIMIZER BFGS
  183.       MAX_FORCE 6.0E-4
  184.    &END GEO_OPT

  186. &END MOTION
  187. #######################  END MOTION SECTION  #########################
复制代码



作者
Author: sobereva    时间: 2022-6-22 20:37
BETA 1.5甭和前面的写到同一行

不熟悉的话就直接用Multiwfn产生输入文件，又省事格式又肯定没错误
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587（http://bbs.keinsci.com/thread-21668-1-1.html）
作者
Author: wangyueda    时间: 2022-6-22 20:47
sobereva 发表于 2022-6-22 20:37
BETA 1.5甭和前面的写到同一行

不熟悉的话就直接用Multiwfn产生输入文件，又省事格式又肯定没错误

好的卢老师 这就用起来



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3