24小时热门版块排行榜

返回列表

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

haowenping

木虫 (著名写手)

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

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

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

回复此楼

» 猜你喜欢

投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复
自荐读博已经有7人回复
自然科学基金委宣布启动申请书“瘦身提质”行动已经有4人回复
求个博导看看已经有18人回复

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

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

查看全部散金贴

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

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

haowenping

木虫 (著名写手)

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

引用回帖:

Originally posted by 渔火江枫 at 2011-04-15 15:23:18:
如果是想得到表面吸附物种与基底原子的成键，需要三个CHGCAR：CHGCAR1 是吸附物种和基底优化后做单点得到，CHGCAR2是把CONCAR中基底原子删掉，并相应修改INCAR，POTCAR，只算单点得到，CHGCAR3，把CONTCAR中吸附 ...

你好，谢谢你的回答，我还想问下：也就是说，我先总体的做了一次静态自洽计算。得到了CHGCAR1，然后呢，我把吸附原子去掉只剩下基底原子了，此时，我是要改变POTCAR的，但是INCAR中是要设置NGXF，NGYF,NGZF参数的，这些参数应该怎么设置？有什么标准没？如果把他设置成为NGX，NGY,NGZ的2倍行不行？这样得到了CHGCAR2文件。第三步就是把基底原子都删除，只做吸附原子的静态计算，（如果我只吸附了一个原子，是不是就只用这一个原子的静态计算），此时的INCAR中，我是不是直接用第二步的INCAR文件中设置的参数？这样我就得到了CHGCAR3。我还想问下，应该怎么将这三个电荷作差呢？怎么才能得到你说的那个CHGDIFF文件呢？希望继续指教，谢谢！