Dear YC Cheng, I am afraid it will take a bit longer to you to understand how to calculate infrared and Raman intensities, than just reading the QE users' manual. Please, make sure you undrestand first how IR and Raman intensities are given in terms of the Born effective charges and derivatives of the polarizability tensor, respectively. This information is NOT contained in the dynamical matrix. Any standard advanced texbook on lattice dynamics would do. You may want to browse for references in our RMP paper first: see: S. Baroni et al. Rev. Mod. Phys. 73, 515 (2001).
Hope this helps Stefano Baroni On Apr 10, 2008, at 4:00 PM, ??? wrote: > Dear PWscf users: > I am trying to test the IR and Raman properties of silicon. > I followed the instruction of the user's manual(scf, phonon at > gamma, dynmat.x) and got diagonalizing the dynamical matrix. I am > shame that I can't understand the dynamical matrix. The matrix is > list as follows: > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > ====================================================================== > q = 0.0000 0.0000 0.0000 > > ************************************************************************** > omega( 1) = 0.000000 [THz] = 0.000005 [cm-1] > ( 0.000000 0.000000 -0.707107 0.000000 0.000000 > 0.000000 ) > ( 0.000000 0.000000 -0.707107 0.000000 0.000000 > 0.000000 ) > omega( 2) = 0.000000 [THz] = 0.000005 [cm-1] > ( -0.707107 0.000000 0.000000 0.000000 0.000000 > 0.000000 ) > ( -0.707107 0.000000 0.000000 0.000000 0.000000 > 0.000000 ) > omega( 3) = 0.000000 [THz] = 0.000005 [cm-1] > ( 0.000000 0.000000 0.000000 0.000000 -0.707107 > 0.000000 ) > ( 0.000000 0.000000 0.000000 0.000000 -0.707107 > 0.000000 ) > omega( 4) = 15.294056 [THz] = 510.158172 [cm-1] > ( 0.000000 0.000000 0.707107 0.000000 0.000000 > 0.000000 ) > ( 0.000000 0.000000 -0.707107 0.000000 0.000000 > 0.000000 ) > omega( 5) = 15.294056 [THz] = 510.158172 [cm-1] > ( 0.000000 0.000000 0.000000 0.000000 -0.707107 > 0.000000 ) > ( 0.000000 0.000000 0.000000 0.000000 0.707107 > 0.000000 ) > omega( 6) = 15.294056 [THz] = 510.158172 [cm-1] > ( -0.707107 0.000000 0.000000 0.000000 0.000000 > 0.000000 ) > ( 0.707107 0.000000 0.000000 0.000000 0.000000 > 0.000000 ) > > ************************************************************************** > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > = > ====================================================================== > How can I get the IR and Raman mode and the intensity? > Thanks in advance! > > Best regards. > > -- > Y. C. Cheng > Department of Phyics > Nanjing University > Nanjing 210093 > P. R. China > Tel: 86-25-83592907 > Email: yccheng.nju at gmail.com > _______________________________________________ > Pw_forum mailing list > Pw_forum at pwscf.org > http://www.democritos.it/mailman/listinfo/pw_forum --- Stefano Baroni - SISSA & DEMOCRITOS National Simulation Center - Trieste [+39] 040 3787 406 (tel) -528 (fax) / stefanobaroni (skype) La morale est une logique de l'action comme la logique est une morale de la pens?e - Jean Piaget Please, if possible, don't send me MS Word or PowerPoint attachments Why? See: http://www.gnu.org/philosophy/no-word-attachments.html -------------- next part -------------- An HTML attachment was scrubbed... URL: http://www.democritos.it/pipermail/pw_forum/attachments/20080410/85b4f642/attachment.htm
