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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 石墨烯声子谱Gamma点附近虚频，怎么处理
北京石油化工学院2026年研究生招生接收调剂公告
61/1返回列表
查看: 4767  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

justinzhouke

新虫 (初入文坛)

[交流] 石墨烯声子谱Gamma点附近虚频，怎么处理 已有4人参与

计算石墨烯的声子谱，计算发现Gamma点附近出现虚频（如下)。不知道怎么处理。

  1   0.000000       0.0000    0.0000    0.0000  879.5294 1560.4160 1560.4160
  2   0.022223      -3.7499   41.7646   65.5168  879.2473 1559.8164 1561.6153
  3   0.044447      -6.5222   83.4849  130.6515  878.3968 1557.9969 1565.0765
  4   0.066670      -6.7089  125.1109  195.0346  876.9649 1554.9026 1570.4069
  5   0.088893       6.6043  166.5984  258.3457  874.9310 1550.4615 1577.0073
  6   0.111117      17.6058  207.9030  320.3191  872.2675 1544.6098 1584.1435
  7   0.133340      29.2947  248.9746  380.7432  868.9412 1537.3218 1591.0273

基本设置包括PP，是按照论文里面来的（PHYSICAL REVIEW B 80, 033406 2009）。
我做过ecut和q-mesh的收敛性测试，q-mesh在16的时候就收敛了。Ecut是55Ry收敛。


In the present analysis, we take advantage of the PWSCF
package of the QUANTUM-ESPRESSO distribution13 to model
the characteristic graphene lattice dynamics. Specifically,
an ultrasoft pseudopotential in the local-density
approximation14 is used with a plane-wave expansion up to
55 Ry and 32x32x1 Monkhorst-Pack meshes adopted for
Brillouin-zone sampling with 0.02 Ry Fermi-Dirac smearing
in the electronic occupation. The calculation also employs
the theoretical lattice parameter a=2.46 Å, which is in excellent
agreement with the experimental value at 300 K. For
bilayer graphene, the Bernal stacking is considered with the
interlayer distance set to the theoretical minimum of
3.36 Å.

输入文件
Gra.scf.in:  
1  &control
  2         calculation = 'scf',
  3         restart_mode= 'from_scratch',
  4         pseudo_dir  = '/home/xuzp/bin/lib/pwscf_pot',
  5         outdir      = './tmp',
  6         prefix='graphene',
  7         !disk_io     = 'none',
  8         tstress = .true.
  9         tprnfor = .true.
10         wf_collect=.true.
11 /
12 &system
13     ibrav           = 4
14     celldm(1)       = 4.64873
15     celldm(3)       = 8.13
16     nat             = 2
17     ntyp            = 1
18     ecutwfc         = 55
19     ecutrho         = 550
20     occupations     = 'smearing'
21     smearing        = 'fd'
22     degauss         = 0.02
23 /
24 &electrons
25     diagonalization = 'david'
26     mixing_beta     = 0.7
27     conv_thr        = 1.0d-10
28 /
29 ATOMIC_SPECIES
30   C  12.01078 C.pbe-rrkjus.UPF  !C.pz-vbc.UPF
31 ATOMIC_POSITIONS crystal
32   C   0.6666666666666667   0.3333333333333333   0.5000000000000000
33   C   0.3333333333333333   0.6666666666666667   0.5000000000000000
34 K_POINTS automatic
35 32 32 1 0 0 0

Gra.ph.in:
  1 phonon for graphene
  2 &inputph
  3         tr2_ph          =1.0d-18,
  4         prefix          ='graphene',
  5         epsil           =.false.,
  6         amass(1)        =12.01078,
  7         outdir          = './tmp',
  8         fildyn='graphene.dynG',
  9         !elph=.false.,
10         !trans=.true.,
11         ldisp=.true.
12         nq1=6, nq2=6, nq3=1
13 /
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
whoisyourdady
1楼 2015-05-05 14:11:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
声子计算总是需要先结构优化，达到足够高的精度，否则由于计算条件略有区别（赝势等）会产生虚频。
赞一下 回复此楼
2楼2015-05-05 14:32:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lchchl

