<div dir="ltr"><div><div>Dear Henry,<br></div> In the PAW method, the atomic pseudo basis functions are not designed to be orthogonal to each other, but there is a "dual relationship" <p_i|\tilde{\psi}_j>=\delta_{ij} with the projector function p_i. One can use the projector function to estimate the charge associated with an atomic site within the augmentation sphere about that site. (Some notes are enclosed about this in case you are interested. The main equations are correct, but the programming details are no longer true.) Paolo Giannozzi programmed this into quantum espresso in version 5.3.0. Perhaps this might suite your needs?? In order to use it you need to use &projwfc:<br><br> &projwfc<br> pawproj= .true. , << new option<br> outdir='$outd/', << old options; change as appropriate<br> prefix='$label',<br> filpdos='pdos',<br> Emin=-25.0, Emax=25.0, DeltaE=0.01,<br> ngauss=0, degauss=0.01<br>/<br><br></div>Sincerely, Natalie Holzwarth<br><div><br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature">N. A. W. Holzwarth email: <a href="mailto:natalie@wfu.edu" target="_blank">natalie@wfu.edu</a><div>Department of Physics web: <a href="http://www.wfu.edu/~natalie" target="_blank">http://www.wfu.edu/~natalie</a></div><div>Wake Forest University phone: 1-336-758-5510 </div><div>Winston-Salem, NC 27109 USA office: Rm. 300 Olin Physical Lab</div></div></div>
<br><div class="gmail_quote">On Fri, Apr 15, 2016 at 8:02 PM, Henry J Seeley <span dir="ltr"><<a href="mailto:hseeley@uoregon.edu" target="_blank">hseeley@uoregon.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello all,<br>
<br>
I'm currently trying to use Quantum Espresso to generate the eigenfunctions for bulk and slab systems of PbS, which I then plan on projecting onto one another to determine the relative surface or bulk character of specific slab states. I've generated my sample systems and have run the scf/nscf calculations, but I've hit some trouble with the eigenvectors obtained from 'pw_export.x'. I have a couple questions that I hope some of you may be able to answer. Of course all my input files will be attached.<br>
<br>
To my surprise I found that different eigenstates of the same system are NOT orthonormal (<psi_i|psi_i> ~= 0.75, abs(<psi_i|psi_j>) ~= 0.15). I've done some searching and determined this may have to do with the PAW psuedo-potentials I am using. Is the PAW method the problem here?<br>
<br>
In trying to troubleshoot this, I found that eigenstate orthonormality may be obtained by including the 'uspp_spsi = .TRUE.' option in the 'pw_export.x' input file, which produces S|psi>, which is orthonormal by <psi_i|S|psi_j>. Is this correct? If so, how can I project the two different systems on one another (bulk/slab), i.e. which system's eigenvector gets the 'S'?<br>
<br>
Finally (and sorry for all the questions at once), I've tried using the 'uspp_spsi = .TRUE.' option for both the bulk and slab systems, but I receive an error message for the bulk (but not the slab!):<br>
<br>
*** glibc detected *** pw_export.x: malloc(): memory corruption: 0x0000000002b3a650 ***<br>
<br>
This occurs right as pw_export.x tries to produce the S|psi> files.<br>
I'm at a loss, and cannot proceed to troubleshoot this by myself any further. Does anyone have any insight into:<br>
(A) Whether I'm even taking the right approach, or (B) Why I could be getting this error for the bulk system but not the slab?<br>
<br>
<br>
Thank you all for your time,<br>
Hank Seeley<br>
University of Oregon, Chemistry Dept.<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></div>