<div dir="ltr">Thanks for all of your replies<div><br></div><div style>they are really helpful!</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/4/4  <span dir="ltr"><<a href="mailto:wannier-request@quantum-espresso.org" target="_blank">wannier-request@quantum-espresso.org</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send Wannier mailing list submissions to<br>
        <a href="mailto:wannier@quantum-espresso.org">wannier@quantum-espresso.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://www.democritos.it/mailman/listinfo/wannier" target="_blank">http://www.democritos.it/mailman/listinfo/wannier</a><br>
or, via email, send a message with subject or body 'help' to<br>
        <a href="mailto:wannier-request@quantum-espresso.org">wannier-request@quantum-espresso.org</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:wannier-owner@quantum-espresso.org">wannier-owner@quantum-espresso.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of Wannier digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
   1. occupation of Wannier function (???)<br>
   2. Re: occupation of Wannier function (Ivo Souza)<br>
   3. Re: occupation of Wannier function (Arash Mostofi)<br>
   4. Re: occupation of Wannier function (Nicolas Lorente)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Wed, 3 Apr 2013 21:09:41 +0800<br>
From: ??? <<a href="mailto:creaman1986@gmail.com">creaman1986@gmail.com</a>><br>
Subject: [Wannier] occupation of Wannier function<br>
To: <a href="mailto:wannier@quantum-espresso.org">wannier@quantum-espresso.org</a><br>
Message-ID:<br>
        <<a href="mailto:CAE35B%2B8Y04xWB3S7_Q2WSKWnqDN-Lg1xEo8x7hRVG6gSd6ZX2A@mail.gmail.com">CAE35B+8Y04xWB3S7_Q2WSKWnqDN-Lg1xEo8x7hRVG6gSd6ZX2A@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Hi guys!<br>
<br>
I have a question about Wannier function.<br>
<br>
Can we get the occupation number of wannier function through this Wannier90<br>
code?<br>
<br>
thanks!<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://www.democritos.it/pipermail/wannier/attachments/20130403/e3476e2c/attachment.html" target="_blank">http://www.democritos.it/pipermail/wannier/attachments/20130403/e3476e2c/attachment.html</a><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Wed, 3 Apr 2013 17:49:42 +0000 (UTC)<br>
From: Ivo Souza <<a href="mailto:ivo_souza@ehu.es">ivo_souza@ehu.es</a>><br>
Subject: Re: [Wannier] occupation of Wannier function<br>
To: =?GB2312?B?6sXosOqQ?= <<a href="mailto:creaman1986@gmail.com">creaman1986@gmail.com</a>><br>
Cc: <a href="mailto:wannier@quantum-espresso.org">wannier@quantum-espresso.org</a><br>
Message-ID: <alpine.DEB.2.02.1304031737110.9167@levelho><br>
Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed<br>
<br>
> Can we get the occupation number of wannier function through this Wannier90 code?<br>
<br>
At present the code does not provide that information. I can think of a<br>
couple of ways in which this could be done:<br>
<br>
(i) Calculate the DOS projected onto the desired WF, and integrate up<br>
to the Fermi level. The calculation of the orbital-projected DOS will be<br>
included in the next release of wannier90.<br>
<br>
(ii) Calculate the BZ average of the occupation matrix in the Wannier<br>
gauge [see Eqs. 16 and 23 of PRB 85, 014435 (2012) for the definition].<br>
The orbital occupations are given by the diagonal elements.<br>
<br>
   Best,<br>
<br>
   Ivo<br>
<br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Wed, 03 Apr 2013 18:54:52 +0100<br>
From: Arash Mostofi <<a href="mailto:a.mostofi@imperial.ac.uk">a.mostofi@imperial.ac.uk</a>><br>
Subject: Re: [Wannier] occupation of Wannier function<br>
To: <a href="mailto:wannier@quantum-espresso.org">wannier@quantum-espresso.org</a><br>
Message-ID: <<a href="mailto:515C6CEC.4080407@imperial.ac.uk">515C6CEC.4080407@imperial.ac.uk</a>><br>
Content-Type: text/plain; charset=UTF-8; format=flowed<br>
<br>
Hello there,<br>
<br>
Adding to Ivo's reply:<br>
<br>
Without disentanglement the occupation numbers are always 1 (or 2 in the<br>
case of spin degeneracy). With disentanglement, and just considering the<br>
Gamma-point of the BZ for simplicity, the occupation "number" becomes an<br>
occupation matrix f_{ij}. See, for example, the discussion leading up to<br>
Eq (9) for f_{ij} in<br>
<br>
Andrinopoulos et al, J Chem Phys 135, 154105 (2011);<br>
<a href="http://link.aip.org/link/doi/10.1063/1.3647912" target="_blank">http://link.aip.org/link/doi/10.1063/1.3647912</a><br>
<br>
Computing f_{ij} is not implemented in v1.2 of Wannier90, but will be in<br>
the next release. It is fairly straightforward to code and would look<br>
something like this:<br>
<br>
! aam: calculate V = U_opt . U<br>
do s=1,num_wann<br>
  do k=1,ndim<br>
   do m=1,num_wann<br>
     v_matrix(k,s) = v_matrix(k,s) + u_matrix_opt(k,m,1)*u_matrix(m,s,1)<br>
   enddo<br>
  enddo<br>
