| 查看: 1624 | 回复: 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朵）
朱全喜» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有12人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有23人回复
-
面上提前没消息，有中的吗
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有13人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有16人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有14人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有17人回复
-
2026年国自然面上资助率
已经有21人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有10人回复
-
2027年申博
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 急需会对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