| 查看: 1788 | 回复: 6 | |||||
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 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
dxcharlary
专家顾问 (职业作家)
-
专家经验: +83 - 1ST强帖: 2
- 应助: 251 (大学生)
- 金币: 10505.9
- 散金: 155
- 红花: 65
- 帖子: 3141
- 在线: 758.7小时
- 虫号: 936718
- 注册: 2010-01-04
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
c_chi: 金币+50, ★★★很有帮助 2016-07-11 21:55:08
漫天飘雪: 金币+2, 应助指数+1, 谢谢交流 2016-07-12 10:17:09
感谢参与，应助指数 +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
c_chi
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2449.7
- 散金: 228
- 红花: 2
- 帖子: 551
- 在线: 493.6小时
- 虫号: 1253512
- 注册: 2011-04-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2016-07-11 21:54:50
dxcharlary
专家顾问 (职业作家)
-
专家经验: +83 - 1ST强帖: 2
- 应助: 251 (大学生)
- 金币: 10505.9
- 散金: 155
- 红花: 65
- 帖子: 3141
- 在线: 758.7小时
- 虫号: 936718
- 注册: 2010-01-04
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
4楼2016-07-11 22:26:49
c_chi
木虫 (正式写手)
- 应助: 17 (小学生)
- 金币: 2449.7
- 散金: 228
- 红花: 2
- 帖子: 551
- 在线: 493.6小时
- 虫号: 1253512
- 注册: 2011-04-02
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
5楼2016-07-11 23:11:43
6楼2016-07-14 09:01:58
dxcharlary
专家顾问 (职业作家)
-
专家经验: +83 - 1ST强帖: 2
- 应助: 251 (大学生)
- 金币: 10505.9
- 散金: 155
- 红花: 65
- 帖子: 3141
- 在线: 758.7小时
- 虫号: 936718
- 注册: 2010-01-04
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
7楼2016-07-20 16:32:46
20