Skip to Content.
Sympa Menu

forum - Re: [abinit-forum] How to deal with this warning ?

forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

Re: [abinit-forum] How to deal with this warning ?


Chronological Thread 
  • From: lan haiping <lanhaiping@gmail.com>
  • To: forum@abinit.org
  • Subject: Re: [abinit-forum] How to deal with this warning ?
  • Date: Sat, 25 Apr 2009 03:11:09 +0800
  • Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:date:message-id:subject:from:to :content-type; b=qO12HIaNDYKoMyv00lI1QTAn+uA7HfcmnvRnzlNBVKxN4k1HXZo6SgcNNXneMwjFl6 DbStVYQRlUob8liyZ76f5rmiExTDp3jlB4ljZfrTJqSUvif24dhQmklsyJc34w7rvmSt tYWR21AO98Y3ia8ivVxq+hRHvOSMGIegmStt0=

Dear Matteo,

Thank you !


On Sat, Apr 25, 2009 at 2:41 AM, Matteo Giantomassi <Matteo.Giantomassi@uclouvain.be> wrote:
> Dear All,
> I try to use closure trick for empty states in GW calculations, but in
> several test calculations, i always came to WARNING below, and it seems
> that
> i cannot just tune "gwencomp" to avoid this WARNING.
> Would you please do me a favor for this problem ?

 The warning is not due to the value used for gwencomp.
 It's related to the fact that the closure trick requires
 the Fourier components of the product u^*(r) u(r) at G1-G2,
 where G1 and G2 are two plane waves used to describe the
 irreducible polarizability.

 Now, it might happen that G1-G2 falls outside the FFT box, in this
 case the contribution given by this term is neglected.
 This should not affect the accuracy of the results
 provided that the number of G1-G2 pairs is relatively small (as in your
case).

 To get rid of the warning there are several possibilities.
 Either you increase the value of ecutwfn  or you decrease ecuteps.
 The former is the preferred choice, usually I use all the G-vectors
 stored in the KSS file to avoid possible inaccuracies due to
 the truncation of the plane wave representation (the CPU effort
 scales almost linearly wrt npwwfn)
other persons also suggested me leaving ecutwfn as ecut, but i found  physical memory is exhausted with default value for ecutwfn when performing screening calculations, especially for q(0.00,0.00, 0.00).  Therefore i should do testings on ecutwfn and ecuteps to figure out the proper converged values. Do you have any tricks to tune physical memory load for screening calculation ?


Best wishes,

Hai-Ping


 


 Hope it helps
 Regards
 Matteo Giantomassi



>
>  completechi0_deltapart : WARNING  Number of G1-G2 pairs outside the
> G-sphere for Wfns =    72
>
>  completechi0_deltapart : WARNING  Number of G1-G2 pairs outside the
> G-sphere for Wfns =    72




--
Hai-Ping Lan
Department of Electronics ,
Peking University , Bejing, 100871
lanhaiping@gmail.com, hplan@pku.edu.cn



Archive powered by MHonArc 2.6.15.

Top of Page