Ni 100 surface relaxation: Difference between revisions

From VASP Wiki
Line 70: Line 70:


== Calculation ==
== Calculation ==
*The sample output for the forces in the {{TAG|OUTCAR}} file should look like this (first and last step):
**First step:
POSITION                                      TOTAL-FORCE  (eV/Angst)
-----------------------------------------------------------------------------------
      0.00000      0.00000      0.00000        0.000000      0.000000      0.391352
      0.00000      1.76500      1.76500        0.000000      0.000000    -0.397024
      0.00000      0.00000      3.53000        0.000000      0.000000      0.005117
      0.00000      1.76500      5.29500        0.000000      0.000000      0.391161
      0.00000      0.00000      7.06000        0.000000      0.000000    -0.390607
-----------------------------------------------------------------------------------
    total drift:                                0.000000      0.000000      0.016391
**Last step:
POSITION                                      TOTAL-FORCE  (eV/Angst)
-----------------------------------------------------------------------------------
      0.00000      0.00000      0.00000        0.000000      0.000000      0.399012
      0.00000      1.76500      1.76500        0.000000      0.000000    -0.377003
      0.00000      0.00000      3.53000        0.000000      0.000000      0.105799
      0.00000      1.76500      5.32685        0.000000      0.000000    -0.062054
      0.00000      0.00000      7.02377        0.000000      0.000000    -0.065753


== Download ==
== Download ==

Revision as of 10:18, 15 May 2017

Task

Relaxation of the first two layers of a Ni (100) surface.

Input

POSCAR

fcc (100) surface
 3.53
   .50000   .50000   .00000
  -.50000   .50000   .00000
   .00000   .00000  5.00000
  5
Selective Dynamics
Kartesisch
   .00000   .00000   .00000 F F F
   .00000   .50000   .50000 F F F
   .00000   .00000  1.00000 F F F
   .00000   .50000  1.50000 T T T
   .00000   .00000  2.00000 T T T
  • Ni lattice constant of 3.53.
  • 1 atom per layer: p(1x1) cell.
  • 5 nickel layers.
  • First two layers (of one side) relaxed.
  • vacuum.

INCAR

  ISTART = 0; ICHARG = 2
    
general:
  SYSTEM = clean Ni(100) surface
  ENCUT = 270 
  ISMEAR = 2 ; SIGMA = 0.2
  ALGO = Fast
  EDIFF = 1E-6
    
spin:
  ISPIN=2
  MAGMOM = 5*1
    
dynamic:
  NSW = 100
  POTIM = 0.8
  IBRION = 1
  • Initial charge-density in startjob from overlapping atoms.
  • Default energy cut-off of 270 eV.
  • MP-smearing (metal).
  • Spin-polarized calculation with initial moment of 1.
  • Ionic relaxation used.

KPOINTS

k-points
0
Monkhorst-Pack
9 9 1
0 0 0
  • Equally spaced mesh.
  • Odd mesh, centered on .
  • 15 k points in irreducible Brillouin zone (IBZ).
  • Only one k point in z-direction for surface.

Calculation

  • The sample output for the forces in the OUTCAR file should look like this (first and last step):
    • First step:
POSITION                                       TOTAL-FORCE  (eV/Angst)
-----------------------------------------------------------------------------------
     0.00000      0.00000      0.00000         0.000000      0.000000      0.391352
     0.00000      1.76500      1.76500         0.000000      0.000000     -0.397024
     0.00000      0.00000      3.53000         0.000000      0.000000      0.005117
     0.00000      1.76500      5.29500         0.000000      0.000000      0.391161
     0.00000      0.00000      7.06000         0.000000      0.000000     -0.390607
-----------------------------------------------------------------------------------
   total drift:                                0.000000      0.000000      0.016391
    • Last step:
POSITION                                       TOTAL-FORCE  (eV/Angst)
-----------------------------------------------------------------------------------
     0.00000      0.00000      0.00000         0.000000      0.000000      0.399012
     0.00000      1.76500      1.76500         0.000000      0.000000     -0.377003
     0.00000      0.00000      3.53000         0.000000      0.000000      0.105799
     0.00000      1.76500      5.32685         0.000000      0.000000     -0.062054
     0.00000      0.00000      7.02377         0.000000      0.000000     -0.065753

Download

Ni100clean_rel.tgz


To the list of examples or to the main page