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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf-5.0得到的声子谱很乱是怎么回事？
91/1返回列表
查看: 3255  |  回复: 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 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我不理解！ +19 Edward_pc 2026-08-26 27/1350 2026-09-02 12:08 by Ai19911123
[基金申请] 要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 16/800 2026-09-02 12:07 by Ai19911123
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 10/500 2026-09-02 11:35 by 大不刘6
[基金申请] 科研人应该花精力去思考如何解决问题，而不是去凝练问题 +4 瞬息宇宙 2026-09-01 8/400 2026-09-02 10:35 by 瞬息宇宙
[基金申请] 学科评审组评审是指会评吗？ +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[文学芳草园] 梦想 +5 myrtle 2026-08-26 8/400 2026-09-02 02:42 by Leogzhya
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 面上函评意见出来了，像什么等级？ 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了，不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉？ +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 基金不中，共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知？ +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
信息提示
关闭
请填处理意见
关闭