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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问态密度怎么算
51/1返回列表
查看: 1306  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wd12309812

铜虫 (正式写手)

[求助] 请问态密度怎么算 已有1人参与

我按照教程上面，驰豫完之后用
IBRION = -1; NSW = 0
静态计算，得到DOSCAR文件只有几行东西 ，没有太密度的任何信息信息为什么阿？？

把vasprun.xml用p4v打开，也什么都没有

求解！！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-04-16 15:06:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wd12309812

铜虫 (正式写手)
INCAR如下

#NWRITE    =    2     LPETIM=.TURE.        # write-flag&timer
#INIWAV   =    1                          # INIWAV, default is 1,only useful when Istart=0
ISTART    =    0        #job : 0-new  1-constant cut-off  2-constant cell size/shape.and cutoff read from WAVECAR
ICHARG    =    2        #job : 0-from initial wf  1-from CHGCAR   2-atomic charge density , relax/DOS=2  EK=11
LCHARG    =  .T.
LWAVE     =   .F.
ISPIN     =    1        
#ISYM      =    0        #default 1-US-pp  2-PAW  symmetry stuff ON(1 or 2) or OFF(0)
PREC      =   Normal

#LORBIT    =    11      # for pdos
ENCUT     =   300 eV

Electronic Relaxation
NELMIN    =   4
NELM      =   200       #   the max number of electronic SC
EDIFF     =    1E-05    #   stopping-criterion for ELM
LREAL     =    Auto     #   F-reciprocal  T-real  A-auto

Ionic Relaxation
IBRION    =    -1        #   ionic relaxation:0-MD 1-quasi-Newton 2-CG   DOS/EK =-1
#ISIF      =    2
NSW       =    0    #   DOS/EK=0  relax = 400
EDIFFG    =    -0.01
#POTIM     =    0.01

LVDW = .TRUE.

#NEB related tags
#SPRING=-5
#IMAGES=3
#ICHAIN=0
#LCLIMB=T

DOS related values
ISMEAR    =    0
SIGMA     =    0.05
NEDOS     =    1200   # for more dos
LORBIT = 11

ALgorithm
ALGO      =    FAST

#NPAR=4

#spin-orbital interaction
#LSORBIT   = .TRUE.
#SAXIS     = 0 0 1
#MAGMOM    =
#GGA_COMPAT = .FALSE.
#LMAXMIX    =  4
#ADDGRID    = .TRUE.
赞一下 回复此楼
2楼2015-04-16 15:16:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

未岸亡灵

铁杆木虫 (著名写手)

暂时离开小木虫

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
wd12309812: 金币+20, ★★★★★最佳答案 2015-04-16 21:43:08
算dos的时候，先自洽，自洽得到电荷密度文件CHGCAR，在令ICHARG =11，即初始电荷密度从CHGCAR中读取。在算出结果。另外建议算dos时，令ISMEAR=-5更准确，这是官方建议。

最后，就是检查算完你的OUTCAR有没有正常结束
赞一下 回复此楼
3楼2015-04-16 18:19:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wd12309812

铜虫 (正式写手)
引用回帖:
3楼: Originally posted by 未岸亡灵 at 2015-04-16 18:19:45
算dos的时候，先自洽，自洽得到电荷密度文件CHGCAR，在令ICHARG =11，即初始电荷密度从CHGCAR中读取。在算出结果。另外建议算dos时，令ISMEAR=-5更准确，这是官方建议。

最后，就是检查算完你的OUTCAR有没有正常 ...

非常感谢，检查了以前的错误原因发现是K点取得过多节点down掉了，超胞大概10*10*20，K点取15*15*3
K点很多吗 ？可是看说明书说应该足够大才是啊？

改到3*3*1已可以正常运行
现在又遇到一个问题，自恰得到CHGCAR，但是ICHARGE=11时任务出错，output显示
"ERROR: charge density could not be read from file CHGCAR for ICHARG>10"
请问这是什么原因？
赞一下 回复此楼
4楼2015-04-16 21:42:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

未岸亡灵

铁杆木虫 (著名写手)

暂时离开小木虫
引用回帖:
4楼: Originally posted by wd12309812 at 2015-04-16 21:42:58
非常感谢，检查了以前的错误原因发现是K点取得过多节点down掉了，超胞大概10*10*20，K点取15*15*3
K点很多吗 ？可是看说明书说应该足够大才是啊？

改到3*3*1已可以正常运行
现在又遇到一个问题，自恰得到CHGC ...

k点网格大小，这个不能定论，得根据原子数和对称性衡量，原子数少可以多取，反之，少取。对称性高，可以多取，反之，少取。最后实际取到的不可约k点坐标应该在IBZKPT这个文件中能看到。

对于后一个错误，显然，你应该保证在自洽的时候LCHARG=T以便输出电荷密度文件，检查其是否为空。若不为空，将其一同放入计算dos的文件夹中，用于读取
赞一下 回复此楼
5楼2015-04-17 07:56:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wd12309812 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