Skip to Content.
Sympa Menu

forum - Re: [abinit-forum] about nband convergence in GW calculation

forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

Re: [abinit-forum] about nband convergence in GW calculation


Chronological Thread 
  • From: lan haiping <lanhaiping@gmail.com>
  • To: forum@abinit.org
  • Subject: Re: [abinit-forum] about nband convergence in GW calculation
  • Date: Thu, 19 Mar 2009 16:16:42 +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=tLX8pN4Yff67qbLQ80UHvUQfWwW98lyGSUu0XBLAYJGqCsIOUnEg6Z/hc2oPIlPXfe FYlOP86HPnBOEnF+Zo9aF0xEmTnUuEnJFpO6WTlo5hhFdbzcvc9Q7XVMIAeNNy8/qA+E hlVUJrmCHKKmsIiJZwnNRC1sYZPjw7gB31sD8=

Thanks, Thanusit.
It seems I need increase the number bands in KSS calculation.
But the tutorial for GW does not have this problem. In fact , the nband for KSS calculation is set to 9.  Does kssform affect
nband for proceeding calculations ? since i found the only difference between my test and tutorial was that my kssform was set to 3.

Regards,

2009/3/19 ธนูสิทธิ์ บุ <thanusit@kku.ac.th>

Dear lan haiping

I'm not sure I get your question correctly. However, it seems that you wonder why you input the nband parameter from 200, 250, ...400, but in the output the number of bands in the self-energy calculation appears to be 100. If this is the case, please check your KSS file. It should be generated at least with nbandkss = 400 for the whole range of your convergence study to be done. If you use the kssform=3 to generate the kss file, the nband and nbandkss set for the kss file generation should be the same number. That is, you may use

      nband=400

      nbandkss=400

or with the "nbdbuf" imposed your may set

     nband=410

     nbandkss=400

     ndbbuf=10
 

Hope this helps.

Regards,

Thanusit Burinprakhon


On March 19, 2009 2:01:04 PM ICT, "lan haiping" <lanhaiping@gmail.com>, "lan haiping" <lanhaiping@gmail.com> wrote:

Dear  list:
   I did a test calculation on bands convergence of GW calculation with
abinit 5.6.3. But it seems that  the "nband" tag  didnot control over the number of bands in self-energy calculation, and the value for "mband" was ok which agreed with my setting for nband.
Below is my setting for GW parts and related output , any comment or advice  is appreciated.

GW  settings:
 ndtset 5
 # GW corrections
  optdriver  4  # driver options for self energy
  getkss 2   # KSS file used
  getscr 1   # SCR file used
  nband: 200  # bands for self energy
  nband+ 50
  ecutwfn  24 Ha #
  ecutsigx  24.0 Ha
  nkptgw    1    #
  kptgw
     0.000 0.500  0.500
  bdgw  54  55

Output
SIGMA fundamental parameters:
 PLASMON POLE MODEL  1
 number of plane-waves for SigmaX                 5961
 number of plane-waves for SigmaC and W           1139
 number of plane-waves for wavefunctions          5961
 number of bands                                   100
 number of independent spin polarizations            1
 number of k-points in IBZ                          13
 number of q-points in IBZ                          13
 number of symmetry operations                       6
 number of k-points in BZ                           64
 number of q-points in BZ                           64
 number of omega for Sigma on real axis              9
 deltae [eV]                                      0.25
 number of omega for Sigma on real axis              0
 max omega for Sigma on real axis  [eV]           0.00


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



Dr Thanusit Burinprakhon
Physics Department, Faculty of Science
Khon Kaen university
Khon Kaen, Thailand
40002

tel 00 66 4320-2222 ext. 2248
mobile  00 66 86-8620510



--
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