Fcc Ni: Difference between revisions

From VASP Wiki
No edit summary
Line 59: Line 59:


== Download ==
== Download ==
[http://www.vasp.at/vasp-workshop/examples/2_7_fccNi.tgz 2_7_fccNi.tgz]
[http://www.vasp.at/vasp-workshop/examples/fccNi.tgz fccNi.tgz]


----
----

Revision as of 17:12, 23 August 2016

Description: fcc Ni (a spinpolarized metal). The bash-script loop.sh runs fcc Ni at several different lattice constants (3.0-3.9 Å) and collects free energy versus lattice constant into the file SUMMARY.fcc


  • INCAR
 SYSTEM = fcc Ni
 ISTART = 0 ; ICHARG=2
 ENCUT  =    270
 ISMEAR =    1  ; SIGMA = 0.2
 LORBIT=11

 ISPIN=2
 MAGMOM = 1
  • KPOINTS
k-points
 0
Monkhorst Pack
 11 11 11
 0  0  0
  • POSCAR
fcc:
 3.53 
 0.5 0.5 0.0
 0.0 0.5 0.5
 0.5 0.0 0.5
   1
cartesian
0 0 0
  • And a simple bash-script to loop over several lattice constants: loop.sh
#! /bin/bash
BIN=/path/to/your/vasp/executable
rm WAVECAR SUMMARY.fcc
for i in 3.0 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9 ; do
cat >POSCAR <<!
fcc:
   $i
 0.5 0.5 0.0
 0.0 0.5 0.5
 0.5 0.0 0.5
   1
cartesian
0 0 0
!
echo "a= $i" ; mpirun -n 2 $BIN
E=`awk '/F=/ {print $0}' OSZICAR` ; echo $i $E  >>SUMMARY.fcc
done
cat SUMMARY.fcc

Mind: You will have to set the correct path to your VASP executable (i.e., BIN), and invoke VASP with the correct command (e.g., in the above: mpirun -np 2).

  • To make a quick plot of SUMMARY.fcc try:
gnuplot
gnuplot> plot "SUMMARY.fcc" using ($1):($4) w lp

Download

fccNi.tgz


To the list of examples or to the main page