Including the Spin-Orbit Coupling: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
 
(12 intermediate revisions by 4 users not shown)
Line 1: Line 1:
Description: Spin-Orbit Coupling (SOC) included self-consistently
Description: Spin-Orbit Coupling (SOC) included self-consistently


''<u>Exercise :</u>'' Determine the total magnetic moment by adding the orbital moment of the Ni atoms. Change the orientation of spins using the [[SAXIS|SAXIS]]-tag.
The Magnetocrystalline Anisotropy Energy is determined by rotating all spins according to different directions. To modify the orientation of the spins in the crystal, we consider the second approach described [[SAXIS|here]]. For the [[MAGMOM|MAGMOM]]-tag, the total local magnetic moment is written according to the z-direction (necessarily, the x and y-directions are equal to 0). The spin orientation <math>(u,v,w)</math> is defined by the [[SAXIS|SAXIS]]-tag in the Cartesian frame. The Magnetocrystalline Anisotropy Energy is calculated by orientating the spins in different directions and the following equation
----
*INCAR
<pre>
NiO GGA+U SOC
  SYSTEM    = "NiO"


Electronic minimization
::<math>E_{\text{MAE}} = E_{(u,v,w)} - E_0</math>
  ENCUT    = 450
  EDIFF    = 1E-5
  LORBIT    = 11
  LREAL    = .False.
  ISTART    = 0
  ISYM      = -1
  NELMIN    = 6
  LSORBIT  = .True.
  LWAVE    = .False.
  LCHARG  = .False.


DOS
with <math>E_0</math> the energy of the most stable spin orientation.
  ISMEAR    = -5


Magnetism
More details are available in the [[SAXIS|SAXIS]] and [[LSORBIT|LSORBIT]] pages.
  ISPIN    = 2
  MAGMOM    = 0 0 2 0 0 -2 6*0
  SAXIS     = 1 1 1


Orbital Moment
''<u>Exercise :</u>'' Determine the total magnetic moment by adding the orbital moment of the Ni atoms. Calculate the Magnetocrystalline Anisotropy Energy of NiO by orientating the spins along the following path : (2,2,2) --> (2,2,1) --> (2,2,0) --> ... --> (2,2,-6). Identify the most stable spin orientation according to this path.
  LORBMOM  = T  
----
*{{TAG|INCAR}}
NiO GGA+U SOC
  {{TAGBL|SYSTEM}}    = "NiO"
   
Electronic minimization
  {{TAGBL|ENCUT}}    = 450
  {{TAGBL|EDIFF}}    = 1E-7
  {{TAGBL|LORBIT}}    = 11
  {{TAGBL|LREAL}}    = .False.
  {{TAGBL|ISTART}}    = 0
  {{TAGBL|ISYM}}      = -1
  {{TAGBL|NELMIN}}    = 6
  {{TAGBL|LSORBIT}}  = .True.
  {{TAGBL|LWAVE}}    = .False.
  {{TAGBL|LCHARG}}    = .False.
   
DOS
  {{TAGBL|ISMEAR}}    = -5
     
Magnetism
  {{TAGBL|ISPIN}}    = 2
  {{TAGBL|MAGMOM}}    = 0 0 2 0 0 -2 6*0
  {{TAGBL|SAXIS}}    = 2 2 2
     
Orbital Moment
  {{TAGBL|LORBMOM}}   = T  
     
Mixer
  {{TAGBL|AMIX}}      = 0.2
  {{TAGBL|BMIX}}      = 0.00001
  {{TAGBL|AMIX_MAG}}  = 0.8
  {{TAGBL|BMIX_MAG}}  = 0.00001
   
GGA+U
  {{TAGBL|LDAU}}      = .TRUE.
  {{TAGBL|LDAUTYPE}}  = 2
  {{TAGBL|LDAUL}}    = 2 -1
  {{TAGBL|LDAUU}}    = 5.00 0.00
  {{TAGBL|LDAUJ}}    = 0.00 0.00
  {{TAGBL|LDAUPRINT}} = 1
  {{TAGBL|LMAXMIX}}  = 4


Mixer
*{{TAG|KPOINTS}}
  AMIX      = 0.2
  BMIX      = 0.00001
  AMIX_MAG  = 0.8
  BMIX_MAG  = 0.00001
 
GGA+U
  LDAU      = .TRUE.
  LDAUTYPE  = 2
  LDAUL    = 2 -1
  LDAUU    = 5.00 0.00
  LDAUJ    = 0.00 0.00
  LDAUPRINT = 2
  LMAXMIX  = 4
</pre>
 
*KPOINTS
<pre>
<pre>
k-points
k-points
Line 56: Line 62:
</pre>
</pre>


*POSCAR
*{{TAG|POSCAR}}
<pre>
<pre>
NiO
NiO
Line 72: Line 78:


== Download ==
== Download ==
[http://www.vasp.at/vasp-workshop/examples/4_3_NiO_LSDA+U.tgz 4_3_NiO_LSDA+U.tgz]
[[Media:nio_SOC.tgz| nio_SOC.tgz]]
 
----
----
[[VASP_example_calculations|To the list of examples]] or to the [[The_VASP_Manual|main page]]


[[Category:Examples]]
[[Category:Examples]]

Latest revision as of 08:48, 11 April 2023

Description: Spin-Orbit Coupling (SOC) included self-consistently

The Magnetocrystalline Anisotropy Energy is determined by rotating all spins according to different directions. To modify the orientation of the spins in the crystal, we consider the second approach described here. For the MAGMOM-tag, the total local magnetic moment is written according to the z-direction (necessarily, the x and y-directions are equal to 0). The spin orientation is defined by the SAXIS-tag in the Cartesian frame. The Magnetocrystalline Anisotropy Energy is calculated by orientating the spins in different directions and the following equation

with the energy of the most stable spin orientation.

More details are available in the SAXIS and LSORBIT pages.

Exercise : Determine the total magnetic moment by adding the orbital moment of the Ni atoms. Calculate the Magnetocrystalline Anisotropy Energy of NiO by orientating the spins along the following path : (2,2,2) --> (2,2,1) --> (2,2,0) --> ... --> (2,2,-6). Identify the most stable spin orientation according to this path.


NiO GGA+U SOC
 SYSTEM    = "NiO"
    
Electronic minimization
  ENCUT     = 450
  EDIFF     = 1E-7
  LORBIT    = 11
  LREAL     = .False.
  ISTART    = 0
  ISYM      = -1
  NELMIN    = 6
  LSORBIT   = .True.
  LWAVE     = .False.
  LCHARG    = .False.
    
DOS
  ISMEAR    = -5
     
Magnetism
  ISPIN     = 2
  MAGMOM    = 0 0 2 0 0 -2 6*0
  SAXIS     = 2 2 2
     
Orbital Moment
  LORBMOM   = T 
     
Mixer
  AMIX      = 0.2
  BMIX      = 0.00001
  AMIX_MAG  = 0.8
  BMIX_MAG  = 0.00001
    
GGA+U
  LDAU      = .TRUE.
  LDAUTYPE  = 2
  LDAUL     = 2 -1
  LDAUU     = 5.00 0.00
  LDAUJ     = 0.00 0.00
  LDAUPRINT = 1
  LMAXMIX   = 4 
k-points
 0
gamma
 4  4  4 
 0  0  0
NiO
 4.17
 1.0 0.5 0.5
 0.5 1.0 0.5
 0.5 0.5 1.0
 2 2
Cartesian
 0.0 0.0 0.0
 1.0 1.0 1.0
 0.5 0.5 0.5
 1.5 1.5 1.5

Download

nio_SOC.tgz