±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 5216  |  »Ø¸´: 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µÄ»ØÌû
»ØÌûÖö¥ ( ¹²ÓÐ1¸ö )

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 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ²ÄÁϵ÷¼Á +12 Ò»ÑùYWY 2026-04-02 12/600 2026-04-03 09:46 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] ÉúÎïѧ308·ÖÇóµ÷¼Á£¨Ò»Ö¾Ô¸»ª¶«Ê¦´ó£© +6 ÏàÐűػá¹ââÍòÕ 2026-03-31 7/350 2026-04-02 23:16 by JourneyLucky
[¿¼ÑÐ] 326·ÖÇóµ÷¼Á +3 ÓÚÊǺõÄØ 2026-04-01 3/150 2026-04-02 21:30 by jing_jenna
[¿¼ÑÐ] 318Çóµ÷¼Á +3 óÆÐÐÖÂÔ¶. 2026-03-31 4/200 2026-04-02 15:56 by Jaylen.
[¿¼ÑÐ] ²ÄÁÏ»¯¹¤340Çóµ÷¼Á +5 jhx777 2026-03-30 5/250 2026-04-02 12:45 by smileboy2006
[¿¼ÑÐ] 271Çóµ÷¼Á +15 ÀÕ²¼ÀÊ@ 2026-03-31 20/1000 2026-04-02 11:24 by Sammy2
[¿¼ÑÐ] 367Çóµ÷¼Á +8 ÓóÄàà£à£¡­ 2026-03-28 8/400 2026-04-02 09:20 by olim
[¿¼ÑÐ] 348»·¾³¹¤³Ìµ÷¼Á +3 ÎâÑå׿24k 2026-04-01 3/150 2026-04-02 09:14 by nanaliuyun
[¿¼ÑÐ] 303Çóµ÷¼Á +8 DLkz1314. 2026-03-30 8/400 2026-04-02 09:08 by sanrepian
[¿¼ÑÐ] »¯Ñ§¹¤³Ìר˶324·Ö£¬Ò»Ö¾Ô¸Öйú¿óÒµ´óѧÇóµ÷¼Á +7 ¹¢¹¢1314 2026-04-01 7/350 2026-04-02 07:40 by ÉÐË®¸óÖ÷
[¿¼ÑÐ] Ó¢Ò»ÊýÒ»408£¬×Ü·Ö284£¬¶þÕ½Õæ³ÏÇóµ÷¼Á +12 12.27 2026-03-30 14/700 2026-04-02 00:18 by ÐÀϲ777
[¿¼ÑÐ] Çóµ÷¼Á0703 +5 ÖܼÎÒ¢ 2026-03-31 8/400 2026-04-01 20:32 by ltltkkk
[¿¼ÑÐ] ÉúÎïÓëÒ½Ò©¿¼Ñе÷¼Á +5 Ìúº©º©123425 2026-03-31 5/250 2026-04-01 18:01 by syh9288
[¿¼ÑÐ] 350Çóµ÷¼Á +7 °¢¼Ñ¡« 2026-03-31 7/350 2026-04-01 16:12 by yanflower7133
[¿¼ÑÐ] »·¾³¹¤³Ìµ÷¼Á +9 hyzzzzzzz. 2026-04-01 9/450 2026-04-01 14:20 by salamander`
[¿¼ÑÐ] 085600£¬321·ÖÇóµ÷¼Á +13 ´ó²öС×Ó 2026-03-31 13/650 2026-04-01 12:35 by chemdavid
[¿¼ÑÐ] 086000ÉúÎïÓëÒ½Ò©298µ÷¼ÁÇóÖú +4 ÔªÔªÇàÇà 2026-03-31 6/300 2026-04-01 11:13 by syh9288
[¿¼ÑÐ] ÎïÀíѧµ÷¼Á +4 СÑò36 2026-03-30 4/200 2026-03-31 16:16 by lishahe
[¿¼ÑÐ] 277¹òÇóµ÷¼Á +8 1915668 2026-03-27 13/650 2026-03-31 14:58 by ÍõÁÁ_´óÁ¬Ò½¿Æ´ó
[¿¼ÑÐ] 348Çóµ÷¼Á +6 СÀÁ³æ²»ÀÁÁË 2026-03-28 6/300 2026-03-30 10:29 by Evan_Liu
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û