| 查看: 1670 | 回复: 1 | |||
[交流]
chgdiff.pl 处理电荷差分密度 已有1人参与
|
|
本来想用chgdiff.pl处理电荷差分密度, 但处理三个文件时有个小bug的, 浪费半天时间 OMG 写下来供借鉴 1. 因为原子信息保留了file1的信息, 应该倒过来, file1-file2 2. printf OUT输出前加个空格, 保持与VASP格式一致 修改: 81c81 < $line1[$j] = $line1[$j]-$line2[$j]; --- > $line1[$j] = $line2[$j]-$line1[$j]; 83c83 < printf OUT " %18.11E %18.11E %18.11E %18.11E %18.11E\n",$line1[1],$line1[2],$line1[3],$line1[4],$line1[5]; --- > printf OUT "%18.11E %18.11E %18.11E %18.11E %18.11E\n",$line1[1],$line1[2],$line1[3],$line1[4],$line1[5]; chgdiff.pl CHGCAR1 CHGCAR2 mv CHGCAR_diff tmp chgdiff.pl tmp CHGCAR3 得到差分密度:CHGCAR1-CHGCAR2-CHGCAR3 用网上找的fortran处理得到plane-averaged difference electronic density, 注意删除原子符号那行 齐活! |
回复此楼» 本帖已获得的红花（最新10朵）
朱全喜» 猜你喜欢
浙江师范大学是怎么坑我的
已经有12人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
浙江师范大学是怎么坑我的
已经有4人回复
-
闲聊
已经有6人回复
-
江苏省自然基金 什么时候出结果
已经有7人回复
-
求教各位大神：2026教育部人文社科青年基金项目何时公示呀？
已经有6人回复
-
面上有专家说实验设备不是我们单位的
已经有8人回复
-
科研求助
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 急需会对chgdiff.pl 脚本对电荷密度做差的高手帮忙！
已经有4人回复
- 差分电荷密度图的疑惑
已经有41人回复
- CHGCAR的文件含义
已经有7人回复
- P4vasp画PDOS，请教一下，谢谢。
已经有14人回复
- 如何用xcrysden画电荷密度图
已经有8人回复
- vasp如何计算差分电荷密度
已经有7人回复
- 差分电荷密度的计算
已经有11人回复
朱全喜
新虫 (小有名气)
- 应助: 19 (小学生)
- 金币: 2977.9
- 红花: 12
- 帖子: 205
- 在线: 132.6小时
- 虫号: 1972785
- 注册: 2012-09-04
- 性别: GG
- 专业: 理论和计算化学
2楼2015-05-27 16:28:40