²é¿´: 5721  |  »Ø¸´: 23

lu_wei07

Òø³æ (СÓÐÃûÆø)

[½»Á÷] PWSCF ¾ßÌåʵÀý Ïêϸ½â´ð¸÷ÖÖÎÊÌâ ÒÑÓÐ13È˲ÎÓë

PWSCF ¾ßÌåʵÀý Ïêϸ½â´ð¸÷ÖÖÎÊÌ⣬ժ×ÔÍøÂçÒ»¸öblogblog µØÖ·
http://zhaoyunawuli.blog.163.com ... 976200811147143279/
ÍøÒ׵Ġ ÍøÃû£º ËæÃÎÏë¶ø·É
1,¸½ ´¦ÀíPWSCFÖмÆËãµÃµ½±¾Õ÷ƵÂʵijÌÐò   ×Ô¼ºÐÞ¸ÄÁË
ÓõÄʱºòÐèÒª°Ñ*.freq µÄÎļþ ÊÖ¶¯¸Ä³Éinput.freq
program prog
real, allocatable :: e(:,
real, allocatable :: k(:, ,dk( ,x(
real, dimension(3) ::k0,a
character(len=32):: xx, yy, zz
integer ::   nbands, i, j, n, nk




open(10, file='input.freq', status='old')
read(10,*) xx, yy, nbands, zz, nk
write(*,*)"nk=",nk

allocate(k(nk,3))
allocate(e(nk,nbands))
allocate(dk(0:nk))
allocate(x(1:nk))

open(11,file='phoband.dat')
    do i=1,nk

    read(10,*) (k(i,j),j=1,3)
    read(10,*) (e(i,n),n=1,nbands)
   
        enddo
         
                do n=1,nk
           dk(0)=0
            if (n.eq.1) then
        k0=k(n,
            endif
      a=k(n,-k0
      dk(n)=dk(n-1)+sqrt(dot_product(a,a))
      x(n)=dk(n)     
      k0=k(n,
      enddo

      do i=1,nk
        write(11, 8888) x(i), (e(i,j),j=1,nbands)
      enddo
8888 format (1x, f9.4, 2x ,1000(f12.6,2x))

end program prog
2,´¦ÀíÄÜ´øµÄ³ÌÐòÖÐÔÚ do i=1, nhighk-1 ǰһÐÐ ¼ÓÉÏ
   ntkp=0
ÆäËüµÄÎÒÔÝʱ»¹Ã»Óõ½£¬ÓÐʲôÎÊÌâÏ£Íû¹²Í¬½»Á÷[ Last edited by lu_wei07 on 2013-3-5 at 16:37 ]
»Ø¸´´ËÂ¥

» ±¾Ìû¸½¼þ×ÊÔ´Áбí

  • »¶Ó­¼à¶½ºÍ·´À¡£ºÐ¡Ä¾³æ½öÌṩ½»Á÷ƽ̨£¬²»¶Ô¸ÃÄÚÈݸºÔð¡£
    ±¾ÄÚÈÝÓÉÓû§×ÔÖ÷·¢²¼£¬Èç¹ûÆäÄÚÈÝÉæ¼°µ½ÖªÊ¶²úȨÎÊÌ⣬ÆäÔðÈÎÔÚÓÚÓû§±¾ÈË£¬Èç¶Ô°æÈ¨ÓÐÒìÒ飬ÇëÁªÏµÓÊÏ䣺xiaomuchong@tal.com
  • ¸½¼þ 1 : PWSCF½Ì³Ì£¬´ÓÈëÃŵ½¾«Í¨.docx
  • 2013-03-04 20:04:39, 45.86 K
  • ¸½¼þ 2 : phfreq.exe
  • 2013-03-05 16:37:30, 349 K

» ÊÕ¼±¾ÌûµÄÌÔÌûר¼­ÍƼö

¾«»ªÍøÌûÊÕ¼¯ ¿ÆÑм¼ÇÉ µÚÒ»ÐÔÔ­Àí¼ÆË㸨Öú¹¤¾ß condensed matter physics
¼ÆËãÄ£Äâ µÚÒ»ÐÔÔ­Àí¼ÆËã ¼ÆËãºÍÈí¼þ

» ±¾ÌûÒÑ»ñµÃµÄºì»¨£¨×îÐÂ10¶ä£©

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÕâÖ»ÊÇÒ»³¡ÓÎÏ·
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
»ØÌûÖ§³Ö ( ÏÔʾ֧³Ö¶È×î¸ßµÄǰ 50 Ãû )

fyuewen

½ð³æ (ÕýʽдÊÖ)

¡ï
fzx2008: ½ð±Ò+1, лл»ØÌû£¡ 2013-03-05 16:37:31
fzx2008: »ØÌûÖö¥ 2013-03-05 16:37:35
ͬѧ£¬Çë×¢ÖØÖªÊ¶²úȨ£¬Ä㶼˵ÁËÊDz©¿ÍÀïÃæµÄ£¬ÄǾͰѲ©Ö÷Ìù³öÀ´¡£Ã»¿´ÎҾͲ²âÊǺîÖù·åÀÏʦµÄ£¬´ò¿ª¿´¹ûÈ»ÊÇ¡£¶øÇÒÄã·ÖÏíµÄÕâÖÖ·½Ê½×Åʵ²»ºÃ£¬ÄãºÃ´õ˵һЩ¸ÅÀ¨ÔÚÌù×ÓÀïÃæ¡£Ìù³öÒ»²¿·ÖÔ­ÎÄ¡£¡£¡£ÆäʵÕâ¸ö¶«Î÷Ö»Òª´ó¼ÒÁ¬½ÓµÄÁËÍøÂç¾Í¶¼¿ÉÒÔÔÚÍøÒ³ÉÏÖ±½Ó¿´µÄ¡£
ÎÞÐǼ¶£¬²»ÊÇÆÀÂÛ×ÊÔ´¡£¶øÊǸøÄãÕâ¸öÌû×Ó¡£
Physicsloveukulele!Aloha!
7Â¥2013-03-05 13:34:43
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
ÆÕͨ»ØÌû

0.8

¶¥Ò»Ï£¬¸Ðл·ÖÏí£¡
2Â¥2013-03-04 21:23:36
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

hakuna

ľ³æ (ÖªÃû×÷¼Ò)

0.8

×ܽá·Ç³£È«Ã桢ϸÖ£¬ÖµµÃÍÆ¼ö£¡
3Â¥2013-03-05 06:49:10
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

0.8

лл·ÖÏí£¡£¡£¡£¡£¡£¡£¡£¡£¡
4Â¥2013-03-05 06:56:35
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

0.8

¶¥Ò»Ï£¬¸Ðл·ÖÏí£¡
5Â¥2013-03-05 09:15:32
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

0.8

¶¥Ò»Ï£¬¸Ðл·ÖÏí£¡
6Â¥2013-03-05 09:19:03
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

Ц`Ц

Ìú³æ (³õÈëÎÄ̳)

0.8

Õý°Ù˼²»µÃÆä½â£¬ÏÈÔÞÒ»¸ö°É
ÓÐÅó×ÔÔ¶·½À´~~
8Â¥2013-03-05 15:05:23
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

0.8

¶¥Ò»Ï£¬¸Ðл·ÖÏí£¡
9Â¥2013-03-05 16:02:02
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

huazhorg

Ìú¸Ëľ³æ (ÕýʽдÊÖ)

¡ï
Сľ³æ: ½ð±Ò+0.5, ¸ø¸öºì°ü£¬Ð»Ð»»ØÌû
ÎÞÂÛÈçºÎ£¬»¹ÊÇÓеã·ÖÏí¾«Éñ£¬ÔÛÃÇÕâ¸ö°æ¿é˵ÆðÀ´ÀäÇ壬¾ÍÊÇȱÉÙÕâÖÖÉáµÃ·ÖÏíµÄÈË£¬´ó¼Ò¶¼ÓÐʲô¶¼²ØÔÚ¶Ç×ÓÀҲ²»ÊÇʲôºÃÊ¡£¶¥Ò»¸ö¡£
10Â¥2013-03-05 17:20:12
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ lu_wei07 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[½Ìʦ֮¼Ò] µ¼Ê¦Í²ۣºÎÒÔõô̯ÉÏÁËÕâô¸ö¼«Æ·Ñо¿Éú£¡ +4 ËÕ¶«ÆÂ¶þÊÀ 2026-08-23 4/200 2026-08-27 09:18 by Ú¤Íõƪ
[»ù½ðÉêÇë] ¹ú×ÔÈ»ÃæÉϸ´ÅÌ~»¶Ó­ÌÖÂÛ +7 ÇçÌì¼ÓÓÍ 2026-08-26 8/400 2026-08-27 07:53 by longfie172
[˶²©¼ÒÔ°] ÊÛSCIÒ»ÇøT0PÎÄÕ£¬ÎÒ:8.O.55.1.O.54,¿ÆÄ¿ÆëÈ«,¿É+¼± +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[»ù½ðÉêÇë] ÎÒ²»Àí½â£¡ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[»ù½ðÉêÇë] ·¶½øÖоÙÒ»ÎĵÄÖÐÐÄ˼Ïë +9 Ñ׻ƹóëÐ 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[»ù½ðÉêÇë] 2026¹ú×ÔÈ»º¯ÆÀ·Ñµ½ÕË +22 ÑòÑü°å 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[»ù½ðÉêÇë] 2026Äê8ÔÂ25ÈÕ¹ú×ÔÈ»·Å°ñǰͻȻÊÕµ½ÁÐÈëÆÀÉóר¼ÒÓʼþ£¬ÓйØÏµÂ𣿠+25 ľˮ˼¶¹ 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[»ù½ðÉêÇë] ÄÜ·ñÍ˳ö²ÎÓëµÄÃæÉÏÏîÄ¿½â³ýÏÞÏî +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by ±¦±´³æ×Ó
[»ù½ðÉêÇë] ÄÄλ¸ßÈËÖÐÁË£¬°Ñ²éѯµ½µÄ½ØÍ¼Ìù³öÀ´ÈÃÎÒ¿´¿´£¬ÈÃÎÒ³¤³¤¼ûʶ +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[»ù½ðÉêÇë] Ϊʲô¹ú×ÔÈ»²»ÄÜÖ±½Ó¹«²¼ +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[»ù½ðÉêÇë] ¹úºÏÀïÃæÄÜ¿´µ½ÁË +7 Ò»»³Ü°Çï 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[»ù½ðÉêÇë] ÏîÄ¿ÐÅÏ¢ºÍ¾­·ÑÐÅÏ¢ÔÚϵͳÀï¶¼¿ÉÒÔ¿´µ½ÁË +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[»ù½ðÉêÇë] ¹úºÏ¿É²éÁË +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[»ù½ðÉêÇë] ½ñÌìÎñί»á¿ªÍêÁË£¬Ã÷Ìì³ö½á¹ûÂð +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[Ö°³¡ÈËÉú] ѧÉúײ¼û¸¨µ¼Ô±ËÍÍâÂô£¬µÚ¶þÌìÈ«°à¶¼³ÁĬÁË +3 ¾¨ÓãÈÚ½ð_Õã½­_É 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[»ù½ðÉêÇë] ÔÚ¼á±ù»¹¸Ç×ű±º£µÄʱºò£¬ÎÒ¿´µ½ÁËÅ­·ÅµÄ÷»¨¡£ (½ð±Ò+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[»ù½ðÉêÇë] Ã÷ÌìÓ¦¸Ã¿É²éÁË£¡£¿ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[»ù½ðÉêÇë] ûÓÐÈκÎÏûÏ¢-ÊDz»ÊǾÍÁ¹ÁË +9 ͼÀ²Í¼À² 2026-08-24 10/500 2026-08-25 11:59 by ÄϺ£Ð¡¸ç
[»ù½ðÉêÇë] ¿´À´½ñÌì²»»á·Å°ñÁË£¿ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[»ù½ðÉêÇë] »ù½ð°¡»ù½ð +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û