24小时热门版块排行榜    

查看: 5647  |  回复: 23

lu_wei07

银虫 (小有名气)

[交流] PWSCF 具体实例 详细解答各种问题 已有13人参与

PWSCF 具体实例 详细解答各种问题,摘自网络一个blogblog 地址
http://zhaoyunawuli.blog.163.com ... 976200811147143279/
网易的  网名: 随梦想而飞
1,附 处理PWSCF中计算得到本征频率的程序   自己修改了
用的时候需要把*.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
同学,请注重知识产权,你都说了是博客里面的,那就把博主贴出来。没看我就猜测是侯柱峰老师的,打开看果然是。而且你分享的这种方式着实不好,你好歹说一些概括在贴子里面。贴出一部分原文。。。其实这个东西只要大家连接的了网络就都可以在网页上直接看的。
无星级,不是评论资源。而是给你这个帖子。
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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于代码变化问题,想知道的进来 +15 且听虎啸 2026-08-07 22/1100 2026-08-09 22:13 by 布布和一二
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:19 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[基金申请] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 14:27 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
信息提示
请填处理意见