O atom

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 standard calculation for a single oxygen atom in a box. Getting to know the main input and output files of VASP.

Input

POSCAR

O atom in a box
 1.0          ! universal scaling parameters
 8.0 0.0 0.0  ! lattice vector  a(1)
 0.0 8.0 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

We are using a POSCAR file with a single atom. Sufficiently large lattice parameters are selected so that no (significant) interactions between atoms in neighbouring cells is present.

INCAR

SYSTEM = O atom in a box
ISMEAR = 0  ! Gaussian smearing

KPOINTS

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

For atoms or molecules a single k point is sufficient. When more k points are used only the interaction between atoms (which should be zero) is described more accurately.

Calculation

stdout

running on    8 total cores
distrk:  each k-point on    8 cores,    1 groups
distr:  one band on    1 cores,    8 groups
using from now: INCAR     
vasp.5.4.1 05Feb16 (build Aug 22 2016 16:46:23) complex                        
POSCAR found :  1 types and       1 ions
scaLAPACK will be used
LDA part: xc-table for Pade appr. of Perdew
POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: small aliasing (wrap around) errors must be expected
FFT: planning ...
WAVECAR not read
entering main loop
N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.384469664751E+02    0.38447E+02   -0.96726E+02    16   0.293E+02
DAV:   2     0.345965628955E+01   -0.34987E+02   -0.34942E+02    32   0.450E+01
DAV:   3    -0.244485866931E+00   -0.37041E+01   -0.34307E+01    16   0.308E+01
DAV:   4    -0.312557021227E+00   -0.68071E-01   -0.66914E-01    16   0.508E+00
DAV:   5    -0.313520305300E+00   -0.96328E-03   -0.96311E-03    32   0.506E-01    0.286E-01
DAV:   6    -0.314540466589E+00   -0.10202E-02   -0.17853E-03    16   0.332E-01    0.142E-01
DAV:   7    -0.314637222361E+00   -0.96756E-04   -0.22710E-04    16   0.134E-01
1 F= -.31463722E+00 E0= -.16037490E+00  d E =-.308525E+00
writing wavefunctions

The example output (stdout) for the O atom was obtained using VASP version 5.4.1. The initial charge corresponds to the charge of isolated overlapping atoms (POTCAR file). For the first 4 steps the charge remains fixed, then the charge is updated (rms(c) column)

Short explanation of the symbols in the OSZICAR and stdout file:

N	iteration count
E	total energy
dE	change of total energy
d	eps change of the eiegenvalues (fixed potential)
ncg	number of optimisation steps $\mathrm{H} \psi$
rms	total residual vector $\sum_{nk} w_{k} f_{nk} (\mathrm{H} - \epsilon_{nk}) \psi_{nk}$
rms(c)	charge density residual vector

OUTCAR

The individual parts of the OUTCAR file are separated by lines.

----------------------------------------------------------

The OUTCAR file is divided into the following parts:

Reading INCAR, POTCAR, POSCAR

Nearest neighbor distances and analysis of symmetry

Verbose job information

Information on lattice, k points and positions

Information on the basis set (number of plane waves)

Non-local pseudo potential information

Information for each electronic step (one line in OSZICAR)

Timing and energy information

   POTLOK:  cpu time    0.0878: real time    0.0877
   SETDIJ:  cpu time    0.0015: real time    0.0014
    EDDAV:  cpu time    0.0267: real time    0.0434
      DOS:  cpu time    0.0001: real time    0.0001
   --------------------------------------------
     LOOP:  cpu time    0.1165: real time    0.1346
   
eigenvalue-minimisations :    16
total energy-change (2. order) : 0.3844697E+02  (-0.9672571E+02)
number of electron       6.0000000 magnetization
augmentation part        6.0000000 magnetization
    
Free energy of the ion-electron system (eV)      
 ---------------------------------------------------
 alpha Z        PSCENC =         0.27135287
 Ewald energy   TEWEN  =       -91.92708002
 -Hartree energ DENC   =      -281.84385691
 -exchange      EXHF   =         0.00000000
 -V(xc)+E(xc)   XCENC  =        26.11948841
 PAW double counting   =       245.99840262     -247.84808825
 entropy T*S    EENTRO =        -0.08636665
 eigenvalues    EBANDS =       -44.50008162
 atomic energy  EATOM  =       432.26319604
 Solvation  Ediel_sol  =         0.00000000
 ---------------------------------------------------
 free energy    TOTEN  =        38.44696648 eV
    
 energy without entropy =       38.53333313  energy(sigma->0) =       38.49014980

