<div dir="auto">I've recently been using the spdf pseudopotentials for Ce and Gd from the pslibrary. I can say that:<div dir="auto">1. They require a monster cutoff (200 Ry for Ce, and 100 for Gd) </div><div dir="auto">2. They require +U, especially Ce</div><div dir="auto">3. Results were quite good</div><div dir="auto"><br><div data-smartmail="gmail_signature" dir="auto">-- <br>Lorenzo Paulatto<br>Written on a virtual keyboard with real fingers</div></div></div><br><div class="gmail_quote"><div dir="ltr">On Mon, 1 Oct 2018, 18:28 Nicola Marzari, <<a href="mailto:nicola.marzari@epfl.ch">nicola.marzari@epfl.ch</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
<br>
Dear Martin,<br>
<br>
we tested a couple of pseudo libraries with the f frozen, but apart from <br>
La they do not work that well in general. As mentioned by Christopher, <br>
you want the f.<br>
<br>
Here are a few tests:<br>
<a href="https://www.materialscloud.org/discover/sssp/plot/efficiency/Ho" rel="noreferrer noreferrer" target="_blank">https://www.materialscloud.org/discover/sssp/plot/efficiency/Ho</a><br>
<br>
and if you eg brose the provennce tree for the Ho US pslibrary 1.0<br>
until you get to a PWscf calculation:<br>
<br>
<a href="https://www.materialscloud.org/explore/sssp/details/f2383960-43b6-4bf6-b080-b49e43c6ce12?nodeType=CALCULATION" rel="noreferrer noreferrer" target="_blank">https://www.materialscloud.org/explore/sssp/details/f2383960-43b6-4bf6-b080-b49e43c6ce12?nodeType=CALCULATION</a><br>
<br>
you can see that the pseudo tested is<br>
Ho.pbe-spdfn-rrkjus_psl.1.0.0.UPF<br>
<br>
I would tend to recommend the Wentzocvitch series, but it's true that <br>
all rare earths are rarely tested. In addition, they have strong <br>
multiconfigurational character, so standard LDA/GGA/GGA+U/HSE<br>
DFT doesn't do that well.<br>
<br>
                                nicola<br>
<br>
<br>
<br>
<br>
On 01/10/2018 10:30, Ing. Martin Matas wrote:<br>
> Dear experts,<br>
> <br>
> I appreciate that many new ultrasoft pseudopotentials for rare earths by Andrea Dal Corso have recently appeared on QE homepage. I performed calculations with a 2-atom primitive cell of metal holmium using Ho.pbe-spdn-rrkjus_psl.1.0.0.UPF, with both ferromagnetic and antiferromagnetic starting configurations.<br>
> <br>
> The resulting magnetization is always zero, not depending on calculation parameters (like starting_magnetization, occupation type and degauss, or energy cutoffs). However, exactly the same calculations using a norm-conserving holmium pseudopotential give results with a non-zero magnetization what I consider correct since holmium at zero temperature is ferromagnetic.<br>
> <br>
> Did anyone test these new ultrasoft pseudopotentials in this way and is there a way how to get magnetically correct results using them?<br>
> Thanks a lot.<br>
> <br>
> Martin Matas<br>
> PhD student<br>
> University of West Bohemia<br>
> Czech Republic<br>
> <br>
> _______________________________________________<br>
> users mailing list<br>
> <a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a><br>
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
> <br>
<br>
<br>
-- <br>
----------------------------------------------------------------------<br>
Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL<br>
Director, National Centre for Competence in Research NCCR MARVEL, EPFL<br>
<a href="http://theossrv1.epfl.ch/Main/Contact" rel="noreferrer noreferrer" target="_blank">http://theossrv1.epfl.ch/Main/Contact</a> <a href="http://nccr-marvel.ch/en/project" rel="noreferrer noreferrer" target="_blank">http://nccr-marvel.ch/en/project</a><br>
_______________________________________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org" target="_blank" rel="noreferrer">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
</blockquote></div>