| 查看: 3238 | 回复: 8 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
weiyongkai木虫 (职业作家)
|
[求助]
pwscf-5.0得到的声子谱很乱是怎么回事？
|
|
|
最近在用QE-5.0计算声子谱，得到的图形非常乱，大家看看是什么原因？输入文件如下： #! /bin/sh # # SCF at dense k-mesh, good enough for electronic DOS # cat > al.scf.in <<! &control title='al' calculation = 'scf' restart_mode='from_scratch', prefix='al', pseudo_dir = './ps', outdir='./tmp' tprnfor=.true. / &system ibrav=4, celldm(1)=6.05606, celldm(3)=1.71298, nat=2, ntyp=1, nbnd= 30, ecutwfc=30.0, occupations ='smearing', degauss =0.01 smearing ='mp' / &electrons diagonalization='cg' diago_cg_maxiter= 60 mixing_mode = 'plain' mixing_beta = 0.5 conv_thr = 1.0d-6 / ATOMIC_SPECIES Sc 44.955910 Sc.pz-hgh.UPF ATOMIC_POSITIONS (crystal) Sc 0.3333333333333286 0.6666666666666714 0.2500000000000000 Sc 0.6666666666666714 0.3333333333333286 0.7500000000000000 K_POINTS (automatic) 8 8 6 0 0 0 ! mpirun -np 12 pw.x<al.scf.in>al.scf.out cat > al.elph.in << ! phonon for Sc &inputph tr2_ph=1.0d-10, prefix='al', fildvscf='aldv', amass(1)=44.955910, outdir='./tmp', fildyn='al.dyn', ldisp=.true., trans=.true., nq1=4, nq2=4, nq3=2 / ! mpirun -np 12 ph.x<al.elph.in> al.elph.out # q2r and matdyn # cat > q2r.in << ! &input zasr='simple', fildyn='al.dyn', flfrc='Al444.fc', la2F=.false., / ! mpirun -np 12 q2r.x < q2r.in > q2r.out cat > matdyn.in.freq << ! &input asr='simple', amass(1)=44.955910, flfrc='Al444.fc', flfrq='Al444.freq', la2F=.false., dos=.false. / 121 0.000000 0.000000 0.000000 1.00 0.033333 0.000000 0.000000 1.00 0.066667 0.000000 0.000000 1.00 0.100000 0.000000 0.000000 1.00 0.133333 0.000000 0.000000 1.00 0.166667 0.000000 0.000000 1.00 0.200000 0.000000 0.000000 1.00 ............... ! mpirun -np 12 matdyn.x <matdyn.in.freq> matdyn.out.freq # # # cat > matdyn.in.dos <<! / &input asr='simple', amass(1)=44.955910, flfrc='Al444.fc', flfrq='Al444.freq', la2F=.false., dos=.true. fldos='phonon.dos', nk1=10, nk2=10, nk3=10, ndos=50 / ! mpirun -np 12 matdyn.x <matdyn.in.dos> matdyn.out.dos  无标题.png [ Last edited by weiyongkai on 2013-5-8 at 11:10 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : ps.rar
2013-05-08 11:07:29, 148.32 K
» 猜你喜欢
某些机构，以效率低为荣，以效率低作为存在感
已经有4人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有3人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
-
没有任何消息-是不是就凉了
已经有10人回复
-
人气不行了
已经有11人回复
-
2026国自然函评费到账
已经有22人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有3人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有7人回复
-
今天基金会出结果吗？20260819
已经有18人回复
-
我面上完蛋了
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用pwscf计算锂的声子色散谱
已经有11人回复
- PWSCF在哪里设置lapack？再问一句声子谱对角化出错是lapack的问题吗？
已经有4人回复
- 是否有软件可实现沿PWscf计算得到的某一声子振动方向来调节晶格中的原子位置
已经有5人回复
- 【求助】pwscf的声子态密度图和色散曲线图能用origin画出来吗？
已经有5人回复
- 【求助】关于PWSCF介电常数和声子谱的计算
已经有7人回复
weiyongkai
木虫 (职业作家)
- 应助: 13 (小学生)
- 金币: 3192
- 散金: 209
- 红花: 3
- 帖子: 3186
- 在线: 228.1小时
- 虫号: 1104094
- 注册: 2010-09-20
- 专业: 凝聚态物性I:结构、力学和
送红花一朵 4楼2013-05-08 17:37:10
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
weiyongkai: 金币+50, ★★★很有帮助, 谢谢！ 2013-05-08 17:36:45
liliangfang: 金币+2, 谢谢交流 2013-05-11 09:24:11
感谢参与，应助指数 +1
weiyongkai: 金币+50, ★★★很有帮助, 谢谢！ 2013-05-08 17:36:45
liliangfang: 金币+2, 谢谢交流 2013-05-11 09:24:11
| 楼主在使用matdyn.x计算特殊路径的振动谱后，又使用同一个名字Al444.freq计算了PhDOS，所以Al444.freq中存储的是PhDOS中定义的nk1*nk2*nk3的q点网格，自然从Al444.freq中提取出来的色散谱就是如楼主图中所示的了。其实，只需要把PhDOS计算中的flfrq中的文件名修改得跟之间色散计算的不同就可以了。 |
2楼2013-05-08 12:10:25
gemucai
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 1774.6
- 散金: 157
- 红花: 4
- 帖子: 693
- 在线: 201.7小时
- 虫号: 816478
- 注册: 2009-07-27
- 专业: 凝聚态物性 II ：电子结构
3楼2013-05-08 16:32:13
chuanghua304
禁虫 (职业作家)
|
本帖内容被屏蔽 |
5楼2013-06-29 21:34:14
20