H2O molecular dynamics: Difference between revisions
Vaspmaster (talk | contribs) No edit summary |
No edit summary |
||
| Line 1: | Line 1: | ||
*INCAR | *{{TAG|INCAR}} | ||
PREC = Normal ! standard precision | PREC = Normal ! standard precision | ||
ENMAX = 400 ! cutoff should be set manually | ENMAX = 400 ! cutoff should be set manually | ||
| Line 10: | Line 10: | ||
TEBEG = 2000 ; TEEND = 2000 ! temperature | TEBEG = 2000 ; TEEND = 2000 ! temperature | ||
*KPOINTS | *{{TAG|KPOINTS}} | ||
Gamma-point only | Gamma-point only | ||
0 | 0 | ||
| Line 17: | Line 17: | ||
0 0 0 | 0 0 0 | ||
*POSCAR | *{{TAG|POSCAR}} | ||
H2O _2 | H2O _2 | ||
0.52918 ! scaling parameter | 0.52918 ! scaling parameter | ||
| Line 57: | Line 57: | ||
gnuplot -persist plotfile | gnuplot -persist plotfile | ||
</pre> | </pre> | ||
== Used INCAR Tags == | |||
{{TAG|EDIFF}}, {{TAG|EDIFFG}}, {{TAG|ENMAX}}, {{TAG|IBRION}}, {{TAG|ISMEAR}}, {{TAG|ISYM}}, {{TAG|NSW}}, {{TAG|POTIM}}, {{TAG|PREC}}, {{TAG|SMASS}}, {{TAG|TEBEG}}, {{TAG|TEEND}} | |||
== Download == | == Download == | ||
[http://www.vasp.at/vasp-workshop/examples/H2Omd.tgz H2Omd.tgz] | [http://www.vasp.at/vasp-workshop/examples/H2Omd.tgz H2Omd.tgz] | ||
---- | ---- | ||
[[VASP_example_calculations|To the list of examples]] or to the [[The_VASP_Manual|main page]] | [[VASP_example_calculations|To the list of examples]] or to the [[The_VASP_Manual|main page]] | ||
[[Category:Examples]] | [[Category:Examples]] | ||
Revision as of 10:17, 16 February 2017
PREC = Normal ! standard precision ENMAX = 400 ! cutoff should be set manually ISMEAR = 0 ; SIGMA = 0.1 ISYM = 0 ! strongly recommended for MD IBRION = 0 ! molecular dynamics NSW = 100 ! 100 steps POTIM = 1.0 ! timestep 1 fs SMASS = -3 ! Nose Hoover thermostat TEBEG = 2000 ; TEEND = 2000 ! temperature
Gamma-point only 0 Monkhorst Pack 1 1 1 0 0 0
H2O _2
0.52918 ! scaling parameter
12 0 0
0 12 0
0 0 12
1 2
select
cart
0.00 0.00 0.00 T T F
1.10 -1.43 0.00 T T F
1.10 1.43 0.00 T T F
the pair correlation function can be visualized using e.g. the following script:
- plot_PCDAT
awk <PCDAT >PCDAT.dat '
NR==8 { pcskal=$1}
NR==9 { pcfein=$1}
NR>=13 {
line=line+1
if (line==257) {
print " "
line=0
}
else
print (line-0.5)*pcfein/pcskal,$1
}
'
cat >plotfile<<!
# set term postscript enhanced colour lw 2 "Helvetica" 20
# set output "pair_correlation.eps"
set title "pair-correlation of H2O at 2000 K"
set xlabel "r [Angstrom]"
set ylabel "g(r)"
plot [0:15] "PCDAT.dat" w lines
!
gnuplot -persist plotfile
Used INCAR Tags
EDIFF, EDIFFG, ENMAX, IBRION, ISMEAR, ISYM, NSW, POTIM, PREC, SMASS, TEBEG, TEEND
Download
To the list of examples or to the main page