| 查看: 6099 | 回复: 16 | ||||
[交流]
关于VASP计算得到CHGCAR文件的数据处理以及数据单位的
|
||||
|
大家好： 本人在处理vasp计算得到的CHGCAR文件时发现一个问题。我们从vasp的手册上可以知道，CHGCAR里面的数据是指 p(r)*Vcell 前面的P(r)是指电荷密度，单位是e / A^3, 而后面的Vcell是指自己计算体系的原胞体积，当然Vcell的单位就是A^3了——前面的这一段话，不要质疑，是我通过积分验证了的。其实就是将CHGCAR里面的所有数据除以CHGCAR中网格的数目 （NG(X,y,Z)F —— p(r)*Vcell /(NGXF*NGYF*NGZF) 然后再求和这样就等于体系的价电子数。 那么我想问的是大家在处理CHGCAR的数据时，是直接将CHGCAR里面的数据读出来转换成MS可视化的数据，还是把CHGCAR中的数据除以了一个Vcell 然后再转换成MS可视化的数据？这两种的单位分别是什么呢？我们作图的时候，有总电荷密度，差分电荷密度，自旋电荷密度，那么分别在处理数据时采用前面的那种方法？ 组里面前任编的Fortran，总电荷密度和自旋电荷密度是直接读取然后就用了，但是在差分电荷密度的时候，又除以了原胞体积Vcell, 很不统一。然后去看他们的那些文献，发现当涉及到介绍这些图片时，单位都是很含糊，有的是 eV/（au）^3, 也有说成这样的eV / A^3。 这里面的单位在作图时切一个面，面电荷 和 直接使用体电荷密度的 一致吗？ 不一致的话区别在哪里？ 求高人指点！！！ |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理相关文档 |
» 本帖已获得的红花（最新10朵）
refors» 猜你喜欢
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有4人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
什么时候开奖？
已经有10人回复
-
2026国自然函评费到账
已经有16人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
让我中一个面上吧！
已经有14人回复
-
今天放榜吗？
已经有16人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP中如何创建 INCAR , POSCAR, KPOINTS和POTCAR文件
已经有7人回复
- VASP计算得到的OUTCAR为什么没有给出费米能级？
已经有13人回复
- 求助关于用vasp计算二维石墨烯纳米带所得到的能带结果与报道对不上
已经有7人回复
- vesta 读取vasp CHGCAR 的单位
已经有7人回复
- 关于vasp收敛和优化
已经有24人回复
- VASP计算能带时是否读取WAVECAR
已经有16人回复
- 求助关于vasp磁矩设置的说明书的详细理解
已经有12人回复
- VASP计算：如何处理波函数占用空间太大的问题
已经有5人回复
- vasp中CHGCAR是PS电荷还是AE电荷
已经有6人回复
- 【任务征集】VASP跑分子动力学数据处理方法
已经有20人回复
- 【其他】vasp_band.exe处理的数据怎么最后不一样啊？
已经有15人回复
- 【活动】vasp计算分子动力学讨论合集及专贴（筹划中）
已经有33人回复
- 【求助】关于vasp态密度计算
已经有10人回复
- 【求助】vasp中关于电场单位的问题
已经有12人回复
- 【neweroica个人文集】操作VASP的CHGCAR文件的程序：chgcarman（v0.0）
已经有43人回复
» 抢金币啦！回帖就可以得到:
-
散金祈福
+5/3985
-
散金币
+1/941
-
散金祈福
+3/792
-
每年这时候我都再仔细读一遍《范进中举》
+1/473
-
南方医科大学|生物材料与医学工程团队诚聘科研助理
+2/256
-
五大联赛看好哪支球队
+5/200
-
诚征另一半
+1/172
-
上海交通大学张航课题组招聘博士后（电化学储能）
+1/89
-
中山大学化学工程与技术学院欧阳钢锋/吴鑫课题组博士生（2027年入学）/博士后招聘启事
+1/34
-
瑞典-博后-固态电池-粘结剂-锂离子电池方向
+1/33
-
瑞典-博后-固态电池-粘结剂-锂离子电池方向
+1/27
-
中科院大连化学物理研究所 招聘 催化剂研发方向 科研助理 2名
+1/5
-
南京大学 智能驱动与感知材料实验室 诚招申请考核博士生/科研助理/博士后
+1/4
-
新审查指南施行半年，AI方向的专利申请是不是更看"技术味"了？
+1/2
-
清华大学环境课题组招聘环境工作专业客座生研究生2-3名（应用导向）
+1/2
-
出全新蔡司 Plan?APOCHROMAT 20×/0.8 平场复消色差物镜
+1/2
-
清华大学环境课题组招聘环境工作专业客座生研究生2-3名（应用导向）
+1/1
-
找女友
+1/1
-
【产品测评分享】体外转录产量不理想、dsRNA 副产物高？T7 RNA 聚合酶实测体验
+1/1
-
2026年国自然
+1/1
9楼2012-06-20 22:43:26
10楼2012-06-21 10:35:12
4楼2012-06-20 16:30:15
8楼2012-06-20 17:03:13
11楼2012-06-23 11:28:51
12楼2012-06-27 21:34:13
14楼2014-05-29 19:46:20
★
lsh1986914(金币+1): 谢谢参与
lsh1986914(金币+1): 谢谢参与
|
请问下lz，为什么我对CHGCAR 求和得不到体系的价电子总数？ In CHGCAR file, the format is: WRITE(IU,FORM) (((C(NX,NY,NZ),NX=1,NGXC),NY=1,NGYZ),NZ=1,NGZC) according to the manual. However, I made a summation of the total charge is largely different from the electrons in my system. My script is: c--- read the charge density open(4,file='CHGCAR',status='old') read(4,*) (((C(NX,NY,NZ),NX=1,NGXC),NY=1,NGYC),NZ=1,NGZC) c---- I have opened the CHGCAR file and on each axis, there are 160 grids. do Nx=1,160 do Ny=1,160 do Nz=1,160 totalcharge = totalcharge + C(Nx,Ny,Nz)/volume enddo enddo enddo open(15,file='Charge_total.dat') write(15,*) totalcharge My result is 40959.9999. But in fact, my system only has two atoms. It is a CO module. I am very confused about my result but I cannot figure it out. I just read the data and make a summation. Could someone give me a hint or suggestion? Thanks! |
15楼2017-09-20 08:35:08
16楼2019-01-09 14:50:54
17楼2019-02-26 11:25:07
简单回复
shanly2楼
2012-06-20 16:06
回复
lsh1986914(金币+1): 谢谢参与
2012-06-20 16:27
回复
lsh1986914(金币+1): 谢谢参与
2012-06-20 16:31
回复
lsh1986914(金币+1): 谢谢参与
2012-06-20 16:39
回复
lsh1986914(金币+1): 谢谢参与
hr6307楼
2012-06-20 16:52
回复
lsh1986914(金币+1): 谢谢参与
showtimexh13楼
2013-09-25 21:01
回复
