24小时热门版块排行榜    

查看: 5667  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +9 majunge000 2026-08-11 11/550 2026-08-15 23:20 by botar
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
信息提示
请填处理意见