| 查看: 4396 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
w76990480荣誉版主 (文坛精英)
super-emuch
|
[求助]
正格矢和倒格矢之间关系 已有1人参与
|
|
|
vasp OUTCAR给出的倒格矢是否正确？ 下面给出的是石墨烯的 direct lattice vectors reciprocal lattice vectors 2.137655082 -1.234175737 0.000000000 0.467802317 0.000000000 0.000000000 0.000000000 2.468351474 0.000000000 0.233901158 0.405128690 0.000000000 0.000000000 0.000000000 11.999037758 0.000000000 0.000000000 0.083340016 length of vectors 2.468351474 2.468351474 11.999037758 0.467802317 0.467802317 0.083340016 如果根据计算能带的程序 band-procar.f implicit real*8(a-h,o-z) parameter (nbd = 200) parameter (nkd = 500) parameter (nxd = 300) parameter (natmd = 20) dimension a(3,3),b(3,3),c(3),e(nkd,nbd),sk(nkd,3) dimension xx(nxd) ,wei(nkd) dimension dump(20),oc(nkd,nbd,natmd,4) C open(7,file='PROCAR',form='FORMATTED',status='OLD') pi = 3.141592654 C read(7,100) dump open(8,file='POSCAR',form='FORMATTED',status='OLD') read(8,100) temp c write(6,100) temp 100 format(20a4) read (8,*) aa C WRITE(6,*) aa c c *** read lattice constant from POSCAR** c do i=1,3 read (8,*) (a(i,j),j=1,3) c WRITE(6,500) (a(i,j),j=1,3) 500 format (3f12.8) enddo do i=1,3 do j=1,3 a(i,j)=aa*a(i,j) enddo c WRITE(6,500) (a(i,j),j=1,3) enddo c c *** read lattice vector from POSCAR*** c volume=a(1,1)*a(2,2)*a(3,3)+a(1,2)*a(2,3)*a(3,1) &+a(1,3)*a(2,1)*a(3,2)-a(1,1)*a(2,3)*a(3,2) &-a(1,2)*a(2,1)*a(3,3)-a(1,3)*a(2,2)*a(3,1) do i=1,3 if (i .eq. 1) then j=2 k=3 else if (i .eq. 2) then j=3 k=1 else j=1 k=2 endif c(1)=a(j,2)*a(k,3)-a(j,3)*a(k,2) c(2)=a(j,3)*a(k,1)-a(j,1)*a(k,3) c(3)=a(j,1)*a(k,2)-a(j,2)*a(k,1) do j=1,3 b(i,j)=2*pi*c(j)/volume WRITE (6,*) b(i,j) enddo enddo END 算出来的数据为 reciprocal lattice vectors 2.93928872574 0.0000 0.00000 1.46964436287 2.545498 0.00000 0.00000000000 0.00000 5.23641 不知道那个对？ |
回复此楼» 猜你喜欢
基金系统什么内容也没有
已经有10人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有15人回复
-
我不理解！
已经有26人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有8人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
梦想
已经有8人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
西风月123
木虫 (著名写手)
- 应助: 15 (小学生)
- 金币: 2356.8
- 红花: 9
- 帖子: 1048
- 在线: 392小时
- 虫号: 2864449
- 注册: 2013-12-11
- 专业: 凝聚态物性 II ：电子结构
3楼2016-09-26 11:00:28
漫天飘雪
版主 (知名作家)
认真做事，踏实做人
- 1ST强帖: 1
- 应助: 362 (硕士)
- 贵宾: 0.108
- 金币: 16530.7
- 散金: 13945
- 红花: 392
- 沙发: 2
- 帖子: 5768
- 在线: 2314.3小时
- 虫号: 914852
- 注册: 2009-11-27
- 性别: GG
- 专业: 材料物理化学
- 管辖: 第一性原理
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
w76990480: 金币+100, ★★★★★最佳答案 2016-09-22 11:08:28
蛋蛋小童鞋: 金币+2, 谢谢参与，小木虫鼓励积极回答问题者 2016-09-22 11:43:12
感谢参与，应助指数 +1
w76990480: 金币+100, ★★★★★最佳答案 2016-09-22 11:08:28
蛋蛋小童鞋: 金币+2, 谢谢参与，小木虫鼓励积极回答问题者 2016-09-22 11:43:12
|
lattice vectors in bohr: 4.039582654 -2.332254133 0.000000000 0.000000000 4.664508265 0.000000000 0.000000000 0.000000000 22.674895121 lattice vectors in ang: 2.137655082 -1.234175737 0.000000000 0.000000000 2.468351474 0.000000000 0.000000000 0.000000000 11.999037758 reciprocal vectors in 1/bohr: 1.555404567 0.000000000 0.000000000 0.777702283 1.347019868 -0.000000000 -0.000000000 0.000000000 0.277098759 reciprocal vectors in 1/ang: 2.939288644 0.000000000 0.000000000 1.469644322 2.545498635 -0.000000000 -0.000000000 0.000000000 0.523640765 reciprocal vectors in 2pi/bohr: 0.247550325 0.000000000 0.000000000 0.123775163 0.214384870 -0.000000000 -0.000000000 0.000000000 0.044101637 reciprocal vectors in 2pi/ang: 0.467802317 0.000000000 0.000000000 0.233901158 0.405128690 -0.000000000 -0.000000000 0.000000000 0.083340016 应该是单位的问题。  |
2楼2016-09-22 10:52:16
w76990480
荣誉版主 (文坛精英)
super-emuch
- 应助: 47 (小学生)
- 贵宾: 13.729
- 金币: 109764
- 散金: 29386
- 红花: 337
- 沙发: 50
- 帖子: 23156
- 在线: 3166.8小时
- 虫号: 942733
- 注册: 2010-01-15
- 专业: 固体力学
- 管辖: 文献求助
4楼2016-09-26 17:25:53
30