Dear Jonathan,<br><br>thank you for your reply. <br><br>I think the symmetry of the WFs should be of p character.<br><br>I have however one problem.<br>I have 4 bands.<br>The number of Carbon ions in one molecule is 22 (and there are 2 molecules).<br>
<br>How can i set the flags:<br><br>begin projections<br> <br> end projections<br><br>in such case in order to avoid the code messagge: "param_get_projection: too many projections defined"?<br><br>ciao<br><br>G.<br>
<br><br><div class="gmail_quote">On Tue, Apr 13, 2010 at 11:48 PM, Jonathan Yates <span dir="ltr"><<a href="mailto:jonathan.yates@materials.ox.ac.uk">jonathan.yates@materials.ox.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div style=""><br><div><div class="im"><div>On 13 Apr 2010, at 14:08, Gianluca Giovannetti wrote:</div><br><blockquote type="cite">Dear All,<br><br>i`m calculating Wannier functions for a molecular crystal (sites are molecule with C and H ions within each molecule).<br>
<br>i want to plot Wannier functions of the 4 bands at Fermi level.<br>the system is metallic.<br> these four bands come from HOMO and LUMO molecular states of the two inequivalent molecules in the unit cell.<br></blockquote>
<div><br></div></div><div>Gianluca,</div><div><br></div><div> I'm assuming that these four bands which span the fermi level form an isolated group (so no disentanglement is needed)</div><div class="im"><div><br></div>
<br><blockquote type="cite">My feeling is that the setting in the win file:<br>begin projections<br> random<br> end projections<br>is not really good as starting point.<br>Indeed the CONV is 1.0E-03.<br>HOMO and LUMO are localized on the molecules.<br>
What is the best strategy in this case?<br></blockquote><div><br></div></div><div>I think using randomly placed Gaussian functions is likely to be particularly bad in this case. If the orbitals are already fairly localised, and the cell is pretty large, then the trial functions might have little overlap with the bloch states. That won't make the calculation wrong - just that the starting point could be poor.</div>
<div> So if you think that the states are 2p like on Oxygen then try such projections as a starting guess. Sometimes a two-step procedure might help: look at the WF resulting from the random projections - if they appear to have some symmetry character then try again using a projection of that form.</div>
<div><br></div><div> It might help in this procedure to use a coarser kpoint mesh (4x4x4 or even 2x2x2) - I find that the minimiser can cope better with a poor starting guess in this case. Then when I know what the MLWF should look like, I increase the mesh.</div>
<div class="im"><div><br></div><div><blockquote type="cite"><font color="#000000">However i get Wannier eigenvalues which are not exactly the same as the one i have calculated Quantum espresso (the discrepancy is not big but finite at some k-points).</font></blockquote>
<br></div></div><div> You're obtaining the band-structure by interpolation. At the kpoints which lie in the original (ab-initio) grid the energies will be reproduced exactly. At any other points it is an interpolation. Importantly the quality of the interpolation increases exponentially as the size of kpoint mesh is increased.</div>
<div> However, it is probably a good plan to get well localised wannier functions first - as this may improve the interpolation.</div><div><br></div><div> Jonathan</div><div><br></div><div> </div></div><br><div><div><div>
<div><div><div><br>-- <br>Department of Materials, University of Oxford, Parks Road, Oxford, OX1 3PH, UK<br>tel: +44 (0)1865 612797                <a href="http://users.ox.ac.uk/%7Eoums0549/" target="_blank">http://users.ox.ac.uk/~oums0549/</a><br>
</div></div></div></div></div> </div><br></div><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></blockquote></div><br>