| 查看: 2783 | 回复: 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
» 猜你喜欢
一志愿南京航空航天大学 材料与化工329分求调剂
已经有7人回复
-
一志愿C9的化学工程（085602） 340分，感觉校内调剂无望，求调剂
已经有12人回复
-
312求调剂
已经有4人回复
-
材料专硕322
已经有9人回复
-
求调剂
已经有7人回复
-
085600材料与化工301分求调剂院校
已经有8人回复
-
一志愿 江南大学 085602 化工专硕 338分求调剂
已经有14人回复
-
070300化学学硕311分求调剂
已经有13人回复
-
材料调剂
已经有9人回复
-
一志愿河北工业大学材料工程，初试344求专硕调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用pwscf计算锂的声子色散谱
已经有11人回复
- PWSCF在哪里设置lapack？再问一句声子谱对角化出错是lapack的问题吗？
已经有4人回复
- 是否有软件可实现沿PWscf计算得到的某一声子振动方向来调节晶格中的原子位置
已经有5人回复
- 【求助】pwscf的声子态密度图和色散曲线图能用origin画出来吗？
已经有5人回复
- 【求助】关于PWSCF介电常数和声子谱的计算
已经有7人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +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 ：电子结构
weiyongkai3楼2013-05-08 16:32:13
weiyongkai
木虫 (职业作家)
- 应助: 13 (小学生)
- 金币: 3192
- 散金: 209
- 红花: 3
- 帖子: 3186
- 在线: 228.1小时
- 虫号: 1104094
- 注册: 2010-09-20
- 专业: 凝聚态物性I:结构、力学和
4楼2013-05-08 17:37:10
chuanghua304
禁虫 (职业作家)
|
本帖内容被屏蔽 |
5楼2013-06-29 21:34:14
6楼2013-08-23 23:44:50
chuanghua304
禁虫 (职业作家)
|
本帖内容被屏蔽 |
7楼2013-08-24 09:00:12
111hzs
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 132.8
- 帖子: 63
- 在线: 44.8小时
- 虫号: 1951933
- 注册: 2012-08-22
- 专业: 凝聚态物性 II ：电子结构
|
你好，我在用QE5.0计算的声子和QE4.3.2计算得到的结果不一样。 用QE5.0，计算声子，q网格取4*4*1时，在ph.out里面发现 Dynamical matrices for ( 4, 4, 1) uniform grid of q-points ( 4q-points): N xq(1) xq(2) xq(3) 1 0.000000000 0.000000000 0.000000000 2 0.000000000 0.288675135 0.000000000 3 0.000000000 -0.577350269 0.000000000 4 0.250000000 -0.721687836 0.000000000 而用4.3.2时，里面显示： Dynamical matrices for ( 4, 4, 1,) uniform grid of q-points ( 7q-points): N xq(1) xq(2) xq(3) 1 0.000000000 0.000000000 0.000000000 2 0.000000000 0.288675135 0.000000000 3 0.000000000 -0.577350269 0.000000000 4 0.250000000 -0.144337567 0.000000000 5 0.250000000 -0.721687836 0.000000000 6 -0.500000000 0.288675135 0.000000000 7 -0.500000000 0.577350269 0.000000000 为什么有这样的差别呢？？ |
8楼2013-12-17 21:41:56
9楼2018-05-06 16:33:21
