other questions about the calculation error

chun chun... at gmail.com
Sun Oct 7 21:43:30 UTC 2007


Teo,
    Thank you very much for your reply last time. I have other more
questions. Hopfully, you can enlighten me.
    I used cp2k to calculate one water molecule at the center of a
unit cell without periodic boundary condition (PBC NONE). If I change
the size of the unit cell (ABC) from 9.4 Angstrom cubic to 5.4
Angstrom cubic, the two results are different. After I checked the
output files, I found cp2k used two different sets of ' max spherical
cutoff ' in PW_GRID.

First of all, what is the ' max spherical cutoff ' ? Why is there
spherical cutoff in PW_GRID (planewave grid?) ?

Secondly, why is the ' max spherical cutoff ' changed for the
different size of the unit cell ?

Thirdly, which one is more reliable? the small unit cell or the large
one?


Thanks

Chun-Min



&FORCE_EVAL
  METHOD Quickstep
  &DFT
    BASIS_SET_FILE_NAME ~/cp2k/tests/QS/ALL_BASIS_SETS
    POTENTIAL_FILE_NAME ~/cp2k/tests/QS/POTENTIAL
    &MGRID
      CUTOFF 158
    &END MGRID
    &QS
      METHOD GAPW
      FULL_GAPW
    &END QS
    &SCF
      SCF_GUESS ATOMIC
    &END SCF
    &XC
      &XC_FUNCTIONAL BLYP
      &END XC_FUNCTIONAL
    &END XC
  &END DFT
  &SUBSYS
    &CELL
      ABC 9.4 9.4 9.4
#      ABC 5.4 5.4 5.4
      UNIT ANGSTROM
      PERIODIC NONE
    &END CELL
    &COORD
    O   4.700000    4.700000    4.634413
    H   4.700000    3.942864    5.220545
    H   4.700000    5.457136    5.220545
#    O   2.700000    2.700000    2.634413
#    H   2.700000    1.942864    3.220545
#    H   2.700000    3.457136    3.220545
    &END COORD
    &KIND H
      BASIS_SET 6-31G*
      POTENTIAL ALL
    &END KIND
    &KIND O
      BASIS_SET 6-31G*
      POTENTIAL ALL
    &END KIND
  &END SUBSYS
&END FORCE_EVAL
&GLOBAL
  PROJECT H2O
  PRINT_LEVEL MEDIUM
&END GLOBAL




More information about the CP2K-user mailing list