24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

haowenping

木虫 (著名写手)

应助: 6 (幼儿园)
金币: 2330.2
帖子: 1104
在线: 244.5小时
虫号: 1124850

[交流] 【求助】求助差分电荷密度的脚本！

各位大虾，本人最近一直想找一个处理查分电荷密度的脚本，有的希望大家发一个给我。顺便说一下我的做差分电荷密度的方法不知道这样是不是行：我做的也是吸附，首先我做了一次静态计算，里面没有设置NGFX,NGFY,NGFZ参数。下一步我就是把这四个输入文件放到一个文件夹下面，只是修改了其中的INCAR的ICHARG=12，同样并没有设置NGFX,NGFY,NGFZ参数，这样又得到了一个CHACAR。不知道这样可不可以希望大家给点建议，很是感激了，我也很急。如果我这个方法是对的希望，有人能发一个处理这个的脚本，我的邮箱是hwp22651304@163.com。收到后重金感谢！

回复此楼

» 猜你喜欢

论文终于录用啦！满足毕业条件了已经有16人回复
求个博导看看已经有19人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-04-13 08:51:38

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

渔火江枫

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 707.7
帖子: 239
在线: 74.5小时
虫号: 340411

★
haowenping(金币+10): 谢谢！ 2011-04-15 15:45:26
franch(金币+1): 谢谢指点 2011-04-15 20:45:22

引用回帖:

Originally posted by haowenping at 2011-04-15 14:54:07:
你好，非常感谢你！真的，我做的是表面吸附，请你建议下，我应该用那个做的比较好？谢谢！

如果是想得到表面吸附物种与基底原子的成键，需要三个CHGCAR：CHGCAR1 是吸附物种和基底优化后做单点得到，CHGCAR2是把CONCAR中基底原子删掉，并相应修改INCAR，POTCAR，只算单点得到，CHGCAR3，把CONTCAR中吸附分子删掉，同样修改INCAR，POTCAR，再做单点得到（注意：后面两个的INCAR 中要增加NGFX,NGFY,NGFZ，与第一个计算保持一样）。然后，将三个CHGCAR做差得到CHGDIFF。最后，处理CHGDIFF有多种方法，最简单的是用VESTA（这个没有教程，自己摸索），第二个方法可以用lev00，配合origin做图，lev00教程里有这个方法的说明。第三个方法用emuch上流传的MS读取charge density的方法做（emuch搜索）。