24小时热门版块排行榜    

查看: 7642  |  回复: 16

笑笑愁

金虫 (小有名气)

未知

引用回帖:
9楼: Originally posted by bird007 at 2014-11-12 17:02:30
我的是用VASP算的电荷密度,然后bader得到的
你的是MS不知道和VASP有区别没...

你的VASP的输入文件能传给我一个吗
有志者,事竟成。
11楼2014-11-12 17:16:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)

引用回帖:
3楼: Originally posted by 笑笑愁 at 2014-11-12 09:46:22
I still don't know what to do. Can you be more specific ? I mean which software should be used, what method is needed, what steps could be followed. Many thanks....

自带NaCl例子的INCAR:
SYSTEM = NaCl charge
IBRION=2
POTIM=0.1
ISYM=0
#ISMEAR=-1;SIGMA=.1
ISIF=0
ISPIN=2
NSW=0
EDIFF=0.1E-04
EDIFFG=-0.001
#ENMAX=250
PREC=med
LWAVE=.FALSE.
LCHARG=.TRUE.
IALGO=48
NGXF=120
NGYF=120
NGZF=120
LAECHG = .TRUE.

-------------------------
主要是增加NG(X,Y,Z)F 的值直到total charge is correct。
这时NUMBER OF ELECTRONS 为整数,反之NUMBER OF ELECTRONS 带小数,则说明总电荷没有对,需要增加NG(X,Y,Z)F 的值。
NGXF=120
NGYF=120
NGZF=120
12楼2014-11-13 07:57:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bird007

荣誉版主 (职业作家)

踩单车的季节

引用回帖:
11楼: Originally posted by 笑笑愁 at 2014-11-12 17:16:30
你的VASP的输入文件能传给我一个吗...

和静态计算比 就多了个
LAECHG=.TRUE.
一日,老师在课堂上教偶认字——帅.偶百思不得其解,正在这时,同桌mm递来一面镜子,偶照之,顿时恍然大悟......[此广告位招商!][此广告位招商!]
13楼2014-11-13 22:03:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
笑笑愁: 金币+10, ★★★很有帮助 2014-11-15 10:44:21
NUMBER OF ELECTRONS:         8.6116
可以肯定你的计算是错的。。。。。。这个应该是体系的总电子数,你这么多原子怎么可能还不到9个电子。
另外我建议不要这么换,直接用vasp算比较好。
14楼2014-11-14 21:00:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

笑笑愁

金虫 (小有名气)

未知

引用回帖:
14楼: Originally posted by dxcharlary at 2014-11-14 21:00:56
NUMBER OF ELECTRONS:         8.6116
可以肯定你的计算是错的。。。。。。这个应该是体系的总电子数,你这么多原子怎么可能还不到9个电子。
另外我建议不要这么换,直接用vasp算比较好。

用VASP怎么算,能传个输入文件看看吗
有志者,事竟成。
15楼2014-11-15 10:46:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
笑笑愁: 金币+10, ★★★很有帮助 2014-11-16 10:02:22
引用回帖:
15楼: Originally posted by 笑笑愁 at 2014-11-15 10:46:31
用VASP怎么算,能传个输入文件看看吗...

就是正常的单点能计算,加个LAECHG=.TRUE.
然后用austin的脚本,具体用法见:
http://theory.cm.utexas.edu/henkelman/code/bader/
16楼2014-11-15 14:34:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wd12309812

铜虫 (正式写手)

引用回帖:
12楼: Originally posted by Illusionist at 2014-11-13 07:57:41
自带NaCl例子的INCAR:
SYSTEM = NaCl charge
IBRION=2
POTIM=0.1
ISYM=0
#ISMEAR=-1;SIGMA=.1
ISIF=0
ISPIN=2
NSW=0
EDIFF=0.1E-04
EDIFFG=-0.001
#ENMAX=250
PREC=med
LWAVE=.FALSE.
LCHARG=.TRUE ...

请问,总电子数是从哪里看
bader生成的ACF.dat最后一行,还是OUTCAR里'total charge'呢?
17楼2016-05-24 11:43:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 笑笑愁 的主题更新
信息提示
请填处理意见