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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 硅烯声子谱G点附近虚频
51/1返回列表
查看: 6465  |  回复: 29
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

坏逊平儿

新虫 (初入文坛)

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

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

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

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

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

golddoushi

木虫 (正式写手)

【答案】应助回帖

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

坏逊平儿

新虫 (初入文坛)
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的回帖
查看全部 30 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 6/300 2026-08-31 02:35 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一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 6/300 2026-08-31 02:23 by hZiFeudZyoGR
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +7 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-31 01:15 by hZiFeudZyoGR
[基金申请] 国社科又开始会评了，不知道这次命运如何 +5 雨打竹帘 2026-08-30 9/450 2026-08-30 22:23 by 余韵清
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:30 by hZiFeudZyoGR
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 4FFAWE8HcgUD 2026-08-29 5/250 2026-08-30 20:29 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 19:57 by hZiFeudZyoGR
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:50 by hZiFeudZyoGR
[博后之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 18:28 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-30 17:56 by hZiFeudZyoGR
[论文投稿] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 17:03 by hZiFeudZyoGR
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 06:10 by ZPa0EcMwuECS
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
关闭
请填处理意见
关闭