首页时事民生政务教育文化科技财富体娱健康情感

旅行百科职场楼市企业乐活学术汽车时尚创业美食幽默美体文摘

VASP脚本使用手册 | 电荷密度计算脚本之电荷密度处理——求和、做差

2025-01-05 08:10 广东

【做计算找华算】理论计算助攻顶刊，50000+成功案例，全职海归技术团队、正版商业软件版权！

经费预存选华算科技：新年献礼，重磅加码，高至30%预存增值！

往期推荐

VASP脚本使用手册 | 获取力和能量收敛趋势、一键续算功能

VASP脚本使用手册 | POSCAR处理之等比缩放晶胞、获取中间态结构、输出径向分布

VASP脚本使用手册 | POSCAR处理之对比POSCAR文件的差异

VASP脚本使用手册 | 不同格式结构文件转换

VASP脚本使用手册 | 输出优化过程中原子运动轨迹及速度

VASP脚本使用手册 | NEB计算脚本系列之生成中间态结构、读取极值点能量及受力

VASP脚本使用手册 | NEB计算脚本系列之绘制势垒拟合线、导出轨迹图

VASP脚本使用手册 | NEB计算脚本系列之获取收敛信息、修改原子坐标

VASP脚本使用手册 | 态密度计算脚本之分解DOSCAR、分析频带中心及绘图

本期脚本介绍如何处理VASP计算得到的CHGCAR电荷文件，因此计算前需确认是否设置了LCHARG = .TRUE.。以Si单质自洽计算得到的CHGCAR为例，下文中涉及到的CHGCAR1与CHGCAR2为同一文件。注：该例只用于解释脚本用途及处理效果，不具备实际科研意义。

脚本028：chgsum.pl

介绍：多个CHGCAR文件求和，输出文件为CHGCAR_sum。CHGCAR_sum的计算公式为 (CHGCAR1*fact1+CHGCAR2+fact2)。默认情况下，fact1=fact2=1.0，因此输出是输入电荷密度文件的总和

用法：chgsum.pl (CHGCAR1) (CHGCAR2) (fact1) (fact2)

示例：

第一步：

当前窗口输入chgsum.pl CHGCAR1 CHGCAR2。

第二步：

ls查看。得到CHGCAR_sum文件。

相同isosurface level下电荷密度分布示意图如下：

脚本029：chgdiff.pl

介绍：产生电荷差值文件CHGCAR_diff文件

用法：输入chgdiff.pl (CHGCAR1) (CHGCAR2)

示例：

第一步：

当前工作窗口输入chgdiff.pl (CHGCAR1) (CHGCAR2)

第二步：

ls查看。得到CHGCAR_sum文件。

相同isosurface level电荷密度分布示意图如下：

由于篇幅有限，后续脚本干货将在本公众号（DFT计算参谋长）持续更新，关注公众号可及时获取最新教程。

这是一套与过渡态、电荷密度、态密度等性质相关的脚本包VTST (Transition State Tools for VASP)。

这140个脚本文件具体功能包括过渡态计算的初始路线设置和能垒的查看、电荷密度的求和与做差、态密度文件的拆解与作图，这些功能覆盖几乎所有的DFT计算内容，数据处理效率提升99%！

【做计算找华算】

🏅 我们提供专业的第一性原理、分子动力学、生物模拟、量子化学、机器学习、有限元仿真等代算服务。

🎯我们的理论计算服务，累计助力5️⃣0️⃣0️⃣0️⃣0️⃣➕篇科研成果，计算数据已发表在Nature & Science正刊及大子刊、JACS、Angew、PNAS、AM系列等国际顶刊。 👏👏👏

👉 点击阅读原文加我，探索更多优惠💖

DFT计算参谋长

华算科技朱参谋长公众号，十余年VASP重度使用经验，每天给大家带来期刊速览，学术新闻，欢迎关注！

最新文章

吸附催化计算不会怎么办?

0元算名额告急！DFT计算限时免单，不容错过！

VASP考虑范德瓦尔斯计算结构要点总结

材料结构、电子、吸附性质计算入门

新年福利翻倍！理论计算全场8折、预存高至30%增值、科研日历加赠，数量有限，先到先得！

Python辅助建模、计算模拟、MDAnalysis辅助分析及论文复现！

华算科技福利大放送！理论计算全场8折、预存高至30%增值，加赠科研日历！

0元算名额告急！DFT计算限时免单，不容错过！

模型大了算不动，小了被质疑

钙钛矿能带与态密度计算不会怎么办？

