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

24小时热门版块排行榜    

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

木虫 (正式写手)
引用回帖:
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 个回答

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 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