<div dir="ltr"><div>Dear Joshua,<br></div>        If you send me the full input file I will be glad to test the PAW datasets constructed with ATOMPAW code.    THe functions on our webpage <a href="http://pwpaw.wfu.edu">http://pwpaw.wfu.edu</a> are mostly focused on accuracy in comparison with the all-electron results and may not rapidly converge with ecut.  However, there are several other efforts to generate PAW functions that may converge more rapidly....   Sincerely, Natalie Holzwarth<br><br></div><div class="gmail_extra"><br clear="all"><div>N. A. W. Holzwarth                                       email: <a href="mailto:natalie@wfu.edu" target="_blank">natalie@wfu.edu</a><div>Department of Physics                                  web: <a href="http://www.wfu.edu/~natalie" target="_blank">http://www.wfu.edu/~natalie</a></div><div>Wake Forest University                                 phone: 1-336-758-5510 </div><div>Winston-Salem, NC 27109 USA                     office: Rm. 300 Olin Physical Lab</div></div>
<br><div class="gmail_quote">On Fri, Sep 19, 2014 at 4:16 PM, Joshua Davis <span dir="ltr"><<a href="mailto:davis101@chemistry.msu.edu" target="_blank">davis101@chemistry.msu.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="ltr" style="font-family:arial,sans-serif;font-size:13px">The plot of the convergence does look like it is oscillating and sort of converging to a point but at fairly large energy cut-off. These potentials were off the "ready to use" potentials on the QE website.  Would you have confidence in the potentials?  PAW is the newest method and there are not many of them.  They all use have scalar relativity and not fully relativity.  I figured that would be good enough.<div><br></div><div>Thanks</div><div><br></div></div><div class="gmail_extra"><br clear="all"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><font face="arial, helvetica, sans-serif">----------------------------------------------------------------------------------------------------------------</font></div><div dir="ltr"><font face="arial, helvetica, sans-serif">Joshua D. Davis</font></div><div dir="ltr"><font face="arial, helvetica, sans-serif"><br>Graduate Assistant</font></div><div dir="ltr"><div><span style="background-color:transparent"><font face="arial, helvetica, sans-serif">Department of Chemistry</font></span></div><div><span style="background-color:transparent"><font face="arial, helvetica, sans-serif">Michigan State University</font></span></div><div><span style="background-color:transparent"><p style="margin:0in 0in 0pt"><font face="arial, helvetica, sans-serif"><span lang="EN-US">578 S. Shaw Lane</span><u></u><u></u><span lang="EN-US">, room 432<u></u></span></font></p></span></div><div><font face="arial, helvetica, sans-serif"><span style="background-color:transparent">East Lansing, MI 48824</span><br></font></div><div><span style="background-color:transparent"><font face="arial, helvetica, sans-serif">-----------------------------------------------------------------------------------------------------------------</font></span></div></div></div></div></div>
<br><div class="gmail_quote">On Fri, Sep 19, 2014 at 3:43 PM, Paolo Giannozzi <span dir="ltr"><<a href="mailto:paolo.giannozzi@uniud.it" target="_blank">paolo.giannozzi@uniud.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Maybe the nonlinear core correction is slowly converging,<br>
especially in Sodium which has s and p semicore electrons<br>
and a very localized (pseudo-)core charge. Just guessing.<br>
<br>
P.<br>
<span><br>
On Fri, 2014-09-19 at 14:20 -0400, Joshua Davis wrote:<br>
> I was performing a convergence test for the energy cut-off for the<br>
> density for a set of PAW potentials just to make sure, but I have some<br>
> results that confuse me a bit.<br>
><br>
><br>
> The compound is a NaB5C (space group: Pm-3m) so the potentials I used<br>
> were:<br>
><br>
> B.pbe-n-kjpaw_psl.0.1.UPF<br>
> C.pbe-n-kjpaw_psl.0.1.UPF<br>
> Na.pbe-spn-kjpaw_psl.0.2.UPF<br>
><br>
><br>
> The cell was also completely relaxed each time with cell and ion<br>
> dynamic being "bfgs" algorithm<br>
><br>
><br>
><br>
> Already did a ecutwfc and kpoint convergence test for each and I found<br>
> and ecutwfc of 35 Ry and a kpoint mesh of 4x8x8 to be sufficient (I am<br>
> calculating a cell doubled in one direction).<br>
><br>
><br>
> When I did the ecutwfc convergence test I went from 10 to 50 Ry in<br>
> increments of 5 Ry, and I kept the default ecutrho=4*ecutwfc for each<br>
> run.<br>
><br>
><br>
> so when I did the ecutrho convergence kept ecutwfc at 35 Ry and the<br>
> kpoints at 4x8x8 and started at 140 and went to 400 Ry.<br>
><br>
><br>
> The graph of my convergence is shown below:<br>
</span><div><div>> Inline image 1<br>
><br>
><br>
> and if the fourm scrubs the image here is the raw data:<br>
> ecutrho(Ry),TotalEnergy(ry)<br>
> 140,-370.0979858<br>
> 150,-370.0980023<br>
> 160,-370.0995505<br>
> 170,-370.1019804<br>
> 180,-370.1047998<br>
> 190,-370.1072663<br>
> 200,-370.1086284<br>
> 250,-370.1048218<br>
> 300,-370.1049835<br>
> 400,-370.1060986<br>
><br>
><br>
><br>
><br>
> The convergence shows a density cut-off of 140 Ry and 150 Ry to be<br>
> pretty similar in calculated energy, but after that the calculated<br>
> energy quickly drops until a cut-off of 200 Ry.  After that the<br>
> calculated energy shoots up at a cut-off of 250 Ry by 0.004 then<br>
> tappers down  by 0.002 Ry as it gets to a cut-off of 400 Ry.<br>
><br>
><br>
> I thought the energy cut-off for density would be the easiest after<br>
> ecutwfc and kpoint, but I want to make sure my results would be fine<br>
> based on what I read for he input documentation for ecutrho.  It says<br>
> ecutrho=4*ecutwfc is usually good (which is pretty standard for the<br>
> PAW method), but I wanted to make sure.  I thought the calculated<br>
> energy would not have changed or the calculated energy would have<br>
> converged at a reasonable value.<br>
><br>
><br>
> Am I doing something wrong? Is the ratio of ecutwfc and ecutrho too<br>
> off at the higher energy ecutrho cut-off?  Other posts suggest<br>
> lowering ecutwfc and keeping ecutrho cut-off constant for proper<br>
> convergence, should I have done this?  I figured sneaking up ecutrho<br>
> would be okay.  Am I dumb and this is just an artifact of the PAW<br>
> method at a large ratio of ecutwfc to ecutrho?  Is this an artifact of<br>
> how the FFT mesh is implemented in QE 5.1?  Should I have not gone<br>
> that high in my density cut-off?<br>
><br>
><br>
> Any help would be much appreciated.<br>
><br>
><br>
> Other considerations:<br>
> I used QE 5.1 compiled with ifort 13 and mkl 11.<br>
><br>
><br>
>  &CONTROL<br>
>    calculation = 'vc-relax',<br>
>    etot_conv_thr = 1.0D-5,<br>
>    forc_conv_thr = 1.0D-4<br>
>  /<br>
><br>
><br>
>  &SYSTEM<br>
>    ibrav = 0,<br>
>    nat = 14,<br>
>    ntyp = 3,<br>
>    nosym = .TRUE.,<br>
>    nspin = 1,<br>
>    occupations = 'fixed',<br>
><br>
> ----------------------------------------------------------------------------------------------------------------<br>
> Joshua D. Davis<br>
><br>
> Graduate Assistant<br>
> Department of Chemistry<br>
> Michigan State University<br>
> 578 S. Shaw Lane, room 432<br>
><br>
><br>
> East Lansing, MI 48824<br>
><br>
> -----------------------------------------------------------------------------------------------------------------<br>
><br>
</div></div>> _______________________________________________<br>
> Pw_forum mailing list<br>
> <a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
> <a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
<br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</blockquote></div><br></div></div>
<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></div>