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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 金属性二位材料声子计算（拉曼）中虚频无法消除
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1658  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

c_chi

木虫 (正式写手)

[求助] 金属性二位材料声子计算（拉曼）中虚频无法消除 已有1人参与

各位虫子好，我现在在用PWscf计算二位材料V2C的拉曼，在计算声子的过程中总会出现一点虚频，如下所示。
第一个问题：我尝试调高晶格优化、电子自洽的精度、ecutwfc 、KPOINTS，以及降低degauss，但是都不能消除这个低频，请问有什么解决办法吗？输入见后面。而且SCF找到的对称差性一直是C_1 (1)，这里有什么问题没有，如何解决？另外，看到有些帖子提到二位材料声子计算的虚频很难消除，但如果足够小可以忽略，请问这里的虚频可以忽略吗？

第二个问题：我用pwscf计算拉曼的方法来自于帖子http://muchong.com/bbs/viewthread.php?tid=1999091&fpage=1&target=blank
但我的体系是金属性的，所以在声子就算时epsil和lraman参数是不能打开的，但愿贴中特意提到这两个参数，请问这样计算出来的声子频率可以最后用来拉曼的计算吗？
或者有什么解决的办法没有？

谢谢！

############频率输出#########
**************************************************************************
     freq (    1) =      -0.685249 [THz] =     -22.857456 [cm-1]
     freq (    2) =       0.477902 [THz] =      15.941095 [cm-1]
     freq (    3) =       0.877071 [THz] =      29.255946 [cm-1]
     freq (    4) =       6.089501 [THz] =     203.123876 [cm-1]
     freq (    5) =       6.104889 [THz] =     203.637164 [cm-1]
     freq (    6) =      10.123378 [THz] =     337.679526 [cm-1]
     freq (    7) =      18.108704 [THz] =     604.041362 [cm-1]
     freq (    8) =      20.018180 [THz] =     667.734614 [cm-1]
     freq (    9) =      20.078524 [THz] =     669.747470 [cm-1]
**************************************************************************

     Mode symmetry, C_1 (1)     point group:

     freq (  1 -  1) =        -22.9  [cm-1]   --> A               I+R
     freq (  2 -  2) =         15.9  [cm-1]   --> A               I+R
     freq (  3 -  3) =         29.3  [cm-1]   --> A               I+R
     freq (  4 -  4) =        203.1  [cm-1]   --> A               I+R
     freq (  5 -  5) =        203.6  [cm-1]   --> A               I+R
     freq (  6 -  6) =        337.7  [cm-1]   --> A               I+R
     freq (  7 -  7) =        604.0  [cm-1]   --> A               I+R
     freq (  8 -  8) =        667.7  [cm-1]   --> A               I+R
     freq (  9 -  9) =        669.7  [cm-1]   --> A               I+R


########自洽计算输入##########
&control
    calculation='scf',
    restart_mode='from_scratch',
    pseudo_dir = './',
    outdir= './',
    prefix='V2C',
    tprnfor = .true. ,
    tstress = .true. ,
    etot_conv_thr = 1.0D-5,    !Default: 1.0D-4  (a.u)
    forc_conv_thr = 1.0D-4,    !Default: 1.0D-3  (a.u)
    nstep = 300,  !number of ionic + electronic steps
    verbosity = 'high' ,
/
&system
    ibrav = 4,  celldm(1)=5.4819,     celldm(3)=6.58615,        !celldm(1) = 1/0.529177,
    nat = 3, ntyp=2,
    ecutwfc = 80.0,              ! kinetic energy cutoff (Ry) for wavefunctions
!   ecutrho = 500.0,             ! kinetic energy cutoff (Ry) for charge density and potential;Default:4 * ecutwfc
    occupations='smearing',      ! 'tetrahedra' especially suited for calculation of DOS,Not suitable for force/optimization/dynamics calculations
    smearing='mp', degauss=0.01,  ! 'fermi-dirac'
    nosym = .false. ,
/
&electrons
    electron_maxstep=500,
    diagonalization = 'david'  ! cg
    mixing_mode     = 'local-TF'  ! TF
    mixing_beta     = 0.3,
    conv_thr        = 1.D-12,
/
ATOMIC_SPECIES
V    23   V.pw-mt_fhi.UPF
C     6   C.pw-mt_fhi.UPF
ATOMIC_POSITIONS {crystal}    !alat | bohr | angstrom | crystal ; 0, fixed
V        0.666818248   0.333332706   0.563114857
V        0.333485971   0.666664798   0.447286012
C        0.000159931  -0.000001312   0.505200379
K_POINTS {automatic}
27 27 1 0 0 0

########声子计算输入##########
phonon calculation
&inputph
  tr2_ph=1.0d-12,
  prefix='V2C',
! epsil=.true.,
  trans=.true.,
! lraman=.true.,
  amass(1)=50.9415,
  amass(2)=12.011,
  outdir='./',
  fildyn='V2C.dynG',
  fildrho='V2C.drho',
