<div dir="ltr">Dear researchers<div><br></div><div>I have been trying to calculate the dielectric constant of multilayer TMDC materials. I have used the ph.x package of quantum espresso for this. From simulation I have got the following values:</div><div>         </div><div>  Dielectric constant in cartesian axis </div><div><br></div><div>          (       6.990898310       0.000000000       0.000000000 )</div><div>          (       0.000000000       6.990898310       0.000000000 )</div><div>          (       0.000000000       0.000000000       1.520000669 )</div><div><br></div><div>The values were obtained for Mos2 trilayer. I did the scf calculation before using the phonon package. The following code was used for phonon:</div><div><br></div><div><div>&inputph</div><div>  tr2_ph=1.0d-14,</div><div>  prefix='MoS2AAA187',</div><div>  epsil=.true.</div><div>  outdir='D:\Quantum ESPRESSO 64-bit 5.1.2-mpich2\MoS2AAA187\outdir',</div><div>  !outdir='./'</div><div>  fildyn='MoS2AAA187.dynG',</div><div>  </div><div> /</div><div> 0.0 0.0 0.0</div></div><div><br></div><div>For scf calculation I used pbe-hgh pseudopotentials. </div><div>The question is, why the z axis value of dielectric constant is so low? Values found in the literature is typically close. Is there something that must be added in phonon package?</div><div><br></div><div>Thanks in advance.</div><div><br></div><div>Sincerely yours</div><div>Kanak</div><div>EEE, BUET</div></div>