李云海
分享一个算平均电势和面内积分电荷密度的程序
2019-5-7 11:08
阅读:7215

下载源码后解压,直接编译即可。


运行时需要输入这个几个参数:

图片.png

第一个是包含电荷密度或者静电势原始数据的文件名称。第二个是输出文件名称。第三个是针对原始数据的缩放因子,因为cube格式使用原子单位制,需要转换为国际单位制(立方埃或eV)则指定该因子,不需转换时则设置为1。最后一个是积分算法,1为梯形法,2为简单求和,取1即可。


输出文件中有三列数据:第一列是z轴坐标,单位为bohr;第二列是积分后的格点数据;第三列为平均后的格点数据。做面内积分电荷密度采用第一、第二列数据,做平均电势取第一、第三列数据。运行效果如下(以黑磷加电场后,面内积分电荷密度为例):

rho_3d.jpg

rho.png

程序目前只支持cube格式的数据,且只能在xOy平面内积分。VASP CHGCAR可以用VTST附带的chg2cube.pl转换成cube格式。


potslice.tar.bz2

转载本文请联系原作者获取授权,同时请注明本文来自李云海科学网博客。

链接地址:https://wap.sciencenet.cn/blog-2909108-1177586.html?mobile=1

收藏

分享到:

当前推荐数:4
推荐到博客首页
网友评论3 条评论
确定删除指定的回复吗?
确定删除本博文吗?