24小时热门版块排行榜    

Znn3bq.jpeg
查看: 6048  |  回复: 29

坏逊平儿

新虫 (初入文坛)

[求助] 硅烯声子谱G点附近虚频 已有4人参与

请教各位大神,小弟入手pwscf不久,计算硅烯的声子谱,可发现两个问题实在难以解决:1、G点附近总有虚频,该怎么消除?2、在一个高对称点处两条光学支本该简并的,但是我计算出来的结果却没有,是什么原因?如下图

(硅烯的结构是组里已经优化好的结构,算声子谱直接做自洽计算吧?!)
跪求大神指导,万分感谢。

硅烯声子谱G点附近虚频
搜狗截图20160802150307.png
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

liqizuiyang

木虫 (著名写手)

【答案】应助回帖

★ ★
蛋蛋小童鞋: 金币+2, 应助指数+1, 鼓励回答问题详细者 2016-08-03 14:34:12
重新优化结构,力收敛判据调至0.0001eV/A。

改用MP型展宽,做严格收敛测试。

K点改成18*18*1.

diagonalization改成david。cg只在空带数很大的时候才有优势,速度上并不占优。
6楼2016-08-03 09:37:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
蛋蛋小童鞋: 金币+1, 谢谢参与 2016-08-11 08:37:38
漫天飘雪: 回帖置顶 2018-06-24 13:12:01
按照如下文章里的参数,肯定能得到好的声子谱:
http://journals.aps.org/prb/abstract/10.1103/PhysRevB.91.205433

尤其你的conv_thr还太大,建议调到d-12。nq应该还太小,一般在8以上。
先进使役材料计算(https://aidme.nimte.ac.cn)
14楼2016-08-04 00:21:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

【答案】应助回帖


漫天飘雪: 回帖置顶 2018-06-24 13:12:05
当你nq是3的倍数或足够密的的时候,你自然会得到简并了。

» 本帖已获得的红花(最新10朵)

先进使役材料计算(https://aidme.nimte.ac.cn)
16楼2016-08-04 00:26:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c19940408

银虫 (初入文坛)

【答案】应助回帖


蛋蛋小童鞋: 金币+1, 谢谢回答 2016-08-11 08:37:15
建议使用和优化时相近的参数来计算。。比如相同的能量截断,k点等等
20楼2016-08-08 02:14:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

坏逊平儿

新虫 (初入文坛)

scf.in参数设置


&control
calculation='scf',
restart_mode='from_scratch',
prefix='si'
pseudo_dir = '/home/bxu/zjp/2sixi/4/si/',
outdir='/home/bxu/zjp/2sixi/4/si/tmp/'
/
&system
ibrav = 4, celldm(1) =7.3052690583, celldm(3) =3.0593725672, nat= 2, ntyp= 1,
ecutwfc = 50,ecutrho=450
occupations='smearing',degauss=0.02,smearing='gaussian'
/
&electrons
electron_maxstep=200
mixing_beta = 0.5
conv_thr = 1.0d-10
diagonalization='cg'
/
ATOMIC_SPECIES
Si 28.086 Si.pbe-n-rrkjus_psl.1.0.0.UPF
ATOMIC_POSITIONS {crystal}
Si 0.333330005 0.666660011 0.480724089
Si 0.666660011 0.333330005 0.519265927
K_POINTS {automatic}
20 20 1 0 0 0
2楼2016-08-02 15:31:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

坏逊平儿

新虫 (初入文坛)

ph.in参数设置


&inputph
tr2_ph=1.0d-16,
prefix='si',
amass(1)=28.086,
outdir='/home/bxu/zjp/2sixi/4/si/tmp',
fildyn='si.dyn',
elph=.false.,
trans=.true.,
ldisp=.true.
nq1=4, nq2=4,nq3=1
/
3楼2016-08-02 15:32:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

坏逊平儿

新虫 (初入文坛)

q2r.in 参数设置


&input
zasr='simple',fildyn='si.dyn',flfrc='si.fc',la2F=.false.
4楼2016-08-02 15:32:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

golddoushi

木虫 (正式写手)

【答案】应助回帖

蛋蛋小童鞋: 应助指数+1, 谢谢参与 2016-08-03 08:42:11
gamma 点本来就有虚频,对应长波极限。一般直接忽略不计
5楼2016-08-03 07:04:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

坏逊平儿

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by golddoushi at 2016-08-03 07:04:41
gamma 点本来就有虚频,对应长波极限。一般直接忽略不计

大神,你应该搞错了。。。问题不是在gamma点,是在gamma点附近的出现较大虚频,这跟文献上明显不对应啊。其次,在另一个好对称点处两光学支没有简并。是什么原因呢
7楼2016-08-03 11:01:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

坏逊平儿

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by liqizuiyang at 2016-08-03 09:37:36
重新优化结构,力收敛判据调至0.0001eV/A。

改用MP型展宽,做严格收敛测试。

K点改成18*18*1.

diagonalization改成david。cg只在空带数很大的时候才有优势,速度上并不占优。

请问改用mp展宽,做测试,能具体讲下细节吗,刚入门,诸多不太懂。。

发自小木虫Android客户端
8楼2016-08-03 11:27:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)

引用回帖:
8楼: Originally posted by 坏逊平儿 at 2016-08-03 11:27:34
请问改用mp展宽,做测试,能具体讲下细节吗,刚入门,诸多不太懂。。
...

可以参考vasp手册。测试流程就是不断增大展开,同时平均到每个原子上的TS项低于某个数值。

» 本帖已获得的红花(最新10朵)

9楼2016-08-03 14:21:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

坏逊平儿

新虫 (初入文坛)

送红花一朵
引用回帖:
9楼: Originally posted by liqizuiyang at 2016-08-03 14:21:59
可以参考vasp手册。测试流程就是不断增大展开,同时平均到每个原子上的TS项低于某个数值。...

谢谢耐心指导
10楼2016-08-03 16:23:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 坏逊平儿 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 327求调剂 +27 Xxjc1107. 2026-04-13 30/1500 2026-04-19 08:22 by cuisz
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 291求调剂 +10 关忆北. 2026-04-14 10/500 2026-04-18 23:32 by 路病情
[考研] 生物学调剂 +10 纸扇zhishan 2026-04-13 10/500 2026-04-18 21:19 by zqndavala
[考研] 297,工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 22408 312求调剂 +24 门路摸摸 2026-04-14 26/1300 2026-04-18 13:04 by wunaiy88
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考研] 化工学硕294分,求导师收留 +33 yzyzx 2026-04-12 37/1850 2026-04-17 23:00 by wunaiy88
[考研] 304求调剂 +7 castLight 2026-04-16 7/350 2026-04-17 20:05 by 关一盏灯cd
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[考研] 297,工科调剂?河南农业大学本科 +14 河南农业大学-能 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 求调剂学校 +14 不会吃肉 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 各位老师好,求调剂,本科211,一志愿天津大学生物与医药学硕,差两名录取。 +11 路六六jjj 2026-04-13 11/550 2026-04-14 16:01 by zs92450
[考研] 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 试管破裂
[考研] 考研英一数一338分 +9 长江大学东校区 2026-04-13 10/500 2026-04-14 00:41 by 王珺璞
[考研] 339求调剂 +4 hanwudada 2026-04-12 4/200 2026-04-13 12:03 by 蓝云思雨
信息提示
请填处理意见