<div dir="ltr">i have run the calculation for Raman tensor and have obtained the output file as given below. But raman tensor given in this output file is not the same as we see in papers. so how to extract raman tensor for raman active modes from this output file.<div><br></div><div><br>     Program PHONON v.6.4.1 starts on 30Aug2020 at  6: 8:38 <br><br>     This program is part of the open-source Quantum ESPRESSO suite<br>     for quantum simulation of materials; please cite<br>         "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);<br>         "P. Giannozzi et al., J. Phys.:Condens. Matter 29 465901 (2017);<br>          URL <a href="http://www.quantum-espresso.org">http://www.quantum-espresso.org</a>", <br>     in publications or presentations arising from this work. More details at<br>     <a href="http://www.quantum-espresso.org/quote">http://www.quantum-espresso.org/quote</a><br><br>     Parallel version (MPI & OpenMP), running on      40 processor cores<br>     Number of MPI processes:                40<br>     Threads/MPI process:                     1<br><br>     MPI processes distributed on     2 nodes<br>     R & G space division:  proc/nbgrp/npool/nimage =      40<br><br>     Reading data from directory:<br>     ./MoS2.save/<br><br>     IMPORTANT: XC functional enforced from input :<br>     Exchange-correlation      = PW ( 1  4  0  0 0 0)<br>     Any further DFT definition will be discarded<br>     Please, verify this is what you really want<br><br> <br>     Parallelization info<br>     --------------------<br>     sticks:   dense  smooth     PW     G-vecs:    dense   smooth      PW<br>     Min          16      16      4                 3658     3658     560<br>     Max          17      17      5                 3735     3735     625<br>     Sum         673     673    199               148251   148251   24293<br> <br>        1 /   1 q-points for this run, from  1 to  1:<br>       N       xq(1)         xq(2)         xq(3) <br>       1   0.000000000   0.000000000   0.000000000<br> <br><br>     Calculation of q =    0.0000000   0.0000000   0.0000000<br><br>     Restart after Phonon calculation<br><br>     phonons of MoS2 at Gamma                                                   <br><br>     bravais-lattice index     =            4<br>     lattice parameter (alat)  =       5.9337  a.u.<br>     unit-cell volume          =    1874.4360 (a.u.)^3<br>     number of atoms/cell      =            3<br>     number of atomic types    =            2<br>     kinetic-energy cut-off    =      70.0000  Ry<br>     charge density cut-off    =     280.0000  Ry<br>     convergence threshold     =      1.0E-12<br>     beta                      =       0.7000<br>     number of iterations used =            4<br>     Exchange-correlation      = PW ( 1  4  0  0 0 0)<br><br><br>     celldm(1)=    5.93374  celldm(2)=    0.00000  celldm(3)=   10.35987<br>     celldm(4)=    0.00000  celldm(5)=    0.00000  celldm(6)=    0.00000<br><br>     crystal axes: (cart. coord. in units of alat)<br>               a(1) = (  1.0000  0.0000  0.0000 )  <br>               a(2) = ( -0.5000  0.8660  0.0000 )  <br>               a(3) = (  0.0000  0.0000 10.3599 )  <br><br>     reciprocal axes: (cart. coord. in units 2 pi/alat)<br>               b(1) = (  1.0000  0.5774  0.0000 )  <br>               b(2) = (  0.0000  1.1547  0.0000 )  <br>               b(3) = (  0.0000  0.0000  0.0965 )  <br><br><br>     Atoms inside the unit cell: <br><br>     Cartesian axes<br><br>     site n.  atom      mass           positions (alat units)<br>        1     Mo  95.9500   tau(    1) = (   -0.00005    0.57738    4.18232  )<br>        2     S   32.0600   tau(    2) = (    0.50005    0.28865    4.68113  )<br>        3     S   32.0600   tau(    3) = (    0.50005    0.28865    3.68351  )<br><br>     Computing dynamical matrix for <br>                    q = (   0.0000000   0.0000000   0.0000000 )<br> <br>      3 Sym.Ops. (with q -> -q+G )<br><br><br>     G cutoff =  249.7212  (   3660 G-vectors)     FFT grid: ( 32, 32,360)<br>     number of k points=    14<br>                       cart. coord. in units 2pi/alat<br>        k(    1) = (   0.1250000   0.2165064  -0.0482631), wk =   0.0833333<br>        k(    2) = (   0.1250000   0.5051815  -0.0482631), wk =   0.1666667<br>        k(    3) = (   0.1250000  -0.3608439  -0.0482631), wk =   0.1666667<br>        k(    4) = (   0.1250000  -0.0721688  -0.0482631), wk =   0.0833333<br>        k(    5) = (   0.3750000   0.6495191  -0.0482631), wk =   0.0833333<br>        k(    6) = (   0.3750000  -0.2165064  -0.0482631), wk =   0.0833333<br>        k(    7) = (  -0.1250000   0.2165064   0.0482631), wk =   0.1666667<br>        k(    8) = (  -0.1250000   0.5051815   0.0482631), wk =   0.1666667<br>        k(    9) = (  -0.3750000   0.3608439  -0.0482631), wk =   0.1666667<br>        k(   10) = (  -0.1250000  -0.3608439   0.0482631), wk =   0.1666667<br>        k(   11) = (   0.3750000  -0.0721688  -0.0482631), wk =   0.1666667<br>        k(   12) = (  -0.1250000  -0.0721688   0.0482631), wk =   0.1666667<br>        k(   13) = (  -0.3750000   0.6495191   0.0482631), wk =   0.1666667<br>        k(   14) = (  -0.3750000  -0.2165064   0.0482631), wk =   0.1666667<br><br>     PseudoPot. # 1 for Mo read from file:<br>     /home/ankitasingh.rs.mst18.iitbhu/Tensor/Mo.upf<br>     MD5 check sum: 757ed1ca8eaf384e7a6ffdf32476f2ff<br>     Pseudo is Norm-conserving + core correction, Zval = 14.0<br>     Generated using ONCVPSP code by D. R. Hamann<br>     Using radial grid of 1602 points,  6 beta functions with: <br>                l(1) =   0<br>                l(2) =   0<br>                l(3) =   1<br>                l(4) =   1<br>                l(5) =   2<br>                l(6) =   2<br><br>     PseudoPot. # 2 for S  read from file:<br>     /home/ankitasingh.rs.mst18.iitbhu/Tensor/S.upf<br>     MD5 check sum: cb156da9049275a86b4f507c9016d0ed<br>     Pseudo is Norm-conserving + core correction, Zval =  6.0<br>     Generated using ONCVPSP code by D. R. Hamann<br>     Using radial grid of 1132 points,  6 beta functions with: <br>                l(1) =   0<br>                l(2) =   0<br>                l(3) =   1<br>                l(4) =   1<br>                l(5) =   2<br>                l(6) =   2<br><br>     Mode symmetry, C_s (m)     point group:<br><br><br>     Electric field:<br>     Dielectric constant<br>     Born effective charges in two ways <br>     Raman tensor<br>     Electro-optic tensor<br><br><br>     Atomic displacements:<br>     There are   9 irreducible representations<br><br>     Representation     1      1 modes -A'  Done<br><br>     Representation     2      1 modes -A'  Done<br><br>     Representation     3      1 modes -A'  Done<br><br>     Representation     4      1 modes -A'  Done<br><br>     Representation     5      1 modes -A'  Done<br><br>     Representation     6      1 modes -A'  Done<br><br>     Representation     7      1 modes -A''  Done<br><br>     Representation     8      1 modes -A''  Done<br><br>     Representation     9      1 modes -A''  Done<br><br><br>     PHONON       :      1.89s CPU      6.28s WALL<br><br><br>          Dielectric constant in cartesian axis <br><br>          (       3.790423733       0.000809259      -0.000223671 )<br>          (       0.000809259       3.791358185       0.000129136 )<br>          (      -0.000223671       0.000129136       1.188378071 )<br><br>          Effective charges (d Force / dE) in cartesian axis<br><br>           atom      1   Mo <br>      Ex  (       -1.21765        0.00046        0.00111 )<br>      Ey  (        0.00046       -1.21712       -0.00064 )<br>      Ez  (       -0.00000        0.00000       -0.08891 )<br>           atom      2   S  <br>      Ex  (        0.48972       -0.00025       -0.00004 )<br>      Ey  (       -0.00025        0.48943        0.00002 )<br>      Ez  (        0.00001       -0.00001        0.04238 )<br>           atom      3   S  <br>      Ex  (        0.50559       -0.00025        0.00051 )<br>      Ey  (       -0.00025        0.50530       -0.00030 )<br>      Ez  (       -0.00001        0.00001        0.05387 )<br><br>              Electro-optic tensor is defined as <br>            the derivative of the dielectric tensor <br>              with respect to one electric field <br>                 units are Rydberg a.u. <br><br>            to obtain the static chi^2 multiply by 1/2  <br><br>            to convert to pm/Volt multiply per 2.7502  <br><br><br>          Electro-optic tensor in cartesian axis: <br><br>          (      -4.417946244      -7.617122054       1.871368586 )<br>          (      -7.617122054      -1.472648748       3.040246832 )<br>          (       1.871368586       3.040246832      -0.037270747 )<br><br>          (      -7.617122054      -1.472648748       3.040246832 )<br>          (      -1.472648748      11.018058658       5.381943240 )<br>          (       3.040246832       5.381943240       0.021518276 )<br><br>          (       1.871368586       3.040246832      -0.037270747 )<br>          (       3.040246832       5.381943240       0.021518276 )<br>          (      -0.037270747       0.021518276       0.016833272 )<br><br><br>     Raman tensor (A^2)<br><br>     atom #    1    pol.  1<br>     -0.329668500707E+00     -0.142014670420E+02      0.257292263558E-01<br>     -0.142014670420E+02     -0.341151793663E+00      0.700276792529E-01<br>      0.257292263558E-01      0.700276792529E-01     -0.569106019847E-03<br>     atom #    1    pol.  2<br>     -0.140011880210E+02      0.574164647847E-02      0.700276792529E-01<br>      0.574164647847E-02      0.143884862985E+02      0.106590225291E+00<br>      0.700276792529E-01      0.106590225291E+00      0.328573513756E-03<br>     atom #    1    pol.  3<br>     -0.460515743770E+00     -0.101466465886E+01     -0.222065709668E-02<br>     -0.101466465886E+01     -0.163214957163E+01      0.128209697255E-02<br>     -0.222065709668E-02      0.128209697255E-02      0.130524522359E-01<br>     atom #    2    pol.  1<br>      0.140581302092E+01      0.707398985511E+01      0.930516667954E-01<br>      0.707398985511E+01      0.139791465548E+01     -0.724986307385E-02<br>      0.930516667954E-01     -0.724986307385E-02      0.774338492435E-02<br>     atom #    2    pol.  2<br>      0.626918351417E+01      0.394918272009E-02     -0.724986307385E-02<br>      0.394918272009E-02     -0.788791644287E+01      0.846802460008E-01<br>     -0.724986307385E-02      0.846802460008E-01     -0.447064537051E-02<br>     atom #    2    pol.  3<br>      0.121000444708E+02      0.567513288288E+00      0.715323559367E-02<br>      0.567513288288E+00      0.127553523703E+02     -0.412992249558E-02<br>      0.715323559367E-02     -0.412992249558E-02      0.120848475018E+01<br>     atom #    3    pol.  1<br>     -0.125370383234E+01      0.404592247444E+01     -0.152464706768E+00<br>      0.404592247444E+01     -0.126543444413E+01     -0.375870307386E-01<br>     -0.152464706768E+00     -0.375870307386E-01      0.787368048152E-05<br>     atom #    3    pol.  2<br>      0.477990772733E+01      0.586530589243E-02     -0.375870307386E-01<br>      0.586530589243E-02     -0.332548256529E+01     -0.195866471398E+00<br>     -0.375870307386E-01     -0.195866471398E+00     -0.454587154552E-05<br>     atom #    3    pol.  3<br>     -0.114151241036E+02      0.601448588700E+00      0.396186167582E-02<br>      0.601448588700E+00     -0.107206310944E+02     -0.228738190503E-02<br>      0.396186167582E-02     -0.228738190503E-02     -0.113597287830E+01<br> <br>     Number of q in the star =    1<br>     List of q in the star:<br>          1   0.000000000   0.000000000   0.000000000<br><br>          Dielectric constant in cartesian axis <br><br>          (       3.790423733       0.000809259      -0.000223671 )<br>          (       0.000809259       3.791358185       0.000129136 )<br>          (      -0.000223671       0.000129136       1.188378071 )<br><br>          Effective charges (d Force / dE) in cartesian axis<br><br>           atom      1   Mo <br>      Ex  (       -1.21765        0.00046        0.00111 )<br>      Ey  (        0.00046       -1.21712       -0.00064 )<br>      Ez  (       -0.00000        0.00000       -0.08891 )<br>           atom      2   S  <br>      Ex  (        0.48972       -0.00025       -0.00004 )<br>      Ey  (       -0.00025        0.48943        0.00002 )<br>      Ez  (        0.00001       -0.00001        0.04238 )<br>           atom      3   S  <br>      Ex  (        0.50559       -0.00025        0.00051 )<br>      Ey  (       -0.00025        0.50530       -0.00030 )<br>      Ez  (       -0.00001        0.00001        0.05387 )<br><br>          Effective charges (d P / du) in cartesian axis <br><br>           atom      1   Mo <br>      Px  (       -1.21566       -0.00141        0.00008 )<br>      Py  (       -0.00141       -1.21728       -0.00005 )<br>      Pz  (        0.00033       -0.00019       -0.09196 )<br>           atom      2   S  <br>      Px  (        0.48710       -0.00525       -0.00018 )<br>      Py  (       -0.00525        0.48104        0.00010 )<br>      Pz  (        0.00134       -0.00078        0.04350 )<br>           atom      3   S  <br>      Px  (        0.50354       -0.00552        0.00021 )<br>      Py  (       -0.00552        0.49717       -0.00012 )<br>      Pz  (        0.00168       -0.00097        0.05486 )<br><br>     Diagonalizing the dynamical matrix<br><br>     q = (    0.000000000   0.000000000   0.000000000 ) <br><br> **************************************************************************<br>     freq (    1) =      -2.829149 [THz] =     -94.370238 [cm-1]<br>     freq (    2) =       0.120485 [THz] =       4.018943 [cm-1]<br>     freq (    3) =       0.994796 [THz] =      33.182814 [cm-1]<br>     freq (    4) =       8.096315 [THz] =     270.064006 [cm-1]<br>     freq (    5) =       8.426042 [THz] =     281.062493 [cm-1]<br>     freq (    6) =      11.370676 [THz] =     379.284933 [cm-1]<br>     freq (    7) =      11.749743 [THz] =     391.929246 [cm-1]<br>     freq (    8) =      11.973769 [THz] =     399.401955 [cm-1]<br>     freq (    9) =      13.954583 [THz] =     465.474780 [cm-1]<br> **************************************************************************<br><br>     Mode symmetry, C_s (m)     point group:<br><br>     freq (  1 -  1) =        -94.4  [cm-1]   --> A''             I+R<br>     freq (  2 -  2) =          4.0  [cm-1]   --> A'              I+R<br>     freq (  3 -  3) =         33.2  [cm-1]   --> A'              I+R<br>     freq (  4 -  4) =        270.1  [cm-1]   --> A''             I+R<br>     freq (  5 -  5) =        281.1  [cm-1]   --> A'              I+R<br>     freq (  6 -  6) =        379.3  [cm-1]   --> A'              I+R<br>     freq (  7 -  7) =        391.9  [cm-1]   --> A''             I+R<br>     freq (  8 -  8) =        399.4  [cm-1]   --> A'              I+R<br>     freq (  9 -  9) =        465.5  [cm-1]   --> A'              I+R<br> <br>     PHONON       :      1.91s CPU      6.71s WALL<br><br>     INITIALIZATION: <br>     phq_setup    :      0.01s CPU      0.21s WALL (       1 calls)<br> <br>     init_vloc    :      0.01s CPU      0.02s WALL (       1 calls)<br>     init_us_1    :      0.01s CPU      0.04s WALL (       1 calls)<br> <br>     DIELECTRIC CONSTANT AND EFFECTIVE CHARGES:<br> <br>     RAMAN COEFFICIENTS, THIRD-ORDER CHI:<br> <br>     DYNAMICAL MATRIX:<br>     phqscf       :      0.00s CPU      0.02s WALL (       1 calls)<br>     dynmatrix    :      0.01s CPU      0.18s WALL (       1 calls)<br> <br>     phqscf       :      0.00s CPU      0.02s WALL (       1 calls)<br> <br> <br> <br>     phqscf       :      0.00s CPU      0.02s WALL (       1 calls)<br> <br> <br> <br> <br> <br> <br> <br> <br>      General routines<br>     fft          :      0.00s CPU      0.22s WALL (       5 calls)<br>     davcio       :      0.00s CPU      0.13s WALL (      14 calls)<br> <br> <br>     PHONON       :      1.91s CPU      6.71s WALL<br><br> <br>   This run was terminated on:   6: 8:45  30Aug2020            <br><br>=------------------------------------------------------------------------------=<br>   JOB DONE.<br>=------------------------------------------------------------------------------=<br></div></div>