Re: [Wien] Problem with structeditor for slabs

2014-02-21 Thread Michael Sluydts

Thanks Robert,

All slabs have now generated succesfully without any holes.


Kind regards,

Michael Sluydts


Robert Laskowski schreef op 21/02/2014 3:18:


Hi,

cp attached file to $WIENROOT/SRC_structeditor/bin.

thanks

Robert

On 20 February 2014 PM 10:17:55 Michael Sluydts wrote:

Hello everyone,

I'm trying to generate a series of (110) surface slabs of varying 
thickness for a zincblende crystal. I'm doing this using the 
structeditor functions in an octave script (and have done this with 
both wien2k (11), 12 and 13).


At certain thicknesses of the slab a problem occurs, which I've seen 
in the past for other systems, where it appears to make a diagonal cut 
at the bottom of the slab along a {001} direction. Does anyone by any 
chance have a patch to resolve this? I know ways around it myself but 
being able to use the provided functions directly would be nice.


I've added the bulk struct and octave script needed to reproduce the 
problem below.



Thanks,

Michael Sluydts

bulk.struct
blebleble
F   LATTICE,NONEQUIV.ATOMS   2  216 F-43m
MODE OF CALC=RELA unit=bohr
 10.847217 10.847217 10.847217 90.00 90.00 90.00
ATOM  -1: X=0. Y=0. Z=0.
  MULT= 1  ISPLIT=15
Zn NPT=  781  R0=.5 RMT=   2.0   Z: 30.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.2500 Y=0.2500 Z=0.2500
  MULT= 1  ISPLIT=15
Se NPT=  781  R0=.5 RMT=   2.0   Z: 34.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

script.m
bulk = makeconventional(loadstruct('bulk.struct'));
surf = makesurface(bulk,[1 1 0],2,(bulk.a(3)/2*sin(pi/4))*10+0.01,37.8);
savestruct(makeprimitive(surf),['bug.struct'],0);



--

==

Dr. Robert Laskowski

Senior Scientist, Materials Science  Engineering Department

Institute of High Performance Computing, A*STAR

1 Fusionopolis Way, #16-16, Connexis, Singapore 138632

Tel(Off): +65. 64191493 Fax: +65. 64632536

=


IHPC Values :: Impact :: Honesty :: Performance :: Co-operation
This email is confidential and may be privileged. If you are not the 
intended recipient, please delete it and notify us immediately. Please 
do not copy or use it for any purpose, or disclose its contents to any 
other person. Thank you.



___
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


[Wien] steucture phases stability

2014-02-21 Thread MAHDI SALMANI HIRMAND
Dear Wien2k users,
How can i find structure phases stability of a compound in two symmetries?
How can i find the pressure related to it?
Is it necessary to use Gibbs energy for finding the pressure which symmetry of 
compound will change?

Please give me some advices.
With best,
mahdi
___
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] steucture phases stability

2014-02-21 Thread Oleg Rubel
Please, check this text:
http://onlinelibrary.wiley.com/book/10.1002/9780470447710
Chapter 7 Equilibrium Phase Diagrams from Ab Initio Thermodynamics (pages
163-177) can be a good starting point.

Oleg


On Fri, Feb 21, 2014 at 7:51 AM, MAHDI SALMANI HIRMAND 
mahdisa_ri...@yahoo.com wrote:

 Dear Wien2k users,
 How can i find structure phases stability of a compound in two symmetries?
 How can i find the pressure related to it?
 Is it necessary to use Gibbs energy for finding the pressure which
symmetry of compound will change?
 Please give me some advices.
 With best,
 mahdi

 ___
 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


[Wien] Error in lapw2 during DOS

2014-02-21 Thread Lawal Mohammed
Dear Developers and Users,

We are trying to improve our calculation by using MBJ potential. Electron 
charge density finished successfully,as we run x lapw2 -qtl for DOS this error 
appears.

forrtl: severe (24): end-of-file during read, unit 10, file 
/root/WIEN2k/mbj/./mbj.vector
Image  PCRoutineLineSource  
   
lapw2  0053BE1A  Unknown   Unknown  Unknown
lapw2  0053A995  Unknown   Unknown  Unknown
lapw2  004DF6B6  Unknown   Unknown  Unknown
lapw2  0049C476  Unknown   Unknown  Unknown
lapw2  0049BBE9  Unknown   Unknown  Unknown
lapw2  004B7DD4  Unknown   Unknown  Unknown
lapw2  0045E0D5  l2main_   618  
l2main_tmp_.F
lapw2  00474CAC  MAIN__564  lapw2_tmp_.F
lapw2  00403CEC  Unknown   Unknown  Unknown
libc.so.6  00346442169D  Unknown   Unknown  Unknown
lapw2  00403BE9  Unknown   Unknown  Unknown
0.441u 0.025s 0:00.46 100.0%0+0k 0+3424io 0pf+0w
error: command   /home/mohammed/WIEN2K/lapw2 lapw2.def   failed 
We need help from you to overcome this problem.

With best regards
Lawal Mohammed
___
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