DFT计算限时免单，不容错过！

入门必备！344页PPT《基于DFT的第一性原理计算方法简介》

新年福利翻倍！理论计算全场8折、预存高至30%增值、科研日历加赠，数量有限，先到先得！

大家都在学的DFT建模与计算课程！

超权威VASP官方Workshop：维也纳大学开发组Dr. Martijn Marsman亲授！

0元算名额告急！DFT计算限时免单，不容错过！

磁性材料不会计算怎么办！

天降福利！6套VASP计算课程：态密度、能带、过渡态、扩散、电池、表面吸附、半导体缺陷、LAMMPS等！价值过万，免费拿走！

福利翻倍！理论计算全场8折、预存高至30%增值、科研日历加赠，数量有限，先到先得！

定个小目标！春节前拿下编程语言的新王者！

最新版！晶体材料结构文件CIF的查找和下载方法！

炸裂！居然免费算！DFT计算免单，名额有限，快抢！

各种过渡态计算原来这么简单！

VASP脚本使用手册 | 电荷密度计算脚本之电荷密度处理——求和、做差

福利翻倍！理论计算全场8折、预存高至30%增值、科研日历加赠，数量有限，先到先得！

磁性材料计算怎么学！

快冲！华算科技年终超级福利来袭！理论计算全场8折、高至30%预存增值，更有科研日历加赠！

理论计算答疑合集027：如果计算非晶的能带没有意义，那是不是计算非晶的态密度也没有意义？

炸裂！居然免费算！DFT计算免单，名额有限，快抢！

二维材料弄不明白怎么办！

这些功能居然能覆盖几乎所有的DFT计算内容！数据处理效率提升99%！

预存高至30%增值，科研福利超级重磅加码，全年无忧，狂发顶刊！

立刻学会钙钛矿能带与态密度计算！

理论计算答疑合集026：VASP结构优化是怎么判断收敛的？

炸裂！居然免费算！DFT计算免单，名额有限，快抢！

分子催化剂、单原子、双原子催化剂模型构建

理论计算答疑合集025：算2000个原子到5000个原子需要多少个服务器做集群？

DFT计算限时免单！不容错过！

正负极材料及异质结构电子结构怎么分析

VASP实例计算课程，送送送

DFT计算限时免单！不容错过！

分子催化剂、单原子、双原子催化剂模型构建渡态

理论计算答疑合集024：不同结构的优化计算，K点如何选取？

科研人的年终特惠！华算科技理论计算全场8折，还有科研日历免费领！

燃料电池怎么入门？

理论计算答疑合集023：结构优化的时候，面对能量收敛但是力不收敛的情况，这应该怎么办？

科研人的年终特惠！华算科技理论计算全场8折，还有科研日历免费领！

VASP磁性计算：非线性磁矩、GGA+U、磁矩调控、自旋轨道耦合、四态法、蒙特卡洛、紧束缚模型、贝里曲率、量子反常霍尔效应

朱老师讲VASP | 理论计算答疑合集022：请问在选择赝势时，有的元素有好几个赝势，应该怎么选择呢？

DFT计算限时免单！不容错过！

分类

原创标签

时事社会财经军事教育体育科技汽车科学房产搞笑综艺明星音乐动漫游戏时尚健康旅游美食生活摄影宠物职场育儿情感小说曲艺文化历史三农文学娱乐电影视频图片新闻宗教电视剧纪录片广告创意壁纸头像心灵鸡汤星座命理教育培训艺术文化金融财经健康医疗美妆时尚餐饮美食母婴育儿社会新闻工业农业时事政治星座占卜幽默笑话独立短篇连载作品文化历史科技互联网

发布位置

广东北京山东江苏河南浙江山西福建河北上海四川陕西湖南安徽湖北内蒙古江西云南广西甘肃辽宁黑龙江贵州新疆重庆吉林天津海南青海宁夏西藏香港澳门台湾美国加拿大澳大利亚日本新加坡英国西班牙新西兰韩国泰国法国德国意大利缅甸菲律宾马来西亚越南荷兰柬埔寨俄罗斯巴西智利卢森堡芬兰瑞典比利时瑞士土耳其斐济挪威朝鲜尼日利亚阿根廷匈牙利爱尔兰印度老挝葡萄牙乌克兰印度尼西亚哈萨克斯坦塔吉克斯坦希腊南非蒙古奥地利肯尼亚加纳丹麦津巴布韦埃及坦桑尼亚捷克阿联酋安哥拉