<div dir="ltr"><div><div><div><div><div><div><div>Dear QE users and developers <br></div>I try to obtain Raman intensity using ph.x and dynmat.x. I run example 15 in the phonon examples (espresso 5.4.0) and I get only IR intensity. The plot of gnuplot give the Raman intensity but it is just the frequency in Thz and not really the intensity:<br><br>plot 'plotdata_zno.dat' u (\$2):(\$3) title ' ZnO-RAMAN' w i lw 2<br><br></div>the file 'plotdata_zno.dat' contain:<br>   1     -0.00   -0.0000    0.0000<br>    2     -0.00   -0.0000    0.0000<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      0.00    0.0000    0.0000<br>    7     50.57    1.5161    0.0000<br>    8    107.07    3.2100    0.0000<br>    9    158.38    4.7482    0.0000<br>   10    215.37    6.4565    0.0000<br>   11    276.67    8.2943    0.0000<br>   12    337.61   10.1212    0.0000<br></div>(\$2):(\$3) This is the frequency in cm^-1 vs the frequency in THz. I do not understand where is the real intensity. <br></div>After that I try to calculate the intensity with the following input:<br></div>For ph.x:<br>phonon  for NaCl<br> &inputph<br>  tr2_ph=1.0d-15,<br>  prefix='NaCl',<br>  fildvscf='aldv',<br>  amass(1)=22.9898,<br>  amass(2)=35.453,<br>  outdir='/home/argaman/tmp/tmp17116/',<br>  fildyn='NaCl.dyn',<br>  trans=.true.,<br>  lraman=.true.,<br>  epsil=.true.,<br> /<br>0 0 0<br><br></div>For dynmat.x:<br> &input fildyn='NaCl.dyn', asr='simple'/<br><br></div>With norm-conserving PP and LDA. I still do not obtain the intensity for Raman in the file NaCl.dm.out, only IR:<br><br>     Reading Dynamical Matrix from file NaCl.dyn<br>     ...Force constants read<br>     ...epsilon and Z* read<br>     ...Raman cross sections read<br>     A direction for q was not specified:TO-LO splitting will be absent<br><br>     Polarizability (A^3 units)<br>     multiply by 0.655304 for Clausius-Mossotti correction<br>         5.147271    0.000000    0.000000<br>         0.000000    5.147271    0.000000<br>         0.000000    0.000000    5.147271<br><br>     IR activities are in (D/A)^2/amu units<br><br># mode   [cm-1]    [THz]      IR<br>    1      0.00    0.0000    0.0000<br>    2      0.00    0.0000    0.0000<br>    3      0.00    0.0000    0.0000<br>    4    183.55    5.5027    1.8492<br>    5    183.55    5.5027    1.8492<br>    6    183.55    5.5027    1.8492<br> <br>Electronic dielectric permittivity tensor (F/m units)<br>         2.578030    0.000000    0.000000<br>         0.000000    2.578030    0.000000<br>         0.000000    0.000000    2.578030<br> <br> ... with zone-center polar mode contributions<br>         5.433926    0.000000    0.000000<br>         0.000000    5.433926    0.000000<br>         0.000000    0.000000    5.433926<br><br><div><div><div><div><div><div>What did I do wrong? <br><br></div><div>Uri Argaman<br></div><div>Ben-Gurion University<br></div><div>Israel<br></div></div></div></div></div></div></div>