24小时热门版块排行榜    

查看: 5156  |  回复: 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的回帖

0.8

谢谢分享!!!!!!!!!
4楼2013-03-05 06:56:35
已阅   回复此楼   关注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

顶一下,感谢分享!
5楼2013-03-05 09:15:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +4 @taotao 2026-03-20 4/200 2026-03-20 10:37 by EBSD
[考研] 08工学调剂 +4 用户573181 2026-03-20 4/200 2026-03-20 10:32 by 朗月清风1
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考研] 279分求调剂 一志愿211 +7 chaojifeixia 2026-03-19 7/350 2026-03-20 09:00 by sunny300
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +8 步川酷紫123 2026-03-13 8/400 2026-03-20 08:47 by xingguangj
[考研] 281求调剂(0805) +11 烟汐忆海 2026-03-16 21/1050 2026-03-20 08:45 by 无际的草原
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 317求调剂 +3 申子申申 2026-03-19 6/300 2026-03-19 14:16 by 申子申申
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考研] 278求调剂 +5 烟火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
信息提示
请填处理意见