[Wannier] Re: Too many projections defined error

Brad Malone brad.malone@gmail.com
Sat, 26 May 2007 11:57:30 -0700


------=_Part_147304_20622608.1180205850595
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit
Content-Disposition: inline

Thanks for the response Arash.  Setting the projections to random (they were
set to sp3) did indeed fix that issue.  I ran the wannier code and generated
the .nnkp file.  However, in the nnkp file it lists my excluded bands as
being 25-32 (which I had specified in the .win input file; perhaps
unnecessary because my window excludes them anyway) but also it lists band
#8 as being excluded, and I did not want or specify that.  Is there a reason
for this?  I'm worried that it may be causing me an error when I run the
Wannier code again because I get an error of "param_read: mismatch in
silicon.eig".  Inspecting the silicon.eig file it looks like it should be
okay, unless the code thinks that there are 23 bands (having excluded band
8), as there are 24 bands (for my 64 points) listed in the silicon.eig file.

I'm also interested in the bottom of the conduction band.  Do you have a
feel for how robust this is to the initial guess for the projections?  And
is it possible for me to specify sp3 projections without having 48 wannier
bands?  I'd like to use fewer bands if possible because it seems that I'm
running into memory allocation issues that prevent me from starting with a
denser k-grid. Thanks again for the help.

Brad



> Message: 1
> Date: Fri, 25 May 2007 12:58:31 -0700
> From: "Brad Malone" <brad.malone@gmail.com>
> To: wannier@quantum-espresso.org
> Subject: [Wannier] Too many projections defined error
>
> ------=_Part_139432_23760487.1180123111998
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
> Content-Transfer-Encoding: 7bit
> Content-Disposition: inline
>
> Hi, I'm having an issue that I'm very confused about.  I'm running a
> calculation on silicon with 12 atoms per unit cell.  There is a gap in my
> system (pwSCF tells me that 7.3 is the highest occupied state and 8.4 is
> the
> lowest unoccupied state).  I just care about the valence bands and so I
> set
> my upper energy window to 7.8 eV and set num_wann=24 (num_bands=32).  My
> frozen energy window is the same as my larger window and it covers all 24
> valence bands.  However, when I run the wannier code to generate the .nnkp
> file it tells me the following error:
>
> param_get_projection: too many projections defined
>
>
>
> Any help in understanding why this is would be greatly appreciated!
>
> Best,
> Brad
>
> ------=_Part_139432_23760487.1180123111998
> Content-Type: text/html; charset=ISO-8859-1
> Content-Transfer-Encoding: 7bit
> Content-Disposition: inline
>
> Hi, I&#39;m having an issue that I&#39;m very confused about.&nbsp;
> I&#39;m running a calculation on silicon with 12 atoms per unit cell.&nbsp;
> There is a gap in my system (pwSCF tells me that 7.3 is the highest
> occupied state and 8.4
> is the lowest unoccupied state).&nbsp; I just care about the valence bands
> and so I set my upper energy window to 7.8 eV and set num_wann=24
> (num_bands=32).&nbsp; My frozen energy window is the same as my larger
> window and it covers all 24 valence bands.&nbsp; However, when I run the
> wannier code to generate the .nnkp file it tells me the following error:
> <br><br><blockquote style="border-left: 1px solid rgb(204, 204, 204);
> margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"
> class="gmail_quote">&nbsp;param_get_projection: too many projections
> defined</blockquote><div><br><br>Any help in understanding why this is would
> be greatly appreciated!
> <br><br>Best,<br>Brad<br></div><br><br>
>
> ------=_Part_139432_23760487.1180123111998--
>
> --__--__--
>
> Message: 2
> Date: Fri, 25 May 2007 16:22:28 -0400 (EDT)
> From: Arash A Mostofi <mostofi@MIT.EDU>
> To: Wannier List <wannier@quantum-espresso.org>
> Subject: Re: [Wannier] Too many projections defined error
>
> Dear Brad
>
> I think the problem is in the "projections" block of the input file, not
> the energy windows. As you say, a 12 atom unit cell of silicon will need
> 24
> maximally localised Wannier functions (MLWFs) to span the valence manifold
> and this is the number of projections you should define in the projections
> block.
>
> What we've found is that the generation of MLWFs for the valence manifold
> is very robust to the initial guess, so even specifying randomly centred
> s-orbital projections is good enough. This can be done by setting
>
> Begin Projections
> random
> End Projections
>
> in the input file.
>
> Cheers
>
> Arash
>
> -----------------------------------------
> Dr Arash A Mostofi
> Dept of Materials Science and Engineering
> 13-4084 MIT
> 77 Massachusetts Avenue
> Cambridge MA 02139-4307 USA
> Tel: 617.2536026
> -----------------------------------------
>
>
>

------=_Part_147304_20622608.1180205850595
Content-Type: text/html; charset=ISO-8859-1
Content-Transfer-Encoding: 7bit
Content-Disposition: inline

