²é¿´: 11246  |  »Ø¸´: 274
µ±Ç°Ö÷ÌâÒѾ­´æµµ¡£

wuli8

ÈÙÓþ°æÖ÷ (ÖªÃû×÷¼Ò)

¡­¡­¡­¡­

ÓÅÐã°æÖ÷ÓÅÐã°æÖ÷

¡ï ¡ï ¡ï
csfn(½ð±Ò+3,VIP+0):many thanks, ÐÁ¿àÁË
Ask:
ÒýÓûØÌû:
I want to plot Band structure by use
Case.Spaghetti_ene file but I don't know, what are the five columns in this
file? Could you please let me know, how can I plot Band structure by
Case.Spaghrtti_ene file and or other files?

Answer
ÒýÓûØÌû:
The first three columns correspond to coordinates of the particular k-point,
the 4th column is the total distance of your path in the k-space, and the
last column is the actual energy corresponding to this k-point and energy
band. Bands are listed below each other.

To plot the band structure you need only the 4th and 5th columns. The first
three are there mainly for you orientation and to allow you to label
important points in your band structure correctly.

¡­¡­¡­¡­
121Â¥2008-12-28 12:04:57
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wysong

½ð³æ (СÓÐÃûÆø)

¡ï ¡ï ¡ï ¡ï
wuli8(½ð±Ò+4,VIP+0):¸ÐлÄúµÄ½»Á÷£¡(*^__^*) ¡­¡­
ÒýÓûØÌû:
Originally posted by zhangyunli at 2008-12-25 10:43:

ºÃÏñ²»ÐÐ

ллÄãµÄ´ð¸´ ÎÒÊÇÕâÑù×öµÄ£ºÊ×ÏÈÓÃFSM³õ²½´ïµ½Ò»¸ö·´Ìú´ÅµÄ×ÔÇ¢½á¹û È»ºóÔËÐÐÕý³£µÄGGA ÔÚÕâ¸ö½á¹ûµÄ»ù´¡ÉÏÃæÎҾͿÉÒÔ¼Ó¹ìµÀñîºÏÁË ºÇºÇ
122Â¥2008-12-28 14:56:20
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wuli8

ÈÙÓþ°æÖ÷ (ÖªÃû×÷¼Ò)

¡­¡­¡­¡­

ÓÅÐã°æÖ÷ÓÅÐã°æÖ÷

ÒýÓûØÌû:
Originally posted by wysong at 2008-12-28 14:56:



ллÄãµÄ´ð¸´ ÎÒÊÇÕâÑù×öµÄ£ºÊ×ÏÈÓÃFSM³õ²½´ïµ½Ò»¸ö·´Ìú´ÅµÄ×ÔÇ¢½á¹û È»ºóÔËÐÐÕý³£µÄGGA ÔÚÕâ¸ö½á¹ûµÄ»ù´¡ÉÏÃæÎҾͿÉÒÔ¼Ó¹ìµÀñîºÏÁË ºÇºÇ

¸Ðл½»Á÷
¡­¡­¡­¡­
123Â¥2008-12-28 15:05:19
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wuli8

ÈÙÓþ°æÖ÷ (ÖªÃû×÷¼Ò)

¡­¡­¡­¡­

ÓÅÐã°æÖ÷ÓÅÐã°æÖ÷

¡ï ¡ï ¡ï
csfn(½ð±Ò+3,VIP+0):many thanks, ÐÁ¿àÁË
Ask
ÒýÓûØÌû:
I read userguide and plot my band structure according to that.
> The process is
> after a successful scf
>   making case.klist_band
>    x lapw1 -band
>    x lapw2 -qtl -band
>    editing case.insp
>    x spaghetti
>    plot ...
>
> my questions:
>
> 1- I saw case.energy but I want to know from which file I can
> calculate band gap and band width exactly  and how.
>
> 2. Is there any difference between  bandstructure of one case with the
> same case.klist_band but different case.klist ?
>

Answer
ÒýÓûØÌû:
1.A possible (probably not the fastest) route is as follows :
   Plot the spaghetti, to start with.
   Get the K0-value at which the gap occurs (eg. GAMMA point)  
   Get EFERMI from case.outputt.
   Examine (eigen)energy values at K0 in case.output1 : The energy
difference between EFERMI and the
   next E value is the required energy gap.
   Similarly, knowing the k-values where band extrema (min,max) occur (
see spaghetti plot),
   you can obtain the bandwidth from examination of E(K) values in
case.output1.

2. YES ! Increasing the case.klist set, in principle, improves accuracy
of SC potential.
    See UG  pg. 19 - "x kgen"

