24小时热门版块排行榜    

查看: 1756  |  回复: 6

lorna639

木虫 (著名写手)

[求助] bader电荷分析问题

请问我用bader分析电荷,为什么得不到AECCAR0,AECCAR1和AECCAR2,只能得到ACF.dat, AVF.dat和BCF.dat?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
ben_ladeng: 金币+2, 辛苦了 2016-07-23 08:04:12
AECCAR0 2是vasp自身产生的。http://theory.cm.utexas.edu/henkelman/code/bader/

“By adding the LAECHG=.TRUE. to the INCAR file, the core charge is written to AECCAR0 and the valance charge to AECCAR2. These two charge density files can be summed using the chgsum.pl script:

  chgsum.pl AECCAR0 AECCAR2”
不一定挂在论坛,计算问题问题欢迎留言。
2楼2016-07-12 09:21:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lorna639

木虫 (著名写手)

引用回帖:
2楼: Originally posted by 卡开发发 at 2016-07-12 09:21:32
AECCAR0 2是vasp自身产生的。http://theory.cm.utexas.edu/henkelman/code/bader/

“By adding the LAECHG=.TRUE. to the INCAR file, the core charge is written to AECCAR0 and the valance charge to AECCAR ...

我也看到这段文字了,不过我在INCAR里面加入LAECHG=.TRUE.,结果输出错误,去掉这个关键词后一切正常,不知道怎么回事啊,我用的是vasp5
3楼2016-07-12 09:59:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
3楼: Originally posted by lorna639 at 2016-07-12 09:59:25
我也看到这段文字了,不过我在INCAR里面加入LAECHG=.TRUE.,结果输出错误,去掉这个关键词后一切正常,不知道怎么回事啊,我用的是vasp5...

INCAR和跳错信息都发来。
不一定挂在论坛,计算问题问题欢迎留言。
4楼2016-07-12 10:25:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lorna639

木虫 (著名写手)

引用回帖:
4楼: Originally posted by 卡开发发 at 2016-07-12 10:25:03
INCAR和跳错信息都发来。...

System = Cr2AlB2 P mmm
ENCUT  = 700
ISTART = 0; ICHARG = 2
IBRION = -1; ISIF = 2; NSW=0;
ISMEAR = 0; SIGMA = 0.2;
PREC   = Accurate
LAECHG = .TRUE.
LCHARG = .TRUE.
这是输入文件,在计算时,在Scratch文件夹中可以看到AECCAR0和  AECCAR1,但是没有AECCAR2。计算结束时,解压EXTRAFILE,发现里面AECCAR0  AECCAR1也没有了,还请你予以指导,多谢!
5楼2016-07-13 06:11:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
5楼: Originally posted by lorna639 at 2016-07-13 06:11:23
System = Cr2AlB2 P mmm
ENCUT  = 700
ISTART = 0; ICHARG = 2
IBRION = -1; ISIF = 2; NSW=0;
ISMEAR = 0; SIGMA = 0.2;
PREC   = Accurate
LAECHG = .TRUE.
LCHARG = .TRUE.
这是输入文件,在计算 ...

我不知道为啥你会输出错误,还是看看输出文件,还有你得确定版本是否符合以及是不是采用了PAW。

[ 发自手机版 http://muchong.com/3g ]
不一定挂在论坛,计算问题问题欢迎留言。
6楼2016-07-13 16:21:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
5楼: Originally posted by lorna639 at 2016-07-13 06:11:23
System = Cr2AlB2 P mmm
ENCUT  = 700
ISTART = 0; ICHARG = 2
IBRION = -1; ISIF = 2; NSW=0;
ISMEAR = 0; SIGMA = 0.2;
PREC   = Accurate
LAECHG = .TRUE.
LCHARG = .TRUE.
这是输入文件,在计算 ...

我不知道你们计算怎么提交的,为啥作业结束会EXTRAFILE,是否有可能没有把这三个文件打包在内?不妨把提交作业的脚本也发来看看。
不一定挂在论坛,计算问题问题欢迎留言。
7楼2016-07-13 16:28:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lorna639 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 4/200 2026-11-02 00:50 by xzcv0Tehd8lQ
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-02 00:47 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 hKe7U1kad5Sg 2026-10-31 8/400 2026-11-02 00:37 by xzcv0Tehd8lQ
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-02 00:11 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 19/950 2026-11-02 00:02 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 23:55 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 15/750 2026-11-01 23:45 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-01 21:52 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 21:45 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 8/400 2026-11-01 21:35 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 21:26 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 4/200 2026-11-01 11:01 by BCdvBAEpt3zm
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
请填处理意见