Thanks for the response Arash.&nbsp; Setting the projections to random (they were set to sp3) did indeed fix that issue.&nbsp; I ran the wannier code and generated the .nnkp file.&nbsp; However, in the nnkp file it lists my excluded bands as being 25-32 (which I had specified in the .win input file; perhaps unnecessary because my window excludes them anyway) but also it lists band #8 as being excluded, and I did not want or specify that.&nbsp; Is there a reason for this?&nbsp; I&#39;m worried that it may be causing me an error when I run the Wannier code again because I get an error of &quot;param_read: mismatch in 
silicon.eig&quot;.&nbsp; Inspecting the silicon.eig file it looks like it should be okay, unless the code thinks that there are 23 bands (having excluded band 8), as there are 24 bands (for my 64 points) listed in the silicon.eig
 file.<br><br>I&#39;m also interested in the bottom of the conduction band.&nbsp; Do you have a feel for how robust this is to the initial guess for the projections?&nbsp; And is it possible for me to specify sp3 projections without having 48 wannier bands?&nbsp; I&#39;d like to use fewer bands if possible because it seems that I&#39;m running into memory allocation issues that prevent me from starting with a denser k-grid. Thanks again for the help.
<br><br>Brad<br><br><div><span class="gmail_quote"><br></span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><br>Message: 1<br>Date: Fri, 25 May 2007 12:58:31 -0700
<br>From: &quot;Brad Malone&quot; &lt;<a href="mailto:brad.malone@gmail.com">brad.malone@gmail.com</a>&gt;<br>To: <a href="mailto:wannier@quantum-espresso.org">wannier@quantum-espresso.org</a><br>Subject: [Wannier] Too many projections defined error
<br><br>------=_Part_139432_23760487.1180123111998<br>Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>Content-Transfer-Encoding: 7bit<br>Content-Disposition: inline<br><br>Hi, I&#39;m having an issue that I&#39;m very confused about.&nbsp;&nbsp;I&#39;m running a
<br>calculation on silicon with 12 atoms per unit cell.&nbsp;&nbsp;There is a gap in my<br>system (pwSCF tells me that 7.3 is the highest occupied state and 8.4 is the<br>lowest unoccupied state).&nbsp;&nbsp;I just care about the valence bands and so I set
<br>my upper energy window to 7.8 eV and set num_wann=24 (num_bands=32).&nbsp;&nbsp;My<br>frozen energy window is the same as my larger window and it covers all 24<br>valence bands.&nbsp;&nbsp;However, when I run the wannier code to generate the .nnkp
<br>file it tells me the following error:<br><br> param_get_projection: too many projections defined<br><br><br><br>Any help in understanding why this is would be greatly appreciated!<br><br>Best,<br>Brad<br><br>------=_Part_139432_23760487.1180123111998
<br>Content-Type: text/html; charset=ISO-8859-1<br>Content-Transfer-Encoding: 7bit<br>Content-Disposition: inline<br><br>Hi, I&amp;#39;m having an issue that I&amp;#39;m very confused about.&amp;nbsp; I&amp;#39;m running a calculation on silicon with 12 atoms per unit cell.&amp;nbsp; There is a gap in my system (pwSCF tells me that 
7.3 is the highest occupied state and 8.4<br> is the lowest unoccupied state).&amp;nbsp; I just care about the valence bands and so I set my upper energy window to 7.8 eV and set num_wann=24 (num_bands=32).&amp;nbsp; My frozen energy window is the same as my larger window and it covers all 24 valence bands.&amp;nbsp; However, when I run the wannier code to generate the .nnkp file it tells me the following error:
<br>&lt;br&gt;&lt;br&gt;&lt;blockquote style=&quot;border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;&quot; class=&quot;gmail_quote&quot;&gt;&amp;nbsp;param_get_projection: too many projections defined&lt;/blockquote&gt;&lt;div&gt;&lt;br&gt;&lt;br&gt;Any help in understanding why this is would be greatly appreciated!
<br>&lt;br&gt;&lt;br&gt;Best,&lt;br&gt;Brad&lt;br&gt;&lt;/div&gt;&lt;br&gt;&lt;br&gt;<br><br>------=_Part_139432_23760487.1180123111998--<br><br>--__--__--<br><br>Message: 2<br>Date: Fri, 25 May 2007 16:22:28 -0400 (EDT)<br>
From: Arash A Mostofi &lt;<a href="mailto:mostofi@MIT.EDU">mostofi@MIT.EDU</a>&gt;<br>To: Wannier List &lt;<a href="mailto:wannier@quantum-espresso.org">wannier@quantum-espresso.org</a>&gt;<br>Subject: Re: [Wannier] Too many projections defined error
<br><br>Dear Brad<br><br>I think the problem is in the &quot;projections&quot; block of the input file, not<br>the energy windows. As you say, a 12 atom unit cell of silicon will need 24<br>maximally localised Wannier functions (MLWFs) to span the valence manifold
<br>and this is the number of projections you should define in the projections<br>block.<br><br>What we&#39;ve found is that the generation of MLWFs for the valence manifold<br>is very robust to the initial guess, so even specifying randomly centred
<br>s-orbital projections is good enough. This can be done by setting<br><br>Begin Projections<br>random<br>End Projections<br><br>in the input file.<br><br>Cheers<br><br>Arash<br><br>-----------------------------------------
<br>Dr Arash A Mostofi<br>Dept of Materials Science and Engineering<br>13-4084 MIT<br>77 Massachusetts Avenue<br>Cambridge MA 02139-4307 USA<br>Tel: 617.2536026<br>-----------------------------------------<br><br><br></blockquote>
</div><br>

------=_Part_147304_20622608.1180205850595--