²é¿´: 5165  |  »Ø¸´: 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µÄ»ØÌû

0.8

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

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µÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ²ÄÁÏѧ˶297ÒѹýËÄÁù¼¶Çóµ÷¼ÁÍÆ¼ö +6 adaie 2026-03-19 6/300 2026-03-20 14:12 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] 290Çóµ÷¼Á +4 ^O^Ø¿ 2026-03-19 4/200 2026-03-20 11:51 by ѧԱ8dgXkO
[¿¼ÑÐ] ²ÄÁÏѧ˶318Çóµ÷¼Á +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á +10 ÄÝÄÝninicgb 2026-03-15 14/700 2026-03-19 22:59 by ѧԱ8dgXkO
[¿¼²©] É격26Äê +3 °Ë6°Ë68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[¿¼ÑÐ] ¸´ÊÔµ÷¼Á +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[¿¼ÑÐ] Ò»Ö¾Ô¸Î÷ÄϽ»´ó£¬Çóµ÷¼Á +4 ²Ä»¯ÖðÃÎÈË 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[¿¼ÑÐ] ÊÕ¸´ÊÔµ÷¼ÁÉú +4 ÓêºóÇïºÉ 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[¿¼ÑÐ] 302Çóµ÷¼Á +10 ºôºôºô¡£¡£¡£¡£ 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[¿¼ÑÐ] 280Çóµ÷¼Á +6 ¹¾ààÏþÏþ 2026-03-18 7/350 2026-03-18 11:25 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 303Çóµ÷¼Á +4 î£08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[¿¼ÑÐ] 334Çóµ÷¼Á +3 Ö¾´æ¸ßÔ¶ÒâÔÚ»úÐ 2026-03-16 3/150 2026-03-18 08:34 by lm4875102
[¿¼ÑÐ] 293Çóµ÷¼Á +11 zjlµÄºÅ 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[¿¼ÑÐ] 268Çóµ÷¼Á +8 Ò»¶¨ÓÐѧÉÏ- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[˶²©¼ÒÔ°] ºþ±±¹¤Òµ´óѧ ÉúÃü¿ÆÑ§Ó뽡¿µÑ§Ôº-¿ÎÌâ×éÕÐÊÕ2026¼¶Ê³Æ·/ÉúÎï·½Ïò˶ʿ +3 1ϲ´º8 2026-03-17 5/250 2026-03-17 17:18 by ber´¨cool×Ó
[¿¼ÑÐ] ÓÐûÓеÀÌú/ÍÁľµÄÏëµ÷¼ÁÄÏÁÖ£¬¸ø×Ô¼ºÕÐʦµÜÖС« +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[¿¼ÑÐ] Ò»Ö¾Ô¸ÄϾ©´óѧ£¬080500²ÄÁÏ¿ÆÑ§Ó빤³Ì£¬µ÷¼Á +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[¿¼ÑÐ] 283Çóµ÷¼Á +3 Ìý·ç¾ÍÊÇÓꣻ 2026-03-16 3/150 2026-03-17 07:41 by ÈÈÇéɳĮ
[»ù½ðÉêÇë] ½ñÄêµÄ¹ú»ù½ðÊÇ´ò·ÖÖÆÂ𣿠50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by ±±¾©À³ÒðÈóÉ«
[¿¼ÑÐ] 297Çóµ÷¼Á +4 ѧº£Æ¯²´ 2026-03-13 4/200 2026-03-14 11:51 by ÈÈÇéɳĮ
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û