Page 1 of 1

Phonon Calc at Gamma

Posted: Sun Dec 04, 2005 9:27 pm
by Tiger-paw
1.
To calc Phonon Freq in a crystal supercell at Gamma, I use IBRION=50, NFREE = 2, but how to choose a reasonable POTIM? Will the calculated phonon freq be sensitive to POTIM?

2.
In my calc, all freq I obtained for Gamma are nonzero. Am I supposed to see 3 zero freq and 3N-3 nonzero normal mode freq? What can I do for this?

Many thanks for your help,

tiger-paw

Phonon Calc at Gamma

Posted: Mon Dec 05, 2005 9:59 am
by admin
1) The Hessian Matrix is constructed via displacements from the equilibrium positons as set by POTIM, therefore the valuehas to be small enough to ensure that you are in the harmonic regime (a few hundereths of Ã…).
btw: IBRION =5 (not 50)
2) please choose a reasonable input in POSCAR (selective dynamics, T/F): if all positions are free to be displaced, the number of degrees of freedom is 3N, including the translation of the whole species as well.

Phonon Calc at Gamma

Posted: Sat Feb 11, 2006 7:40 pm
by dfcox
Admin-

Could you be a little more specific about point 1) when IBRION=5? What value of POTIM is appropriate to give a displacement of a few hundreths of an Angstom? Does POTIM represent the magnitude of the displacement in Angstroms?

Thanks

Phonon Calc at Gamma

Posted: Mon Feb 13, 2006 12:46 pm
by admin
1) please read the manual (http://cms.mpi.univie.ac.at/vasp/vasp/node101.html) concerning the
valid settings of IBRION.
2) POTIM is the displacement in Ã… you are right

Phonon Calc at Gamma

Posted: Fri Sep 15, 2006 4:10 pm
by Tiger-paw
Admin-

Could the direct/supercell phonon code be used to obtain all normal mode freqs for a surface slab, with vacuum and even adsorbent on both sides, and hence total free energy? (assuming all atoms are free to be displaced)

Is Gamma-point-only enough (or exact) for this C_1 sysmetry surface slab?

Many thanks!

Phonon Calc at Gamma

Posted: Fri Sep 15, 2006 4:15 pm
by Tiger-paw
Sorry, "adsorbent" should be "adsorbates", i.e. adsorbed atoms or moleculers.

Phonon Calc at Gamma

Posted: Sat Sep 16, 2006 6:51 am
by bandy
why you want to study adsorbates on both sides(surfaces) of slab? Is there any specific reason? because the calculation will be very expensive?

Phonon Calc at Gamma

Posted: Thu Sep 21, 2006 3:48 pm
by Tiger-paw
surface dipoles on two sides can be cancelled.

Phonon Calc at Gamma

Posted: Tue Sep 26, 2006 11:35 am
by admin
to correct for the dipol moments, you can also set IDIPOL=3 if it helps to spare computing time
in principle, you can leave all atoms free to relax and calculate all vibration modes, however, this will be very time-consuming