<DIV>Dear Paolo:</DIV>
<DIV>       I have calculated the band structure with lattice constant <FONT color=#ff0000>5.64613 <STRONG>(A)</STRONG></FONT><FONT color=#000000> , and set  calculation = 'relax'. Then I get the band gap of the germanium 0.2269 eV. It's  also too large.</FONT></DIV>
<DIV>      I tried to optimize band structure, just like CO in example3.</DIV>
<DIV>input file :</DIV>
<DIV>&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&</DIV>
<DIV><FONT color=#000000>       &control<BR>    calculation = 'relax'<BR>    restart_mode='from_scratch',<BR>    prefix='germaniun',<BR>    tstress = .true.<BR>    tprnfor = .true.<BR>    pseudo_dir = '/usr/people/lty/code/pseudo//',<BR>    outdir='/usr/people/lty/tmp//'<BR> /<BR> &system<BR>    ibrav=  0, nat=  2, ntyp= 1,<BR>    ecutwfc =30.0, nbnd = 8<BR> /<BR> &electrons<BR>    diagonalization='cg'<BR>    mixing_mode = 'plain'<BR>    mixing_beta = 0.7<BR>    conv_thr =  1.0d-8<BR> /<BR> &ions<BR>    upscale = 10.0<BR> /<BR>ATOMIC_SPECIES<BR> Ge  72.59  Ge.pz-bhs.UPF<BR>ATOMIC_POSITIONS<BR> Ge 0.00 0.00 0.00 1 1 1<BR> Ge 0.25 0.25 0.25 1 1
 1<BR>CELL_PARAMETERS<BR> 10.0 00.0 00.0<BR> 00.0 10.0 00.0<BR> 00.0 00.0 10.0<BR>K_POINTS<BR>  6<BR>   0.0  0.0  0.0   1.00<BR>   0.1  0.1  0.1   1.00<BR>   0.2  0.2  0.2   1.00<BR>   0.3  0.3  0.3   1.00<BR>   0.4  0.4  0.4   1.00<BR>   0.5  0.5  0.5   1.00<BR>&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&</FONT></DIV>
<DIV><FONT color=#000000> Then I get the band gap :</FONT></DIV>
<DIV><FONT color=#000000>***********************************************************************</FONT></DIV><FONT color=#000000>
<DIV><BR>          k = 0.0000 0.0000 0.0000 ( 2777 PWs)   bands (ev):</DIV>
<DIV>   -11.1883  -7.1028  -2.6627  -1.6724  -1.6724  -0.1296  -0.1296   4.2633</DIV>
<DIV>          k = 0.1000 0.1000 0.1000 ( 2785 PWs)   bands (ev):</DIV>
<DIV>   -11.1650  -7.1625  -2.5315  -1.8738  -1.8738   0.1408   0.1408   2.0930</DIV>
<DIV>          k = 0.2000 0.2000 0.2000 ( 2785 PWs)   bands (ev):</DIV>
<DIV>   -11.1051  -7.2871  -2.3407  -2.2155  -2.2155   0.5997   0.5997   1.8475</DIV>
<DIV>          k = 0.3000 0.3000 0.3000 ( 2785 PWs)   bands (ev):</DIV>
<DIV>   -11.0338  -7.3931  -2.4812  -2.4812  -2.3644   0.7896   0.7896   2.4967</DIV>
<DIV>          k = 0.4000 0.4000 0.4000 ( 2764 PWs)   bands (ev):</DIV>
<DIV>   -10.9785  -7.4498  -2.6302  -2.6302  -2.4865   0.6643   0.6643   3.4208</DIV>
<DIV>          k = 0.5000 0.5000 0.5000 ( 2800 PWs)   bands (ev):</DIV>
<DIV>   -10.9579  -7.4662  -2.6758  -2.6758  -2.5383   0.5696   0.5696   3.0471</DIV>
<DIV>&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&</DIV>
<DIV>The band-gap is -0.8659 eV. I think the band-gap is error. Which must be adjusted in my input-file, then I get the right band-gap ?</DIV>
<DIV> </DIV>
<DIV>Best wishes</DIV>
<DIV> </DIV>
<DIV>tylv<BR><BR><B><I>Paolo Giannozzi <giannozz@nest.sns.it></I></B> wrote:</FONT></DIV>
<BLOCKQUOTE class=replbq style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #1010ff 2px solid">On Saturday 31 July 2004 11:04, Bachir Bouhafs wrote:<BR><BR>> [...] I think that you have calculated your band structure at the<BR>> experimental lattice parameter [...] which is greater than the LDA <BR>> value. Which means, that you have applied a negative pressure <BR>> to your crystal, and consequently the band gap increases.<BR><BR>[ please do not post HTML ]<BR><BR>and this is likely to be the main reason of the discrepancy. Other <BR>factors, such as the pseudopotential, cutoff, number of k-points, <BR>may contribute to the difference. For instance: in GaAs, the use <BR>of a nonrelativistic pseudopotential for Ga increases the band gap<BR>by approx. 0.4 eV.<BR><BR>By the way: the tests and examples are NOT meant to give any<BR>realistic description of any realistic systems. <BR><BR>Paolo<BR>-- <BR>Paolo Giannozzi e-mail: giannozz@nest.sns.it<BR>Scuola Normale
 Superiore Phone: +39/050-509876, Fax:-563513 <BR>Piazza dei Cavalieri 7 I-56126 Pisa, Italy<BR>_______________________________________________<BR>Pw_forum mailing list<BR>Pw_forum@pwscf.org<BR>http://www.democritos.it/mailman/listinfo/pw_forum<BR></BLOCKQUOTE><p><br><hr size=1><b>Do You Yahoo!?</b><br>
<a href="http://music.yisou.com" target=blank>150万曲MP3疯狂搜，带您闯入音乐殿堂</a><br><a href="http://image.yisou.com" target=blank>美女明星应有尽有，搜遍美图、艳图和酷图</a><br>
<a href="http://cn.rd.yahoo.com/mail_cn/tag/1g/*http://cn.mail.yahoo.com/event/mail_1g/" target=blank>1G就是1000兆，雅虎电邮自助扩容！</a>