<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">

<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>

</head>

<body dir="ltr">

<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">

<p><span id="ms-rterangepaste-start"></span></p>

<p style="margin-bottom:0in; line-height:100%">Dear WIEN2K users and developers,</p>

<p style="margin-bottom:0in; line-height:100%"><br>

</p>

<p style="margin-bottom:0in; line-height:100%">I am currently playing with the wavefunctions in WIEN2K in order to compute material-dependent properties. I must then extract them !

</p>

<p style="margin-bottom:0in; line-height:100%"><br>

</p>

<p style="margin-bottom:0in; line-height:100%">I managed to extract the relevant quantities (C_nk(G), the coefficents A_lm, B_lm, C_lm and the radial functions) in order to build the wavefunction. When I try to compute the norm of the wavefunction, I find that

 it is quite well normalized BUT only for states which are below the Fermi energy Е For states above, discrepancies arise and I find numbers like 0.8 or 1.3 Е</p>

<p style="margin-bottom:0in; line-height:100%"><br>

</p>

<p style="margin-bottom:0in; line-height:100%">To output the A_lm, B_lm of unoccupied states with "x lapw2 -alm", I modified slightly the l2main.F routine. The following has been done :

</p>

<p style="margin-bottom:0in; line-height:100%">IF(MODUS.EQ.'ALM ') then</p>

<p style="margin-bottom:0in; line-height:100%">!LEANDRO INCLUDE EMPTY BANDS (START)</p>

<p style="margin-bottom:0in; line-height:100%">NEMAX_SAVE=NEMAX</p>

<p style="margin-bottom:0in; line-height:100%">NEMAX=NE</p>

<p style="margin-bottom:0in; line-height:100%">!LEANDRO INCLUDE EMPTY BANDS (END)</p>

<p style="margin-bottom:0in; line-height:100%">WRITE(24,2055) s_kvec,t_kvec,z_kvec,n,ne,bname</p>

<p style="margin-bottom:0in; line-height:100%">write(24,*) jatom,nemin,nemax,' jatom,nemin,nemax'</p>

<p style="margin-bottom:0in; line-height:100%">endif</p>

<p style="margin-bottom:0in; line-height:100%"><br>

</p>

<p style="margin-bottom:0in; line-height:100%">where my modification is in between the У!LEANDRO ЕФ. What I did it just saying that the NEMAX (which normally refers to the number of occupied bands) should go up to the highest computed state and thus, can go

 above the Fermi level.</p>

<p style="margin-bottom:0in; line-height:100%"><br>

</p>

<p style="margin-bottom:0in; line-height:100%">I was wondering that may be, the states are automatically normalized but only for those below the Fermi level. If this is the case, then I can compute the norm and divide by the sqrt. Am I right or am I missing

 something ? Since the wavefunction is the basic mathematical description of the material, extracting wrong quantities would be quite problematic ...</p>

<p style="margin-bottom:0in; line-height:100%">I am of course taking into account the local orbitals in the process !<br>

<br>

Are the wavefunctions only normalized below E_F ? If yes, do you know in which part the normalization is done ?</p>

<p style="margin-bottom:0in; line-height:100%"><br>

</p>

<p style="margin-bottom:0in; line-height:100%">If anyone has experience with this topic or any suggestion, I would be please to hear ! I have already checked throughout the mailing list and the user guide ...</p>

<p style="margin-bottom:0in; line-height:100%"><br>

</p>

<p style="margin-bottom:0in; line-height:100%">Thank you,</p>

<p style="margin-bottom:0in; line-height:100%"><br>

</p>

<p style="margin-bottom:0in; line-height:100%">Leandro Salemi</p>

<span id="ms-rterangepaste-end"></span><br>

<p></p>

</div>

</body>

</html>