| 查看: 4299 | 回复: 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 不知道那个对？ |
回复此楼» 猜你喜欢
河北省自然科学基金
已经有7人回复
-
揭秘青基评审内幕：几个A才能顺利中标
已经有4人回复
-
投MDPI的三个期刊都被桌拒
已经有5人回复
-
26/27申博自荐
已经有6人回复
-
青B发送上会通知了吗
已经有7人回复
-
博士申请
已经有3人回复
-
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！
已经有4人回复
-
今年审到国自然15份，谈谈感受
已经有28人回复
-
投稿求助，期刊
已经有8人回复
-
献血感触
已经有13人回复
» 本主题相关商家推荐: (我也要在这里推广)
漫天飘雪
版主 (知名作家)
认真做事，踏实做人
- 1ST强帖: 1
- 应助: 362 (硕士)
- 贵宾: 0.108
- 金币: 16449.7
- 散金: 13541
- 红花: 392
- 沙发: 2
- 帖子: 5768
- 在线: 2310.9小时
- 虫号: 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
西风月123
木虫 (著名写手)
- 应助: 15 (小学生)
- 金币: 2356.8
- 红花: 9
- 帖子: 1048
- 在线: 392小时
- 虫号: 2864449
- 注册: 2013-12-11
- 专业: 凝聚态物性 II ：电子结构
3楼2016-09-26 11:00:28
w76990480
荣誉版主 (文坛精英)
super-emuch
- 应助: 47 (小学生)
- 贵宾: 13.729
- 金币: 109757
- 散金: 29386
- 红花: 337
- 沙发: 50
- 帖子: 23155
- 在线: 3166.8小时
- 虫号: 942733
- 注册: 2010-01-15
- 专业: 固体力学
- 管辖: 文献求助
4楼2016-09-26 17:25:53
10