[Wien] Projection of Wavefunctions on bessel plane

2013-03-21 Thread محمد ارشد فرحان
Dear Wien users/developers,
hi,

i'm trying to figure out bands arising from surface states in the band
structure of a slab geometry.

one method reported in
http://iopscience.iop.org/0953-8984/24/3/035502using VASP code is by
projecting the wave function on top/bottom
plane/layer of the unit cell (slab)

can any body suggest me how can i do it in Wien2k.

Regards,

-- 
*M.Arshad Farhan*
Dept of Chemistry,
Pohang Univ of Sci & Tech
Pohang, Republic of Korea
-- next part --
An HTML attachment was scrubbed...
URL: 
<http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20130321/18a05a0a/attachment.htm>


[Wien] Tr : elastic

2013-03-21 Thread Jameson Maibam
Use the elastic calculation packages available in the wien2k unsupported goodies

--
Sent from my Nokia phone

--Original message--
From: Mourad Karima 
To: "wien at zeus.theochem.tuwien.ac.at" 
List-Post: wien@zeus.theochem.tuwien.ac.at
Date: Wednesday, March 20, 2013 7:14:35 PM UTC
Subject: [Wien] Tr : elastic








Respected Prof.PBalha sir and all wien2k users



how to install the program genetempl to calculate the properties? elastic


-
thanks

___
Wien mailing list
Wien at zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien



[Wien] QTL-B message in scf2 after "x lapw2 -qtl"

2013-03-21 Thread Sanae Fujita
Dear Prof.Blaha: Thank you very much for your kind help.

> Always modify the energy parameters of the atom and l, where qtl-bs occur.

