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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf-5.0得到的声子谱很乱是怎么回事？
91/1返回列表
查看: 2809  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 297，工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 一志愿沪9，326求生物学调剂 +12 刘墨墨 2026-04-13 12/600 2026-04-18 23:31 by 路病情
[考研] 320求调剂 +5 深郊akm 2026-04-17 5/250 2026-04-18 19:52 by 王珺璞
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 297，工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 一志愿中科大材料与化工，353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[考研] 295分求调剂 +5 ?要上岸? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[考研] 一志愿沪9，生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 297，工科调剂?河南农业大学本科 +14 河南农业大学-能 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[考研] 一志愿A区211，22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
信息提示
关闭
请填处理意见
关闭