Spin-orbit coupling in a Ni monolayer

From VASP Wiki
The printable version is no longer supported and may have rendering errors. Please update your browser bookmarks and please use the default browser print function instead.

Task

Spin-orbit coupling (SOC) in a freestanding Ni monolayer.

Input

POSCAR

fcc Ni 100 surface
 3.53
   .50000   .50000   .00000
  -.50000   .50000   .00000
   .00000   .00000  5.00000
  1
Cartesian
   .00000   .00000   .00000

INCAR

SYSTEM        = Ni (100) monolayer
ISTART        = 0
ENCUT         = 270.00
LNONCOLLINEAR = .TRUE.
MAGMOM        = 0.0 0.0 1.0
VOSKOWN       = 1
LSORBIT       = .TRUE.
    
LMAXMIX       = 4
  • Initialization of moment along z-direction (out of plane) (MAGMOM = 0.0 0.0 1.0).
  • Spin-orbit interaction switched on (LSORBIT=.TRUE.).
  • For the second calculation, switch to in-plane magnetization by setting MAGMOM = 1.0 0.0 0.0 in the INCAR file.

KPOINTS

k-points
 0
Monkhorst-Pack
9 9 1
0 0 0

Calculation

  • The sample output for the total energy using out of plane magnetization is given in the OSZICAR file:
...
DAV:  20    -0.371322930070E+01    0.15852E-03   -0.11632E-03   636   0.235E-01    0.225E-02
DAV:  21    -0.371323204989E+01   -0.27492E-05   -0.13047E-05   500   0.184E-02
   1 F= -.37132320E+01 E0= -.37139803E+01  d E =0.224478E-02  mag=  0.0000  0.0000  0.9035
  • The sample output for the total energy using in plane magnetization looks like the following:
...
DAV:  19    -0.371443443024E+01   -0.80757E-04   -0.35822E-03  1084   0.323E-01    0.119E-02
DAV:  20    -0.371446032472E+01   -0.25894E-04   -0.42423E-05   916   0.263E-02
   1 F= -.37144603E+01 E0= -.37150300E+01  d E =0.170900E-02  mag=  0.9049  0.0000  0.0000
  • From the energy difference of these calculations we see that the easy axis lies in plane:

Download

4_4_SOI_Ni.tgz

Retrieved from "https://www.vasp.at/wiki/index.php?title=Spin-orbit_coupling_in_a_Ni_monolayer&oldid=10462"