Iron surface

marcello rosini.... at gmail.com
Tue Mar 8 16:28:58 UTC 2011


Hallo everybody...
I am new to cp2k so forgive if I my questions may sound silly...
After some succesfull calculations on smaller systems I am trying to
setup the simulation of a Fe surface, but my calculation does not
converge. Actually the scf cycle does not converge

Q1- The basis for Fe is the only one I found on the CP2k source
tree... is it correct to mix PADE basis with PBE potential and XC?
Q2- Are the parameters controlling the scf correct? I cannot manage
them to obtain convergence...
Q3- is the dimension of the surface too small for a metal?
Q4- does anyone have some suggestion/basis-set/potential for my
system?

Thanks in advance to everybody
Marcello

----------------
this is my input:

&GLOBAL
  RUN_TYPE GEO_OPT
&END GLOBAL
!
&FORCE_EVAL
  METHOD Quickstep
  &DFT
    BASIS_SET_FILE_NAME MIOBASIS
    POTENTIAL_FILE_NAME MIOPOT
    LSD
    &MGRID
      CUTOFF 240
    &END MGRID
    &QS
      EPS_DEFAULT 1.0E-14
      EXTRAPOLATION PS
    &END QS
    &SCF
      EPS_SCF 4.0E-7
      SCF_GUESS atomic
      MAX_SCF 500
      CHOLESKY INVERSE
      ADDED_MOS 120 120
      &DIAGONALIZATION on
        ALGORITHM standard
      &END DIAGONALIZATION
      &MIXING on
        METHOD broyden_mixing
        ALPHA 0.05
        BETA 1.5
        NBROYDEN 8
      &END MIXING
      &SMEAR on
         METHOD energy_window
         WINDOW_SIZE [Ry] 0.02
      &END SMEAR
    &END SCF
    &XC
      &XC_FUNCTIONAL pbe
      &END XC_FUNCTIONAL
    &END XC
  &END DFT
  &SUBSYS
    &CELL
      A 8.06666835 -5.70399589   0.0
      B 8.06666835  5.70399589   0.0
      C 0.0         0.0         28.23
      PERIODIC XYZ
    &END CELL
    &COORD
   Fe   12.10000000   -0.00000315    8.05240000
   Fe   10.08330000    1.42600000    8.05240000
   Fe   10.08330000   -1.42600000    8.05240000
   Fe    8.06667000    0.00000088    8.05240000
   Fe    8.06667000   -2.85200000    8.05240000
   Fe    8.06667000    2.85199000    8.05240000
   Fe    6.05000000    4.27800000    8.05240000
   Fe    6.05000000   -4.27800000    8.05240000
   Fe    6.05000000    1.42600000    8.05240000
   Fe    6.05000000   -1.42600000    8.05240000
   Fe    4.03333000   -0.00000158    8.05240000
   Fe    4.03333000    2.85200000    8.05240000
   Fe    4.03333000   -2.85200000    8.05240000
   Fe    2.01667000    1.42600000    8.05240000
   Fe    2.01667000   -1.42600000    8.05240000
   Fe    0.00000000    0.00000000    8.05240000
   Fe   14.11670000   -0.00000315    6.05311000
   Fe   12.10000000    1.42600000    6.05311000
   Fe   12.10000000   -1.42600000    6.05311000
   Fe   10.08330000    0.00000088    6.05311000
   Fe   10.08330000   -2.85200000    6.05311000
   Fe   10.08330000    2.85199000    6.05311000
   Fe    8.06667000    4.27800000    6.05311000
   Fe    8.06667000   -4.27800000    6.05311000
   Fe    8.06667000    1.42600000    6.05311000
   Fe    8.06667000   -1.42600000    6.05311000
   Fe    6.05000000   -0.00000158    6.05311000
   Fe    6.05000000    2.85200000    6.05311000
   Fe    6.05000000   -2.85200000    6.05311000
   Fe    4.03333000    1.42600000    6.05311000
   Fe    4.03333000   -1.42600000    6.05311000
   Fe    2.01667000    0.00000000    6.05311000
   Fe   12.10000000   -0.00000315    4.03452000   BOT
   Fe   10.08330000    1.42600000    4.03452000   BOT
   Fe   10.08330000   -1.42600000    4.03452000   BOT
   Fe    8.06667000    0.00000088    4.03452000   BOT
   Fe    8.06667000   -2.85200000    4.03452000   BOT
   Fe    8.06667000    2.85199000    4.03452000   BOT
   Fe    6.05000000    4.27800000    4.03452000   BOT
   Fe    6.05000000   -4.27800000    4.03452000   BOT
   Fe    6.05000000    1.42600000    4.03452000   BOT
   Fe    6.05000000   -1.42600000    4.03452000   BOT
   Fe    4.03333000   -0.00000158    4.03452000   BOT
   Fe    4.03333000    2.85200000    4.03452000   BOT
   Fe    4.03333000   -2.85200000    4.03452000   BOT
   Fe    2.01667000    1.42600000    4.03452000   BOT
   Fe    2.01667000   -1.42600000    4.03452000   BOT
   Fe    0.00000000    0.00000000    4.03452000   BOT
    &END COORD
    &KIND Fe
      BASIS_SET DZV-GTH-PADE
      POTENTIAL GTH-PBE-q16
    &END KIND
  &END SUBSYS
&END FORCE_EVAL
!
&MOTION
  &GEO_OPT
    OPTIMIZER BFGS
    MAX_ITER 200
  &END GEO_OPT
  &CONSTRAINT
    &FIXED_ATOMS
      COMPONENTS_TO_FIX XYZ
      MOLNAME BOT
    &END FIXED_ATOMS
  &END CONSTRAINT
&END MOTION


More information about the CP2K-user mailing list