<p>Dear pwscf users:<br>&nbsp;&nbsp;&nbsp; I am trying to use the v3.2.3 to study Raman properties. In the &quot;INPUT_PH&quot;, <br>it mentions literature &quot;M. Lazzeri and F. Mauri, Phys. Rev. Lett. 90, 036401 (2003)&quot;.<br>In this paper, Lazzeri computed the raman intensities form the second order derivative<br>
of the electronic density matrix with respect to a uniform electric field. Moreover, <br>another paper &quot;J. Sun and et. al. , Phys. Rev. B 74, 193101&quot; discussed the raman spectrum <br>via Lazzeri&#39;s method. I think this method, which has been implented int pwscf code, is <br>
different from the linear response method (1.S. Baroni et al. Rev. Mod. Phys. 73, 515 (2001),<br>2.P. Umari and et. al. Phys. Rev. B 63, 094305), which get the raman intensities by <br>polarizability tensors.<br>&nbsp;&nbsp;&nbsp; However, I cannot understant the results of example15(AlAs)which study the raman tensor of AlAs.<br>
The results of Raman tensor are showed as below:<br>=======================================================================================================<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Raman tensor (A^2)</p>
<p>&nbsp;&nbsp;&nbsp;&nbsp; atom #&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; pol.&nbsp; 1<br>&nbsp;&nbsp;&nbsp;&nbsp; -0.201910099100E-14&nbsp;&nbsp;&nbsp;&nbsp; -0.540966680608E-15&nbsp;&nbsp;&nbsp;&nbsp; -0.190481225566E-15<br>&nbsp;&nbsp;&nbsp;&nbsp; -0.145527656333E-14&nbsp;&nbsp;&nbsp;&nbsp; -0.739067155197E-15&nbsp;&nbsp;&nbsp;&nbsp; -0.517081932644E+01<br>&nbsp;&nbsp;&nbsp;&nbsp; -0.556205178654E-15&nbsp;&nbsp;&nbsp;&nbsp; -0.517081932644E+01&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.906690633696E-15<br>
&nbsp;&nbsp;&nbsp;&nbsp; atom #&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; pol.&nbsp; 2<br>&nbsp;&nbsp;&nbsp;&nbsp; ...<br>&nbsp;&nbsp;&nbsp;&nbsp; atom #&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; pol.&nbsp; 3<br>&nbsp;&nbsp;&nbsp;&nbsp; ...<br>&nbsp;&nbsp;&nbsp;&nbsp; atom #&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; pol.&nbsp; 1<br>&nbsp;&nbsp;&nbsp;&nbsp; ...<br>&nbsp;&nbsp;&nbsp;&nbsp; atom #&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; pol.&nbsp; 2<br>&nbsp;&nbsp;&nbsp;&nbsp; ...<br>&nbsp;&nbsp;&nbsp;&nbsp; atom #&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; pol.&nbsp; 3<br>&nbsp;&nbsp;&nbsp;&nbsp; ...<br>=======================================================================================================<br>
My questions: what does&nbsp; A^2 mean? what does pol. * mean?</p>
<p>I have another question about the alas.dynG. <br>1.Do omega 1, 2, 3 stand for the acoustic mode?<br>2.How can I get the Raman intensities?(the last but not the least!)<br>&nbsp;**************************************************************************<br>
&nbsp;&nbsp;&nbsp;&nbsp; omega( 1) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.043635 [THz] =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.455501 [cm-1]<br>&nbsp;( -0.499893&nbsp; 0.000000&nbsp; 0.432764&nbsp; 0.000000 -0.250591&nbsp; 0.000000 )<br>&nbsp;( -0.499920&nbsp; 0.000000&nbsp; 0.432787&nbsp; 0.000000 -0.250604&nbsp; 0.000000 )<br>&nbsp;&nbsp;&nbsp;&nbsp; omega( 2) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.043635 [THz] =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.455501 [cm-1]<br>
&nbsp;(&nbsp; 0.189458&nbsp; 0.000000 -0.164016&nbsp; 0.000000 -0.661194&nbsp; 0.000000 )<br>&nbsp;(&nbsp; 0.189468&nbsp; 0.000000 -0.164025&nbsp; 0.000000 -0.661229&nbsp; 0.000000 )<br>&nbsp;&nbsp;&nbsp;&nbsp; omega( 3) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.043635 [THz] =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.455501 [cm-1]<br>&nbsp;( -0.462802&nbsp; 0.000000 -0.534591&nbsp; 0.000000&nbsp; 0.000000&nbsp; 0.000000 )<br>
&nbsp;( -0.462827&nbsp; 0.000000 -0.534620&nbsp; 0.000000&nbsp; 0.000000&nbsp; 0.000000 )<br>&nbsp;&nbsp;&nbsp;&nbsp; omega( 4) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 10.686462 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 356.464369 [cm-1]<br>&nbsp;( -0.105298&nbsp; 0.000000&nbsp; 0.234486&nbsp; 0.000000&nbsp; 0.253351&nbsp; 0.000000 )<br>&nbsp;(&nbsp; 0.272091&nbsp; 0.000000 -0.605913&nbsp; 0.000000 -0.654659&nbsp; 0.000000 )<br>
&nbsp;&nbsp;&nbsp;&nbsp; omega( 5) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 10.686462 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 356.464369 [cm-1]<br>&nbsp;(&nbsp; 0.329240&nbsp; 0.000000&nbsp; 0.147848&nbsp; 0.000000&nbsp; 0.000000&nbsp; 0.000000 )<br>&nbsp;( -0.850757&nbsp; 0.000000 -0.382040&nbsp; 0.000000&nbsp; 0.000000&nbsp; 0.000000 )<br>&nbsp;&nbsp;&nbsp;&nbsp; omega( 6) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 10.686462 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 356.464369 [cm-1]<br>
&nbsp;( -0.103785&nbsp; 0.000000&nbsp; 0.231118&nbsp; 0.000000 -0.257044&nbsp; 0.000000 )<br>&nbsp;(&nbsp; 0.268182&nbsp; 0.000000 -0.597208&nbsp; 0.000000&nbsp; 0.664201&nbsp; 0.000000 )<br>&nbsp;**************************************************************************</p><br clear="all">
<br>-- <br>Y. C. Cheng<br>Department of Phyics<br>Nanjing University<br>Nanjing 210093<br>P. R. China<br>Tel: 86-25-83592907<br>Email: <a href="mailto:yccheng.nju@gmail.com">yccheng.nju@gmail.com</a>