| 查看: 5644 | 回复: 23 | ||||||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||||||||
[交流]
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
- 附件 2 : phfreq.exe
2013-03-04 20:04:39, 45.86 K
2013-03-05 16:37:30, 349 K
» 收录本帖的淘帖专辑推荐
 精华网帖收集 | 科研技巧 | 第一性原理计算辅助工具 | condensed matter physics |
| 计算模拟 | 第一性原理计算 | 计算和软件 |
» 本帖已获得的红花（最新10朵）
renren123123» 猜你喜欢
面上项目filecode邪修
已经有4人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有13人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 我的PWSCF一直报错，貌似是输入问题，但不知错在哪，新手求指导。
已经有5人回复
- 请大牛解答一个概念性问题，表观活化能与绝对能垒，有实例分析
已经有9人回复
- 以石墨烯为例Pwscf计算态密度的详细步骤和要点
已经有11人回复
- 虫友结构解析修正问题解答实例
已经有70人回复
hakuna
木虫 (知名作家)
- 1ST强帖: 3
- 应助: 82 (初中生)
- 金币: 8735.9
- 散金: 158
- 红花: 126
- 帖子: 5599
- 在线: 1367.2小时
- 虫号: 167622
- 注册: 2006-01-15
- 专业: 考古理论
3楼2013-03-05 06:49:10
0.8
| 顶一下，感谢分享！ |
2楼2013-03-04 21:23:36
0.8
| 谢谢分享！！！！！！！！！ |
4楼2013-03-05 06:56:35
0.8
| 顶一下，感谢分享！ |
5楼2013-03-05 09:15:32
20