Information on the Eigenvalues

E-fermi :  -8.8431     XC(G=0):  -0.8043     alpha+bet : -0.1463
k-point     1 :       0.0000    0.0000    0.0000
band No.  band energies     occupation
1     -23.8439      2.00000
2      -8.9040      1.33333
3      -8.9040      1.33333
4      -8.9040      1.33333
5      -0.4676      0.00000
6       1.8633      0.00000
7       1.8633      0.00000
8       1.8633      0.00000

Information on stress tensor

The O atom (Example: Oatom)
FORCE on cell =-STRESS in cart. coord.  units (eV):
Direction    XX          YY          ZZ          XY          YZ          ZX
--------------------------------------------------------------------------------------
Alpha Z     0.27135     0.27135     0.27135
Ewald     -30.64236   -30.64236   -30.64236     0.00000     0.00000     0.00000
Hartree    93.90244    93.90244    93.90244    -0.00000    -0.00000    -0.00000
E(xc)     -27.93035   -27.93035   -27.93035    -0.00000    -0.00000    -0.00000
Local    -147.86211  -147.86211  -147.86211     0.00000     0.00000     0.00000
n-local   -20.54942   -20.54942   -20.54942    -0.00000    -0.00000    -0.00000
augment     5.55366     5.55366     5.55366     0.00000    -0.00000     0.00000
Kinetic   126.50998   126.50998   126.50997    -0.00000     0.00000    -0.00000
Fock        0.00000     0.00000     0.00000     0.00000     0.00000     0.00000
-------------------------------------------------------------------------------------
Total      -0.74681    -0.74681    -0.74681     0.00000    -0.00000    -0.00000
in kB      -2.33695    -2.33695    -2.33695     0.00000    -0.00000    -0.00000
external pressure =       -2.34 kB  Pullay stress =        0.00 kB

Information on the energy

FREE ENERGIE OF THE ION-ELECTRON SYSTEM (eV)
---------------------------------------------------
free  energy   TOTEN  =        -0.31463722 eV
energy  without entropy=       -0.00611258  energy(sigma->0) =       -0.16037490

The relevant energy for molecules and atoms is energy without entropy.

energy  without entropy=     -0.00611258  energy(sigma->0) =     -0.16037490

Three degenerate p orbitals are occupied by 2/3 electrons causing an unphysical electronic entropy

entropy T*S    EENTRO =      -0.30852464

A tiny value of SIGMA=0.01 would reduce the entropy but might slow convergence (default is SIGMA=0.2). SIGMA controls the electronic temperature, which is not a very meaningful quantity for molecules and atoms.

The total energy is found to be essentially zero. VASP subtracts from any calculated energy the energy of the atom in the configuration for which the pseudo potential was generated. All pseudo potentials were generated using non spin-polarized reference atoms.

Restart of the calculation

When VASP is restarted the WAVECAR file is read and the run is continued from the previous wave functions (converging rapidly).

running on    8 total cores
distrk:  each k-point on    8 cores,    1 groups
distr:  one band on    1 cores,    8 groups
using from now: INCAR
vasp.5.4.1 05Feb16 (build Aug 22 2016 16:46:23) complex
  
POSCAR found :  1 types and       1 ions
scaLAPACK will be used
LDA part: xc-table for Pade appr. of Perdew
found WAVECAR, reading the header
POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: small aliasing (wrap around) errors must be expected
FFT: planning ...
reading WAVECAR
the WAVECAR file was read successfully
initial charge from wavefunction
entering main loop
      N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.314680766875E+00   -0.31468E+00   -0.83090E-05    16   0.564E-02    0.107E-02
DAV:   2    -0.314677281013E+00    0.34859E-05   -0.10030E-05    16   0.198E-02
  1 F= -.31467728E+00 E0= -.16041496E+00 d E =-.308525E+00
writing wavefunctions

Download

Oatom.tgz

To the list of tutorials or to the main page