About convergence

Concerns issues with computing quasiparticle corrections to the DFT eigenvalues - i.e., the self-energy within the GW approximation (-g n), or considering the Hartree-Fock exchange only (-x)

Moderators: Davide Sangalli, andrea.ferretti, myrta gruning, andrea marini, Daniele Varsano

Post Reply
lesheng
Posts: 26
Joined: Thu Feb 07, 2013 11:39 pm

About convergence

Post by lesheng » Wed Mar 13, 2013 4:07 pm

Dear all,

I am a new user of Yambo, and after I followed the tutorial on Si to perform the GW calculation, I have a general question about the convergence test. For each system I am dealing with, I have to carry out a set of calculations with different variables (such as EXXRLvcs, BndsRnX, GbndRnge etc) in order to determine the "best" value for these variables, so I can obtain an accurate result (with convergence outcomes). And such a strategy will be absolutely "time expensive". Is that correct?

Thanks,
PhD candidate
209 South Rd, Chapel Hill, NC 27599
Department of Chemistry
The University of North Carolina at Chapel Hill
Email: leshengl@live.unc.edu

User avatar
Daniele Varsano
Posts: 3816
Joined: Tue Mar 17, 2009 2:23 pm
Contact:

Re: About convergence

Post by Daniele Varsano » Wed Mar 13, 2013 5:24 pm

Dear Lesheng,
the heaviness of the convergence studies depends totally on your system size (number of electrons, k-points etc..).
Most of them are independent so you can push a variable keeping the others low (but reasonable) in order to assess the each converged
value and finally run the production calculations.

Cheers,

Daniele
Dr. Daniele Varsano
S3-CNR Institute of Nanoscience and MaX Center, Italy
MaX - Materials design at the Exascale
http://www.nano.cnr.it
http://www.max-centre.eu/

lesheng
Posts: 26
Joined: Thu Feb 07, 2013 11:39 pm

Re: About convergence

Post by lesheng » Wed Mar 13, 2013 6:49 pm

Daniele Varsano wrote:Dear Lesheng,
the heaviness of the convergence studies depends totally on your system size (number of electrons, k-points etc..).
Most of them are independent so you can push a variable keeping the others low (but reasonable) in order to assess the each converged
value and finally run the production calculations.

Cheers,

Daniele
Thank you Daniele, and others in this forum, I find very helpful information during reading these posts.

Thanks,

Lesheng Li
PhD candidate
209 South Rd, Chapel Hill, NC 27599
Department of Chemistry
The University of North Carolina at Chapel Hill
Email: leshengl@live.unc.edu

User avatar
claudio
Posts: 456
Joined: Tue Mar 31, 2009 11:33 pm
Location: Marseille
Contact:

Re: About convergence

Post by claudio » Wed Mar 13, 2013 7:54 pm

Dear Lesheng,

moreover I advise you to use mHa as unit for the different parameters instead of RL (reciprocal lattice vectors = numer of g-vectors).
The mHa unit is in some sense less system dependent respect to vacuuum etc..

Cla
Claudio Attaccalite
[CNRS/ Aix-Marseille Université/ CINaM laborarory / TSN department
Campus de Luminy – Case 913
13288 MARSEILLE Cedex 09
web site: http://www.attaccalite.com

lesheng
Posts: 26
Joined: Thu Feb 07, 2013 11:39 pm

Re: About convergence

Post by lesheng » Tue Mar 26, 2013 2:35 am

claudio wrote:Dear Lesheng,

moreover I advise you to use mHa as unit for the different parameters instead of RL (reciprocal lattice vectors = numer of g-vectors).
The mHa unit is in some sense less system dependent respect to vacuuum etc..

Cla
Dear Cla,

Sorry I just notice your reply, and thank you for your advice, I'll try this way and compare the difference. :-)

Thank you again.

Best,
lesheng
PhD candidate
209 South Rd, Chapel Hill, NC 27599
Department of Chemistry
The University of North Carolina at Chapel Hill
Email: leshengl@live.unc.edu

Post Reply