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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 硅烯声子谱G点附近虚频
301/3返回列表123下一页
查看: 6039  |  回复: 29
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

坏逊平儿

新虫 (初入文坛)

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

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

（硅烯的结构是组里已经优化好的结构，算声子谱直接做自洽计算吧？！）
跪求大神指导，万分感谢。

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

» 猜你喜欢
1楼 2016-08-02 15:28:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有3个 )

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以上。
赞一下(1人) 回复此楼
先进使役材料计算（https://aidme.nimte.ac.cn）
14楼2016-08-04 00:21:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

【答案】应助回帖

漫天飘雪: 回帖置顶 2018-06-24 13:12:03
另外，声子计算用什么参数，结构优化也最好用一样严格的参数，另外，力的收敛最好达到d-4。
赞一下 回复此楼
先进使役材料计算（https://aidme.nimte.ac.cn）
15楼2016-08-04 00:24:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

【答案】应助回帖


漫天飘雪: 回帖置顶 2018-06-24 13:12:05
当你nq是3的倍数或足够密的的时候，你自然会得到简并了。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

坏逊平儿
先进使役材料计算（https://aidme.nimte.ac.cn）
16楼2016-08-04 00:26:00
已阅   回复此楼   关注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只在空带数很大的时候才有优势，速度上并不占优。
赞一下(1人) 回复此楼
6楼2016-08-03 09:37:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c19940408

银虫 (初入文坛)

【答案】应助回帖

★
蛋蛋小童鞋: 金币+1, 谢谢回答 2016-08-11 08:37:15
建议使用和优化时相近的参数来计算。。比如相同的能量截断，k点等等
赞一下(1人) 回复此楼
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的回帖
相关版块跳转 我要订阅楼主 坏逊平儿 的主题更新
301/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 还有化工二轮调剂的学校吗 5+13 化工人999 2026-04-09 47/2350 2026-04-11 00:20 by 骑牛渡寒江
[考研] 生物学308分求调剂（一志愿华东师大）做过分子实验 +8 相信必会光芒万� 2026-04-07 9/450 2026-04-10 21:03 by zhouxiaoyu
[考研] 材料调剂 +13 一样YWY 2026-04-04 13/650 2026-04-10 11:07 by mattzhming
[考研] 江苏大学 工科调剂 捡漏 +3 Evan_Liu 2026-04-09 5/250 2026-04-10 10:22 by Evan_Liu
[考研] 一志愿中科大070300化学，314分求调剂 +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
[考研] 一志愿华工085600 331分 +6 天下ww 2026-04-09 6/300 2026-04-09 18:59 by l_paradox
[考研] 085501机械专硕 302分 不挑专业求调剂 +5 汪某. 2026-04-09 5/250 2026-04-09 15:38 by 蒋皓禹
[考研] 353求调剂 +8 晴空万里air 2026-04-07 8/400 2026-04-09 00:18 by GouQ
[考研] 0703化学调剂 348分 +14 唉我超真没招了 2026-04-06 15/750 2026-04-08 19:16 by 我减肥1
[考研] 一志愿吉大化学327求调剂 +12 王王白石 2026-04-06 13/650 2026-04-08 16:05 by luoyongfeng
[考研] 304求调剂 +16 c297914 2026-04-05 17/850 2026-04-08 13:00 by grayjzr
[考研] 22408 266求调剂 +11 masss11222 2026-04-07 14/700 2026-04-08 11:06 by yulian1987
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 材料调剂 +13 汉123456 2026-04-07 14/700 2026-04-07 22:53 by 来看流星雨10
[考研] 312求调剂 +4 LR6 2026-04-06 4/200 2026-04-07 08:42 by jp9609
[考研] 一志愿太原理工大学计算机技术专硕348，求调剂指导 +3 nexious 2026-04-05 3/150 2026-04-07 08:19 by jp9609
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
[考研] 专硕304找调剂，一线城市最好 +3 李lsl李 2026-04-05 3/150 2026-04-06 12:16 by ffffjjjj
[考研] 一志愿南航，数一英一学硕317求调剂！！ +6 Acaciad 2026-04-04 6/300 2026-04-06 12:13 by 考研学校招点人
信息提示
关闭
请填处理意见
关闭