²é¿´: 5158  |  »Ø¸´: 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 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 26µ÷¼Á/²ÄÁÏ/Ó¢Ò»Êý¶þ/×Ü·Ö289/ÒѹýAÇøÏß +8 ²½´¨¿á×Ï123 2026-03-13 8/400 2026-03-20 08:47 by xingguangj
[¿¼ÑÐ] Ò»Ö¾Ô¸ËÕÖÝ´óѧ²ÄÁÏÇóµ÷¼Á£¬×Ü·Ö315£¨Ó¢Ò»£© +3 sbdksD 2026-03-19 3/150 2026-03-19 23:21 by fmesaito
[¿¼ÑÐ] 320Çóµ÷¼Á0856 +3 ²»ÏëÆðÃû×Ö112 2026-03-19 3/150 2026-03-19 22:53 by ѧԱ8dgXkO
[¿¼ÑÐ] Ò»Ö¾Ô¸ÄϾ©Àí¹¤´óѧ085701×ÊÔ´Óë»·¾³302·ÖÇóµ÷¼Á +3 ¿ûè÷ÎÀ¶Ó 2026-03-18 5/250 2026-03-19 19:35 by ¸øÄãÄã×¢ÒâÐÝÏ¢
[¿¼ÑÐ] 286Çóµ÷¼Á +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[¿¼ÑÐ] 328Çóµ÷¼Á£¬Ó¢ÓïÁù¼¶551£¬ÓпÆÑо­Àú +4 ÉúÎ﹤³Ìµ÷¼Á 2026-03-16 12/600 2026-03-19 11:10 by ÉúÎ﹤³Ìµ÷¼Á
[¿¼ÑÐ] ±¾¿ÆÖ£ÖÝ´óѧÎïÀíѧԺ£¬Ò»Ö¾Ô¸»ª¿Æ070200ѧ˶£¬346Çóµ÷¼Á +4 ÎÒ²»ÊÇÒ»¸ù´Ð 2026-03-18 4/200 2026-03-19 09:11 by ¸¡ÔÆ166
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤ +5 °²È«Éϰ¶£¡ 2026-03-16 5/250 2026-03-18 15:33 by cmz0325
[¿¼ÑÐ] 297Çóµ÷¼Á +8 Ï·¾«µ¤µ¤µ¤ 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[¿¼ÑÐ] 331Çóµ÷¼Á£¨0703Óлú»¯Ñ§ +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[¿¼ÑÐ] 312Çóµ÷¼Á +8 İå·Ï£ 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[¿¼ÑÐ] 0703»¯Ñ§Çóµ÷¼Á ×Ü·Ö331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á +3 ÄÝÄÝninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[¿¼ÑÐ] ²ÄÁÏ£¬·ÄÖ¯£¬ÉúÎ0856¡¢0710£©£¬»¯Ñ§ÕÐÉúÀ² +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[»ù½ðÉêÇë] ±»ÎÒÑÔÖУºÐÂÄ£°å²»Ç¿µ÷¸ñʽÁË£¬¼Ùר¼Ò¿ªÊ¼¹Ü¸ñʽÁË +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by »ÆÄñÓÚ·ÉChao
[¿¼ÑÐ] ¿¼ÑÐÇóµ÷¼Á +3 éÙËÌ. 2026-03-17 4/200 2026-03-17 21:43 by ÓÐÖ»ÀêÅ«
[¿¼ÑÐ] Ò»Ö¾Ô¸ËÕÖÝ´óѧ²ÄÁϹ¤³Ì£¨085601£©×¨Ë¶ÓпÆÑо­ÀúÈýÏî¹ú½±Á½¸öʵÓÃÐÍרÀûÒ»ÏîÊ¡¼¶Á¢Ïî +6 ´ó»ðɽС»ðɽ 2026-03-16 8/400 2026-03-17 15:05 by ÎÞи¿É»÷111
[¿¼ÑÐ] 0854¿ØÖƹ¤³Ì 359Çóµ÷¼Á ¿É¿çרҵ +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[¿¼ÑÐ] 070300»¯Ñ§Ñ§Ë¶Çóµ÷¼Á +6 Ì«Ïë½ø²½ÁË0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[¿¼ÑÐ] 0703 ÎïÀí»¯Ñ§µ÷¼Á +3 ÎÒ¿ÉÒÔÉϰ¶µÄ¶Ô 2026-03-13 5/250 2026-03-16 10:50 by ÎÒ¿ÉÒÔÉϰ¶µÄ¶ÔÂ
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û