24小时热门版块排行榜

返回列表

w76990480

荣誉版主 (文坛精英)

super-emuch

应助: 47 (小学生)
贵宾: 13.729
金币: 109737
散金: 29386
红花: 337
沙发: 50
帖子: 23151
在线: 3166.8小时
虫号: 942733
注册: 2010-01-15
专业: 固体力学
管辖: 文献求助

[求助] 正格矢和倒格矢之间关系已有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

不知道那个对？

回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

专业文献求助

1楼 2016-09-22 10:27:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

漫天飘雪

版主 (知名作家)

认真做事，踏实做人

1ST强帖: 1
应助: 362 (硕士)
贵宾: 0.108
金币: 16799.9
散金: 13410
红花: 392
沙发: 2
帖子: 5766
在线: 2302.5小时
虫号: 914852
注册: 2009-11-27
性别: GG
专业: 材料物理化学
管辖: 第一性原理

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +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

应该是单位的问题。