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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 金属性二位材料声子计算（拉曼）中虚频无法消除
51/1返回列表
查看: 1802  |  回复: 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的回帖

司马星空

银虫 (初入文坛)
引用回帖:
5楼: Originally posted by c_chi at 2016-07-11 23:11:43
还有个问题，我用的是QE官网提供的由ab inito的模守恒赝势转换过来的赝势，您知道这个的可靠性如何吗？...

可以尝试下超软赝势
赞一下 回复此楼
6楼2016-07-14 09:01:58
已阅   回复此楼   关注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的回帖

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的回帖

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的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 关于Filecode分析方法 +8 majunge000 2026-08-10 9/450 2026-08-12 22:17 by 虫友是什么虫
[基金申请] FileCode能看出啥？ +8 要乐观耀哥 2026-08-10 21/1050 2026-08-12 16:08 by 云上清扬
[基金申请] 长年满屏的广告，版主太不责任了。基金也等的急 +4 gltch 2026-08-12 4/200 2026-08-12 16:08 by eulota
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果？ +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
信息提示
关闭
请填处理意见
关闭