O atom spinpolarized low symmetry: Difference between revisions

From VASP Wiki
No edit summary
 
(19 intermediate revisions by 2 users not shown)
Line 1: Line 1:
<div style="background:#FFFFFF; border:solid 1px #AAAAAA; margin-bottom:1em; padding:0.4em 0.6em;">[[O atom tutorial|Overview]]&nbsp;<span style="color:#AAAAAA;">&gt;</span> [[O atom|O atom]]&nbsp;<span style="color:#AAAAAA;">&gt;</span> [[O atom spinpolarized|O atom spinpolarized]]&nbsp;<span style="color:#AAAAAA;">&gt;</span> [[O atom spinpolarized low symmetry|O atom spinpolarized low symmetry]]&nbsp;<span style="color:#AAAAAA;">&gt;</span> [[/4|Videotutorial]]</div>
{{Template:At_and_mol - Tutorial}}


*{{TAG|INCAR}}
== Task ==
 
Performing a spin-polarized low symmetry calculation of a single oxygen atom in a non cubic box to get the right energy minimum.
 
== Input ==
 
=== {{FILE|POSCAR}} ===
O atom in a box
  1.0          ! universal scaling parameters
  7.0 0.0 0.0  ! lattice vector  a(1)
  0.0 7.5 0.0  ! lattice vector  a(2)
  0.0 0.0 8.0  ! lattice vector  a(3)
1            ! number of atoms
cart          ! positions in cartesian coordinates
  0 0 0
 
 
 
=== {{FILE|INCAR}} ===
  {{TAGBL|SYSTEM}} = O atom in a box
  {{TAGBL|SYSTEM}} = O atom in a box
  {{TAGBL|ISMEAR}} = 0  ! Gaussian smearing
  {{TAGBL|ISMEAR}} = 0  ! Gaussian smearing
Line 7: Line 25:
  {{TAGBL|ISPIN}} =  2  ! spin polarized calculation
  {{TAGBL|ISPIN}} =  2  ! spin polarized calculation


*{{TAG|KPOINTS}}
=== {{FILE|KPOINTS}} ===
  Gamma-point only
  Gamma-point only
   0
   0
Line 14: Line 32:
   0 0 0
   0 0 0


*{{TAG|POSCAR}}
== Calculation ==
  O atom in a box
 
   1.0          ! universal scaling parameters
*In the GGA most atoms are characterized by a symmetry broken solution. VASP, however, symmetrizes the charge density according to the determined symmetry of the cell. Check the {{TAG|OUTCAR}} file, to see what symmetry VASP is using.
   7.0 0.0 0.0  ! lattice vector a(1)
 
   0.0 7.5 0.0  ! lattice vector a(2)
*To lower the symmetry, simply change the lattice parameters to 7.0, 7.5 and 8.0 in the {{TAG|POSCAR}} file (see the example file above) and reduce {{TAG|SIGMA}} to {{TAG|SIGMA}}=0.01 in the {{TAG|INCAR}} file.
   0.0 0.0 8.0 ! lattice vector a(3)
 
1            ! number of atoms
*By rerunning VASP one finds a much lower energy:
cart          ! positions in cartesian coordinates
 
  0 0 0
  vasp.5.4.1 05Feb16 (build Aug 22 2016 16:46:23) complex
...   ...   . ..
DAV:  15    -0.189071145737E+01  -0.29321E-03  -0.39183E-05    48  0.478E-02    0.995E-03
  DAV: 16    -0.189071145737E+01   -0.27775E-03  -0.39294E-05    40  0.290E-02    0.541E-03
  DAV: 17    -0.189104076616E+01   -0.51555E-04  -0.34087E-06    48  0.132E-02   
    1 F= -.18910408E+01 E0= -.18910408E+01 d E =-.309633E-20 mag=    1.9998
 
== Further things to try ==
 
*How does the energy change when one decreases {{TAG|SIGMA}} to {{TAG|SIGMA}}=0.01 in the {{TAG|INCAR}} file? Why?


== Download ==
== Download ==
[http://www.vasp.at/vasp-workshop/examples/Oatomspinlowsym.tgz Oatomspinlowsym.tgz]
[[Media:Oatomspinlowsym.tgz|Oatomspinlowsym.tgz]]
----
[[VASP_example_calculations|To the list of examples]] or to the [[The_VASP_Manual|main page]]
{{Template:At_and_mol}}


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

Latest revision as of 14:16, 14 November 2019

Task

Performing a spin-polarized low symmetry calculation of a single oxygen atom in a non cubic box to get the right energy minimum.

Input

POSCAR

O atom in a box
 1.0          ! universal scaling parameters
 7.0 0.0 0.0  ! lattice vector  a(1)
 0.0 7.5 0.0  ! lattice vector  a(2)
 0.0 0.0 8.0  ! lattice vector  a(3)
1             ! number of atoms
cart          ! positions in cartesian coordinates
 0 0 0


INCAR

SYSTEM = O atom in a box
ISMEAR = 0  ! Gaussian smearing
SIGMA  = 0.01
ISPIN =  2  ! spin polarized calculation

KPOINTS

Gamma-point only
 0
Monkhorst Pack
 1 1 1
 0 0 0

Calculation

  • In the GGA most atoms are characterized by a symmetry broken solution. VASP, however, symmetrizes the charge density according to the determined symmetry of the cell. Check the OUTCAR file, to see what symmetry VASP is using.
  • To lower the symmetry, simply change the lattice parameters to 7.0, 7.5 and 8.0 in the POSCAR file (see the example file above) and reduce SIGMA to SIGMA=0.01 in the INCAR file.
  • By rerunning VASP one finds a much lower energy:
vasp.5.4.1 05Feb16 (build Aug 22 2016 16:46:23) complex
...   ...   . ..
DAV:  15    -0.189071145737E+01   -0.29321E-03   -0.39183E-05    48   0.478E-02    0.995E-03
DAV:  16    -0.189071145737E+01   -0.27775E-03   -0.39294E-05    40   0.290E-02    0.541E-03
DAV:  17    -0.189104076616E+01   -0.51555E-04   -0.34087E-06    48   0.132E-02    
   1 F= -.18910408E+01 E0= -.18910408E+01  d E =-.309633E-20  mag=     1.9998

Further things to try

  • How does the energy change when one decreases SIGMA to SIGMA=0.01 in the INCAR file? Why?

Download

Oatomspinlowsym.tgz