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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf-5.0得到的声子谱很乱是怎么回事？
91/1返回列表
查看: 3249  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-31 03:18 by hZiFeudZyoGR
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 zICmwzsBXjbN 2026-08-29 6/300 2026-08-31 03:09 by hZiFeudZyoGR
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 6/300 2026-08-31 02:45 by hZiFeudZyoGR
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-31 02:34 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-31 02:34 by hZiFeudZyoGR
[找工作] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-31 01:59 by hZiFeudZyoGR
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +6 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-31 01:16 by hZiFeudZyoGR
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 22:24 by xiangy672
[基金申请] 为什么到现在没收到通知？ +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
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 5/250 2026-08-30 20:41 by hZiFeudZyoGR
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:21 by hZiFeudZyoGR
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 基金未中，这种答复是模板吗？ +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
关闭
请填处理意见
关闭