应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】怎么根据PROCAR作DOS和PDOS图
51/1返回列表
查看: 2572  |  回复: 21
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

huangyc

木虫 (正式写手)

[交流] 【求助】怎么根据PROCAR作DOS和PDOS图

大家好，计算DOS时得到PROCAR，准备学习画DOS和PDOS，共126个K点，每个K点100 bands，用lev00处理时提示缺少PROCAR或PROCAR错误，用p4vasp也打不开，请帮我看看，谢谢
PROCAR new format
# of k-points:  126         # of bands: 100         # of ions:  24

k-point    1 :    0.04166667 0.04166667 0.04166667     weight = 0.00462963

band   1 # energy  -29.57997136 # occ.  2.00000000

ion      s      p      d      f    tot
  1  0.233  0.000  0.000  0.000  0.233
  2  0.233  0.000  0.000  0.000  0.233
  3  0.233  0.000  0.000  0.000  0.233
  4  0.233  0.000  0.000  0.000  0.233
  5  0.000  0.000  0.000  0.000  0.000
  6  0.000  0.000  0.000  0.000  0.000
  7  0.000  0.000  0.000  0.000  0.000
  8  0.000  0.000  0.000  0.000  0.000
  9  0.000  0.000  0.000  0.000  0.000
10  0.000  0.000  0.000  0.000  0.000
11  0.000  0.000  0.000  0.000  0.000
12  0.000  0.000  0.000  0.000  0.000
13  0.000  0.000  0.000  0.000  0.000
14  0.000  0.000  0.000  0.000  0.000
15  0.000  0.000  0.000  0.000  0.000
16  0.000  0.000  0.000  0.000  0.000
17  0.000  0.000  0.000  0.000  0.000
18  0.000  0.000  0.000  0.000  0.000
19  0.000  0.000  0.000  0.000  0.000
20  0.000  0.000  0.000  0.000  0.000
21  0.000  0.000  0.000  0.000  0.000
22  0.000  0.000  0.000  0.000  0.000
23  0.000  0.000  0.000  0.000  0.000
24  0.000  0.000  0.000  0.000  0.000
tot  0.937  0.001  0.000  0.000  0.938
这是其中一部分，下面还有band2-100，接着就是k-points2-126，很长的文件
不知道大家都怎么处理DOS？学习VASP时间不长，还希望虫友们指点

[ Last edited by SHY31 on 2008-12-4 at 20:41 ]
回复此楼

» 收录本帖的淘帖专辑推荐

VASP

» 猜你喜欢
1楼 2008-11-27 22:17:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wyl04

金虫 (小有名气)
★ ★ ★
wuchenwf(金币+3,VIP+0):xiexie
implicit real*8(a-h,o-z)
      parameter (nkd =900)
      parameter (nbd = 50 )
      parameter (natmd = 20)
      parameter (ned = 1001)
      dimension dump(20),oc(nkd,nbd,natmd,4),eig_ev(nkd,nbd),wt(nkd)
      dimension ee(ned),gpdos(ned),gpdost(ned)
      dimension gpdos_s(ned),gpdos_p(ned),gpdos_d(ned)
      open(7,file='PROCAR')
      pi = 3.141592654

      write(*,*) 'Spin polarized calculation? (no=1,yes=2):'
      read (*,*) ispin
      if ((ispin.ne.1).and.(ispin.ne.2)) then
      write(*,*) ' INPUT ERROR, ispin must equal to 1 or 2 '
      stop
      endif
      if (ispin.eq.1) then
      open(31,file='ldos.dat')
      open(32,file='dos-tot.dat')
      elseif (ispin.eq.2) then
      open(31,file='ldos-up.dat')
      open(32,file='dos-tot-up.dat')
      open(41,file='ldos-dn.dat')
      open(42,file='dos-tot-dn.dat')
      endif

      write(6,*) ' Enter fermi energy: '
      read(5,*) fermi
c     write(6,102) fermi
      read(7,103) dump                    
c     write(6,103) dump                    

      write(6,*) 'which atom (na) you want to plot LDOS:'
      read(5,*) na
      write(6,*) 'enter the gaussian smearing factor:'
      read(5,*) gaussian

      do 9000 isp=1,ispin
      read(7,104) nk,nband,nion
c     write(6,104) nk,nband,nion
      if (nk .gt. nkd) stop ' nk too large '
      emin = 1000.0
      emax = -1000.0
      do 1000 k = 1,nk
      read(7,103) dump
      read(7,105) kp,pt1,pt2,pt3,wt(k)
c     write(6,105) kp,pt1,pt2,pt3,wt(k)
      read(7,103) dump
      do  nb = 1,nband
      read(7,106) nb1,eig_ev(k,nb),occ
c     write(6,106) nb1,eig_ev(k,nb),occ
      eig_ev(k,nb) = eig_ev(k,nb)-fermi
      if (eig_ev(k,nb) .gt. emax) emax = eig_ev(k,nb)
      if (eig_ev(k,nb) .lt. emin) emin = eig_ev(k,nb)
      read(7,103) dump
      read(7,103) dump
c     write(6,*) 'nion=',nion
      niont = nion +1
      if (nion .eq. 1) niont = 1
      do  ion = 1,niont
      read(7,107) (oc(k,nb,ion,j),j=1,4)
