<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Dear all,<br>
<br>
I have few questions related to the above mentioned PP's.<br>
<br>
(1) After generating the PAW and USPP versions of these PP (pslib1.0.0) using ld1 code provided with q.e.5.1.2, I applied the test script. The output seems to be fine in the case of USPP but for PAW extra few lines are printed. For example, in the case of
LDA scalar-rel Fe, I get:<br>
<br>
Making Fe.pz-n-kjpaw_psl.0.2.4.test.in ... Done !<br>
3 2 3D 1( 7.00) -0.41920 -0.41398 -0.00522 !<br>
Making Fe.pz-n-kjpaw_psl.1.0.0.test.in ... Done<br>
Making Fe.pz-n-rrkjus_psl.0.2.4.test.in ... Done<br>
Making Fe.pz-n-rrkjus_psl.1.0.0.test.in ... Done<br>
Making Fe.pz-spn-kjpaw_psl.1.0.0.test.in ... Done<br>
Making Fe.pz-spn-rrkjus_psl.1.0.0.test.in ... Done<br>
<br>
And in the case of Fe-rel I get:<br>
<br>
Making Fe.rel-pz-n-kjpaw_psl.0.2.4.test.in ... Done !<br>
3 2 1.5 3D 1( 4.00) -0.42153 -0.41623 -0.00530 !<br>
3 2 2.5 3D 1( 3.00) -0.41137 -0.40618 -0.00519 !<br>
Making Fe.rel-pz-n-kjpaw_psl.1.0.0.test.in ... Done<br>
Making Fe.rel-pz-n-rrkjus_psl.0.2.4.test.in ... Done<br>
Making Fe.rel-pz-n-rrkjus_psl.1.0.0.test.in ... Done<br>
Making Fe.rel-pz-spn-kjpaw_psl.1.0.0.test.in ... Done<br>
Making Fe.rel-pz-spn-rrkjus_psl.1.0.0.test.in ... Done<br>
<br>
Similar output obtained for Co PAW. Are these extra printed lines indicative of some error in the generation?
<br>
<br>
<br>
(2) The charge density cutoff recommended (in the UPF files) for these PP is about 12
<a name="ecutrho">*ecutwfc</a>. This large cutoff for the charge density was also employed in a recent paper using Fe (pslib0.31) (http://arxiv.org/pdf/1502.01534.pdf) . However, when I did convergence test for the total energy on a unit cell of SrFeO3 , and
SrCO3, and it looks like 4* <a name="idp81536"></a><a name="ecutrho">ecutwfc</a> is sufficient to achieve convergence. The convergence in the case of SrFeO3 (pslib0.3, USPP) can be found here (http://web.mit.edu/myoussef/www/Fe_16e_pslib.pdf). The pink bound
indicates that convergence was achieved within 25 mev/unit cell (5 meV/atom) at
<a name="ecutrho">ecutwfc=85 Ry and </a><a name="idp81536"></a><a name="ecutrho">ecutrho =4 *</a><a name="ecutrho">ecutwfc</a>. It seems to me that the recommend
<a name="ecutrho">ecutrho it too strict. </a> Are there other tests that one needs to perform to check the convergence of charge density cutoff?
<br>
<br>
(3) The recommended <a name="ecutrho">ecutwfc and </a><a name="ecutrho">ecutrho for this set of PP in Pslib do not change from PAW to USPP and do not depend on the functional. I did preliminary tests </a> and it seems so far that the values work for GGA also
work for LDA consistent with the recommendation. I wonder if one could test convergence of
<a name="ecutrho">ecutwfc and </a><a name="ecutrho">ecutrho for LDA only and apply the outcome to all other semilocal functionals, or are there cases where one encounters different convergence depending on the functional.<br>
<br>
<br>
Thank you in advance for reading and help!<br>
Mostafa Youssef<br>
MIT<br>
</a><br>
<br>
</div>
</body>
</html>