I understand. Thank you.  (The purpose of my strange procedure was to check
it by changing atom2,l=0's energy parameter 2.5-4Ry..)



> And small values like 2.x for states at high energy are probably ok 

After my iterative process, The lowest QTL-B value I got is 2.94021 in
1.22167Ry of atom2 l=1. Its Fermi energy is about 0.29Ry.
Is 0.93Ry above Fermi energy high enough?
If so, I can get the satisfaction of this result!



> Why case.in1_st ??? It must be case.in1 

I thought that rewiting case.in1_st was necessary in order to change
energy parameter for SCF.

But with your advice, I checked and knew that rewriting case.in1 before SCF
is enough. Thank you.


Best regards,
S.Fujita
-- next part --
An HTML attachment was scrubbed...
URL: 
<http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20130321/bc2ed93f/attachment.htm>


[Wien] IBM AIX error

2013-03-21 Thread Peter Blaha
I don't think you can use ksh or bsh, but have to use  csh (or tcsh)

I suggest to set in x_lapw

#!/bin/csh -xf

or

#!/bin/tcsh -xf  (works only, if you have a tcsh)

On some systems csh and tcsh are different, on others the same.

The -x   will give you an even longer output when runningx lapw2
and you should find the values of eferm, efold, ... used in the crucial 
line:

 > set eferm=`grep ':FER  :' $file.scf2$sc$updn|cut -f 2 -d=`
 > set efold=`head -1 $file.in1$cmplx${sc}|cut -c11-17 | sed -e
 > "s/[a-zA-Z(,=]//g"`
 > if ( $#efold == 0 ) set efold=0.5

 > set eferm=`echo " if($eferm > $efold) {if ($eferm > $efold + 0.2)
 > {$efold + 0.2} else {$eferm} } else {if($eferm < $efold - 0.2) 
{$efold - 0.2} else ($eferm*0.75+$efold*0.25) }" | bc -l`

Once you know that $eferm,.. are ok (reasonable numbers) you may want to 
test it in a terminal window:

/bin/csh  (or tcsh) (to define the csh
set eferm=0.55
set efold=0.45
echo " if($eferm > $efold) {if ($eferm > $efold + 0.2)
 > {$efold + 0.2} else {$eferm} } else {if($eferm < $efold - 0.2) 
{$efold - 0.2} else ($eferm*0.75+$efold*0.25) }" | bc -l

This should give you the same error (> syntax error on line 1 stdin)

then remove the | bc -l

any error ? If not, the problem is in the bc program  (maybe this 
version does not understand the nested "if" statements ? Then you would 
have to program this in csh.


On 03/20/2013 11:18 PM, Oliver Albertini wrote:
> Here is the relevant info from adding the verbose flag to the x_lapw file:
>
>
> set lapw2test=`echo $exe | cut -c 1-5`
> if ( "$lapw2test" == 'lapw2' ) then
> if ( ! $?in1orig & ! $?band & ! $?emin1 & ! $?emax1 ) then
> if ( -e $file.scf2$sc$updn && ! -z $file.scf2$sc$updn ) then
> if ( $cmplx == 'c' ) then
> set eferm=`grep ':FER  :' $file.scf2$sc$updn|cut -f 2 -d=`
> set efold=`head -1 $file.in1$cmplx${sc}|cut -c11-17 | sed -e
> "s/[a-zA-Z(,=]//g"`
> if ( $#efold == 0 ) set efold=0.5
> set eferm=`echo " if($eferm > $efold) {if ($eferm > $efold + 0.2)
> {$efold + 0.2} else {$eferm} } else {if($eferm < $efold - 0.2) {$efold -
> 0.2} else ($eferm*0
> .75+$efold*0.25) }" | bc -l`
> syntax error on line 1 stdin
> set modus=`head -1 $file.in1$cmplx${sc} | cut -c1-5`
> echo "$modus  EF=$eferm   (WFFIL, WFPRI, ENFIL, SUPWF) " >
> $file.in1$cmplx${sc}_tmp
> tail -n +2 $file.in1$cmplx${sc} >> $file.in1$cmplx${sc}_tmp
> mv $file.in1$cmplx${sc}_tmp $file.in1$cmplx${sc}
> endif
> endif
> endif
>
>
>
> On Wed, Mar 20, 2013 at 3:13 PM, Oliver Albertini  > wrote:
>
> So after making a fresh initialization, the error remains the same.
> The file uplapw2.error is indeed empty. The default on AIX is ksh,
> so I tried both #!/bin/ksh -vp
> and #!/bin/bsh -vf. The -p for Kourne Shell means to ignore .profile.
>
>
> On Wed, Mar 20, 2013 at 12:58 PM, Laurence Marks
> mailto:L-marks at northwestern.edu>> wrote:
>
> Did you run from a fresh initialization or use the old files? It
> could
> be that something went wrong before so perhaps everything is OK.
>
> Is the file uplapw2.error empty? If it is then it is something after
> that in lapw2para. You can edit the first line to be
> #!/bin/csh -vf
>
> (assuming that your csh is standard) and then do x lapw2 -up.
> You will
> get a mess of output but the bottom will give you an idea where
> it is
> failing. (You might want to do something like "x lapw2 -up >> junk"
> and then look at the file junk, and/or a different piping of the
> output which will depend upon what shell you are using.)
>
> On Wed, Mar 20, 2013 at 2:49 PM, Oliver Albertini
> mailto:ora at georgetown.edu>> wrote:
>  > Thanks to you both,
>  >
>  > I am currently using a standard lapack library. You are
> right, the routines
>  > dggglm and dgelsy were not present in the SRC_lib version.
> Removing essl
>  > resulted in even more missing routines, so I tried replacing
> it with
>  > standard 64 bit BLAS. It compiled, and I tried with the new
> mixer, but the
>  > same messages appeared as before.
>  >
>  > I am not sure what the meaning of 'syntax error on line 1
> stdin' is. Yes it
>  > still happens when running x lapw2 -up.
>  >
>  >
>  >
>  >
>  > On Wed, Mar 20, 2013 at 1:07 AM, Peter Blaha
> mailto:pblaha at 
> theochem.tuwien.ac.at>>
>  > wrote:
>  >>
>  >> Unfortunately I do not have access to an AIX system, so
> cannot do anything
>  >> myself.
>  >>
>  >>
>  >>> I recently compiled 12.1 on AIX (v 6.1) pwr6. Like Luis, I
> also had to
>  >>> make some changes to SRC's in order to finish the
> compilation. These
>  >>> were mostly issues with xlf like syntax

[Wien] Projection of Wavefunctions on bessel plane

2013-03-21 Thread Peter Blaha
In wien2k this is very simple.  You plot bandstructures with "fat bands" 
and use the partial charges (qtl) to find out the character of each band 
(state).

Once you found in interesting looking state you can look directly in the 
partial charges of this state (case.qtl file).

A surface state is a state which has its charge predominantly at the 
surface-layer atoms and contributions from subsurface layers decay 
rapidly.  (A bulk state will have contributions from "bulk-like" atoms, 
often with oscillatory behavior as this corresponds to a certain 
k-perpendicular state.

On 03/20/2013 06:30 PM,   ? wrote:
> Dear Wien users/developers,
> hi,
> i'm trying to figure out bands arising from surface states in the band
> structure of a slab geometry.
> one method reported in http://iopscience.iop.org/0953-8984/24/3/035502
> using VASP code is by projecting the wave function on top/bottom
> plane/layer of the unit cell (slab)
> can any body suggest me how can i do it in Wien2k.
> Regards,
>
> --
> *_M.Arshad Farhan_*
> Dept of Chemistry,
> Pohang Univ of Sci & Tech
> Pohang, Republic of Korea
>
>
> ___
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
>

-- 

   P.Blaha
--
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
Phone: +43-1-58801-165300 FAX: +43-1-58801-165982
Email: blaha at theochem.tuwien.ac.atWWW: 
http://info.tuwien.ac.at/theochem/
--