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>
<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>