c     write(6,107) (oc(k,nb,ion,j),j=1,4)
      enddo
      read(7,103) dump
c     write(6,103) dump
      enddo
1000 continue
      weight = 0.0
      do k = 1, nk
      weight = weight + wt(k)
      enddo
      do k = 1,nk
      wt(k) =  wt(k) / weight
      enddo
c     estart_ev = - 5.0
c     eend_ev = 2.0
      estart_ev = int(emin -1.0)
      eend_ev = int(emax +1.0)
      netot =  701
      de_ev = (eend_ev - estart_ev)/(netot-1)
      do 2000 ne = 1, netot
      ee(ne) = estart_ev + (ne-1) * de_ev
      gpdos(ne) = 0.0d0
      gpdos_s(ne) = 0.0d0
      gpdos_p(ne) = 0.0d0
      gpdos_d(ne) = 0.0d0
      gpdost(ne) = 0.0d0
2000 continue
      ascal = de_ev   / ( gaussian * sqrt(pi) )
      ascal = 1.0     / ( gaussian * sqrt(pi) )
      scal_spin = 2.0
      do 5000 k=1,nk
      do 4000 nb=1,nband
      ddos = scal_spin * oc(k,nb,na,4)*wt(k)
      ddos_s = scal_spin * oc(k,nb,na,1)*wt(k)
      ddos_p = scal_spin * oc(k,nb,na,2)*wt(k)
      ddos_d = scal_spin * oc(k,nb,na,3)*wt(k)
      ddost = scal_spin *wt(k)
      do ne = 1,netot
      dij = ( eig_ev(k,nb) - ee(ne) )**2 / (gaussian**2)
      gpdos(ne) = gpdos(ne) + ascal * ddos * exp(-dij)
      gpdos_s(ne) = gpdos_s(ne) + ascal * ddos_s * exp(-dij)
      gpdos_p(ne) = gpdos_p(ne) + ascal * ddos_p * exp(-dij)
      gpdos_d(ne) = gpdos_d(ne) + ascal * ddos_d * exp(-dij)
      gpdost(ne) = gpdost(ne) + ascal * ddost * exp(-dij)
      enddo
4000 continue
5000 continue
      do 6000 ne = 1,netot
      if (gpdos(ne) .lt. 1.0E-29) gpdos(ne)=0.0
      if (gpdos_s(ne) .lt. 1.0E-29) gpdos_s(ne)=0.0
      if (gpdos_p(ne) .lt. 1.0E-29) gpdos_p(ne)=0.0
      if (gpdos_d(ne) .lt. 1.0E-29) gpdos_d(ne)=0.0
      if (gpdost(ne) .lt. 1.0E-29) gpdost(ne)=0.0
      write(21+10*isp,108) ee(ne),gpdos(ne),gpdos_s(ne),gpdos_p(ne)
     &       ,gpdos_d(ne)  
      write(22+10*isp,108) ee(ne),gpdost(ne)
6000 continue
9000 continue
101  format(10x,f9.5)
102  format(f10.5)
103  format(20a4)
104  format(16x,i3,20x,i5,19x,i4)
105  format(10x,i3,5x,3f11.8,13x,f11.8)
106  format(4x,i4,9x,f14.8,7x,f12.8)
107  format(3x,4f7.3)
108  format(f10.5,4e12.4)
      stop
      end
赞一下(10人) 回复此楼
7楼2008-12-05 18:57:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 22 个回答

intro987

铁杆木虫 (小有名气)
★ ★
jghe(金币+2,VIP+0):感谢你参与讨论！
我使用DOSCAR处理的，应该和PROCAR是一起出来的结果，里面很多列数据，第一列是能量，后面对应s, p_y, p_z, p_x, d_xy, d_yz, d_z2, d_xz, d_x2-y2轨道（VASP出来的结果应该顺序都是一样的）。数据上下顺序对应POSCAR里面的原子顺序
lev00是处理电荷密度的，p4vasp我只用来打开过POSCAR，没试过PROCAR，可能不支持这种类型吧~
赞一下(10人) 回复此楼
2楼2008-12-04 18:43:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hym5701

木虫 (小有名气)
★ ★
jghe(金币+2,VIP+0):谢谢你提供的资源！
引用回帖:
Originally posted by huangyc at 2008-11-27 22:17:
大家好，计算DOS时得到PROCAR，准备学习画DOS和PDOS，共126个K点，每个K点100 bands，用lev00处理时提示缺少PROCAR或PROCAR错误，用p4vasp也打不开，请帮我看看，谢谢
PROCAR new format
# of k-points:  126   ...

网上不是用vasp的一些tools吗，里面就有一个分析PROCAR的文件，文件名是dos-procar.f，你现将这个文件编译一下，然后就运行这个文件，他会提示你要干什么，运行完后，就可以得到totdos。dat 和ldos。dat，接着用origin等画图就可以了。
赞一下(10人) 回复此楼
3楼2008-12-04 20:22:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangyc

木虫 (正式写手)
楼上的这个东东编译不起来啊，有语法错误
哪位再指点，谢谢
赞一下 回复此楼
4楼2008-12-05 11:15:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 22 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