<div dir="ltr"><div><div><div><div><div>Dear developers of QE, <br><br></div>I noticed a backward compatibility issue on the electron-phonon matrix element between QE 4.0.3 and 5.1.2. <br><br></div>I was wondering if it was a bug and if not, what was the reason.<br><br></div><div>I added the following line in the two versions of the code in PHonon/PH/elphon.f90<br><br>if (ik==215) then<br>write(*,*)'elphmat(:,:,:)**2',SUM((REAL(REAL(elphmat(:,:,:))))**2)+SUM((REAL(AIMAG(elphmat(:,:,:))))**2)<br></div><div><br></div><div>This should be a physically relevant quantity that should not change between version (the order of the mode changed but it is not an issue).<br></div><div><br></div>Test: <br></div>6x6x6 q-grid same pseudo, tol and inputs parameters <br><br></div>There is 16 irr q-points:<br><br><div>     (  16q-points):<br>       N         xq(1)         xq(2)         xq(3)<br>       1   0.000000000   0.000000000   0.000000000<br>       2  -0.166666667   0.166666667  -0.166666667<br>       3  -0.333333333   0.333333333  -0.333333333<br>       4   0.500000000  -0.500000000   0.500000000<br>       5   0.000000000   0.333333333   0.000000000<br>       6  -0.166666667   0.500000000  -0.166666667<br>       7   0.666666667  -0.333333333   0.666666667<br>       8   0.500000000  -0.166666667   0.500000000<br>       9   0.333333333   0.000000000   0.333333333<br>      10   0.000000000   0.666666667   0.000000000<br>      11   0.833333333  -0.166666667   0.833333333<br>      12   0.666666667  -0.000000000   0.666666667<br>      13   0.000000000  -1.000000000   0.000000000<br>      14   0.666666667  -0.333333333   1.000000000<br>      15   0.500000000  -0.166666667   0.833333333<br>      16  -0.333333333  -1.000000000   0.000000000<br><br></div><div>All the elphmat**2 are the same (up to 4 digits) for 13 of the 16 q-points. For example<br></div><div><br>QE 5.1.2  q = (   -0.166666667   0.166666667  -0.166666667 )</div><div>elphmat(:,:,:)**2  0.5446267<br>elphmat(:,:,:)**2  0.5446253<br>elphmat(:,:,:)**2  0.9506719<br>elphmat(:,:,:)**2  0.9506719<br><br>     freq (    1) =       9.137536 [THz] =     304.795384 [cm-1]<br>     freq (    2) =       9.137536 [THz] =     304.795384 [cm-1]<br>     freq (    3) =      15.052865 [THz] =     502.109513 [cm-1]<br>     freq (    4) =      39.653912 [THz] =    1322.712114 [cm-1]<br>     freq (    5) =      39.653912 [THz] =    1322.712114 [cm-1]<br>     freq (    6) =      40.741956 [THz] =    1359.005377 [cm-1]<br><br>QE 4.0.3 q = (   -0.166666667   0.166666667  -0.166666667 )<br>elphmat(:,:,:)**2  0.5446267<br>elphmat(:,:,:)**2  0.9506719<br>elphmat(:,:,:)**2  0.5446266<br>elphmat(:,:,:)**2  0.9506719<br><br>     omega( 1) =       9.137611 [THz] =     304.799915 [cm-1]<br>     omega( 2) =       9.137611 [THz] =     304.799915 [cm-1]<br>     omega( 3) =      15.052641 [THz] =     502.105399 [cm-1]<br>     omega( 4) =      39.653770 [THz] =    1322.716173 [cm-1]<br>     omega( 5) =      39.653770 [THz] =    1322.716173 [cm-1]<br>     omega( 6) =      40.741800 [THz] =    1359.009209 [cm-1]<br><br><br></div><div>But the three irr-q<b> 5, 10</b> and <b>16</b> are NOT the same:<br>QE 5.1.2 q = (    0.000000000   0.333333333   0.000000000 )<br>elphmat(:,:,:)**2  0.1192376<br>elphmat(:,:,:)**2  0.9861287<br>elphmat(:,:,:)**2   1.348592<br>elphmat(:,:,:)**2  0.3134981<br><br></div><div>     freq (    1) =      11.633039 [THz] =     388.036416 [cm-1]<br>     freq (    2) =      11.633039 [THz] =     388.036416 [cm-1]<br>     freq (    3) =      16.188917 [THz] =     540.004156 [cm-1]<br>     freq (    4) =      39.110566 [THz] =    1304.588071 [cm-1]<br>     freq (    5) =      39.110566 [THz] =    1304.588071 [cm-1]<br>     freq (    6) =      40.947507 [THz] =    1365.861806 [cm-1]<br><br></div><div>QE 4.0.3 q = (    0.000000000   0.333333333   0.000000000 <br></div><div>elphmat(:,:,:)**2  0.3075966<br>elphmat(:,:,:)**2   1.354492<br>elphmat(:,:,:)**2  0.1192375<br>elphmat(:,:,:)**2  0.9861288<br><br>     omega( 1) =      11.633152 [THz] =     388.042751 [cm-1]<br>     omega( 2) =      11.633152 [THz] =     388.042751 [cm-1]<br>     omega( 3) =      16.188922 [THz] =     540.007888 [cm-1]<br>     omega( 4) =      39.110425 [THz] =    1304.592025 [cm-1]<br>     omega( 5) =      39.110425 [THz] =    1304.592025 [cm-1]<br>     omega( 6) =      40.947359 [THz] =    1365.865939 [cm-1]<br><br></div><div>As you can see the difference for the 0,1/3,0 is not huge but significant. As a matter of fact when<br></div><div>I'm using the dvscf produced by QE 4 or 5 for that irr-q it can lead to physically very different results. Therefore I would like to know which dvscf file is correct (QE 4 or 5)?<br><br></div><div>Thank you, <br><br></div><div>Best Regards, <br><br></div><div>Samuel Ponce <br>Department of Materials, University of Oxford<br></div><div><br><br><br></div></div>