I have a question about the k-point that is given by vasp.
I found some k-points would not include the Brillouin zone.
For example, a hexagonal crystal system, kpoints 4x4x4, and the reciprocal crystal length (b1, b2, b3), got 20 kpoints(obtained in OUTCAR) through symmetry operator.
{this crystal character: b1=b2, angle between b1 and b2 is 60 degree}
Regardless of z axis , this k-points on x-y plane spread in a parallelogram b1xb2 (only consider xy plane), but BZ of this crystal shows hexagon shape. (3/8,3/8, 1/8) is one of this case that don't locate in Brullouin zone.
I want to confirm whether these kpoints out of BZ would fold into BZ during vasp calculation or not.
If they fold into BZ, the k+G numbers will be wrong. that would results in the mistake result in my calculation.
If not, they keep the same value, and proceed the calculation, k+G number will be right.
Does anyone know the answer of this question?
Thanks a lot.
K point & Brillouin zone
Moderators: Global Moderator, Moderator
-
- Newbie
- Posts: 2
- Joined: Thu Sep 03, 2009 9:10 am
K point & Brillouin zone
Last edited by chenmy on Thu May 27, 2010 8:09 pm, edited 1 time in total.
-
- Sr. Member
- Posts: 339
- Joined: Mon Apr 24, 2006 9:07 am
- License Nr.: 173
- Location: Gothenburg, Sweden
K point & Brillouin zone
I'm not sure I understand your question. Nevertheless, please notice that an even mesh breaks the full symmetry of a hexagonal cell. The subsequent symmetrization and reduction to the irreducible Brillouin zone then provides the k-point grid that is used in the calculation. Please compare with page 30 in
http://cms.mpi.univie.ac.at/vasp-worksh ... points.pdf
A gamma centered grid however preserves the full symmetry of the cell and also converges significantly faster than the even grid.
http://cms.mpi.univie.ac.at/vasp/vasp/h ... tices.html
The recommendation is therefore to use a gamma centered grid for a hexagonal cell.
Cheers,
/Dan
<span class='smallblacktext'>[ Edited Thu May 27 2010, 10:48PM ]</span>
http://cms.mpi.univie.ac.at/vasp-worksh ... points.pdf
A gamma centered grid however preserves the full symmetry of the cell and also converges significantly faster than the even grid.
http://cms.mpi.univie.ac.at/vasp/vasp/h ... tices.html
The recommendation is therefore to use a gamma centered grid for a hexagonal cell.
Cheers,
/Dan
<span class='smallblacktext'>[ Edited Thu May 27 2010, 10:48PM ]</span>
Last edited by forsdan on Thu May 27, 2010 8:33 pm, edited 1 time in total.
- oamnehc
- Newbie
- Posts: 5
- Joined: Fri May 21, 2010 5:45 am
- License Nr.: 5-269
- Location: Taipei
K point & Brillouin zone
I'd like to re-state the question raised by chenmy, for I have the similar problem.
Indeed it is natural and correct that by using a 4x4x4 k-mesh input (which is already a gamma-centered grid) for a hexagonal cell (the angle between a1 and a2 is pi/3), one gets 20 k points after symmetrization.
The problem is, why, among the 20 k points, are there 2 points described in VASP OUTCAR as (3/8,3/8,1/8) and (3/8,3/8,3/8)? These two points are OUT OF the conventionally defined Brillouin zone, and actually SHOULD be described by (3/8,-5/8,1/8) and (3/8,-5/8,3/8) instead. They should sit on the BZ boundary.
This is not just a matter of different descriptions. When it comes to any follow-up calculations using k+G, one would wonder if the confusing description of the k points is on purpose or not.
Best,
oamnehc
Indeed it is natural and correct that by using a 4x4x4 k-mesh input (which is already a gamma-centered grid) for a hexagonal cell (the angle between a1 and a2 is pi/3), one gets 20 k points after symmetrization.
The problem is, why, among the 20 k points, are there 2 points described in VASP OUTCAR as (3/8,3/8,1/8) and (3/8,3/8,3/8)? These two points are OUT OF the conventionally defined Brillouin zone, and actually SHOULD be described by (3/8,-5/8,1/8) and (3/8,-5/8,3/8) instead. They should sit on the BZ boundary.
This is not just a matter of different descriptions. When it comes to any follow-up calculations using k+G, one would wonder if the confusing description of the k points is on purpose or not.
Best,
oamnehc
Last edited by oamnehc on Fri Jun 04, 2010 3:43 am, edited 1 time in total.