<div dir="ltr"><div>Dear colleagues, <br><br></div>I have found strange differences in results for the calculation of pressure using pwscf in versions 4.3 and 5.0.N (N=1,2).  We are calculating wurtzite CdSe, using the following UPF files. <br>
<br clear="all"><div><div>PAW calculation<br>   Cd  112.40000  Cd.pbe-dn-kjpaw_psl.0.2.UPF <br>   Se   78.96000  Se.pbe-n-kjpaw_psl.0.2.UPF <br></div><div><br><br>Ultrasoft PP calculation<br>   Cd  112.40000  Cd.pbe-dn-rrkjus_psl.0.2.UPF <br>
   Se   78.96000  Se.pbe-n-rrkjus_psl.0.2.UPF <br><br><br><br></div><div>For the PAW calculations, the convergence of pressure with the density cutoff (euctrho) shows jumps  (see the attached  figure paw-press.eps) and no convergence in the case of version 5.0.N, while there is convergence with version 4.3. For both versions the total energies are the same. <br>
</div><div>The jumps are associated with differences in the number of symmetry operations detected. The message <br><br> "         (note:  6 additional sym.ops. were found but ignored<br>           their fractional translations are incommensurate with FFT grid)"<br>
</div><div>is associated with higher pressures (33 kbar). Lower pressures ~ 26 kbar) are found when the FFT grid is commensurate and all 12 sym. ops. are used. <br>In the calculations with version 4.3 the pressures are somewhat lower, but apparently convergent. Also, there are less symmetry operations detected or used than for version 5.0.N. <br>
</div><div><br></div><div>In the ultrasoft calculations, (see figure uspp-press), both versions 4.3 and 5.0.N produce convergent pressures vs ecutrho, but the pressures differ almost by 10 kbar. (The total energies are almost equal, as expected). <br>
</div><div><br></div><div>Thanks for any comment. <br></div><div><br>-- <br>
<div>Eduardo Menendez Proupin</div><div>Departamento de Fisica, Facultad de Ciencias, Universidad de Chile <br>URL: <a href="http://fisica.ciencias.uchile.cl/%7Eemenendez" target="_blank">http://fisica.ciencias.uchile.cl/~emenendez</a><div>
<span style="font-family:sans-serif;font-size:13px;line-height:19px"><font><span style="font-family:sans-serif;line-height:19px"><i><span style="font-family:"Arial","sans-serif"">"Lograremos el anhelado desarrollo cuando sea la ciudadanía, y no los científicos,
la que pida más ciencia para Chile." J. Babul</span></i></span></font><br></span></div></div>
</div></div></div>