<div dir="ltr"><div><div><div>Dear Users,<br><br></div>I am trying to perform phonon calculations for BaTiO3. After performing the the calculations at gamma point, with the help of 'ph.x' executable, I get following results for phonon frequencies,<br><br>------------------------------------------------------------------<br>Mode symmetry, C_4v (4mm)  point group:<br><br>     freq (  1 -  2) =       -190.6  [cm-1]   --> E    G_5 D_5    I+R<br>     freq (  3 -  3) =        -45.6  [cm-1]   --> A_1  G_1 D_1    I+R<br>     freq (  4 -  5) =         26.2  [cm-1]   --> E    G_5 D_5    I+R<br>     freq (  6 -  6) =        115.4  [cm-1]   --> A_1  G_1 D_1    I+R<br>     freq (  7 -  8) =        146.5  [cm-1]   --> E    G_5 D_5    I+R<br>     freq (  9 -  9) =        245.7  [cm-1]   --> B_1  G_3 D_2    R  <br>     freq ( 10 - 10) =        268.0  [cm-1]   --> A_1  G_1 D_1    I+R<br>     freq ( 11 - 12) =        277.1  [cm-1]   --> E    G_5 D_5    I+R<br>     freq ( 13 - 14) =        457.1  [cm-1]   --> E    G_5 D_5    I+R<br>     freq ( 15 - 15) =        496.7  [cm-1]   --> A_1  G_1 D_1    I+R<br>---------------------------------------------------------------------<br><br></div>I could not set lraman=.true. as I was using ultra-soft pseudopotentials. <br><br></div>Following were the results for the phonon frequencies upon the application of asr=simple with 'dynmat.x' executable. <br><br>--------------------------------------------------------------------<br> <br><div>     IR activities are in (D/A)^2/amu units<br><br># mode   [cm-1]    [THz]      IR<br>    1   -189.03   -5.6670   78.9556<br>    2   -189.03   -5.6670   78.9556<br>    3      0.00    0.0000    0.0000<br>    4      0.00    0.0000    0.0000<br>    5      0.00    0.0000    0.0000<br>    6    129.14    3.8716    8.7326<br>    7    140.75    4.2197    0.7690<br>    8    140.75    4.2197    0.7690<br>    9    260.86    7.8203    0.0292<br>   10    260.86    7.8203    0.0292<br>   11    263.08    7.8868    0.0000<br>   12    273.06    8.1862   42.5449<br>   13    456.98   13.6999    2.7989<br>   14    456.98   13.6999    2.7989<br>   15    502.29   15.0582   12.3293<br><br>-----------------------------------------------------------------------<br><div><br></div><div>These results are only for IR-modes. So my questions are as follows,<br><br>1. Does the "I+R" tag at the end of the output file given by ph.x calculations, specify that the modes are both IR as well as Raman active?<br><br>2. Which are the correct values of IR frequencies? Those given by ph.x calculations or by dynmat.x calculations?<br><br></div><div>Thanks in advance.<br><br></div><div>Best,<br></div><div>DS.<br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div>Deepashri Saraf<br></div><div>PhD Student<br></div><div>Department of Physics, SPPU<br></div><div>Pune (India)<br></div></div></div>
</div></div></div>