Skip to Content.
Sympa Menu

forum - Re: [abinit-forum] Why my SCF calculation converges so slow?

forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

Re: [abinit-forum] Why my SCF calculation converges so slow?


Chronological Thread 
  • From: "Anglade Pierre-Matthieu" <anglade@gmail.com>
  • To: forum@abinit.org
  • Subject: Re: [abinit-forum] Why my SCF calculation converges so slow?
  • Date: Tue, 6 Jun 2006 11:26:53 +0200
  • Domainkey-signature: a=rsa-sha1; q=dns; c=nofws; s=beta; d=gmail.com; h=received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:references; b=YQwqzpCmLL2jkrm7EMAB4aAaFcBLVBrnJU0x15rY4PAVh3GGcOfHCgyO85krZA4oW6T9uH9/fTyAHAVLVSWDkgP9vfhgsMtanG0t9+mlVOaFKU0FIcBSn6POzr2ozq7AE3cuXPIh/0++THToyoprnbXhxl3Qo6oCaTxfSxpKpls=

Hi,

It might have been useful to join you input file... Our comment would have been more accurate then.
Nevertheless here are a few basic recommandation:
- try changing the precondtioning shcheme your using. For very difficult system
"iprcel 49" may work best.
-try changing the mixing scheme: even simple mixing can sometimes out perform conjugate gradient. In abinit version 5.x.x you have "iscf 7" which is usualy very efficient
-try using "nline " and "nnsclo" they'll help improving the true convergence of your caclulation
-and yes tolvrs 1e-18 is almost mandatory if you wan't an nice phonon calculation.

regards

PMA

On 6/6/06, zhangtingPKU <zhangting@pku.edu.cn> wrote:
Dear ABINIT users:

    I'm calculating the (5,5) CNT's phonon spectrum with ABINIT code, with the input file like in the lesson_rf2. But I found that the first database, which is the initio SCF calculation, converges very slow. With the accuracy in the lesson, tolvrs=1E-18, the calculation has already run for about 600 SCF-steps, but has not reached the convergence condition yet. In fact it oscillates at about deltavrs = 1E-17 ~ 1E-16. My question is that in the first step in a phonon calculation, does it really necessary to set a converge condition as telvrs=1E-18 ? Or am I got something wrong within the input file, such that it can't reach the converge condition?(My input files is attached) Thanks!





礼!


zhangtingPKU
zhangting@pku.edu.cn
2006-06-06






--
Pierre-Matthieu Anglade



Archive powered by MHonArc 2.6.16.

Top of Page