<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></p>

<div class="msgHead" style="font-family:helvetica,arial,sans-serif; font-size:14px">

<p class="darkgray font13" style="font-size:13px; color:rgb(102,102,102)"></p>

<p style="margin-bottom:0in; line-height:100%"><tt class="western"><font color="#000000"><font face="courier, courier new, monospace"><font style="font-size:10pt">Dear Pr. Blaha,<br>

<br>

Thank you for your quick answer. You are totally right with the case.in2 file Е</font></font></font></tt></p>

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

</p>

<p style="margin-bottom:0in; line-height:100%"><tt class="western"><font color="#000000"><font face="courier, courier new, monospace"><font style="font-size:10pt">For the wavefunction, thatТs how I do :</font></font></font></tt></p>

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

</p>

<p style="margin-bottom:0in; line-height:100%"><tt class="western"><font color="#000000"><font face="courier, courier new, monospace"><font style="font-size:10pt">For the interstitial, I integrate the product of PW over the whole space and subtract the part

 from the spheres using a Rayleigh-expansion of a PW in terms of spherical harmonics.</font></font></font></tt></p>

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

</p>

<p style="margin-bottom:0in; line-height:100%"><tt class="western"><font color="#000000"><font face="courier, courier new, monospace"><font style="font-size:10pt">In the spheres, I compute terms like Alm* Blm \int u_l(r) dot(u)_l(r) r*r dr for the spheres.

 These coefficients are found in case.almblm and are actually Alm(k) = SUM_G C_G Alm(k+G) (this I got from a previous mailing thread). I have also terms like Alm* Clm \int u_l(r) u_lo__l(r) r*r dr to take into account the APW+lo and LO ! I guess that the ClmТs

 in case.almclm are also the ones multiplied by their coefficient C_G(LO) (which are not really part of the PW expansion), just like for Alm and Blm. I sum those terms over the different indexes l,m and over the spheres.</font></font></font></tt></p>

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

</p>

<p style="margin-bottom:0in; line-height:100%"><tt class="western"><font color="#000000"><font face="courier, courier new, monospace"><font style="font-size:10pt">Actually, I follow more less what is presented in the paper of C. Ambrosch-Draxl and J. O. Sofo

 (Computer Physics Communications 175 (2006) 1Ц14). They computed <Psi_mk| p |Psi_nk> so I of course adapted the formalism to my need.</font></font></font></tt></p>

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

</p>

<p style="margin-bottom:0in; line-height:100%"><tt class="western"><font color="#000000"><font face="courier, courier new, monospace"><font style="font-size:10pt">Since you said that wavefunctions where all normalized, even above E_F, I must have missed something

 somewhere Е If you see any mistake in what I have written there or have any references, works or source code I should look at, I would be pleased to hear that ! My aim is to extract exactly the wavefunctions. I will then manipulate them.</font></font></font></tt></p>

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

</p>

<p style="margin-bottom:0in; line-height:100%"><tt class="western"><font color="#000000"><font face="courier, courier new, monospace"><font style="font-size:10pt">Thank you very much,</font></font></font></tt></p>

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

</p>

<p style="margin-bottom:0in; line-height:100%"><tt class="western"><font color="#000000"><font face="courier, courier new, monospace"><font style="font-size:10pt">Best regards,</font></font></font></tt></p>

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

</p>

<p style="margin-bottom:0in; line-height:100%"><tt class="western"><font color="#000000"><font face="courier, courier new, monospace"><font style="font-size:10pt">Leandro Salemi</font></font></font></tt></p>

<p style="margin-bottom:0in; line-height:100%"><tt class="western"><font color="#000000"><font face="courier, courier new, monospace"><font style="font-size:10pt"><br>

</font></font></font></tt></p>

<p style="margin-bottom:0in; line-height:100%"><tt class="western"><font color="#000000"><font face="courier, courier new, monospace"><font style="font-size:10pt">P.S. : I might have answered in a not-proper way since I had not received the mail (I did not

 pay attention that my account was disabled for the reception of the mail ...). Sorry for the inconvenience (now it is enabled).</font></font></font></tt></p>

<br>

<p></p>

<p class="darkgray font13" style="font-size:13px; color:rgb(102,102,102)"><br>

</p>

<p class="darkgray font13" style="font-size:13px; color:rgb(102,102,102)"><br>

</p>

<p class="darkgray font13" style="font-size:13px; color:rgb(102,102,102)"><br>

</p>

