Dear Nicola, Hanghui and Trinh,<br>      <br>         I'm sorry for interupting but I've been trying to find someone who have had experiences with generating an ultrasoft pseudopotential for QE, and I saw your email back and forth so I just want to give it a shot to ask you couple of questions. Please accept my apology for any inconvenience. <br>
         By the way, I've tried to generate an US PP for Cs and have not get much luck with it. To be honest I really don't have any background on PP, so it was really tough for me to be able to understand what each variable means and how to determine them (even though I have read the PP generation notes for QE users, and several papers or information online about PP). So would you spare me some time and give me some guidance on how to build a meaningful input file for PP generation? <br>
        For Cs (Z=55), I want to include the 5s, 5p and 6s state for my valence and below is the input file I have build so far (it is not running correctly though). However, there are several things that I'm still unclear about.<br>
&input<br>     title='Cs',<br>     zed=55.0,<br>     prefix='Cs',<br>     iswitch=3,<br>     rel=1,<br>     dft='PBE',<br>     config='[Kr] 4d10 5s2 5p6 6s1 ',<br><font color="#666600">till this part, I understand that is to run the all-electron (AE) calculation </font><br>
<br>     rlderiv=4.20, <font color="#ff6600">(here I'm not quite sure how one can come up with the value for rlderiv, in this case I used the radius at most outpeak, 6s, from the AE wfc plot), so is it a right approach to obtain rlderiv? and what effect does it have on pseudizing the wfc process? is this the same as the cut of radius?</font> <br>
     eminld=-8.0,<br>     emaxld=3.0,<br>     deld=0.01d0,<br>     nld=3, <font color="#ff6600">For this value, I assume that I want to include 3 energy state, 5s 5p 6s, so I need to have 3 number of logarithmic derivative. is that correct?</font><br>
/<br> &inputp<br>     lloc=0, <font color="#ff6600">here, from the ld1-input.txt, i know the lloc has to be equal to the last lls # in the card list below (in this case =0), but what does it mean by the "local chanel"? (what is a local chanel?) if I have lloc=0 here, does it mean that my local chanel is the 5s ?</font><br>
     pseudotype=3,<br>     rho0=0.01,<br>     file_pseudopw='Cs.uspp-tb.UPF'<br>/<br>4<br>6S 1 0 1.00 0.00 <span style="color:rgb(0,102,0)">4.20 4.50</span><br>6S 1 0 0.00<span style="color:rgb(0,0,153)"> 1.00</span> <span style="color:rgb(0,102,0)">4.00 4.50</span><br>
5P 2 1 6.00 0.00 <span style="color:rgb(0,102,0)">1.60 1.60</span><br>5S 1 0 2.00 0.00 <span style="color:rgb(0,102,0)">1.40 1.40</span><br><font color="#cc0000"> For the wfc configuration above,<br>- what does it indicate when one have difference (or the same) rcuts for norm-conversing and US PP? </font><br>
<font color="#cc0000">- when rcuts of US and NC are difference, why do one need two energies/channel?</font><br><span style><span style="background-color:rgb(204,0,0)"></span></span><font color="#cc0000">- In the case of when the second row is needed when rcuts are different, based on what one can determine the value for the "ener" (energy used to pseudize correspoding state)? I only know what it ener is zero for bound state, but most of the examples i looked at, everyone was using a nonzero value when they have another row for the same wfc( e.g. 6s in this case), does it mean the second row is the unbound state?</font><br>
<br>I'm really sorry for the long email and the color code i used if that is bothering you. I'm doing this all alone without any help of any expert around here. So things are very confusing to me especially when it comes to quantum physics in which I have very limited understanding. I really appreciate any input you have for me.<br>
<br>Best Regards,<br>Tram Bui<br><br><div class="gmail_quote">On Wed, Apr 27, 2011 at 12:48 PM, Vo, Trinh (388C) <span dir="ltr"><<a href="mailto:Trinh.Vo@jpl.nasa.gov">Trinh.Vo@jpl.nasa.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Dear Nicola,<br>
<br>
I am working on generating pseudo potentials for a few elements of<br>
Lanthanides too, but so far I am not successful yet (it is still bad).  I<br>
appreciate if you could share with me your parameter input.<br>
<br>
Thank you very much in advance,<br>
<br>
Trinh<br>
<br>
<br>
On 4/27/11 10:30 AM, "Nicola Marzari" <<a href="mailto:nicola.marzari@materials.ox.ac.uk">nicola.marzari@materials.ox.ac.uk</a>><br>
wrote:<br>
<div class="im HOEnZb"><br>
> On 4/27/11 5:45 PM, hanghui chen wrote:<br>
>> Dear PWscf community,<br>
>>       I am wondering if any one has succeeded in generating the<br>
>> ultrasoft pseudopotential of Eu, using uspp-736. Now in the literature,<br>
>> for most of the work that involves Eu, the calculation is done with VASP<br>
>> where the pseudopotential is provided in its own library.<br>
>>       Thank you very much indeed.<br>
>><br>
>> Hanghui Chen<br>
>> Department of Physics,<br>
>> Yale University<br>
><br>
><br>
> Dear Hanghui,<br>
><br>
> we have a test one normconserving, with f electrons in the valence, and<br>
> an ultrasoft one, with f electrons in the core. This is part of a<br>
> project developing and testing lanthanides pseudos, to be released later<br>
> on, but if you want to give them a try, I can send them to you (would<br>
> be great to have some test feedback).<br>
><br>
> nicola<br>
><br>
<br>
</div><div class="HOEnZb"><div class="h5">_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Tram Bui<br><br>M.S. Materials Science & Engineering<br><a href="mailto:trambui@u.boisestate.edu" target="_blank">trambui@u.boisestate.edu</a><br><br>