No subject


Wed Jul 6 14:13:09 CEST 2005


(for the primitives, ignoring ky for simplicity)
0 < g1 < 1 ; 0 < g3 < 1, with g1, g3 the primitive vectors. However, this
is not the irreproducible part of the BZ, which should be:
0 < g1 < 1 ; -g1 < g3 < g1
unless I'm being stupid (possible).

N.B., after the line "NUMBER OF POINT GROUP OPERATION" something looks
fishy, since I'm getting:
    SYMMETRY MATRIX NR.  1   SYMMETRY MATRIX NR.  2
    ***************          ***************
    ***************          ***************
    ***************          ***************

On Tue, 29 Jun 2004, Jorissen Kevin wrote:

> To come back to this question :
>
> kgen is a bit of a mess, made up of Peter Blöchl's programs and parts added to implement his code into some old version of the wien code.  I definitely agree that eg. the symmetry operations are taken from case.struct, and not from the symmetry dedicated routines in kgen.  However, I thought that the reciprocal unit cell (which is to be sampled by a uniform mesh and then reduced by symmetry) is calculated by some of the old routines, and these indeed reduce some systems to primitive cells, don't they?
>
> Let's go into main.f, and take the example of an F lattice with basis a * (1 0 0), (0 1 0), (0 0 1).  The lattice is read from case.struct.  It is fed into bravai.f, which redefines the basis as a * (1/2, 1/2, 0) etc.  THis is in the variable rbas, which is later used (in main and in arbmsh) by the routine gbas to calculate the reciprocal basis.  This basis defines the reciprocal cell which will be sampled by basdiv, zuord ...
>
>
> This is the reason why I referred to the routine bravai.f.
>
> Please correct me if I am mistaken.
>
> Kevin.
>
>
> 	-----Oorspronkelijk bericht-----
> 	Van: Peter Blaha [mailto:pblaha at zeus.theochem.tuwien.ac.at]
> 	Verzonden: ma 6/28/2004 4:14
> 	Aan: wien at zeus.theochem.tuwien.ac.at
> 	CC:
> 	Onderwerp: RE: [Wien] k-point sampling
>
>
>
> 	> Kgen may indeed change your lattice.  Check out SRC_kgen/bravai.f to get an idea of what it's doing in your case.
>
> 	I havn't checked this now, but what I remember is that bravai.f is still
> 	present (and also called), but actually more or less NOT used !!!
> 	(most variables (eg. symmetry operations,..) are reset to values from
> 	case.struct,....)
>
> 	                                      P.Blaha
> 	--------------------------------------------------------------------------
> 	Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
> 	Phone: +43-1-58801-15671             FAX: +43-1-58801-15698
> 	Email: blaha at theochem.tuwien.ac.at    WWW: http://info.tuwien.ac.at/theochem/
> 	--------------------------------------------------------------------------
>
> 	_______________________________________________
> 	Wien mailing list
> 	Wien at zeus.theochem.tuwien.ac.at
> 	http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
>
>
>
>

-----------------------------------------------
Laurence Marks
Department of Materials Science and Engineering
MSE Rm 2036 Cook Hall
2225 N Campus Drive
Northwestern University
Evanston, IL 60201, USA
Tel: (847) 491-3996 Fax: (847) 491-7820
mailto:L-marks at northwestern.edu
http://www.numis.northwestern.edu
-----------------------------------------------





More information about the Wien mailing list