24小时热门版块排行榜    

查看: 4881  |  回复: 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 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见