应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf-5.0得到的声子谱很乱是怎么回事？
【调剂】北京石油化工学院2024年16个专业接受调剂
91/1返回列表
查看: 2203  |  回复: 8
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com
  • 附件 1 : ps.rar
    • 2013-05-08 11:07:29, 148.32 K

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
凌云
1楼2013-05-08 11:07:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

DFPT

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +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中的文件名修改得跟之间色散计算的不同就可以了。
赞一下(3人) 回复此楼
010-8393-9164
2楼2013-05-08 12:10:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
引用回帖:
2楼: Originally posted by DFPT at 2013-05-08 12:10:25
楼主在使用matdyn.x计算特殊路径的振动谱后，又使用同一个名字Al444.freq计算了PhDOS，所以Al444.freq中存储的是PhDOS中定义的nk1*nk2*nk3的q点网格，自然从Al444.freq中提取出来的色散谱就是如楼主图中所示的了。其 ...

额，真仔细，我还没看到这个错误呢。确实重名了。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

weiyongkai
3楼2013-05-08 16:32:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiyongkai

木虫 (职业作家)
送红花一朵
引用回帖:
3楼: Originally posted by gemucai at 2013-05-08 16:32:13
额，真仔细，我还没看到这个错误呢。确实重名了。...

谢谢谢谢！
回复此楼
凌云
4楼2013-05-08 17:37:10
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuanghua304

禁虫 (职业作家)
本帖内容被屏蔽
5楼2013-06-29 21:34:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsdwrk

新虫 (小有名气)
引用回帖:
5楼: Originally posted by chuanghua304 at 2013-06-29 21:34:14
请问计算声子色散的高对称点，为什么没有了关键字"K_POINTS"，那又怎么设置后面的tpiba（默认）， crystal ，tpiba_b 等等，难道只能用默认的嘛？

应该指的是分数坐标！
赞一下 回复此楼
6楼2013-08-23 23:44:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuanghua304

禁虫 (职业作家)
本帖内容被屏蔽
7楼2013-08-24 09:00:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

111hzs

铜虫 (小有名气)
你好，我在用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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

AryaLu

新虫 (初入文坛)
请问问题解决了么
赞一下 回复此楼
9楼2018-05-06 16:33:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 weiyongkai 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 哎 +5 s150535912 2024-04-18 6/300 2024-04-19 20:33 by s150535912
[论文投稿] 一审一个审稿人，小修，会怎么样呀？ +9 林师妹 2024-04-18 9/450 2024-04-19 20:00 by nono2009
[论文投稿] 最近遇到这样一个问题 2+4 asd123gfa689 2024-04-18 10/500 2024-04-19 18:57 by asd123gfa689
[育儿交流] 小孩7岁，上一年级，不自信，学习较差 10+7 ljaxida 2024-04-16 34/1700 2024-04-19 18:49 by 水可南阳
[基金申请] 国自然青年基金只有一篇能第一标注，其他都是第二标注能结题吗？ 10+4 lice_1987 2024-04-18 6/300 2024-04-19 08:29 by kudofaye
[论文投稿] The total environment格式 25+4 lyyyyds0 2024-04-13 9/450 2024-04-18 09:33 by 我为人人研
[无机/物化] H3PO2在130度附近分解为磷化氢（剧毒），怎样防护，求指点。 1(金币+5) +3 大力2010 2024-04-13 9/450 2024-04-17 21:55 by 大力2010
[考研] 工学322，英语六级求调剂 +3 domax 2024-04-17 10/500 2024-04-17 21:25 by 东方花园小万
[公派出国] CSC在职申请求问 +5 Amber0919 2024-04-16 5/250 2024-04-17 16:15 by liaojiayan
[考研] 广州大学光电信息工程专业调剂，招收物理学专业学生 +5 txhx4010 2024-04-14 7/350 2024-04-16 10:52 by domax
[考博] 博士申请 +5 啦啦啦3737 2024-04-13 9/450 2024-04-16 07:07 by 1147wzzzz
[考研] 材料专硕329调剂遗留难民 +9 Kaylawander 2024-04-13 9/450 2024-04-15 19:21 by mumin1990
[考研] 273求调剂 +5 Late婉安 2024-04-15 7/350 2024-04-15 13:01 by Late婉安
[考研] 278求调剂 +4 月亮就蒜 2024-04-13 5/250 2024-04-14 23:03 by 永字号
[考研] 300求调剂 +5 @Xqy668800 2024-04-14 6/300 2024-04-14 23:03 by lincunhui
[考研] 一志愿厦门大学329分材料工程专硕求调剂 +10 Kaylawander 2024-04-13 11/550 2024-04-14 22:39 by liu823948201
[考研] 290求调剂 +3 杨yhr 2024-04-14 5/250 2024-04-14 21:50 by coco1981
[考研] 334求调剂 +4 学药救人 2024-04-13 6/300 2024-04-13 20:27 by 献世的王
[考研] 086000生医没考数学304分，一志愿西农求调剂 +3 王有田有学上 2024-04-13 5/250 2024-04-13 16:18 by coco1981
[考博] 现在申博很晚了吗 +3 熙风慕寒 2024-04-13 3/150 2024-04-13 11:08 by Nanamiwww
信息提示
关闭
请填处理意见
关闭