[Pw_forum] Help for reading data
Paolo Giannozzi
giannozz at nest.sns.it
Wed Feb 23 11:48:03 CET 2005
On Tuesday 22 February 2005 19:14, wei luo wrote:
> Now I get the *.dnyG file. But I don't know how to read it and
> identify the Raman vibration mode.
=======================================
Please post from the same e-mail you used to register:
posts appearing to come from non-registered e-mails will
be automatically deleted (effective NOW: this post is the last
one I rescue)
=======================================
"make tools" compiles code "dynmat.x" in the pwtools/ directory.
It reads the dynamical matrix file at q=0 and can
- apply the acoustic sum rule, so that you get omega=0 for
translation modes at q=0
- add the non-analytical terms that are responsible for the TO-LO
splitting, calculate IR cross sections (provided that the dynamical
matrix file contains the effective charges, option "epsil=.true." in ph.x)
- produce a file with eigendisplacements that can be directly plotted
by "molden"
- Raman cross sections can also be calculated but it is much trickier
It shouldn't be difficult to classify phonon modes in a 8-atom system.
In any case, the bunch of numbers following each frequency are, for
each atom: Re(u_x) Im(u_x) Re(u_y) Im(u_y) Re(u_z) Im(u_z), where u
is the normalized eigendisplacement. The imaginary terms are obviously
zero for q=0
Paolo
--
Paolo Giannozzi e-mail: giannozz at nest.sns.it
Scuola Normale Superiore Phone: +39/050-509876, Fax:-563513
Piazza dei Cavalieri 7 I-56126 Pisa, Italy
More information about the Pw_forum
mailing list