H2O vibration: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
Line 1: Line 1:
*{{TAG|INCAR}}
*{{TAG|INCAR}}
<pre>
{{TAGBL|SYSTEM}} = H2O vibration
SYSTEM = H2O vibration
{{TAGBL|PREC}} = A
PREC = A
# {{TAGBL|IBRION}} = 1 ; {{TAGBL|NSW}} = 10 ; {{TAGBL|NFREE}} = 2 ; {{TAGBL|EDIFFG}} = -1E-4
# IBRION = 1 ; NSW = 10 ; NFREE = 2 ; EDIFFG = -1E-4
{{TAGBL|ENMAX}} = 400
ENMAX = 400
{{TAGBL|ISMEAR}} = 0    # Gaussian smearing
ISMEAR = 0    # Gaussian smearing
{{TAGBL|IBRION}} = 6    # finite differences with symmetry
IBRION = 6    # finite differences with symmetry
{{TAGBL|NFREE}} = 2    # central differences (default)
NFREE = 2    # central differences (default)
{{TAGBL|POTIM}} = 0.015 # default as well
POTIM = 0.015 # default as well
{{TAGBL|EDIFF}} = 1E-8  
EDIFF = 1E-8  
{{TAGBL|NSW}} = 1      # ionic steps > 0
NSW = 1      # ionic steps > 0
</pre>


*{{TAG|KPOINTS}}
*{{TAG|KPOINTS}}
Line 42: Line 40:
cutoff {{TAG|ENCUT}} to 800 eV. The important and physically meaningful
cutoff {{TAG|ENCUT}} to 800 eV. The important and physically meaningful
frequencies are, however, insensitive to the choice of the cutoff.
frequencies are, however, insensitive to the choice of the cutoff.
== Used INCAR Tags ==
{{TAG|EDIFF}}, {{TAG|EDIFFG}}, {{TAG|ENMAX}}, {{TAG|IBRION}}, {{TAG|ISMEAR}}, {{TAG|NFREE}}, {{TAG|NSW}}, {{TAG|PREC}}, {{TAG|SYSTEM}}


== Download ==
== Download ==

Revision as of 10:48, 15 March 2017

SYSTEM = H2O vibration
PREC = A
# IBRION = 1 ; NSW = 10 ; NFREE = 2 ; EDIFFG = -1E-4
ENMAX = 400
ISMEAR = 0    # Gaussian smearing
IBRION = 6    # finite differences with symmetry
NFREE = 2     # central differences (default)
POTIM = 0.015 # default as well
EDIFF = 1E-8 
NSW = 1       # ionic steps > 0
Gamma-point only
 0
Monkhorst Pack
 1 1 1
 0 0 0
H2O _2                                                                          
1.0000000
  8.0000000   0.0000000   0.0000000
  0.0000000   8.0000000   0.0000000
  0.0000000   0.0000000   8.0000000
   1    2
cart
  0.0000000   0.0000000   0.0000000
  0.5960812  -0.7677068   0.0000000
  0.5960812   0.7677068   0.0000000

How many zero frequency modes should be observed and why? Try to use the linear response code (IBRION=8 and EDIFF=1E-8) to obtain reference results. For finite differences, are the results sensitive to the step width POTIM. In this specific case, the drift in the forces is too large to obtain the zero frequency modes "exactly", and it is simplest to increase the cutoff ENCUT to 800 eV. The important and physically meaningful frequencies are, however, insensitive to the choice of the cutoff.

Download

H2Ovib.tgz


To the list of examples or to the main page