
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN">

<html><body style='font-family: Verdana,Geneva,sans-serif'>

<p>Dear Yogesh, <br /><br />I assume you are trying to use PSLibrary ultrasoft PPs. <br />The one you've listed is spdn, meaning no 4f in valence. <br /><br />PSLibrary's lanthanides with 4f (spdfn) require crazy high plane wave cutoffs. <br />As long as you are using +U anyway, I recommend Topsakal's PAWs. <br />http://www.mineralscloud.com/resources/repaw/index.shtml<br />Please see the respective references to get the values of +U.<br /> <br />However, note that they've changed the QE code to use +U on both d and f shells. <br />I'd use these PPs with +U at 4f. <br /><br />You can use USPPs for the other atoms, or find other PAWs for them.<br /><br />Best regards. <br />Andrii Shyichuk, University of Wrocław.<br /><br />W dniu 2020-10-28 05:59, Yogesh Kumar napisał:</p>

<blockquote type="cite" style="padding-left:5px; border-left:#1010ff 2px solid; margin-left:5px"><!-- html ignored --><!-- head ignored --><!-- meta ignored -->

<div dir="ltr">Dear QE experts,

<div> </div>

<div>I was trying to calculate PDOS of CePO4 using qe-6.2.1. After completing all the steps, I found that No file having Ce(4f) contribution was there. So I looked at Ce PP and found that Ce 4f is not considered in valence electrons (I am not sure if I am right!). Then I also checked the PP of few more Lanthanides and found that no 4f electrons are considered in valence electrons. How can I then calculate Ce 4f contribution in CePO4 or other compounds having Ce or other lanthanides? Pseudopotentials used by me are as following:</div>

<div> </div>

<div>Ce.pbe-spdn-kjpaw_psl.1.0.0.UPF</div>

<div>O.pbe-n-kjpaw_psl.0.1.UPF</div>

<div>P.pbe-n-kjpaw_psl.0.1.UPF</div>

<div> </div>

<div>Thanks and Regards,</div>

<div> </div>

<div>Yogesh Kumar</div>

<div>Scientific Officer</div>

<div>BARC, Mumbai</div>

</div>

<br />

<pre>_______________________________________________

Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu">www.max-centre.eu</a>)

users mailing list <a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>

<a href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a>

</pre>

</blockquote>

<p> </p>

<div> </div>

</body></html>