24小时热门版块排行榜    

查看: 1562  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lgxyz

铜虫 (正式写手)

[求助] 请教差分电荷

各位,我想做这样的一个差分电荷:Δρ=ρ总-ρ1-ρ2-ρ3,
ρ总是总电荷;ρ1,ρ2,ρ3分别是金属表面,吸附物1和吸附物2的电荷
我用chgdiff.pl脚本来做,chgdiff.pl只能求两个电荷之差得到CHGCAR_diff,而CHGCAR_diff的格式和CHGCAR不一样,就不能用chgdiff.pl脚本处理了。
不知道大家有什么好的求多个电荷差的脚本或现成的程序,谢谢赐教!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)

引用回帖:
3楼: Originally posted by fanchen021 at 2011-12-20 16:57:09:
CHGCAR_diff的格式和CHGCAR的格式是完全一样的

谢谢!
好像每个CHGCAR里还有下面这个东西,不知啥意思,CHGCAR_diff里没有这些,请指教下:
augmentation occupancies   1  15
  0.4793829E+00 -0.3411239E+00  0.9496256E-01 -0.1884780E+00 -0.5617611E-01
  0.8677252E-01 -0.3696490E-01  0.7850851E-01  0.2418509E-01  0.3642196E-01
-0.2488471E-02 -0.1165034E-01  0.1220837E-01  0.8700472E-02 -0.2518154E-02
augmentation occupancies   2  15
  0.4691445E+00 -0.3404651E+00 -0.6024096E-01  0.1146301E+00 -0.1735776E+00
  0.8738012E-01  0.2632801E-01 -0.4503226E-01  0.7270577E-01  0.3668362E-01
  0.9227440E-02 -0.3974533E-02 -0.1543592E-02 -0.1452878E-01  0.9781519E-02
augmentation occupancies   3  15
  0.4772481E+00 -0.3406651E+00 -0.1733294E+00 -0.3171649E-01  0.1301968E+00
  0.8690734E-01  0.7226215E-01  0.1472136E-01 -0.5184649E-01  0.3654300E-01
-0.1519754E-01  0.5406170E-02 -0.9723277E-02 -0.2009987E-02 -0.5033813E-02
augmentation occupancies   4  15
...................
...............
..................
4楼2011-12-21 11:33:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

fanchen021

铁杆木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
fzx2008(金币+1): 谢谢交流 2011-12-20 21:23:23
lgxyz(金币+5): 2011-12-21 15:02:14
CHGCAR_diff的格式和CHGCAR的格式是完全一样的
3楼2011-12-20 16:57:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)

【答案】应助回帖


fzx2008(金币+1): 谢谢指教 2011-12-21 14:43:40
那部分内容不影响你做差分。 那个脚本仍然可以用。
5楼2011-12-21 14:35:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)

引用回帖:
5楼: Originally posted by fanchen021 at 2011-12-21 14:35:37:
那部分内容不影响你做差分。 那个脚本仍然可以用。

谢谢!还有个问题:
我们用chgdiff.pl CHGCAR1 CHGCAR2,是想用CHGCAR1减去CHGCAR2,可是chgdiff.pl里写的是用CHGCAR2减去CHGCAR1,请给些建议。
6楼2011-12-21 15:01:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302材料工程求调剂 +4 Doleres 2026-03-01 5/250 2026-03-01 11:52 by liqiongjy
[考研] 272求调剂 +5 材紫有化 2026-02-28 5/250 2026-03-01 11:51 by gaoxiaoniuma
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +3 虫虫虫虫虫7 2026-03-01 6/300 2026-03-01 11:50 by gaoxiaoniuma
[考研] 290求调剂 +7 材料专硕调剂; 2026-02-28 8/400 2026-03-01 11:45 by 王伟要上岸啊
[考研] 材料化工调剂 +6 今夏不夏 2026-03-01 6/300 2026-03-01 11:28 by gaoxiaoniuma
[考研] 高分子化学与物理调剂 +5 好好好1233 2026-02-28 9/450 2026-03-01 10:59 by fengyu211
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 10/500 2026-03-01 10:02 by 科研狗111
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[考研] 284求调剂 +5 天下熯 2026-02-28 5/250 2026-03-01 08:49 by ms629
[考研] 298求调剂 +5 axyz3 2026-02-28 5/250 2026-03-01 06:45 by 刘兵
[考研] 材料调剂 +4 爱擦汗的可乐冰 2026-02-28 4/200 2026-03-01 00:38 by 猫猫球alter
[基金申请] 面上模板改不了页边距吧? +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
[考研] 295求调剂 +5 19171856320 2026-02-28 5/250 2026-02-28 21:39 by gaoxiaoniuma
[考博] 26申博 +4 想申博! 2026-02-26 4/200 2026-02-28 21:37 by limorning
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[基金申请] 面上可以超过30页吧? +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
信息提示
请填处理意见