铜虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
Gamma点附近的一点虚频，也许改一下声学声子求和规则就可以消除了。
赞一下 回复此楼
为什么小木虫每天都是扣四个金币！！！！！！
3楼2015-05-18 15:10:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magicmonk

至尊木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
3楼: Originally posted by lchchl at 2015-05-18 15:10:06
Gamma点附近的一点虚频，也许改一下声学声子求和规则就可以消除了。

此话怎讲，求详细说明求和规则
赞一下 回复此楼
我们最深的恐惧不是我们能力不够，我们最深的恐惧是能力超越了极限。
4楼2015-05-18 15:33:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lchchl

铜虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
4楼: Originally posted by magicmonk at 2015-05-18 15:33:50
此话怎讲，求详细说明求和规则...

有个参数叫zasr，这个是控制求和规则的，印象中可以设置的求和规则包括，no，simple，crystal等。
赞一下 回复此楼
为什么小木虫每天都是扣四个金币！！！！！！
5楼2015-05-18 19:44:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

三叶虫的脚印

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主你好。请问你Gamma点的虚频消掉了吗？我最近计算声子谱时也出现了相同的问题。Gamma点处总是有3个小虚频。。
赞一下(1人) 回复此楼
6楼2016-12-22 20:17:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 justinzhouke 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085701求调剂初试286分 +4 secret0328 2026-03-28 4/200 2026-03-28 21:09 by 15366876211
[考研] 生物学学硕，一志愿湖南大学，初试成绩338 +6 YYYYYNNNNN 2026-03-26 7/350 2026-03-28 20:52 by 唐沐儿
[考研] 311（085601）求调剂 +4 liziyeyeye 2026-03-28 4/200 2026-03-28 18:50 by 535743368
[考研] 291求调剂 +15 hhhhxn.. 2026-03-23 21/1050 2026-03-28 11:26 by self2008
[考研] 086502化学工程342求调剂 +6 阿姨复古不过 2026-03-27 6/300 2026-03-28 07:06 by wangy0907
[考研] 复试调剂，一志愿南农083200食品科学与工程 +5 XQTJZ 2026-03-26 5/250 2026-03-27 14:49 by 狂炫麦当当
[考研] 一志愿华东理工大学081700，初试分数271 +6 kotoko_ik 2026-03-23 7/350 2026-03-27 12:29 by 惠州彭于晏
[考研] 292求调剂 +4 求求了收下我吧� 2026-03-26 4/200 2026-03-27 10:37 by zhshch
[考研] 324求调剂 +5 hanamiko 2026-03-26 5/250 2026-03-27 10:33 by wangjy2002
[考研] 【双一流院校新能源、环境材料，材料加工与模拟招收大量调剂】 +4 Higraduate 2026-03-22 8/400 2026-03-26 20:34 by Higraduate
[考研] 总分322求生物学/生化与分子/生物信息学相关调剂 +5 星沉uu 2026-03-26 6/300 2026-03-26 19:02 by macy2011
[考研] 材料考研求调剂 +3 Dendel 2026-03-23 6/300 2026-03-26 17:51 by fmesaito
[考研] 一志愿河工大 081700 276求调剂 +4 地球绕着太阳转 2026-03-23 4/200 2026-03-26 14:27 by zzll406
[考研] 309求调剂 +4 gajsj 2026-03-25 5/250 2026-03-26 00:27 by Dyhoer
[考研] 0854电子信息求调剂 324 +4 Promise-jyl 2026-03-23 4/200 2026-03-25 11:36 by Sugarlight
[考研] 289材料与化工（085600）B区求调剂 +4 这么名字咋样 2026-03-22 5/250 2026-03-25 08:20 by mx.yue
[有机交流] 有机合成求助 20+3 FENGSHUJEI 2026-03-23 5/250 2026-03-24 19:31 by 88817753
[考研] 一志愿南航材料专317分求调剂 +5 炸呀炸呀炸薯条 2026-03-23 5/250 2026-03-24 16:52 by 星空星月
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[考研] 石河子大学（211、双一流）硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
信息提示
关闭
请填处理意见
关闭