应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】读CHGCAR里的数据
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 1132  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

yunml

主管区长

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！

[交流] 【求助】读CHGCAR里的数据

虽然用vesta能很好的看电荷密度，但是不能得到几个相减的电荷密度差。网上有小软件做这个，但我的老是出问题，所以特别想自己拿最原始的数据，自己动手在origin做。

问题是chgcar里面的数据根本读不同，有没有哪位高人能指点下？ 不胜感激。
回复此楼

» 猜你喜欢
一以贯之
1楼 2009-12-04 10:47:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TIGERYZZ

兑换贵宾

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+3,VIP+0):感谢指点 12-4 13:46
chgcar首先输出的是结构信息，然后输出FFT-grid的维度，即NGX，NGY，NGZ，然后输出网格中每一个点的电荷密度，输出的格式在vasp手册里有提到：
WRITE(IU,FORM) (((C(NX,NY,NZ),NX=1,NGXC),NY=1,NGYZ),NZ=1,NGZC)
赞一下(11人) 回复此楼
2楼2009-12-04 13:24:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

灯塔守望者

主管区长

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+1,VIP+0):鼓励交流 12-5 16:37
OK，下一步我将在vaspkit中增加提取差分电荷密度数据！
赞一下(11人) 回复此楼
3楼2009-12-04 17:50:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

专家顾问

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
自己动手编程吧
回复此楼
4楼2009-12-04 17:51:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yunml

版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
Originally posted by TIGERYZZ at 2009-12-4 13:24:
chgcar首先输出的是结构信息，然后输出FFT-grid的维度，即NGX，NGY，NGZ，然后输出网格中每一个点的电荷密度，输出的格式在vasp手册里有提到：
WRITE(IU,FORM) (((C(NX,NY,NZ),NX=1,NGXC),NY=1,NGYZ),NZ=1,NGZC)

能不能解释一下"WRITE(IU,FORM) (((C(NX,NY,NZ),NX=1,NGXC),NY=1,NGYZ),NZ=1,NGZC) "?  因为我看chgcar里面电荷密度部分每行有五个数据，按照（x,y,z,density）理解好像不对啊。

还有，上面的语句是Fortran写的吗？  多谢了！

CHGCAR 如下：
Fe-CO                                    
   1.00000000000000     
     0.000000    5.732800    0.000000
     0.000000    0.000000    5.732800
    17.166000    0.000000    0.000000
   1   1  24
Direct
  0.500000  0.464705  0.451881
  0.500000  0.641670  0.501367
  0.000000  0.000000  0.000000
  0.000000  0.500000  0.000000
  0.500000  0.000000  0.000000
...

   48   48  144
0.11710375410E+05 0.11082654420E+05 0.93465338201E+04 0.69185564327E+04 0.43845496481E+04
0.23036430807E+04 0.98892296509E+03 0.40174174331E+03 0.23369168341E+03 0.16858018791E+03
0.13671233690E+03 0.12337598318E+03 0.11990845078E+03 0.12335148597E+03 0.13665683923E+03
0.16846542258E+03 0.23345974810E+03 0.40125016805E+03 0.98777940755E+03
...

augmentation occupancies  25 138
-0.8339124E-02  0.1781598E-03 -0.2282241E-03 -0.1983454E-03 -0.3674109E-04
  0.4917278E-03 -0.1912815E-01  0.8767631E-03 -0.7149171E-03 -0.8271628E-03
-0.2565947E-03  0.2019992E-02  0.9772558E-02 -0.5522550E-02  0.3923657E-01
...
augmentation occupancies  26 138
-0.9033073E-02  0.1112419E-03  0.3308886E-03 -0.2789434E-03  0.9789807E-04
  0.5204789E-03 -0.1939311E-01  0.4184915E-03  0.1347132E-02 -0.1239935E-02
  0.2499834E-03  0.2262646E-02  0.5939990E-02  0.3275988E-02  0.3513974E-01
  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
...
赞一下 回复此楼
一以贯之
5楼2009-12-05 09:24:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yunml

专家顾问

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
Originally posted by 灯塔守望者 at 2009-12-4 17:50:
OK，下一步我将在vaspkit中增加提取差分电荷密度数据！

问下，我想用Fe吸附CO后的减去单独优化的Fe表面的，再减去单独CO的可以吗？

网上那个小软件似乎必须得结构参数一样才能减。

非常感谢！
赞一下 回复此楼
一以贯之
6楼2009-12-05 09:30:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

灯塔守望者

