<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hello,<div><br><div><div>Il giorno 27/mag/2011, alle ore 20.05, Tram Bui ha scritto:</div><br class="Apple-interchange-newline"><blockquote type="cite">Dear Gabriele,<br>    Thank you for your respond. Just to make sure that I'm understanding your answer correctly, below is a shortcut on the output file of my PP testing: <br><br>---------------------- Testing the pseudopotential ----------------------<br>
<br><br>     scalar relativistic calculation<br><br>     atomic number is 55.00   valence charge is  1.00<br>     dft =PBE   lsd =0 sic =0 latt =0  beta=0.20 tr2=1.0E-14<br>     mesh =1251 r(mesh) = 101.25750 xmin = -7.00 dx = 0.01250<br>
<br>     n l     nl             e AE (Ry)        e PS (Ry)   <span style="color: rgb(0, 0, 0); background-color: rgb(51, 51, 51);"> De AE-PS (Ry)</span><span style=""><span style="background-color: rgb(51, 51, 51);"></span></span><br>
     1 0     6S   1( 0.50)       -0.17728       -0.17238       -0.00490<br>     2 1     6P   1( 0.50)       -0.07657       -0.07687        0.00030<br>     3 2     5D   1( 0.00)       -0.05580       -0.05081       -0.00499<br>
<br>     eps = 8.5E-15  iter = 22<br><br>     Etot =  -15573.817276 Ry,   -7786.908638 Ha, -211892.560039 eV<br>     Etotps =    -0.217507 Ry,      -0.108753 Ha,      -2.959329 eV<br>     dEtot_ae =       0.049938 Ry<br>     dEtot_ps =       0.048795 Ry,  <span style="background-color: rgb(102, 102, 102);"> </span><span style="color: rgb(255, 153, 0); background-color: rgb(102, 102, 102);">Delta E=       0.001144 Ry</span><br>
<br>     Ekin =       0.080184 Ry,       0.040092 Ha,       1.090963 eV<br>So I should look at the "De AE-PS (Ry)", which is the difference in the eigenvalues? or the "Delta E"? I assume it is the Delta E based on your respond. But am I correct?<br></blockquote><div><br></div><div>I would say that you should look at both, but what Paolo was referring to in the guide is the second thing.</div><div>The Delta E looks fair to me, that's what you can expect from standard USPP. The difference in the eigenvalues seems instead a bit too large, maybe you can still improve a bit the quality of the PP.</div><div><br></div><div><br></div><div>HTH</div><div><br></div><div>GS</div><br><blockquote type="cite">
<br>Regards,<br>Tram<br><br><div class="gmail_quote">On Thu, May 26, 2011 at 2:16 AM, Gabriele Sclauzero <span dir="ltr"><<a href="mailto:sclauzer@sissa.it">sclauzer@sissa.it</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;">
<div style="word-wrap: break-word;"><br><div><div>Il giorno 26/mag/2011, alle ore 01.00, Tram Bui ha scritto:</div><div class="im"><br><blockquote type="cite">Dear Paul,<br>       In the file that I have attached here (Paul's note on PP generation), under section 3.10: checking for transferability, you have mentioned about the:"The error on total energy differences<br>
between PP and AE results gives a feeling on how good the PP is." Do you mean the difference in AE and PP eigenvalues? </blockquote><div><br></div></div><div>No, I think here Paolo "Paul" Giannozzi really meant what you understood first, but you didn't get it in the correct way.</div>
<div>You should really compare total energy differences, which means:</div><div>E_2(AE) - E_1(AE)    versus   E_2(PS) - E_1(PS)</div><div>where 1 and 2 label two different electronic configurations, and NOT</div><div>E_1(PS)  versus  E_1(AE)    </div>
<div>as you are probably trying to do.</div><div class="im"><br><blockquote type="cite">because I saw the Etot and Etotps are huge difference </blockquote><div><br></div></div><div>This is often the case, given the high degree of freedom in the choice of PP parameters (especially for the local potential, I suppose).</div>
<div class="im"><div><br></div><blockquote type="cite">but their eigenvalue (e, Ry) are closed.<br></blockquote><div><br></div></div><div>Indeed this is how it should be, since the modern PP are designed to reproduce the scattering properties of the AE atomic potential as best as possible.</div>
<div class="im"><br><blockquote type="cite">
<br>Dear QE users,<br>     I'm trying to test the Cs' US PP with reference configuration of :[Xe] 6s1 5d0 6p0. To test the PP, I'm using these five configuration (just by guessing): <br>     configts(1)= '6s1 6p0 5d0'<br>

     configts(2)= '6s0.5 6p0.5 5d0'<br>     configts(3)= '6s0.5 6p0 5d0.5'<br>     configts(4)= '6s0 6p1 5d0'<br>     configts(5)= '6s0 6p0 5d1'<br>My questions are: <br> - how to choose the configurations for PP testing based on the reference one? is there a strategy or some sort?<br>
</blockquote><div><br></div></div><div>Use some reasonable configurations (remove or add one electron, or a fraction of it, move one or more electrons between different orbitals, ...).</div><div>The choice may also be driven by the class of materials that you are going to study (e.g. ionic crystals).</div>
<div class="im"><br><blockquote type="cite">
 - how many configurations are enough for the test?<br></blockquote><div><br></div></div><div>That's up to you... 5 look enough for me. Don't forget to do some testing with the solid state code (e.g. bond distance of dimers, lattice constant of simple solids...).</div>
<div><br></div><div><br></div><div>HTH</div><div><br></div><div><br></div><div>GS</div><br><blockquote type="cite"><div class="im"><br>Thank you very much,<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>
</div><span><pseudo-gen.pdf></span>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org" target="_blank">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>
</blockquote></div><font color="#888888"><br><div>
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div>
<span style="color: rgb(126, 126, 126); font-size: 16px; font-style: italic;"><br>§ Gabriele Sclauzero, EPFL SB ITP CSEA</span></div><div><font color="#7e7e7e"><i>   PH H2 462, Station 3, CH-1015 Lausanne</i></font></div>
</span>
</div>
<br></font></div><br>_______________________________________________<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>
<br></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>
_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>http://www.democritos.it/mailman/listinfo/pw_forum<br></blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div><span class="Apple-style-span" style="color: rgb(126, 126, 126); font-size: 16px; font-style: italic; "><br class="Apple-interchange-newline">§ Gabriele Sclauzero, EPFL SB ITP CSEA</span></div><div><font class="Apple-style-span" color="#7E7E7E"><i>   PH H2 462, Station 3, CH-1015 Lausanne</i></font></div></span>
</div>
<br></div></body></html>