[ Last edited by wuli8 on 2008-12-28 at 17:05 ]
¡­¡­¡­¡­
124Â¥2008-12-28 17:02:37
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wuli8

ÈÙÓþ°æÖ÷ (ÖªÃû×÷¼Ò)

¡­¡­¡­¡­

ÓÅÐã°æÖ÷ÓÅÐã°æÖ÷

¡ï ¡ï ¡ï
csfn(½ð±Ò+3,VIP+0):many thanks, ÐÁ¿àÁË
Ask
ÒýÓûØÌû:
Do you mean I should use the imaginary and real part values in the file
case.epsilon to calculate the reflectivity? Today I tried to calculate the
reflectivity of anatase TiO2, but the obtained values of reflectivity are
less than 35% in the range of 200-900nm. While my experimental values are
more than 80%. What's the matter?

Answer
ÒýÓûØÌû:
> Do you mean I should use the imaginary and real part values
> in the file case.epsilon to calculate the reflectivity?

I mean only what is written at the top of page 131 of the
usersguide.  That is, "All dielectric tensor components can
be found in file case.epsilon[up|dn]."

> Today I tried to calculate the reflectivity of anatase TiO2, but
> the obtained values of reflectivity are less than 35% in the
> range of 200-900nm. While my experimental values are
> more than 80%. What's the matter?

The reflectivity comes after a lot of post-processing.
Why not go back to the more "raw" theoretical results.

I suggest you do the following:
1) find previous DFT-based calculations on anatase TiO2.
One that seems particulary relevant is PRB 51, 13023 (1995).
It even mentions reflectivity calculations.
2) Compare your band structure, DOS and computed dieletric
functions to the ones in the above paper.
3) Check your scissor-shift.  Did you use a scissor shift?  Note
that (for normal incidence) the reflectivity is
        r = (n-1)^2/(n+1)^2.
A scissor shift lowers the value of epsilon(0), which lowers
the theoretical value of the reflectivity.

¡­¡­¡­¡­
125Â¥2008-12-28 17:23:49
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wysong

½ð³æ (СÓÐÃûÆø)

ÇóÖú

Ê×ÏÈÏò°æÖ÷µÀǸ ÔÚÕâÀïÎÊÒ»¸öÓëʹÓÃwien2kÎ޹صÄÎÊÌ⠺ǺÇ

ÎÒÏëÈø÷ΪÓÃwien2kµÄÕ½ÓÑÃǰïæÁ˽âһϠÓÐÔÚÀ¼ÖݵÄͬѧÂð À¼ÖÝÓÐÓÃwien2k»òÕßvaspÀ´×ö²ÄÁÏÉè¼ÆµÄÑо¿µ¥Î»Âð  Âé·³´ó¼ÒÁË ºÇºÇ
126Â¥2008-12-28 19:15:38
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wuli8

ÈÙÓþ°æÖ÷ (ÖªÃû×÷¼Ò)

ÒýÓûØÌû:
Originally posted by wysong at 2008-12-28 19:15:
Ê×ÏÈÏò°æÖ÷µÀǸ ÔÚÕâÀïÎÊÒ»¸öÓëʹÓÃwien2kÎ޹صÄÎÊÌ⠺ǺÇ

ÎÒÏëÈø÷ΪÓÃwien2kµÄÕ½ÓÑÃǰïæÁ˽âһϠÓÐÔÚÀ¼ÖݵÄͬѧÂð À¼ÖÝÓÐÓÃwien2k»òÕßvaspÀ´×ö²ÄÁÏÉè¼ÆµÄÑо¿µ¥Î»Âð  Âé·³´ó¼ÒÁË ºÇºÇ

127Â¥2008-12-28 22:28:56
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wysong

½ð³æ (СÓÐÃûÆø)

ÇóÖú

¸÷λ´óÏÀ ÎÒÏëÓÃxcrsdenÀ´×öµç×ÓÃܶÈͼ ¿´¿´ÎÒµÄÌåϵÊÇ·ñÓйìµÀÓÐÐòÏÖÏ󣬵«Êǵõ½µÄͼһֱÏÔʾÊÇÇòÐÎ ºÜÃ÷ÏÔÊÇÓÐÎÊÌâ ÎÒÒѾ­ÔÚin2ÎļþÖа´ÕÕDOSÉÏÏÔʾµÄÄÜÁ¿Çø¼äÑ¡ÔñÁËÎÒÏëÒªÏÔʾԭ×ÓµÄÇø¼ä ²»ÖªµÀ¸ÃÔõô°ìÁË ºÇºÇ ÓÐÄÇλ¸ßÊÖ°ï°ïæ°¡ ллÁË Õâ¸ö¶«Î÷ºÜÂé·³ ÄÜÓÃQQÖ¸µãÒ»ÏÂÂð ÎÒµÄÊÇ292222016 ·Ç³£¸Ðл£¡£¡
128Â¥2008-12-29 15:36:25
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wax54

