[Wien] The calculation principles of X-ray absorption spectra by wien2k

2008-09-10 Thread Lily Anh
Dear all,
I know that the x-ray absorption spectra are calculated for the dipole allowed 
transitions, generating matrix elements, which are multipliedwith a radial 
transition probability and the partial densities of states as explained in the 
userguide.

But I want to know more detailed principles employed in Wien2k. However, the 
references in userguide seem too old to find out.
 
Here are the references mentioned in userguide:
Schwarz K., Neckel A and Nordgren J, J.Phys.F:Metal Phys. 9, 2509 (1979)
Schwarz K., andWimmer E, J.Phys.F:Metal Phys. 10, 1001 (1980)
Neckel A., Schwarz K., Eibler R. and Rastl P. 1975 Microchim.Acta, Suppl.6, 257.

 Can anybody help me about the references above or recommend some other 
references about the principles used in Wien2k?

I will highly appreciate your help.
Best regards. 
LiLy
lily_physics at yahoo.com.sg


  New Email addresses available on Yahoo!
Get the Email name you've always wanted on the new @ymail and @rocketmail. 
Hurry before someone else does!
http://mail.promotions.yahoo.com/newdomains/sg/
-- next part --
An HTML attachment was scrubbed...
URL: 
http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20080910/2793be43/attachment.html


[Wien] Subscript out of range warning

2008-09-10 Thread Yurko Natanzon
Dear Prof. Blaha,
Thank you for your explanation. The problem was that the forces were
too large. I performed the force minimization with min_lapw and the
warning disappeared.

regards,
Yurko

2008/8/11 Peter Blaha :
> atest   is used in  run_lapw to find out whether a certain atom has
> forces or not.
>
> It might point to another incompatibility of Linux/csh/tcsh versions,
> which I cannot test.
>
> If you are somehow familiar with Linux, you may edit   $WIENROOT/run_lapw
> search for  "atest":
> set atest=`head -$itest $file.inM |tail -1`
> set itest=`echo " $atest[1] + $atest[2] + $atest[3]"|bc`
>
> and print some variables (echo $itest   before the two lines and
> echo $atestafter the first line)
>
> This may give a hint to solve the problem.
>
>
>
> Yurko Natanzon schrieb:
>> Dear wien2k users,
>> I'm running scf calculation of TiO2 (rutile) with GGA-PBE functional
>> and other parameters as in the example section of the User manual.
>> After each lapw0 cykle I get a warning:
>>
>> atest: Subscript out of range
>>
>> However, the calculation continues and finishes normally, and no
>> .error files are created. This warning appears twice after succecfull
>> LAPW0 end in each SCF cycle.
>>
>> Could you explain the meaning of this warning and what can be done to
>> avoid it? Seems that nobody have posted such question before, the
>> maillist search gave no results.
>>
>>
>
> --
> -
> Peter Blaha
> Inst. Materials Chemistry, TU Vienna
> Getreidemarkt 9, A-1060 Vienna, Austria
> Tel: +43-1-5880115671
> Fax: +43-1-5880115698
> email: pblaha at theochem.tuwien.ac.at
> -
> ___
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
>



-- 
Yurko (aka Yuriy, Iurii, Jurij etc) Natanzon
PhD Student
Henryk Niewodniczan`ski Institute of Nuclear Physics
Polish Academy of Sciences
ul. Radzikowskiego 152,
31-342 Krako`w, Poland
Email: Yurii.Natanzon at ifj.edu.pl, yurko.natanzon at gmail.com


[Wien] d5/2 and d3/2 states of tungsten

2008-09-10 Thread Md. Fhokrul Islam

Dear Wien2k users,
 
 I am trying to calculate the density of states of d5/2 and d3/2 states of 
tungsten using qtl package of wien2k following the instructions manual and have 
run 
 
run_lapw -p -i 200 -ec 0.1 -cc 0.001 -so
x qtl -p -so
 
It generates case.inq file as follows
 
-9.0   3.0   Emin  Emax   1   number of atoms   1  -1  
1  0   iatom,qsplit,symmetrize,locrot3   0  1  2  nL, l-values
 
but resulting case.qtl file gives only zeros (a part of it is shown below). I 
am not sure where I am going wrong. Any suggestions will be appreciated.
 
Thanks,
fislam
 
 
 LATTICE CONST.=  5.9814  5.9814  5.9814   FERMI ENERGY=   1.02990  102 < NMAT 
<  118   SPIN=1   NAT=  1  SO 0 KLmax 11 JATOM  1  MULT= 1  ISPLIT=-1  
tot,s,p,p1/2(-1/2),p1/2(1/2),p3/2(-3/2),,,p3/2(3/2),d,d3/2(-3/2),,,d3/2(3/2),d5/2(-5/2),d5/2(5/2),
 
  -0.95643  1 0.00.0 0.0 0.0 0.0 0.0 0.0 
0.0 0.0 0.0 0.0 0.0  -0.95643  2 1.0  -0.95645  1 
0.00.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 
0.0 0.0 0.0  -0.95645  2 1.0 BAND:  23   0.33609  1 0.0
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 
0.0   0.33609  2 1.0   0.36760  1 0.00.0 0.0 0.0 
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0   0.36760  2 
1.0
 
_
Connect to the next generation of MSN Messenger?
http://imagine-msn.com/messenger/launch80/default.aspx?locale=en-us&source=wlmailtagline
-- next part --
An HTML attachment was scrubbed...
URL: 
http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20080910/8ee32479/attachment.html