Hi all,
I am wondering which unit should I use during the calculation. Since in the input file I generated, it uses RL as the unit for EXXRLvcs and NGsBlkXp, but in the first tutorial (Basic concepts of the GW approximation) in Yambo website, it says use Ry instead. So what is the difference?
Thanks for all your attention.
Best,
lesheng
About the unit in Exxrlvcs
Moderators: Davide Sangalli, andrea.ferretti, myrta gruning, andrea marini, Daniele Varsano
-
- Posts: 26
- Joined: Thu Feb 07, 2013 11:39 pm
About the unit in Exxrlvcs
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
209 South Rd, Chapel Hill, NC 27599
Department of Chemistry
The University of North Carolina at Chapel Hill
Email: leshengl@live.unc.edu
- Daniele Varsano
- Posts: 3816
- Joined: Tue Mar 17, 2009 2:23 pm
- Contact:
Re: About the unit in Exxrlvcs
Dear Lesheng,
The G vectors can be set in both ways. In RL unit (Reciprocal Lattice unit) you put in your input the numbers of Gvectors you want to use in your calculation. In this way have a control on the "size" of our calculation, while using energy units Ry, all the Gvectors sucht that G^2/2 < E(Ry) are included. In this way the number will depend on the cell size (for instances on the vacuum included in a non 3d periodic calculations), so it is useful when looking at volume effects convergences. You can use both ways, as you feel more confortable. Please note that in the report file you have the correspondence of the shell included for the first and last G-vectors:
Moreover when using energy units, you will also have the corresponding number of Gvectors in the generated databases showed in the report files.
Hope it helps,
Daniele
The G vectors can be set in both ways. In RL unit (Reciprocal Lattice unit) you put in your input the numbers of Gvectors you want to use in your calculation. In this way have a control on the "size" of our calculation, while using energy units Ry, all the Gvectors sucht that G^2/2 < E(Ry) are included. In this way the number will depend on the cell size (for instances on the vacuum included in a non 3d periodic calculations), so it is useful when looking at volume effects convergences. You can use both ways, as you feel more confortable. Please note that in the report file you have the correspondence of the shell included for the first and last G-vectors:
Code: Select all
[02.03] RL shells
=================
Shells, format: [S#] G_RL(mHa)
Hope it helps,
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/
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/