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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】读CHGCAR里的数据
91/1返回列表
查看: 1229  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 jzg3xk1po8 2026-07-31 12/600 2026-08-01 17:50 by Kddn1pB6pAzX
[基金申请] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 jzg3xk1po8 2026-07-31 12/600 2026-08-01 17:50 by Kddn1pB6pAzX
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 jzg3xk1po8 2026-07-31 9/450 2026-08-01 17:43 by Kddn1pB6pAzX
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 jzg3xk1po8 2026-07-31 7/350 2026-08-01 17:42 by Kddn1pB6pAzX
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 jzg3xk1po8 2026-07-31 9/450 2026-08-01 17:40 by Kddn1pB6pAzX
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +16 低垂的野花 2026-07-31 22/1100 2026-08-01 16:39 by Cathyau
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 jzg3xk1po8 2026-07-31 9/450 2026-08-01 15:50 by y4XuShGyllvT
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 jzg3xk1po8 2026-07-31 6/300 2026-08-01 15:46 by y4XuShGyllvT
[基金申请] 微信指数没变化，科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[教师之家] 基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 8/400 2026-07-31 18:13 by beefly
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了，估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了？ +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了，我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
关闭
请填处理意见
关闭