<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br><div><div>On Jan 15, 2008, at 9:00 AM, Min Seung Kyu wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">Dear users and developers,<br><br>I have a question about eigenfunctions in PWSCF codes.<br><br>I looked at PWSCF source files and I found that the 'evc' variable is <br>KS eigenvectors in G-space, right?<br><br>Because Hamiltonian matrix H(G,G') is hermitian, 'evc' should be real.</blockquote><div><br class="webkit-block-placeholder"></div>Who said this?&nbsp;</div><div><br class="webkit-block-placeholder"></div><div>The eigenvalues of a Hermitean matrix are real, the eigenvectors not necessarily so. If in addition to being Hermitean the matrix is also real (which implies it is symmetric), the the eigenvectors CAN BE CHOSEN to be real. If you feed a real matrix to a complex eigenvalue solver, then the eigenvectors will most probably come out complex, but can they be easily converted to real (if needed) by an appropriate choice of phases (again, if the original matrix is real symmetric, not merely Hermitean)</div><div><br class="webkit-block-placeholder"></div><div>Hope this helps</div><div>Stefano</div><div><br><blockquote type="cite"><br><br>When I print out 'evc', however, 'evc' contains imaginary parts.<br><br>For example, when I run 'co.rx.in' (example03 folder), I got these values:<br><br> &nbsp;DBLE(evc(1:4,1)) &nbsp;&nbsp;&nbsp;&nbsp;0.186884 &nbsp;&nbsp;&nbsp;&nbsp;0.154551 &nbsp;&nbsp;&nbsp;&nbsp;0.154551 &nbsp;&nbsp;&nbsp;&nbsp;0.142085<br>AIMAG(evc(1:4,1)) &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;-0.048560<br> &nbsp;DBLE(evc(1:4,2)) &nbsp;&nbsp;&nbsp;-0.024799 &nbsp;&nbsp;&nbsp;-0.021751 &nbsp;&nbsp;&nbsp;-0.021751 &nbsp;&nbsp;&nbsp;-0.068645<br>AIMAG(evc(1:4,2)) &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;-0.145963<br> &nbsp;DBLE(evc(1:4,3)) &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;&nbsp;0.000000<br>AIMAG(evc(1:4,3)) &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;-0.005668 &nbsp;&nbsp;&nbsp;&nbsp;0.152298 &nbsp;&nbsp;&nbsp;&nbsp;0.000000<br> &nbsp;DBLE(evc(1:4,4)) &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;&nbsp;0.000000<br>AIMAG(evc(1:4,4)) &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;-0.152298 &nbsp;&nbsp;&nbsp;-0.005668 &nbsp;&nbsp;&nbsp;&nbsp;0.000000<br> &nbsp;DBLE(evc(1:4,5)) &nbsp;&nbsp;&nbsp;-0.299154 &nbsp;&nbsp;&nbsp;-0.185970 &nbsp;&nbsp;&nbsp;-0.185970 &nbsp;&nbsp;&nbsp;&nbsp;0.011252<br>AIMAG(evc(1:4,5)) &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;&nbsp;0.000000 &nbsp;&nbsp;&nbsp;&nbsp;0.168998<br><br>( I added 'WRITE(*,*) evc(1:4,ibnd)' in 'sum_band.f90'. )<br><br>What is the problem?<br><br>Thank you in advance.<br><br>Sincerely, Min.<br><br>-------------------------------Sig.---------------------------<br>Center for Superfunctional Materials, Department of Chemistry,<br>Pohang University of Science and Technology (POSTECH)<br><br>Mr. Seung Kyu, Min<br><br>Ph.D. Candidate<br><br>E-Mail: <a href="mailto:min0220@postech.ac.kr">min0220@postech.ac.kr</a><br>Tel. &nbsp;: 82-54-279-5858<br>----------------------------------------------<br><br><br><br><br><br>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>http://www.democritos.it/mailman/listinfo/pw_forum<br></blockquote></div><br><div> <span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div style="word-wrap: break-word; -khtml-nbsp-mode: space; -khtml-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px;"><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; ">---</span></span></span></span></span></span></span></span></span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px;"><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; ">Stefano Baroni - SISSA</span></span></span></span></span></span></span></span></span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px;"><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; ">&nbsp; </span></span></span></span></span></span></span></span></span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px;"><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; ">&amp;</span></span></span></span></span></span></span></span></span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px;"><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; ">&nbsp; </span></span></span></span></span></span></span></span></span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px;"><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; ">DEMOCRITOS National Simulation Center - Trieste</span></span></span></span></span></span></span></span></span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px;"><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; ">[+39] 040 3787 406 (tel) -528 (fax) / stefanobaroni (skype)</span></span></span></span></span></span></span></span></span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px;"><br class="khtml-block-placeholder"></span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px;"><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; ">Please, if possible, don't</span></span></span></span></span></span></span></span></span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px;"><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; ">&nbsp; </span></span></span></span></span></span></span></span></span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px;"><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; ">send me MS Word or PowerPoint attachments</span></span></span></span></span></span></span></span></span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px;"><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; ">Why? See:</span></span></span></span></span></span></span></span></span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px;"><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; "><span class="Apple-style-span" style="font-size: 12px; ">&nbsp; </span></span></span></span></span></span></span></span></span></font><a href="http://www.gnu.org/philosophy/no-word-attachments.html"><font class="Apple-style-span" color="#0023E9" size="3"><span class="Apple-style-span" style="font-size: 12px;; color: rgb(0, 35, 233); -khtml-text-decorations-in-effect: underline; "><span class="Apple-style-span" style="color: rgb(0, 35, 233); font-size: 12px; -khtml-text-decorations-in-effect: underline; "><span class="Apple-style-span" style="color: rgb(0, 35, 233); font-size: 12px; -khtml-text-decorations-in-effect: underline; "><span class="Apple-style-span" style="color: rgb(0, 35, 233); font-size: 12px; -khtml-text-decorations-in-effect: underline; "><span class="Apple-style-span" style="color: rgb(0, 35, 233); font-size: 12px; -khtml-text-decorations-in-effect: underline; "><span class="Apple-style-span" style="color: rgb(0, 35, 233); font-size: 12px; -khtml-text-decorations-in-effect: underline; "><span class="Apple-style-span" style="color: rgb(0, 35, 233); font-size: 12px; -khtml-text-decorations-in-effect: underline; "><span class="Apple-style-span" style="color: rgb(0, 35, 233); font-size: 12px; -khtml-text-decorations-in-effect: underline; "><span class="Apple-style-span" style="color: rgb(0, 35, 233); font-size: 12px; -khtml-text-decorations-in-effect: underline; ">http://www.gnu.org/philosophy/no-word-attachments.html</span></span></span></span></span></span></span></span></span></font></a></div><div><br class="khtml-block-placeholder"></div><br class="Apple-interchange-newline"></span></span></span></span></span></span></span></span></div></span> </div><br></body></html>