O atom spinpolarized low symmetry: Difference between revisions
No edit summary |
|||
(13 intermediate revisions by 2 users not shown) | |||
Line 1: | Line 1: | ||
{{Template:At_and_mol - Tutorial}} | |||
== Task == | == Task == | ||
Line 7: | Line 7: | ||
== Input == | == Input == | ||
=== {{ | === {{FILE|POSCAR}} === | ||
O atom in a box | O atom in a box | ||
1.0 ! universal scaling parameters | 1.0 ! universal scaling parameters | ||
Line 19: | Line 19: | ||
=== {{ | === {{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 25: | Line 25: | ||
{{TAGBL|ISPIN}} = 2 ! spin polarized calculation | {{TAGBL|ISPIN}} = 2 ! spin polarized calculation | ||
=== {{ | === {{FILE|KPOINTS}} === | ||
Gamma-point only | Gamma-point only | ||
0 | 0 | ||
Line 34: | Line 34: | ||
== Calculation == | == Calculation == | ||
*In the GGA most atoms are characterized by a symmetry broken solution. VASP, however, | *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. | ||
*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. | *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. | ||
Line 46: | Line 46: | ||
DAV: 17 -0.189104076616E+01 -0.51555E-04 -0.34087E-06 48 0.132E-02 | 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 | 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 == | ||
[ | [[Media:Oatomspinlowsym.tgz|Oatomspinlowsym.tgz]] | ||
{{Template:At_and_mol}} | |||
[[Category:Examples]] | [[Category:Examples]] |
Latest revision as of 14:16, 14 November 2019
Overview > O atom > O atom spinpolarized > O atom spinpolarized low symmetry > O dimer > CO > CO vibration > CO partial DOS > H2O >
H2O vibration > H2O molecular dynamics > Further things to try > List of tutorials
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