yhli0906的个人博客分享 http://blog.sciencenet.cn/u/yhli0906

博文

分享一个算平均电势和面内积分电荷密度的程序

已有 6316 次阅读 2019-5-7 11:08 |个人分类:计算机软件|系统分类:科研笔记

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


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

图片.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://blog.sciencenet.cn/blog-2909108-1177586.html

上一篇:巧用宏处理器生成输入文件
下一篇:分享一个无耻的抢资源脚本
收藏 IP: 121.249.15.*| 热度|

4 张宇 刘鹏飞 王安良 唐刚

该博文允许注册用户评论 请点击登录 评论 (3 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-3-29 04:24

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部