<div dir="ltr">Dear Davide,<div><font color="#333333" face="normal arial, sans-serif"><span style="line-height:20px"><font size="3">Last 2-3 days, I also trying to run some Meta GGA calculations with TPSS functional (Very small molecule, like NH3, CH4).</font></span></font></div>
<div><font color="#333333" face="normal arial, sans-serif"><font size="3"><span style="line-height:20px">When the box length was less than 10 angstrom, its works, at least I dont face the convergence failure problem.</span></font></font></div>
<div><font color="#333333" face="normal arial, sans-serif"><font size="3"><span style="line-height:20px">When I increase the box length(15-20 angstrom) i got the similar kind of error with QE-5.0.3 version.</span></font></font></div>
<div><font color="#333333" face="normal arial, sans-serif"><font size="3"><span style="line-height:20px">But its work for QE-4.2.1 version.  May be you can try it with QE-4.2.1 .</span></font></font></div><div><font color="#333333" face="normal arial, sans-serif"><font size="3"><span style="line-height:20px">I dont understand, why the QE-5.0.0 to QE-5.0.3 or trunk version is not working.</span></font></font></div>
<div><span style="line-height:20px;font-size:medium;color:rgb(51,51,51);font-family:'normal arial',sans-serif">Probable thay have changes something in the new Meta GGA code.</span></div><div><span style="line-height:20px;font-size:medium;color:rgb(51,51,51);font-family:'normal arial',sans-serif">Thanks</span></div>
<div><span style="line-height:20px;font-size:medium;color:rgb(51,51,51);font-family:'normal arial',sans-serif">Ananya  </span></div><div><font color="#333333" face="normal arial, sans-serif"><font size="3"><span style="line-height:20px"><br>
</span></font></font></div><div><font color="#333333" face="normal arial, sans-serif"><span style="line-height:20px"><font size="3"> </font><br></span></font><div><span style="color:rgb(51,51,51);font-family:'normal arial',sans-serif;font-size:16px;line-height:20px"><br>
</span></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Apr 25, 2014 at 9:10 AM, Mike Marchywka <span dir="ltr"><<a href="mailto:marchywka@hotmail.com" target="_blank">marchywka@hotmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
<br>
<br>
<br>
<br>
----------------------------------------<br>
> From: <a href="mailto:paolo.giannozzi@uniud.it">paolo.giannozzi@uniud.it</a><br>
> To: <a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a><br>
> Date: Fri, 25 Apr 2014 03:42:06 +0200<br>
> Subject: Re: [Pw_forum] tpss: too many bands are not converged<br>
<div class="">><br>
> On Thu, 2014-04-24 at 09:03 -0400, Mike Marchywka wrote:<br>
><br>
>> By free volume, you mean places where n and likely grad everything are numerical noise?<br>
><br>
> exactly: the regions of space where the charge has decayed to noise.<br>
> Apparently (the current implementation of) meta-GGA doesn't like the<br>
> presence of such large regions (and it doesn't seem to like H atoms<br>
> as well). Even plain GGA diverges in those regions: it is less of a<br>
> problem, but still, the exact energy values depend upon how the<br>
<br>
</div>Does it actually diverge or just reach some kind of limit cycle with<br>
amplitude in energy greater than convergence criteria?  <br>
<div class=""><br>
> potential is cut off. The divergence comes from gradient corrections<br>
> to exchange: there is no divergence in atoms with the correct large-r<br>
<br>
</div>so this is related to PP and lack of exponential decay giving more weight<br>
to "free volume" than it should physically have? I had been doing some<br>
tests with analytical hydrogen since it was easy to check :) That is not noise<br>
alone then but a modelling issue. <br>
<div class=""><br>
> limit of the charge density, but this is not of much help in plane-wave<br>
> calculations of condensed-matter systems.<br>
<br>
</div>Molecular crystals better or worse than say metals? <br>
<br>
>From what I can tell of the TPSS "z" anyway, trying to calculate it<br>
as a ratio of intermediate real quantities seems to allow for many errors<br>
such as aliasing that may occur by taking magnitudes.  I did some<br>
algebra it looks like the magnitude of the grad of the wf angle is more direct but I have<br>
not checked the math or tested it much although IIRC some related results<br>
vs analytical hydrogen did not appear obviously wrong but this was hardly a complete test.<br>
<div class=""><br>
<br>
<br>
><br>
> (disclaimer: this is what I remember. I looked at this problem several<br>
> years ago for meta-GGA, many years ago for GGA)<br>
<br>
<br>
</div>Thanks, I'm just trying to understand a few things and see if I can play with<br>
the problematic terms and evaluate some alternatives.<br>
<div class="HOEnZb"><div class="h5"><br>
><br>
> Paolo<br>
>> AFAICT, the converge of SCF is largely empirical<br>
>> as there does not seem to be much theory here. In any case, you would like to think that the areas lacking "physics" should not be that<br>
>> big a factor in the overall calculation and if that is not the case there may be something interesting here to explore :)<br>
>><br>
>>> You can try to decrease the size of the unit-cell in order to minimize the free space it contains. Of course, this will be at the expense of the inter-polyemer interaction you probably want to avoid.<br>
>>><br>
>>> Nevertheless, this will probably not cure your problem, which to my knowledge has no other solution than changing the meta-GGA functional to a GGA one.<br>
>><br>
>> I took a look at the "z" parameter which being a ratio could have some issue with small denominators. With some initial work,<br>
>> it looks like there are better ways to calculate it than directly from real quantities such as n and grad mags,<br>
>> but I am still trying to test what may be trivially wrong math :)<br>
>> I think my latest result was something to the effect that you could reduce z to A/(A+B)( both non-negative)<br>
>> if you calculate everything from psi and the grad of psi expressed in polar form, not much different from some formula<br>
>> for current densities, and maybe find a limit when A and B are both zero.<br>
>><br>
>> I had originally hoped to try to test it for trivial issues in JDFTX and c++ or c with libxc using an expanded interface scheme<br>
>> to let me pass z or grad psi in polar form but it will require a bit of effort. Curious if anyone else here has examined<br>
>> this or has literature references they care to share. Seems disconcerting that something with physical meaning would have these<br>
>> kinds of implementation issues - perhaps something else can be learned from investigating.<br>
>><br>
>> Thanks.<br>
>><br>
>>> Best regards,<br>
>>><br>
>>> Yves<br>
>>><br>
>>><br>
>>><br>
>>> Le 24 avr. 2014 à 13:36, Davide Tiana a écrit :<br>
>>><br>
>>>><br>
>>>> Dear all,<br>
>>>> I've been trying to calculate a polymer using tpss. I already fully<br>
>>>> optimised it with different functional (pbesol, pw91, hse) and<br>
>>>> pseudopotential (nc, paw). For some reason (apologise my ignorance if<br>
>>>> it should be well known) when I try to calculate using tpss nc the<br>
>>>> system simply does not converge crashing with this error:<br>
>>>><br>
>>>> Error in routine c_bands (1):<br>
>>>> too many bands are not converged<br>
>>>><br>
>>>> After some reaserch in old pw-forum mails, I tried decreasing<br>
>>>> mixing_beta=0.2, increasing mixing_ndim=12, as well as mixing mode TF<br>
>>>> and local-TF.<br>
>>>> everythings failed.<br>
>>>><br>
>>>> Does anyone have any tips or idea of why the system can't converge?<br>
>>>><br>
>>>> Thanks a lot,<br>
>>>> Davide<br>
>>>><br>
>>>><br>
>>>><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>
>>><br>
>>><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>
>><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>
><br>
> --<br>
> Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
> Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
> Phone <a href="tel:%2B39-0432-558216" value="+390432558216">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222">+39-0432-558222</a><br>
><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>
<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>
</div></div></blockquote></div><br></div>