| 查看: 1526 | 回复: 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朵）
朱全喜» 猜你喜欢
球磨粉体时遇到了大的问题，请指教！
已经有15人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有5人回复
-
情人节自我反思：在爱情中有过遗憾吗？
已经有5人回复
-
江汉大学解明教授课题组招博士研究生/博士后
已经有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