24小时热门版块排行榜    

查看: 2151  |  回复: 6

xiaoluoj

金虫 (正式写手)

[求助] vasp手册中关于CHGCAR的说明怎么理解

This file contains the lattice vectors, atomic coordinates, the total charge density multiplied by the volume ρ(r)∗Vcell
on the fine FFT-grid (NG(X,Y,Z)F)
, and the PAW one-center occupancies. CHGCAR can be used to restart VASP
from an existing charge density, for visualisation the CHG file should be used, since the PAW-one centre occupancies
are difficult to parse.
  1.红色字体部分怎么理解? 假设对于一个280个原子的超胞,体积为V,FFT为200*200*120,价电子数1280。那么 ρ(r)乘以的是超胞体积还是每个格点所占的体积(即V/200*200*120)?
2. ρ(r)∗体积不就是电子数了吗?但是CHG中的是如此之大,最大的都达到二千多了,但体系的总的电子数目才1280个,这如何理解?
————————————————————————————
我的部分CHG
200  200  192
  20.506      20.610      21.071      22.262      24.571      28.168      32.917      38.583      45.244      53.605   
  64.975      80.907      102.61      130.44      163.92      202.79      248.41      304.75      378.11      475.60   
  602.94      762.33      951.61      1165.3      1396.7      1640.2      1889.7      2137.6      2371.6      2575.2   
  2729.5      2817.6      2828.7      2761.3      2623.7      2431.6      2204.1      1958.6      1708.3      1462.0   
  1226.1      1007.4      813.33      650.06      520.47      423.09      353.14      304.76      272.90      254.46   
  248.15      253.72      271.39      302.30      349.30      417.14      511.50      637.14      795.77      984.82   
  1198.2      1428.2      1667.9      1910.3      2146.7      2364.6      2547.6      2678.3      2742.1      2731.2   
  2647.2      2500.0      2305.4      2080.8      1841.8      1599.9      1363.0      1137.3      928.95      744.73   
  589.95      466.50      371.94      300.54      245.37      200.43      162.01      128.93      101.53      80.216   
  64.588      53.405      45.128      38.485      32.818      28.076      24.502      22.224      21.058      20.608   
  20.506      20.610      21.071      22.262      24.571      28.168      32.917      38.583      45.244      53.605   
  64.975      80.907      102.61      130.44      163.92      202.79      248.41      304.75      378.11      475.60   
  602.94      762.33      951.61      1165.3      1396.7      1640.2      1889.7      2137.6      2371.6      2575.2   
  2729.5      2817.6      2828.7      2761.3      2623.7      2431.6      2204.1      1958.6      1708.3      1462.0   
  1226.1      1007.4      813.33      650.06      520.47      423.09      353.14      304.76      272.90      254.46   
  248.15      253.72      271.39      302.30      349.30      417.14      511.50      637.14      795.77      984.82   
  1198.2      1428.2      1667.9      1910.3      2146.7      2364.6      2547.6      2678.3      2742.1      2731.2   
  2647.2      2500.0      2305.4      2080.8      1841.8      1599.9      1363.0      1137.3      928.95      744.73
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
2楼2012-11-19 08:59:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoluoj

金虫 (正式写手)

我顶,千万别沉了
3楼2012-11-19 11:48:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiaozaoququ

木虫 (著名写手)

还真没注意手册上的这句话。期待大牛出现解释
再牛逼的选择也需要更傻逼的坚持
4楼2012-11-20 21:01:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)

我觉得ρ(r)是超胞内向量为r处的电荷密度,所以乘以体积不会是总得电子
timeflies..
5楼2012-11-21 15:23:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoluoj

金虫 (正式写手)

6楼2013-07-31 14:40:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)

★ ★
franch: 金币+2, 谢谢回帖交流,, 2013-08-01 00:40:34
the total charge density multiplied by the volume ρ(r)∗Vcell on the fine FFT-grid (NG(X,Y,Z)F)

===============================================================

总电荷密度n(r)是个函数,存储的是这个函数的傅立叶变换系数。

ρ(r)∗体积不就是电子数了吗?
================================================================
应该是n(r)在全空间的积分(可能还要乘以体积)等于电子数。
7楼2013-07-31 15:48:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaoluoj 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上复盘~欢迎讨论 +4 晴天加油 2026-08-26 5/250 2026-08-26 21:27 by 木风9012
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
信息提示
请填处理意见