/
0.0 0.0 0.0
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2016-07-11 04:55:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c_chi

木虫 (正式写手)
引用回帖:
2楼: Originally posted by dxcharlary at 2016-07-11 21:33:20
如果是gamma点处的虚频，那其实不用消。正常声学支在gamma点处为零。如果硬要消的话，加ecutrho和ph的收敛精度。
普通拉曼的计算方法实际上对金属不适用。同时理论上，金属对拉曼红外全吸收，没有拉曼和红外信号。 ...

谢谢您的回复。
加ecutrho收敛精度貌似效果不明显，我再试试ph收敛精度。
这里occupations改成fixed意思是强制把金属性材料作为半导体来计算吗？
赞一下 回复此楼
3楼2016-07-11 21:54:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
c_chi: 金币+50, ★★★很有帮助 2016-07-11 21:55:08
漫天飘雪: 金币+2, 应助指数+1, 谢谢交流 2016-07-12 10:17:09
如果是gamma点处的虚频，那其实不用消。正常声学支在gamma点处为零。如果硬要消的话，加ecutrho和ph的收敛精度。
普通拉曼的计算方法实际上对金属不适用。同时理论上，金属对拉曼红外全吸收，没有拉曼和红外信号。当然实验上金属可以做表面增强。当然硬要算的话，把occupations改成fixed，对于大多数金属性不是很好的材料是适用的。
另外提醒一句V体系可能需要加强关联。
赞一下 回复此楼
2楼2016-07-11 21:33:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

引用回帖:
3楼: Originally posted by c_chi at 2016-07-11 21:54:50
谢谢您的回复。
加ecutrho收敛精度貌似效果不明显，我再试试ph收敛精度。
这里occupations改成fixed意思是强制把金属性材料作为半导体来计算吗？...

是的，你可以这么理解。
赞一下 回复此楼
4楼2016-07-11 22:26:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c_chi

木虫 (正式写手)
引用回帖:
4楼: Originally posted by dxcharlary at 2016-07-11 22:26:49
是的，你可以这么理解。...

还有个问题，我用的是QE官网提供的由ab inito的模守恒赝势转换过来的赝势，您知道这个的可靠性如何吗？
赞一下 回复此楼
5楼2016-07-11 23:11:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学308分求调剂（一志愿华东师大） +6 相信必会光芒万� 2026-03-31 7/350 2026-04-02 23:16 by JourneyLucky
[考研] 298分 070300求调剂 +8 zwen03 2026-04-02 8/400 2026-04-02 21:55 by barlinike
[考研] 266求调剂 +12 阳阳哇塞 2026-03-27 13/650 2026-04-02 21:32 by zllcz
[考研] 296求调剂 +4 sdhu 2026-04-02 4/200 2026-04-02 21:29 by baoball
[考研] 调剂求助 +8 想换手机不想解� 2026-04-02 8/400 2026-04-02 20:47 by dongzh2009
[考研] 085602化工求调剂（331分） +9 111@127 2026-03-30 9/450 2026-04-02 20:00 by dick_runner
[考研] 413求调剂 +3 柯某某 2026-03-31 3/150 2026-04-02 16:59 by zzsw+
[考研] 346求调剂 +5 郑诚乐 2026-04-02 5/250 2026-04-02 16:38 by SZW_UJN
[考研] 07生物学求调剂 一志愿同济大学359分 +3 LAMC. 2026-03-30 3/150 2026-04-02 10:26 by 18828373951
[考研] 372求调剂 +3 jj涌77 2026-04-02 3/150 2026-04-02 09:57 by olim
[考研] 生物学327，求调剂 +5 书上的梅子 2026-04-01 6/300 2026-04-02 06:47 by ilovexiaobin
[考研] 070300一志愿211，312分求调剂院校 +14 小黄鸭宝 2026-03-30 14/700 2026-04-01 20:19 by 赖春艳
[考研] 320分，材料与化工专业，求调剂 +14 一定上岸aaa 2026-03-27 18/900 2026-04-01 20:10 by 积极调剂的小学�
[考研] 材料科学与工程339求调剂 +11 hyz0119 2026-03-31 12/600 2026-04-01 18:40 by 伟大河北
[考研] 267求调剂 +13 uiybh 2026-03-31 13/650 2026-04-01 10:25 by 探123
[考研] 329求调剂，一志愿西北工业大学，材料工程（085601） +6 小小机灵虫 2026-03-29 12/600 2026-03-31 16:58 by 记事本2026
[考研] 一志愿南昌大学324求调剂 +6 hanamiko 2026-03-29 6/300 2026-03-31 16:35 by hypershenger
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
[考研] 298求调剂 +3 种圣赐 2026-03-29 3/150 2026-03-29 12:06 by longlotian
[考研] 298求调剂 +4 种圣赐 2026-03-28 4/200 2026-03-29 08:42 by q1092522407
信息提示
关闭
请填处理意见
关闭