| 查看: 4687 | 回复: 16 | ||||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 lsh1986914 的 2 个金币 ，回帖就立即获得 1 个金币，每人有 1 次机会 | ||||
[交流]
关于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» 猜你喜欢
申博雅思求助
已经有7人回复
-
专科院校申请国自然
已经有18人回复
-
厅级项目出校却没中
已经有20人回复
-
博后网上办公系统维护
已经有6人回复
-
基金申请书名称有变化
已经有7人回复
-
关于读博感觉自己很抓马
已经有21人回复
-
高校辞职，要求赔偿，这到底合不合理
已经有34人回复
-
考博英语
已经有8人回复
-
博士该不该读？
已经有5人回复
-
申博好难
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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/318
-
东莞理工学院金属功能材料团队招科研助理1名
+1/88
-
西安“石油地质研究实验室”18791013812田老师
+1/83
-
哈工大（深圳）材料学院张嘉恒课题组（青千）—— 招聘科研助理
+1/39
-
中山大学孙逸仙纪念医院彭丽课题组实验和生信岗位招聘公告
+1/32
-
深圳技术大学招聘化学专业副研究员
+1/32
-
感慨
+1/29
-
某大型新能源企业诚招【正极材料研发专家/总监】！！
+1/23
-
四川大学生物材料方向招聘科研助理
+1/14
-
中科院纳米生物效应与安全性实验室2024年招收2-3名联合培养/客座研究生/博士后
+1/14
-
四川大学华西医院罕见病研究院招聘启事
+1/12
-
江苏省农科院植保所小麦病害防控团队诚招博士后
+1/9
-
湖南大学半导体学院招收直博生
+1/7
-
中科院半导体所二维量子材料方向招聘科研助理和博士后
+1/6
-
92年离异男，真诚找对象，坐标北京、河北
+1/6
-
纳米氧化镍 热敏性材料
+1/6
-
天津大学浙江研究院（宁波）诚聘青年研究员/博士后/研发工程师
+1/5
-
西工大新能源智能电驱动团队招收2024级博士后、博士、硕士
+1/3
-
知名化工外企招聘胶黏剂研发博士！高分子，聚合，化学方向均可！可看应届！
+1/3
-
NUS_化学科学硕士项目 MSc in Chemical Sciences, 2025年8月入学申请现已启动!
+1/2
+5/318
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
回复
