Re: [Wien] error dstart in sm-structure in hexagonal axes

2014-02-09 Thread Peter Blaha

Wien2k will by default always use the smallest possible unit cell.

Eg. for a fcc Cu cell (F -attice), you have only one atom (0,0,0), but if you
want to run this in a simple cubic lattice (P lattice, identical lattice 
parameters),
you would have 4 atoms (+(0.5,0.5,0) and permutations).
However, w2k does NOT allow such a calculation with 4 EQUIVALENT atoms.
If you still want to use a P cell with FOUR atoms, you have to make them
NON_EQUIVALENT. You can do this eg. by naming the atoms  Cu1, Cu2, Cu3 and Cu4.

Now the same happens for your problem. The R cell is 3 times smaller than the
H setting. If you want to use the H lattice, you must make some/all atoms
non-equivalent by the method explained above.

PS: Do you know why you want the H setting ?
It makes only sense for E-tot comparison with different phases (in real 
H-cells)
or when you want to do impurities, vacancies,... (although than the cell is 
still
too small).

Am 08.02.2014 12:14, schrieb Stefaan Cottenier:



I want to run sm metal in hexagonal axes, as Prof. Cottenier said, the
file struct should be as below;


You asked for the hexagonal information -- which you got. This does define the 
correct lattice, and wien2k can read/understand this, but cannot work with such 
a structure
file as such. As said before (see usersguide), for rhombohedral lattices you 
have to give the hexagonal lattice constant but rhombohedral coordinates. 
During init_lapw,
sgroup will make this conversion for you. If you accept its suggestion (which 
has only 3 atoms), dstart will run fine.

Stefaan




H   LATTICE,NONEQUIV.ATOMS   2  166 R-3m
MODE OF CALC=RELA unit=bohr
   6.856115  6.856115 49.573186 90.00 90.00120.00
ATOM  -1: X=0. Y=0. Z=0.
   MULT= 3  ISPLIT=15
   -1: X=0.6667 Y=0. Z=0.
   -1: X=0. Y=0.6667 Z=0.6667
Sm NPT=  781  R0=.1 RMT=   2.5   Z:  62.0
LOCAL ROT MATRIX:1.000 0.000 0.000
  0.000 1.000 0.000
  0.000 0.000 1.000
ATOM  -2: X=0. Y=0. Z=0.
   MULT= 6  ISPLIT=15
   -2: X=0.6667 Y=0. Z=0.5556
   -2: X=0. Y=0.6667 Z=0.8889
   -2: X=0. Y=0. Z=0.7778
   -2: X=0.6667 Y=0. Z=0.
   -2: X=0. Y=0.6667 Z=0.
Sm NPT=  781  R0=.1 RMT=   2.5   Z:  62.0
LOCAL ROT MATRIX:1.000 0.000 0.000
  0.000 1.000 0.000
  0.000 0.000 1.000
0  NUMBER OF SYMMETRY OPERATIONS


I make this structure, but when i initialized it, the new error in
dstart is obsearved as
Error in DSTART
  'ROTDEF' - no symmetry operation found.
  'ROTDEF' - for jatom, index   1   2
  'ROTDEF' - atomposition of jatom   0.000   0.000   0.000
  'ROTDEF' - atomposition of index   0.667   0.333   0.333
Please let me know, how can i solve it?


___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html



___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html


--
-
Peter Blaha
Inst. Materials Chemistry, TU Vienna
Getreidemarkt 9, A-1060 Vienna, Austria
Tel: +43-1-5880115671
Fax: +43-1-5880115698
email: pbl...@theochem.tuwien.ac.at
-
___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html


Re: [Wien] error dstart in sm-structure in hexagonal axes

2014-02-08 Thread Stefaan Cottenier



I want to run sm metal in hexagonal axes, as Prof. Cottenier said, the
file struct should be as below;


You asked for the hexagonal information -- which you got. This does 
define the correct lattice, and wien2k can read/understand this, but 
cannot work with such a structure file as such. As said before (see 
usersguide), for rhombohedral lattices you have to give the hexagonal 
lattice constant but rhombohedral coordinates. During init_lapw, sgroup 
will make this conversion for you. If you accept its suggestion (which 
has only 3 atoms), dstart will run fine.


Stefaan




H   LATTICE,NONEQUIV.ATOMS   2  166 R-3m
MODE OF CALC=RELA unit=bohr
   6.856115  6.856115 49.573186 90.00 90.00120.00
ATOM  -1: X=0. Y=0. Z=0.
   MULT= 3  ISPLIT=15
   -1: X=0.6667 Y=0. Z=0.
   -1: X=0. Y=0.6667 Z=0.6667
Sm NPT=  781  R0=.1 RMT=   2.5   Z:  62.0
LOCAL ROT MATRIX:1.000 0.000 0.000
  0.000 1.000 0.000
  0.000 0.000 1.000
ATOM  -2: X=0. Y=0. Z=0.
   MULT= 6  ISPLIT=15
   -2: X=0.6667 Y=0. Z=0.5556
   -2: X=0. Y=0.6667 Z=0.8889
   -2: X=0. Y=0. Z=0.7778
   -2: X=0.6667 Y=0. Z=0.
   -2: X=0. Y=0.6667 Z=0.
Sm NPT=  781  R0=.1 RMT=   2.5   Z:  62.0
LOCAL ROT MATRIX:1.000 0.000 0.000
  0.000 1.000 0.000
  0.000 0.000 1.000
0  NUMBER OF SYMMETRY OPERATIONS


I make this structure, but when i initialized it, the new error in
dstart is obsearved as
Error in DSTART
  'ROTDEF' - no symmetry operation found.
  'ROTDEF' - for jatom, index   1   2
  'ROTDEF' - atomposition of jatom   0.000   0.000   0.000
  'ROTDEF' - atomposition of index   0.667   0.333   0.333
Please let me know, how can i solve it?


___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html



___
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:  
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html