enddo<br>
<br>
! aam: calculate f = V^dagger . V<br>
do r=1,num_wann<br>
  do s=1,num_wann<br>
   do k=1,ndim<br>
     f_w(r,s) = f_w(r,s) + v_matrix(k,s) * conjg(v_matrix(k,r))<br>
   enddo<br>
  enddo<br>
enddo<br>
<br>
Best wishes,<br>
<br>
Arash<br>
<br>
--<br>
Dr Arash A Mostofi<br>
Departments of Materials and Physics<br>
Deputy Director, CDT on Theory and Simulation of Materials<br>
Imperial College London, London SW7 2AZ, UK<br>
+44 (0)207 594 8154 | <a href="http://www.cmth.ph.ic.ac.uk/people/a.mostofi" target="_blank">www.cmth.ph.ic.ac.uk/people/a.mostofi</a><br>
<br>
On 03/04/2013 14:09, ??? wrote:<br>
> Hi guys!<br>
><br>
> I have a question about Wannier function.<br>
><br>
> Can we get the occupation number of wannier function through this<br>
> Wannier90 code?<br>
><br>
> thanks!<br>
><br>
><br>
><br>
><br>
> _______________________________________________<br>
> Wannier mailing list<br>
> <a href="mailto:Wannier@quantum-espresso.org">Wannier@quantum-espresso.org</a><br>
> <a href="http://www.democritos.it/mailman/listinfo/wannier" target="_blank">http://www.democritos.it/mailman/listinfo/wannier</a><br>
><br>
<br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Thu, 04 Apr 2013 08:15:19 +0200<br>
From: Nicolas Lorente <<a href="mailto:nicolas.lorente@cin2.es">nicolas.lorente@cin2.es</a>><br>
Subject: Re: [Wannier] occupation of Wannier function<br>
To: <a href="mailto:wannier@quantum-espresso.org">wannier@quantum-espresso.org</a><br>
Message-ID: <<a href="mailto:515D1A77.7050401@cin2.es">515D1A77.7050401@cin2.es</a>><br>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
<br>
On 03/04/13 19:49, Ivo Souza wrote:<br>
>> Can we get the occupation number of wannier function through this Wannier90 code?<br>
> At present the code does not provide that information. I can think of a<br>
> couple of ways in which this could be done:<br>
><br>
> (i) Calculate the DOS projected onto the desired WF, and integrate up<br>
> to the Fermi level. The calculation of the orbital-projected DOS will be<br>
> included in the next release of wannier90.<br>
><br>
> (ii) Calculate the BZ average of the occupation matrix in the Wannier<br>
> gauge [see Eqs. 16 and 23 of PRB 85, 014435 (2012) for the definition].<br>
> The orbital occupations are given by the diagonal elements.<br>
><br>
>     Best,<br>
><br>
>     Ivo<br>
> _______________________________________________<br>
> Wannier mailing list<br>
> <a href="mailto:Wannier@quantum-espresso.org">Wannier@quantum-espresso.org</a><br>
> <a href="http://www.democritos.it/mailman/listinfo/wannier" target="_blank">http://www.democritos.it/mailman/listinfo/wannier</a><br>
<br>
To follow on idea (i), I sent our wpdos.f90 code to compute the PDOS<br>
using wannier90 in this thread:<br>
<br>
<a href="http://www.democritos.it/pipermail/wannier/2012-November/000603.html" target="_blank">http://www.democritos.it/pipermail/wannier/2012-November/000603.html</a><br>
<br>
Best,<br>
<br>
Nicol?s<br>
<br>
<br>
<br>
--<br>
Nicol?s Lorente<br>
CIN2 (CSIC)<br>
CIN2 building, UAB<br>
08193 Bellaterra, Spain<br>
Tel. <a href="tel:%2B34-937373622" value="+34937373622">+34-937373622</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
Wannier mailing list<br>
<a href="mailto:Wannier@quantum-espresso.org">Wannier@quantum-espresso.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/wannier" target="_blank">http://www.democritos.it/mailman/listinfo/wannier</a><br>
<br>
<br>
End of Wannier Digest, Vol 63, Issue 1<br>
**************************************<br>
</blockquote></div><br></div>