管理员

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+3,VIP+0):谢谢指点 12-5 16:37
可以，只要保证前后晶胞结构参数和K点不变，所产生的电荷密度维数是相同的，直接减就行。
WRITE(IU,FORM) (((C(NX,NY,NZ),NX=1,NGXC),NY=1,NGYZ),NZ=1,NGZC)等价于
for NZ=1,NGZC
    for Ny=1,NGYZ
         for NX,=1, NGXC
         print   ***

看看fortran书就知道了，有这时间还不如自己查资料搞定。
赞一下(12人) 回复此楼
7楼2009-12-05 10:50:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yunml

兑换贵宾

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
引用回帖:
Originally posted by 灯塔守望者 at 2009-12-5 10:50:
可以，只要保证前后晶胞结构参数和K点不变，所产生的电荷密度维数是相同的，直接减就行。
WRITE(IU,FORM) (((C(NX,NY,NZ),NX=1,NGXC),NY=1,NGYZ),NZ=1,NGZC)等价于
for NZ=1,NGZC
    for Ny=1,NGYZ
         ...

谢谢回复. 不过你的解释也没看懂. 以前没接触过Fortran,马上就开始学.

我不明白你说的"晶胞结构参数不变"的意思. 我前面说的(slab+CO)-CO,减号前后构型是一样的吗？
赞一下 回复此楼
一以贯之
8楼2009-12-05 12:34:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

灯塔守望者

实习版主

优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！优秀！！有木有！！！
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+1,VIP+0):谢谢指点 12-5 16:38
就是晶胞大小不变。计算CO是把铁原子去掉！
赞一下(11人) 回复此楼
9楼2009-12-05 13:30:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yunml 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 285求调剂 +4 AZMK 2026-04-01 4/200 2026-04-01 12:15 by liumengping
[考研] 生物与医药考研调剂 +3 铁憨憨123425 2026-03-31 3/150 2026-04-01 10:53 by judy105
[考研] 301求调剂 +8 axibli 2026-04-01 8/400 2026-04-01 09:51 by 我的船我的海
[考研] 350求调剂 +6 阿佳～ 2026-03-31 6/300 2026-04-01 08:49 by 不吃魚的貓
[考研] 考研生物与医药调剂 +7 铁憨憨123425 2026-03-31 7/350 2026-04-01 08:45 by JourneyLucky
[考研] 080200学硕，机械工程专业277分，求带走！ +4 瓶子PZ 2026-03-31 4/200 2026-03-31 20:16 by vgtyfty
[考研] 375求调剂 +7 雨夏整夜 2026-03-29 7/350 2026-03-31 18:52 by xhai2011
[考研] 张芳铭-中国农业大学-环境工程专硕-298 +9 手机用户 2026-03-26 9/450 2026-03-31 18:09 by 544594351
[考研] 求调剂 生物学 377分 +6 zzll03 2026-03-31 6/300 2026-03-31 17:33 by 唐沐儿
[考研] 330分求调剂 +6 qzenlc 2026-03-29 6/300 2026-03-31 16:44 by Wang200018
[考研] 一志愿南昌大学324求调剂 +6 hanamiko 2026-03-29 6/300 2026-03-31 16:35 by hypershenger
[基金申请] 面上5B能上会吗？ +8 redcom 2026-03-29 8/400 2026-03-31 15:53 by niuailing
[考研] 323分 食品与营养调剂 +3 嘿ooo 2026-03-31 3/150 2026-03-31 09:38 by longlotian
[考研] 一志愿华中师范化学332分求调剂 +3 Lyy930824@ 2026-03-29 3/150 2026-03-30 20:15 by DHUSHUAI
[考研] 考研调剂 +5 Sanmu-124 2026-03-26 5/250 2026-03-30 13:36 by chemdavid
[考研] 一志愿北京理工大学本科211材料工程294求调剂 +8 mikasa的围巾 2026-03-28 8/400 2026-03-29 12:48 by 无际的草原
[考研] 一志愿211院校 344分 东北农业大学生物学学硕，求调剂 +5 丶风雪夜归人丶 2026-03-26 8/400 2026-03-27 19:22 by 丶风雪夜归人丶
[考研] 一志愿上海理工能源动力（085800）310分求调剂 +3 zhangmingc 2026-03-27 4/200 2026-03-27 19:01 by 给你你注意休息
[考研] 292求调剂 +4 求求了收下我吧� 2026-03-26 4/200 2026-03-27 10:37 by zhshch
[考研] 材料与化工304求B区调剂 +3 邱gl 2026-03-25 3/150 2026-03-25 19:03 by Ainin_
信息提示
关闭
请填处理意见
关闭