Hello,
Given a hexagonal cell and a kpoint set that doesn't have the full symmetry of the lattice I see the vasp code applying the symmetry operations of the full symmetry group to each kpoint before reducing it with the symmetry of the lattice. ( in function IBZKPT in symmetry.F). This leads to a non-equidistant kpoint set with significantly more kpoints in the irreducible zone than maximally expected from the choice of the monkhorst pack parameters.
Despite the explicit discouragement of this choice of kpoint set, I wonder why the above stated is necessary, particularly because only the lattice symmetry is used in the symmetrization of the charge density.
What consequences would not enforcing the full symmetry group have?
Thanks in advance
kpoints and symmetrization
Moderators: Global Moderator, Moderator
-
- Newbie
- Posts: 1
- Joined: Fri May 23, 2008 11:50 pm
- License Nr.: 516
kpoints and symmetrization
Last edited by atheb on Wed Feb 18, 2009 4:17 am, edited 1 time in total.
-
- Administrator
- Posts: 2921
- Joined: Tue Aug 03, 2004 8:18 am
- License Nr.: 458
kpoints and symmetrization
the consequence is a (significant) increase in computation time (depending on how large the fraction of the IBZ is of the full BZ).
especially for k-meshes of non-rectangular lattices the symmetrization of the k-mesh is sometimes not easily done by vasp if the mesh does not explicitely include Gamma (typically for even meshes, please have a look at the k-points.pdf tutorial file of the vasp-workshop). Please try if the behavior pesists with an odd k-mesh as well
especially for k-meshes of non-rectangular lattices the symmetrization of the k-mesh is sometimes not easily done by vasp if the mesh does not explicitely include Gamma (typically for even meshes, please have a look at the k-points.pdf tutorial file of the vasp-workshop). Please try if the behavior pesists with an odd k-mesh as well
Last edited by admin on Thu Feb 19, 2009 3:47 pm, edited 1 time in total.