24小时热门版块排行榜    

查看: 1690  |  回复: 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 XlxB7q11Bc3x 2026-09-01 3/150 2026-09-01 21:52 by Y3OMSsg8u9zj
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
信息提示
请填处理意见