<div dir="ltr"><span style="font-size:12.8000001907349px">I confused about choice of pseudo potential for system about this atoms (C,N,O,Fe)</span><br><div><span style="font-size:12.8000001907349px"><br></span></div><div><span style="font-size:12.8000001907349px">first I used ultra soft pseudo potential </span></div><div><div style="font-size:12.8000001907349px">C  12.0107     C.pbe-van_ak.UPF</div><div style="font-size:12.8000001907349px">Fe 55.8450     Fe.pbe-sp-van_ak.UPF</div><div style="font-size:12.8000001907349px">N  14.0067     N.pbe-van_ak.UPF</div></div><div style="font-size:12.8000001907349px"><img src="cid:ii_14f4bd1788fd5667" alt="Inline image 1" width="313" height="167" style="margin-right: 0px;"><br></div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">Ecut optimum:   ecutwfc   = 60.0,  ecutrho   = 600.0,</span><br></div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px"><br></span></div><div><span style="font-size:12.8000001907349px">This Ecut is too high for ultra soft pseudo potential.</span></div><div><span style="font-size:12.8000001907349px"><br></span></div><div><span style="font-size:12.8000001907349px">I use another </span><span style="font-size:12.8000001907349px">ultra soft pseudo potential</span></div><div><div style="font-size:12.8000001907349px"><div>C  12.0107     C.pbe-n-rrkjus_psl.0.1.UPF</div><div>Fe 55.8450     Fe.pbe-nd-rrkjus.UPF</div><div>N  14.0067     N.pbe-n-rrkjus_psl.0.1.UPF</div></div><div style="font-size:12.8000001907349px"><br></div><div style="font-size:12.8000001907349px">Ecut optimum:   ecutwfc   = 40.0,  ecutrho   = 400.0</div></div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px"><img src="cid:ii_14f4bd3a739c2fd6" alt="Inline image 2" width="312" height="187" style="margin-right: 0px;"><br></span></div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px"><br></span></div><div><span style="font-size:12.8000001907349px">I use from this comment Prof. Stefano  in this question from forum.</span></div><div><br></div><div><font face="times new roman, serif">"<span style="color:black;text-align:justify">My recommended strategy to fix ecutwfc and
ecutrho is to perform total </span><span style="color:black;text-align:justify"> </span><span style="color:black;text-align:justify">energy (and
possibly, force and stress) convergence test increasing ecutwfc keeping ecutrho
at its default vaule (=4*ecutwfc)</span><span style="color:black;text-align:justify">  </span><span style="color:black;text-align:justify">until
satisfactory stability is reached (typically ~</span><span style="color:black;text-align:justify;background-image:initial;background-color:yellow;background-repeat:initial">1 mry/atom in the energy, 1.d-4 ry/au in the forces, a
fraction of a KBar in the stress)</span><span style="color:black;text-align:justify"> ...</span><span style="color:black;text-align:justify"> 
</span><span style="color:black;text-align:justify">this fixes the converged value of ecutrho to 4 times the resulting
ecutwfc.</span></font></div>

<p class="MsoNormal" style="margin-bottom:0.0001pt;text-align:justify"><span style="color:black"><font face="times new roman, serif">Now keeping this value for ecutrho one can try
to reduce ecutwfc and see how much this can be done without deteriorating the
convergence."</font><span style="font-family:'Times New Roman',serif;font-size:13pt"></span></span></p><div><span style="font-size:12.8000001907349px">I use PAW pseudo potential . </span></div><div><span style="font-size:12.8000001907349px"><br></span></div><div><span style="font-size:12.8000001907349px">first I use Ecutrho in default (4*Ecut) but don't reach difference of energy from mRy.</span></div><div><br></div><div><span style="font-size:12.8000001907349px"><img src="cid:ii_14f4bd537611abd3" alt="Inline image 3" width="286" height="172" style="margin-right: 0px;"><br></span></div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px"><br></span></div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">second I select Ecutrho=256 and fix Ecut rho and changed Ecut wfc .</span></div><div><span style="font-size:12.8000001907349px">and Ecutwfc=52 and Ecut rho=256 is nearly converged.</span></div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px"><br></span></div><div style="font-size:12.8000001907349px">in higher Ecut wfc oscillation of energy is seen.</div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px"><img src="cid:ii_14f4bdf552fbe0c8" alt="Inline image 4" width="285" height="175" style="margin-right: 0px;"><br></span></div><div style="font-size:12.8000001907349px"><pre style="color:rgb(0,0,0)">Any help would be appreciated.</pre></div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px">Elham</span></div></div>