Í­³æ (СÓÐÃûÆø)

¡ï ¡ï ¡ï ¡ï ¡ï
wuli8(½ð±Ò+5,VIP+0):¸ÐлÄúµÄ½»Á÷£¡(*^__^*) ¡­¡­
ʹÓÃwien2kʱºòÐèÒªµÄ¿´Í¼¹¤¾ßÊDz»ÊÇÒ»¹²ÓÐÈý¸ö°¡£¬°üÀ¨
xcrysden£ºÓÃÓڲ鿴¾§¸ñ½á¹¹ºÍµçºÉÃܶÈ
gnuplot£ºÓÃÓڲ鿴DOSºÍµçºÉÃܶÈ
ghostview£ºÓÃÓڲ鿴ÄÜ´ø½á¹¹
°ßÖñ£¬ÎÒµÄÀí½âÊÇ·ñÕùÈ¡£¿
129Â¥2008-12-30 19:39:43
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wax54

Í­³æ (СÓÐÃûÆø)

´ò´íÁË£¬ÎÒÊÇÏë˵ÊÇ·ñÕýÈ·£¿
130Â¥2008-12-30 19:40:30
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ wuli8 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 085700×ÊÔ´Óë»·¾³308Çóµ÷¼Á +3 īīĮ 2026-03-18 3/150 2026-03-18 22:35 by bingxueer79
[¿¼ÑÐ] ²ÄÁÏרҵÇóµ÷¼Á +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] 295Çóµ÷¼Á +3 Ò»Ö¾Ô¸¾©Çø211 2026-03-18 5/250 2026-03-18 17:03 by zhaoqian0518
[¿¼ÑÐ] 085601ר˶£¬×Ü·Ö342Çóµ÷¼Á£¬µØÇø²»ÏÞ +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤Ò»Ö¾Ô¸Äϲý´óѧ327Çóµ÷¼ÁÍÆ¼ö +8 Ncdx123456 2026-03-13 9/450 2026-03-18 14:40 by haxia
[¿¼ÑÐ] 299Çóµ÷¼Á +5 ¡÷С͸Ã÷* 2026-03-17 5/250 2026-03-18 11:49 by ¾¡Ë´Ò¢1
[¿¼ÑÐ] 303Çóµ÷¼Á +4 î£08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[¿¼ÑÐ] 290Çóµ÷¼Á +6 ¿×Ö¾ºÆ 2026-03-12 11/550 2026-03-17 14:41 by ÖÜÖÛÖÛ77
[¿¼ÑÐ] 275Çóµ÷¼Á +4 Ì«Ñô»¨ÌìÌ쿪ÐÄ 2026-03-16 4/200 2026-03-17 10:53 by ¹¦·ò·è¿ñ
[¿¼ÑÐ] 11408 Ò»Ö¾Ô¸Î÷µç£¬277·ÖÇóµ÷¼Á +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[¿¼ÑÐ] 0854¿ØÖƹ¤³Ì 359Çóµ÷¼Á ¿É¿çרҵ +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[¿¼ÑÐ] 326Çóµ÷¼Á +4 ŵ±´¶û»¯Ñ§½±êéê 2026-03-15 7/350 2026-03-16 17:11 by ŵ±´¶û»¯Ñ§½±êéê
[¿¼ÑÐ] 0703Ò»Ö¾Ô¸211 285·ÖÇóµ÷¼Á +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by ŶŶ123
[¿¼ÑÐ] 283Çóµ÷¼Á +10 С¥¡£ 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[¿¼ÑÐ] 327Çóµ÷¼Á +6 ʰ¹âÈÎȾ 2026-03-15 11/550 2026-03-15 22:47 by ʰ¹âÈÎȾ
[¿¼ÑÐ] 080500£¬²ÄÁÏѧ˶302·ÖÇóµ÷¼ÁѧУ +4 ³õʶ¿ÉÀÖ 2026-03-14 5/250 2026-03-14 21:08 by peike
[¿¼ÑÐ] 297Ò»Ö¾Ô¸ÉϽ»085600Çóµ÷¼Á +5 Ö¸¼â°ËǧÀï 2026-03-14 5/250 2026-03-14 17:26 by a²»Ò×
[¿¼ÑÐ] 0817»¯Ñ§¹¤³ÌÓë¼¼Êõ¿¼ÑÐ312·Öµ÷¼Á +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[¿¼ÑÐ] 290Çóµ÷¼Á +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[¿¼ÑÐ] 333Çóµ÷¼Á +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û