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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf-5.0得到的声子谱很乱是怎么回事？
91/1返回列表
查看: 3230  |  回复: 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一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 8/400 2026-08-25 05:12 by BZKMTicpDhFj
[基金申请] 明天应该可查了！？ +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] filecode，4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤？
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