<p class="darkgray font13" style="font-size:13px; color:rgb(102,102,102)"><span class="sender pipe" style="margin-right:0.5em; border-right:1px solid rgb(221,221,221); padding-right:0.8em"><span itemprop="author" itemscope="" itemtype="http://schema.org/Person" style="color:rgb(102,102,102)"><span itemprop="name"><a href="https://www.mail-archive.com/search?l=wien@zeus.theochem.tuwien.ac.at&q=from:%22Peter+Blaha%22" rel="nofollow" style="color:rgb(102,102,102)">Peter

 Blaha</a></span></span></span> <span class="date"><a href="https://www.mail-archive.com/search?l=wien@zeus.theochem.tuwien.ac.at&q=date:20171204" rel="nofollow" style="color:rgb(102,102,102)">Mon, 04 Dec 2017 06:13:01 -0800</a></span></p>

<p class="darkgray font13" style="font-size:13px; color:rgb(102,102,102)"><br>

</p>

</div>

<div itemprop="articleBody" class="msgBody" style="font-family:helvetica,arial,sans-serif; font-size:14px">

<tt style="margin:0px; font-family:courier,"courier new",monospace; font-size:1em">Of course all the wave functions are normalized. It comes automatically </tt><tt style="margin:0px; font-family:courier,"courier new",monospace; font-size:1em">after solving

 the generalized eigenvalue problem in lapw1. </tt><tt style="margin:0px; font-family:courier,"courier new",monospace; font-size:1em">Without properly normalized wave functions, how should we calculate an </tt><tt style="margin:0px; font-family:courier,"courier new",monospace; font-size:1em">electron

 density which yields the desired number of electrons ??</tt>

<pre style="font-family:courier,"courier new",monospace; font-size:1em; white-space:pre-wrap; word-wrap:break-word; margin:0em">Why did you modify the code ?

It would be much simpler to set NE in case.in2 to a larger number ....</pre>

<pre style="font-family:courier,"courier new",monospace; font-size:1em; white-space:pre-wrap; word-wrap:break-word; margin-top:0px; margin-bottom:0px"></pre>

<tt style="margin:0px; font-family:courier,"courier new",monospace; font-size:1em">I don't know how you actually calculated the norm of a wave function, </tt><tt style="margin:0px; font-family:courier,"courier new",monospace; font-size:1em">but due to the dual

 representation, it is not so straightforward.</tt></div>

<div itemprop="articleBody" class="msgBody" style="font-family:helvetica,arial,sans-serif; font-size:14px">

<br>

</div>

<div itemprop="articleBody" class="msgBody" style="font-size:14px"><tt style="font-family:courier,"courier new",monospace; margin:0px; font-size:1em"></tt>

<pre style="font-family:courier,"courier new",monospace; font-size:1em; white-space:pre-wrap; word-wrap:break-word; margin:0em">On 12/04/2017 02:55 PM, Leandro Salemi wrote:

</pre>

<blockquote style="font-family:helvetica,arial,sans-serif; margin:0em; padding:0px 0px 0px 0.85em; color:rgb(160,30,30); border-left:0.2em solid rgb(85,85,238)">

<pre style="font-family:courier,"courier new",monospace; font-size:1em; white-space:pre-wrap; word-wrap:break-word; margin:0em">Dear WIEN2K users and developers,

I am currently playing with the wavefunctions in WIEN2K in order to

compute material-dependent properties. I must then extract them !

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 Е

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 :

IF(MODUS.EQ.'ALM ') then

!LEANDRO INCLUDE EMPTY BANDS (START)

NEMAX_SAVE=NEMAX

NEMAX=NE

!LEANDRO INCLUDE EMPTY BANDS (END)

WRITE(24,2055) s_kvec,t_kvec,z_kvec,n,ne,bname

write(24,*) jatom,nemin,nemax,' jatom,nemin,nemax'

endif

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.

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 ...

I am of course taking into account the local orbitals in the process !

Are the wavefunctions only normalized below E_F ? If yes, do you know in

which part the normalization is done ?

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 ...

Thank you,

Leandro Salemi

_______________________________________________

Wien mailing list

Wien@zeus.theochem.tuwien.ac.at

<a rel="nofollow" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" style="color:rgb(160,30,30)">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a>

SEARCH the MAILING-LIST at:  

<a rel="nofollow" href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" style="color:rgb(160,30,30)">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a>

</pre>

</blockquote>

</div>

<br>

<p></p>

</div>

</body>

</html>