24小时热门版块排行榜    

查看: 1351  |  回复: 2

twjdpnf287

铜虫 (初入文坛)

[求助] 磷单质的DOS计算方法 已有1人参与

各位大神,请指教一下磷单质的DOS计算方法。
我采用的是超软贗势的方法,我的INCAR设置为:
SYSTEM=P               ! Name of the system
ISTART=1                   ! startjo: 0-new 1-const 2-samecut
ISPIN=2                    ! spin polarized calculation (2-yes 1-no)
ICHARG=11                   !charge:1-file 2-atom 10-const
GGA=91
EDIFF=10E-6                ! stopping criterion for electr updt
EDIFFG=10E-10              ! stopping criterion for + energy - forces
ENCUT=400                   ! max cut off introduced based on previous calc.
#NSW=600                      ! number of steps for ionic update
#IBRION=2                 ! ionic relaxation:0-MD 1-quasi new 2-CG
#ISIF=2                    ! switch to control if calculate stress tensor etc. and what to relax?
ISMEAR=-5                   ! for semiconductor, and using 1 or 2 k points, otherwise <0,=1 for metals
SIGMA=0.2                  ! broadening in eV -4 tet -1 fermi 0-gaus
VOSKOWN=1                  ! use voski, Wilk , Nusair interpolation! create WAVECAR, CHGCAR,LOCPOT
RWIGS  =    1.233  

但是系统一直报错,out报错如下:
LDA part: xc-table for Ceperly-Alder, Vosko type interpolation para-ferro
WARNING: stress and forces are not correct
POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...
WAVECAR not read
charge-density read from file: B2                                      
magnetization density of overlapping atoms calculated
internal ERROR RSPHER:running out of buffer            0           0
           8           1           0
nonlr.F:Out of buffer RSPHER
internal ERROR RSPHER:running out of buffer            0           0
           8           1           0
nonlr.F:Out of buffer RSPHER
internal ERROR RSPHER:running out of buffer            0           0
           8           1           0
.......
求各位大神指点迷津,在此感激涕零。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

冯亚娟

金虫 (正式写手)

【答案】应助回帖

尝试跟换ISMEAR=0与SIGMA=0.05或者SIGMA=0.1,看看是不是会有变化?我最近也在计算DOS 的时候出现问题,这样更改后就对了不知道适不适合你的体系。
2楼2016-04-08 11:03:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scuedu2

金虫 (小有名气)

为什么用超软魇势
3楼2016-04-08 16:11:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 twjdpnf287 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:47 by oEVWOejN9taj
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[基金申请] 关于代码变化问题,想知道的进来 +11 且听虎啸 2026-08-07 14/700 2026-08-08 12:50 by 布布和一二
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&amp;fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见