[Wien] Charge Convergence is not achieved

shamik chakrabarti shamikphy at gmail.com
Thu Jul 7 16:18:13 CEST 2011


Dear Dr. Peter Blaha Sir,

                        Thank you very much for all your responses. I have
to think over all the discussions calmly and it should definitely help us to
understand the basics of this convergence problem.

Thank you very much Sir.

with best regards,

On Thu, Jul 7, 2011 at 6:58 PM, Peter Blaha <pblaha at theochem.tuwien.ac.at>wrote:

> The problem of "functional" could for instance happen for 4f-compounds.
>
> In LDA/GGA all 4f bands will be around EF (which is physically wrong) and
> convergence is naturally very different (14 extremely FLAT bands/atom !!).
>
> Using "open core" or LDA+U or Hybrid-DFT you remove the failure of GGA and
> immediately also convergence should improve (since you have only a FEW
> WIDE s,d-bands at EF).
>
> Am 07.07.2011 15:05, schrieb shamik chakrabarti:
>
>> Dear Laurence Marks Sir,
>>
>>               Thank you very much for your reply.....yes the question may
>> not have a very good answer....if in any calculation we are not getting the
>> desired convergence, it may
>> happen that:
>> (1) Our chosen functional is not appropriate for the system
>>
>> (2) The system (say while trying to predict a new material!) may not be
>> feasible at all..
>>
>> Probably all we can say that if we are able to achieve desired convergence
>> (say 0.0001) we can say that we have used the appropriate functional for the
>> system and the system is
>> (may be) feasible (at least theoretically!! ).
>>
>> Sir please correct me if I am wrong in my concept.
>>
>> with best regards,
>>
>> On Thu, Jul 7, 2011 at 6:12 PM, Laurence Marks <L-marks at northwestern.edu<mailto:
>> L-marks at northwestern.**edu <L-marks at northwestern.edu>>> wrote:
>>
>>    2011/7/7 Shamik Chakrabarti <shamikiitkgp at gmail.com <mailto:
>> shamikiitkgp at gmail.com**>>:
>>
>>     > Dear Peter Blaha Sir,
>>     >                        Indeed by increasing number of K points we
>> got the
>>     > convergence. Sir I have now some basic queries on this topic. You
>> have said
>>     > that
>>     > "sometimes you cannot reach (easily) arbitrary
>>     > convergence"
>>     > why in some cases we can not reach convergence up to our desired
>> limit?...is
>>     > it the limitation of DFT?....or it means that the feasibility of the
>>     > solution is only up to the achieved convergence?
>>
>>    This is in fact a deep, and very good question, at least in my opinion.
>>
>>    Unfortunately that does not mean that there is a good answer to it!
>>
>>    With the perfect functional convergence should (I believe, others may
>>    disagree) always be good. With a very imperfect functional it is quite
>>    possible that a DFT calculation will not converge, i.e. it is
>>    unfeasible. Empirically many (but not all) metals do not converge well
>>    with small numbers of k-points, but some others do. Why....I do not
>>    understand as I cannot write down a mathematical analysis to explain
>>    this and do not believe that there is a formal analysis in the
>>    literature, it is just empirical knowledge (folklore).
>>
>>
>>    --
>>    Laurence Marks
>>    Department of Materials Science and Engineering
>>    MSE Rm 2036 Cook Hall
>>    2220 N Campus Drive
>>    Northwestern University
>>    Evanston, IL 60208, USA
>>    Tel: (847) 491-3996 Fax: (847) 491-7820
>>    email: L-marks at northwestern dot edu
>>    Web: www.numis.northwestern.edu <http://www.numis.**northwestern.edu<http://www.numis.northwestern.edu>
>> >
>>
>>    Chair, Commission on Electron Crystallography of IUCR
>>    www.numis.northwestern.edu/ <http://www.numis.**northwestern.edu/<http://www.numis.northwestern.edu/>
>> >
>>
>>    Research is to see what everybody else has seen, and to think what
>>    nobody else has thought
>>    Albert Szent-Gyorgi
>>    ______________________________**_________________
>>    Wien mailing list
>>    Wien at zeus.theochem.tuwien.ac.**at <Wien at zeus.theochem.tuwien.ac.at><mailto:
>> Wien at zeus.theochem.**tuwien.ac.at <Wien at zeus.theochem.tuwien.ac.at>>
>>
>>    http://zeus.theochem.tuwien.**ac.at/mailman/listinfo/wien<http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien>
>>
>>
>>
>>
>> --
>> Shamik Chakrabarti
>> Research Scholar
>> Dept. of Physics & Meteorology
>> Material Processing & Solid State Ionics Lab
>> IIT Kharagpur
>> Kharagpur 721302
>> INDIA
>>
>>
>>
>> ______________________________**_________________
>> Wien mailing list
>> Wien at zeus.theochem.tuwien.ac.**at <Wien at zeus.theochem.tuwien.ac.at>
>> http://zeus.theochem.tuwien.**ac.at/mailman/listinfo/wien<http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien>
>>
>
> --
>
>                                      P.Blaha
> ------------------------------**------------------------------**
> --------------
> Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
> Phone: +43-1-58801-15671             FAX: +43-1-58801-15698
> Email: blaha at theochem.tuwien.ac.at    WWW: http://info.tuwien.ac.at/**
> theochem/ <http://info.tuwien.ac.at/theochem/>
> ------------------------------**------------------------------**
> --------------
> ______________________________**_________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.**at <Wien at zeus.theochem.tuwien.ac.at>
> http://zeus.theochem.tuwien.**ac.at/mailman/listinfo/wien<http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien>
>



-- 
Shamik Chakrabarti
Research Scholar
Dept. of Physics & Meteorology
Material Processing & Solid State Ionics Lab
IIT Kharagpur
Kharagpur 721302
INDIA
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20110707/3d89e0a2/attachment.htm>


More information about the Wien mailing list