²é¿´: 5159  |  »Ø¸´: 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

¶¥Ò»Ï£¬¸Ðл·ÖÏí£¡
9Â¥2013-03-05 16:02:02
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢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µÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 286·ÖÈ˹¤ÖÇÄÜרҵÇëÇóµ÷¼ÁÔ¸Òâ¿ç¿¼£¡ +3 lemonzzn 2026-03-17 4/200 2026-03-20 11:04 by lemonzzn
[¿¼ÑÐ] ²ÄÁÏѧ˶297ÒѹýËÄÁù¼¶Çóµ÷¼ÁÍÆ¼ö +5 adaie 2026-03-19 5/250 2026-03-20 10:44 by EBSD
[¿¼ÑÐ] Ò»Ö¾Ô¸Î人Àí¹¤²ÄÁϹ¤³Ìר˶µ÷¼Á +7 Doleres 2026-03-19 7/350 2026-03-20 10:39 by EBSD
[¿¼ÑÐ] 304Çóµ÷¼Á +5 ÂüÊâ2266 2026-03-18 5/250 2026-03-20 09:00 by ZHANG0tao
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤ Çóµ÷¼Á +17 enenenhui 2026-03-13 18/900 2026-03-20 08:55 by xingguangj
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á +10 ÄÝÄÝninicgb 2026-03-15 14/700 2026-03-19 22:59 by ѧԱ8dgXkO
[¿¼ÑÐ] 320Çóµ÷¼Á0856 +3 ²»ÏëÆðÃû×Ö112 2026-03-19 3/150 2026-03-19 22:53 by ѧԱ8dgXkO
[¿¼²©] ¶«»ªÀí¹¤´óѧ»¯²Äרҵ26½ì˶ʿ²©Ê¿ÉêÇë +8 zlingli 2026-03-13 8/400 2026-03-19 16:32 by ÇáËɲ»ÉÙËæ
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[¿¼ÑÐ] 346Çóµ÷¼Á[0856] +3 WayneLim327 2026-03-16 6/300 2026-03-19 11:21 by WayneLim327
[¿¼ÑÐ] 330Çóµ÷¼Á +3 С²Ä»¯±¾¿Æ 2026-03-18 3/150 2026-03-18 21:55 by ÎÞи¿É»÷111
[¿¼ÑÐ] 334Çóµ÷¼Á +3 Ö¾´æ¸ßÔ¶ÒâÔÚ»úÐ 2026-03-16 3/150 2026-03-18 08:34 by lm4875102
[¿¼ÑÐ] 277µ÷¼Á +5 ×ÔÓɼå±ý¹û×Ó 2026-03-16 6/300 2026-03-17 19:26 by Àîleezz
[˶²©¼ÒÔ°] ºþ±±¹¤Òµ´óѧ ÉúÃü¿ÆÑ§Ó뽡¿µÑ§Ôº-¿ÎÌâ×éÕÐÊÕ2026¼¶Ê³Æ·/ÉúÎï·½Ïò˶ʿ +3 1ϲ´º8 2026-03-17 5/250 2026-03-17 17:18 by ber´¨cool×Ó
[¿¼ÑÐ] Ò»Ö¾Ô¸ËÕÖÝ´óѧ²ÄÁϹ¤³Ì£¨085601£©×¨Ë¶ÓпÆÑо­ÀúÈýÏî¹ú½±Á½¸öʵÓÃÐÍרÀûÒ»ÏîÊ¡¼¶Á¢Ïî +6 ´ó»ðɽС»ðɽ 2026-03-16 8/400 2026-03-17 15:05 by ÎÞи¿É»÷111
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤×¨Ë¶µ÷¼Á +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by Ó¸ÒÌ«¼àÍõ¹«¹«
[¿¼ÑÐ] »úеר˶325£¬Ñ°ÕÒµ÷¼ÁԺУ +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[¿¼ÑÐ] 070300»¯Ñ§Ñ§Ë¶Çóµ÷¼Á +6 Ì«Ïë½ø²½ÁË0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[¿¼ÑÐ] ÖпÆÔº²ÄÁÏ273Çóµ÷¼Á +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[¿¼ÑÐ] 0856ר˶279Çóµ÷¼Á +5 ¼ÓÓͼÓÓÍ£¡? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û