| 查看: 21614 | 回复: 142 | ||||||||||||||||||||
| 【奖励】 本帖被评价116次，作者gaottgougou增加金币 90.6501 个 | ||||||||||||||||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||||||||||||||||||
[资源]
【分享】将ELFCAR、CHGCAR转化成MS阅读可读的fortran代码
|
||||||||||||||||||||
|
在程序中只要修改读取格式就能转换不同的vasp生成文件（ELFCAR或CHGCAR） http://d.namipan.com/d/chgcar2gr ... 23622d46c678c080000 [ Last edited by zzy870720z on 2010-9-9 at 20:32 ] youzhizhe版主补充本地下载[ Last edited by youzhizhe on 2011-10-6 at 22:49 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 纳米技术与能源及模拟 | 精华网帖收集 | 第一性原理计算资源总汇 | VASP and MS |
| 研究生资料集 | VASP | VASP | VASP使用技巧 |
| 第一性原理计算-tg | vasp | first principle | VASP算例 |
| VASP资料 | @个人收集@ | VASP计算 | 爪机党 |
| Gaussian软件 | VASP相关 |
» 本帖已获得的红花（最新10朵）
hookah» 猜你喜欢
建议基金发布提前给出明确的时间点
已经有14人回复
-
范进中举一文的中心思想
已经有5人回复
-
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- GC-MS定量测定脂肪酸时标品的问题
已经有10人回复
- 纳米薄膜太阳能电池转化效率达8.1%
已经有11人回复
- Intel Visual Fortran应用程序开发
已经有259人回复
- 有关用HPLC-MS检测食物中农药残留可以投的英文杂志（影响因子低于1 一下）
已经有17人回复
- 【活动】Show 脚本，赢奖励
已经有23人回复
- 大家来帮我看看这两个INCAR 文件有什么区别哦。。为啥计算结果不一样呢
已经有4人回复
- 【原创】CHGCAR转化成MS可读文件修正
已经有18人回复
|
实际上ELFCAR和CHGCAR是有差别的。在两点上： （1）网格划分的控制参数不同。CHGCAR是NGXF*NGYF*NGZF;ELGCAR是NGX*NGY*NGZ. (2)每个数字的格式不同。 另外，你的程序中这段： IF(NSPIN==1) THEN ALLOCATE(work1(nx,ny,nz)) READ(16,vasp_form) (((work1(i,j,k),i=1,nx),j=1,ny),k=1,nz) WRITE(15,dmol_form) (((work1(i,j,k),i=1,nx),j=1,ny),k=1,nz) DEALLOCATE(work1) ELSE ALLOCATE(work1(nx,ny,nz),work2(nx,ny,nz)) READ(16,vasp_form) (((work1(i,j,k),i=1,nx),j=1,ny),k=1,nz) READ(16,*) READ(16,vasp_form) (((work2(i,j,k),i=1,nx),j=1,ny),k=1,nz) work1=work1-work2 WRITE(15,dmol_form) (((work1(i,j,k),i=1,nx),j=1,ny),k=1,nz) DEALLOCATE(work1,work2) END IF 从红色部分看，如果是考虑自旋（NSPIN=2）的话,你只给出了部分结果。 请问ELFCAR或者CHGCAR中自旋结果是不是这样给出的： 第一部分：total(spin up + spin down) 第二部分：spin up or spin down. 你用减法来处理是得到单纯上旋结果？ [ Last edited by cenwanglai on 2010-9-12 at 11:22 ] |
9楼2010-09-12 10:45:01
5楼2010-09-10 17:23:33
6楼2010-09-10 19:26:01
7楼2010-09-11 13:07:59
简单回复
2010-09-10 08:53
回复
顶 。
