<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
Dear Dr Giuseppe Mattioli:<BR>
Many thanks!<BR>
&nbsp;<BR>
You mean I should adopt the&nbsp;&nbsp;e = -5.55656 eV, because the # 408 wavefunction is the&nbsp;largest part of&nbsp;wavefunction psi?<BR>
&nbsp;<BR>
The eigenvalue is slightly different in different k points, should I average them and get the mean value as the He 1S eigenvalue?<BR>
&nbsp;<BR>
<BR>My pdos.out is as follows:<BR>
<SPAN lang=EN>......<BR>
state # 397: atom 62 (Cd ), wfc 3 (l=1 m= 2)<BR>
state # 398: atom 62 (Cd ), wfc 3 (l=1 m= 3)<BR>
state # 399: atom 63 (Cd ), wfc 1 (l=2 m= 1)<BR>
state # 400: atom 63 (Cd ), wfc 1 (l=2 m= 2)<BR>
state # 401: atom 63 (Cd ), wfc 1 (l=2 m= 3)<BR>
state # 402: atom 63 (Cd ), wfc 1 (l=2 m= 4)<BR>
state # 403: atom 63 (Cd ), wfc 1 (l=2 m= 5)<BR>
state # 404: atom 63 (Cd ), wfc 2 (l=0 m= 1)<BR>
state # 405: atom 63 (Cd ), wfc 3 (l=1 m= 1)<BR>
state # 406: atom 63 (Cd ), wfc 3 (l=1 m= 2)<BR>
state # 407: atom 63 (Cd ), wfc 3 (l=1 m= 3)<BR>
state # 408: atom 64 (He ), wfc 1 (l=0 m= 1)<BR>
.....<BR><SPAN lang=EN>
k = 0.0000000000 0.0000000000 0.0000000000<BR>
.....<BR></SPAN></SPAN><SPAN lang=EN>
e = -5.55656 eV<BR>
psi = <FONT color=#ff0000>0.696*[# 408]+</FONT>0.052*[# 105]+0.051*[# 5]+0.050*[# 77]+0.040*[# 33]+<BR>
+0.003*[# 203]+0.003*[# 223]+0.004*[# 229]+0.003*[# 157]+0.003*[# 239]+<BR>
+0.003*[# 358]+0.003*[# 6]+0.003*[# 108]+0.003*[# 7]+0.003*[# 8]+<BR></SPAN>
......<BR><SPAN lang=EN>
k = 0.0000000000 0.0000000000 -0.5000000000<BR></SPAN>
......<BR><SPAN lang=EN>
Lowdin Charges:<BR></SPAN>
......<BR>
&nbsp;<BR>
&nbsp;<BR>
&nbsp;<BR>
Thanks !<BR>
&nbsp;<BR>
&nbsp;<BR>
Rgds<BR>Jianchun Wu<BR>Department of physics<BR>Soochow university China.215006<BR><BR>
&nbsp;<BR>
&nbsp;<BR>
&gt; From: giuseppe.mattioli@mlib.ism.cnr.it<BR>&gt; To: pw_forum@pwscf.org<BR>&gt; Date: Mon, 3 May 2010 13:12:17 +0200<BR>&gt; Subject: Re: [Pw_forum] the potential of core electron<BR>&gt; <BR>&gt; &gt; Shall I do the NSCF between the scf and projwfc.x as the DOS calculation?<BR>&gt; <BR>&gt; No, but if you need it for different purposes.<BR>&gt; <BR>&gt; &gt; Could you tell me where I can find He 1s eigenvalue ?<BR>&gt; <BR>&gt; state #651: atom 192 (H ), wfc 1 (l=0 m= 1)<BR>&gt; state #652: atom 193 (H ), wfc 1 (l=0 m= 1)<BR>&gt; state #653: atom 194 (H ), wfc 1 (l=0 m= 1)<BR>&gt; state #654: atom 195 (H ), wfc 1 (l=0 m= 1)<BR>&gt; state #655: atom 196 (H ), wfc 1 (l=0 m= 1)<BR>&gt; state #656: atom 197 (H ), wfc 1 (l=0 m= 1)<BR>&gt; state #657: atom 198 (H ), wfc 1 (l=0 m= 1)<BR>&gt; state #658: atom 199 (He ), wfc 1 (l=0 m= 1)<BR>&gt; <BR>&gt; look for the # XXX wavefunction corresponding to your He atom at the beginning <BR>&gt; of the pdos.out file, and match it with a given eigenvalue in the bulk of the <BR>&gt; same file.<BR>&gt; <BR>&gt; Giuseppe<BR>&gt; <BR>&gt; <BR>&gt; <BR>&gt; On Friday 30 April 2010 12:19:15 wujianchun wrote:<BR>&gt; &gt; Dear Dr Giuseppe Mattioli:<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt; Thank you again!<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt; Shall I do the NSCF between the scf and projwfc.x as the DOS calculation?<BR>&gt; &gt;<BR>&gt; &gt; I think I should get the PDOS of He.<BR>&gt; &gt;<BR>&gt; &gt; Could you tell me where I can find He 1s eigenvalue ?<BR>&gt; &gt;<BR>&gt; &gt; There are some "e" values after each k points in the outfile:pdos.out<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt; "k = 0.0000000000 0.0000000000 0.0000000000<BR>&gt; &gt;<BR>&gt; &gt; e = -7.05254 eV<BR>&gt; &gt;<BR>&gt; &gt; psi = 0.022*[# 1]+0.022*[# 5]+0.022*[# 9]+0.022*[# 13]+0.022*[# 17]+<BR>&gt; &gt;<BR>&gt; &gt; +0.022*[# 21]+0.022*[# 25]+0.022*[# 29]+0.022*[# 33]+0.022*[# 37]+<BR>&gt; &gt;<BR>&gt; &gt; +0.022*[# 41]+0.022*[# 45]+0.022*[# 49]+0.022*[# 53]+0.022*[# 57]+<BR>&gt; &gt;<BR>&gt; &gt; +0.022*[# 61]+0.022*[# 65]+0.022*[# 69]+0.022*[# 73]+0.022*[# 77]+<BR>&gt; &gt;<BR>&gt; &gt; +0.022*[# 81]+0.022*[# 85]+0.022*[# 89]+0.022*[# 93]+0.022*[# 97]+<BR>&gt; &gt;<BR>&gt; &gt; +0.022*[# 101]+0.022*[# 105]+0.022*[# 109]+0.022*[# 113]+0.022*[# 117]+<BR>&gt; &gt;<BR>&gt; &gt; +0.022*[# 121]+0.022*[# 125]+0.009*[# 134]+0.009*[# 143]+0.009*[# 152]+<BR>&gt; &gt;<BR>&gt; &gt; +0.009*[# 161]+0.009*[# 170]+0.009*[# 179]+0.009*[# 188]+0.009*[# 197]+<BR>&gt; &gt;<BR>&gt; &gt; +0.009*[# 206]+0.009*[# 215]+0.009*[# 224]+0.009*[# 233]+0.009*[# 242]+<BR>&gt; &gt;<BR>&gt; &gt; +0.009*[# 251]+0.009*[# 260]+0.009*[# 269]+0.009*[# 278]+0.009*[# 287]+<BR>&gt; &gt;<BR>&gt; &gt; +0.009*[# 296]+0.009*[# 305]+0.009*[# 314]+0.009*[# 323]+0.009*[# 332]+<BR>&gt; &gt;<BR>&gt; &gt; +0.009*[# 341]+0.009*[# 350]+0.009*[# 359]+0.009*[# 368]+0.009*[# 377]+<BR>&gt; &gt;<BR>&gt; &gt; +0.009*[# 386]+0.009*[# 395]+0.009*[# 404]+0.009*[# 413]+<BR>&gt; &gt;<BR>&gt; &gt; |psi|^2 = 0.998<BR>&gt; &gt;<BR>&gt; &gt; e = -6.79505 eV"<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt; I guess they are the eigenvalues, and the lowest is He 1s eigenvalue in our<BR>&gt; &gt; CdS system,<BR>&gt; &gt;<BR>&gt; &gt; am I right?<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt; RgdS<BR>&gt; &gt;<BR>&gt; &gt; Jianchun<BR>&gt; &gt; Department of physics<BR>&gt; &gt; Soochow university China.215006<BR>&gt; &gt;<BR>&gt; &gt; &gt; From: giuseppe.mattioli@mlib.ism.cnr.it<BR>&gt; &gt; &gt; To: pw_forum@pwscf.org<BR>&gt; &gt; &gt; Date: Fri, 30 Apr 2010 12:43:04 +0200<BR>&gt; &gt; &gt; Subject: Re: [Pw_forum] the potential of core electron<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; And, sorry, make a projwfc.x calculation after the scf one to find the He<BR>&gt; &gt; &gt; 1s eigenvalue if it does not stand alone by comparing scf and scf+He<BR>&gt; &gt; &gt; calculations<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; Giuseppe<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; On Friday 30 April 2010 03:17:41 wujianchun wrote:<BR>&gt; &gt; &gt; &gt; Dear Dr Giuseppe Mattioli:<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; Could you tell me the details on using the 1s level of a He atom to<BR>&gt; &gt; &gt; &gt; align the potential.<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; My system is 64 atoms CdS supercell. Need we replace the Cd with He at<BR>&gt; &gt; &gt; &gt; the corner?How to get the 1S level of He in the output file?<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; Thanks a lot!<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; Regdards<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; Jianchun<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; Department of physics<BR>&gt; &gt; &gt; &gt; Soochow university China.215006<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; From: jcwu.suda@hotmail.com<BR>&gt; &gt; &gt; &gt; To: pw_forum@pwscf.org<BR>&gt; &gt; &gt; &gt; Date: Thu, 29 Apr 2010 18:57:26 +0800<BR>&gt; &gt; &gt; &gt; Subject: Re: [Pw_forum] the potential of core electron<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; Dear Dr Giuseppe Mattioli:<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; Many thanks!<BR>&gt; &gt; &gt; &gt; I'll try the method.<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; Rgds<BR>&gt; &gt; &gt; &gt; Jianchun<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; To: pw_forum@pwscf.org<BR>&gt; &gt; &gt; &gt; &gt; From: giuseppe.mattioli@mlib.ism.cnr.it<BR>&gt; &gt; &gt; &gt; &gt; Date: Thu, 29 Apr 2010 13:17:38 +0200<BR>&gt; &gt; &gt; &gt; &gt; Subject: Re: [Pw_forum] the potential of core electron<BR>&gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; Dear Jianchun<BR>&gt; &gt; &gt; &gt; &gt; You cannot get the core levels in pseudopotential calculations. At<BR>&gt; &gt; &gt; &gt; &gt; least, you cannot get such levels in a straightforward way. And you<BR>&gt; &gt; &gt; &gt; &gt; don't need! If your host matrix is not close packed, try to use the<BR>&gt; &gt; &gt; &gt; &gt; 1s level of a He atom to align the potential. The results can be<BR>&gt; &gt; &gt; &gt; &gt; often quite good...<BR>&gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; Giuseppe<BR>&gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; On Wednesday 28 April 2010 03:29:25 wujianchun wrote:<BR>&gt; &gt; &gt; &gt; &gt; &gt; Dear All,<BR>&gt; &gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; &gt; In order to align the potential between the defect supercell and<BR>&gt; &gt; &gt; &gt; &gt; &gt; perfect supercell.<BR>&gt; &gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; &gt; I want to know the potential of core electrons(core-level).<BR>&gt; &gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; &gt; Could you tell me whether the PWSCF can provide the value of the<BR>&gt; &gt; &gt; &gt; &gt; &gt; potential of core electrons?<BR>&gt; &gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; &gt; Who knows how to get the core level of each atom after the scf<BR>&gt; &gt; &gt; &gt; &gt; &gt; calculation?<BR>&gt; &gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; &gt; Thanks in advance!<BR>&gt; &gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; &gt; Rgds<BR>&gt; &gt; &gt; &gt; &gt; &gt; Jianchun<BR>&gt; &gt; &gt; &gt; &gt; &gt; Department of physics<BR>&gt; &gt; &gt; &gt; &gt; &gt; Soochow university China.215006<BR>&gt; &gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; &gt; _________________________________________________________________<BR>&gt; &gt; &gt; &gt; &gt; &gt; MSNÊ®Äê»ØÀ¡£¬Ã¿Î»Óû§¿ÉÃâ·Ñ»ñµÃ¼ÛÖµ25ÔªµÄ¿¨°Í˹»ù·´²¡¶¾Èí¼þ2010¼¤»îÂ룬¿ìÀ´ÁìÈ¡£¡<BR>&gt; &gt; &gt; &gt; &gt; &gt; http://kaba.msn.com.cn/?k=1<BR>&gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; --<BR>&gt; &gt; &gt; &gt; &gt; ********************************************************<BR>&gt; &gt; &gt; &gt; &gt; - Article premier - Les hommes naissent et demeurent<BR>&gt; &gt; &gt; &gt; &gt; libres et ¨¨gaux en droits. Les distinctions sociales<BR>&gt; &gt; &gt; &gt; &gt; ne peuvent ¨ºtre fond¨¨es que sur l'utilit¨¨ commune<BR>&gt; &gt; &gt; &gt; &gt; - Article 2 - Le but de toute association politique<BR>&gt; &gt; &gt; &gt; &gt; est la conservation des droits naturels et<BR>&gt; &gt; &gt; &gt; &gt; imprescriptibles de l'homme. Ces droits sont la libert¨¨,<BR>&gt; &gt; &gt; &gt; &gt; la propri¨¨t¨¨, la s¨´ret¨¨ et la r¨¨sistance ¨¤ l'oppression.<BR>&gt; &gt; &gt; &gt; &gt; ********************************************************<BR>&gt; &gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; &gt; Giuseppe Mattioli<BR>&gt; &gt; &gt; &gt; &gt; CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<BR>&gt; &gt; &gt; &gt; &gt; v. Salaria Km 29,300 - C.P. 10<BR>&gt; &gt; &gt; &gt; &gt; I 00016 - Monterotondo Stazione (RM)<BR>&gt; &gt; &gt; &gt; &gt; Tel + 39 06 90672836 - Fax +39 06 90672316<BR>&gt; &gt; &gt; &gt; &gt; E-mail: &lt;giuseppe.mattioli@ism.cnr.it&gt;<BR>&gt; &gt; &gt; &gt; &gt; _______________________________________________<BR>&gt; &gt; &gt; &gt; &gt; Pw_forum mailing list<BR>&gt; &gt; &gt; &gt; &gt; Pw_forum@pwscf.org<BR>&gt; &gt; &gt; &gt; &gt; http://www.democritos.it/mailman/listinfo/pw_forum<BR>&gt; &gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; ËÑË÷±¾Ó¦ÊDzÊÉ«µÄ,¿ìÀ´ÌåÑéÐÂÒ»´úËÑË÷ÒýÇæ-±ØÓ¦,¾«ÃÀͼƬÿÌ컻Ŷ! Á¢¼´ÊÔÓã¡<BR>&gt; &gt; &gt; &gt; _________________________________________________________________<BR>&gt; &gt; &gt; &gt; SkyDriveµç×Ó»­²á£¬´øÄãÁìÂÔ¾«²ÊÕÕÆ¬£¬·ÖÏí¡°ÃÀ¡±Ê±¡°ÃÀ¡±¿Ì!<BR>&gt; &gt; &gt; &gt; http://www.windowslive.cn/campaigns/e-magazine/ngmchina/?a=c<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; --<BR>&gt; &gt; &gt; ********************************************************<BR>&gt; &gt; &gt; - Article premier - Les hommes naissent et demeurent<BR>&gt; &gt; &gt; libres et ¨¨gaux en droits. Les distinctions sociales<BR>&gt; &gt; &gt; ne peuvent ¨ºtre fond¨¨es que sur l'utilit¨¨ commune<BR>&gt; &gt; &gt; - Article 2 - Le but de toute association politique<BR>&gt; &gt; &gt; est la conservation des droits naturels et<BR>&gt; &gt; &gt; imprescriptibles de l'homme. Ces droits sont la libert¨¨,<BR>&gt; &gt; &gt; la propri¨¨t¨¨, la s¨´ret¨¨ et la r¨¨sistance ¨¤ l'oppression.<BR>&gt; &gt; &gt; ********************************************************<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; Giuseppe Mattioli<BR>&gt; &gt; &gt; CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<BR>&gt; &gt; &gt; v. Salaria Km 29,300 - C.P. 10<BR>&gt; &gt; &gt; I 00016 - Monterotondo Stazione (RM)<BR>&gt; &gt; &gt; Tel + 39 06 90672836 - Fax +39 06 90672316<BR>&gt; &gt; &gt; E-mail: &lt;giuseppe.mattioli@ism.cnr.it&gt;<BR>&gt; &gt; &gt; _______________________________________________<BR>&gt; &gt; &gt; Pw_forum mailing list<BR>&gt; &gt; &gt; Pw_forum@pwscf.org<BR>&gt; &gt; &gt; http://www.democritos.it/mailman/listinfo/pw_forum<BR>&gt; &gt;<BR>&gt; &gt; _________________________________________________________________<BR>&gt; &gt; Ô¼»á˵²»Ç嵨·½£¿À´ÊÔÊÔ΢ÈíµØÍ¼×îÐÂmsn»¥¶¯¹¦ÄÜ£¡<BR>&gt; &gt; http://ditu.live.com/?form=TL&amp;swm=1<BR>&gt; <BR>&gt; <BR>&gt; <BR>&gt; -- <BR>&gt; ********************************************************<BR>&gt; - Article premier - Les hommes naissent et demeurent<BR>&gt; libres et ¨¨gaux en droits. Les distinctions sociales<BR>&gt; ne peuvent ¨ºtre fond¨¨es que sur l'utilit¨¨ commune<BR>&gt; - Article 2 - Le but de toute association politique<BR>&gt; est la conservation des droits naturels et <BR>&gt; imprescriptibles de l'homme. Ces droits sont la libert¨¨,<BR>&gt; la propri¨¨t¨¨, la s¨´ret¨¨ et la r¨¨sistance ¨¤ l'oppression.<BR>&gt; ********************************************************<BR>&gt; <BR>&gt; &nbsp; &nbsp;Giuseppe Mattioli &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<BR>&gt; &nbsp; &nbsp;CNR - ISTITUTO DI STRUTTURA DELLA MATERIA &nbsp; <BR>&gt; &nbsp; &nbsp;v. Salaria Km 29,300 - C.P. 10 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<BR>&gt; &nbsp; &nbsp;I 00016 - Monterotondo Stazione (RM) &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<BR>&gt; &nbsp; &nbsp;Tel + 39 06 90672836 - Fax +39 06 90672316 &nbsp; &nbsp;<BR>&gt; &nbsp; &nbsp;E-mail: &lt;giuseppe.mattioli@ism.cnr.it&gt;<BR>&gt; _______________________________________________<BR>&gt; Pw_forum mailing list<BR>&gt; Pw_forum@pwscf.org<BR>&gt; http://www.democritos.it/mailman/listinfo/pw_forum<BR><BR>                                               <br /><hr />ËÑË÷±¾Ó¦ÊDzÊÉ«µÄ,¿ìÀ´ÌåÑéÐÂÒ»´úËÑË÷ÒýÇæ-±ØÓ¦,¾«ÃÀͼƬÿÌ컻Ŷ! <a href='http://cn.bing.com/?form=CRMADS ' target='_new'>Á¢¼´ÊÔÓã¡</a></body>
</html>